CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

Discretization of transient stokes equation using finite element method

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 25, 2021, 16:05
Default Discretization of transient stokes equation using finite element method
  #1
New Member
 
Join Date: Dec 2021
Posts: 1
Rep Power: 0
cfdplus is on a distinguished road
Hello,
Before everything I need to apologize for not writing the equation in a suitable format, with a quick search I came to know that I have to use Latex but unfortunately I have no experience with Latex.
I am trying to solve the transient stokes equation using finite element method. I discretized the equation using triangular element and ended up with below equation.

Stokes equation:

∂u/∂t=ϑ∆u-∂P/ρ

velocity x-direction:
(LNN+dt*nu*DNDN)*(u_x at time step n+1)=LNN*(u_x at time step n)+dt/rho(transpose of NDNx)*(Pressure at time step n+1)

velocity y-direction:
(LNN+dt*nu*DNDN)*(u_y at time step n+1)=LNN*(u_y at time step n)+dt/rho(transpose of NDNy)*(Pressure at time step n+1)+g*LNN

pressure equation:

DNDN*(pressure at time step n+1)=-(rho/dt)*(transpose of NDNx* (u_x at time step n+1)+transpose NDNy*(u_y at time step n+1)) + (Transpose NDNx* gradient P in x direction)+(transpose NDNy*gradient P in y direction))

I used the dt=0.001, nu= 1e-3, density=1000 and gravity=10. The problem is I get very large value for pressure and I'm notsure whether the discretized equations are correct ot not. The initial value for ux=2, uy=0, pressure=0. Imagine the domain is between -1<x<1, -1<y<1. The boundary condition is ux=2 for x=-1, uy=0 everywhere and p=1 at x=1.



LNN: lumped matrix of mass matrix
dt: time
nu:kinematic viscosity
rho:density
g:gravity
N : finite element shape function
Nx: derivation of shape function with respect to x
Ny: derivation of shape function with respect to y
NN:Mass Matrix
NDNx=Matrix of NNx in finite element
NDNy=Matrix of NNy in finite element
DNDN: Laplacian Matrix =-(NxNx+NyNy)
cfdplus is offline   Reply With Quote

Reply

Tags
finite element method, fluid dynamics

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
What is C.V. based finite element method C-H Kuo Main CFD Forum 4 September 19, 2022 15:06
Domain Reference Pressure and mass flow inlet boundary AdidaKK CFX 75 August 20, 2018 06:37
[Other] Mesh Importing Problem cuteapathy ANSYS Meshing & Geometry 2 June 24, 2017 06:29
The finite element method - “displacement” or “stiffness” method 3Lin7 ANSYS 0 February 15, 2013 19:36
Finite element method - number of elements poornima Main CFD Forum 6 April 14, 2006 08:23


All times are GMT -4. The time now is 08:24.