English

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

Computational Physics 2026-03-03 v3 Numerical Analysis Numerical Analysis Plasma Physics

Abstract

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 10410^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.

Keywords

Cite

@article{arxiv.2501.08482,
  title  = {Surrogate-based multilevel Monte Carlo methods for uncertainty quantification in the Grad-Shafranov free boundary problem},
  author = {Howard Elman and Jiaxing Liang and Tonatiuh Sánchez-Vizuet},
  journal= {arXiv preprint arXiv:2501.08482},
  year   = {2026}
}

Comments

In memory of Prof. Howard C. Elman, esteemed friend, mentor, and outstanding computational scientist

R2 v1 2026-06-28T21:06:37.274Z