Finite element and implicit Runge–Kutta implementation of an acoustics–convection upstream resolution algorithm for the time‐dependent two‐dimensional Euler equations

The second of a two‐paper series, this paper details a solver for the characteristics‐bias system from the acoustics–convection upstream resolution algorithm for the Euler and Navier–Stokes equations. An integral formulation leads to several surface integrals that allow effective enforcement of boundary conditions. Also presented is a new multi‐dimensional procedure to enforce a pressure boundary condition at a subsonic outlet, a procedure that remains accurate and stable. A classical finite element Galerkin discretization of the integral formulation on any prescribed grid directly yields an optimal discretely conservative upstream approximation for the Euler and Navier–Stokes equations, an approximation that remains multi‐dimensional independently of the orientation of the reference axes and computational cells. The time‐dependent discrete equations are then integrated in time via an implicit Runge–Kutta procedure that in this paper is proven to remain absolutely non‐linearly stable for the spatially‐discrete Euler and Navier–Stokes equations and shown to converge rapidly to steady states, with maximum Courant number exceeding 100 for the linearized version. Even on relatively coarse grids, the acoustics–convection upstream resolution algorithm generates essentially non‐oscillatory solutions for subsonic, transonic and supersonic flows, encompassing oblique‐ and interacting‐shock fields that converge within 40 time steps and reflect reference exact solutions. Copyright © 2005 John Wiley & Sons, Ltd.     

Derivation and characteristics analysis of an acoustics–convection upstream resolution algorithm for the two‐dimensional Euler and Navier–Stokes equations

The first of a two‐paper series, this paper introduces a new decomposition not of the hyperbolic flux vector but of the flux vector Jacobian. The paper then details for the Euler and Navier–Stokes equations an intrinsically infinite directional upstream‐bias formulation that rests on the mathematics and physics of multi‐dimensional acoustics and convection. Based upon characteristic velocities, this formulation introduces the upstream bias directly at the differential equation level, before the spatial discretization, within a characteristics‐bias governing system. Through a decomposition of the Euler flux divergence into multi‐dimensional acoustics and convection–acoustics components, this characteristics‐bias system induces consistent upstream bias along all directions of spatial wave propagation, with anisotropic variable‐strength upstreaming that correlates with the spatial distribution of characteristic velocities. Copyright © 2005 John Wiley & Sons, Ltd.     

A Mach‐uniform preconditioner for incompressible and subsonic flows

In this study, a novel Mach‐uniform preconditioning method is developed for the solution of Euler equations at low subsonic and incompressible flow conditions. In contrast to the methods developed earlier in which the conservation of mass equation is preconditioned, in the present method, the conservation of energy equation is preconditioned, which enforces the divergence free constraint on the velocity field even at the limiting case of incompressible, zero Mach number flows. Despite most preconditioners, the proposed Mach‐uniform preconditioning method does not have a singularity point at zero Mach number. The preconditioned system of equations preserves the strong conservation form of Euler equations for compressible flows and recovers the artificial compressibility equations in the case of zero Mach number. A two‐dimensional Euler solver is developed for validation and performance evaluation of the present formulation for a wide range of Mach number flows. The validation cases studied show the convergence acceleration, stability, and accuracy of the present Mach‐uniform preconditioner in comparison to the non‐preconditioned compressible flow solutions. The convergence acceleration obtained with the present formulation is similar to those of the well‐known preconditioned system of equations for low subsonic flows and to those of the artificial compressibility method for incompressible flows. Copyright © 2013 John Wiley & Sons, Ltd.     

A high‐order discontinuous Galerkin solver for low Mach number flows

In this work, we present a high‐order discontinuous Galerkin method (DGM) for simulating variable density flows at low Mach numbers. The corresponding low Mach number equations are an approximation of the compressible Navier–Stokes equations in the limit of zero Mach number. To the best of the authors'y knowledge, it is the first time that the DGM is applied to the low Mach number equations. The mixed‐order formulation is applied for spatial discretization. For steady cases, we apply the semi‐implicit method for pressure‐linked equation (SIMPLE) algorithm to solve the non‐linear system in a segregated manner. For unsteady cases, the solver is implicit in time using backward differentiation formulae, and the SIMPLE algorithm is applied to solve the non‐linear system in each time step. Numerical results for the following three test cases are shown: Couette flow with a vertical temperature gradient, natural convection in a square cavity, and unsteady natural convection in a tall cavity. Considering a fixed number of degrees of freedom, the results demonstrate the benefits of using higher approximation orders. Copyright © 2015 John Wiley & Sons, Ltd.     

