Device specific plasma classes

[RemDelaporteMathurin]

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?

[shimwell]

We have some of the inputs over here
https://github.com/fusion-energy/paramak/blob/main/examples/example_parametric_components/make_plasmas.ipynb

[RemDelaporteMathurin]

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.

[tokamaster]

@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:

• Major radius: 9.0
• Minor radius: 9.0/3.1
• Elongation (at 95%): 1.65
• Triangularity (at 95%): 0.33

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.

[RemDelaporteMathurin]

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
    )

[shimwell]

I'm thinking we could have a new class, something like EuDemo(TokamakSource) which inherits from TokamakSource and defines defaults

[tokamaster]

@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]

@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.... 😄