A CFD Euler solver from a physical acoustics–convection flux Jacobian decomposition

This paper introduces a continuum, i.e. non‐discrete, upstream‐bias formulation that rests on the physics and mathematics of acoustics and convection. The formulation induces the upstream‐bias at the differential equation level, within a characteristics‐bias system associated with the Euler equations with general equilibrium equations of state. For low subsonic Mach numbers, this formulation returns a consistent upstream‐bias approximation for the non‐linear acoustics equations. For supersonic Mach numbers, the formulation smoothly becomes an upstream‐bias approximation of the entire Euler flux. With the objective of minimizing induced artificial diffusion, the formulation non‐linearly induces upstream‐bias, essentially locally, in regions of solution discontinuities, whereas it decreases the upstream‐bias in regions of solution smoothness. The discrete equations originate from a finite element discretization of the characteristic‐bias system and are integrated in time within a compact block tridiagonal matrix statement by way of an implicit non‐linearly stable Runge–Kutta algorithm for stiff systems. As documented by several computational results that reflect available exact solutions, the acoustics–convection solver induces low artificial diffusion and generates essentially non‐oscillatory solutions that automatically preserve a constant enthalpy, as well as smoothness of both enthalpy and mass flux across normal shocks. Copyright © 1999 John Wiley & Sons, Ltd.     

A momentum coupling method for the unsteady incompressible Navier–Stokes equations on the staggered grid

A new numerical method is developed to efficiently solve the unsteady incompressible Navier–Stokes equations with second-order accuracy in time and space. In contrast to the SIMPLE algorithms, the present formulation directly solves the discrete x- and y-momentum equations in a coupled form. It is found that the present implicit formulation retrieves some cross convection terms overlooked by the conventional iterative methods, which contribute to accuracy and fast convergence. The finite volume method is applied on the fully staggered grid to solve the vector-form momentum equations. The preconditioned conjugate gradient squared method (PCGS) has proved very efficient in solving the associate linearized large, sparse block-matrix system. Comparison with the SIMPLE algorithm has indicated that the present momentum coupling method is fast and robust in solving unsteady as well as steady viscous flow problems. © 1998 John Wiley & Sons, Ltd.     

Temporal accuracy analysis of phase change convection simulations using the JFNK‐SIMPLE algorithm

The incompressible Navier–Stokes and energy conservation equations with phase change effects are applied to two benchmark problems: (1) non‐dimensional freezing with convection; and (2) pure gallium melting. Using a Jacobian‐free Newton–Krylov (JFNK) fully implicit solution method preconditioned with the SIMPLE (Numerical Heat Transfer and Fluid Flow. Hemisphere: New York, 1980) algorithm using centred discretization in space and three‐level discretization in time converges with second‐order accuracy for these problems. In the case of non‐dimensional freezing, the temporal accuracy is sensitive to the choice of velocity attenuation parameter. By comparing to solutions with first‐order backward Euler discretization in time, it is shown that the second‐order accuracy in time is required to resolve the fine‐scale convection structure during early gallium melting. Qualitative discrepancies develop over time for both the first‐order temporal discretized simulation using the JFNK‐SIMPLE algorithm that converges the nonlinearities and a SIMPLE‐based algorithm that converges to a more common mass balance condition. The discrepancies in the JFNK‐SIMPLE simulations using only first‐order rather than second‐order accurate temporal discretization for a given time step size appear to be offset in time. Copyright © 2007 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.     

二维对流扩散方程的高精度全隐式多重网格方法

提出了数值求解二维非定常变系数对流扩散方程的一种时间二阶、空间四阶精度的三层全隐紧致差分格式。为了加快迭代求解隐格式时在每一个时间步上的收敛速度，采用多重网格加速技术，建立了适用于本文高精度金隐紧致格式的多重网格算法。数值实验结果验证了本文方法的精确性、稳定性和对高网格雷诺数问题的强适应性。     

Accelerated convergence of the numerical simulation of incompressible flow in general curvilinear co‐ordinates by discretizations on the double‐staggered grids

