Explicit time integration of the stiff chemical Langevin equations using computational singular perturbation

Abstract

A stable explicit time-scale splitting algorithm for stiff chemical Langevin equations (CLEs) is developed, based on the concept of computational singular perturbation. The drift term of the CLE is projected onto basis vectors that span the fast and slow subdomains. The corresponding fast modes exhaust quickly, in the mean sense, and the system state then evolves, with a mean drift controlled by slow modes, on a random manifold. The drift-driven time evolution of the state due to fast exhausted modes is modeled algebraically as an exponential decay process, while that due to slow drift modes and diffusional processes is integrated explicitly. This allows time integration step sizes much larger than those required by typical explicit numerical methods for stiff stochastic differential equations. The algorithm is motivated and discussed, and extensive numerical experiments are conducted to illustrate its accuracy and stability with a number of model systems.

Localization of Noise in Biochemical Networks

“…It greatly reduces the computational cost compared to SSA, but it can be less accurate . Some implementations that shift between SSA and τ-leaping have been developed to prevent most of the loss of accuracy. , In terms of speed, a less accurate, but faster approach than τ-leaping is the chemical Langevin equation. − With this approach, the error can be large for some systems, especially if the step size is not properly chosen. , …”

Optimal explicit stabilized postprocessed $τ$-leap method for the simulation of chemical kinetics

“…Roughly speaking, most of such schemes divide the reactions (or species) into fast and slow ones. Then, the fast dynamics are resolved by making use of a quasiequilibrium assumption and the slow terms are integrated employing larger step sizes -see [13,18,26,27,28,30,37]. In this paper we do not assume that the system is clearly separable into fast and slow dynamics, therefore multirate methods are not discussed in what follows.…”

Computational Singular Perturbation Method and Tangential Stretching Rate Analysis of Large Scale Simulations of Reactive Flows: Feature Tracking, Time Scale Characterization, and Cause/Effect Identification. Part 2, Analyses of Ignition Systems, Laminar and Turbulent Flames

“…It must be noted that the CSP method provided the theoretical framework for designing solvers for deterministic sti problems (the CSP solver [48] and the G-Scheme solver [31,32]), for sti uncertain problems [49], for sti stochastic problems [50,51].…”