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.     

Solution to transient Navier–Stokes equations by the coupling of differential quadrature time integration scheme with dual reciprocity boundary element method

The two‐dimensional time‐dependent Navier–Stokes equations in terms of the vorticity and the stream function are solved numerically by using the coupling of the dual reciprocity boundary element method (DRBEM) in space with the differential quadrature method (DQM) in time. In DRBEM application, the convective and the time derivative terms in the vorticity transport equation are considered as the nonhomogeneity in the equation and are approximated by radial basis functions. The solution to the Poisson equation, which links stream function and vorticity with an initial vorticity guess, produces velocity components in turn for the solution to vorticity transport equation. The DRBEM formulation of the vorticity transport equation results in an initial value problem represented by a system of first‐order ordinary differential equations in time. When the DQM discretizes this system in time direction, we obtain a system of linear algebraic equations, which gives the solution vector for vorticity at any required time level. The procedure outlined here is also applied to solve the problem of two‐dimensional natural convection in a cavity by utilizing an iteration among the stream function, the vorticity transport and the energy equations as well. The test problems include two‐dimensional flow in a cavity when a force is present, the lid‐driven cavity and the natural convection in a square cavity. The numerical results are visualized in terms of stream function, vorticity and temperature contours for several values of Reynolds (Re) and Rayleigh (Ra) numbers. Copyright © 2008 John Wiley & Sons, Ltd.     

Vorticity–velocity formulation of the 3D Navier–Stokes equations in cylindrical co‐ordinates

A finite difference method is presented for solving the 3D Navier–Stokes equations in vorticity–velocity form. The method involves solving the vorticity transport equations in ‘curl‐form’ along with a set of Cauchy–Riemann type equations for the velocity. The equations are formulated in cylindrical co‐ordinates and discretized using a staggered grid arrangement. The discretized Cauchy–Riemann type equations are overdetermined and their solution is accomplished by employing a conjugate gradient method on the normal equations. The vorticity transport equations are solved in time using a semi‐implicit Crank–Nicolson/Adams–Bashforth scheme combined with a second‐order accurate spatial discretization scheme. Special emphasis is put on the treatment of the polar singularity. Numerical results of axisymmetric as well as non‐axisymmetric flows in a pipe and in a closed cylinder are presented. Comparison with measurements are carried out for the axisymmetric flow cases. Copyright © 2003 John Wiley & Sons, Ltd.     

A fourth‐order compact finite difference scheme for the steady stream function–vorticity formulation of the Navier–Stokes/Boussinesq equations

A fourth‐order compact finite difference scheme on the nine‐point 2D stencil is formulated for solving the steady‐state Navier–Stokes/Boussinesq equations for two‐dimensional, incompressible fluid flow and heat transfer using the stream function–vorticity formulation. The main feature of the new fourth‐order compact scheme is that it allows point‐successive overrelaxation (SOR) or point‐successive underrelaxation iteration for all Rayleigh numbers Ra of physical interest and all Prandtl numbers Pr attempted. Numerical solutions are obtained for the model problem of natural convection in a square cavity with benchmark solutions and compared with some of the accurate results available in the literature. Copyright © 2003 John Wiley & Sons, Ltd.     

Effect of boundary vorticity discretization on explicit stream‐function vorticity calculations

The numerical solution of the time‐dependent Navier–Stokes equations in terms of the vorticity and a stream function is a well tested process to describe two‐dimensional incompressible flows, both for fluid mixing applications and for studies in theoretical fluid mechanics. In this paper, we consider the interaction between the unsteady advection–diffusion equation for the vorticity, the Poisson equation linking vorticity and stream function and the approximation of the boundary vorticity, examining from a practical viewpoint, global iteration stability and error. Our results show that most schemes have very similar global stability constraints although there may be small stability gains from the choice of method to determine boundary vorticity. Concerning accuracy, for one model problem we observe that there were cases where the boundary vorticity discretization did not propagate to the interior, but for the usual cavity flow all the schemes tested had error close to second order. Copyright © 2005 John Wiley & Sons, Ltd.     

Fourth‐ and tenth‐order compact finite difference solutions of perturbed circular vortex flows

