Finite element solution of the Navier–Stokes equations by a velocity–vorticity method

A velocity–vorticity formulation of the Navier–Stokes equations is presented as an alternative to the primitive variables approach. The velocity components and the vorticity are solved for in a fully coupled manner using a Newton method. No artificial viscosity is required in this formulation. The pressure is updated by a method allowing natural imposition of boundary conditions. Incompressible and subsonic results are presented for two-dimensional laminar internal flows up to high Reynolds numbers.     

Three-dimensional incompressible Navier–Stokes equations on non-orthogonal staggered grids using the velocity–vorticity formulation

This paper is concerned with the numerical resolution of the incompressible Navier–Stokes equations in the velocity–vorticity form on non-orthogonal structured grids. The discretization is performed in such a way, that the discrete operators mimic the properties of the continuous ones. This allows the discrete equivalence between the primitive and velocity–vorticity formulations to be proved. This last formulation can thus be seen as a particular technique for solving the primitive equations. The difficulty associated with non-simply connected computational domains and with the implementation of the boundary conditions are discussed. One of the main drawback of the velocity–vorticity formulation, relative to the additional computational work required for solving the additional unknowns, is alleviated. Two- and three-dimensional numerical test cases validate the proposed method. © 1998 John Wiley & Sons, Ltd.     

A new velocity–vorticity boundary integral formulation for Navier–Stokes equations

In the present work, an indirect boundary integral method for the numerical solution of Navier–Stokes equations formulated in velocity–vorticity dependent variables is proposed. This wholly integral approach, based on Helmholtz's decomposition, deals directly with the vorticity field and gives emphasis to the establishment of appropriate boundary conditions for the vorticity transport equation. The coupling between the vorticity and the vortical velocity fields is expressed by an iterative procedure. The present analysis shows the usefulness of an integral formulation not only in providing a potentially more efficient computational tool, but also in giving a better understanding to the physics of the phenomenon. Copyright © 2000 John Wiley & Sons, Ltd.     

A comparison of primitive variables and streamfunction–vorticity for fem depth-averaged modelling of steady flow

The governing equations for depth-averaged turbulent flow are presented in both the primitive variable and streamfunction–vorticity forms. Finite element formulations are presented, with special emphasis on the handling of bottom stress terms and spatially varying eddy viscosity. The primitive variable formulation is found to be preferable because of its flexibility in handling spatial variation in viscosity, variability in water surface elevations, and inflow and outflow boundaries. The substantial reduction in computational effort afforded by the streamfunction–vorticity formulation is found not to be sufficient to recommend its use for general depth-averaged flows. For those flows in which the surface can be approximated as a fixed level surface, the streamfunction–vorticity form can produce results equivalent to the primitive variable form as long as turbulent viscosity can be estimated as a constant.     

A finite element formulation for laminar flow of a fluid with microstructure

A finite element formulation for the steady laminar flow of an incompressible fluid with microstructure has been developed. The particular fluids considered are commonly known as micropolar fluids, in which case suspended particulate microstructures are modelled by an ‘extended’ continuum formulation. The particle microspin is a new kinematic variable which is independent of the classical vorticity vector and thereby allows relative rotation between particles and the surrounding fluid. This formulation also gives rise to couple stresses in addition to classical force or traction stresses. The finite element formulation utilizes a variational approach and imposes conservation of mass through a penalty function. A general boundary condition for microspin has been incorporated whereby microspin at a solid boundary is constrained to be proportional to the fluid vorticity. The proportionality constant in this case can vary from zero to unity. Sample solutions are presented for fully developed flow through a straight tube and compared with an analytical solution. Results are also generated for flow through a constricted tube and compared with a Newtonian fluid solution.     

A stream function–vorticity formulation‐based immersed boundary method and its applications

A new stream function–vorticity formulation‐based immersed boundary method is presented in this paper. Different from the conventional immersed boundary method, the main feature of the present model is to accurately satisfy both governing equations and boundary conditions through velocity correction and vorticity correction procedures. The velocity correction process is performed implicitly based on the requirement that velocity at the immersed boundary interpolated from the corrected velocity field accurately satisfies the nonslip boundary condition. The vorticity correction is made through the stream function formulation rather than the vorticity transport equation. It is evaluated from the firstorder derivatives of velocity correction. Two simple and efficient ways are presented for approximation of velocity‐correction derivatives. One is based on finite difference approximation, while the other is based on derivative expressions of Dirac delta function and velocity correction. It was found that both ways can work very well. The main advantage of the proposed method lies in its simple concept, easy implementation, and robustness in stability. Numerical experiments for both stationary and moving boundary problems were conducted to validate the capability and efficiency of the present method. Good agreements with available data in the literature were achieved. Copyright © 2011 John Wiley & Sons, Ltd.     

