Skip to content

structure

Structure

Class containing support structure calculations

This class contains routines for calculating the parameters of the support structure for a fusion power plant.

Source code in process/models/structure.py 
 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
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
class Structure:
    """Class containing support structure calculations

    This class contains routines for calculating the
    parameters of the support structure for a
    fusion power plant.
    """

    def __init__(self):
        self.outfile = constants.NOUT  # output file unit

    def run(self, output: bool = False):
        """Structure calculation caller

        This subroutine calls the support structure mass calculations.

        Parameters
        ----------
        output :
            indicate whether output should be written to the output file, or not
        """

        # Total weight of the PF coil conductor and its structure
        total_weight_pf = pfv.m_pf_coil_conductor_total + pfv.m_pf_coil_structure_total

        (
            stv.fncmass,
            stv.aintmass,
            stv.clgsmass,
            stv.coldmass,
            stv.gsmass,
        ) = self.structure(
            pv.plasma_current,
            pv.rmajor,
            pv.rminor,
            pv.kappa,
            pv.b_plasma_toroidal_on_axis,
            tfv.i_tf_sup,
            pfv.i_pf_conductor,
            bv.dr_tf_inner_bore + bv.dr_tf_outboard + bv.dr_tf_inboard,
            bv.z_tf_inside_half,
            fwbsv.whtshld,
            divv.m_div_plate,
            total_weight_pf,
            tfv.m_tf_coils_total,
            fwbsv.m_fw_total,
            fwbsv.m_blkt_total,
            fwbsv.m_fw_blkt_div_coolant_total,
            fwbsv.dewmkg,
            output=output,
        )

    def structure(
        self,
        ai,
        r0,
        a,
        akappa,
        b0,
        i_tf_sup,
        i_pf_conductor,
        tf_h_width,
        tfhmax,
        shldmass,
        dvrtmass,
        pfmass,
        tfmass,
        m_fw_total,
        blmass,
        m_fw_blkt_div_coolant_total,
        dewmass,
        output,
    ):
        """Method to calculate mass of support structure

        Programmed to match the ITER (i.e. B. Spears) rules.

        Parameters
        ----------
        ai : float
            plasma current (max design value) (A)
        r0 : float
            plasma major radius (m)
        a : float
            plasma minor radius (m)
        akappa : float
            plasma elongation
        b0 : float
            axial B-field (T)
        i_tf_sup :
            switch denoting whether TF coils are superconducting
        i_pf_conductor : integer
            switch denoting whether PF & CS coils are resistive
        tf_h_width : float
            TF coil horizontal dr_bore (m)
        tfhmax : float
            TF coil max height (m)
        shldmass : float
            total mass of shield (kg)
        dvrtmass : float
            total mass of divertor and assoc. structure (kg)
        pfmass : float
            total mass of PF coils plus cases (kg)
        tfmass : float
            total mass of TF coils plus cases (kg)
        blmass : float
            blanket mass (kg)
        m_fw_total : float
            first wall mass (kg)
        m_fw_blkt_div_coolant_total : float
            total water coolant mass (kg)
        dewmass : float
            vacuum vessel + cryostat mass (kg)
        output :
            indicate whether output should be written to the output file, or not


        Returns
        -------
        type
            fncmass (`float`) mass of outer pf coil support fence (kg)
            - aintmass (`float`) mass of intercoil support (kg)
            - clgsmass (`float`) coil gravity support mass (kg)
            - coldmass (`float`) total mass of cryogenic temp. stuff (kg)
            - gsm (`float`) gravity support for magnets, and shield/blanket (kg)
        """

        #  Outer PF coil fence (1990 ITER fit)
        fncmass = 2.1e-11 * ai * ai * r0 * akappa * a

        #  Intercoil support between TF coils to react overturning moment
        #  (scaled to 1990 ITER fit)
        aintmass = 1.4e6 * (ai / 2.2e7) * b0 / 4.85e0 * tf_h_width**2 / 50.0e0
        try:
            assert aintmass < np.inf
        except AssertionError:
            logger.exception("aintmass is inf. Kludging to 1e10.")
            aintmass = 1e10

        #  Total mass of coils plus support plus vacuum vessel + cryostat
        coilmass = tfmass + pfmass + aintmass + dewmass

        #  Total mass of cooled components
        coldmass = 0.0e0
        if i_tf_sup == 1:
            coldmass = coldmass + tfmass + aintmass + dewmass
        if i_pf_conductor != 1:
            coldmass = coldmass + pfmass

        #  Coil gravity support mass
        #  Set density (kg/m3) and allowable stress (Pa)

        dens = 7.8e3
        sigal = 2.5e7
        clgsmass = coilmass * (r0 / 6.0e0) * 9.1e0 * 9.807e0 * dens / sigal

        #  Gravity support masses scaled from Spears algorithms (9/90) :

        #  Torus leg support

        ws1 = m_fw_blkt_div_coolant_total + m_fw_total + blmass + shldmass + dvrtmass
        gsm1 = 5.0e0 * 9.807e0 * ws1 * dens / sigal

        #  Ring beam

        ws2 = ws1 + tfmass + pfmass + aintmass + clgsmass
        gsm2 = 1.0e-3 * 34.77e0 * (r0 + 1.0e0) * math.sqrt(0.001e0 * ws2) * dens

        #  Ring legs (JG: this term may be too big, need to check)

        gsm3 = 1.0e-6 * 0.3e0 * (tfhmax + 2.0e0) * ws2 * dens

        gsm = gsm1 + gsm2 + gsm3

        #  Output section

        if output:
            po.oheadr(self.outfile, "Support Structure")
            po.ovarre(
                self.outfile,
                "Outer PF coil fence mass (kg)",
                "(fncmass)",
                fncmass,
                "OP ",
            )
            po.ovarre(
                self.outfile,
                "Intercoil support structure mass (kg)",
                "(aintmass)",
                aintmass,
                "OP ",
            )
            po.ovarre(
                self.outfile,
                "Mass of cooled components (kg)",
                "(coldmass)",
                coldmass,
                "OP ",
            )
            po.ovarre(
                self.outfile,
                "Gravity support structure mass (kg)",
                "(clgsmass)",
                clgsmass,
                "OP ",
            )
            po.ovarre(self.outfile, "Torus leg support mass (kg)", "(gsm1)", gsm1, "OP ")
            po.ovarre(self.outfile, "Ring beam mass (kg)", "(gsm2)", gsm2, "OP ")
            po.ovarre(self.outfile, "Ring legs mass (kg)", "(gsm3)", gsm3, "OP ")

        return fncmass, aintmass, clgsmass, coldmass, gsm