In this study, high‐order compact finite difference calculations are reported for 2D unsteady incompressible circular vortex flow in primitive variable formulation. The fourth‐order Runge–Kutta temporal discretization is used together with fourth‐ or tenth‐order compact spatial discretization. Dependent on the perturbation initially imposed, the solutions display a tripole, triangular or square vortex. The comparison of the predictions with the detailed spectral calculations of Kloosterziel and Carnevale (J. Fluid Mech. 1999; 388 :217–257) shows that the vorticity fields are very well captured. The spectral resolution of the present method was quantified from the decomposition of the vorticity distribution in its azimuthal components and compared with reported spectral results. Using identical grid resolution to the reference results yields negligible differences in the main features of the flow. The perturbation amplitude and its first harmonic are virtually identical to the reference results for both fourth‐ or tenth‐order spatial discretization, as theoretically expected but seldom a posteriori verified. The differences between the two spatial discretizations appear only for coarser grids, favouring the tenth‐order discretization. Copyright © 2005 John Wiley & Sons, Ltd.     

The advantages of using high‐order finite differences schemes in laminar–turbulent transition studies

This paper presents various finite difference schemes and compare their ability to simulate instability waves in a given flow field. The governing equations for two‐dimensional, incompressible flows were solved in vorticity–velocity formulation. Four different space discretization schemes were tested, namely, a second‐order central differences, a fourth‐order central differences, a fourth‐order compact scheme and a sixth‐order compact scheme. A classic fourth‐order Runge–Kutta scheme was used in time. The influence of grid refinement in the streamwise and wall normal directions were evaluated. The results were compared with linear stability theory for the evolution of small‐amplitude Tollmien–Schlichting waves in a plane Poiseuille flow. Both the amplification rate and the wavenumber were considered as verification parameters, showing the degree of dissipation and dispersion introduced by the different numerical schemes. The results confirmed that high‐order schemes are necessary for studying hydrodynamic instability problems by direct numerical simulation. Copyright © 2005 John Wiley & Sons, Ltd.     

A non‐conforming least‐squares finite element method for incompressible fluid flow problems

In this paper, we develop least‐squares finite element methods (LSFEMs) for incompressible fluid flows with improved mass conservation. Specifically, we formulate a new locally conservative LSFEM for the velocity–vorticity–pressure Stokes system, which uses a piecewise divergence‐free basis for the velocity and standard C⁰ elements for the vorticity and the pressure. The new method, which we term dV‐VP improves upon our previous discontinuous stream‐function formulation in several ways. The use of a velocity basis, instead of a stream function, simplifies the imposition and implementation of the velocity boundary condition, and eliminates second‐order terms from the least‐squares functional. Moreover, the size of the resulting discrete problem is reduced because the piecewise solenoidal velocity element is approximately one‐half of the dimension of a stream‐function element of equal accuracy. In two dimensions, the discontinuous stream‐function LSFEM [1] motivates modification of our functional, which further improves the conservation of mass. We briefly discuss the extension of this modification to three dimensions. Computational studies demonstrate that the new formulation achieves optimal convergence rates and yields high conservation of mass. We also propose a simple diagonal preconditioner for the dV‐VP formulation, which significantly reduces the condition number of the LSFEM problem. Published 2012. This article is a US Government work and is in the public domain in the USA.     

Numerical solution of three‐dimensional velocity–vorticity Navier–Stokes equations by finite difference method

This paper describes the finite difference numerical procedure for solving velocity–vorticity form of the Navier–Stokes equations in three dimensions. The velocity Poisson equations are made parabolic using the false‐transient technique and are solved along with the vorticity transport equations. The parabolic velocity Poisson equations are advanced in time using the alternating direction implicit (ADI) procedure and are solved along with the continuity equation for velocities, thus ensuring a divergence‐free velocity field. The vorticity transport equations in conservative form are solved using the second‐order accurate Adams–Bashforth central difference scheme in order to assure divergence‐free vorticity field in three dimensions. The velocity and vorticity Cartesian components are discretized using a central difference scheme on a staggered grid for accuracy reasons. The application of the ADI procedure for the parabolic velocity Poisson equations along with the continuity equation results in diagonally dominant tri‐diagonal matrix equations. Thus the explicit method for the vorticity equations and the tri‐diagonal matrix algorithm for the Poisson equations combine to give a simplified numerical scheme for solving three‐dimensional problems, which otherwise requires enormous computational effort. For three‐dimensional‐driven cavity flow predictions, the present method is found to be efficient and accurate for the Reynolds number range 100?Re?2000. Copyright © 2004 John Wiley & Sons, Ltd.     

