4422 articles – 2353 Notices  [english version]
HAL : hal-00649170, version 2

Voir la fiche détaillée  BibTeX,EndNote,...
Versions disponibles
Simulating diffusion processes in discontinuous media: a numerical scheme with constant time steps
Antoine Lejay 1, 2, Géraldine Pichot 3
TOSCA ; SAGE Collaboration(s)
(02/12/2011)

In this article, we propose new Monte Carlo techniques for moving a diffusive particle in a discontinuous media. In this framework, we characterize the stochastic process that governs the positions of the particle. The key tool is the reduction of the process to a Skew Brownian Motion (SBM). In a zone where the coefficients are locally constant on each side of the discontinuity, the new position of the particle after a constant time step is sampled from the exact distribution of the SBM process at the considered time. To do so, we propose two different but equivalent algorithms: a two-steps simulation with a stop at the discontinuity and a one-step direct simulation of the SBM dynamic. Some benchmark tests illustrate their effectiveness.
1 :  Institut Elie Cartan Nancy (IECN)
CNRS : UMR7502 – INRIA – Université Henri Poincaré - Nancy I – Université Nancy II – Institut National Polytechnique de Lorraine (INPL)
2 :  TOSCA (INRIA Sophia Antipolis / INRIA Lorraine / IECN)
INRIA – CNRS : UMR7502 – Université Henri Poincaré - Nancy I – Université Nancy II – Institut National Polytechnique de Lorraine (INPL)
3 :  SAGE (INRIA - IRISA)
CNRS : UMR6074 – INRIA – Université de Rennes 1
Probabilités et statistiques
Mathématiques/Probabilités

Physique/Physique/Géophysique

Planète et Univers/Sciences de la Terre/Géophysique

Sciences de l'environnement/Milieux et Changements globaux
divergence form operators – stochastic differential equation – skew Brownian motion – Monte Carlo simulation – Euler scheme – geophysics – diffusive media with interfaces
Liste des fichiers attachés à ce document :
PDF
lejay_pichot_2011_v2.pdf(708.1 KB)