A split random time-stepping method for stiff and nonstiff detonation capturing

Abstract

In this paper, a new fractional step method is proposed for simulating stiff and nonstiff chemically reacting flows. In stiff cases, a well-known spurious numerical phenomenon, i.e. the incorrect propagation speed of discontinuities, may be produced by general fractional step methods due to the under-resolved discretization in both space and time. The previous random projection method has been successfully applied for stiff detonation capturing in under-resolved conditions. Not to randomly project the intermediate state into two presumed equilibrium states (completely burnt or unburnt) as in the random projection method, the present study is to randomly choose the time-dependent advance or stop of a reaction process. Each one-way reaction has been decoupled from the multi-reaction kinetics using operator splitting and the local smeared temperature due to numerical dissipation of shock-capturing schemes is compared with a random one within two limited temperatures corresponding to the advance and its inverse states, respectively, to control the random reaction. The random activation or deactivation in the reaction step is thus promising to correct the deterministic accumulative error of the propagation of discontinuities. Extensive numerical experiments, including model problems and realistic reacting flows in one and two dimensions, demonstrate this expectation as well as the effectiveness and robustness of the method. Meanwhile, for nonstiff problems when spatial and temporal resolutions are fine, the proposed random method recovers the results as general fractional step methods, owing to the increasing possibility of activation with diminishing randomness by adding a shift term.

Review on Smoothed Particle Hydrodynamics: Methodology development and recent achievement

“…This limitation may lead excessive computational efforts especially when large damping ratio and high resolution ratio are applied. To release this limitation, an operator splitting scheme [168,22] is first applied to decouple the momentum conservation equation Eq. ( 88) into the original momentum part and the damping part.…”

An integrative smoothed particle hydrodynamics framework for modeling cardiac function

“…In this work, we employ a reaction-by-reaction splitting method proposed by Wang et al [24]. The multi-reaction system can be decoupled, e.g.…”

“…This consideration also leads to a much larger time step size compared to the simple forward Euler method [2]. Then, we introduce a splitting reaction-by-reaction method [24] combined with quasi-steady-state (QSS) solver to capture the stiff dynamics of the transmembrane potential and the gating variables of the ionic model governed by nonlinear ordinary differential equations (ODEs). Furthermore, an anisotropic SPH discretization for diffusion equation derived by Tran-Duc et al [25] is modified by introducing a linear operator to improve the computational efficiency and using a correction kernel matrix to improve the numerical accuracy.…”

SPHinXsys: an open-source multi-physics and multi-resolution library based on smoothed particle hydrodynamics

“…In SPHinXsys, a reaction-by-reaction splitting method [29] is introduced for solving the system of ODEs defined by Eq. ( 47), which is generally stiff and induces numerical instability when the integration time step is not sufficiently small.…”