A high-order flux reconstruction method with adaptive mesh refinement and artificial diffusivity on unstructured moving/deforming mesh for shock capturing

Abstract This article presents new progress of shock capturing on unstructured dynamic grids by using adaptive mesh refinement (AMR) method in conjunction with artificial diffusivity (AD) under the framework of flux reconstruction (FR)/correction procedure via reconstruction (CPR) method. The proposed AMR algorithm is applicable even when the grid undergoes dynamic motion. It features an innovative data structure that is capable of direct addressing in managing the system of cells on multiple levels of refinement. A conservative mortar method is used to handle non-conforming interfaces resulted from cell subdivision. The AD is added in regions of strong dilatation to counter the Gibbs oscillation that is disruptive to high-order solutions when shock discontinuities exist. Several test cases are used to verify that this FR-AMR-AD framework can achieve high-order accuracy for smooth flows, and obtain stable solution with the presence of shock discontinuity. By further limiting the shock within subdivided adjacent cells, a sharp shock solution can be obtained with reduced amount of AD addition, hence errors from adding dissipation, on dynamically refined grids.

[1]  Hui Hu,et al.  A high-order spectral difference method for unstructured dynamic grids , 2011 .

[2]  A. Jameson.  ANALYSIS AND DESIGN OF NUMERICAL SCHEMES FOR GAS DYNAMICS, 1: ARTIFICIAL DIFFUSION, UPWIND BIASING, LIMITERS AND THEIR EFFECT ON ACCURACY AND MULTIGRID CONVERGENCE , 1995 .

[3]  Marsha Berger,et al.  Three-Dimensional Adaptive Mesh Refinement for Hyperbolic Conservation Laws , 1994, SIAM J. Sci. Comput..

[4]  John H. Kolias,et al.  A CONSERVATIVE STAGGERED-GRID CHEBYSHEV MULTIDOMAIN METHOD FOR COMPRESSIBLE FLOWS , 1995 .

[5]  Steven J. Ruuth,et al.  A New Class of Optimal High-Order Strong-Stability-Preserving Time Discretization Methods , 2002, SIAM J. Numer. Anal..

[6]  R. D. Richtmyer,et al.  A Method for the Numerical Calculation of Hydrodynamic Shocks , 1950 .

[7]  P. Colella,et al.  Local adaptive mesh refinement for shock hydrodynamics , 1989 .

[8]  Bin Zhang,et al.  An efficient correction procedure via reconstruction for simulation of viscous flow on moving and deforming domains , 2014, J. Comput. Phys..

[9]  H. T. Huynh,et al.  High-Order Methods for Computational Fluid Dynamics: A Brief Review of Compact Differential Formulations on Unstructured Grids , 2013 .

[10]  R. Anderson,et al.  An arbitrary Lagrangian-Eulerian method with adaptive mesh refinement for the solution of the Euler equations , 2004 .

[11]  Friedrich L. Bauer,et al.  Supercritical Wing Sections II , 1974 .

[12]  Francis X. Giraldo,et al.  Analysis of adaptive mesh refinement for IMEX discontinuous Galerkin solutions of the compressible Euler equations with application to atmospheric simulations , 2014, J. Comput. Phys..

[13]  Soshi Kawai,et al.  Localized artificial diffusivity scheme for discontinuity capturing on curvilinear meshes , 2008, J. Comput. Phys..

[14]  Philip L. Roe,et al.  An Entropy Adjoint Approach to Mesh Refinement , 2010, SIAM J. Sci. Comput..

[15]  Lei Shi,et al.  Adjoint-based error estimation and mesh adaptation for the correction procedure via reconstruction method , 2015, J. Comput. Phys..

[16]  Zhi Jian Wang,et al.  A unifying lifting collocation penalty formulation including the discontinuous Galerkin, spectral volume/difference methods for conservation laws on mixed grids , 2009, J. Comput. Phys..

[17]  H. Tijdeman,et al.  Transonic Flow Past Oscillating Airfoils , 1980 .

[18]  Chunlei Liang,et al.  A comparison of computational efficiencies of spectral difference method and correction procedure via reconstruction , 2013, J. Comput. Phys..

[19]  R J Magnus.  Some Numerical Solutions of Inviscid, Unsteady, Transonic Flows over the NLR 7301 Airfoil. , 1978 .

[20]  J. B. Mcdevitt,et al.  Static and dynamic pressure measurements on a NACA 0012 airfoil in the Ames High Reynolds Number Facility , 1985 .

[21]  Soshi Kawai,et al.  Toward Accurate Simulation of Shockwave-Turbulence Interaction on Unstructured Meshes: A Coupling of High- Order FR and LAD Schemes , 2013 .

[22]  A. Jameson,et al.  Numerical solution of the Euler equations by finite volume methods using Runge Kutta time stepping schemes , 1981 .

[23]  Taku Nonomura,et al.  Fully-conservative High-order FR Scheme on Moving and Deforming Grids , 2015 .

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

[25]  David A. Kopriva,et al.  Discontinuous Galerkin spectral element approximations on moving meshes , 2011, J. Comput. Phys..

[26]  Chunlei Liang,et al.  Computation of flows with shocks using the Spectral Difference method with artificial viscosity, II: Modified formulation with local mesh refinement , 2014 .

[27]  M. Berger,et al.  Adaptive mesh refinement for hyperbolic partial differential equations , 1982 .

[28]  Sanjiva K. Lele,et al.  A modified artificial viscosity approach for compressible turbulence simulations , 2008, J. Comput. Phys..

[29]  D. Kopriva.  A Conservative Staggered-Grid Chebyshev Multidomain Method for Compressible Flows. II. A Semi-Structured Method , 1996 .

[30]  P. Roe.  Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

[31]  Bin Zhang,et al.  A simple, efficient, and high-order accurate curved sliding-mesh interface approach to spectral difference method on coupled rotating and stationary domains , 2015, J. Comput. Phys..

[32]  A. Jameson.  Solution of the Euler equations for two dimensional transonic flow by a multigrid method , 1983 .

[33]  W. Cabot,et al.  A high-wavenumber viscosity for high-resolution numerical methods , 2004 .

[34]  H. T. Huynh,et al.  A Flux Reconstruction Approach to High-Order Schemes Including Discontinuous Galerkin Methods , 2007 .

[35]  Marco Luciano Savini,et al.  Discontinuous Galerkin solution of the Reynolds-averaged Navier–Stokes and k–ω turbulence model equations , 2005 .

[36]  John B. Bell,et al.  A Cartesian Grid Projection Method for the Incompressible Euler Equations in Complex Geometries , 1994, SIAM J. Sci. Comput..

[37]  Q.F. Stout,et al.  Adaptive Blocks: A High Performance Data Structure , 1997, ACM/IEEE SC 1997 Conference (SC'97).

[38]  David S. Dolling,et al.  Fifty Years of Shock-Wave/Boundary-Layer Interaction Research: What Next? , 2001 .