Simulating diffusion processes in discontinuous media: Benchmark tests

JOURNAL OF COMPUTATIONAL PHYSICS

卷:314

Article

Lejay, Antoine^1,2,3,5  Pichot, Geraldine⁴ 

[1] Univ Lorraine, IECL, UMR 7502, F-54500 Vandoeuvre Les Nancy, France

[2] CNRS, IECL, UMR 7502, F-54500 Vandoeuvre Les Nancy, France

[3] Inria, F-54600 Villers Les Nancy, France

[4] Inria Paris, 2 Rue Simone Iff,CS 42112, F-75589 Paris 12, France

[5] IECL, BP 70238, F-54506 Vandoeuvre Les Nancy, France

关键词: Monte Carlo methods for discontinuous media;    Fick's law;    Breakthrough curve;   

DOI  :  10.1016/j.jcp.2016.03.003

来源: Elsevier

【 摘 要 】

We present several benchmark tests for Monte Carlo methods simulating diffusion in one-dimensional discontinuous media. These benchmark tests aim at studying the potential bias of the schemes and their impact on the estimation of micro- or macroscopic quantities (repartition of masses, fluxes, mean residence time,...). These benchmark tests are backed by a statistical analysis to filter out the bias from the unavoidable Monte Carlo error. We apply them on four different algorithms. The results of the numerical tests give a valuable insight into the fine behavior of these schemes, as well as rules to choose between them. (C) 2016 Elsevier Inc. All rights reserved.

【 授权许可】

Free