Abstract
We consider boundary-value problems for elliptic equations with constant coefficients and apply Monte Carlo methods to solving these equations. To take into account boundary conditions involving solution’s normal derivative, we apply the new mean-value relation written down at boundary point. This integral relation is exact and provides a possibility to get rid of the bias caused by usually used finite-difference approximation. We consider Neumann and mixed boundary-value problems, and also the problem with continuity boundary conditions, which involve fluxes. Randomization of the mean-value relation makes it possible to continue simulating walk-on-spheres trajectory after it hits the boundary. We prove the convergence of the algorithm and determine its rate. In conclusion, we present the results of some model computations.
Supported in part by NATO Linkage Grant and Grant for Leading Scientific Schools.
This is a preview of subscription content, log in via an institution to check access.
Access this chapter
Tax calculation will be finalised at checkout
Purchases are for personal use only
Preview
Unable to display preview. Download preview PDF.
References
Müller, M.E.: Some continuous Monte Carlo methods for the Dirichlet problem. Ann. Math. Statistics 27(3), 569–589 (1956)
Ermakov, S.M., Mikhailov, G.A.: Statisticheskoe modelirovanie [Statistical simulation] (in Russian). Nauka, Moscow (1982)
Elepov, B.S., et al.: Reshenie kraevyh zadach metodom Monte-Carlo [Solution of boundary value problems by the Monte Carlo method] (in Russian). Nauka, Novosibirsk (1980)
Sabelfeld, K.K., Simonov, N.A.: Random Walks on Boundary for solving PDEs. VSP, Utrecht (1994)
Haji-Sheikh, A., Sparrow, E.M.: The floating random walk and its application to Monte Carlo solutions of heat equations. SIAM J. Appl. Math. 14(2), 570–589 (1966)
Kronberg, A.A.: On algorithms for statistical simulation of the solution of boundary value problems of elliptic type (in Russian). Zh. Vychisl. Mat. i Mat. Phyz. 84(10), 1531–1537 (1984)
Mikhailov, G.A., Makarov, R.N.: Solution of boundary value problems of the second and third kind by Monte Carlo methods (in Russian). Sibirsk. Mat. Zh. 38(3), 603–614 (1997)
Mascagni, M., Simonov, N.A.: Monte Carlo methods for calculating some physical properties of large molecules. SIAM J. Sci. Comp. 26(1), 339–357 (2004)
Karaivanova, A., Mascagni, M., Simonov, N.A.: Solving BVPs using quasirandom walks on the boundary. In: Lirkov, I., et al. (eds.) LSSC 2003. LNCS, vol. 2907, pp. 162–169. Springer, Heidelberg (2004)
Karaivanova, A., Mascagni, M., Simonov, N.A.: Parallel quasi-random walks on the boundary. Monte Carlo Methods and Applications 10(3-4), 311–320 (2004)
Miranda, C.: Partial differential equations of elliptic type. Springer, New York (1970)
Simonov, N.A.: A random walk algorithm for solving boundary value problems with partition into subdomains (in Russian). In: Metody i algoritmy statisticheskogo modelirovanija [Methods and algorithms for statistical modelling], pp. 48–58. Akad. Nauk SSSR Sibirsk. Otdel., Vychisl. Tsentr, Novosibirsk (1983)
Author information
Authors and Affiliations
Editor information
Rights and permissions
Copyright information
© 2007 Springer Berlin Heidelberg
About this paper
Cite this paper
Simonov, N.A. (2007). Random Walks for Solving Boundary-Value Problems with Flux Conditions. In: Boyanov, T., Dimova, S., Georgiev, K., Nikolov, G. (eds) Numerical Methods and Applications. NMA 2006. Lecture Notes in Computer Science, vol 4310. Springer, Berlin, Heidelberg. https://doi.org/10.1007/978-3-540-70942-8_21
Download citation
DOI: https://doi.org/10.1007/978-3-540-70942-8_21
Publisher Name: Springer, Berlin, Heidelberg
Print ISBN: 978-3-540-70940-4
Online ISBN: 978-3-540-70942-8
eBook Packages: Computer ScienceComputer Science (R0)