Fluid and kinetic studies of tokamak disruptions using Bayesian optimization

<p>The codes and the data in this directory corresponds to the code and results used in the paper [I. Ekmark et al (2024) J. Plasma Phys., Fluid and kinetic studies of tokamak disruptions using Bayesian optimization, http://arxiv.org/abs/2402.05843]. References to figures below refer to this publication.&nbsp;</p> <p>The optimizations have been performed using the Python package by Fernando Nogueira [https://github.com/bayesian-optimization/BayesianOptimization] and the simulations are performed using the disruptions simulation simulation tool DREAM [https://github.com/chalmersplasmatheory/DREAM, git hash: 0d786e859f6228185ef68b6b3639747e8d96172d], for more information on the latter code visit https://ft.nephy.chalmers.se/dream/.</p> <p>The codes:<br>&nbsp;- BayesianOptimization.py: Runs the optimization, first in fluid and then in isotropic mode. For activated simulations, use the flag "-A".<br>&nbsp; &nbsp; &nbsp;- BlackBox.py: Contains the functions that are run in the optimizations and sets up the simulations.<br>&nbsp; &nbsp; &nbsp;- utils.py: Contains the settings for the simulations as well as some other functions needed in BlackBox.py<br>&nbsp; &nbsp; &nbsp;- ITER.py: Contains all the ITER specific settings.<br>&nbsp; &nbsp; &nbsp;- Exceptions.py: Contains exceptions needed during the simulations.&nbsp;<br>&nbsp; &nbsp; &nbsp;- CostFunction.py: Contains the functions used for evaluating the cost function value for specified values of the representative runaway current, final Ohmic current, current quench time and transported heat fraction.<br>&nbsp;- RunCases.py: Sets up simulations for the cases of table 1 in the paper, as well as for all the optima found.</p> <p>The data:<br>&nbsp;- Optimization results:<br>&nbsp;<br>&nbsp; &nbsp; - Data/OptimizationResults/optresult_fluid.json: Contains the optimization data for the non-activated case using the fluid model. Used to produce figure 1.a.&nbsp;<br>&nbsp; &nbsp; - Data/OptimizationResults/optresult_isotropic.json: Contains the optimization data for the non-activated case using the isotropic model. Used to produce figure 1.b.&nbsp;<br>&nbsp; &nbsp; - Data/OptimizationResults/optresult_fluid_activated.json: Contains the optimization data for the activated case using the fluid model. Used to produce figure 5.a.<br>&nbsp; &nbsp; - Data/OptimizationResults/optresult_isotropic_activated.json: Contains the optimization data for the activated case using the isotropic model. Used to produce figure 5.b.<br>&nbsp; &nbsp;&nbsp;<br>&nbsp; &nbsp; - Data/OptimizationResults/components_fluid.json: Contains the cost function components for each sample from the optimization of the non-activated case using the fluid model. Used to produce figure 2.a.&nbsp;<br>&nbsp; &nbsp; - Data/OptimizationResults/components_isotropic.json: Contains the cost function components for each sample from the optimization of the non-activated case using the isotropic model. Used to produce figure 2.b.&nbsp;<br>&nbsp; &nbsp; - Data/OptimizationResults/components_fluid_activated.json: Contains the cost function components for each sample from the optimization of the activated case using the fluid model. Used to produce figure 6.a.<br>&nbsp; &nbsp; - Data/OptimizationResults/components_isotropic_activated.json: Contains the cost function components for each sample from the optimization of the activated case using the isotropic model. Used to produce figure 6.b.<br>&nbsp; &nbsp;&nbsp;<br>&nbsp;- Cases:<br>&nbsp; &nbsp; - Data/Cases/nonActivatedOpts/fluidOpt/: Contains outputfiles for fluid and isotropic simulations of the optimal case found for the non-activated scenario using the fluid model.<br>&nbsp; &nbsp; - Data/Cases/nonActivatedOpts/isoOpt/: Contains outputfiles for fluid and isotropic simulations of the optimal case found for the non-activated scenario using the isotropic model.<br>&nbsp; &nbsp; - Data/Cases/activatedOpts/fluidOpt/: Contains outputfiles for fluid and isotropic simulations of the optimal case found for the activated scenario using the fluid model.<br>&nbsp; &nbsp; - Data/Cases/activatedOpts/isoOpt/: Contains outputfiles for fluid and isotropic simulations of the optimal case found for the activated scenario using the isotropic model.<br>&nbsp; &nbsp; - Data/Cases/[circle, cross, square, triangle]: Contains outputfiles for fluid and isotropic simulations corresponding to the cases presented in table 1.</p>