A finite volume scheme preserving extremum principle for convection–diffusion equations on polygonal meshes

We propose a nonlinear finite volume scheme for convection–diffusion equation on polygonal meshes and prove that the discrete solution of the scheme satisfies the discrete extremum principle. The approximation of diffusive flux is based on an adaptive approach of choosing stencil in the construction of discrete normal flux, and the approximation of convection flux is based on the second‐order upwind method with proper slope limiter. Our scheme is locally conservative and has only cell‐centered unknowns. Numerical results show that our scheme can preserve discrete extremum principle and has almost second‐order accuracy. Copyright © 2017 John Wiley & Sons, Ltd.     

The third‐order polynomial method for two‐dimensional convection and diffusion

Using the upstream polynomial approximation a series of accurate two‐dimensional explicit numerical schemes is developed for the solution of the convection–diffusion equation. A third‐order polynomial approximation (TOP) of the convection term and a consistent second‐order approximation of the diffusion term are combined in a single‐step flux‐difference algorithm. Stability analysis confirms that the TOP‐12 scheme satisfies the CFL condition for two dimensions. Using smaller and narrower flux stencils compared to algorithms of similar accuracy, the TOP‐12 scheme is more efficient in terms of computations per single node. Numerical tests and comparison with other well‐known algorithms show a high performance of the developed schemes. Copyright © 2003 John Wiley & Sons, Ltd.     

An Enthalpy Control Volume Method for Transient Mass and Heat Transport with Solidification

A single domain enthalpy control volume method is developed for solving the coupled fluid flow and heat transfer with solidification problem arising from the continuous casting process. The governing equations consist of the continuity equation, the Navier–Stokes equations and the convection–diffusion equation. The formulation of the method is cast into the framework of the Petrov–Galerkin finite element method with a step test function across the control volume and locally constant approximation to the fluxes of heat and fluid. The use of the step test function and the constant flux approximation leads to the derivation of the exponential interpolating functions for the velocity and temperature fields within each control volume. The exponential fitting makes it possible to capture the sharp boundary layers around the solidification front. The method is then applied to investigate the effect of various casting parameters on the solidification profile and flow pattern of fluids in the casting process.     

Improvements to compressible Euler methods for low‐Mach number flows

In the present study improvements to numerical algorithms for the solution of the compressible Euler equations at low Mach numbers are investigated. To solve flow problems for a wide range of Mach numbers, from the incompressible limit to supersonic speeds, preconditioning techniques are frequently employed. On the other hand, one can achieve the same aim by using a suitably modified acoustic damping method. The solution algorithm presently under consideration is based on Roe's approximate Riemann solver [Roe PL. Approximate Riemann solvers, parameter vectors and difference schemes. Journal of Computational Physics 1981; 43 : 357–372] for non‐structured meshes. The numerical flux functions are modified by using Turkel's preconditioning technique proposed by Viozat [Implicit upwind schemes for low Mach number compressible flows. INRIA, Rapport de Recherche No. 3084, January 1997] for compressible Euler equations and by using a modified acoustic damping of the stabilization term proposed in the present study. These methods allow the compressible Euler equations at low‐Mach number flows to be solved, and they are consistent in time. The efficiency and accuracy of the proposed modifications have been assessed by comparison with experimental data and other numerical results in the literature. Copyright © 2000 John Wiley & Sons, Ltd.     

A multiscale control volume finite element method for advection–diffusion equations

We present a new stabilized method for advection–diffusion equations, which combines a control volume FEM formulation of the governing equations with a novel multiscale approximation of the total flux. The latter incorporates information about the exact solution that cannot be represented on the mesh. To define this flux, we solve the governing equations along suitable mesh segments under the assumption that the flux varies linearly along these segments. This procedure yields second‐order accurate fluxes on the edges of the mesh. Then, we use curl‐conforming elements of the same order to lift these edge fluxes into the mesh elements. In so doing, we obtain a stabilized control volume FEM formulation that is second‐order accurate and does not require mesh‐dependent stabilization parameters. Numerical convergence studies on uniform and nonuniform grids along with several standard advection tests illustrate the computational properties of the new method. Published 2015. This article is a U.S. Government work and is in the public domain in the USA.     

