A multiblock operator‐splitting algorithm for unsteady flows at all speeds in complex geometries

This paper describes a non‐iterative operator‐splitting algorithm for computing all‐speed flows in complex geometries. A pressure‐based algorithm is adopted as the base, in which pressure, instead of density, is a primary variable, thus allowing for a unified formulation for all Mach numbers. The focus is on adapting the method for (a) flows at all speeds, and (b) multiblock, non‐orthogonal, body‐fitted grids for very complex geometries. Key features of the formulation include special treatment of mass fluxes at control volume interfaces to avoid pressure–velocity decoupling for incompressible (low Mach number limit) flows and to provide robust pressure–velocity–density coupling for compressible (high‐speed) flows. The method is shown to be robust for all Mach number regimes for both steady and unsteady flows; it is found to be stable for CFL numbers of order ten, allowing large time steps to be taken for steady flows. Enhancements to the method which allow for stable solutions to be obtained on non‐orthogonal grids are also discussed. The method is found to be very reliable even in complex engineering applications such as unsteady rotor–stator interactions in turbulent, all‐speed turbomachinery flows. Copyright © 2004 John Wiley & Sons, Ltd.     

GPU‐accelerated direct numerical simulations of decaying compressible turbulence employing a GKM‐based solver

Gas Kinetic Method‐based flow solvers have become popular in recent years owing to their robustness in simulating high Mach number compressible flows. We evaluate the performance of the newly developed analytical gas kinetic method (AGKM) by Xuan et al. in performing direct numerical simulation of canonical compressible turbulent flow on graphical processing unit (GPU)s. We find that for a range of turbulent Mach numbers, AGKM results shows excellent agreement with high order accurate results obtained with traditional Navier–Stokes solvers in terms of key turbulence statistics. Further, AGKM is found to be more efficient as compared with the traditional gas kinetic method for GPU implementation. We present a brief overview of the optimizations performed on NVIDIA K20 GPU and show that GPU optimizations boost the speedup up‐to 40x as compared with single core CPU computations. Hence, AGKM can be used as an efficient method for performing fast and accurate direct numerical simulations of compressible turbulent flows on simple GPU‐based workstations. Copyright © 2016 John Wiley & Sons, Ltd.     

Accuracy of high‐order density‐based compressible methods in low Mach vortical flows

At low Mach numbers, Godunov‐type approaches, based on the method of lines, suffer from an accuracy problem. This paper shows the importance of using the low Mach number correction in Godunov‐type methods for simulations involving low Mach numbers by utilising a new, well‐posed, two‐dimensional, two‐mode Kelvin–Helmholtz test case. Four independent codes have been used, enabling the examination of several numerical schemes. The second‐order and fifth‐order accurate Godunov‐type methods show that the vortex‐pairing process can be captured on a low resolution with the low Mach number correction applied down to 0.002. The results are compared without the low Mach number correction and also three other methods, a Lagrange‐remap method, a fifth‐order accurate in space and time finite difference type method based on the wave propagation algorithm, and fifth‐order spatial and third‐order temporal accurate finite volume Monotone Upwind Scheme for Conservation Laws (MUSCL) approach based on the Godunov method and Simple Low Dissipation Advection Upstream Splitting Method (SLAU) numerical flux with low Mach capture property. The ability of the compressible flow solver of the commercial software, ANSYS FLUENT , in solving low Mach flows is also demonstrated for the two time‐stepping methods provided in the compressible flow solver, implicit and explicit. Results demonstrate clearly that a low Mach correction is required for all algorithms except the Lagrange‐remap approach, where dissipation is independent of Mach number. © 2013 Crown copyright. International Journal for Numerical Methods in Fluids. © 2013 John Wiley & Sons, Ltd.     

Spectral methods for compressible reactive flows

High order simulations are necessary in order to capture fine details in resolving supersonic reactive flows. However, high Mach number compressible flows feature sharp gradients and discontinuities, which present a challenge to successful simulations using high order methods. Spectral methods have proven a powerful tool in simulation of incompressible turbulent flows, and recent advances allow the application of spectral methods to compressible reactive flows. We review the recent advances in the theory and application of spectral methods which allow stable computations of discontinuous phenomena, and the recovery of high order information via postprocessing, and present applications of high Mach number reactive flows. To cite this article: D. Gottlieb, S. Gottlieb, C. R. Mecanique 333 (2005).     

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.     

