A gas‐kinetic scheme for the simulation of turbulent flows on unstructured meshes

This paper presents a numerical method for simulating turbulent flows via coupling the Boltzmann BGK equation with Spalart–Allmaras one equation turbulence model. Both the Boltzmann BGK equation and the turbulence model equation are carried out using the finite volume method on unstructured meshes, which is different from previous works on structured grid. The application of the gas‐kinetic scheme is extended to the simulation of turbulent flows with arbitrary geometries. The adaptive mesh refinement technique is also adopted to reduce the computational cost and improve the efficiency of meshes. To organize the unstructured mesh data structure efficiently, a non‐manifold hybrid mesh data structure is extended for polygonal cells. Numerical experiments are performed on incompressible flow over a smooth flat plate and compressible turbulent flows around a NACA 0012 airfoil using unstructured hybrid meshes. These numerical results are found to be in good agreement with experimental data and/or other numerical solutions, demonstrating the applicability of the proposed method to simulate both subsonic and transonic turbulent flows. Copyright © 2016 John Wiley & Sons, Ltd.     

Simulation of multiple shock–shock interference using implicit anti‐diffusive WENO schemes

Accurate computations of two‐dimensional turbulent hypersonic shock–shock interactions that arise when single and dual shocks impinge on the bow shock in front of a cylinder are presented. The simulation methods used are a class of lower–upper symmetric‐Gauss–Seidel implicit anti‐diffusive weighted essentially non‐oscillatory (WENO) schemes for solving the compressible Navier–Stokes equations with Spalart–Allmaras one‐equation turbulence model. A numerical flux of WENO scheme with anti‐diffusive flux correction is adopted, which consists of first‐order and high‐order fluxes and allows for a more flexible choice of first‐order dissipative methods. Experimental flow fields of type IV shock–shock interactions with single and dual incident shocks by Wieting are computed. By using the WENO scheme with anti‐diffusive flux corrections, the present solution indicates that good accuracy is maintained and contact discontinuities are sharpened markedly as compared with the original WENO schemes on the same meshes. Computed surface pressure distribution and heat transfer rate are also compared with experimental data and other computational results and good agreement is found. Copyright © 2009 John Wiley & Sons, Ltd.     

Applications of a higher‐order bounded numerical scheme to turbulent flows

This paper applies the higher‐order bounded numerical scheme Weighted Average Coefficients Ensuring Boundedness (WACEB) to simulate two‐ and three‐dimensional turbulent flows. In the scheme, a weighted average formulation is used for interpolating the variables at cell faces and the weighted average coefficients are determined from a normalized variable formulation and total variation diminishing (TVD) constraints to ensure the boundedness of the solution. The scheme is applied to two turbulent flow problems: (1) two‐dimensional turbulent flow around a blunt plate; and (2) three‐dimensional turbulent flow inside a mildly curved U‐bend. In the present study, turbulence is evaluated by using a low‐Reynolds number version of the k–ω model. For the flow simulation, the QUICK scheme is applied to the momentum equations while either the WACEB scheme (Method 1) or the UPWIND scheme (Method 2) is used for the turbulence equations. The present study shows that the WACEB scheme has at least second‐order accuracy while ensuring boundedness of the solutions. The present numerical study for a pure convection problem shows that the ‘TVD’ slope ranges from 2 to 4. For the turbulent recirculating flow, two different mixed procedures (Method 1 and Method 2) produce a substantial difference for the mean velocities as well as for the turbulence kinetic energy. Method 1 predicts better results than Method 2 does, comparing the analytical solution and the experimental data. For the turbulent flow inside the mildly curved U‐bend, although the predictions of velocity distributions with two procedures are very close, a noticeable difference of turbulence kinetic energy is exhibited. It is noticed that the discrepancy exists between numerical results and the experimental data. The reason is the limit of the two‐equation turbulence model to such complex turbulent flows with extra strain‐rates. Copyright © 2001 John Wiley & Sons, Ltd.     

Computation of Compressible Flows using a Two-equation Turbulence Model on Unstructured Grids

