An adaptive high-order hybrid scheme for compressive, viscous flows with detailed chemistry
A hybrid weighted essentially non-oscillatory (WENO)/centered-difference numerical method, with low numerical dissipation, high-order shock-capturing, and structured adaptive mesh refinement (SAMR), has been developed for the direct numerical simulation of the multicomponent, compressible, reactive Navier–Stokes equations. The method enables accurate resolution of diffusive processes within reaction zones. The approach combines time-split reactive source terms with a high-order, shock-capturing scheme specifically designed for diffusive flows. A description of the order-optimized, symmetric, finite difference, flux-based, hybrid WENO/centered-difference scheme is given, along with its implementation in a high-order SAMR framework. The implementation of new techniques for discontinuity flagging, scheme-switching, and high-order prolongation and restriction is described. In particular, the refined methodology does not require upwinded WENO at grid refinement interfaces for stability, allowing high-order prolongation and thereby eliminating a significant source of numerical diffusion within the overall code performance. A series of one-and two-dimensional test problems is used to verify the implementation, specifically the high-order accuracy of the diffusion terms. One-dimensional benchmarks include a viscous shock wave and a laminar flame. In two-space dimensions, a Lamb–Oseen vortex and an unstable diffusive detonation are considered, for which quantitative convergence is demonstrated. Further, a two-dimensional high-resolution simulation of a reactive Mach reflection phenomenon with diffusive multi-species mixing is presented.
© 2011 Elsevier Inc. Received 18 June 2010; revised 16 June 2011; Accepted 16 June 2011. Available online 25 June 2011. Jack Ziegler is supported by the Department of Energy Computational Science Graduate Fellowship program (DOE CSGF). This research used resources of the National Energy Research Scientific Computing (NERSC) Center. D. Pullin and R. Deiterding were partially supported by the Department of Energy Advanced Scientific and Computing (ASC) program under subcontract B341492 of DOE contract W-7405-ENG-48.
Accepted Version - ziegler_jcp_2011.pdf