On the use of reverse Brownian motion to accelerate hybrid simulations

JOURNAL OF COMPUTATIONAL PHYSICS

卷:334

Article

Bakarji, Joseph¹  Tartakovsky, Daniel M.¹ 

[1] Stanford Univ, Dept Energy Resources Engn, 367 Panama St, Stanford, CA 94305 USA

关键词: Hybrid simulation;    Local Monte Carlo;    Reverse Brownian motion;   

DOI  :  10.1016/j.jcp.2016.12.032

来源: Elsevier

【 摘 要 】

Multiscale and multiphysics simulations are two rapidly developing fields of scientific computing. Efficient coupling of continuum (deterministic or stochastic) constitutive solvers with their discrete (stochastic, particle-based) counterparts is a common challenge in both kinds of simulations. We focus on interfacial, tightly coupled simulations of diffusion that combine continuum and particle -based solvers. The latter employs the reverse Brownian motion (rBm), a Monte Carlo approach that allows one to enforce inhomogeneous Dirichlet, Neumann, or Robin boundary conditions and is trivially parallelizable. We discuss numerical approaches for improving the accuracy of rBm in the presence of inhomogeneous Neumann boundary conditions and alternative strategies for coupling the rBm solver with its continuum counterpart. Numerical experiments are used to investigate the convergence, stability, and computational efficiency of the proposed hybrid algorithm. (C) 2016 Elsevier Inc. All rights reserved.

【 授权许可】

Free