Global Sensitivity Analysis on Numerical Solver Parameters of Particle-In-Cell Models in Particle Accelerator Systems

Abstract

Every computer model depends on numerical input parameters that are chosen according to mostly conservative but rigorous numerical or empirical estimates. These parameters could for example be the step size for time integrators, a seed for pseudo-random number generators, a threshold or the number of grid points to discretize a computational domain. In case a numerical model is enhanced with new algorithms and modelling techniques the numerical influence on the quantities of interest, the running time as well as the accuracy is often initially unknown. Usually parameters are chosen on a trial-and-error basis neglecting the computational cost versus accuracy aspects. As a consequence the cost per simulation might be unnecessarily high which wastes computing resources. Hence, it is essential to identify the most critical numerical parameters and to analyze systematically their effect on the result in order to minimize the time-to-solution without losing significantly on accuracy. Relevant parameters are identified by global sensitivity studies where Sobol' indices are common measures. These sensitivities are obtained by surrogate models based on polynomial chaos expansion. In this paper, we first introduce the general methods for uncertainty quantification. We then demonstrate their use on numerical solver parameters to reduce the computational costs and discuss further model improvements based on the sensitivity analysis. The sensitivities are evaluated for neighbouring bunch simulations of the existing PSI Injector II and PSI Ring as well as the proposed Daedalus Injector cyclotron and simulations of the rf electron gun of the Argonne Wakefield Accelerator.

References in corpus (10)

• Scikit-learn: Machine Learning in Python
• API design for machine learning software: experiences from the scikit-learn project
• An active-learning algorithm that combines sparse polynomial chaos expansions and bootstrap for structural reliability analysis
• Polynomial chaos expansions for dependent random variables
• Beam Dynamics in High Intensity Cyclotrons Including Neighboring Bunch Effects: Model, Implementation and Application
• OPAL a Versatile Tool for Charged Particle Accelerator Simulations
• Intensity limits of the PSI Injector II cyclotron
• High Current H2+ Cyclotrons for Neutrino Physics: The IsoDAR and DAE δALUS Projects
• Matching of turn pattern measurements for cyclotrons using multi-objective optimization
• Modelling challenges of the high power cyclotrons for the DAE$\delta$ALUS project

Cited by in corpus (2)

• Order-of-Magnitude Beam Current Improvement in Compact Cyclotrons
• Computational Models for High-Power Cyclotrons and FFAs