This paper presents a numerical method for solving compressible turbulent flows using a k - l turbulence model on unstructured meshes. The flow equations and turbulence equations are solved in a loosely coupled manner. The flow equations are advanced in time using a multi-stage Runge-Kutta time stepping scheme, while the turbulence equations are advanced using a multi-stage point-implicit scheme. The positivity of turbulence variables is achieved using a simple change of dependent variables. The developed method is used to compute a variety of turbulent flow problems. The results obtained are in good agreement with theoretical and experimental data, indicating that the present method provides a viable and robust algorithm for computing turbulent flows on unstructured meshes.     

A high‐order finite difference method for numerical simulations of supersonic turbulent flows

This paper describes a new class of three‐dimensional finite difference schemes for high‐speed turbulent flows in complex geometries based on the high‐order monotonicity‐preserving (MP) method. Simulations conducted for various 1D, 2D, and 3D problems indicate that the new high‐order MP schemes can preserve sharp changes in the flow variables without spurious oscillations and are able to capture the turbulence at the smallest computed scales. Our results also indicate that the MP method has less numerical dissipation and faster grid convergence than the weighted essentially non‐oscillatory method. However, both of these methods are computationally more demanding than the COMP method and are only used for the inviscid fluxes. To reduce the computational cost for reacting flows, the scalar equations are solved by the COMP method, which is shown to yield similar results to those obtained by the MP in supersonic turbulent flows with strong shock waves. Copyright © 2011 John Wiley & Sons, Ltd.     

Calculation of three‐dimensional turbulent flow with a finite volume multigrid method

A numerical method for the efficient calculation of three‐dimensional incompressible turbulent flow in curvilinear co‐ordinates is presented. The mathematical model consists of the Reynolds averaged Navier–Stokes equations and the k–ε turbulence model. The numerical method is based on the SIMPLE pressure‐correction algorithm with finite volume discretization in curvilinear co‐ordinates. To accelerate the convergence of the solution method a full approximation scheme‐full multigrid (FAS‐FMG) method is utilized. The solution of the k–ε transport equations is embedded in the multigrid iteration. The improved convergence characteristic of the multigrid method is demonstrated by means of several calculations of three‐dimensional flow cases. Copyright © 1999 John Wiley & Sons, Ltd.     

Evaluation of one‐ and two‐equation low‐Re turbulence models. Part II—Vortex‐generator jet and diffusing S‐duct flows

This second segment of the two‐part paper systematically examines several turbulence models in the context of two flows, namely a vortex flow created by an inclined jet in crossflow, and the flow field in a diffusing S‐shaped duct. The test cases are chosen on the basis of availability of high‐quality and detailed experimental data. The tested turbulence models are integrated to solid surfaces and consist of: Rodi's two‐layer k–ε model, Wilcox's k–ω model, Menter's two‐equation shear–stress‐transport model, and the one‐equation model of Spalart and Allmaras. The objective of the study is to establish the prediction accuracy of these turbulence models with respect to three‐dimensional separated flows with streamline curvature. At the same time, the study establishes the minimum spatial resolution requirements for each of these turbulence closures, and identifies the proper low‐Mach‐number preconditioning and artificial diffusion settings of a Reynolds‐averaged Navier–Stokes algorithm for optimum rate of convergence and minimum adverse impact on prediction accuracy. Copyright © 2003 John Wiley & Sons, Ltd.     

A stress‐based least‐squares finite‐element model for incompressible Navier–Stokes equations

In this paper we present a stress‐based least‐squares finite‐element formulation for the solution of the Navier–Stokes equations governing flows of viscous incompressible fluids. Stress components are introduced as independent variables to make the system first order. Continuity equation becomes an algebraic equation and is eliminated from the system with suitable modifications. The h and p convergence are verified using the exact solution of Kovasznay flow. Steady flow past a large circular cylinder in a channel is solved to test mass conservation. Transient flow over a backward‐facing step problem is solved on several meshes. Results are compared with that obtained using vorticity‐based first‐order formulation for both benchmark problems. Copyright © 2007 John Wiley & Sons, Ltd.     