The convergence rate of a methodology for solving incompressible flow in general curvilinear co‐ordinates is analyzed. Double‐staggered grids (DSGs), each defined by the same boundaries as the physical domain, are used for discretization. Both grids are MAC quadrilateral meshes with scalar variables (pressure, temperature, etc.) arranged at the center and the Cartesian velocity components at the middle of the sides of the mesh cells. The problem was checked against benchmark solutions of natural convection in a squeezed cavity, heat transfer in concentric horizontal cylindrical annuli, and a hot cylinder in a duct. Poisson's pressure‐correction equations that arise from the SIMPLE‐like procedure are solved by several methods: successive overrelaxation, symmetric overrelaxation, modified incomplete factorization preconditioner, conjugate gradient (CG), and CG with preconditioner. A genetic algorithm was developed to solve problems of numerical optimization of SIMPLE‐like calculation time in a space of iteration numbers and relaxation parameters. The application provides a means of making an unbiased comparison between the DSGs method and the widely used interpolation method. Furthermore, the convergence rate was demonstrated by application to the calculation of natural convection heat transfer in concentric horizontal cylindrical annuli. Calculation times when DSGs were used were 2–10 times shorter than those achieved by interpolation. With the DSGs method, calculation time increases slightly with increasing non‐orthogonality of the grids, whereas an interpolation method calls for very small iteration parameters that lead to unacceptable calculation times. Copyright © 2007 John Wiley & Sons, Ltd.     

A comparison of time integration methods in an unsteady low‐Reynolds‐number flow

This paper describes three different time integration methods for unsteady incompressible Navier–Stokes equations. Explicit Euler and fractional‐step Adams–Bashford methods are compared with an implicit three‐level method based on a steady‐state SIMPLE method. The implicit solver employs a dual time stepping and an iteration within the time step. The spatial discretization is based on a co‐located finite‐volume technique. The influence of the convergence limits and the time‐step size on the accuracy of the predictions are studied. The efficiency of the different solvers is compared in a vortex‐shedding flow over a cylinder in the Reynolds number range of 100–1600. A high‐Reynolds‐number flow over a biconvex airfoil profile is also computed. The computations are performed in two dimensions. At the low‐Reynolds‐number range the explicit methods appear to be faster by a factor from 5 to 10. In the high‐Reynolds‐number case, the explicit Adams–Bashford method and the implicit method appear to be approximately equally fast while yielding similar results. Copyright © 2002 John Wiley & Sons, Ltd.     

A comparative study of high‐order variable‐property segregated algorithms for unsteady low Mach number flows

In this paper, five different algorithms are presented for the simulation of low Mach flows with large temperature variations, based on second‐order central‐difference or fourth‐order compact spatial discretization and a pressure projection‐type method. A semi‐implicit three‐step Runge–Kutta/Crank–Nicolson or second‐order iterative scheme is used for time integration. The different algorithms solve the coupled set of governing scalar equations in a decoupled segregate manner. In the first algorithm, a temperature equation is solved and density is calculated from the equation of state, while the second algorithm advances the density using the differential form of the equation of state. The third algorithm solves the continuity equation and the fourth algorithm solves both the continuity and enthalpy equation in conservative form. An iterative decoupled algorithm is also proposed, which allows the computation of the fully coupled solution. All five algorithms solve the momentum equation in conservative form and use a constant‐ or variable‐coefficient Poisson equation for the pressure. The efficiency of the fourth‐order compact scheme and the performances of the decoupling algorithms are demonstrated in three flow problems with large temperature variations: non‐Boussinesq natural convection, channel flow instability, flame–vortex interaction. Copyright © 2010 John Wiley & Sons, Ltd.     

A time‐accurate pseudo‐compressibility approach based on an unstructured hybrid finite volume technique applied to unsteady turbulent premixed flame propagation

