An assessment of atypical mesh topologies for low-Mach large-eddy simulation

Abstract An implicit, low-dissipation, low-Mach, variable density control volume finite element formulation is used to explore foundational understanding of numerical accuracy for large-eddy simulation applications on hybrid meshes. Detailed simulation comparisons are made between low-order hexahedral, tetrahedral, pyramid, and wedge/prism topologies against a third-order, unstructured hexahedral topology. Using smooth analytical and manufactured low-Mach solutions, design-order convergence is established for the hexahedral, tetrahedral, pyramid, and wedge element topologies using a new open boundary condition based on energy-stable methodologies previously deployed within a finite-difference context. A wide range of simulations demonstrate that low-order hexahedral- and wedge-based element topologies behave nearly identically in both computed numerical errors and overall simulation timings. Moreover, low-order tetrahedral and pyramid element topologies also display nearly the same numerical characteristics. Although the superiority of the hexahedral-based topology is clearly demonstrated for trivial laminar, principally-aligned flows, e.g., a 1x2x10 channel flow with specified pressure drop, this advantage is reduced for non-aligned, turbulent flows including the Taylor–Green Vortex, turbulent plane channel flow (Reτ395), and buoyant flow past a heated cylinder. With the order of accuracy demonstrated for both homogenous and hybrid meshes, it is shown that solution verification for the selected complex flows can be established for all topology types. Although the number of elements in a mesh of like spacing comprised of tetrahedral, wedge, or pyramid elements increases as compared to the hexahedral counterpart, for wall-resolved large-eddy simulation, the increased assembly and residual evaluation computational time for non-hexahedral is offset by more efficient linear solver times. Finally, most simulation results indicate that modest polynomial promotion provides a significant increase in solution accuracy.

[1]  Alison L. Marsden,et al.  Construction of commutative filters for LES on unstructured meshes , 2000 .

[2]  Christer Fureby,et al.  Large-Eddy Simulation: Current Capabilities, Recommended Practices, and Future Research , 2009 .

[3]  R. Verzicco,et al.  Large Eddy Simulation in Complex Geometric Configurations Using Boundary Body Forces , 2000 .

[4]  P. Moin,et al.  Unstructured Large Eddy Simulation for Prediction of Noise Issued from Turbulent Jets in Various Configurations , 2011 .

[5]  Sourabh V. Apte,et al.  Large-Eddy Simulation of Realistic Gas Turbine Combustors , 2004 .

[6]  Christopher C. Pain,et al.  A finite element LES methodology for anisotropic inhomogeneous meshes , 2012 .

[7]  Stefan P. Domino,et al.  A Framework for Characterizing Structural Uncertainty in Large-Eddy Simulation Closures , 2018 .

[8]  Antony Jameson,et al.  High-Order Flux Reconstruction Schemes for LES on Tetrahedral Meshes , 2015 .

[9]  Magnus Svärd,et al.  A stable high-order finite difference scheme for the compressible Navier-Stokes equations: No-slip wall boundary conditions , 2008, J. Comput. Phys..

[10]  Kenneth E. Jansen,et al.  Large-Eddy Simulation of Flow Around a NACA 4412 Airfoil Using Unstructured Grids , 1996 .

[11]  G. Taylor,et al.  Mechanism of the production of small eddies from large ones , 1937 .

[12]  G. Stipcich,et al.  Discontinuous control-volume/finite-element method for advection–diffusion problems , 2011 .

[13]  Stefan P. Domino,et al.  Design-order, non-conformal low-Mach fluid algorithms using a hybrid CVFEM/DG approach , 2018, J. Comput. Phys..

[14]  P. Quemere,et al.  Verification and validation considerations regarding the qualification of numerical schemes for LES for dilution problems , 2010 .

[15]  John Kim,et al.  DIRECT NUMERICAL SIMULATION OF TURBULENT CHANNEL FLOWS UP TO RE=590 , 1999 .

[16]  P. Roache.  Code Verification by the Method of Manufactured Solutions , 2002 .

[17]  Thomas Höhne,et al.  Modeling of a buoyancy-driven flow experiment at the ROCOM test facility using the CFD codes CFX-5 and Trio_U , 2006 .

[18]  R. Codina.  Pressure Stability in Fractional Step Finite Element Methods for Incompressible Flows , 2001 .

[19]  Thomas W. Grasser,et al.  Experimental investigation of a cylinder in turbulent thermal convection with an imposed shear flow. , 2005 .

[20]  Thomas S. Lund,et al.  Kinetic energy conservation issues associated with the collocated mesh scheme for incompressible flow , 2006, J. Comput. Phys..

[21]  Todd Simons,et al.  Large eddy simulation of turbulent flows using unstructured grids , 1998 .

[22]  Onkar Sahni,et al.  Finite element based large eddy simulation of flow over bluff bodies , 2017 .

[23]  Arif Masud,et al.  Residual-based variational multiscale turbulence models for unstructured tetrahedral meshes , 2013 .

[24]  Koen Hillewaert,et al.  Assessment of a discontinuous Galerkin method for the simulation of vortical flows at high Reynolds number , 2014 .

[25]  Gianluca Iaccarino,et al.  DNS of buoyancy-dominated turbulent flows on a bluff body using the immersed boundary method , 2009, J. Comput. Phys..

[26]  James C. Sutherland,et al.  Verification for multi-mechanics applications , 2007 .

[27]  F. Nicoud,et al.  Subgrid-Scale Stress Modelling Based on the Square of the Velocity Gradient Tensor , 1999 .

[28]  P. Moin,et al.  A numerical method for large-eddy simulation in complex geometries , 2004 .

[29]  Gerry E. Schneider,et al.  A SKEWED, POSITIVE INFLUENCE COEFFICIENT UPWINDING PROCEDURE FOR CONTROL-VOLUME-BASED FINITE-ELEMENT CONVECTION-DIFFUSION COMPUTATION , 1986 .

[30]  R. Koch,et al.  Compressible large eddy simulation of turbulent combustion in complex geometry on unstructured meshes , 2004 .

[31]  Sean P. Kearney,et al.  Mixed convection heat transfer to and from a horizontal cylinder in cross-flow with heating from below. , 2006 .

[32]  Simon P. Walker,et al.  Polynomial basis functions on pyramidal elements , 2007 .

[33]  Parviz Moin,et al.  ADVANCES IN LARGE EDDY SIMULATION METHODOLOGY FOR COMPLEX FLOWS , 2002, Proceeding of Second Symposium on Turbulence and Shear Flow Phenomena.

[34]  Malte Braack,et al.  Finite elements with local projection stabilization for incompressible flow problems , 2009 .

[35]  Oriol Lehmkuhl,et al.  Unsteady forces on a circular cylinder at critical Reynolds numbers , 2014 .

[36]  C. Rhie,et al.  Numerical Study of the Turbulent Flow Past an Airfoil with Trailing Edge Separation , 1983 .

[37]  Fotini Katopodes Chow,et al.  The effect of numerical errors and turbulence models in large-eddy simulations of channel flow, with and without explicit filtering , 2003, Journal of Fluid Mechanics.