UNIFIED COMPUTATION OF FLOW WITH COMPRESSIBLE AND INCOMPRESSIBLE FLUID BASED ON ROE＇S SCHEME

A unified numerical scheme for the solutions of the compressible and incompressible Navier-Stokes equations is investigated based on a time-derivative preconditioning algorithm. The primitive variables are pressure, velocities and temperature. The time integration scheme is used in conjunction with a finite volume discretization. The preconditioning is coupled with a high order implicit upwind scheme based on the definition of a Roe＇s type matrix. Computational capabilities are demonstrated through computations of high Mach number, middle Mach number, very low Mach number, and incompressible flow. It has also been demonstrated that the discontinuous surface in flow field can be captured for the implementation Roe＇s scheme.     

Low-diffusion preconditioning scheme for numerical simulation of low-speed flows past airfoil

The preconditioning technique can address the stiffness of a low Mach number flow, while its stability is poor. Based on the conventional preconditioning method of Roe's scheme, a low-diffusion scheme is proposed. An adjustable parameter is introduced to control numerical dissipation, especially over the dissipation in the boundary layer and extremely in a low speed region. Numerical simulations of the low Mach number and low Reynolds number flows past a cylinder and the low Mach number and high Reynolds number flows past NACA0012 and S809 airfoils are performed to validate the new scheme. Results of the three tests well agree with experimental data, showing the applicability of the proposed scheme to low Mach number flow simulations.     

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 accurate discontinuous Galerkin finite element method for laminar low Mach number flows

In this paper we present a discontinuous Galerkin (DG) method designed to improve the accuracy and efficiency of laminar flow simulations at low Mach numbers using an implicit scheme. The algorithm is based on the flux preconditioning approach, which modifies only the dissipative terms of the numerical flux. This formulation is quite simple to implement in existing implicit DG codes, it overcomes the time‐stepping restrictions of explicit multistage algorithms, is consistent in time and thus applicable to unsteady flows. The performance of the method is demonstrated by solving the flow around a NACA0012 airfoil and on a flat plate, at different low Mach numbers using various degrees of polynomial approximations. Computations with and without flux preconditioning are performed on different grid topologies to analyze the influence of the spatial discretization on the accuracy of the DG solutions at low Mach numbers. The time accurate solution of unsteady flow is also demonstrated by solving the vortex shedding behind a circular cylinder at the Reynolds number of 100. Copyright © 2012 John Wiley & Sons, Ltd.     

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.     

Continuous and discrete adjoint approaches for aerodynamic shape optimization with low Mach number preconditioning

Discrete and continuous adjoint approaches for use in aerodynamic shape optimization problems at all flow speeds are developed and assessed. They are based on the Navier–Stokes equations with low Mach number preconditioning. By alleviating the large disparity between acoustic waves and fluid speeds, the preconditioned flow and adjoint equations are numerically solved with affordable CPU cost, even at the so‐called incompressible flow conditions. Either by employing the adjoint to the preconditioned flow equations or by preconditioning the adjoint to the ‘standard’ flow equations (under certain conditions the two formulations become equivalent, as proved in this paper), efficient optimization methods with reasonable cost per optimization cycle, even at very low Mach numbers, are derived. During the mathematical development, a couple of assumptions are made which are proved to be harmless to the accuracy in the computed gradients and the effectiveness of the optimization method. The proposed approaches are validated in inviscid and viscous flows in external aerodynamics and turbomachinery flows at various Mach numbers. Copyright © 2007 John Wiley & Sons, Ltd.     

Numerical simulation of a pulsating flow arising over an axisymmetric spiked blunt body at Mach 2.21 and Mach 6.00

The present paper reports a numerical simulation of the supersonic/hypersonic unsteady flow over a spiked blunt body. Axisymmetric compressible Navier-Stokes equations are solved using a high-resolution unfactored implicit upwind Roe's scheme and a time-accurate pseudo-time method is employed for advancing in time. Unsteady flows arising at Mach 2.21 and Mach 6.00 around a spiked cylinder are simulated and the computational results are compared with measurements. The simulated results are used to increase understanding of the mechanisms of the flow. Received 28 September 1999 / Accepted 26 July 2000     

Role of the momentum interpolation mechanism of the Roe scheme in shock instability