An implicit mixed finite-volume–finite-element method for solving 3D turbulent compressible flows

The development of new aeronautic projects require accurate and efficient simulations of compressible flows in complex geometries. It is well known that most flows of interest are at least locally turbulent and that the modelling of this turbulence is critical for the reliability of the computations. A turbulence closure model which is both cheap and reasonably accurate is an essential part of a compressible code. An implicit algorithm to solve the 2D and 3D compressible Navier–Stokes equations on unstructured triangular/tetrahedral grids has been extended to turbulent flows. This numerical scheme is based on second-order finite element–finite volume discretization: the diffusive and source terms of the Navier–Stokes equations are computed using a finite element method, while the other terms are computed with a finite volume method. Finite volume cells are built around each node by means of the medians. The convective fluxes are evaluated with the approximate Riemann solver of Roe coupled with the van Albada limiter. The standard k–ϵ model has been introduced to take into account turbulence. Implicit integration schemes with efficient numerical methods (CGS, GMRES and various preconditioning techniques) have also been implemented. Our interest is to present the whole method and to demonstrate its limitations on some well-known test cases in three-dimensional geometries. © 1997 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.     

On higher‐order mixed FEM for low Mach number flows: application to a natural convection benchmark problem

We consider higher‐order mixed finite elements with continuous pressures for the computation of stationary compressible flows at low Mach number. The proposed approach is based on a fully coupled treatment of the governing equations and therefore, for steady‐state calculations, does not rely on time‐stepping techniques. The non‐linear problem is solved by means of a quasi‐Newton iteration. The strongly coupled system resulting from higher‐order discretization of the linearized equations requires adequate solvers. We propose a new scheme based on multigrid methods with varying FEM ansatz orders on the grid hierarchy as well as multiplicative smoothers based on blocking techniques. Computational results are described for a benchmark configuration including a flow with heat transfer in the low Mach number regime. Furthermore, the issue of anisotropic grids is addressed in that context. Copyright © 2003 John Wiley & Sons, Ltd.     

Cover Image

A simple, robust, and accurate HLLC-type Riemann solver for the compressible Euler equations at various Mach numbers is built. To cure shock instability of the HLLC solver at strong shocks, a pressure-control technique, which plays a role in limiting the propagation of erroneous pressure perturbation, is proposed. With an all Mach correction method for the compressible Euler system, the proposed method is further extended to compute flow problems at low Mach numbers. The proposed all Mach HLLC-type scheme has been implemented and used to compute a variety of flow problems ranging from hypersonic compressible to low Mach incompressible flow regimes. Various numerical results demonstrate that the obtained all Mach HLLC-type scheme is both accurate and stable for all speed ranges.     

A high‐order discontinuous Galerkin method for all‐speed flows

In this article, we present a discontinuous Galerkin (DG) method designed to improve the accuracy and efficiency of steady solutions of the compressible fully coupled Reynolds‐averaged Navier–Stokes and k ? ω turbulence model equations for solving all‐speed flows. The system of equations is iterated to steady state by means of an implicit scheme. The DG solution is extended to the incompressible limit by implementing a low Mach number preconditioning technique. A full preconditioning approach is adopted, which modifies both the unsteady terms of the governing equations and the dissipative term of the numerical flux function by means of a new preconditioner, on the basis of a modified version of Turkel's preconditioning matrix. At sonic speed the preconditioner reduces to the identity matrix thus recovering the non‐preconditioned DG discretization. An artificial viscosity term is added to the DG discretized equations to stabilize the solution in the presence of shocks when piecewise approximations of order of accuracy higher than 1 are used. Moreover, several rescaling techniques are implemented in order to overcome ill‐conditioning problems that, in addition to the low Mach number stiffness, can limit the performance of the flow solver. These approaches, through a proper manipulation of the governing equations, reduce unbalances between residuals as a result of the dependence on the size of elements in the computational mesh and because of the inherent differences between turbulent and mean‐flow variables, influencing both the evolution of the Courant Friedrichs Lewy (CFL) number and the inexact solution of the linear systems. The performance of the method is demonstrated by solving three turbulent aerodynamic test cases: the flat plate, the L1T2 high‐lift configuration and the RAE2822 airfoil (Case 9). The computations are performed at different Mach numbers using various degrees of polynomial approximations to analyze the influence of the proposed numerical strategies on the accuracy, efficiency and robustness of a high‐order DG solver at different flow regimes. Copyright © 2014 John Wiley & Sons, Ltd.     

