Adaptive hybrid (prismatic–tetrahedral) grids for incompressible flows

Hybrid grids consisting of prisms and tetrahedra are employed for the solution of the 3-D Navier–Stokes equations of incompressible flow. A pressure correction scheme is employed with a finite volume–finite element spatial discretization. The traditional staggered grid formulation has been substituted with a collocated mesh approach which uses fourth-order artificial dissipation. The hybrid grid is refined adaptively in local regions of appreciable flow variations. The scheme operations are performed on an edge-wise basis which unifies treatment of both types of grid elements. The adaptive method is employed for incompressible flows in both single and multiply-connected domains. © 1998 John Wiley & Sons, Ltd.     

An overlapping control volume method for the Navier–Stokes equations on non‐staggered grids

An algorithm, based on the overlapping control volume (OCV) method, for the solution of the steady and unsteady two‐dimensional incompressible Navier–Stokes equations in complex geometry is presented. The primitive variable formulation is solved on a non‐staggered grid arrangement. The problem of pressure–velocity decoupling is circumvented by using momentum interpolation. The accuracy and effectiveness of the method is established by solving five steady state and one unsteady test problems. The numerical solutions obtained using the technique are in good agreement with the analytical and benchmark solutions available in the literature. On uniform grids, the method gives second‐order accuracy for both diffusion‐ and convection‐dominated flows. There is little loss of accuracy on grids that are moderately non‐orthogonal. Copyright © 1999 John Wiley & Sons, Ltd.     

A simple and efficient error analysis for multi‐step solution of the Navier–Stokes equations

A simple error analysis is used within the context of segregated finite element solution scheme to solve incompressible fluid flow. An error indicator is defined based on the difference between a numerical solution on an original mesh and an approximated solution on a related mesh. This error indicator is based on satisfying the steady‐state momentum equations. The advantages of this error indicator are, simplicity of implementation (post‐processing step), ability to show regions of high and/or low error, and as the indicator approaches zero the solution approaches convergence. Two examples are chosen for solution; first, the lid‐driven cavity problem, followed by the solution of flow over a backward facing step. The solutions are compared to previously published data for validation purposes. It is shown that this rather simple error estimate, when used as a re‐meshing guide, can be very effective in obtaining accurate numerical solutions. Copyright © 2002 John Wiley & Sons, Ltd.     

Stretched Cartesian grids for solution of the incompressible Navier–Stokes equations

Two Cartesian grid stretching functions are investigated for solving the unsteady incompressible Navier–Stokes equations using the pressure–velocity formulation. The first function is developed for the Fourier method and is a generalization of earlier work. This function concentrates more points at the centre of the computational box while allowing the box to remain finite. The second stretching function is for the second‐order central finite difference scheme, which uses a staggered grid in the computational domain. This function is derived to allow a direct discretization of the Laplacian operator in the pressure equation while preserving the consistent behaviour exhibited by the uniform grid scheme. Both functions are analysed for their effects on the matrix of the discretized pressure equation. It is shown that while the second function does not spoil the matrix diagonal dominance, the first one can. Limits to stretching of the first method are derived for the cases of mappings in one and two directions. A limit is also derived for the second function in order to prevent a strong distortion of a sine wave. The performances of the two types of stretching are examined in simulations of periodic co‐flowing jets and a time developing boundary layer. Copyright © 2000 John Wiley & Sons, Ltd.     

On numerical errors and turbulence modeling in tip vortex flow prediction

The accuracy of tip vortex flow prediction in the near‐field region is investigated numerically by attempting to quantify the shortcomings of the turbulence models and the flow solver. In particular, some turbulence models can produce a ‘numerical diffusion’ that artificially smears the vortex core. Low‐order finite differencing techniques of the convective and pressure terms of the Navier–Stokes equations and inadequate grid density and distribution can also produce the same adverse effect. The flow over a wing and the near‐wake with the wind tunnel walls included was simulated using 2.5 million grid points. Two subset problems, one using a steady, three‐dimensional analytical vortex, and the other, a vortex obtained from experiment and propagated downstream, were also devised in order to make the study of vortex preservation more tractable. The method of artificial compressibility is used to solve the steady, three‐dimensional, incompressible Navier–Stokes equations. Two one‐equation turbulence models (Baldwin–Barth and Spalart–Allmaras turbulence models), have been used with the production term modified to account for the stabilizing effect of the nearly solid body rotation in the vortex core. Finally, a comparison between the computed results and experiment is presented. Published in 1999 by John Wiley & Sons, Ltd.     

Colocated schemes for the incompressible Navier–Stokes equations on non‐smooth grids for two‐dimensional problems