A high‐order accurate method for two‐dimensional incompressible viscous flows

A high‐order accurate solution method for complex geometries is developed for two‐dimensional flows using the stream function–vorticity formulation. High‐order accurate spectrally optimized compact schemes along with appropriate boundary schemes are used for spatial discretization while a two‐level backward Euler implicit scheme is used for the time integration. The linear system of equations for stream function and vorticity are solved by an inner iteration while contravariant velocities constitute outer iterations. The effect of curvilinear grids on the solution accuracy is studied. The method is used to compute Cartesian and inclined driven cavity, flow in a triangular cavity and viscous flow in constricted channel. Benchmark‐like accuracy is obtained in all the problems with fewer grid points compared to reported studies. Copyright © 2006 John Wiley & Sons, Ltd.     

Combined compact difference method for solving the incompressible Navier–Stokes equations

This paper presents a numerical method for solving the two‐dimensional unsteady incompressible Navier–Stokes equations in a vorticity–velocity formulation. The method is applicable for simulating the nonlinear wave interaction in a two‐dimensional boundary layer flow. It is based on combined compact difference schemes of up to 12th order for discretization of the spatial derivatives on equidistant grids and a fourth‐order five‐ to six‐alternating‐stage Runge–Kutta method for temporal integration. The spatial and temporal schemes are optimized together for the first derivative in a downstream direction to achieve a better spectral resolution. In this method, the dispersion and dissipation errors have been minimized to simulate physical waves accurately. At the same time, the schemes can efficiently suppress numerical grid‐mesh oscillations. The results of test calculations on coarse grids are in good agreement with the linear stability theory and comparable with other works. The accuracy and the efficiency of the current code indicate its potential to be extended to three‐dimensional cases in which full boundary layer transition happens. Copyright © 2011 John Wiley & Sons, Ltd.     

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.     

Two‐dimensional compact finite difference immersed boundary method

We present a compact finite differences method for the calculation of two‐dimensional viscous flows in biological fluid dynamics applications. This is achieved by using body‐forces that allow for the imposition of boundary conditions in an immersed moving boundary that does not coincide with the computational grid. The unsteady, incompressible Navier–Stokes equations are solved in a Cartesian staggered grid with fourth‐order Runge–Kutta temporal discretization and fourth‐order compact schemes for spatial discretization, used to achieve highly accurate calculations. Special attention is given to the interpolation schemes on the boundary of the immersed body. The accuracy of the immersed boundary solver is verified through grid convergence studies. Validation of the method is done by comparison with reference experimental results. In order to demonstrate the application of the method, 2D small insect hovering flight is calculated and compared with available experimental and computational results. Copyright © 2009 John Wiley & Sons, Ltd.     

Gas‐assisted fluid displacement in a circular tube and a rectangular channel

In this paper the amount of liquid left inside of a circular tube and a rectangular channel when displaced by another immiscible fluid are determined by solving the full creeping‐motion equations. The exact continuity of stress on the free surface is employed with a finite difference method. In order to solve the equations, the steady‐state shape of the interface is guessed and the normal stress boundary condition is dropped. The equations based on a stream function‐vorticity formulation are solved with the aid of elliptic grid generation. The computed results are compared with the experimental results of Taylor (J. Fluid Mech. 1961; 10: 161), the theoretical results of Reinelt and Saffman (SIAM J. Sci. Stat. Comput. 1985; 6: 542) and our experimental data. The computed results are in close agreement with our experimental data and those of previous investigators. Copyright © 2002 John Wiley & Sons, Ltd.     

High‐order compact numerical schemes for non‐hydrostatic free surface flows

The development of a numerical scheme for non‐hydrostatic free surface flows is described with the objective of improving the resolution characteristics of existing solution methods. The model uses a high‐order compact finite difference method for spatial discretization on a collocated grid and the standard, explicit, single step, four‐stage, fourth‐order Runge–Kutta method for temporal discretization. The Cartesian coordinate system was used. The model requires the solution of two Poisson equations at each time‐step and tridiagonal matrices for each derivative at each of the four stages in a time‐step. Third‐ and fourth‐order accurate boundaries for the flow variables have been developed including the top non‐hydrostatic pressure boundary. The results demonstrate that numerical dissipation which has been a problem with many similar models that are second‐order accurate is practically eliminated. A high accuracy is obtained for the flow variables including the non‐hydrostatic pressure. The accuracy of the model has been tested in numerical experiments. In all cases where analytical solutions are available, both phase errors and amplitude errors are very small. Copyright © 2006 John Wiley & Sons, Ltd.     

