Parallel computation of unsteady three-dimensional incompressible viscous flow using an unstructured multigrid method

The development and validation of a parallel unstructured tetrahedral non-nested multigrid (MG) method for simulation of unsteady 3D incompressible viscous flow is presented. The Navier-Stokes solver is based on the artificial compressibility method (ACM) and a higher-order characteristics-based finite-volume scheme on unstructured MG. Unsteady flow is calculated with an implicit dual time stepping scheme. The parallelization of the solver is achieved by a MG domain decomposition approach (MG-DD), using the Single Program Multiple Data (SPMD) programming paradigm. The Message-Passing Interface (MPI) Library is used for communication of data and loop arrays are decomposed using the OpenMP standard. The parallel codes using single grid and MG are used to simulate steady and unsteady incompressible viscous flows for a 3D lid-driven cavity flow for validation and performance evaluation purposes. The speedups and efficiencies obtained by both the parallel single grid and MG solvers are reasonably good for all test cases, using up to 32 processors on the SGI Origin 3400. The parallel results obtained agree well with those of serial solvers and with numerical solutions obtained by other researchers, as well as experimental measurements.     

An implicit high-order discontinuous Galerkin method for steady and unsteady incompressible flows

This paper presents the latest developments of a discontinuous Galerkin (DG) method for incompressible flows introduced in [Bassi F, Crivellini A, Di Pietro DA, Rebay S. An artificial compressibility flux for the discontinuous Galerkin solution of the incompressible Navier–Stokes equations. J Comput Phys 2006;218(2):794–815] for the steady Navier–Stokes equations and extended in [Bassi F, Crivellini A. A high-order discontinuous Galerkin method for natural convection problems. In: Wesseling P, Oñate E, Periaux J, editors. Electronic proceedings of the ECCOMAS CFD 2006 conference, Egmond aan Zee, The Netherlands, September 5–8; 2006. TU Delft] to the coupled Navier–Stokes and energy equations governing natural convection flows.

The method is fully implicit and applies to the governing equations in primitive variable form. Its distinguishing feature is the formulation of the inviscid interface flux, which is based on the solution of local Riemann problems associated with the artificial compressibility perturbation of the Euler equations. The tight coupling between pressure and velocity so introduced stabilizes the method and allows using equal-order approximation spaces for both pressure and velocity. Since, independently of the amount of artificial compressibility added, the interface flux reduces to the physical one for vanishing interface jumps, the resulting method is strongly consistent.

In this paper, we present a review of the method together with two recently developed issues: (i) the high-order DG discretization of the incompressible Euler equations; (ii) the high-order implicit time integration of unsteady flows. The accuracy and versatility of the method are demonstrated by a suite of computations of steady and unsteady, inviscid and viscous incompressible flows.     

A Legendre spectral element method for simulation of unsteady incompressible viscous free-surface flows

A new Legendre spectral element method is presented for the solution of viscous incompressible free-surface flows. It is based on the following extensions of the fixed-domain spectral element method: use of the full viscous stress tensor for natural imposition of traction (surface tension) boundary conditions; use of arbitrary-Lagrangian-Eulerian methods for accurate representation of moving boundaries; and use of semi-implicit time-stepping procedures to partially decouple the free-surface evolution and interior Navier-Stokes equations. For purposes of analysis and clarity of presentation, attention is focused on the stability of falling films. Analysis of the spectrum of the linear stability problem (Orr-Sommerfeld equation) associated with film flow reveals physical effects that limit the stability of semi-implicit schemes and suggests optimal formulas for temporal discretization of the spectral element equations. Detailed results are presented for the spectral element simulation of the film flow problem.     

Pseudospectral multi-domain method for incompressible viscous flow computation

We present a multi-domain pseudospectral method for the calculation of incompressible viscous flow. Governing equations are written in primitive variables formulation. Velocity components and pressure are discretized on the same grid of collocation points. A coupling algorithm for the Stokes problem is investigated and preliminary results are presented in the two-dimensional case.     

A computational scheme in generalized coordinates for viscous incompressible flows

