Open RemDelaporteMathurin opened 2 years ago
We have some of the inputs over here https://github.com/fusion-energy/paramak/blob/main/examples/example_parametric_components/make_plasmas.ipynb
Ah yes good point. These are the geometrical properties but we also need some more like "ion temperature at the pedestal", "ion density at the separatrix", etc.
@RemDelaporteMathurin @Shimwell Hi. From https://www.sciencedirect.com/science/article/pii/S0920379620301514 I could gather the following parameters for the EU-DEMO 2018:
From Table 1:
Taking information from figures may be a bit subjective, but for example from Figure 3 (or 1) we could set the pedestal radius at around 0.86 so pedestal_radius=0.86 * 9.0/3.1 I believe temperatures and densities (at core or pedestal) could be taken from Figure 1 too.
Hi @tokamaster
Thanks for taking the time to gather the info. and thanks for the ref
I think DEMO's minor radius is 3.1m (9.0/3.1 would be way too large 😄 ).
I think the parameters in the example were taken for DEMO already but this paper will maybe refine them :-)
my_demo_plasma = TokamakSource(
elongation=1.557,
ion_density_centre=1.09e20,
ion_density_peaking_factor=1,
ion_density_pedestal=1.09e20,
ion_density_separatrix=3e19,
ion_temperature_centre=45.9,
ion_temperature_peaking_factor=8.06,
ion_temperature_pedestal=6.09,
ion_temperature_separatrix=0.1,
major_radius=9.0,
minor_radius=3.1,
pedestal_radius=0.8 * 2.92258,
mode="H",
shafranov_factor=0.44789,
triangularity=0.270,
ion_temperature_beta=6
)
I'm thinking we could have a new class, something like EuDemo(TokamakSource)
which inherits from TokamakSource
and defines defaults
@RemDelaporteMathurin If I'm not mistaken 3.1 is the aspect ratio, so minor radius would be 9.0/3.1 = 2.90322580645 m.
@RemDelaporteMathurin If I'm not mistaken 3.1 is the aspect ratio, so minor radius would be 9.0/3.1 = 2.90322580645 m.
@tokamaster Yes you're right, my bad, I probably hadn't had my coffee yet.... 😄
As @pshriwise suggested, it would be more convenient for users to have device specific plasma classes.
I can think of DEMO, ITER, MAST-U, SPARC (?)
Then again some parameters for these machines can change, but I guess the geometrical parameters can be fixed.
@Shimwell do you know where we could find the info for these devices?