Stable Robin solid wall boundary conditions for the Navier–Stokes equations

In this paper we prove stability of Robin solid wall boundary conditions for the compressible Navier–Stokes equations. Applications include the no-slip boundary conditions with prescribed temperature or temperature gradient and the first order slip-flow boundary conditions. The formulation is uniform and the transitions between different boundary conditions are done by a change of parameters. We give different sharp energy estimates depending on the choice of parameters.     

Absorbing boundary conditions for the Euler and Navier–Stokes equations with the spectral difference method

Two absorbing boundary conditions, the absorbing sponge zone and the perfectly matched layer, are developed and implemented for the spectral difference method discretizing the Euler and Navier–Stokes equations on unstructured grids. The performance of both boundary conditions is evaluated and compared with the characteristic boundary condition for a variety of benchmark problems including vortex and acoustic wave propagations. The applications of the perfectly matched layer technique in the numerical simulations of unsteady problems with complex geometries are also presented to demonstrate its capability.     

A Fourier spectral method for the Navier–Stokes equations with volume penalization for moving solid obstacles

An efficient numerical scheme to compute flows past rigid solid bodies moving through viscous incompressible fluid is presented. Solid obstacles of arbitrary shape are taken into account using the volume penalization method to impose no-slip boundary condition. The 2D Navier–Stokes equations, written in the vorticity-streamfunction formulation, are discretized using a Fourier pseudo-spectral scheme. Four different time discretization schemes of the penalization term are proposed and compared. The originality of the present work lies in the implementation of time-dependent penalization, which makes the above method capable of solving problems where the obstacle follows an arbitrary motion. Fluid–solid coupling for freely falling bodies is also implemented. The numerical method is validated for different test cases: the flow past a cylinder, Couette flow between rotating cylinders, sedimentation of a cylinder and a falling leaf with elliptical shape.     

A fourth-order adaptive mesh refinement algorithm for the multicomponent,reacting compressible Navier–Stokes equations

In this paper, we present a fourth-order in space and time block-structured adaptive mesh refinement algorithm for the compressible multicomponent reacting Navier–Stokes equations. The algorithm uses a finite-volume approach that incorporates a fourth-order discretisation of the convective terms. The time-stepping algorithm is based on a multi-level spectral deferred corrections method that enables explicit treatment of advection and diffusion coupled with an implicit treatment of reactions. The temporal scheme is embedded in a block-structured adaptive mesh refinement algorithm that includes subcycling in time with spectral deferred correction sweeps applied on levels. Here we present the details of the multi-level scheme paying particular attention to the treatment of coarse–fine boundaries required to maintain fourth-order accuracy in time. We then demonstrate the convergence properties of the algorithm on several test cases including both non-reacting and reacting flows. Finally we present simulations of a vitiated dimethyl ether jet in 2D and a turbulent hydrogen jet in 3D, both with detailed kinetics and transport.     

Third order residual distribution schemes for the Navier–Stokes equations

We construct a third order multidimensional upwind residual distribution scheme for the system of the Navier–Stokes equations. The underlying approximation is obtained using standard P² Lagrange finite elements. To discretise the inviscid component of the equations, each element is divided in sub-elements over which we compute a high order residual defined as the integral of the inviscid fluxes on the boundary of the sub-element. The residuals are distributed to the nodes of each sub-element in a multi-dimensional upwind way. To obtain a discretisation of the viscous terms consistent with this multi-dimensional upwind approach, we make use of a Petrov–Galerkin analogy. The analogy allows to find a family of test functions which can be used to obtain a weak approximation of the viscous terms. The performance of this high-order method is tested on flows with high and low Reynolds number.     

Output-based space–time mesh adaptation for the compressible Navier–Stokes equations

This paper presents an output-based adaptive algorithm for unsteady simulations of convection-dominated flows. A space–time discontinuous Galerkin discretization is used in which the spatial meshes remain static in both position and resolution, and in which all elements advance by the same time step. Error estimates are computed using an adjoint-weighted residual, where the discrete adjoint is computed on a finer space obtained by order enrichment of the primal space. An iterative method based on an approximate factorization is used to solve both the forward and adjoint problems. The output error estimate drives a fixed-growth adaptive strategy that employs hanging-node refinement in the spatial domain and slab bisection in the temporal domain. Detection of space–time anisotropy in the localization of the output error is found to be important for efficiency of the adaptive algorithm, and two anisotropy measures are presented: one based on inter-element solution jumps, and one based on projection of the adjoint. Adaptive results are shown for several two-dimensional convection-dominated flows, including the compressible Navier–Stokes equations. For sufficiently-low accuracy levels, output-based adaptation is shown to be advantageous in terms of degrees of freedom when compared to uniform refinement and to adaptive indicators based on approximation error and the unweighted residual. Time integral quantities are used for the outputs of interest, but entire time histories of the integrands are also compared and found to converge rapidly under the proposed scheme. In addition, the final output-adapted space–time meshes are shown to be relatively insensitive to the starting mesh.     

