TY - JOUR
T1 - Numerical solution of the incompressible Navier-Stokes equations in primitive variables and velocity-vorticity formulation
AU - Loukopoulos, V. C.
AU - Messaris, G. T.
AU - Bourantas, G. C.
PY - 2013
Y1 - 2013
N2 - A finite-difference method is presented for the numerical solution of the Navier-Stokes equations of motion of a viscous incompressible fluid in two dimensions in primitive-variables and velocity-vorticity formulation. For the case of primitive-variables, using a staggered grid and introducing an auxiliary function of the coordinate system and considering the form of the initial equation on lines passing through the nodal point (x0, y0) and parallel to the coordinate axes, we can separate it into two parts that are finally reduced to ordinary linear differential equations, one for each dimension. Discretization of these equations leads to a system of linear equations in n-unknowns which is solved by an iterative technique and the method converges rapidly giving satisfactory results. For the pressure variable we consider a pressure Poisson equation with suitable Neumann type boundary conditions. In case of velocity-vorticity formulation similar procedure is used, for a collocated grid. Numerical results up to Reynolds number 5000, confirming the accuracy of the proposed method, are presented for configurations of interest, such as Poiseuille flow, lid-driven cavity flow in square domain and flow in a backward-facing step in the presence of a pressure gradient.
AB - A finite-difference method is presented for the numerical solution of the Navier-Stokes equations of motion of a viscous incompressible fluid in two dimensions in primitive-variables and velocity-vorticity formulation. For the case of primitive-variables, using a staggered grid and introducing an auxiliary function of the coordinate system and considering the form of the initial equation on lines passing through the nodal point (x0, y0) and parallel to the coordinate axes, we can separate it into two parts that are finally reduced to ordinary linear differential equations, one for each dimension. Discretization of these equations leads to a system of linear equations in n-unknowns which is solved by an iterative technique and the method converges rapidly giving satisfactory results. For the pressure variable we consider a pressure Poisson equation with suitable Neumann type boundary conditions. In case of velocity-vorticity formulation similar procedure is used, for a collocated grid. Numerical results up to Reynolds number 5000, confirming the accuracy of the proposed method, are presented for configurations of interest, such as Poiseuille flow, lid-driven cavity flow in square domain and flow in a backward-facing step in the presence of a pressure gradient.
KW - LFDM
KW - Navier-Stokes equations
KW - Numerical solution
KW - Primitive variables
KW - Staggered grid
KW - Velocity-Vorticity
UR - http://www.scopus.com/inward/record.url?scp=84883179525&partnerID=8YFLogxK
U2 - 10.1016/j.amc.2013.08.002
DO - 10.1016/j.amc.2013.08.002
M3 - Article
AN - SCOPUS:84883179525
SN - 0096-3003
VL - 222
SP - 575
EP - 588
JO - Applied Mathematics and Computation
JF - Applied Mathematics and Computation
ER -