Quasi-Monte Carlo algorithms for diffusion equations in high dimensions


Dateien zu dieser Ressource

Dateien Größe Format Anzeige

Zu diesem Dokument gibt es keine Dateien.

VENKITESWARAN, Gopalakrishnan, Michael JUNK, 2005. Quasi-Monte Carlo algorithms for diffusion equations in high dimensions. In: Mathematics and Computers in Simulation. 68(1), pp. 23-41. ISSN 0378-4754. eISSN 1872-7166

@article{Venkiteswaran2005Quasi-25405, title={Quasi-Monte Carlo algorithms for diffusion equations in high dimensions}, year={2005}, doi={10.1016/j.matcom.2004.09.003}, number={1}, volume={68}, issn={0378-4754}, journal={Mathematics and Computers in Simulation}, pages={23--41}, author={Venkiteswaran, Gopalakrishnan and Junk, Michael} }

eng Venkiteswaran, Gopalakrishnan deposit-license 2013-12-13T10:56:08Z Venkiteswaran, Gopalakrishnan 2005 Quasi-Monte Carlo algorithms for diffusion equations in high dimensions Mathematics and Computers in Simulation ; 68 (2005), 1. - S. 23-41 2013-12-13T10:56:08Z Junk, Michael Diffusion equation posed on a high dimensional space may occur as a sub-problem in advection-diffusion problems (see [G. Venkiteswaran, M. Junk, A QMC approach for high dimensional Fokker–Planck equations modelling polymeric liquids, Math. Comput. Simul. 68 (2005) 43–56.] for a specific application). Although the transport part can be dealt with the method of characteristics, the efficient simulation of diffusion in high dimensions is a challenging task. The traditional Monte Carlo method (MC) applied to diffusion problems converges and is N<sup>−1/2</sup> accurate, where N is the number of particles. It is well known that for integration, quasi-Monte Carlo (QMC) outperforms Monte Carlo in the sense that one can achieve N<sup>−1</sup> convergence, up to a logarithmic factor. This is our starting point to develop methods based on Lécot’s approach [C. Lécot, F.E. Khettabi, Quasi-Monte Carlo simulation of diffusion, Journal of Complexity 15 (1999) 342–359.], which are applicable in high dimensions, with a hope to achieve better speed of convergence. Through a number of numerical experiments we observe that some of the QMC methods not only generalize to high dimensions but also show faster convergence in the results and thus, slightly outperform standard MC. Junk, Michael

Das Dokument erscheint in:

KOPS Suche


Mein Benutzerkonto