This paper presents a computational scheme suitable for analyzing viscous incompressible flows in generalized curvilinear coordinate system. The scheme is based on finite volume algorithm with an overlapping staggered grid. The pseudo-diffusive terms arising from the coordinate transformation are treated as source terms. The system of nonlinear algebraic equations is solved by a semi-implicit procedure based upon line-relaxation and a generalization of Patankar's pressure correction algorithm. Examples of the application of the algorithm to flow in convergent channels, developing flow in a U-bend, and flow past backward facing step, are given. In addition, the case of flow past backward facing step is analyzed in detail, and the computed flowfields are found to be in close agreement with previous experimental and numerical results for expansion ratio (defined as the ratio of step height to channel height) of 0.5. The results are summarized in the form of a correlation relating the primary separation length, Reynolds number and expansion ratio.     

An iterative scheme for numerical solution of steady incompressible viscous flows

An iterative solution scheme is proposed for application to steady incompressible viscous flows in simple and complex geometries. The iterative scheme solves the vorticity-stream function form of the Navier-Stokes equations in generalized curvilinear coordinates. The flow system of equations are cast into a Newton's iterative form which are solved using the modified strongly implicit procedure. The solution scheme is benchmarked using two test cases, namely: a shear-driven steady laminar flow in a square cavity; and a simple laminar flow in a complex expanding channel. The iterative process to steady-state convergence in both test cases is highly stable and the convergence rate is without spurious oscillations. At convergence, the flow solutions are second-order accurate.     

A SVD-GFD scheme for computing 3D incompressible viscous fluid flows

A generalized finite difference (GFD) scheme for the simulation of three-dimensional (3D) incompressible viscous fluid flows in primitive variables is described in this paper. Numerical discretization is carried out on a hybrid Cartesian cum meshfree grid, with derivative approximation on non-Cartesian grids being carried out by a singular value decomposition (SVD) based GFD procedure. The Navier-Stokes equations are integrated by a time-splitting pressure correction scheme with second-order Crank-Nicolson and second-order discretization of time and spatial derivatives respectively. Axisymmetric and asymmetric 3D flows past a sphere with Reynolds numbers of up to 300 are simulated and compared with the results of Johnson and Patel [Johnson TA, Patel VC. Flow past a sphere up to a Reynolds number of 300. J Fluid Mech 1999;378:19-70] and others. Flows past toroidal rings are also simulated to illustrate the ability of the scheme to deal with more complex body geometry. The current method can also deal with flow past 3D bodies with sharp edges and corners, which is shown by a simple 3D case.     

Algorithmic and theoretical results on computation of incompressible viscous flows by finite element methods

Primitive variable as well as streamfunction-vorticity and pure streamfunction formulations are discussed. For the primitive variable case alternative choices of the viscous stress term are shown to produce natural boundary conditions which are well suited for matching to various far field conditions. For the other cases recent analytical results, including error estimates are described, and an optical algorithm for pressure recovery as well as treatment for multiply connected domains are given.     

Computational challenges of viscous incompressible flows

Over the past 30 years, numerical methods and simulation tools for incompressible flows have been advanced as a subset of the computational fluid dynamics (CFD) discipline. Although incompressible flows are encountered in many areas of engineering, simulation of compressible flow has been the major driver for developing computational algorithms and tools. This is probably due to the rather stringent requirements for predicting aerodynamic performance characteristics of flight vehicles, while flow devices involving low-speed or incompressible flow could be reasonably well designed without resorting to accurate numerical simulations. As flow devices are required to be more sophisticated and highly efficient, CFD tools become increasingly important in fluid engineering for incompressible and low-speed flow. This paper reviews some of the successes made possible by advances in computational technologies during the same period, and discusses some of the current challenges faced in computing incompressible flows.     

A finite element method for diffusion dominated unsteady viscous flows

A general conforming finite element scheme for computing viscous flows is presented which is of second-order accuracy in space and time. Viscous terms are treated implicitly and advection terms are treated explicitly in the time marching segment of the algorithm. A method for solving the algebraic equations at each time step is given. The method is demonstrated on two test problems, one of them being a plane vortex flow for which asymptotic methods are used to obtain suitable numerical boundary conditions at each time step.     

An implicit matrix-free Discontinuous Galerkin solver for viscous and turbulent aerodynamic simulations

This paper presents some recent advancements of the computational efficiency of a Discontinuous Galerkin (DG) solver for the Navier–Stokes (NS) and Reynolds Averaged Navier Stokes (RANS) equations. The implementation and the performance of a Newton–Krylov matrix-free (MF) method is presented and compared with the matrix based (MB) counterpart. Moreover two solution strategies, developed in order to increase the solver efficiency, are discussed and experimented. Numerical results of some test cases proposed within the EU ADIGMA (Adaptive Higher-Order Variational Methods for Aerodynamic Applications in Industry) project demonstrate the capabilities of the method.     