Multiple pressure variables methods for fluid flow at all Mach numbers

In this paper we present a method for the simulation of incompressible as well as compressible unsteady flows. At first we discuss three different forms, i.e. a primitive‐, conservative‐ and a semi‐conservative form of the governing equations. We use a semi‐implicit time integration in such a fashion that the stability is guaranteed independently of the speed of sound and the resulting method is independent of the Mach number range. Moreover, with the application of the so‐called multiple pressure variables (MPV) approach the difficulties with the pressure term can be circumvented as in the incompressible limit the hydrodynamic pressure decouples from the equation of state. Increasing approximation errors in the low Mach number regime are avoided. As a result, the proposed algorithm can also simulate incompressible flows as limit for zero Mach number. Copyright © 2005 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.     

Development and validation of a SUPG finite element scheme for the compressible Navier–Stokes equations using a modified inviscid flux discretization

This paper considers the streamline‐upwind Petrov–Galerkin (SUPG) method applied to the unsteady compressible 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, and the second‐order accurate time discretization are presented. The numerical method is discussed in detail. The performance of the algorithm is then investigated by considering inviscid flow past a circular cylinder. Validation of the finite element formulation via comparisons with experimental data for high‐Mach number perfect gas laminar flows is presented, with a specific focus on comparisons with experimentally measured skin friction and convective heat transfer on a 15° compression ramp. Copyright © 2007 John Wiley & Sons, Ltd.     

Recent developments in the computation of compressible low Mach number flows

The computation of low speed flows can usualy be performed by incompressible models or various Low Mach number approximations (Boussinesq, Anelastic, etc). However, there is a large number of flows where although the velocities are small, compressible effects cannot be ignored and the use of full compressible models is required.Unfortunately, standard finite element or finite volume method experiences various difficulties in the computation of these low Mach number flows. This talk will explain the origin of these difficulties and describe some techniques to circumvent them.     

亚、跨、超音速及不可压流动的数值分析方法的研究

为了对亚、跨、超音速及不可压无粘流动进行数值模拟,将LU－SGS方法与预处理方法结合,给出了PLU－SGS方法。方程离散基于有限体积法,采用高阶精度AUSMPW格式。方程求解采用了特征边界条件。通过典型算例的数值试验对比分析,表明PLU－SGS方法可以有效地对亚、跨、超音速及不可压流动进行数值模拟,并具有较高的计算精度和收敛速度。     

Numerical simulation of cavitating flow in 2D and 3D inducer geometries

A computational method is proposed to simulate 3D unsteady cavitating flows in spatial turbopump inducers. It is based on the code FineTurbo, adapted to take into account two‐phase flow phenomena. The initial model is a time‐marching algorithm devoted to compressible flow, associated with a low‐speed preconditioner to treat low Mach number flows. The presented work covers the 3D implementation of a physical model developed in LEGI for several years to simulate 2D unsteady cavitating flows. It is based on a barotropic state law that relates the fluid density to the pressure variations. A modification of the preconditioner is proposed to treat efficiently as well highly compressible two‐phase flow areas as weakly compressible single‐phase flow conditions. The numerical model is applied to time‐accurate simulations of cavitating flow in spatial turbopump inducers. The first geometry is a 2D Venturi type section designed to simulate an inducer blade suction side. Results obtained with this simple test case, including the study of its general cavitating behaviour, numerical tests, and precise comparisons with previous experimental measurements inside the cavity, lead to a satisfactory validation of the model. A complete three‐dimensional rotating inducer geometry is then considered, and its quasi‐static behaviour in cavitating conditions is investigated. Numerical results are compared to experimental measurements and visualizations, and a promising agreement is obtained. Copyright © 2004 John Wiley & Sons, Ltd.     

Simulation of thermoacoustic waves by a pressure-based algorithm for compressible flows

