On the performance of WENO/TENO schemes to resolve turbulence in DNS/LES of high‐speed compressible flows

High‐speed compressible turbulent flows typically contain discontinuities and have been widely modeled using Weighted Essentially Non‐Oscillatory (WENO) schemes due to their high‐order accuracy and sharp shock capturing capability. However, such schemes may damp the small scales of turbulence and result in inaccurate solutions in the context of turbulence‐resolving simulations. In this connection, the recently developed Targeted Essentially Non‐Oscillatory (TENO) schemes, including adaptive variants, may offer significant improvements. The present study aims to quantify the potential of these new schemes for a fully turbulent supersonic flow. Specifically, DNS of a compressible turbulent channel flow with M = 1.5 and Re_τ = 222 is conducted using OpenSBLI, a high‐order finite difference computational fluid dynamics framework. This flow configuration is chosen to decouple the effect of flow discontinuities and turbulence and focus on the capability of the aforementioned high‐order schemes to resolve turbulent structures. The effect of the spatial resolution in different directions and coarse grid implicit LES are also evaluated against the WALE LES model. The TENO schemes are found to exhibit significant performance improvements over the WENO schemes in terms of the accuracy of the statistics and the resolution of the three‐dimensional vortical structures. The sixth‐order adaptive TENO scheme is found to produce comparable results to those obtained with nondissipative fourth‐ and sixth‐order central schemes and reference data obtained with spectral methods. Although the most computationally expensive scheme, it is shown that this adaptive scheme can produce satisfactory results if used as an implicit LES model.     

Improved third‐order weighted essentially nonoscillatory schemes with new smoothness indicators

In this article, we present two improved third‐order weighted essentially nonoscillatory (WENO) schemes for recovering their design‐order near first‐order critical points. The schemes are constructed in the framework of third‐order WENO‐Z scheme. Two new global smoothness indicators, τ_L3 and τ_L4, are devised by a nonlinear combination of local smoothness indicators (IS_k) and reference values (IS_G) based on Lagrangian interpolation polynomial. The performances of the proposed schemes are evaluated on several numerical tests governed by one‐dimensional linear advection equation or one‐ and two‐dimensional Euler equations. Numerical results indicate that the presented schemes provide less dissipation and higher resolution than the original WENO3‐JS and subsequent WENO3‐N scheme.     

Application of central differencing and low‐dissipation weights in a weighted compact nonlinear scheme

This paper proposes WCNS‐CU‐Z, a weighted compact nonlinear scheme, that incorporates adapted central difference and low‐dissipative weights together with concepts of the adaptive central‐upwind sixth‐order weighted essentially non‐oscillatory scheme (WENO‐CU) and WENO‐Z schemes. The newly developed WCNS‐CU‐Z is a high‐resolution scheme, because interpolation of this scheme employs a central stencil constructed by upwind and downwind stencils. The smoothness indicator of the downwind stencil is calculated using the entire central stencil, and the downwind stencil is stopped around the discontinuity for stability. Moreover, interpolation of the sixth‐order WCNS‐CU‐Z exhibits sufficient accuracy in the smooth region through use of low‐dissipative weights. The sixth‐order WCNS‐CU‐Zs are implemented with a robust linear difference formulation (R‐WCNS‐CU6‐Z), and the resolution and robustness of this scheme were evaluated. These evaluations showed that R‐WCNS‐CU6‐Z is capable of achieving a higher resolution than the seventh‐order classical robust weighted compact nonlinear scheme and can provide a crisp result in terms of discontinuity. Among the schemes tested, R‐WCNS‐CU6‐Z has been shown to be robust, and variable interpolation type R‐WCNS‐CU6‐Z (R‐WCNS‐CU6‐Z‐V) provides a stable computation by modifying the first‐order interpolation when negative density or negative pressure arises after nonlinear interpolation. Copyright © 2016 John Wiley & Sons, Ltd.     

Convergence issues in using high‐resolution schemes and lower–upper symmetric Gauss–Seidel method for steady shock‐induced combustion problems