A class of high-resolution algorithms for incompressible flows

We present a class of a high-resolution Godunov-type algorithms for solving flow problems governed by the incompressible Navier-Stokes equations. The algorithms use high-resolution finite volume methods developed in LeVeque (SIAM J Numer Anal 1996;33:627-665) for the advective terms and finite difference methods for the diffusion and the Poisson pressure equation. The high-resolution algorithm advects the cell-centered velocities using the divergence-free cell-edge velocities. The resulting cell-centered velocity is then updated by the solution of the Poisson equation. The algorithms are proven to be robust for constant-density flows at high Reynolds numbers via an example of lid-driven cavity flow. With a slight modification for the projection operator in the constant-density solvers, the algorithms also solve incompressible flows with finite-amplitude density variation. The strength of such algorithms is illustrated through problems like Rayleigh-Taylor instability and the Boussinesq equations for Rayleigh-Bénard convection. Numerical studies of the convergence and order of accuracy for the velocity field are provided. While simulations for two-dimensional regular-geometry problems are presented in this study, in principle, extension of the algorithms to three dimensions with complex geometry is feasible.     

Applications of stencil-adaptive finite difference method to incompressible viscous flows with curved boundary

This paper investigates the applicability of the stencil-adaptive finite difference method for the simulation of two-dimensional unsteady incompressible viscous flows with curved boundary. The adaptive stencil refinement algorithm has been proven to be able to continuously adapt the stencil resolution according to the gradient of flow parameter of interest [Ding H, Shu C. A stencil adaptive algorithm for finite difference solution of incompressible viscous flows. J Comput Phys 2006;214:397-420], which facilitates the saving of the computational efforts. On the other hand, the capability of the domain-free discretization technique in dealing with the curved boundary provides a great flexibility for the finite difference scheme on the Cartesian grid. Here, we show that their combination makes it possible to simulate the unsteady incompressible flow with curved boundary on a dynamically changed grid. The methods are validated by simulating steady and unsteady incompressible viscous flows over a stationary circular cylinder.     

Computation of complex turbulent flow using matrix-free implicit dual time-stepping scheme and LRN turbulence model on unstructured grids

In this study, a matrix-free implicit dual time-stepping method has been developed. It is implemented, together with a low-Reynolds-number q-ω turbulence model, in a high-order upwind finite-volume solver on unstructured grids. Semi-implicit treatment of the source terms of the q and ω equations is also introduced to further stabilize the numerical solution. It has been found that these techniques provide strong stabilization in the computation of a supersonic flow with complex shock-boundary-layer interactions in a channel with a backward-facing step. The proposed method has a low-memory overhead, similar to an explicit scheme, while it shows good stability and computational efficiency as an implicit scheme. The method developed has been validated by comparing the computed results with the corresponding experimental measurements and other calculated results, which shows good agreement. Research is being done to extend the method to calculate unsteady turbulent flows.     

Optimal aerodynamic design of airfoils in unsteady viscous flows

A continuous adjoint formulation is used to determine optimal airfoil shapes in unsteady viscous flows at Re = 1 × 10⁴. The Reynolds number is based on the free-stream speed and the chord length of the airfoil. A finite element method based on streamline-upwind Petrov/Galerkin (SUPG) and pressure-stabilized Petrov/Galerkin (PSPG) stabilizations is used to solve both the flow and adjoint equations. The airfoil is parametrized via a Non-Uniform Rational B-Splines (NURBS) curve. Three different objective functions are used to obtain optimal shapes: maximize lift, minimize drag and minimize ratio of drag to lift. The objective functions are formulated on the basis of time-averaged aerodynamic coefficients. The three objective functions result in diverse airfoil geometries. The resulting airfoils are thin, with the largest thickness to chord ratio being only 5.4%. The shapes obtained are further investigated for their aerodynamic performance. Maximization of time-averaged lift leads to an airfoil that produces more than six times more lift compared to the NACA 0012 airfoil. The excess lift is a consequence of the large peak and extended region of high suction on the upper surface and high pressure on the lower surface. Minimization of drag results in an airfoil with a sharp leading edge. The flow remains attached for close to 70% of the chord length. Minimization of the ratio of drag to lift results in an airfoil with a shallow dimple on the upper surface. It leads to a fairly large value of the time-averaged ratio of lift to drag (～ 17.8). The high value is mostly achieved by a 447% increase in lift and 16% reduction in drag, compared to a NACA 0012 airfoil. Imposition of volume constraint, for the cases studied, is found to result in airfoils that have lower aerodynamic performance.     