outfile = constants.NOUT instance-attribute

run(output=False)

Structure calculation caller

This subroutine calls the support structure mass calculations.

Parameters:

Name Type Description Default
output bool

indicate whether output should be written to the output file, or not

 False
Source code in process/models/structure.py 
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
def run(self, output: bool = False):
    """Structure calculation caller

    This subroutine calls the support structure mass calculations.

    Parameters
    ----------
    output :
        indicate whether output should be written to the output file, or not
    """

    # Total weight of the PF coil conductor and its structure
    total_weight_pf = pfv.m_pf_coil_conductor_total + pfv.m_pf_coil_structure_total

    (
        stv.fncmass,
        stv.aintmass,
        stv.clgsmass,
        stv.coldmass,
        stv.gsmass,
    ) = self.structure(
        pv.plasma_current,
        pv.rmajor,
        pv.rminor,
        pv.kappa,
        pv.b_plasma_toroidal_on_axis,
        tfv.i_tf_sup,
        pfv.i_pf_conductor,
        bv.dr_tf_inner_bore + bv.dr_tf_outboard + bv.dr_tf_inboard,
        bv.z_tf_inside_half,
        fwbsv.whtshld,
        divv.m_div_plate,
        total_weight_pf,
        tfv.m_tf_coils_total,
        fwbsv.m_fw_total,
        fwbsv.m_blkt_total,
        fwbsv.m_fw_blkt_div_coolant_total,
        fwbsv.dewmkg,
        output=output,
    )

structure(ai, r0, a, akappa, b0, i_tf_sup, i_pf_conductor, tf_h_width, tfhmax, shldmass, dvrtmass, pfmass, tfmass, m_fw_total, blmass, m_fw_blkt_div_coolant_total, dewmass, output)

Method to calculate mass of support structure

Programmed to match the ITER (i.e. B. Spears) rules.

Parameters:

Name Type Description Default
ai float

plasma current (max design value) (A)

 required
r0 float

plasma major radius (m)

 required
a float

plasma minor radius (m)

 required
akappa float

plasma elongation

 required
b0 float

axial B-field (T)

 required
i_tf_sup

switch denoting whether TF coils are superconducting

 required
i_pf_conductor integer

switch denoting whether PF & CS coils are resistive

 required
tf_h_width float

TF coil horizontal dr_bore (m)

 required
tfhmax float

TF coil max height (m)

 required
shldmass float

total mass of shield (kg)

 required
dvrtmass float

total mass of divertor and assoc. structure (kg)

 required
pfmass float

total mass of PF coils plus cases (kg)

 required
tfmass float

total mass of TF coils plus cases (kg)

 required
blmass float

blanket mass (kg)

 required
m_fw_total float

first wall mass (kg)

 required
m_fw_blkt_div_coolant_total float

total water coolant mass (kg)

 required
dewmass float

vacuum vessel + cryostat mass (kg)

 required
output

indicate whether output should be written to the output file, or not

 required

Returns:

Type Description
type

fncmass (float) mass of outer pf coil support fence (kg) - aintmass (float) mass of intercoil support (kg) - clgsmass (float) coil gravity support mass (kg) - coldmass (float) total mass of cryogenic temp. stuff (kg) - gsm (float) gravity support for magnets, and shield/blanket (kg)
Source code in process/models/structure.py 
 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
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
def structure(
    self,
    ai,
    r0,
    a,
    akappa,
    b0,
    i_tf_sup,
    i_pf_conductor,
    tf_h_width,
    tfhmax,
    shldmass,
    dvrtmass,
    pfmass,
    tfmass,
    m_fw_total,
    blmass,
    m_fw_blkt_div_coolant_total,
    dewmass,
    output,
):
    """Method to calculate mass of support structure

    Programmed to match the ITER (i.e. B. Spears) rules.

    Parameters
    ----------
    ai : float
        plasma current (max design value) (A)
    r0 : float
        plasma major radius (m)
    a : float
        plasma minor radius (m)
    akappa : float
        plasma elongation
    b0 : float
        axial B-field (T)
    i_tf_sup :
        switch denoting whether TF coils are superconducting
    i_pf_conductor : integer
        switch denoting whether PF & CS coils are resistive
    tf_h_width : float
        TF coil horizontal dr_bore (m)
    tfhmax : float
        TF coil max height (m)
    shldmass : float
        total mass of shield (kg)
    dvrtmass : float
        total mass of divertor and assoc. structure (kg)
    pfmass : float
        total mass of PF coils plus cases (kg)
    tfmass : float
        total mass of TF coils plus cases (kg)
    blmass : float
        blanket mass (kg)
    m_fw_total : float
        first wall mass (kg)
    m_fw_blkt_div_coolant_total : float
        total water coolant mass (kg)
    dewmass : float
        vacuum vessel + cryostat mass (kg)
    output :
        indicate whether output should be written to the output file, or not


    Returns
    -------
    type
        fncmass (`float`) mass of outer pf coil support fence (kg)
        - aintmass (`float`) mass of intercoil support (kg)
        - clgsmass (`float`) coil gravity support mass (kg)
        - coldmass (`float`) total mass of cryogenic temp. stuff (kg)
        - gsm (`float`) gravity support for magnets, and shield/blanket (kg)
    """

    #  Outer PF coil fence (1990 ITER fit)
    fncmass = 2.1e-11 * ai * ai * r0 * akappa * a

    #  Intercoil support between TF coils to react overturning moment
    #  (scaled to 1990 ITER fit)
    aintmass = 1.4e6 * (ai / 2.2e7) * b0 / 4.85e0 * tf_h_width**2 / 50.0e0
    try:
        assert aintmass < np.inf
    except AssertionError:
        logger.exception("aintmass is inf. Kludging to 1e10.")
        aintmass = 1e10

    #  Total mass of coils plus support plus vacuum vessel + cryostat
    coilmass = tfmass + pfmass + aintmass + dewmass

    #  Total mass of cooled components
    coldmass = 0.0e0
    if i_tf_sup == 1:
        coldmass = coldmass + tfmass + aintmass + dewmass
    if i_pf_conductor != 1:
        coldmass = coldmass + pfmass

    #  Coil gravity support mass
    #  Set density (kg/m3) and allowable stress (Pa)

    dens = 7.8e3
    sigal = 2.5e7
    clgsmass = coilmass * (r0 / 6.0e0) * 9.1e0 * 9.807e0 * dens / sigal

    #  Gravity support masses scaled from Spears algorithms (9/90) :

    #  Torus leg support

    ws1 = m_fw_blkt_div_coolant_total + m_fw_total + blmass + shldmass + dvrtmass
    gsm1 = 5.0e0 * 9.807e0 * ws1 * dens / sigal

    #  Ring beam

    ws2 = ws1 + tfmass + pfmass + aintmass + clgsmass
    gsm2 = 1.0e-3 * 34.77e0 * (r0 + 1.0e0) * math.sqrt(0.001e0 * ws2) * dens

    #  Ring legs (JG: this term may be too big, need to check)

    gsm3 = 1.0e-6 * 0.3e0 * (tfhmax + 2.0e0) * ws2 * dens

    gsm = gsm1 + gsm2 + gsm3

    #  Output section

    if output:
        po.oheadr(self.outfile, "Support Structure")
        po.ovarre(
            self.outfile,
            "Outer PF coil fence mass (kg)",
            "(fncmass)",
            fncmass,
            "OP ",
        )
        po.ovarre(
            self.outfile,
            "Intercoil support structure mass (kg)",
            "(aintmass)",
            aintmass,
            "OP ",
        )
        po.ovarre(
            self.outfile,
            "Mass of cooled components (kg)",
            "(coldmass)",
            coldmass,
            "OP ",
        )
        po.ovarre(
            self.outfile,
            "Gravity support structure mass (kg)",
            "(clgsmass)",
            clgsmass,
            "OP ",
        )
        po.ovarre(self.outfile, "Torus leg support mass (kg)", "(gsm1)", gsm1, "OP ")
        po.ovarre(self.outfile, "Ring beam mass (kg)", "(gsm2)", gsm2, "OP ")
        po.ovarre(self.outfile, "Ring legs mass (kg)", "(gsm3)", gsm3, "OP ")

    return fncmass, aintmass, clgsmass, coldmass, gsm