A preconditioning approach based on the artificial compressibility formulation is extended to solve the governing equations for unsteady turbulent reactive flows with heat release, at low Mach numbers, on an unstructured hybrid grid context. Premixed reactants are considered and a flamelet approach for combustion modelling is adopted using a continuous quenched mean reaction rate. An overlapped cell‐vertex finite volume method is adopted as a discretisation scheme. Artificial dissipation terms for hybrid grids are explicitly added to ensure a stable, discretised set of equations. A second‐order, explicit, hybrid Runge–Kutta scheme is applied for the time marching in pseudo‐time. A time derivative of the dependent variable is added to recover the time accuracy of the preconditioned set of equations. This derivative is discretised by an implicit, second‐order scheme. The resulting scheme is applied to the calculation of an infinite planar (one‐dimensional) turbulent premixed flame propagating freely in reactants whose turbulence is supposed to be frozen, homogeneous and isotropic. The accuracy of the results obtained with the proposed method proves to be excellent when compared to the data available in the literature. Copyright © 2004 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,     

Three‐dimensional numerical modelling of free surface flows with non‐hydrostatic pressure

A three‐dimensional numerical model is developed for incompressible free surface flows. The model is based on the unsteady Reynolds‐averaged Navier–Stokes equations with a non‐hydrostatic pressure distribution being incorporated in the model. The governing equations are solved in the conventional sigma co‐ordinate system, with a semi‐implicit time discretization. A fractional step method is used to enable the pressure to be decomposed into its hydrostatic and hydrodynamic components. At every time step one five‐diagonal system of equations is solved to compute the water elevations and then the hydrodynamic pressure is determined from a pressure Poisson equation. The model is applied to three examples to simulate unsteady free surface flows where non‐hydrostatic pressures have a considerable effect on the velocity field. Emphasis is focused on applying the model to wave problems. Two of the examples are about modelling small amplitude waves where the hydrostatic approximation and long wave theory are not valid. The other example is the wind‐induced circulation in a closed basin. The numerical solutions are compared with the available analytical solutions for small amplitude wave theory and very good agreement is obtained. Copyright © 2002 John Wiley & Sons, Ltd.     

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.     

Two‐dimensional prediction of time dependent,turbulent flow around a square cylinder confined in a channel

This paper presents two‐dimensional and unsteady RANS computations of time dependent, periodic, turbulent flow around a square block. Two turbulence models are used: the Launder–Sharma low‐Reynolds number k–ε model and a non‐linear extension sensitive to the anisotropy of turbulence. The Reynolds number based on the free stream velocity and obstacle side is Re=2.2×10⁴. The present numerical results have been obtained using a finite volume code that solves the governing equations in a vertical plane, located at the lateral mid‐point of the channel. The pressure field is obtained with the SIMPLE algorithm. A bounded version of the third‐order QUICK scheme is used for the convective terms. Comparisons of the numerical results with the experimental data indicate that a preliminary steady solution of the governing equations using the linear k–ε does not lead to correct flow field predictions in the wake region downstream of the square cylinder. Consequently, the time derivatives of dependent variables are included in the transport equations and are discretized using the second‐order Crank–Nicolson scheme. The unsteady computations using the linear and non‐linear k–ε models significantly improve the velocity field predictions. However, the linear k–ε shows a number of predictive deficiencies, even in unsteady flow computations, especially in the prediction of the turbulence field. The introduction of a non‐linear k–ε model brings the two‐dimensional unsteady predictions of the time‐averaged velocity and turbulence fields and also the predicted values of the global parameters such as the Strouhal number and the drag coefficient to close agreement with the data. Copyright © 2009 John Wiley & Sons, Ltd.     

A Fast,Matrix-free Implicit Method for Computing Low Mach Number Flows on Unstructured Grids

A fast, matrix-free implicit method has been developed to solve low Mach number flow problems on unstructured grids. The preconditioned compressible Euler and Navier-Stokes equations are integrated in time using a linearized implicit scheme. A newly developed fast, matrix-free implicit method, GMRES + LU?SGS, is then applied to solve the resultant system of linear equations. A variety of computations has been made for a wide range of flow conditions, for both in viscid and viscous flows, in both 2D and 3D to validate the developed method and to evaluate the effectiveness of the GMRES + LU?SGS method. The numerical results obtained indicate that the use of the GMRES + LU?SGS method leads to a significant increase in performance over the LU?SGS method, while maintaining memory requirements similar to its explicit counterpart. An overall speedup factor from one to more than two order of magnitude for all test cases in comparison with the explicit method is demonstrated.     

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.     

2D thermal/isothermal incompressible viscous flows