Numerical simulations of viscous rotating flows using a new pressure-based method

This paper describes the development and application of a parabolic finite-volume procedure for the computation of viscous rotating flows having complex geometries. A parabolic numerical scheme has been developed to include the influence of the Coriolis force on pressure corrections. Two-equation turbulence models were employed to account for the effect of turbulence on the flows. Numerical tests for 2-D flows have been conducted, and the results show that computation with the new pressure-based method gives reasonably good agreement with the experiments.     

Design of airfoils in incompressible viscous flows by numerical optimization

A method is outlined for the design of airfoils in incompressible viscous flows by numerical optimization wherein a reduced number of design coordinates are used to define the airfoil shape. The optimization problem is formulated as a nongradient search in a finite constrained parameter space. The approach is to define the airfoil as a linear combination of basic shapes which may be analytically or numerically defined. The design problem is to determine the participation of each of these basic shapes in defining the optimum airfoil. The aerodynamic analysis program is specially developed to fit the requirements of the optimization program and is based on the vortex singularity method for inviscid flow analysis and the momentum integral method for boundary layer analysis. Four examples have been worked out to illustrate the proposed design method. In these, modifications to four different airfoil geometries are made to achieve either a minimum drag coefficient or a minimum pitching moment coefficient under prescribed constraints. The results show that significant drag or pitching moment reduction is possible through shape manipulation alone.     

Numerical modelling of interfacial soil erosion with viscous incompressible flows

The aim of this study is to develop a numerical model for simulating surface erosion occurring at a fluid/soil interface subject to a flow process. Balance equations with jump relations are used. A penalization procedure including a fictitious domain method is used to compute the Stokes flow around obstacles, in order to avoid body-fitted unstructured meshes and instead use fast and efficient finite volume approximations on Cartesian meshes. The evolution of the water/soil interface is described by using a level set function. The ability of the model to predict the interfacial erosion of soils is confirmed by several numerical simulations.     

Manifold method coupled velocity and pressure for Navier-Stokes equations and direct numerical solution of unsteady incompressible viscous flow

Numerical manifold method (NMM) application to direct numerical solution for unsteady incompressible viscous flow Navier-Stokes (N-S) equations was discussed in this paper, and numerical manifold schemes for N-S equations were derived based on Galerkin weighted residuals method as well. Mixed covers with linear polynomial function for velocity and constant function for pressure was employed in finite element cover system. The patch test demonstrated that mixed covers manifold elements meet the stability conditions and can be applied to solve N-S equations coupled velocity and pressure variables directly. The numerical schemes with mixed covers have also been proved to be unconditionally stable. As applications, mixed cover 4-node rectangular manifold element has been used to simulate the unsteady incompressible viscous flow in typical driven cavity and flow around a square cylinder in a horizontal channel. High accurate results obtained from much less calculational variables and very large time steps are in very good agreement with the compact finite difference solutions from very fine element meshes and very less time steps in references. Numerical tests illustrate that NMM is an effective and high order accurate numerical method for unsteady incompressible viscous flow N-S equations.     

Assessment of grid interface treatments for multi-block incompressible viscous flow computation

In the multi-block computation of the Navier-Stokes equations, the interface treatment is a key issue. In the present work, we investigate this issue in the context of a pressure-based method using a non-orthogonal grid. For the momentum equations, a straightforward bilinear interpolation seems satisfactory as the interface treatment; on the other hand, because the pressure field depends on the satisfaction of the mass continuity equation, a conservative interface treatment has been found necessary for the pressure-correction equation. Two alternative interface treatments for the pressure-correction equation, one employing the Neumann boundary condition in both grid blocks, based on explicit local, cell-by-cell mass flux conservation, and the other utilizing Neumann-Dirichlet boundary conditions, allowing the interface condition in one block to be derived by interpolating the pressure field from the adjacent block, are assessed in the present work. To evaluate these interface schemes, the laminar flow inside a lid-driven cavity flow, and the turbulent flow around cascades of multiple airfoils have been investigated. For the case tested, both interface treatments give comparable accuracy. The finding that more than one type of interface treatment can work well allows one to devise a flexible multi-block strategy for complex flow computations.