An Adaptive delta f Monte Carlo Method
2010 (English)In: IEEE Transactions on Plasma Science, ISSN 0093-3813, E-ISSN 1939-9375, Vol. 38, no 9, 2190-2197 p.Article in journal (Refereed) Published
A new adaptive delta f Monte Carlo method is presented with an application to radio frequency heating and transport in fusion plasmas. The method is suitable when an initial zeroth-order approximation of the distribution function is known. The difference between our method and earlier delta f methods is that we model the source term, obtained from the delta f ansatz, by adding particles. The rate of particles is defined by the inhomogeneous term in the Fokker-Planck equation. We develop an adaptive scheme for modifying the unperturbed part G(x) such that the number of particles used in the simulation for a fixed weight is minimized. This implicitly reduces the variance and improves computational efficiency. The method is tested on a one-dimensional Fokker-Planck model for RF-heating and compared against the analytical stationary solution.
Place, publisher, year, edition, pages
2010. Vol. 38, no 9, 2190-2197 p.
Control-variate, Monte Carlo (MC), variance reduction, delta f method
Fusion, Plasma and Space Physics
IdentifiersURN: urn:nbn:se:kth:diva-27685DOI: 10.1109/TPS.2010.2051686ISI: 000283252500015ScopusID: 2-s2.0-77956619121OAI: oai:DiVA.org:kth-27685DiVA: diva2:380409
QC 201012212010-12-212010-12-202013-04-05Bibliographically approved