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.     

NUMERICAL SOLUTIONS OF OPTIMAL CONTROL FOR THERMALLY CONVECTIVE FLOWS

We study the numerical solution of optimal control problems associated with two-dimensional viscous incompressible thermally convective flows. Although the techniques apply to more general settings, the presentation is confined to the objectives of minimizing the vorticity in the steady state case and tracking the velocity field in the non-stationary case with boundary temperature controls. In the steady state case we develop a systematic way to use the Lagrange multiplier rules to derive an optimality system of equations from which an optimal solution can be computed; finite element methods are used to find approximate solutions for the optimality system of equations. In the time-dependent case a piecewise-in-time optimal control approach is proposed and the fully discrete approximation algorithm for solving the piecewise optimal control problem is defined. Numerical results are presented for both the steady state and time-dependent optimal control problems. © 1997 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.     

An experimental and numerical investigation of container filling with viscous liquids

This work is concerned with a study of container filling, with particular reference to the food industry. A computer code was developed and an experimental rig was built, the main purpose being to validate the software. The computational fluid dynamic (CFD) code, called GENSMAC, was specifically designed for relatively slow viscous flow and was capable of capturing multiple free surfaces. This paper focuses on the design of the experimental rig and how it functions. The visual output of the code is then compared with high‐speed photographic shots of glucose syrup being jetted into a tub for a selected number of flow regimes. Copyright © 1999 John Wiley & Sons, Ltd.     

A stabilized finite element method for the incompressible Navier–Stokes equations using a hierarchical basis

Stabilized finite element methods have been shown to yield robust, accurate numerical solutions to both the compressible and incompressible Navier–Stokes equations for laminar and turbulent flows. The present work focuses on the application of higher‐order, hierarchical basis functions to the incompressible Navier–Stokes equations using a stabilized finite element method. It is shown on a variety of problems that the most cost‐effective simulations (in terms of CPU time, memory, and disk storage) can be obtained using higher‐order basis functions when compared with the traditional linear basis. In addition, algorithms will be presented for the efficient implementation of these methods within the traditional finite element data structures. Copyright © 2001 John Wiley & Sons, Ltd.     

Analysis of super compact finite difference method and application to simulation of vortex–shock interaction

Turbulence and aeroacoustic noise high‐order accurate schemes are required, and preferred, for solving complex flow fields with multi‐scale structures. In this paper a super compact finite difference method (SCFDM) is presented, the accuracy is analysed and the method is compared with a sixth‐order traditional and compact finite difference approximation. The comparison shows that the sixth‐order accurate super compact method has higher resolving efficiency. The sixth‐order super compact method, with a three‐stage Runge–Kutta method for approximation of the compressible Navier–Stokes equations, is used to solve the complex flow structures induced by vortex–shock interactions. The basic nature of the near‐field sound generated by interaction is studied. Copyright © 2001 John Wiley & Sons, Ltd.     

PEGASE: A NAVIER–STOKES SOLVER FOR DIRECT NUMERICAL SIMULATION OF INCOMPRESSIBLE FLOWS

A hybrid conservative finite difference/finite element scheme is proposed for the solution of the unsteady incompressible Navier–Stokes equations. Using velocity–pressure variables on a non-staggeredgrid system, the solution is obtained with a projection method basedon the resolution of a pressure Poisson equation. The new proposed scheme is derived from the finite element spatial discretization using the Galerkin method with piecewise bilinear polynomial basis functions defined on quadrilateral elements. It is applied to the pressure gradient term and to the non-linear convection term as in the so-called group finite element method. It ensures strong coupling between spatial directions, inhibiting the development of oscillations during long-term computations, as demonstrated by the validation studies. Two- and three-dimensional unsteady separated flows with open boundaries have been simulated with the proposed method using Cartesian uniform mesh grids. Several examples of calculations on the backward-facing step configuration are reported and the results obtained are compared with those given by other methods. © 1997 by John Wiley & Sons, Ltd. Int. j. numer. methods fluids 24: 833–861, 1997.     

Numerical solution of the incompressible Navier–Stokes equations with an upwind compact difference scheme

A new finite difference method for the discretization of the incompressible Navier–Stokes equations is presented. The scheme is constructed on a staggered‐mesh grid system. The convection terms are discretized with a fifth‐order‐accurate upwind compact difference approximation, the viscous terms are discretized with a sixth‐order symmetrical compact difference approximation, the continuity equation and the pressure gradient in the momentum equations are discretized with a fourth‐order difference approximation on a cell‐centered mesh. Time advancement uses a three‐stage Runge–Kutta method. The Poisson equation for computing the pressure is solved with preconditioning. Accuracy analysis shows that the new method has high resolving efficiency. Validation of the method by computation of Taylor's vortex array is presented. Copyright © 1999 John Wiley & Sons, Ltd.     

Comment on ‘generalized potential flow theory and direct calculation of velocities applied to the numerical solution of the Navier–Stokes and the Boussinesq equations’

In a recent paper a generalized potential flow theory and its application to the solution of the Navier–Stokes equation are developed.¹ The purpose of this comment is to show that the analysis presented in that paper is in general not correct. We note that the theoretical development of Reference 1 is in fact an extension—although not cited—of some work first done by Hawthorne for steady inviscid flow.² Hawthorne's solution is correct, and his analysis, which we briefly describe, provides a useful introduction to this note.     

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.