2D thermal and isothermal time‐dependent incompressible viscous flows are presented in rectangular domains governed by the Boussinesq approximation and Navier–Stokes equations in the stream function–vorticity formulation. The results are obtained with a simple numerical scheme based on a fixed point iterative process applied to the non‐linear elliptic systems that result after a second‐order time discretization. The iterative process leads to the solution of uncoupled, well‐conditioned, symmetric linear elliptic problems. Thermal and isothermal examples are associated with the unregularized, driven cavity problem and correspond to several aspect ratios of the cavity. Some results are presented as validation examples and others, to the best of our knowledge, are reported for the first time. The parameters involved in the numerical experiments are the Reynolds number Re, the Grashof number Gr and the aspect ratio. All the results shown correspond to steady state flows obtained from the unsteady problem. Copyright © 2005 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.     

A novel implicit numerical treatment for non‐linear turbulence models using high and low Reynolds number formulations

Non‐linear turbulence models can be seen as an improvement of the classical eddy‐viscosity concept due to their better capacity to simulate characteristics of important flows. However, application of non‐linear models demand robustness of the numerical method applied, requiring a stable discretization scheme for convergence of all variables involved. Usually, non‐linear terms are handled in an explicit manner leading to possible numerical instabilities. Thus, the present work shows the steps taken to adapt a general non‐linear constitutive equation using a new semi‐implicit numerical treatment for the non‐linear diffusion terms. The objective is to increase the degree of implicitness of the solution algorithm to enhance convergence characteristics. Flow over a backward‐facing step was computed using the control volume method applied to a boundary‐fitted coordinate system. The SIMPLE algorithm was used to relax the algebraic equations. Classical wall function and a low Reynolds number model were employed to describe the flow near the wall. The results showed that for certain combination of relaxation parameters, the semi‐implicit treatment proposed here was the sole successful treatment in order to achieve solution convergence. Also, application of the implicit method described here shows that the stability of the solution either increases (high Reynolds with non‐orthogonal mesh) or preserves the same (low Reynolds number applications). Additional advantages of the procedure proposed here lie in the possibility of testing different non‐linear expressions if one considers the enhanced robustness and stability obtained for the entire numerical algorithm. Copyright © 2010 John Wiley & Sons, Ltd.     

Viscous dissipative fluid flow past a semi‐infinite vertical plate with variable surface temperature

An analysis of the effect of viscous dissipative heat on two‐dimensional viscous incompressible fluid flow past a semi‐infinite vertical plate with variable surface temperature is carried out. The dimensionless governing equations are unsteady, two‐dimensional, coupled, and non‐linear governing equations. A most accurate, unconditionally stable and fast converging implicit finite‐difference scheme is used to solve the non‐dimensional governing equations. Velocity and temperature of the flow have been presented graphically for various parameters occurring in the problem. The local and average skin friction and Nusselt number are also shown graphically. It is observed that greater viscous dissipative heat causes a rise in the temperature. Copyright © 2007 John Wiley & Sons, Ltd.     

An implicit high-order preconditioned flux reconstruction method for low-Mach-number flow simulation with dynamic meshes

A fully implicit high-order preconditioned flux reconstruction/correction procedure via reconstruction (FR/CPR) method is developed to solve the compressible Navier-Stokes equations at low Mach numbers. A dual-time stepping approach with the second-order backward differentiation formula (BDF2) is employed to ensure temporal accuracy for unsteady flow simulation. When dynamic meshes are used to handle moving/deforming domains, the geometric conservation law is implicitly enforced to eliminate errors due to the resolution discrepancy between BDF2 and the spatial FR/CPR discretization. The large linear system resulted from the spatial and temporal discretizations is tackled with the restarted generalized minimal residual solver in the PETSc (portable, extensible toolkit for scientific computation) library. Through several benchmark steady and unsteady numerical tests, the preconditioned FR/CPR methods have demonstrated good convergence and accuracy for simulating flows at low Mach numbers. The new flow solver is then used to study the effects of Mach number on unsteady force generation over a plunging airfoil when operating in low-Mach-number flows. It is observed that weak compressibility has a significant impact on thrust generation but has a negligible effect on lift generation of an oscillating airfoil.