Volume 8, Issue 4
An Efficient Algorithm to Simulate a Brownian Motion Over Irregular Domains

S. Zein, A. Lejay & M. Deaconu

Commun. Comput. Phys., 8 (2010), pp. 901-916.

Published online: 2010-08

Preview Full PDF 95 553
Export citation
  • Abstract

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.

  • Keywords

  • AMS Subject Headings

  • Copyright

COPYRIGHT: © Global Science Press

  • Email address
  • References
  • Hide All
    View All

@Article{CiCP-8-901, author = {S. Zein, A. Lejay and M. Deaconu}, title = {An Efficient Algorithm to Simulate a Brownian Motion Over Irregular Domains}, journal = {Communications in Computational Physics}, year = {2010}, volume = {8}, number = {4}, pages = {901--916}, abstract = {

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.

}, issn = {1991-7120}, doi = {https://doi.org/10.4208/cicp.240209.031209a}, url = {http://global-sci.org/intro/article_detail/cicp/7601.html} }
Copy to clipboard
The citation has been copied to your clipboard