Skip to content

forces

Module for coil force calculations in stellarators.

calculate_max_force_density(a_tf_wp_no_insulation, data)

Calculate the maximum force density in the TF coil winding pack from scaling. [MN/m3]

Parameters:

Name Type Description Default
a_tf_wp_no_insulation
required
data DataStructure

data structure object

 required
Source code in process/models/stellarator/coils/forces.py 
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
def calculate_max_force_density(a_tf_wp_no_insulation, data: DataStructure):
    """Calculate the maximum force density in the TF coil winding pack from scaling.
     [MN/m3]

    Parameters
    ----------
    a_tf_wp_no_insulation :

    data: DataStructure
        data structure object
    """
    data.tfcoil.max_force_density = (
        data.stellarator_config.stella_config_max_force_density
        * data.stellarator.f_st_i_total
        / data.stellarator.f_st_n_coils
        * data.tfcoil.b_tf_inboard_peak_symmetric
        / data.stellarator_config.stella_config_wp_bmax
        * data.stellarator_config.stella_config_wp_area
        / a_tf_wp_no_insulation
    )

calculate_max_force_density_mnm(data)

Calculate the maximum force per meter in the TF coil winding pack from scaling. [MN/m]

Source code in process/models/stellarator/coils/forces.py 
28
29
30
31
32
33
34
35
36
37
38
def calculate_max_force_density_mnm(data: DataStructure):
    """Calculate the maximum force per meter in the TF coil winding pack from scaling.
    [MN/m]
    """
    return (
        data.stellarator_config.stella_config_max_force_density_mnm
        * data.stellarator.f_st_i_total
        / data.stellarator.f_st_n_coils
        * data.tfcoil.b_tf_inboard_peak_symmetric
        / data.stellarator_config.stella_config_wp_bmax
    )

calculate_maximum_stress(data)

Approximate, very simple maximum stress (needed for limitation of icc 32), in Pa

Source code in process/models/stellarator/coils/forces.py 
41
42
43
44
45
46
47
def calculate_maximum_stress(data):
    """Approximate, very simple maximum stress (needed for limitation of icc 32),
    in Pa
    """
    data.tfcoil.sig_tf_wp = (
        data.tfcoil.max_force_density * data.tfcoil.dr_tf_wp_with_insulation * 1.0e6
    )

calculate_max_lateral_force_density(a_tf_wp_no_insulation, data)

Calculate the maximum lateral force density in the TF coil winding pack from scaling. [MN/m3]

Parameters:

Name Type Description Default
a_tf_wp_no_insulation
required
data DataStructure

data structure object

 required
Source code in process/models/stellarator/coils/forces.py 
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
def calculate_max_lateral_force_density(a_tf_wp_no_insulation, data: DataStructure):
    """Calculate the maximum lateral force density in the TF coil winding pack
     from scaling. [MN/m3]

    Parameters
    ----------
    a_tf_wp_no_insulation :

    data: DataStructure
        data structure object
    """
    return (
        data.stellarator_config.stella_config_max_lateral_force_density
        * data.stellarator.f_st_i_total
        / data.stellarator.f_st_n_coils
        * data.tfcoil.b_tf_inboard_peak_symmetric
        / data.stellarator_config.stella_config_wp_bmax
        * data.stellarator_config.stella_config_wp_area
        / a_tf_wp_no_insulation
    )

calculate_max_radial_force_density(a_tf_wp_no_insulation, data)

Calculate the maximum radial force density in the TF coil winding pack from scaling. [MN/m3]

Parameters:

Name Type Description Default
a_tf_wp_no_insulation
required
data

data structure object

 required
Source code in process/models/stellarator/coils/forces.py 
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
def calculate_max_radial_force_density(a_tf_wp_no_insulation, data):
    """Calculate the maximum radial force density in the TF coil winding pack
    from scaling. [MN/m3]

    Parameters
    ----------
    a_tf_wp_no_insulation :

    data: DataStructure
        data structure object
    """
    return (
        data.stellarator_config.stella_config_max_radial_force_density
        * data.stellarator.f_st_i_total
        / data.stellarator.f_st_n_coils
        * data.tfcoil.b_tf_inboard_peak_symmetric
        / data.stellarator_config.stella_config_wp_bmax
        * data.stellarator_config.stella_config_wp_area
        / a_tf_wp_no_insulation
    )

calculate_centering_force_max_mn(data)

Calculate the maximum centering force in the TF coils from scaling. [MN]

Source code in process/models/stellarator/coils/forces.py 
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
def calculate_centering_force_max_mn(data: DataStructure):
    """Calculate the maximum centering force in the TF coils from scaling. [MN]"""
    return (
        data.stellarator_config.stella_config_centering_force_max_mn
        * data.stellarator.f_st_i_total
        / data.stellarator.f_st_n_coils
        * data.tfcoil.b_tf_inboard_peak_symmetric
        / data.stellarator_config.stella_config_wp_bmax
        * data.stellarator_config.stella_config_coillength
        / data.tfcoil.n_tf_coils
        / data.tfcoil.len_tf_coil
    )

calculate_centering_force_min_mn(data)

Calculate the minimum centering force in the TF coils from scaling. [MN]

Source code in process/models/stellarator/coils/forces.py 
108
109
110
111
112
113
114
115
116
117
118
119
def calculate_centering_force_min_mn(data: DataStructure):
    """Calculate the minimum centering force in the TF coils from scaling. [MN]"""
    return (
        data.stellarator_config.stella_config_centering_force_min_mn
        * data.stellarator.f_st_i_total
        / data.stellarator.f_st_n_coils
        * data.tfcoil.b_tf_inboard_peak_symmetric
        / data.stellarator_config.stella_config_wp_bmax
        * data.stellarator_config.stella_config_coillength
        / data.tfcoil.n_tf_coils
        / data.tfcoil.len_tf_coil
    )

calculate_centering_force_avg_mn(data)

Calculate the average centering force in the TF coils from scaling. [MN]

Source code in process/models/stellarator/coils/forces.py 
122
123
124
125
126
127
128
129
130
131
132
133
def calculate_centering_force_avg_mn(data: DataStructure):
    """Calculate the average centering force in the TF coils from scaling. [MN]"""
    return (
        data.stellarator_config.stella_config_centering_force_avg_mn
        * data.stellarator.f_st_i_total
        / data.stellarator.f_st_n_coils
        * data.tfcoil.b_tf_inboard_peak_symmetric
        / data.stellarator_config.stella_config_wp_bmax
        * data.stellarator_config.stella_config_coillength
        / data.tfcoil.n_tf_coils
        / data.tfcoil.len_tf_coil
    )