Numerical simulation and experimental visualization of the influence of the deformation frequency of a radially deforming circular cylinder impulsively started on cylinder wake

A periodic superimposed motion may notably influence the flow structure and the development of the convective heat transfer relative to non‐deformable case. In particular, a radial deformation of a circular cylinder, may cause a possible synchronization with the cylinder wake, which is itself periodic when Vortex Street takes place. This synchronization phenomenon, often called ‘lock‐in’, may cause undesirable effects, but may also constitute a way of controlling the wake development. Body deformability may be used as wake control device that would favourably affect the interplay of primary and secondary vorticities, thus reducing the drag coefficient. These numerical and experimental studies are done herein for a Reynolds number equal to 23500. The problem is resolved by using the Navier–Stokes equations in the vorticity‐stream function form. The vorticity transport equation is solved by a second‐order finite difference method in both directions of the domains. The Poisson equation for the stream‐function is solved by a SOR method. The advance in time is achieved by a second‐order Adams–Bashforth scheme. The effect of turbulence is represented by eddy viscosity ν_t, which is determined by a sub‐grid‐scale model. In the present study, we use a Smagorinsky model. Copyright © 2003 John Wiley & Sons, Ltd.     

A parallel adaptive numerical method with generalized curvilinear coordinate transformation for compressible Navier–Stokes equations

A fourth‐order finite‐volume method for solving the Navier–Stokes equations on a mapped grid with adaptive mesh refinement is proposed, implemented, and demonstrated for the prediction of unsteady compressible viscous flows. The method employs fourth‐order quadrature rules for evaluating face‐averaged fluxes. Our approach is freestream preserving, guaranteed by the way of computing the averages of the metric terms on the faces of cells. The standard Runge–Kutta marching method is used for time discretization. Solutions of a smooth flow are obtained in order to verify that the method is formally fourth‐order accurate when applying the nonlinear viscous operators on mapped grids. Solutions of a shock tube problem are obtained to demonstrate the effectiveness of adaptive mesh refinement in resolving discontinuities. A Mach reflection problem is solved to demonstrate the mapped algorithm on a non‐rectangular physical domain. The simulation is compared against experimental results. Future work will consider mapped multiblock grids for practical engineering geometries. Copyright © 2016 John Wiley & Sons, Ltd.     

Vortex particle‐in‐cell method for three‐dimensional viscous unbounded flow computations

A new vortex particle‐in‐cell (PIC) method is developed for the computation of three‐dimensional unsteady, incompressible viscous flow in an unbounded domain. The method combines the advantages of the Lagrangian particle methods for convection and the use of an Eulerian grid to compute the diffusion and vortex stretching. The velocity boundary conditions used in the method are of Dirichlet‐type, and can be calculated using the vorticity field on the grid by the Biot–Savart equation. The present results for the propagation speed of the single vortex ring are in good agreement with the Saffman's model. The applications of the method to the head‐on and head‐off collisions of the two vortex rings show good agreement with the experimental and numerical literature. Copyright © 2000 John Wiley & Sons, Ltd.     

Least‐squares meshfree method for incompressible Navier–Stokes problems

A least‐squares meshfree method based on the first‐order velocity–pressure–vorticity formulation for two‐dimensional incompressible Navier–Stokes problem is presented. The convective term is linearized by successive substitution or Newton's method. The discretization of all governing equations is implemented by the least‐squares method. Equal‐order moving least‐squares approximation is employed with Gauss quadrature in the background cells. The boundary conditions are enforced by the penalty method. The matrix‐free element‐by‐element Jacobi preconditioned conjugate method is applied to solve the discretized linear systems. Cavity flow for steady Navier–Stokes problem and the flow over a square obstacle for time‐dependent Navier–Stokes problem are investigated for the presented least‐squares meshfree method. The effects of inaccurate integration on the accuracy of the solution are investigated. Copyright © 2004 John Wiley & Sons, Ltd.     

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.