Library Subscription: Guest
Begell Digital Portal Begell Digital Library eBooks Journals References & Proceedings Research Collections
International Journal for Uncertainty Quantification
IF: 4.911 5-Year IF: 3.179 SJR: 1.008 SNIP: 0.983 CiteScore™: 5.2

ISSN Print: 2152-5080
ISSN Online: 2152-5099

Open Access

International Journal for Uncertainty Quantification

DOI: 10.1615/Int.J.UncertaintyQuantification.2015006492
pages 327-339

UNCERTAINTY QUANTIFICATION OF THE GEM CHALLENGE MAGNETIC RECONNECTION PROBLEM USING THE MULTILEVEL MONTE CARLO METHOD

Eder M. Sousa
Aerospace and Energetics Research Program, University of Washington, Seattle, Washington 98195, USA
Guang Lin
Computational Science & Mathematics Division, Pacific Northwest National Laboratory, Richland, Washington 99352; Department of Mathematics, School of Mechanical Engineering, Purdue University, West Lafayette, Indiana, USA
Uri Shumlak
Aerospace and Energetics Research Program, University of Washington, Seattle, Washington 98195, USA

ABSTRACT

Plasma modelers often change the ion-to-electron mass ratio and speed of light to Alfven speed ratio to decrease computational cost. Changing these parameters may affect the outcome of simulation and uncertainty in the results. This work aims to quantify the uncertainty of varying the ion-to-electron mass ratio, speed of light to Alfven speed ratio, and the initial magnetic flux perturbation on the reconnected flux to provide a confidence limit. In this study, the multilevel Monte Carlo (MMC) method is employed to estimate the mean and variance, and the results are compared with the standard Monte Carlo (MC) and the probabilistic collocation (PC) methods. The plasma model used here is the two-fluid plasma where ions and electrons are treated as two separate fluids. Numerical simulations are presented showing the effectiveness of the MMC method when applied to the quasi-neutral ion cyclotron waves and the Geospace Environment Modeling (GEM) magnetic reconnection challenge problems. The mean reconnected flux with error bars provides a reconnection flux variation envelope, which can help numerical modelers to evaluate whether their reconnection flux lies inside the envelope for different plasma models. The results of the MMC mean and variance are comparable to the MC method but at a much lower computational cost.


Articles with similar content:

UTILIZING ADJOINT-BASED ERROR ESTIMATES FOR SURROGATE MODELS TO ACCURATELY PREDICT PROBABILITIES OF EVENTS
International Journal for Uncertainty Quantification, Vol.8, 2018, issue 2
Timothy Wildey, Troy Butler
CONTROLLING LIGHT TRAPPING IN DISORDERED STRUCTURES WITH SHORT-RANGE ORDER
International Heat Transfer Conference 16, Vol.17, 2018, issue
Changying Zhao, X.W. Chen, B. X. Wang, M.Q. Liu
Numerical Assessment of Theoretical Error Estimates in Coarse-Grained Kinetic Monte Carlo Simulations: Application to Surface Diffusion
International Journal for Multiscale Computational Engineering, Vol.3, 2005, issue 1
Markos A. Katsoulakis, Dionisios G. Vlachos, Abhijit Chatterjee
GRADIENT-BASED STOCHASTIC OPTIMIZATION METHODS IN BAYESIAN EXPERIMENTAL DESIGN
International Journal for Uncertainty Quantification, Vol.4, 2014, issue 6
Youssef Marzouk, Xun Huan
EXPERIMENTAL STUDY ON EVAPORATION CHARACTERISTICS OF MENISCUS
International Heat Transfer Conference 16, Vol.4, 2018, issue
Shu Soma, Tomoaki Kunugi, Zensaku Kawara, Takehiko Yokomine