High-Order Hybrid WCNS-CPR Scheme for Shock Capturing of Conservation Laws

By introducing hybrid technique into high-order CPR (correction procedure via reconstruction) scheme, a novel hybrid WCNS-CPR scheme is developed for efficient supersonic simulations. Firstly, a shock detector based on nonlinear weights is used to identify grid cells with high gradients or discontinuities throughout the whole flow field. Then, WCNS (weighted compact nonlinear scheme) is adopted to capture shocks in these areas, while the smooth area is calculated by CPR. A strategy to treat the interfaces of the two schemes is developed, which maintains high-order accuracy. Convergent order of accuracy and shock-capturing ability are tested in several numerical experiments; the results of which show that this hybrid scheme achieves expected high-order accuracy and high resolution, is robust in shock capturing, and has less computational cost compared to the WCNS.

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

[2]  Mengping Zhang,et al.  A simple weighted essentially non-oscillatory limiter for the correction procedure via reconstruction (CPR) framework on unstructured meshes , 2015 .

[3]  Mengping Zhang,et al.  A simple weighted essentially non-oscillatory limiter for the correction procedure via reconstruction (CPR) framework , 2015 .

[4]  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..

[5]  Chi-Wang Shu,et al.  TVB Runge-Kutta local projection discontinuous galerkin finite element method for conservation laws. II: General framework , 1989 .

[6]  Tiegang Liu,et al.  A GENERAL HIGH-ORDER MULTI-DOMAIN HYBRID DG/WENO-FD METHOD FOR HYPERBOLIC CONSERVATION LAWS * , 2016 .

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

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

[9]  Tiegang Liu,et al.  Multidomain Hybrid RKDG and WENO Methods for Hyperbolic Conservation Laws , 2013, SIAM J. Sci. Comput..

[10]  Feng Qu,et al.  A new all-speed flux scheme for the Euler equations , 2019, Comput. Math. Appl..

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

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

[13]  G. Sod.  A survey of several finite difference methods for systems of nonlinear hyperbolic conservation laws , 1978 .

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

[15]  Taku Nonomura,et al.  Conservative high-order flux-reconstruction schemes on moving and deforming grids , 2016 .

[16]  Minwei Wu,et al.  Direct numerical simulation of supersonic turbulent boundary layer over a compression ramp , 2007 .

[17]  Chunlei Liang,et al.  A high-order flux reconstruction method with adaptive mesh refinement and artificial diffusivity on unstructured moving/deforming mesh for shock capturing , 2016 .

[18]  Junqiang Bai,et al.  Self-similar structures based genuinely two-dimensional Riemann solvers in curvilinear coordinates , 2020, J. Comput. Phys..

[19]  Feng Qu,et al.  An efficient adaptive high-order scheme based on the WENO process , 2016 .

[20]  Xiaogang Deng,et al.  Osher Flux with Entropy Fix for Two-Dimensional Euler Equations , 2016 .

[21]  Wai-Sun Don,et al.  An improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws , 2008, J. Comput. Phys..

[22]  J. Remacle,et al.  Shock detection and limiting with discontinuous Galerkin methods for hyperbolic conservation laws , 2004 .

[23]  Wanai Li,et al.  A p-weighted limiter for the discontinuous Galerkin method on one-dimensional and two-dimensional triangular grids , 2020, J. Comput. Phys..

[24]  Zhi J. Wang,et al.  Towards industrial large eddy simulation using the FR/CPR method , 2017 .

[25]  Yifeng Zhang,et al.  High-Order and High Accurate CFD Methods and Their Applications for Complex Grid Problems , 2012 .

[26]  Xiaogang Deng,et al.  Developing high-order weighted compact nonlinear schemes , 2000 .

[27]  Zhenguo Yan,et al.  High-Order Hybrid WCNS-CPR Schemes on Hybrid Meshes with Curved Edges for Conservation Laws I: Spatial Accuracy and Geometric Conservation Laws , 2018 .

[28]  P. Woodward,et al.  The numerical simulation of two-dimensional fluid flow with strong shocks , 1984 .

[29]  P. Lax.  Weak solutions of nonlinear hyperbolic equations and their numerical computation , 1954 .

[30]  Chaowei Hu,et al.  No . 98-32 Weighted Essentially Non-Oscillatory Schemes on Triangular Meshes , 1998 .

[31]  Chao Zhang,et al.  A third-order gas-kinetic CPR method for the Euler and Navier-Stokes equations on triangular meshes , 2018, J. Comput. Phys..

[32]  Lu Qilin,et al.  A Parameter-Free Gradient-Based Limiter for the FR/CPR method on Mixed Unstructured Meshes , 2019 .

[33]  Michael Dumbser,et al.  A posteriori subcell limiting of the discontinuous Galerkin finite element method for hyperbolic conservation laws , 2014, J. Comput. Phys..