This paper reports numerical convergence study for simulations of steady shock‐induced combustion problems with high‐resolution shock‐capturing schemes. Five typical schemes are used: the Roe flux‐based monotone upstream‐centered scheme for conservation laws (MUSCL) and weighted essentially non‐oscillatory (WENO) schemes, the Lax–Friedrichs splitting‐based non‐oscillatory no‐free parameter dissipative (NND) and WENO schemes, and the Harten–Yee upwind total variation diminishing (TVD) scheme. These schemes are implemented with the finite volume discretization on structured quadrilateral meshes in dimension‐by‐dimension way and the lower–upper symmetric Gauss–Seidel (LU–SGS) relaxation method for solving the axisymmetric multispecies reactive Navier–Stokes equations. Comparison of iterative convergence between different schemes has been made using supersonic combustion flows around a spherical projectile with Mach numbers M = 3.55 and 6.46 and a ram accelerator with M = 6.7. These test cases were regarded as steady combustion problems in literature. Calculations on gradually refined meshes show that the second‐order NND, MUSCL, and TVD schemes can converge well to steady states from coarse through fine meshes for M = 3.55 case in which shock and combustion fronts are separate, whereas the (nominally) fifth‐order WENO schemes can only converge to some residual level. More interestingly, the numerical results show that all the schemes do not converge to steady‐state solutions for M = 6.46 in the spherical projectile and M = 6.7 in the ram accelerator cases on fine meshes although they all converge on coarser meshes or on fine meshes without chemical reactions. The result is based on the particular preconditioner of LU–SGS scheme. Possible reasons for the nonconvergence in reactive flow simulation are discussed.Copyright © 2012 John Wiley & Sons, Ltd.     

A new well‐balanced Hermite weighted essentially non‐oscillatory scheme for shallow water equations

Hermite weighted essentially non‐oscillatory (HWENO) methods were introduced in the literature, in the context of Euler equations for gas dynamics, to obtain high‐order accuracy schemes characterized by high compactness (e.g. Qiu and Shu, J. Comput. Phys. 2003; 193 :115). For example, classical fifth‐order weighted essentially non‐oscillatory (WENO) reconstructions are based on a five‐cell stencil whereas the corresponding HWENO reconstructions are based on a narrower three‐cell stencil. The compactness of the schemes allows easier treatment of the boundary conditions and of the internal interfaces. To obtain this compactness in HWENO schemes both the conservative variables and their first derivatives are evolved in time, whereas in the original WENO schemes only the conservative variables are evolved. In this work, an HWENO method is applied for the first time to the shallow water equations (SWEs), including the source term due to the bottom slope, to obtain a fourth‐order accurate well‐balanced compact scheme. Time integration is performed by a strong stability preserving the Runge–Kutta method, which is a five‐step and fourth‐order accurate method. Besides the classical SWE, the non‐homogeneous equations describing the time and space evolution of the conservative variable derivatives are considered here. An original, well‐balanced treatment of the source term involved in such equations is developed and tested. Several standard one‐dimensional test cases are used to verify the high‐order accuracy, the C‐property and the good resolution properties of the model. Copyright © 2010 John Wiley & Sons, Ltd.     

A high‐resolution scheme for compressible multicomponent flows with shock waves

A simple methodology for a high‐resolution scheme to be applied to compressible multicomponent flows with shock waves is investigated. The method is intended for use with direct numerical simulation or large eddy simulation of compressible multicomponent flows. The method dynamically adds non‐linear artificial diffusivity locally in space to capture different types of discontinuities such as a shock wave, contact surface or material interface while a high‐order compact differencing scheme resolves a broad range of scales in flows. The method is successfully applied to several one‐dimensional and two‐dimensional compressible multicomponent flow problems with shock waves. The results are in good agreement with experiments and earlier computations qualitatively and quantitatively. The method captures unsteady shock and material discontinuities without significant spurious oscillations if initial start‐up errors are properly avoided. Comparisons between the present numerical scheme and high‐order weighted essentially non‐oscillatory (WENO) schemes illustrate the advantage of the present method for resolving a broad range of scales of turbulence while capturing shock waves and material interfaces. Also the present method is expected to require less computational cost than popular high‐order upwind‐biased schemes such as WENO schemes. The mass conservation for each species is satisfied due to the strong conservation form of governing equations employed in the method. Copyright © 2010 John Wiley & Sons, Ltd.     

High‐order ADER‐WENO ALE schemes on unstructured triangular meshes—application of several node solvers to hydrodynamics and magnetohydrodynamics