The accuracy of colocated finite volume schemes for the incompressible Navier–Stokes equations on non‐smooth curvilinear grids is investigated. A frequently used scheme is found to be quite inaccurate on non‐smooth grids. In an attempt to improve the accuracy on such grids, three other schemes are described and tested. Two of these are found to give satisfactory results. Copyright © 2000 John Wiley & Sons, Ltd.     

SINGLE-LAYER FINITE ELEMENT APPROXIMATION OF 3D NAVIER–STOKES EQUATIONS WITH FREE BOUNDARIES FOR SHALLOW WATER FLOWS

In order to simulate flows in the shallow water limit, the full incompressible Navier–Stokes equations with free boundaries are solved using a single layer of finite elements. This implies a polynomial approximation of the velocity profile in the vertical direction, which in turn distorts the wave speed. This fact is verified by numerical results: the wave speed depends on the vertical discretization. When at least two layers of finite elements are used, the boundary layer at the bottom can be simulated and the correct solution for the shallow water limit is recovered. Then this algorithm is applied to the prediction of Tsunami event.     

On the implicit time integration of semi‐discrete viscous fluxes on unstructured dynamic meshes

Numerical simulations of viscous flow problems with complex moving and/or deforming boundaries commonly require the solution of the corresponding fluid equations of motion on unstructured dynamic meshes. In this paper, a systematic investigation of the importance of the choice of the mesh configuration for evaluating the viscous fluxes is performed when the semi‐discrete Navier–Stokes equations are time‐integrated using the popular second‐order implicit backward difference algorithm. The findings are illustrated with the simulation of a laminar viscous flow problem around an oscillating airfoil. Copyright © 1999 John Wiley & Sons, Ltd.     

Implicit weighted essentially non‐oscillatory schemes for the incompressible Navier–Stokes equations

A class of lower–upper/approximate factorization (LUAF) implicit weighted essentially non‐oscillatory (ENO; WENO) schemes for solving the two‐dimensional incompressible Navier–Stokes equations in a generalized co‐ordinate system is presented. The algorithm is based on the artificial compressibility formulation, and symmetric Gauss–Seidel relaxation is used for computing steady state solutions while symmetric successive overrelaxation is used for treating time‐dependent flows. WENO spatial operators are employed for inviscid fluxes and central differencing for viscous fluxes. Internal and external viscous flow test problems are presented to verify the numerical schemes. The use of a WENO spatial operator not only enhances the accuracy of solutions but also improves the convergence rate for the steady state computation as compared with using the ENO counterpart. It is found that the present solutions compare well with exact solutions, experimental data and other numerical results. Copyright © 1999 John Wiley & Sons, Ltd.     

A method for obtaining benchmark Navier–Stokes solutions

A refinement to an established method for obtaining benchmark Navier–Stokes solutions is presented. Pressure and body forces are derived explicitly such that the momentum equations are satisfied. The problem is reduced to determining a streamfunction in separation of variables form that describes a desired flow pattern. Examples based upon the well‐known shear flow cavity are presented. Copyright © 1999 John Wiley & Sons, Ltd.     

Comparison of pressure correction smoothers for multigrid solution of incompressible flow

We compare the performance of different pressure correction algorithms used as basic solvers in a multigrid method for the solution of the incompressible Navier–Stokes equations on non-staggered grids. Numerical tests were performed on several cases of lid-driven cavity flow using four different pressure correction schemes, including the traditional SIMPLE and SIMPLEC methods as well as novel variants, and varying combinations of underrelaxation parameters. The results show that three of the four algorithms tested are robust smoothers for the multigrid solver and that one of the new methods converges fastest in most of the tests. © 1997 John Wiley & Sons, Ltd.     

A COMPARATIVE STUDY OF TIME-STEPPING TECHNIQUES FOR THE INCOMPRESSIBLE NAVIER-STOKES EQUATIONS: FROM FULLY IMPLICIT NON-LINEAR SCHEMES TO SEMI-IMPLICIT PROJECTION METHODS

We present a numerical comparison of some time-stepping schemes for the discretization and solution of the non-stationary incompressible Navier– Stokes equations. The spatial discretization is by non-conforming quadrilateral finite elements which satisfy the LBB condition. The major focus is on the differences in accuracy and efficiency between the backward Euler, Crank–Nicolson and fractional-step Θ schemes used in discretizing the momentum equations. Further, the differences between fully coupled solvers and operator-splitting techniques (projection methods) and the influence of the treatment of the nonlinear advection term are considered. The combination of both discrete projection schemesand non-conforming finite elementsallows the comparison of schemes which are representative for many methods used in practice. On Cartesian grids this approach encompasses some well-known staggered grid finite difference discretizations too. The results which are obtained for several typical flow problems are thought to be representative and should be helpful for a fair rating of solution schemes, particularly in long-time simulations.