Kinematic Laplacian Equation Method: A Velocity-Vorticity Formulation for the Navier-Stokes Equations

In this work, a novel procedure to solve the Navier-Stokes equations in the vorticity-velocity formulation is presented. The vorticity transport equation is solved as an ordinary differential equation (ODE) problem on each node of the spatial discretization. Evaluation of the right-hand side of the ODE system is computed from the spatial solution for the velocity field provided by a new partial differential equation expression called the kinematic Laplacian equation (KLE). This complete decoupling of the two variables in a vorticity-in-time/velocity-in-space split algorithm reduces the number of unknowns to solve in the time-integration process and also favors the use of advanced ODE algorithms, enhancing the efficiency and robustness of time integration. The issue of the imposition of vorticity boundary conditions is addressed, and details of the implementation of the KLE by isoparametric finite element discretization are given. Validation results of the KLE method applied to the study of the classical case of a circular cylinder in impulsive-started pure-translational steady motion are presented. The problem is solved at several Reynolds numbers in the range 5 < Re < 180 comparing numerical results with experimental measurements and flow visualization plates. Finally, a recent result from a study on periodic vortex-array structures produced in the wake of forced-oscillating cylinders is included.

[1]  Sadatoshi Taneda,et al.  Experimental Investigation of the Wakes behind Cylinders and Plates at Low Reynolds Numbers , 1956 .

[2]  Alexandre J. Chorin,et al.  Vortex sheet approximation of boundary layers , 1978 .

[3]  C. Davies,et al.  A novel velocity-vorticity formulation of the Navier-Stokes equations with applications to boundary layer disturbance evolution , 2001 .

[4]  C. Williamson.  Oblique and parallel modes of vortex shedding in the wake of a circular cylinder at low Reynolds numbers , 1989, Journal of Fluid Mechanics.

[5]  A. Chorin.  Numerical study of slightly viscous flow , 1973, Journal of Fluid Mechanics.

[6]  Charles G. Speziale,et al.  On the advantages of the vorticity-velocity formulations of the equations of fluid dynamics , 1986 .

[7]  F. Stella,et al.  A vorticity-velocity method for the numerical solution of 3D incompressible flows , 1993 .

[8]  William H. Press,et al.  Numerical recipes in C , 2002 .

[9]  L. Quartapelle,et al.  Projection conditions on the vorticity in viscous incompressible flows , 1981 .

[10]  Derek B. Ingham,et al.  Finite-difference methods for calculating steady incompressible flows in three dimensions , 1979 .

[11]  G. Pascazio,et al.  A numerical method for the vorticity-velocity Navier-Stokes equations in two and three dimensions , 1991 .

[12]  Fernando L. Ponta,et al.  Constant-curl Laplacian equation: a new approach for the analysis of flows around bodies , 2003 .

[13]  Hjh Herman Clercx,et al.  A Spectral Solver for the Navier-Stokes Equations in the Velocity-Vorticity Formulation for Flows with Two Nonperiodic Directions , 1997 .

[14]  Fernando L. Ponta,et al.  A vortex model for Darrieus turbine using finite element techniques , 2001 .

[15]  D. L. Young,et al.  Arbitrary Lagrangian-Eulerian finite element analysis of free surface flow using a velocity-vorticity formulation , 2004 .

[16]  M. V. Dyke,et al.  An Album of Fluid Motion , 1982 .

[17]  G. Batchelor,et al.  An Introduction to Fluid Dynamics , 1968 .

[18]  L. Quartapelle,et al.  Numerical solution of the incompressible Navier-Stokes equations , 1993, International series of numerical mathematics.

[19]  L. Quartapelle,et al.  Vorticity conditioning in the computation of two-dimensional viscous flows , 1981 .

[20]  Wagdi G. Habashi,et al.  Finite Element Solution of the 3D Compressible Navier-Stokes Equations by a Velocity-Vorticity Method , 1993 .

[21]  Thomas B. Gatski,et al.  The numerical solution of the Navier-Stokes equations for 3-dimensional, unsteady, incompressible flows by compact schemes , 1989 .

[22]  K. Bathe.  Finite Element Procedures , 1995 .

[23]  Hermann F. Fasel,et al.  Investigation of the stability of boundary layers by a finite-difference model of the Navier—Stokes equations , 1976, Journal of Fluid Mechanics.