In this paper, we present a class of high‐order accurate cell‐centered arbitrary Lagrangian–Eulerian (ALE) one‐step ADER weighted essentially non‐oscillatory (WENO) finite volume schemes for the solution of nonlinear hyperbolic conservation laws on two‐dimensional unstructured triangular meshes. High order of accuracy in space is achieved by a WENO reconstruction algorithm, while a local space–time Galerkin predictor allows the schemes to be high order accurate also in time by using an element‐local weak formulation of the governing PDE on moving meshes. The mesh motion can be computed by choosing among three different node solvers, which are for the first time compared with each other in this article: the node velocity may be obtained either (i) as an arithmetic average among the states surrounding the node, as suggested by Cheng and Shu, or (ii) as a solution of multiple one‐dimensional half‐Riemann problems around a vertex, as suggested by Maire, or (iii) by solving approximately a multidimensional Riemann problem around each vertex of the mesh using the genuinely multidimensional Harten–Lax–van Leer Riemann solver recently proposed by Balsara et al. Once the vertex velocity and thus the new node location have been determined by the node solver, the local mesh motion is then constructed by straight edges connecting the vertex positions at the old time level tⁿ with the new ones at the next time level t^{n + 1}. If necessary, a rezoning step can be introduced here to overcome mesh tangling or highly deformed elements. The final ALE finite volume scheme is based directly on a space–time conservation formulation of the governing PDE system, which therefore makes an additional remapping stage unnecessary, as the ALE fluxes already properly take into account the rezoned geometry. In this sense, our scheme falls into the category of direct ALE methods. Furthermore, the geometric conservation law is satisfied by the scheme by construction. We apply the high‐order algorithm presented in this paper to the Euler equations of compressible gas dynamics as well as to the ideal classical and relativistic magnetohydrodynamic equations. We show numerical convergence results up to fifth order of accuracy in space and time together with some classical numerical test problems for each hyperbolic system under consideration. Copyright © 2014 John Wiley & Sons, Ltd.     

High‐order ENO and WENO schemes for unstructured grids

This work describes the implementation and analysis of high‐order accurate schemes applied to high‐speed flows on unstructured grids. The class of essentially non‐oscillatory schemes (ENO), that includes weighted ENO schemes (WENO), is discussed in the paper with regard to the implementation of third‐ and fourth‐order accurate methods. The entire reconstruction process of ENO and WENO schemes is described with emphasis on the stencil selection algorithms. The stencils can be composed by control volumes with any number of edges, e.g. triangles, quadrilaterals and hybrid meshes. In the paper, ENO and WENO schemes are implemented for the solution of the dimensionless, 2‐D Euler equations in a cell centred finite volume context. High‐order flux integration is achieved using Gaussian quadratures. An approximate Riemann solver is used to evaluate the fluxes on the interfaces of the control volumes and a TVD Runge–Kutta scheme provides the time integration of the equations. Such a coupling of all these numerical tools, together with the high‐order interpolation of primitive variables provided by ENO and WENO schemes, leads to the desired order of accuracy expected in the solutions. An adaptive mesh refinement technique provides better resolution in regions with strong flowfield gradients. Results for high‐speed flow simulations are presented with the objective of assessing the implemented capability. Copyright © 2007 John Wiley & Sons, Ltd.     

DNS of backward‐facing step flow with fully turbulent inflow

Direct numerical simulation (DNS) has been performed to study the channel flow over a backward‐facing step at a Reynolds number Re_b=5600 based on the step height h and the inflow bulk velocity U_b. A dynamic method has been used in order to generate realistic turbulent inflow conditions. The results upstream of the step compared well with the fully developed channel flow. Downstream of the step our results show excellent agreement with experimental data. Copyright © 2009 John Wiley & Sons, Ltd.     

Reynolds‐stress modelling of M = 2.25 shock‐wave/turbulent boundary‐layer interaction

M = 2.25 shock‐wave/turbulent‐boundary‐layer interactions over a compression ramp for several angles (8, 13 and 18°) at Reynolds‐number Re=7 × 10³ were simulated with three low‐Reynolds second‐moment closures and a linear low‐Reynolds standard k–ε model. A detailed assessment of the turbulence closures by comparison with both mean‐flow and turbulent experimental quantities is presented. The Reynolds‐stress model which is wall‐topology free and which uses an optimized redistribution closure, is in good agreement with experimental data both for wall‐pressure and mean‐velocity profiles. Detailed analysis of three components of the Reynolds‐stress tensor (comparison with measurements and transport‐equation budgets) provides a critical evaluation of full Reynolds‐stress models for the separated supersonic compression ramp. The discrepancy observed in the shock‐wave foot region, between computations and measurements for the Reynolds‐stresses profiles, could be explained by considering the experimental shock‐wave oscillation and directions for future modelling work are indicated. Copyright © 2007 John Wiley & Sons, Ltd.