A kinetic energy preserving nodal discontinuous Galerkin spectral element method

In this work, we discuss the construction of a skew‐symmetric discontinuous Galerkin (DG) collocation spectral element approximation for the compressible Euler equations. Starting from the skew‐symmetric formulation of Morinishi, we mimic the continuous derivations on a discrete level to find a formulation for the conserved variables. In contrast to finite difference methods, DG formulations naturally have inter‐domain surface flux contributions due to the discontinuous nature of the approximation space. Thus, throughout the derivations we accurately track the influence of the surface fluxes to arrive at a consistent formulation also for the surface terms. The resulting novel skew‐symmetric method differs from the standard DG scheme by additional volume terms. Those volume terms have a special structure and basically represent the discretization error of the different product rules. We use the summation‐by‐parts (SBP) property of the Gauss–Lobatto‐based DG operator and show that the novel formulation is exactly conservative for the mass, momentum, and energy. Finally, an analysis of the kinetic energy balance of the standard DG discretization shows that because of aliasing errors, a nonzero transport source term in the evolution of the discrete kinetic energy mean value may lead to an inconsistent increase or decrease in contrast to the skew‐symmetric formulation. Furthermore, we derive a suitable interface flux that guarantees kinetic energy preservation in combination with the skew‐symmetric DG formulation. As all derivations require only the SBP property of the Gauss–Lobatto‐based DG collocation spectral element method operator and that the mass matrix is diagonal, all results for the surface terms can be directly applied in the context of multi‐domain diagonal norm SBP finite difference methods. Numerical experiments are conducted to demonstrate the theoretical findings. Copyright © 2014 John Wiley & Sons, Ltd.     

Method of relaxed streamline upwinding for hyperbolic conservation laws

In this work a new finite element based Method of Relaxed Streamline Upwinding is proposed to solve hyperbolic conservation laws. Formulation of the proposed scheme is based on relaxation system which replaces hyperbolic conservation laws by semi-linear system with stiff source term also called as relaxation term. The advantage of the semi-linear system is that the nonlinearity in the convection term is pushed towards the source term on right hand side which can be handled with ease. Six symmetric discrete velocity models are introduced in two dimensions which symmetrically spread foot of the characteristics in all four quadrants thereby taking information symmetrically from all directions. Proposed scheme gives exact diffusion vectors which are very simple. Moreover, the formulation is easily extendable from scalar to vector conservation laws. Various test cases are solved for Burgers equation (with convex and non-convex flux functions), Euler equations and shallow water equations in one and two dimensions which demonstrate the robustness and accuracy of the proposed scheme. New test cases are proposed for Burgers equation, Euler and shallow water equations. Exact solution is given for two-dimensional Burgers test case which involves normal discontinuity and series of oblique discontinuities. Error analysis of the proposed scheme shows optimal convergence rate. Moreover, spectral stability analysis gives implicit expression of critical time step.     

Numerical experiments with several variant WENO schemes for the Euler equations

Numerical experiments with several variants of the original weighted essentially non‐oscillatory (WENO) schemes (J. Comput. Phys. 1996; 126 :202–228) including anti‐diffusive flux corrections, the mapped WENO scheme, and modified smoothness indicator are tested for the Euler equations. The TVD Runge–Kutta explicit time‐integrating scheme is adopted for unsteady flow computations and lower–upper symmetric‐Gauss–Seidel (LU‐SGS) implicit method is employed for the computation of steady‐state solutions. A numerical flux of the variant WENO scheme in flux limiter form is presented, which consists of first‐order and high‐order fluxes and allows for a more flexible choice of low‐order schemes. Computations of unsteady oblique shock wave diffraction over a wedge and steady transonic flows over NACA 0012 and RAE 2822 airfoils are presented to test and compare the methods. Various aspects of the variant WENO methods including contact discontinuity sharpening and steady‐state convergence rate are examined. By using the WENO scheme with anti‐diffusive flux corrections, the present solutions indicate that good convergence rate can be achieved and high‐order accuracy is maintained and contact discontinuities are sharpened markedly as compared with the original WENO schemes on the same meshes. Copyright © 2008 John Wiley & Sons, Ltd.     

