Calculating Effective Diffusivities in the Limit of Vanishing Molecular Diffusion

Abstract

In this paper we study the problem of the numerical calculation (by Monte Carlo Methods) of the effective diffusivity for a particle moving in a periodic divergent-free velocity filed, in the limit of vanishing molecular diffusion. In this limit traditional numerical methods typically fail, since they do not represent accurately the geometry of the underlying deterministic dynamics. We propose a stochastic splitting method that takes into account the volume preserving property of the equations motion in the absence of noise, and when inertial effects can be neglected. An extension of the method is then proposed for the cases where the noise has a non trivial time-correlation structure and when inertial effects cannot be neglected. Modified equations are used to perform backward error analysis. The new stochastic geometric integrators are shown to outperform standard Euler-based integrators. Various asymptotic limits of physical interest are investigated by means of numerical experiments, using the new integrators.

36 pages, 21 figures, removed the labeling from the equations

References in corpus (3)

• Particles and fields in fluid turbulence
• Periodic Homogenization for Inertial Particles
• Stochastic suspensions of heavy particles

Cited by in corpus (9)

• Long-Run Accuracy of Variational Integrators in the Stochastic Context
• Pathwise Accuracy and Ergodicity of Metropolized Integrators for SDEs
• Dispersion of inertial particles in cellular flows in the small-Stokes, large-Péclet regime
• Sharp error estimates on a stochastic structure-preserving scheme in computing effective diffusivity of 3D chaotic flows
• Homogenization for advection-diffusion in a perforated domain
• Strang Splitting for Parametric Inference in Second-order Stochastic Differential Equations
• On Metropolis Integrators for Molecular Dynamics
• Convergence of stochastic structure-preserving schemes for computing effective diffusivity in random flows
• Convergence analysis of a Lagrangian numerical scheme in computing effective diffusivity of 3D time-dependent flows