Solution of 2-D Incompressible Navier Stokes Equations with Artificial Compressibility Method Using Ftcs Scheme
Department of Mechanical Engineering College of EME
National University of Science and Technology
Abstract— The paper deals with the 2-D lid-driven cavity flow governed by the non dimensional incompressible Navier-Stokes theorem in the rectangular domain. Specific boundary conditions for this case study have been defined and the flow characteristics pertaining to the scenario have been coded in MATLAB using artificial compressibility method and FTCS scheme. The results are compared successfully with an authentic research paper by Ghia, Ghia & Shin.
Keywords: Navier …show more content…
pi,jn+1-pi,jn∆t+1/β2[ ui+1,jn-ui-1,jn2∆x+vi,j+1n-vi,j-1n2∆y ]= 0
Once the pressure component is computed for the time level n+1, the velocity equations can be solved by using the pressure value calculated above as
Discretized U equation:
=(1/Re)[ ui+1,jn-2ui,j+ui-1,jn∆x2+ui+1,jn-2ui,j+ui-1,jn∆y2 ]
Discretized V equation:
=(1/Re)[ vi+1,jn-2vi,j+vi-1,jn∆x2+vi+1,jn-2vi,j+vi-1,jn∆y2 ]
Boundary conditions for Pressure:
The pressure at the boundary conditions will come from the pressure at the internal points by using the equation approximated as Neuman boundary condition on all four walls. For pressure, Neuman boundary conditions are applied at top, bottom and left, right walls. The normal derivative ∂p∂n is given by two points having boundary in the middle .
Pi,j+1-Pi,jdy=0 => Pi,j+1=Pi,j
As the point Pi,j+1 is out of domain so it is approximated as Pi,j-1.
The vorticity at a fluid point is defined as twice the angular velocity and is mathematically expressed as
For two-dimensional flow, it reduced to
In descritized form
The vorticity equation is solved in Matlab in order to see the rotational effect of flow particles.