Runaway Electrons in Stellarators: Unlikely or Unavoidable?

The codes and the data in this directory corresponds to the code and results used in the paper [Runaway Electrons in Stellarators: Unlikely or Unavoidable?, https://arxiv.org/abs/2607.29523]. References to figures below refer to this publication. The simulations are performed using the disruptions simulation simulation tool DREAM [https://github.com/chalmersplasmatheory/DREAM, git hash: 8a0f4ca4aeb4ed89f9fdc834a73141448297dc3f], for more information visit https://ft.nephy.chalmers.se/dream/. The DESC equilibrium solver python module is needed, for more information visit https://desc-docs.readthedocs.io/en/stable/index.html. The stellarator equilibrium files, provided by Stefan Buller, and needs to be downloaded and can be found at https://sb0095.mycpanel.princeton.edu/QA/notes.html The codes: ITER.py: Settings for ITER-like configuration runSimulation.py: For running simulations scan.py: For running scans, calls runSimulation.py stellarator.py: For stellarator, circular or elongated configuration Figure 4: Data in "Data/Scans/IpTend" folder Use method "scanIpAndTfin" in "scan.py" to produce data Figure 5: Data in "Data/Scans/Tend" folder Use method "scanTfin" in "scan.py" to produce data Figure 6: Data in "Data/Scans/Tau" folder Use method "scanTau" in "scan.py" to produce data Figure 7: Outputfiles "Data/avalancheMetric/output_Ip10.0_A5.h5" "Data/avalancheMetric/output_Ip10.0_circ.h5" "Data/avalancheMetric/output_Ip10.0_elong.h5" Use method "runSimulation" in "runSimulation.py" to produce data with inputs Ip=10e6 seed=True stellarator=True/False elongated=True/False fb=1 Figure 8: Data in "Data/Scans/fb" folder Use method "scanFb" in "scan.py" to produce data Figure 9 and 10: Data in "Data/Scans/profile" folder Use method "scanProfile" in "scan.py" to produce data, with "constantIp=False" for constant maximum magnitude of the current density Figure 11: Outputfiles in "Data/avalancheMetric/" Use method "runSimulation" in "runSimulation.py" to produce data with inputs seed=True stellarator=True/False elongated=True/False fb=1 Choose plasma current with input "Ip" If you want to run ITER-like, set "runITER=True"