An Asymptotic-Preserving all-speed scheme for the Euler and Navier–Stokes equations

We present an Asymptotic-Preserving ‘all-speed’ scheme for the simulation of compressible flows valid at all Mach-numbers ranging from very small to order unity. The scheme is based on a semi-implicit discretization which treats the acoustic part implicitly and the convective and diffusive parts explicitly. This discretization, which is the key to the Asymptotic-Preserving property, provides a consistent approximation of both the hyperbolic compressible regime and the elliptic incompressible regime. The divergence-free condition on the velocity in the incompressible regime is respected, and an the pressure is computed via an elliptic equation resulting from a suitable combination of the momentum and energy equations. The implicit treatment of the acoustic part allows the time-step to be independent of the Mach number. The scheme is conservative and applies to steady or unsteady flows and to general equations of state. One and two-dimensional numerical results provide a validation of the Asymptotic-Preserving ‘all-speed’ properties.     

A Relaxed Dimensional Factorization preconditioner for the incompressible Navier–Stokes equations

In this paper we introduce a Relaxed Dimensional Factorization (RDF) preconditioner for saddle point problems. Properties of the preconditioned matrix are analyzed and compared with those of the closely related Dimensional Splitting (DS) preconditioner recently introduced by Benzi and Guo [7]. Numerical results for a variety of finite element discretizations of both steady and unsteady incompressible flow problems indicate very good behavior of the RDF preconditioner with respect to both mesh size and viscosity.     

Preconditioning methods for discontinuous Galerkin solutions of the Navier–Stokes equations

A Newton–Krylov method is developed for the solution of the steady compressible Navier–Stokes equations using a discontinuous Galerkin (DG) discretization on unstructured meshes. Steady-state solutions are obtained using a Newton–Krylov approach where the linear system at each iteration is solved using a restarted GMRES algorithm. Several different preconditioners are examined to achieve fast convergence of the GMRES algorithm. An element Line-Jacobi preconditioner is presented which solves a block-tridiagonal system along lines of maximum coupling in the flow. An incomplete block-LU factorization (Block-ILU(0)) is also presented as a preconditioner, where the factorization is performed using a reordering of elements based upon the lines of maximum coupling. This reordering is shown to be superior to standard reordering techniques (Nested Dissection, One-way Dissection, Quotient Minimum Degree, Reverse Cuthill–Mckee) especially for viscous test cases. The Block-ILU(0) factorization is performed in-place and an algorithm is presented for the application of the linearization which reduces both the memory and CPU time over the traditional dual matrix storage format. Additionally, a linear p-multigrid preconditioner is also considered, where Block-Jacobi, Line-Jacobi and Block-ILU(0) are used as smoothers. The linear multigrid preconditioner is shown to significantly improve convergence in term of number of iterations and CPU time compared to a single-level Block-Jacobi or Line-Jacobi preconditioner. Similarly the linear multigrid preconditioner with Block-ILU smoothing is shown to reduce the number of linear iterations to achieve convergence over a single-level Block-ILU(0) preconditioner, though no appreciable improvement in CPU time is shown.     

On physics-based preconditioning of the Navier–Stokes equations

We develop a fully implicit scheme for the Navier–Stokes equations, in conservative form, for low to intermediate Mach number flows. Simulations in this range of flow regime produce stiff wave systems in which slow dynamical (advective) modes coexist with fast acoustic modes. Viscous and thermal diffusion effects in refined boundary layers can also produce stiffness. Implicit schemes allow one to step over the fast wave phenomena (or unresolved viscous time scales), while resolving advective time scales. In this study we employ the Jacobian-free Newton–Krylov (JFNK) method and develop a new physics-based preconditioner. To aid in overcoming numerical stiffness caused by the disparity between acoustic and advective modes, the governing equations are transformed into the primitive-variable form in a preconditioning step. The physics-based preconditioning incorporates traditional semi-implicit and physics-based splitting approaches without a loss of consistency between the original and preconditioned systems. The resulting algorithm is capable of solving low-speed natural circulation problems (M∼10^-4)$(M\sim {10}^{-4})$ with significant heat flux as well as intermediate speed (M∼1)$(M\sim 1)$ flows efficiently by following dynamical (advective) time scales of the problem.     

Beyond the Navier–Stokes equations: Burnett hydrodynamics

