rp‐adaptation for compressible flows

We present a novel rp‐adaptation strategy for high‐fidelity simulations of compressible inviscid flows with shocks. The mesh resolution in regions of flow discontinuities is increased by using a variational optimizer to r‐adapt the mesh and cluster degrees of freedom there. In regions of smooth flow, we locally increase or decrease the local resolution through increasing or decreasing the polynomial order of the elements, respectively. This dual approach allows us to take advantage of the strengths of both methods for best computational performance, thereby reducing the overall cost of the simulation. The adaptation workflow uses a sensor for both discontinuities and smooth regions that is cheap to calculate, but the framework is general and could be used in conjunction with other feature‐based sensors or error estimators. We demonstrate this proof‐of‐concept using two geometries in transonic and supersonic flow regimes. The method has been implemented in the open‐source spectral/hp element framework Nektar++, and adaptivity is performed by its dedicated high‐order mesh generation tool NekMesh. The results show that the proposed rp‐adaptation methodology is a reasonably cost‐effective way of improving simulation accuracy.

[1]  N. D. Wong,et al.  Experimental investigation of a large scale, two-dimensional, mixed-compression inlet system - Performance at design conditions, M=3.0 , 1970 .

[2]  N. D. Wong,et al.  Experimental investigation of a large scale, two-dimensional, mixed compression inlet system - Performance at supersonic conditions, M equals 1.55 to 3.2 , 1971 .

[3]  J. T. Oden,et al.  An hp adaptive strategy for finite element approximations of the Navier‐Stokes equations , 1995 .

[4]  E. Toro.  Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .

[5]  S. Rebay,et al.  A High-Order Accurate Discontinuous Finite Element Method for the Numerical Solution of the Compressible Navier-Stokes Equations , 1997 .

[6]  Chi-Wang Shu,et al.  The Local Discontinuous Galerkin Method for Time-Dependent Convection-Diffusion Systems , 1998 .

[7]  Sorin Mitran.  A Comparison of Adaptive Mesh Refinement Approaches for Large Eddy Simulation , 2001 .

[8]  Pierre Sagaut,et al.  A dynamic p-adaptive Discontinuous Galerkin method for viscous flow with shocks , 2005 .

[9]  Sanjay Mittal,et al.  Euler flow in a supersonic mixed‐compression inlet , 2006 .

[10]  J. Peraire,et al.  Sub-Cell Shock Capturing for Discontinuous Galerkin Methods , 2006 .

[11]  John C. Vassberg,et al.  In Pursuit of Grid Convergence for Two-Dimensional Euler Solutions , 2009 .

[12]  David L. Darmofal,et al.  Shock capturing with PDE-based artificial viscosity for DGFEM: Part I. Formulation , 2010, J. Comput. Phys..

[13]  Antony Jameson,et al.  Comparison of Adaptive h and p Refinements for Spectral Difference Methods , 2010 .

[14]  Viscous Shock Capturing in a Time-Explicit Discontinuous Galerkin Method , 2011, 1102.3190.

[15]  D. Darmofal,et al.  Review of Output-Based Error Estimation and Mesh Adaptation in Computational Fluid Dynamics , 2011 .

[16]  Robert Michael Kirby,et al.  Nektar++: An open-source spectral/hp element framework , 2015, Comput. Phys. Commun..

[17]  Spencer J. Sherwin,et al.  Automatic generation of 3D unstructured high-order curvilinear meshes , 2016 .

[18]  M. A. O. Ignacio,et al.  How to cite this article , 2016 .

[19]  Spencer J. Sherwin,et al.  A p-adaptation method for compressible flow problems using a goal-based error indicator , 2017 .

[20]  David Moxey,et al.  Towards p-Adaptive Spectral/hp Element Methods for Modelling Industrial Flows , 2017 .

[21]  Neil D. Sandham,et al.  An error indicator for finite difference methods using spectral techniques with application to aerofoil simulation , 2018 .

[22]  Spencer J. Sherwin,et al.  Direct numerical simulation of a compressible boundary-layer flow past an isolated three-dimensional hump in a high-speed subsonic regime , 2018 .

[23]  David Moxey,et al.  Curvilinear mesh generation using a variational framework , 2017, Comput. Aided Des..

[24]  Frédéric Coquel,et al.  A comparison of refinement indicators for p-adaptive simulations of steady and unsteady flows using discontinuous Galerkin methods , 2019, J. Comput. Phys..

[25]  Spencer J. Sherwin,et al.  A variational approach to high-order r-adaptation , 2019, ArXiv.

[26]  Philipp Schlatter,et al.  Adaptive mesh refinement for steady flows in Nek5000 , 2020 .

[27]  Robert Michael Kirby,et al.  Nektar++: enhancing the capability and application of high-fidelity spectral/hp element methods , 2019, Comput. Phys. Commun..