Boundary element method for steady 2D high-Reynolds-number flow

This paper deals with the numerical simulation of fluid dynamics using the boundary–domain integral technique (BEM). The steady 2D diffusion–convection equations are discussed and applied to solve the plane Navier-Stokes equations. A vorticity–velocity formulation has been used. The numerical scheme was tested on the well-known ‘driven cavity’ problem. Results for Re = 1000 and 10,000 are compared with benchmark solutions. There are also results for Re = 15,000 but they have only qualitative value. The purpose was to show the stability and robustness of the method even when the grid is relatively coarse.     

Local vorticity computation approach in double distribution functions based lattice Boltzmann methods for flow and scalar transport

Computation of vorticity, or the skew-symmetric velocity gradient tensor, in conjunction with the strain rate tensor, plays an important role in the flow classification, in vortical structure identification and in the modeling of various complex fluids and flows. For the simulation of flows accompanied by the advection-diffusion transport of a scalar field (e.g., temperature), double distribution functions (DDF) based lattice Boltzmann (LB) methods, involving a pair of LB schemes are commonly used. We present a new local vorticity computation approach by introducing an intensional anisotropy of the scalar flux in the third order, off-diagonal moment equilibria of the LB scheme for the scalar field, and then combining the second order non-equilibrium components of both the LB methods. As such, any pair of lattice sets in the DDF formulation that can independently support the third order off-diagonal moments would enable local determination of the complete flow kinematics, with the LB methods for the fluid motion and the transport of the passive scalar respectively providing the necessary moment relationships to determine the symmetric and skew-symmetric components of the velocity gradient tensor. Since the resulting formulation is completely local and do not rely on any finite difference approximations for velocity derivatives, it is by design naturally suitable for parallel computation. As an illustration of our approach, we formulate a DDF-LB scheme for local vorticity computation using a pair of multiple relaxation times (MRT) based collision approaches on two-dimensional, nine velocity (D2Q9) lattices, where the necessary moment relationships to determine the velocity gradient tensor and the vorticity are established via a Chapman-Enskog analysis. Simulations of various benchmark flows demonstrate good accuracy of the predicted vorticity fields using our approach against available solutions, including numerical results, with a second order convergence. Furthermore, extensions of our formulation for a variety of collision models, including those based on cascaded and non-cascaded central moments, to enable local vorticity computation are presented.     

Finite volume solution of the Navier–Stokes equations in velocity–vorticity formulation

For the incompressible Navier–Stokes equations, vorticity‐based formulations have many attractive features over primitive‐variable velocity–pressure formulations. However, some features interfere with the use of the numerical methods based on the vorticity formulations, one of them being the lack of a boundary conditions on vorticity. In this paper, a novel approach is presented to solve the velocity–vorticity integro‐differential formulations. The general numerical method is based on standard finite volume scheme. The velocities needed at the vertexes of each control volume are calculated by a so‐called generalized Biot–Savart formula combined with a fast summation algorithm, which makes the velocity boundary conditions implicitly satisfied by maintaining the kinematic compatibility of the velocity and vorticity fields. The well‐known fractional step approaches are used to solve the vorticity transport equation. The paper describes in detail how we accurately impose no normal‐flow and no tangential‐flow boundary conditions. We impose a no‐flux boundary condition on solid objects by the introduction of a proper amount of vorticity at wall. The diffusion term in the transport equation is treated implicitly using a conservative finite update. The diffusive fluxes of vorticity into flow domain from solid boundaries are determined by an iterative process in order to satisfy the no tangential‐flow boundary condition. As application examples, the impulsively started flows through a flat plate and a circular cylinder are computed using the method. The present results are compared with the analytical solution and other numerical results and show good agreement. Copyright © 2005 John Wiley & Sons, Ltd.     

Solution of the Navier–Stokes equations in velocity–vorticity form using a Eulerian–Lagrangian boundary element method