This work is mainly concerned with the extension of hydrodynamics beyond the Navier–Stokes equations, a regime known as Burnett hydrodynamics. The derivation of the Burnett equations is considered from several theoretical approaches. In particular we discuss the Chapman–Enskog, Grad’s method, and Truesdell’s approach for solving the Boltzmann equation. Also, their derivation using the macroscopic approach given by extended thermodynamics is mentioned. The problems and successes of these equations are discussed and some alternatives proposed to improve them are mentioned. Comparisons of the predictions coming from the Burnett equations with experiments and/or simulations are given in order to have the necessary elements to give a critical assessment of their validity and usefulness.     

Absorbing boundary conditions for nonlinear Euler and Navier–Stokes equations based on the perfectly matched layer technique

Absorbing boundary conditions for the nonlinear Euler and Navier–Stokes equations in three space dimensions are presented based on the perfectly matched layer (PML) technique. The derivation of equations follows a three-step method recently developed for the PML of linearized Euler equations. To increase the efficiency of the PML, a pseudo mean flow is introduced in the formulation of absorption equations. The proposed PML equations will absorb exponentially the difference between the nonlinear fluctuation and the prescribed pseudo mean flow. With the nonlinearity in flux vectors, the proposed nonlinear absorbing equations are not formally perfectly matched to the governing equations as their linear counter-parts are. However, numerical examples show satisfactory results. Furthermore, the nonlinear PML reduces automatically to the linear PML upon linearization about the pseudo mean flow. The validity and efficiency of proposed equations as absorbing boundary conditions for nonlinear Euler and Navier–Stokes equations are demonstrated by numerical examples.     

A finite element method for the numerical solution of the coupled Cahn–Hilliard and Navier–Stokes system for moving contact line problems

In this paper, a semi-implicit finite element method is presented for the coupled Cahn–Hilliard and Navier–Stokes equations with the generalized Navier boundary condition for the moving contact line problems. In our method, the system is solved in a decoupled way. For the Cahn–Hilliard equations, a convex splitting scheme is used along with a P1-P1 finite element discretization. The scheme is unconditionally stable. A linearized semi-implicit P2-P0 mixed finite element method is employed to solve the Navier–Stokes equations. With our method, the generalized Navier boundary condition is extended to handle the moving contact line problems with complex boundary in a very natural way. The efficiency and capacity of the present method are well demonstrated with several numerical examples.     

Adjoint algorithms for the Navier–Stokes equations in the low Mach number limit

This paper describes a derivation of the adjoint low Mach number equations and their implementation and validation within a global mode solver. The advantage of using the low Mach number equations and their adjoints is that they are appropriate for flows with variable density, such as flames, but do not require resolution of acoustic waves. Two versions of the adjoint are implemented and assessed: a discrete-adjoint and a continuous-adjoint. The most unstable global mode calculated with the discrete-adjoint has exactly the same eigenvalue as the corresponding direct global mode but contains numerical artifacts near the inlet. The most unstable global mode calculated with the continuous-adjoint has no numerical artifacts but a slightly different eigenvalue. The eigenvalues converge, however, as the timestep reduces. Apart from the numerical artifacts, the mode shapes are very similar, which supports the expectation that they are otherwise equivalent. The continuous-adjoint requires less resolution and usually converges more quickly than the discrete-adjoint but is more challenging to implement. Finally, the direct and adjoint global modes are combined in order to calculate the wavemaker region of a low density jet.     

Implicit LU-SGS algorithm for high-order methods on unstructured grid with p-multigrid strategy for solving the steady Navier–Stokes equations

The fluid dynamic equations are discretized by a high-order spectral volume (SV) method on unstructured tetrahedral grids. We solve the steady state equations by advancing in time using a backward Euler (BE) scheme. To avoid the inversion of a large matrix we approximate BE by an implicit lower–upper symmetric Gauss–Seidel (LU-SGS) algorithm. The implicit method addresses the stiffness in the discrete Navier–Stokes equations associated with stretched meshes. The LU-SGS algorithm is then used as a smoother for a p-multigrid approach. A Von Neumann stability analysis is applied to the two-dimensional linear advection equation to determine its damping properties. The implicit LU-SGS scheme is used to solve the two-dimensional (2D) compressible laminar Navier–Stokes equations. We compute the solution of a laminar external flow over a cylinder and around an airfoil at low Mach number. We compare the convergence rates with explicit Runge–Kutta (E-RK) schemes employed as a smoother. The effects of the cell aspect ratio and the low Mach number on the convergence are investigated. With the p-multigrid method and the implicit smoother the computational time can be reduced by a factor of up to 5–10 compared with a well tuned E-RK scheme.     

h-Multigrid for space-time discontinuous Galerkin discretizations of the compressible Navier–Stokes equations