The shock instability phenomenon is a well‐known problem for hypersonic flow computation by the shock‐capturing Roe scheme. The pressure checkerboard is another well‐known problem for low‐Mach‐number flow computation. The momentum interpolation method (MIM) is necessary for low‐Mach‐number flows to suppress the pressure checkerboard problem, and the pressure‐difference‐driven modification for cell face velocity can be regarded as a version of the MIM by subdividing the numerical dissipation of the Roe scheme. In this paper, MIM has been discovered through analysis and numerical tests to have the most important function in shock instability. MIM should be completely removed for nonlinear flows. However, the unexpected MIM is activated on the cell face nearly parallel to the flow for the high‐Mach‐number flows or low‐Mach‐number cells in numerical shock. Therefore, MIM should be retained for low‐Mach‐number flows and be completely removed for high‐Mach‐number flows and low‐Mach‐number cells in numerical shock. For such conditions, two coefficients are designed on the basis of the local Mach number and a shock detector. Thereafter, the improved Roe scheme is proposed. This scheme considers the requirement of MIM for incompressible and compressible flows, and is validated for good performance of numerical tests. An acceptable result can also be obtained with only the Mach number coefficient for general practical computation. Therefore, the objective of decreasing rather than increasing numerical dissipation to cure shock instability can be achieved with simple modification. Moreover, the mechanism of shock instability has been profoundly understood, in which MIM plays the most important role, although it is not the only factor. Copyright © 2016 John Wiley & Sons, Ltd.     

Finite Element Computations of Flow around a Wall-mounted Cube

In this paper, we present the application of a finite element scheme to full three-dimensional incompressible flow around a cube mounted on the wall in a channel. This scheme is based on the Petrov-Galerkin weak formulation using exponential weighting functions. The incompressible Navier-Stokes equations are numerically integrated in time by using a fractional step strategy with a second-order accurate Adams-Bashforth scheme. The workability and validity of the present approach are demonstrated through the results of streamlines and pressure coefficients in the flow field up to high Reynolds number regimes.     

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.     

Euler solutions for aerodynamic inverse shape design

Contributions to the aerodynamics development have to be involved to achieve an increase in quality, reducing time and computer costs. Therefore, this work develops an optimization method based on the finite volume explicit Runge–Kutta multi‐stage scheme with central spatial discretization in combination with multigrid and preconditioning. The multigrid approach includes local time‐stepping and residual smoothing. Such a method allows getting the goal of compressible and almost incompressible solution of fluid flows, having a rate of convergence almost independent from the Mach number. Numerical tests are carried out for the NACA 0012 and 0009 airfoils and three‐dimensional wings based on NACA profiles for Mach‐numbers ranging from 0.8 to 0.002 using the Euler equations. These calculations are found to compare favorably with experimental and numerical data available in the literature. Besides, it is worth pointing out that these results build on earlier ones when finding appropriate new three‐dimensional aerodynamical geometries. Copyright © 2004 John Wiley & Sons, Ltd.     

Implicit preconditioned high‐order compact scheme for the simulation of the three‐dimensional incompressible Navier–Stokes equations with pseudo‐compressibility method

This paper combines the pseudo‐compressibility procedure, the preconditioning technique for accelerating the time marching for stiff hyperbolic equations, and high‐order accurate central compact scheme to establish the code for efficiently and accurately solving incompressible flows numerically based on the finite difference discretization. The spatial scheme consists of the sixth‐order compact scheme and 10th‐order numerical filter operator for guaranteeing computational stability. The preconditioned pseudo‐compressible Navier–Stokes equations are marched temporally using the implicit lower–upper symmetric Gauss–Seidel time integration method, and the time accuracy is improved by the dual‐time step method for the unsteady problems. The efficiency and reliability of the present procedure are demonstrated by applications to Taylor decaying vortices phenomena, double periodic shear layer rolling‐up problem, laminar flow over a flat plate, low Reynolds number unsteady flow around a circular cylinder at Re = 200, high Reynolds number turbulence flow past the S809 airfoil, and the three‐dimensional flows through two 90°curved ducts of square and circular cross sections, respectively. It is found that the numerical results of the present algorithm are in good agreement with theoretical solutions or experimental data. Copyright © 2011 John Wiley & Sons, Ltd.     

A preconditioned dual time-stepping procedure coupled with matrix-free LU-SGS scheme for unsteady low speed viscous flows with moving objects