Improving Euler Computations at Low Mach Numbers

Abstract

This paper consists of two parts, both dealing with conditioning techniques for low-Mach-number Euler-flow computations, in which a multigrid technique is applied

In the first part, for subsonic flows and upwind-discretized linearized 1-D Euler equations, the smoothing behavior of multigrid-accelerated point Gauss-Seidel relaxation is investigated. Error decay by convection over domain boundaries is also discussed. A fix to poor convergence rates at low Mach numbers is sought by replacing the point relaxation applied to unconditioned Euler equations by locally implicit “time” stepping applied to preconditioned Euler equations. The locally implicit iteration step is optimized for good damping of high-frequency errors. Numerical inaccuracy at low Mach numbers is also addressed. In the present case it is not necessary to solve this accuracy problem

In the second part, insight is given into the conditions of derivative matrices to be inverted in point-relaxation methods for 1-D and 2-D, upwind-discretized Euler equations. Speed regimes are found where ill-conditioning of these matrices occurs, and 1-D flow equations appear to be less well-conditioned than 2-D flow equations. Fixes to the ill-conditioning follow more or less directly, when thinking of adding regularizing matrices to the ill-conditioned derivative matrices. A smoothing analysis is made of point Gauss-Seidel relaxation applied to discrete Euler equations conditioned by such an additive matrix. The method is successfully applied to a very low-subsonic,     

Multi‐dimensional wave‐oriented upwind schemes with reduced cross‐wind diffusion for flow in porous media

Single‐point upstream weighting schemes remain the reservoir simulation standard for integrating the convective components of the subsurface flow equations. While single‐point upstream weighting is attractive due to simplicity and small stencil, the accompanying first‐order error can contribute to a large amount of grid‐dependent numerical diffusion which smooths the numerical solution. New upwind schemes are introduced for reservoir simulation that significantly reduce cross‐wind diffusion, retain a local flux approximation with local conservation, and are free of spurious oscillations. Copyright © 2010 John Wiley & Sons, Ltd.     

SIMPLE‐type preconditioners for the Oseen problem

In this paper, the steady incompressible Navier–Stokes equations are discretized by the finite element method. The resulting systems of equations are solved by preconditioned Krylov subspace methods. Some new preconditioning strategies, both algebraic and problem dependent are discussed. We emphasize on the approximation of the Schur complement as used in semi implicit method for pressure‐linked equations‐type preconditioners. In the usual formulation, the Schur complement matrix and updates use scaling with the diagonal of the convection–diffusion matrix. We propose a variant of the SIMPLER preconditioner. Instead of using the diagonal of the convection–diffusion matrix, we scale the Schur complement and updates with the diagonal of the velocity mass matrix. This variant is called modified SIMPLER (MSIMPLER). With the new approximation, we observe a drastic improvement in convergence for large problems. MSIMPLER shows better convergence than the well‐known least‐squares commutator preconditioner which is also based on the diagonal of the velocity mass matrix. Copyright © 2008 John Wiley & Sons, Ltd.     

Composite high resolution localized relaxation scheme based on upwinding for hyperbolic conservation laws

In this work we present an upwind‐based high resolution scheme using flux limiters. Based on the direction of flow we choose the smoothness parameter in such a way that it leads to a truly upwind scheme without losing total variation diminishing (TVD) property for hyperbolic linear systems where characteristic values can be of either sign. Here we present and justify the choice of smoothness parameters. The numerical flux function of a high resolution scheme is constructed using wave speed splitting so that it results into a scheme that truly respects the physical hyperbolicity property. Bounds are given for limiter functions to satisfy TVD property. The proposed scheme is extended for non‐linear problems by using the framework of relaxation system that converts a non‐linear conservation law into a system of linear convection equations with a non‐linear source term. The characteristic speed of relaxation system is chosen locally on three point stencil of grid. This obtained relaxation system is solved using composite scheme technique, i.e. using a combination of proposed scheme with the conservative non‐standard finite difference scheme. Presented numerical results show higher resolution near discontinuity without introducing spurious oscillations. Copyright © 2008 John Wiley & Sons, Ltd.     

