Verification methodology for plasma simulations and application to a scrape-off layer turbulence code
Riva, Fabio Matteo
ORCID
Istituto ricerche solari Aldo e Cele Daccò (IRSOL), Faculty of Informatics, Università della Svizzera italiana Switzerland - Centre de Recherches en Physique des Plasmas (CRPP), École Polytechnique Fédérale de Lausanne (EPFL), Switzerland
Ricci, PaoloCentre de Recherches en Physique des Plasmas (CRPP), École Polytechnique Fédérale de Lausanne (EPFL), Switzerland
Halpern, Federico D.Centre de Recherches en Physique des Plasmas (CRPP), École Polytechnique Fédérale de Lausanne (EPFL), Switzerland
Jolliet, SébastienCentre de Recherches en Physique des Plasmas (CRPP), École Polytechnique Fédérale de Lausanne (EPFL), Switzerland
Loizu, JoaquimCentre de Recherches en Physique des Plasmas (CRPP), École Polytechnique Fédérale de Lausanne (EPFL), Switzerland
Mosetto, AnnamariaCentre de Recherches en Physique des Plasmas (CRPP), École Polytechnique Fédérale de Lausanne (EPFL), Switzerland
Physics of Plasmas. - 2014, vol. 21, no. 6, p. 062301
English
Bridging the gap between plasma physics and other scientific domains, in particular, the computational fluid dynamics community, a general, rigorous, and simple-to-apply methodology is presented for both the verification of the correct implementation of the model equations (code verification) and numerical error quantification (solution verification). The proposed code verification procedure consists in using the method of manufactured solutions and executing an order-of-accuracy test, assessing the rate of convergence of the numerical solution to the manufactured one. For the solution verification, the numerical error is quantified by applying the Richardson extrapolation, which provides an approximation of the analytical solution, and by using the grid convergence index to estimate the numerical uncertainty affecting the simulation results. The methodology is applied to verify the correct implementation of the drift-reduced Braginskii equations into the GBS code, and to estimate the numerical error affecting the GBS solutions. The GBS code is successfully verified, and an estimate of the numerical error affecting the simulation results is provided.