This paper describes the Eulerian–Lagrangian boundary element model for the solution of incompressible viscous flow problems using velocity–vorticity variables. A Eulerian–Lagrangian boundary element method (ELBEM) is proposed by the combination of the Eulerian–Lagrangian method and the boundary element method (BEM). ELBEM overcomes the limitation of the traditional BEM, which is incapable of dealing with the arbitrary velocity field in advection‐dominated flow problems. The present ELBEM model involves the solution of the vorticity transport equation for vorticity whose solenoidal vorticity components are obtained iteratively by solving velocity Poisson equations involving the velocity and vorticity components. The velocity Poisson equations are solved using a boundary integral scheme and the vorticity transport equation is solved using the ELBEM. Here the results of two‐dimensional Navier–Stokes problems with low–medium Reynolds numbers in a typical cavity flow are presented and compared with a series solution and other numerical models. The ELBEM model has been found to be feasible and satisfactory. Copyright © 2000 John Wiley & Sons, Ltd.     

The (9,5) HOC formulation for the transient Navier–Stokes equations in primitive variable

We recently proposed an improved (9,5) higher order compact (HOC) scheme for the unsteady two‐dimensional (2‐D) convection–diffusion equations. Because of using only five points at the current time level in the discretization procedure, the scheme was seen to be computationally more efficient than its predecessors. It was also seen to capture very accurately the solution of the unsteady 2‐D Navier–Stokes (N–S) equations for incompressible viscous flows in the stream function–vorticity (ψ – ω) formulation. In this paper, we extend the scope of the scheme for solving the unsteady incompressible N–S equations based on primitive variable formulation on a collocated grid. The parabolic momentum equations are solved for the velocity field by a time‐marching strategy and the pressure is obtained by discretizing the elliptic pressure Poisson equation by the steady‐state form of the (9,5) scheme with the Neumann boundary conditions. In particular, for pressure, we adopt a strategy on the collocated grid in conjunction with ideas borrowed from the staggered grid approach in finite volume. We first apply this extension to a problem having analytical solution and then to the famous lid‐driven square cavity problem. We also apply our formulation to the backward‐facing step problem to see how the method performs for external flow problems. The results are presented and are compared with established numerical results. This new approach is seen to produce excellent comparison in all the cases. Copyright © 2007 John Wiley & Sons, Ltd.     

THE NO-SLIP BOUNDARY CONDITION IN FINITE DIFFERENCE APPROXIMATIONS

A finite difference method for the Navier–Stokes equations in vorticity –streamfunction formulation is proposed to resolve the difficulty of the lack of a vorticity boundary condition at a no-slip boundary. It is particularly suitable for flows in regions with complicated geometries. Convergence with second-order accuracy in vorticity and velocity is established. In numerical experiments the convergence rates agree with theoretical predictions. Test results for the two-dimensional driven cavity problem and for the flow in expansion and contraction channels are given.     

A spectral solver for the Navier–Stokes equations in the velocity–vorticity formulation

A numerical algorithm intended for the study of flows in a cylindrical container under laminar flow conditions is proposed. High resolution of the flow field, governed by the Navier–Stokes equations in velocity–vorticity formulation relative to a cylindrical frame of reference, is achieved through spatial discretisation by means of the spectral method. This method is based on a Fourier expansion in the azimuthal direction and an expansion in Chebyshev polynomials in the (nonperiodic) radial and axial directions. Several regularity constraints are used to take care of the coordinate singularity. These constraints are implemented, together with the boundary conditions at the top, bottom and mantle of the cylinder, via the tau method. The a priori unknown boundary values of the vorticity are evaluated by means of the influence-matrix technique. The compatibility between the mathematical and numerical formulation of the Navier–Stokes equations is established through a tau-correction procedure. The resolved flow field exhibits high-precision satisfaction of the incompressibility constraints for velocity and vorticity and the definition of the vorticity. The performance of the solver is illustrated by resolution of several configurations representative of generic three-dimensional laminar flows.     

Integral vorticity transport method on an overset grid

We present an overset grid method for solution of the integro‐differential vorticity–velocity formulation of the Navier–Stokes equations for two‐dimensional, incompressible flow. The method uses a body‐fitted inner grid, on which vorticity is evolved semi‐implicitly, and a Cartesian outer grid with explicit vorticity evolution. The Biot–Savart integral is solved using an adaptive, optimized multipole acceleration method. The Biot–Savart integration is performed over all inner grid cells, over all ‘active cells’ of the outer grid that lie entirely outside of the inner grid, and over sub‐elements of a set of ‘overhanging’ cells of the outer grid that overlap part of the inner grid. A novel method is developed using a level‐set distance function to rapidly and easily partition the overhanging grid cells, which is essential for the Biot–Savart integration in order to avoid double‐counting vorticity in the overhanging region. A similar decomposition into outer, inner and overhanging cells is used in solving for pressure using a boundary‐element formulation, which requires evaluation of an integral over the vorticity field using a method similar to that used for the Biot–Savart integral. The new overset grid method is applied to flow past stationary and moving bodies in two dimensions and found to agree well with prior experimental and numerical results. Copyright © 2011 John Wiley & Sons, Ltd.     