Being implicit in time, the space-time discontinuous Galerkin discretization of the compressible Navier–Stokes equations requires the solution of a non-linear system of algebraic equations at each time-step. The overall performance, therefore, highly depends on the efficiency of the solver. In this article, we solve the system of algebraic equations with a h-multigrid method using explicit Runge–Kutta relaxation. Two-level Fourier analysis of this method for the scalar advection–diffusion equation shows convergence factors between 0.5 and 0.75. This motivates its application to the 3D compressible Navier–Stokes equations where numerical experiments show that the computational effort is significantly reduced, up to a factor 10 w.r.t. single-grid iterations.     

An efficient method for the incompressible Navier–Stokes equations on irregular domains with no-slip boundary conditions,high order up to the boundary

Common efficient schemes for the incompressible Navier–Stokes equations, such as projection or fractional step methods, have limited temporal accuracy as a result of matrix splitting errors, or introduce errors near the domain boundaries (which destroy uniform convergence to the solution). In this paper we recast the incompressible (constant density) Navier–Stokes equations (with the velocity prescribed at the boundary) as an equivalent system, for the primary variables velocity and pressure. equation for the pressure. The key difference from the usual approaches occurs at the boundaries, where we use boundary conditions that unequivocally allow the pressure to be recovered from knowledge of the velocity at any fixed time. This avoids the common difficulty of an, apparently, over-determined Poisson problem. Since in this alternative formulation the pressure can be accurately and efficiently recovered from the velocity, the recast equations are ideal for numerical marching methods. The new system can be discretized using a variety of methods, including semi-implicit treatments of viscosity, and in principle to any desired order of accuracy. In this work we illustrate the approach with a 2-D second order finite difference scheme on a Cartesian grid, and devise an algorithm to solve the equations on domains with curved (non-conforming) boundaries, including a case with a non-trivial topology (a circular obstruction inside the domain). This algorithm achieves second order accuracy in the L^∞ norm, for both the velocity and the pressure. The scheme has a natural extension to 3-D.     

Weak and strong wall boundary procedures and convergence to steady-state of the Navier–Stokes equations

We study the influence of different implementations of no-slip solid wall boundary conditions on the convergence to steady-state of the Navier–Stokes equations. The various approaches are investigated using the energy method and an eigenvalue analysis. It is shown that the weak implementation is superior and enhances the convergence to steady-state for coarse meshes. It is also demonstrated that all the stable approaches produce the same convergence rate as the mesh size goes to zero. The numerical results obtained by using a fully nonlinear finite volume solver support the theoretical findings from the linear analysis.     

Analogue of the Cole-Hopf transform for the incompressible Navier–Stokes equations and its application

We consider the Navier–Stokes equations written in the stream function in two dimensions and vector potentials in three dimensions, which are critical dependent variables. On this basis, we introduce an analogue of the Cole-Hopf transform, which exactly reduces the Navier–Stokes equations to the heat equations with a potential term (i.e. the nonlinear Schrödinger equation at imaginary times). The following results are obtained. (i) A regularity criterion immediately obtains as the boundedness of condition for the potential term when the equations are recast in a path-integral form by the Feynman-Kac formula. (ii) This in turn gives an additional characterisation of possible singularities for the Navier–Stokes equations. (iii) Some numerical results for the two-dimensional Navier–Stokes equations are presented to demonstrate how the potential term captures near-singular structures. Finally, we extend this formulation to higher dimensions, where the regularity issues are markedly open.     

The Elastoplast Discontinuous Galerkin (EDG) method for the Navier–Stokes equations

The present work details the Elastoplast (this name is a translation from the French “sparadrap”, a concept first applied by Yves Morchoisne for Spectral methods [1]) Discontinuous Galerkin (EDG) method to solve the compressible Navier–Stokes equations. This method was first presented in 2009 at the ICOSAHOM congress with some Cartesian grid applications. We focus here on unstructured grid applications for which the EDG method seems very attractive. As in the Recovery method presented by van Leer and Nomura in 2005 for diffusion, jumps across element boundaries are locally eliminated by recovering the solution on an overlapping cell. In the case of Recovery, this cell is the union of the two neighboring cells and the Galerkin basis is twice as large as the basis used for one element. In our proposed method the solution is rebuilt through an L² projection of the discontinuous interface solution on a small rectangular overlapping interface element, named Elastoplast, with an orthogonal basis of the same order as the one in the neighboring cells. Comparisons on 1D and 2D scalar diffusion problems in terms of accuracy and stability with other viscous DG schemes are first given. Then, 2D results on acoustic problems, vortex problems and boundary layer problems both on Cartesian or unstructured triangular grids illustrate stability, precision and versatility of this method.