Multiple-correction hybrid k-exact schemes for high-order compressible RANS-LES simulations on fully unstructured grids

Abstract A Godunov's type unstructured finite volume method suitable for highly compressible turbulent scale-resolving simulations around complex geometries is constructed by using a successive correction technique. First, a family of k-exact Godunov schemes is developed by recursively correcting the truncation error of the piecewise polynomial representation of the primitive variables. The keystone of the proposed approach is a quasi-Green gradient operator which ensures consistency on general meshes. In addition, a high-order single-point quadrature formula, based on high-order approximations of the successive derivatives of the solution, is developed for flux integration along cell faces. The proposed family of schemes is compact in the algorithmic sense, since it only involves communications between direct neighbors of the mesh cells. The numerical properties of the schemes up to fifth-order are investigated, with focus on their resolvability in terms of number of mesh points required to resolve a given wavelength accurately. Afterwards, in the aim of achieving the best possible trade-off between accuracy, computational cost and robustness in view of industrial flow computations, we focus more specifically on the third-order accurate scheme of the family, and modify locally its numerical flux in order to reduce the amount of numerical dissipation in vortex-dominated regions. This is achieved by switching from the upwind scheme, mostly applied in highly compressible regions, to a fourth-order centered one in vortex-dominated regions. An analytical switch function based on the local grid Reynolds number is adopted in order to warrant numerical stability of the recentering process. Numerical applications demonstrate the accuracy and robustness of the proposed methodology for compressible scale-resolving computations. In particular, supersonic RANS/LES computations of the flow over a cavity are presented to show the capability of the scheme to predict flows with shocks, vortical structures and complex geometries.

[1]  Takeo Nakagawa,et al.  Vortex shedding behind a square cylinder in transonic flows , 1987, Journal of Fluid Mechanics.

[2]  Zhenbing Luo,et al.  Flow Visualization of Supersonic Flow over a Finite Cylinder , 2012 .

[3]  Abrar H. Mohammad.  Large Eddy Simulation of Flow over a Cylinder Using High-Order Spectral Difference Method , 2010 .

[4]  H. J.,et al.  Hydrodynamics , 1924, Nature.

[5]  Steven A. Orszag,et al.  The Taylor-Green vortex and fully developed turbulence , 1984 .

[6]  Jonathan Dudley,et al.  A high-order CFD method using successive differentiation , 2015, J. Comput. Phys..

[7]  A. Prasad,et al.  The instability of the shear layer separating from a bluff body , 1997, Journal of Fluid Mechanics.

[8]  Friedrich Ringleb,et al.  Exakte Lösungen der Differentialgleichungen einer adiabatischen Gasströmung , 1940 .

[9]  Chunlei Liang,et al.  Large Eddy Simulation of Flow over a Cylinder Using High-Order Spectral Difference Method , 2008 .

[10]  Antony Jameson,et al.  On the Development of a High-Order, Multi-GPU Enabled, Compressible Viscous Flow Solver for Mixed Unstructured Grids , 2011 .

[11]  P. Spalart.  Strategies for turbulence modelling and simulations , 2000 .

[12]  Grégoire Pont,et al.  Self adaptive turbulence models for unsteady compressible flows Modèles de turbulence auto-adaptatifs pour la simulation des écoulements compressibles instationnaires , 2015 .

[13]  Paola Cinnella,et al.  Automatic Hybrid RANS/LES Strategy for Industrial CFD , 2015 .

[14]  P. Cinnella,et al.  Investigation of High-Order Methods in Large-Eddy Simulation of Separated Flow in a Channel with Periodic Constrictions , 2018 .

[15]  P. Spalart.  Comments on the feasibility of LES for wings, and on a hybrid RANS/LES approach , 1997 .

[16]  Florian R. Menter,et al.  The Scale-Adaptive Simulation Method for Unsteady Turbulent Flow Predictions. Part 1: Theory and Model Description , 2010 .

[17]  Florian Haider,et al.  Discrétisation en maillage non structuré général et applications LES , 2009 .

[18]  Laslo T. Diosady,et al.  Case 3.3: Taylor-Green Vortex Evolution , 2015 .

[19]  Alain Lerat,et al.  Spectral properties of high-order residual-based compact schemes for unsteady compressible flows , 2013, J. Comput. Phys..

[20]  Clinton P. T. Groth,et al.  High-order solution-adaptive central essentially non-oscillatory (CENO) method for viscous flows , 2011, J. Comput. Phys..

[21]  G. Aubard,et al.  Comparison of Subgrid-scale Viscosity Models and Selective Filtering Strategy for Large-eddy Simulations , 2013 .

[23]  Paola Cinnella,et al.  Assessment of automatic Hybrid RANS/LES Models for Industrial CFD , 2014 .

[24]  Pierre Sagaut,et al.  Multiscale And Multiresolution Approaches In Turbulence , 2006 .

[25]  Pierre Brenner.  Three-dimensional aerodynamics with moving bodies applied to solid propellant , 1991 .

