Parallel adaptive mesh refinement method based on WENO finite difference scheme for the simulation of multi-dimensional detonation

For numerical simulation of detonation, computational cost using uniform meshes is large due to the vast separation in both time and space scales. Adaptive mesh refinement (AMR) is advantageous for problems with vastly different scales. This paper aims to propose an AMR method with high order accuracy for numerical investigation of multi-dimensional detonation. A well-designed AMR method based on finite difference weighted essentially non-oscillatory (WENO) scheme, named as AMR&WENO is proposed. A new cell-based data structure is used to organize the adaptive meshes. The new data structure makes it possible for cells to communicate with each other quickly and easily. In order to develop an AMR method with high order accuracy, high order prolongations in both space and time are utilized in the data prolongation procedure. Based on the message passing interface (MPI) platform, we have developed a workload balancing parallel AMR&WENO code using the Hilbert space-filling curve algorithm. Our numerical experiments with detonation simulations indicate that the AMR&WENO is accurate and has a high resolution. Moreover, we evaluate and compare the performance of the uniform mesh WENO scheme and the parallel AMR&WENO method. The comparison results provide us further insight into the high performance of the parallel AMR&WENO method.

[1]  K. Benkiewicz,et al.  Two-dimensional numerical simulations of multi-headed detonations in oxygen-aluminum mixtures using an adaptive mesh refinement , 2003 .

[2]  Jianxian Qiu,et al.  Simulations of detonation wave propagation in rectangular ducts using a three-dimensional WENO scheme , 2008 .

[3]  Andrew J. Majda,et al.  Theoretical and numerical structure for unstable one-dimensional detonations , 1991 .

[4]  S. Sivier,et al.  Adaptive unstructured finite element method for two-dimensional detonation simulations , 1998 .

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

[6]  Antonio Baeza,et al.  Adaptive mesh refinement techniques for high‐order shock capturing schemes for multi‐dimensional hydrodynamic simulations , 2006 .

[7]  Ralf Deiterding,et al.  An adaptive high-order hybrid scheme for compressive, viscous flows with detailed chemistry , 2011, J. Comput. Phys..

[8]  Ralf Deiterding,et al.  Parallel adaptive simulation of multi-dimensional detonation structures , 2003 .

[9]  Gary J. Sharpe,et al.  Transverse waves in numerical simulations of cellular detonations , 2001, Journal of Fluid Mechanics.

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

[11]  Elaine S. Oran,et al.  Numerical simulation of detonation initiation in a flame brush : The role of hot spots , 1999 .

[12]  Jianguo Ning,et al.  High Resolution WENO Simulation of 3D Detonation Waves , 2013 .

[13]  B C Khoo,et al.  The cellular structure of a two-dimensional H2/O2/Ar detonation wave , 2004 .

[14]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

[15]  A. M. Khokhlov,et al.  Fully Threaded Tree for Adaptive Refinement Fluid Dynamics Simulations , 1998 .

[16]  Thomas Huld,et al.  An adaptive 3-D CFD solver for modeling explosions on large industrial environmental scales , 1999 .

[17]  Alexei M. Khokhlov,et al.  Fully Threaded Tree Algorithms for Adaptive Refinement Fluid Dynamics Simulations , 1997, astro-ph/9701194.

[18]  Hua Ji,et al.  A new adaptive mesh refinement data structure with an application to detonation , 2010, J. Comput. Phys..

[19]  Ronald Fedkiw,et al.  Numerical resolution of pulsating detonation waves , 2000 .

[20]  Elaine S. Oran,et al.  A numerical study of the mechanisms of self-reignition in low-overdrive detonations , 1996 .

[21]  Fue-Sang Lien,et al.  Numerical simulation of detonation using an adaptive Cartesian cut-cell method combined with a cell-merging technique , 2010 .

[22]  Manish Parashar,et al.  An Application-Centric Characterization of Domain-Based SFC Partitioners for Parallel SAMR , 2002, IEEE Trans. Parallel Distributed Syst..

[23]  M. Dumbser,et al.  High order space–time adaptive ADER-WENO finite volume schemes for non-conservative hyperbolic systems , 2013, 1304.5408.

[24]  J. Hyman,et al.  Adaptive Mesh Refinement for Finite Difference WENO Schemes , 2022 .

[25]  William D. Henshaw,et al.  Parallel computation of three-dimensional flows using overlapping grids with adaptive mesh refinement , 2008, J. Comput. Phys..

[26]  Chaopeng Shen,et al.  Adaptive mesh refinement based on high order finite difference WENO scheme for multi-scale simulations , 2011, J. Comput. Phys..

[27]  Chi-Wang Shu,et al.  Multidomain WENO Finite Difference Method with Interpolation at Subdomain Interfaces , 2003, J. Sci. Comput..

[28]  W. Fickett,et al.  Flow Calculations for Pulsating One‐Dimensional Detonations , 1966 .

[29]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .