Skip to content

exhaust

PlasmaExhaust

Class to hold plasma exhaust calculations for plasma processing.

Source code in process/models/physics/exhaust.py 
  8
  9
 10
 11
 12
 13
 14
 15
 16
 17
 18
 19
 20
 21
 22
 23
 24
 25
 26
 27
 28
 29
 30
 31
 32
 33
 34
 35
 36
 37
 38
 39
 40
 41
 42
 43
 44
 45
 46
 47
 48
 49
 50
 51
 52
 53
 54
 55
 56
 57
 58
 59
 60
 61
 62
 63
 64
 65
 66
 67
 68
 69
 70
 71
 72
 73
 74
 75
 76
 77
 78
 79
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
class PlasmaExhaust:
    """Class to hold plasma exhaust calculations for plasma processing."""

    def __init__(self):
        self.outfile = constants.NOUT
        self.mfile = constants.MFILE

    @staticmethod
    def calculate_separatrix_power(
        f_p_alpha_plasma_deposited: float,
        p_alpha_total_mw: float,
        p_non_alpha_charged_mw: float,
        p_hcd_injected_total_mw: float,
        p_plasma_ohmic_mw: float,
        p_plasma_rad_mw: float,
    ) -> float:
        """
        Calculate the power crossing the separatrix (P_sep).

        Parameters
        ----------
        f_p_alpha_plasma_deposited : float
            Fraction of alpha power deposited in plasma.
        p_alpha_total_mw : float
            Total alpha power produced (MW).
        p_non_alpha_charged_mw : float
            Power from non-alpha charged particles (MW).
        p_hcd_injected_total_mw : float
            Total power injected by heating and current drive (MW).
        p_plasma_ohmic_mw : float
            Ohmic heating power (MW).
        p_plasma_rad_mw : float
            Radiated power from plasma (MW).

        Returns
        -------
        float
            Power crossing the separatrix (MW).
        """

        return (
            f_p_alpha_plasma_deposited * p_alpha_total_mw
            + p_non_alpha_charged_mw
            + p_hcd_injected_total_mw
            + p_plasma_ohmic_mw
            - p_plasma_rad_mw
        )

    @staticmethod
    def calculate_psep_over_r_metric(
        p_plasma_separatrix_mw: float, rmajor: float
    ) -> float:
        """
        Calculate the power crossing the separatrix per unit major radius (P_sep/R).

        Parameters
        ----------
        p_plasma_separatrix_mw : float
            Power crossing the separatrix (MW).
        rmajor : float
            Plasma major radius (m).

        Returns
        -------
        float
            Power crossing the separatrix per unit major radius (MW/m).
        """
        return p_plasma_separatrix_mw / rmajor

    @staticmethod
    def calculate_eu_demo_re_attachment_metric(
        p_plasma_separatrix_mw: float,
        b_plasma_toroidal_on_axis: float,
        q95: float,
        aspect: float,
        rmajor: float,
    ) -> float:
        """Calculate the EU DEMO divertor protection re-attachment metric for plasma exhaust.

        Parameters
        ----------
        p_plasma_separatrix_mw : float
            Power crossing the separatrix (MW).
        b_plasma_toroidal_on_axis : float
            Toroidal magnetic field on axis (T).
        q95 : float
            Safety factor at 95% flux surface.
        aspect : float
            Aspect ratio of the plasma.
        rmajor : float
            Plasma major radius (m).

        Returns
        -------
        float
            EU DEMO re-attachment metric (MW T /m).

        References
        ----------
        - M. Siccinio, G. Federici, R. Kembleton, H. Lux, F. Maviglia, and J. Morris,
          "Figure of merit for divertor protection in the preliminary design of the EU-DEMO reactor,"
          Nuclear Fusion, vol. 59, no. 10, pp. 106026-106026, Jul. 2019,
          doi: https://doi.org/10.1088/1741-4326/ab3153.

        - H. Zohm et al.,
          "A stepladder approach to a tokamak fusion power plant,"
          Nuclear Fusion, vol. 57, no. 8, pp. 086002-086002, May 2017,
          doi: https://doi.org/10.1088/1741-4326/aa739e.
        """

        return (p_plasma_separatrix_mw * b_plasma_toroidal_on_axis) / (
            q95 * aspect * rmajor
        )

outfile = constants.NOUT instance-attribute

mfile = constants.MFILE instance-attribute

calculate_separatrix_power(f_p_alpha_plasma_deposited, p_alpha_total_mw, p_non_alpha_charged_mw, p_hcd_injected_total_mw, p_plasma_ohmic_mw, p_plasma_rad_mw) staticmethod

Calculate the power crossing the separatrix (P_sep).

Parameters:

Name Type Description Default
f_p_alpha_plasma_deposited float

Fraction of alpha power deposited in plasma.

 required
p_alpha_total_mw float

Total alpha power produced (MW).

 required
p_non_alpha_charged_mw float

Power from non-alpha charged particles (MW).

 required
p_hcd_injected_total_mw float

Total power injected by heating and current drive (MW).

 required
p_plasma_ohmic_mw float

Ohmic heating power (MW).

 required
p_plasma_rad_mw float

Radiated power from plasma (MW).

 required

Returns:

Type Description
float