[26]  P. Spalart,et al.  A New Version of Detached-eddy Simulation, Resistant to Ambiguous Grid Densities , 2006 .

[27]  Nikolaus A. Adams,et al.  An adaptive local deconvolution method for implicit LES , 2005, J. Comput. Phys..

[28]  Sébastien Deck,et al.  Zonal Detached Eddy Simulation of a Controlled Propulsive Jet , 2007 .

[29]  Z. Wang.  High-order methods for the Euler and Navier–Stokes equations on unstructured grids , 2007 .

[30]  B. Ghosh.  Stresses in a long hollow cylindrically aeolotropic circular cylinder with plane ends subjected to surface forces , 1972 .

[31]  P. Spalart,et al.  Physical and Numerical Upgrades in the Detached-Eddy Simulation of Complex Turbulent Flows , 2002 .

[32]  Richard A. Gentry,et al.  An Eulerian differencing method for unsteady compressible flow problems , 1966 .

[33]  Rémi Abgrall,et al.  High‐order CFD methods: current status and perspective , 2013 .

[34]  Jean-Pierre Croisille,et al.  Efficient Implementation of High Order Reconstruction in Finite Volume Methods , 2011 .

[35]  Erik Dick,et al.  Construction of explicit and implicit dynamic finite difference schemes and application to the large-eddy simulation of the Taylor-Green vortex , 2009, J. Comput. Phys..

[36]  S. Girimaji.  Partially-Averaged Navier-Stokes Model for Turbulence: A Reynolds-Averaged Navier-Stokes to Direct Numerical Simulation Bridging Method , 2006 .

[37]  Bruno Maugars.  Méthodes de volumes finis d'ordre élevé en maillages non coïncidents pour les écoulements dans les turbomachines , 2016 .

[38]  C. Ollivier-Gooch,et al.  A high-order-accurate unstructured mesh finite-volume scheme for the advection-diffusion equation , 2002 .

[39]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

[40]  Michael Dumbser,et al.  Quadrature-free non-oscillatory finite volume schemes on unstructured meshes for nonlinear hyperbolic systems , 2007, J. Comput. Phys..

[41]  焦予秦,et al.  Detached—Eddy Simulation方法模拟不同类型翼型的失速特性 , 2005 .

[42]  Sangsan Lee,et al.  Large eddy simulation of shock turbulence interaction , 1993 .

[43]  Miguel R. Visbal,et al.  High-Order-Accurate Methods for Complex Unsteady Subsonic Flows , 1999 .

[44]  P. Frederickson,et al.  Higher order solution of the Euler equations on unstructured grids using quadratic reconstruction , 1990 .

[45]  Neil D. Sandham,et al.  Large Eddy Simulation of Flow Around an Airfoil Near Stall , 2010 .

[46]  C. Tam,et al.  Dispersion-relation-preserving finite difference schemes for computational acoustics , 1993 .

[47]  E. Koschmieder.  Taylor vortices between eccentric cylinders , 1976 .

[48]  Doru Caraeni,et al.  Efficient 3rd-Order Finite-Volume Discretization Using Iterative Quadratic Data Reconstruction on Unstructured Grids , 2009 .

[49]  Lin He,et al.  Characteristics of the cylinder-induced shock wave and turbulent boundary layer interactions , 2016, J. Vis..

[50]  Chi-Wang Shu.  Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws , 1998 .

[51]  Pierre Brenner,et al.  Hybrid RANS/LES Simulation of a Space Launcher Using a High Order Finite Volume Scheme and Grid Intersections Technique , 2016 .

[52]  Eiji Shima,et al.  Green–Gauss/Weighted-Least-Squares Hybrid Gradient Reconstruction for Arbitrary Polyhedra Unstructured Grids , 2013 .

[53]  J. Blair Perot,et al.  A self-adapting turbulence model for flow simulation at any mesh resolution , 2007 .

[54]  Lawrence L. Takacs,et al.  A Two-Step Scheme for the Advection Equation with Minimized Dissipation and Dispersion Errors , 1985 .

[55]  Timothy J. Barth,et al.  Recent developments in high order K-exact reconstruction on unstructured meshes , 1993 .

[56]  Stephane Catris,et al.  Density corrections for turbulence models , 2000 .

[57]  D. Drikakis.  Advances in turbulent flow computations using high-resolution methods , 2003 .

[58]  E. Dick,et al.  On the performance of relaxation filtering for large-eddy simulation , 2013 .

[59]  Dimitri J. Mavriplis,et al.  Revisiting the Least-squares Procedure for Gradient Reconstruction on Unstructured Meshes , 2003 .

[60]  Paola Cinnella,et al.  High-Order Hybrid RANS/LES Strategy for Industrial Applications , 2018 .

[61]  Sergio Pirozzoli,et al.  On the spectral properties of shock-capturing schemes , 2006, J. Comput. Phys..

[62]  P. Sagaut,et al.  On the Use of Shock-Capturing Schemes for Large-Eddy Simulation , 1999 .

[63]  Herman Deconinck,et al.  Higher order upwind finite volume schemes with ENO-properties for general unstructured meshes , 1992 .