Adiabatic shock capturing in perfect gas hypersonic flows

This paper considers the streamline‐upwind Petrov/Galerkin (SUPG) method applied to the compressible Euler and Navier–Stokes equations in conservation‐variable form. The spatial discretization, including a modified approach for interpolating the inviscid flux terms in the SUPG finite element formulation, is briefly reviewed. Of particular interest is the behavior of the shock‐capturing operator, which is required to regularize the scheme in the presence of strong, shock‐induced gradients. A standard shock‐capturing operator that has been widely used in previous studies by several authors is presented and discussed. Specific modifications are then made to this standard operator that is designed to produce a more physically consistent discretization in the presence of strong shock waves. The actual implementation of the term in a finite‐dimensional approximation is also discussed. The behavior of the standard and modified scheme is then compared for several supersonic/hypersonic flows. The modified shock‐capturing operator is found to preserve enthalpy in the inviscid portion of the flowfield substantially better than the standard operator. Published in 2009 by John Wiley & Sons, Ltd.     

Anisotropic goal‐oriented error analysis for a third‐order accurate CENO Euler discretization

In this paper, a central essentially non‐oscillatory approximation based on a quadratic polynomial reconstruction is considered for solving the unsteady 2D Euler equations. The scheme is third‐order accurate on irregular unstructured meshes. The paper concentrates on a method for a metric‐based goal‐oriented mesh adaptation. For this purpose, an a priori error analysis for this central essentially non‐oscillatory scheme is proposed. It allows us to get an estimate depending on the polynomial reconstruction error. As a third‐order error is not naturally expressed in terms of a metric, we propose a least‐square method to approach a third‐order error by a quadratic term. Then an optimization problem for the best mesh metric is obtained and analytically solved. The resulting mesh optimality system is discretized and solved using a global unsteady fixed‐point algorithm. The method is applied to an acoustic propagation benchmark.     

Generalized lattice‐BGK concept for thermal and chemically reacting flows at low Mach numbers

The lattice‐BGK method has been extended by introducing additional, free parameters in the original formulation of the lattice‐BGK methods. The relationship between these parameters and the macroscopic moment equations is analysed by Taylor series and Chapman–Enskog expansion. The parameters are determined from the macroscopic moment equations by comparisons with the governing equations to be modelled. Extensions are presented for the Navier–Stokes equations at low Mach numbers in Cartesian or axisymmetric coordinates with constant or variable density, for scalar convection–diffusion equations and for equations of Poisson type. The generalized lattice‐BGK concept is demonstrated by two applications of chemical engineering. These are the computation of chemically reacting flow through an axisymmetric reactor and of the transport and deposition of particles to filters under the action of different forces. Copyright © 2006 John Wiley & Sons, Ltd.     

A simple incompressible flux splitting for sharp free surface capturing

This paper first applies a flux vector‐type splitting method based on the numerical speed of sound for computing incompressible single and multifluid flows. Here, a preconditioning matrix based on Chorin's artificial compressibility concept is used to modify the incompressible multifluid Navier–Stokes equations to be hyperbolic and density or volume fraction‐independent. The current approach can reduce eigenvalues disparity induced from density or volume fraction ratios and enhance numerical stability. Also, a simple convection‐pressure flux‐splitting method with high‐order essentially nonoscillatory‐type primitive variable extrapolations coupled with monotone upstream‐centered schemes for conservation laws‐type volume fraction recompressed reconstruction is used to maintain the preservation of sharp interface evolutions in multifluid flow simulations. Benchmark tests including a solid rotation test of a notched two‐dimensional cylinder, the evolution of spiral and rotational shapes of deformable circles, a dam breaking problem, and the Rayleigh–Taylor instability were chosen to validate the current incompressible multifluid methodology. An incompressible driven cavity was also chosen to check the robustness of the proposed method on the computation of single fluid incompressible flow problems. Copyright © 2011 John Wiley & Sons, Ltd.