Power crossing the separatrix (MW).
Source code in process/models/physics/exhaust.py 
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
@staticmethod
def calculate_separatrix_power(
    f_p_alpha_plasma_deposited: float,
    p_alpha_total_mw: float,
    p_non_alpha_charged_mw: float,
    p_hcd_injected_total_mw: float,
    p_plasma_ohmic_mw: float,
    p_plasma_rad_mw: float,
) -> float:
    """
    Calculate the power crossing the separatrix (P_sep).

    Parameters
    ----------
    f_p_alpha_plasma_deposited : float
        Fraction of alpha power deposited in plasma.
    p_alpha_total_mw : float
        Total alpha power produced (MW).
    p_non_alpha_charged_mw : float
        Power from non-alpha charged particles (MW).
    p_hcd_injected_total_mw : float
        Total power injected by heating and current drive (MW).
    p_plasma_ohmic_mw : float
        Ohmic heating power (MW).
    p_plasma_rad_mw : float
        Radiated power from plasma (MW).

    Returns
    -------
    float
        Power crossing the separatrix (MW).
    """

    return (
        f_p_alpha_plasma_deposited * p_alpha_total_mw
        + p_non_alpha_charged_mw
        + p_hcd_injected_total_mw
        + p_plasma_ohmic_mw
        - p_plasma_rad_mw
    )

calculate_psep_over_r_metric(p_plasma_separatrix_mw, rmajor) staticmethod

Calculate the power crossing the separatrix per unit major radius (P_sep/R).

Parameters:

Name Type Description Default
p_plasma_separatrix_mw float

Power crossing the separatrix (MW).

 required
rmajor float

Plasma major radius (m).

 required

Returns:

Type Description
float

Power crossing the separatrix per unit major radius (MW/m).
Source code in process/models/physics/exhaust.py 
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
@staticmethod
def calculate_psep_over_r_metric(
    p_plasma_separatrix_mw: float, rmajor: float
) -> float:
    """
    Calculate the power crossing the separatrix per unit major radius (P_sep/R).

    Parameters
    ----------
    p_plasma_separatrix_mw : float
        Power crossing the separatrix (MW).
    rmajor : float
        Plasma major radius (m).

    Returns
    -------
    float
        Power crossing the separatrix per unit major radius (MW/m).
    """
    return p_plasma_separatrix_mw / rmajor

calculate_eu_demo_re_attachment_metric(p_plasma_separatrix_mw, b_plasma_toroidal_on_axis, q95, aspect, rmajor) staticmethod

Calculate the EU DEMO divertor protection re-attachment metric for plasma exhaust.

Parameters:

Name Type Description Default
p_plasma_separatrix_mw float

Power crossing the separatrix (MW).

 required
b_plasma_toroidal_on_axis float

Toroidal magnetic field on axis (T).

 required
q95 float

Safety factor at 95% flux surface.

 required
aspect float

Aspect ratio of the plasma.

 required
rmajor float

Plasma major radius (m).

 required

Returns:

Type Description
float

EU DEMO re-attachment metric (MW T /m).
References

  • M. Siccinio, G. Federici, R. Kembleton, H. Lux, F. Maviglia, and J. Morris, "Figure of merit for divertor protection in the preliminary design of the EU-DEMO reactor," Nuclear Fusion, vol. 59, no. 10, pp. 106026-106026, Jul. 2019, doi: https://doi.org/10.1088/1741-4326/ab3153.

  • H. Zohm et al., "A stepladder approach to a tokamak fusion power plant," Nuclear Fusion, vol. 57, no. 8, pp. 086002-086002, May 2017, doi: https://doi.org/10.1088/1741-4326/aa739e.

Source code in process/models/physics/exhaust.py 
 77
 78
 79
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
@staticmethod
def calculate_eu_demo_re_attachment_metric(
    p_plasma_separatrix_mw: float,
    b_plasma_toroidal_on_axis: float,
    q95: float,
    aspect: float,
    rmajor: float,
) -> float:
    """Calculate the EU DEMO divertor protection re-attachment metric for plasma exhaust.

    Parameters
    ----------
    p_plasma_separatrix_mw : float
        Power crossing the separatrix (MW).
    b_plasma_toroidal_on_axis : float
        Toroidal magnetic field on axis (T).
    q95 : float
        Safety factor at 95% flux surface.
    aspect : float
        Aspect ratio of the plasma.
    rmajor : float
        Plasma major radius (m).

    Returns
    -------
    float
        EU DEMO re-attachment metric (MW T /m).

    References
    ----------
    - M. Siccinio, G. Federici, R. Kembleton, H. Lux, F. Maviglia, and J. Morris,
      "Figure of merit for divertor protection in the preliminary design of the EU-DEMO reactor,"
      Nuclear Fusion, vol. 59, no. 10, pp. 106026-106026, Jul. 2019,
      doi: https://doi.org/10.1088/1741-4326/ab3153.

    - H. Zohm et al.,
      "A stepladder approach to a tokamak fusion power plant,"
      Nuclear Fusion, vol. 57, no. 8, pp. 086002-086002, May 2017,
      doi: https://doi.org/10.1088/1741-4326/aa739e.
    """

    return (p_plasma_separatrix_mw * b_plasma_toroidal_on_axis) / (
        q95 * aspect * rmajor
    )