Implicit meanflow–multigrid algorithms for Reynolds stress model computation of 3‐D anisotropy‐driven and compressible flows

The present paper investigates the multigrid (MG) acceleration of compressible Reynolds‐averaged Navier–Stokes computations using Reynolds‐stress model 7‐equation turbulence closures, as well as lower‐level 2‐equation models. The basic single‐grid SG algorithm combines upwind‐biased discretization with a subiterative local‐dual‐time‐stepping time‐integration procedure. MG acceleration, using characteristic MG restriction and prolongation operators, is applied on meanflow variables only (MF–MG), turbulence variables being simply injected onto coarser grids. A previously developed non‐time‐consistent (for steady flows) full‐approximation‐multigrid (s–MG) is assessed for 3‐D anisotropy‐driven and/or separated flows, which are dominated by the convergence of turbulence variables. Even for these difficult test cases CPU‐speed‐ups r_CPUSUP∈[3, 5] are obtained. Alternative, potentially time‐consistent approaches (unsteady u–MG), where MG acceleration is applied at each subiteration, are also examined, using different subiterative strategies, MG cycles, and turbulence models. For 2‐D shock wave/turbulent boundary layer interaction, the fastest s–MG approach, with a V(2, 0) sawtooth cycle, systematically yields CPU‐speed‐ups of 5±½, quasi‐independent of the particular turbulence closure used. Copyright © 2008 John Wiley & Sons, Ltd.     

Large eddy simulations of wall‐bounded flows using a simplified immersed boundary method and high‐order compact schemes

This paper presents a solution algorithm based on an immersed boundary (IB) method that can be easily implemented in high‐order codes for incompressible flows. The time integration is performed using a predictor‐corrector approach, and the projection method is used for pressure‐velocity coupling. Spatial discretization is based on compact difference schemes and is performed on half‐staggered meshes. A basic algorithm for body‐fitted meshes using the aforementioned solution method was developed by A. Tyliszczak (see article “A high‐order compact difference algorithm for half‐staggered grids for laminar and turbulent incompressible flows” in Journal of Computational Physics) and proved to be very accurate. In this paper, the formulated algorithm is adapted for use with the IB method in the framework of large eddy simulations. The IB method is implemented using its simplified variant without the interpolation (stepwise approach). The computations are performed for a laminar flow around a 2D cylinder, a turbulent flow in a channel with a wavy wall, and around a sphere. Comparisons with literature data confirm that the proposed method can be successfully applied for complex flow problems. The results are verified using the classical approach with body‐fitted meshes and show very good agreement both in laminar and turbulent regimes. The mean (velocity and turbulent kinetic energy profiles and drag coefficients) and time‐dependent (Strouhal number based on the drag coefficient) quantities are analyzed, and they agree well with reference solutions. Two subfilter models are compared, ie, the model of Vreman (see article “An eddy‐viscosity subgrid‐scale model for turbulent shear flow: algebraic theory and applications” in Physics and Fluids) and σ model (Nicoud et al, see article “Using singular values to build a subgrid‐scale model for large eddy simulations” in Physics and Fluids). The tests did not reveal evident advantages of any of these models, and from the point of view of solution accuracy, the quality of the computational meshes turned out to be much more important than the subfilter modeling.     

On coupling the Reynolds‐averaged Navier–Stokes equations with two‐equation turbulence model equations

Two methods for coupling the Reynolds‐averaged Navier–Stokes equations with the q–ω turbulence model equations on structured grid systems have been studied; namely a loosely coupled method and a strongly coupled method. The loosely coupled method first solves the Navier–Stokes equations with the turbulent viscosity fixed. In a subsequent step, the turbulence model equations are solved with all flow quantities fixed. On the other hand, the strongly coupled method solves the Reynolds‐averaged Navier–Stokes equations and the turbulence model equations simultaneously. In this paper, numerical stabilities of both methods in conjunction with the approximated factorization‐alternative direction implicit method are analysed. The effect of the turbulent kinetic energy terms in the governing equations on the convergence characteristics is also studied. The performance of the two methods is compared for several two‐ and three‐dimensional problems. Copyright © 2005 John Wiley & Sons, Ltd.     

