Randomized quasi-Monte Carlo simulation of fast-ion thermalization
2012 (English)In: Computational Science & Discovery, ISSN 1749-4680, Vol. 5, no 1, 014010- p.Article in journal (Refereed) Published
This work investigates the applicability of the randomized quasi-Monte Carlo method for simulation of fast-ion thermalization processes in fusion plasmas, e.g. for simulation of neutral beam injection and radio frequency heating. In contrast to the standard Monte Carlo method, the quasi-Monte Carlo method uses deterministic numbers instead of pseudo-random numbers and has a statistical weak convergence close to O(N -1), where N is the number of markers. We have compared different quasi-Monte Carlo methods for a neutral beam injection scenario, which is solved by many realizations of the associated stochastic differential equation, discretized with the Euler-Maruyama scheme. The statistical convergence of the methods is measured for time steps up to 2 14.
Place, publisher, year, edition, pages
Institute of Physics (IOP), 2012. Vol. 5, no 1, 014010- p.
Monte Carlo, quasi-Monte Carlo, low discrepancy sequence, variance reduction
Fusion, Plasma and Space Physics
IdentifiersURN: urn:nbn:se:kth:diva-91197DOI: 10.1088/1749-4699/5/1/014010ScopusID: 2-s2.0-84866342521OAI: oai:DiVA.org:kth-91197DiVA: diva2:508790
QC 20121114. Updated from submitted to published.2012-03-132012-03-092013-04-05Bibliographically approved