A low-numerical dissipation, patch-based adaptive-mesh-refinement method for large-eddy simulation of compressible flows

This paper describes a hybrid finite-difference method for the large-eddy simulation of compressible flows with low-numerical dissipation and structured adaptive mesh refinement (SAMR). A conservative flux-based approach is described with an explicit centered scheme used in turbulent flow regions while a weighted essentially non-oscillatory (WENO) scheme is employed to capture shocks. Three-dimensional numerical simulations of a Richtmyer-Meshkov instability are presented

Detailed Simulations of Shock-Bifurcation and Ignition of an Argon-diluted Hydrogen/Oxygen Mixture in a Shock Tube

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/106508/1/AIAA2013-538.pd

Hybrid Spectral Difference/Embedded Finite Volume Method for Conservation Laws

A novel hybrid spectral difference/embedded finite volume method is introduced in order to apply a discontinuous high-order method for large scale engineering applications involving discontinuities in the flows with complex geometries. In the proposed hybrid approach, the finite volume (FV) element, consisting of structured FV subcells, is embedded in the base hexahedral element containing discontinuity, and an FV based high-order shock-capturing scheme is employed to overcome the Gibbs phenomena. Thus, a discontinuity is captured at the resolution of FV subcells within an embedded FV element. In the smooth flow region, the SD element is used in the base hexahedral element. Then, the governing equations are solved by the SD method. The SD method is chosen for its low numerical dissipation and computational efficiency preserving high-order accurate solutions. The coupling between the SD element and the FV element is achieved by the globally conserved mortar method. In this paper, the 5th-order WENO scheme with the characteristic decomposition is employed as the shock-capturing scheme in the embedded FV element, and the 5th-order SD method is used in the smooth flow field. The order of accuracy study and various 1D and 2D test cases are carried out, which involve the discontinuities and vortex flows. Overall, it is shown that the proposed hybrid method results in comparable or better simulation results compared with the standalone WENO scheme when the same number of solution DOF is considered in both SD and FV elements.Comment: 27 pages, 17 figures, 2 tables, Accepted for publication in the Journal of Computational Physics, April 201

Detonation wave diffraction in H₂-O₂-Ar mixtures

In the present study, we have examined the diffraction of detonation in weakly unstable hydrogen–oxygen–argon mixtures. High accuracy and computational efficiency are obtained using a high-order WENO scheme together with adaptive mesh refinement, which enables handling realistic geometries with resolution at the micrometer level. Both detailed chemistry and spectroscopic models of laser induced fluorescence and chemiluminescence were included to enable a direct comparison with experimental data. Agreement was found between the experiments and the simulations in terms of detonation diffraction structure both for sub-critical and super-critical regimes. The predicted wall reflection distance is about 12–14 cell widths, in accordance with previous experimental studies. Computations show that the re-initiation distance is relatively constant, at about 12–15 cell widths, slightly above the experimental value of 11 cell widths. The predicted critical channel height is 10–11 cell widths, which differs from experiments in circular tubes but is consistent with rectangular channel results

An Application of Gaussian Process Modeling for High-order Accurate Adaptive Mesh Refinement Prolongation

We present a new polynomial-free prolongation scheme for Adaptive Mesh Refinement (AMR) simulations of compressible and incompressible computational fluid dynamics. The new method is constructed using a multi-dimensional kernel-based Gaussian Process (GP) prolongation model. The formulation for this scheme was inspired by the GP methods introduced by A. Reyes et al. (A New Class of High-Order Methods for Fluid Dynamics Simulation using Gaussian Process Modeling, Journal of Scientific Computing, 76 (2017), 443-480; A variable high-order shock-capturing finite difference method with GP-WENO, Journal of Computational Physics, 381 (2019), 189-217). In this paper, we extend the previous GP interpolations and reconstructions to a new GP-based AMR prolongation method that delivers a high-order accurate prolongation of data from coarse to fine grids on AMR grid hierarchies. In compressible flow simulations special care is necessary to handle shocks and discontinuities in a stable manner. To meet this, we utilize the shock handling strategy using the GP-based smoothness indicators developed in the previous GP work by A. Reyes et al. We demonstrate the efficacy of the GP-AMR method in a series of testsuite problems using the AMReX library, in which the GP-AMR method has been implemented

Development of a Chemically Reacting Flow Solver on the Graphic Processing Units

The focus of the current research is to develop a numerical framework on the Graphic Processing Units (GPU) capable of modeling chemically reacting flow. The framework incorporates a high-order finite volume method coupled with an implicit solver for the chemical kinetics. Both the fluid solver and the kinetics solver are designed to take advantage of the GPU architecture to achieve high performance. The structure of the numerical framework is shown, detailing different aspects of the optimization implemented on the solver. The mathematical formulation of the core algorithms is presented along with a series of standard test cases, including both nonreactive and reactive flows, in order to validate the capability of the numerical solver. The performance results obtained with the current framework show the parallelization efficiency of the solver and emphasize the capability of the GPU in performing scientific calculations. Distribution A: Approved for public release; distribution unlimited. PA #1117

Spurious behavior of shock-capturing methods by the fractional step approach: Problems containing stiff source terms and discontinuities

The goal of this paper is to relate numerical dissipations that are inherited in high order shock-capturing schemes with the onset of wrong propagation speed of discontinuities. For pointwise evaluation of the source term, previous studies indicated that the phenomenon of wrong propagation speed of discontinuities is connected with the smearing of the discontinuity caused by the discretization of the advection term. The present study focuses only on solving the reactive system by the fractional step method using the Strang splitting. Studies shows that the degree of wrong propagation speed of discontinuities is highly dependent on the accuracy of the numerical method. The manner in which the smearing of discontinuities is contained by the numerical method and the overall amount of numerical dissipation being employed play major roles. Depending on the numerical method, time step and grid spacing, the numerical simulation may lead to (a) the correct solution (within the truncation error of the scheme), (b) a divergent solution, (c) a wrong propagation speed of discontinuities solution or (d) other spurious solutions that are solutions of the discretized counterparts but are not solutions of the governing equations. The findings might shed some light on the reported difficulties in numerical combustion and problems with stiff nonlinear (homogeneous) source terms and discontinuities in general