In this paper an effective method is developed to solve unsteady low speed viscous flow problems with moving objects by using the governing equations of compressible fluids. The method is based on a dual time-stepping scheme, combined with low Mach number preconditioning and an implicit matrix-free Lower-Upper Symmetric Gauss-Seidel iteration on unstructured dynamic meshes. Because preconditioning modifies the governing equations, that induces the change of system's eigenvalues and eigenvectors, characteristic boundary conditions are also modified to suit the preconditioned characteristic system. Several test cases are simulated, including an in-line oscillating cylinder in a fluid at rest, flow over a flapping NACA0014 airfoil and low speed flow past a flapping-wing micro-air vehicle. Compared with experimental results whenever possible, the computed results indicate that this algorithm shows satisfactory improvement of solution efficiency and accuracy for low speed flow problems.     

Application of time preconditioning and high‐order compact discretization method for low Mach number flows

The paper describes a combination of a preconditioning method with a high‐order compact discretization scheme for the purpose of solving the compressible Navier–Stokes equations in moderate and low Mach number regimes. When combined with properly modified characteristic boundary conditions, the proposed approach is very efficient from the point of view of convergence acceleration and accuracy of the results. The computations were performed in typical benchmark cases including the Burggraf flow for which an analytical solution exists, the flow over a backward facing step, and also the flow in 2D and 3D shear‐driven cavities. Depending on the test case, the results were obtained for the Mach number in the range M = 0.001 ? 0.5 and the Reynolds number Re = 1 ? 1000. Copyright © 2012 John Wiley & Sons, Ltd.     

High order numerical simulation of the thermal load on a lobed strut injector for scramjet applications

In this paper, the thermal load on an actively cooled lobed strut injector for scramjet (supersonic combustion ramjet) applications is investigated numerically. This requires coupled simulations of the strut internal and external flow fields together with the heat conduction in the solid injector body. In order to achieve a fast mixing, the lobed strut is positioned at the channel axis to inject hydrogen into the core of a Mach 3 air stream. There it is exposed to the extremely high temperatures of the high speed flow. While the external air and hydrogen flows are supersonic, the strut internal hydrogen flow is mainly subsonic, in some regions at very low Mach numbers. To enable a simulation of the internal flow field which ranges from very low to very high Mach numbers (approximately Mach 2.25 at the nozzle exit), a preconditioning technique is employed. The compressible finite‐volume scheme uses a spatially fourth order multi‐dimensional limiting process discretization, which is used here for a first time to simulate a geometrically and fluid mechanically highly complex problem. It will be demonstrated that besides its high accuracy the multi‐dimensional limiting process scheme is numerically stable even in case of demanding practical applications. The coupled simulation of the lobed strut injector delivers unique insight into the flow phenomena inside and outside the strut, the heat fluxes, the temperature distribution in the solid material, the required hydrogen mass flux with respect to cooling requirements and details concerning the conditions at the exit of the injector. Copyright © 2016 John Wiley & Sons, Ltd.     

Incompressible Limit for the Compressible Flow of Liquid Crystals

The connection between the compressible flow of liquid crystals with low Mach number and the incompressible flow of liquid crystals is studied in a bounded domain. In particular, the convergence of weak solutions of the compressible flow of liquid crystals to the weak solutions of the incompressible flow of liquid crystals is proved when the Mach number approaches zero; that is, the incompressible limit is justified for weak solutions in a bounded domain.     

Comparison of low Mach number models for natural convection problems

 We investigate in this paper two numerical methods for solving low Mach number compressible flows and their application to single-phase natural convection flow problems. The first method is based on an asymptotic model of the Navier–Stokes equations valid for small Mach numbers, whereas the second is based on the full compressible Navier–Stokes equations with particular care given to the discretization at low Mach numbers. These models are more general than the Boussinesq incompressible flow model, in the sense that they are valid even for cases in which the fluid is subjected to large temperature differences, that is when the compressibility of the fluid manifests itself through low Mach number effects. Numerical solutions are computed for a series of test problems with fixed Rayleigh number and increasing temperature differences, as well as for varying Rayleigh number for a given temperature difference. Numerical difficulties associated with low Mach number effects are discussed, as well as the accuracy of the approximations. Received on 17 January 2000     

Numerical schemes for low Mach wave breaking

In this work, we describe a finite volume scheme for the computation of incompressible air–water flows. We use an artificial compressibility approach that permits us to use a completely explicit scheme. We describe successively the low Mach preconditioning of the scheme, the Riemann solver and then the non-conservative approach that is used to suppress velocity-pressure oscillations, the second order extensions and the parallel implementation. Then this is applied to the simulation of the breaking of a wave on a 15% slope.     

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

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