A 3D incompressible Navier–Stokes velocity–vorticity weak form finite element algorithm

The velocity–vorticity formulation is selected to develop a time‐accurate CFD finite element algorithm for the incompressible Navier–Stokes equations in three dimensions.The finite element implementation uses equal order trilinear finite elements on a non‐staggered hexahedral mesh. A second order vorticity kinematic boundary condition is derived for the no slip wall boundary condition which also enforces the incompressibility constraint. A biconjugate gradient stabilized (BiCGSTAB) sparse iterative solver is utilized to solve the fully coupled system of equations as a Newton algorithm. The solver yields an efficient parallel solution algorithm on distributed‐memory machines, such as the IBM SP2. Three dimensional laminar flow solutions for a square channel, a lid‐driven cavity, and a thermal cavity are established and compared with available benchmark solutions. Copyright © 2002 John Wiley & Sons, Ltd.     

Numerical simulation of unsteady interaction of centrifugal impeller with its diffuser using Lagrangian discrete vortex method

In this study, an advanced Lagrangian vortex- boundary element method is applied to simulate the unsteady impeller-diffuser interactions in a diffuser pump not only for design but also for off-design considerations. In velocity calculations based on the Biot-Savart law we do not have to grid large portions of the flow field and the calculation points are concentrated in the regions where vorticity is present. Lagrangian representation of the evolving vorticity field is well suited to moving boundaries. An integral pressure equation shows that the pressure distribution can be estimated directly from the instantaneous velocity and vorticity field. The numerical results are compared with the experimental data and the comparisons show that the method used in this study can provide us insight into the complicated unsteady impeller-diffuser interaction phenomena in a diffuser pump.     

A short note on the reattachment length for BFS problem

The standard backward‐facing step flow problem is solved for steady state laminar case using stream function‐vorticity method. The steady state results are obtained as the asymptotic solution of the transient formulation. The primary reattachment length is studied and the discrepancy in the v velocity is reported. A method for determining appropriate locations for comparison is proposed. The energy equation is solved and found to be in good agreement with benchmark results. Copyright © 2005 John Wiley & Sons, Ltd.     

A numerical‐variational procedure for laminar flow in curved square ducts

A new numerical method is presented for the solution of the Navier–Stokes and continuity equations governing the internal incompressible flows. The method denoted as the CVP method consists in the numerical solution of these equations in conjunction with three additional variational equations for the continuity, the vorticity and the pressure field, using a non‐staggered grid. The method is used for the study of the characteristics of the laminar fully developed flows in curved square ducts. Numerical results are presented for the effects of the flow parameters like the curvature, the Dean number and the stream pressure gradient on the velocity distributions, the friction factor and the appearance of a pair of vortices in addition to those of the familiar secondary flow. The accuracy of the method is discussed and the results are compared with those obtained by us, using a variation of the velocity–pressure linked equation methods denoted as the PLEM method and the results obtained by other methods. Copyright © 2004 John Wiley & Sons, Ltd.     

A fluid solver based on vorticity–helical density equations with application to a natural convection in a cubic cavity

We study numerically a recently introduced formulation of incompressible Newtonian fluid equations in vorticity–helical density and velocity–Bernoulli pressure variables. Unlike most numerical methods based on vorticity equations, the current approach provides discrete solutions with mass conservation, divergence‐free vorticity, and accurate kinetic energy balance in a simple and natural way. The method is applied to compute buoyancy‐driven flows in a differentially heated cubic enclosure in the Boussinesq approximation for Ra ∈ {10⁴,10⁵,10⁶}. The numerical solutions on a finer grid are of benchmark quality. The computed helical density allows quantification of the three‐dimensional nature of the flow. Copyright © 2011 John Wiley & Sons, Ltd.     

A hybrid vortex method for the simulation of three‐dimensional flows

This paper presents an integral vorticity method for solving three‐dimensional Navier–Stokes equations. A finite volume scheme is implemented to solve the vorticity transport equation, which is discretized on a structured hexahedral mesh. A vortex sheet algorithm is used to enforce the no‐slip boundary condition through a vorticity flux at the boundary. The Biot–Savart integral is evaluated to compute the velocity field, in conjunction with a fast algorithm based on multipole expansion. This method is applied to the simulation of uniform flow past a sphere. Copyright © 2007 John Wiley & Sons, Ltd.