| Type de publication : |
 |
Documents sans référence de publication (Preprint) |
|
| Domaine : |
|
|
|
| Titre : |
|
Simulating diffusion processes in discontinuous media: a numerical scheme with constant time steps |
|
| Auteur(s) : |
|
Antoine Lejay ( ,  ) 1, 2, Géraldine Pichot () 3 |
|
| Projet(s) / laboratoire(s) : |
|
|
|
| Équipe de recherche (sauf Inria) : |
|
Probabilités et statistiques |
| Résumé : |
|
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. |
|
| Résumé français : |
|
Dans cet article, nous proposons une nouvelle méthode de Monte Carlo pour déplacer des particules dans un milieu discontinu. Nous caractérisons le processus stochastique spécifiant la dynamique des particules. Le point clef est la réduction de ce problème à l'utilisation d'un mouvement brownien biaisé. Dans une zone où les coefficients sont localement constants de chaque côté d'une discontinuité, la nouvelle position de la particule après un pas de temps constant est tiré selon la loi exacte du mouvement brownien biaisé au temps considéré. Pour cela, nous proposons deux algorithmes : l'un reposant sur une simuation à deux pas avec un arrêt à l'interface, et l'autre reposant sur une simulation directe. Des cas tests illustrent leurs efficacités. |
|
| Classification ACM : |
|
| I.: Computing Methodologies/I.6: SIMULATION AND MODELING/I.6.8: Types of Simulation/I.6.8.6: Monte Carlo |
| G.: Mathematics of Computing/G.3: PROBABILITY AND STATISTICS/G.3.7: Probabilistic algorithms (including Monte Carlo) |
| G.3.15: Stochastic processes |
|
|
| Classification autre : |
|
(AMS 2010) 60J60 ; 65C05 ; 76S05 |
|
| Langue du document : |
|
Anglais |
|
|
| Mots-clés : |
|
divergence form operators – stochastic differential equation – skew Brownian motion – Monte Carlo simulation – Euler scheme – geophysics – diffusive media with interfaces |
|
| Date de rédaction : |
|
02/12/2011 |
|
| Projet ANR : |
|
| Référence du projet |
ANR-10-COSI-SIMUDMRI |
| Année |
2010 |
| Acronyme du projet |
SIMUDMRI |
| Titre du projet |
Simulation du signal d'IRM diffusion dans tissus biologiques |
|
|
| Collaboration(s) : |
|
TOSCA ; SAGE |
|
|
BibTeX,EndNote,...
