An Efficient Algorithm to Simulate a Brownian Motion Over Irregular Domains  

作　　者：S.Zein A.Lejay M.Deaconu 

机构地区：[1]INRIA Centre de Recherche Nancy-Grand Est-IECN,France

出　　处：《Communications in Computational Physics》2010年第9期901-916,共16页计算物理通讯（英文）

摘　　要：In this paper,we present an algorithm to simulate a Brownian motion by coupling two numerical schemes:the Euler scheme with the random walk on the hyper-rectangles.This coupling algorithm has the advantage to be able to compute the exit time and the exit position of a Brownian motion from an irregular bounded domain(with corners at the boundary),and being of order one with respect to the time step of the Euler scheme.The efficiency of the algorithm is studied through some numerical examples by comparing the analytical solution with the Monte Carlo solution of some Poisson problems.The Monte Carlo solution of these PDEs requires simulating Brownian motions of different types(natural,reflected or drifted)over an irregular domain.

关 键 词：Brownian motion Monte Carlo methods partial differential equations Euler scheme random walk on rectangles 

分 类 号：O17[理学—数学]