Multilevel Monte Carlo methods for the Grad-Shafranov free boundary problem
dc.contributor.author | Elman, Howard C. | |
dc.contributor.author | Liang, Jiaxing | |
dc.contributor.author | Sánchez-Vizuet, Tonatiuh | |
dc.date.accessioned | 2024-02-09T20:10:03Z | |
dc.date.available | 2024-02-09T20:10:03Z | |
dc.date.issued | 2024-01-19 | |
dc.identifier.citation | Elman, H. C., Liang, J., & Sánchez-Vizuet, T. (2024). Multilevel Monte Carlo methods for the Grad-Shafranov free boundary problem. Computer Physics Communications, 109099. | en_US |
dc.identifier.issn | 0010-4655 | |
dc.identifier.doi | 10.1016/j.cpc.2024.109099 | |
dc.identifier.uri | http://hdl.handle.net/10150/670933 | |
dc.description.abstract | The equilibrium configuration of a plasma in an axially symmetric reactor is described mathematically by a free boundary problem associated with the celebrated Grad-Shafranov equation. The presence of uncertainty in the model parameters introduces the need to quantify the variability in the predictions. This is often done by computing a large number of model solutions on a computational grid for an ensemble of parameter values and then obtaining estimates for the statistical properties of solutions. In this study, we explore the savings that can be obtained using multilevel Monte Carlo methods, which reduce costs by performing the bulk of the computations on a sequence of spatial grids that are coarser than the one that would typically be used for a simple Monte Carlo simulation. We examine this approach using both a set of uniformly refined grids and a set of adaptively refined grids guided by a discrete error estimator. Numerical experiments show that multilevel methods dramatically reduce the cost of simulation, with cost reductions typically on the order of 60 or more and possibly as large as 200. Adaptive griding results in more accurate computation of geometric quantities such as x-points associated with the model. | en_US |
dc.description.sponsorship | U.S. Department of Energy | en_US |
dc.language.iso | en | en_US |
dc.publisher | Elsevier BV | en_US |
dc.rights | © 2024 Elsevier B.V. All rights reserved. | en_US |
dc.rights.uri | https://rightsstatements.org/vocab/InC/1.0/ | en_US |
dc.subject | General Physics and Astronomy | en_US |
dc.subject | Hardware and Architecture | en_US |
dc.subject | Adaptive finite element discretization | en_US |
dc.subject | Free boundary Grad-Shafranov problem | en_US |
dc.subject | Multilevel Monte Carlo Finite-Element | en_US |
dc.subject | Uncertainty quantification | en_US |
dc.title | Multilevel Monte Carlo methods for the Grad-Shafranov free boundary problem | en_US |
dc.type | Article | en_US |
dc.contributor.department | Department of Mathematics, The University of Arizona | en_US |
dc.identifier.journal | Computer Physics Communications | en_US |
dc.description.note | 24 month embargo; first published 19 January 2024 | en_US |
dc.description.collectioninformation | This item from the UA Faculty Publications collection is made available by the University of Arizona with support from the University of Arizona Libraries. If you have questions, please contact us at repository@u.library.arizona.edu. | en_US |
dc.eprint.version | Final accepted manuscript | en_US |
dc.identifier.pii | S0010465524000225 | |
dc.source.journaltitle | Computer Physics Communications | |
dc.source.volume | 298 | |
dc.source.beginpage | 109099 |