Related papers: Efficient time stepping for reactive turbulent sim…
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…
The goal of this work is to develop a novel splitting approach for the numerical solution of multiscale problems involving the coupling between Stokes equations and ODE systems, as often encountered in blood flow modeling applications. The…
The nonlinear gyrokinetic equations describe plasma turbulence in laboratory and astrophysical plasmas. To solve these equations, massively parallel codes have been developed and run on present-day supercomputers. This paper describes…
This paper introduces an adaptive time splitting technique for the solution of stiff evolutionary PDEs that guarantees an effective error control of the simulation, independent of the fastest physical time scale for highly unsteady…
The efficiency of exact simulation methods for the reaction-diffusion master equation (RDME) is severely limited by the large number of diffusion events if the mesh is fine or if diffusion constants are large. Furthermore, inherent…
Stochastic reaction networks governed by Chemical Langevin Equations (CLE) exhibit pronounced multiscale dynamics spanning fast molecular reactions, intermediate transport, and slow cellular regulation, posing significant challenges for…
Tau-leaping is a family of algorithms for the approximate simulation of the discrete state continuous time Markov chains. Motivation for the development of such methods can be found, for instance, in the fields of chemical kinetics and…
A dual time-stepping strategy of gas-kinetic scheme for the simulation of unsteady flows is introduced in this work. The dual time-stepping strategy is always used in the approaches of unsteady flows, and the ability of dual time-stepping…
Many problems in science and engineering require an efficient numerical approximation of integrals or solutions to differential equations. For systems with rapidly changing dynamics, an equidistant discretization is often inadvisable as it…
A component-splitting method is proposed to improve convergence characteristics for implicit time integration of compressible multicomponent reactive flows. The characteristic decomposition of flux jacobian of multicomponent Navier-Stokes…
A second-order accurate in time, positivity-preserving, and unconditionally energy stable operator splitting numerical scheme is proposed and analyzed for the system of reaction-diffusion equations with detailed balance. The scheme is…
We address the problem of simulation and parameter inference for chemical reaction networks described by the chemical Langevin equation, a stochastic differential equation (SDE) representation of the dynamics of the chemical species. This…
A new approach for integration of motion in many-body systems of interacting polyatomic molecules is proposed. It is based on splitting time propagation of pseudo-variables in a modified phase space, while the real translational and…
We present a detailed convergence analysis for an operator splitting scheme proposed in [C. Liu et al.,J. Comput. Phys., 436, 110253, 2021] for a reaction-diffusion system with detailed balance. The numerical scheme has been constructed…
In addressing the demands of industrial high-fidelity computation, the present study introduces a rapid and accurate customized solver developed on the OpenFOAM platform. To enhance computational efficiency, a novel integrated acceleration…
The computational cost of exact methods for quantum simulation using classical computers grows exponentially with system size. As a consequence, these techniques can only be applied to small systems. By contrast, we demonstrate that quantum…
A fourth-order exponential time differencing (ETD) Runge-Kutta scheme with dimensional splitting is developed to solve multidimensional non-linear systems of reaction-diffusion equations (RDE). By approximating the matrix exponential in the…
Simulating complex astrophysical reacting flows is computationally expensive -- reactions are stiff and typically require implicit integration methods. The reaction update is often the most expensive part of a simulation, which motivates…
A fast and robust Jacobian-free time-integration method - called Minimum-error Adaptation of a Chemical-Kinetic ODE Solver (MACKS) - for solving stiff ODEs pertaining to chemical-kinetics is proposed herein. The MACKS formulation is based…
Models invoking the chemical master equation are used in many areas of science, and, hence, their simulation is of interest to many researchers. The complexity of the problems at hand often requires considerable computational power, so a…