Related papers: A hybrid simulation for a system of singularly per…
We develop symbolic methods of asymptotic approximations for solutions of linear ordinary differential equations and use to them stabilize numerical calculations. Our method follows classical analysis for first-order systems and…
This article describes algorithms for the hybrid parallelization and SIMD vectorization of molecular dynamics simulations with short-range forces. The parallelization method combines domain decomposition with a thread-based parallelization…
We present a new time-stepping algorithm for nonlinear PDEs that exhibit scale separation in time. Our scheme combines asymptotic techniques (which are inexpensive but can have insufficient accuracy) with parallel-in-time methods (which,…
A singularly perturbed parabolic problem of convection-diffusion type with a discontinuous initial condition is examined. An analytic function is identified which matches the discontinuity in the initial condition and also satisfies the…
In this paper, we develop a second-order accurate time-stepping scheme for the tempered time-fractional advection-dispersion equation based on a sum-of-exponentials (SOE) approximation to the convolution kernel involved in the fractional…
Shocks in astrophysical fluids can generate suprathermal particles by first order (or diffusive) Fermi acceleration. In the test particle regime there is a simple relation between the spectrum of the accelerated particles and the jump…
Several relaxation approximations to partial differential equations have been recently proposed. Examples include conservation laws, Hamilton-Jacobi equations, convection-diffusion problems, gas dynamics problems. The present paper focuses…
Different relaxation approximations to partial differential equations, including conservation laws, Hamilton-Jacobi equations, convection-diffusion problems, gas dynamics problems, have been recently proposed. The present paper focuses onto…
We present a test for the problem of decentralized sequential hypothesis testing, which is asymptotically optimum. By selecting a suitable sampling mechanism at each sensor, communication between sensors and fusion center is asynchronous…
The reaction-diffusion model can generate a wide variety of spatial patterns, which has been widely applied in chemistry, biology, and physics, even used to explain self-regulated pattern formation in the developing animal embryo. In this…
The multi-level method for discrete state systems, first introduced by Anderson and Higham [Multiscale Model. Simul. 10:146--179, 2012], is a highly efficient simulation technique that can be used to elucidate statistical characteristics of…
We show that the asymptotic iteration method converges and yields accurate energies for a perturbed Coulomb model. We also discuss alternative perturbation approaches to that model.
Stability and convergence of a time-weighted discrete scheme with nonuniform time steps are established for linear reaction-subdiffusion equations. The Caupto derivative is approximated at an offset point by using linear and quadratic…
A new method is presented for mesoscopic simulations of particle dispersions in nematic liquid crystal solvents. It allows efficient first-principle simulations of the dispersions involving many particles with many-body interactions…
The reversible reactions like A+B <-> C in the many-component diffusive system affect the diffusive properties of the constituents. The effective conjugation of irreversible processes of different dimensionality takes place due to the…
To achieve efficient and accurate long-time integration, we propose a fast, accurate, and stable high-order numerical method for solving fractional-in-space reaction-diffusion equations. The proposed method is explicit in nature and…
The shrinking core model describes the reaction of a spherical solid particle with a surrounding fluid. In this work, we revisit the SCM by deriving it from the underlying physical processes and performing a careful non-dimensionalisation,…
We consider the problem of inference for nonlinear, multivariate diffusion processes, satisfying It\^o stochastic differential equations (SDEs), using data at discrete times that may be incomplete and subject to measurement error. Our…
We introduce the concept of Hybrid Surrogate Models (HSMs) -- combining multivariate polynomials with Heavyside functions -- as approximates of functions with finitely many jump discontinuities. We exploit the HSMs for formulating a…
In this paper, we propose high order numerical methods to solve a 2D advection diffusion equation, in the highly oscillatory regime. We use an integrator strategy that allows the construction of arbitrary high-order schemes {leading} to an…