A modified SIMPLEC method which can solve compressible flows at low Mach number is introduced and used to study thermoacoustic waves induced by a rapid change of temperature at a solid wall and alternatingdirection flows generated by thermoacoustic effects in a tapered resonator.The results indicate that the algorithm adopted in this paper can be used for calculating compressible flows and thermoacoustic waves.It is found that the pressure and velocity in the resonator behave as standing waves,and the tapered resonator can suppress highfrequency harmonic waves as observed in a cylindrical resonator.     

A matrix‐free preconditioned Newton/GMRES method for unsteady Navier–Stokes solutions

The unsteady compressible Reynolds‐averaged Navier–Stokes equations are discretized using the Osher approximate Riemann solver with fully implicit time stepping. The resulting non‐linear system at each time step is solved iteratively using a Newton/GMRES method. In the solution process, the Jacobian matrix–vector products are replaced by directional derivatives so that the evaluation and storage of the Jacobian matrix is removed from the procedure. An effective matrix‐free preconditioner is proposed to fully avoid matrix storage. Convergence rates, computational costs and computer memory requirements of the present method are compared with those of a matrix Newton/GMRES method, a four stage Runge–Kutta explicit method, and an approximate factorization sub‐iteration method. Effects of convergence tolerances for the GMRES linear solver on the convergence and the efficiency of the Newton iteration for the non‐linear system at each time step are analysed for both matrix‐free and matrix methods. Differences in the performance of the matrix‐free method for laminar and turbulent flows are highlighted and analysed. Unsteady turbulent Navier–Stokes solutions of pitching and combined translation–pitching aerofoil oscillations are presented for unsteady shock‐induced separation problems associated with the rotor blade flows of forward flying helicopters. Copyright © 2000 John Wiley & Sons, Ltd.     

Computational study of unsteady turbulent flows around oscillating and ramping aerofoils

The aim of this work is to computationally investigate subsonic and transonic turbulent flows around oscillating and ramping aerofoils under dynamic‐stall conditions. The investigation is based on a high‐resolution Godunov‐type method and several turbulence closures. The Navier–Stokes and turbulence transport equations are solved in a strongly coupled fashion via an implicit‐unfactored scheme. We present results from several computations of flows around oscillating and ramping aerofoils at various conditions in order to (i) assess the accuracy of different turbulence models and (ii) contribute towards a better understanding of dynamic‐stall flows. The results show that the employed non‐linear eddy‐viscosity model generally improves the accuracy of the computations compared to linear models, but at low incidence angles the Spalart–Allmaras one‐equation model was found to provide adequate results. Further, the computations reveal strong similarities between laminar and high‐Reynolds number dynamic‐stall flows as well as between ramping and oscillating aerofoil cases. Investigation of the Mach number effects on dynamic‐stall reveals a delay of the stall angle within a range of Mach numbers. Investigation of the reduced frequency effects suggests the existence of an (almost) linear variation between pitch rate and stall angle, with higher slope at lower pitch rates. The pitch rate affects both the onset of dynamic‐stall as well as the evolution of the associated vortical structures. Copyright © 2003 John Wiley & Sons, Ltd.     

A sharp‐interface immersed boundary framework for simulations of high‐speed inviscid compressible flows

A new finite‐volume flow solver based on the hybrid Cartesian immersed boundary (IB) framework is developed for the solution of high‐speed inviscid compressible flows. The IB method adopts a sharp‐interface approach, wherein the boundary conditions are enforced on the body geometry itself. A key component of the present solver is a novel reconstruction approach, in conjunction with inverse distance weighting, to compute the solutions in the vicinity of the solid‐fluid interface. We show that proposed reconstruction leads to second‐order spatial accuracy while also ensuring that the discrete conservation errors diminish linearly with grid refinement. Investigations of supersonic and hypersonic inviscid flows over different geometries are carried out for an extensive validation of the proposed flow solver. Studies on cylinder lift‐off and shape optimisation in supersonic flows further demonstrate the efficacy of the flow solver for computations with moving and shape‐changing geometries. These studies conclusively highlight the capability of the proposed IB methodology as a promising alternative for robust and accurate computations of compressible fluid flows on nonconformal Cartesian meshes.