Surrogate-based multilevel Monte Carlo methods for uncertainty quantification in the Grad-Shafranov free boundary problem

Howard Elman, Jiaxing Liang, Tonatiuh Sánchez-Vizuet

arXiv:2501.08482·physics.comp-ph·Published 2025-01-14·Updated 2026-03-01

We explore a hybrid technique to quantify the variability in the numerical solutions to a free boundary problem associated with magnetic equilibrium in axisymmetric fusion reactors amidst parameter uncertainties. The method aims at reducing computational costs by integrating a surrogate model into a multilevel Monte Carlo method. The resulting surrogate-enhanced multilevel Monte Carlo methods reduce the cost of simulation by factors as large as $10^4$ compared to standard Monte Carlo simulations involving direct numerical solutions of the associated Grad-Shafranov partial differential equation. Accuracy assessments also show that surrogate-based sampling closely aligns with the results of direct computation, confirming its effectiveness in capturing the behavior of plasma boundary and geometric descriptors.

TopicsProcess Modeling & System Identification, Uncertainty Quantification & Bayesian Methods

Tagssurrogate-modeling

arXiv categoriesphysics.comp-ph, math.NA, physics.plasm-ph

arXiv abstract pagePDF