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: 4d32f454cc4ea33181b389ff8544b739bceeacab], 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, 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/elongation/output_Ip10.0_tau4.00_Tfin1_s0.00_fb1.00.h5" - "Data/elongation/output_Ip10.0_tau4.00_Tfin1_s0.00_fb1.00_circ.h5" - "Data/elongation/output_Ip10.0_tau4.00_Tfin1_s0.00_fb1.00_elong.h5"- Use method "runSimulation" in "runSimulation.py" to produce data with inputs - Ip=10e6 - tau=4e-3 - Tfin=1 - 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"