Evaluation of one‐ and two‐equation low‐Re turbulence models. Part I—Axisymmetric separating and swirling flows

This first segment of the two‐part paper systematically examines several turbulence models in the context of three flows, namely a simple flat‐plate turbulent boundary layer, an axisymmetric separating flow, and a swirling flow. The test cases are chosen on the basis of availability of high‐quality and detailed experimental data. The tested turbulence models are integrated to solid surfaces and consist of: Rodi's two‐layer k–ε model, Chien's low‐Reynolds number k–ε model, Wilcox's k–ω model, Menter's two‐equation shear‐stress‐transport model, and the one‐equation model of Spalart and Allmaras. The objective of the study is to establish the prediction accuracy of these turbulence models with respect to axisymmetric separating flows, and flows of high streamline curvature. At the same time, the study establishes the minimum spatial resolution requirements for each of these turbulence closures, and identifies the proper low‐Mach‐number preconditioning and artificial diffusion settings of a Reynolds‐averaged Navier–Stokes algorithm for optimum rate of convergence and minimum adverse impact on prediction accuracy. Copyright © 2003 John Wiley & Sons, Ltd.     

A marker‐and‐cell approach to free surface 2‐D multiphase flows

This work describes a methodology to simulate free surface incompressible multiphase flows. This novel methodology allows the simulation of multiphase flows with an arbitrary number of phases, each of them having different densities and viscosities. Surface and interfacial tension effects are also included. The numerical technique is based on the GENSMAC front‐tracking method. The velocity field is computed using a finite‐difference discretization of a modification of the Navier–Stokes equations. These equations together with the continuity equation are solved for the two‐dimensional multiphase flows, with different densities and viscosities in the different phases. The governing equations are solved on a regular Eulerian grid, and a Lagrangian mesh is employed to track free surfaces and interfaces. The method is validated by comparing numerical with analytic results for a number of simple problems; it was also employed to simulate complex problems for which no analytic solutions are available. The method presented in this paper has been shown to be robust and computationally efficient. Copyright © 2012 John Wiley & Sons, Ltd.     

Direct simulation of two‐dimensional turbulent flow over a surface‐mounted obstacle

Two‐dimensional turbulent flow over a surface‐mounted obstacle is studied as a numerical experiment that takes place in a wind tunnel. The transient Navier–Stokes equations are solved directly with Galerkin finite elements. The Reynolds number defined with respect to the height of the wind tunnel is 12 518. Instantaneous streamline patterns are shown that give a complete picture of the flow phenomena. Energy and enstrophy spectra yield the dual cascade of two‐dimensional turbulence and the ?1 power law decay of enstrophy. Mean values of velocities and root mean square fluctuations are compared with the available experimental results. Other statistical characteristics of turbulence such as Eulerian autocorrelation coefficients, longitudinal and lateral coefficients are also computed. Finally, oscillation diagrams of computed velocity fluctuations yield the chaotic behaviour of turbulence. Copyright © 2007 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.     

Momentum/continuity coupling with large non‐isotropic momentum source terms

Pressure‐based methods such as the SIMPLE algorithm are frequently used to determine a coupled solution between the component momentum equations and the continuity equation. This paper presents a colocated variable pressure correction algorithm for control volumes of polyhedral/polygonal cell topologies. The correction method is presented independent of spatial approximation. The presence of non‐isotropic momentum source terms is included in the proposed algorithm to ensure its applicability to multi‐physics applications such as gas and particulate flows. Two classic validation test cases are included along with a newly proposed test case specific to multiphase flows. The classic validation test cases demonstrate the application of the proposed algorithm on truly arbitrary polygonal/polyhedral cell meshes. A comparison between the current algorithm and commercially available software is made to demonstrate that the proposed algorithm is competitively efficient. The newly proposed test case demonstrates the benefits of the current algorithm when applied to a multiphase flow situation. The numerical results from this case show that the proposed algorithm is more robust than other methods previously proposed. Copyright © 2009 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.