def _make_pf_coils(self):
list_of_components = []
if (self.pf_coil_vertical_thicknesses,
self.pf_coil_radial_thicknesses,
self.pf_coil_to_rear_blanket_radial_gap) != (None, None, None):
self._pf_coil = paramak.PoloidalFieldCoilSet(
heights=self.pf_coil_vertical_thicknesses,
widths=self.pf_coil_radial_thicknesses,
center_points=self._pf_coils_xy_values,
rotation_angle=self.rotation_angle,
stp_filename='pf_coils.stp',
stl_filename='pf_coils.stl',
name="pf_coil",
material_tag="pf_coil_mat",
)
self._pf_coils_casing = paramak.PoloidalFieldCoilCaseSetFC(
pf_coils=self._pf_coil,
casing_thicknesses=self.pf_coil_case_thickness,
rotation_angle=self.rotation_angle,
stp_filename='pf_coil_cases.stp',
stl_filename='pf_coil_cases.stl',
name="pf_coil_case",
material_tag="pf_coil_case_mat",
)
list_of_components = [self._pf_coils_casing, self._pf_coil]
return list_of_components
def test_PoloidalFieldCoilCaseSetFC_with_number_thickness(self):
"""Creates a set of PF coil cases from a list of PF coils with a
single numerical thicknesses."""
pf_coils_1 = paramak.PoloidalFieldCoil(height=10,
width=10,
center_point=(100, 100))
pf_coils_2 = paramak.PoloidalFieldCoil(height=10,
width=10,
center_point=(100, 150))
pf_coils_3 = paramak.PoloidalFieldCoil(height=20,
width=20,
center_point=(50, 200))
pf_coils_4 = paramak.PoloidalFieldCoil(height=20,
width=40,
center_point=(50, 50))
test_shape = paramak.PoloidalFieldCoilCaseSetFC(
pf_coils=[pf_coils_1, pf_coils_2, pf_coils_3, pf_coils_4],
casing_thicknesses=10,
)
assert test_shape.casing_thicknesses == 10
assert test_shape.solid is not None
assert len(test_shape.solid.Solids()) == 4
def _make_coils(self):
list_of_components = []
self._inboard_tf_coils = paramak.CenterColumnShieldCylinder(
height=self._blanket_rear_wall_end_height * 2,
inner_radius=self._inboard_tf_coils_start_radius,
outer_radius=self._inboard_tf_coils_end_radius,
rotation_angle=self.rotation_angle,
stp_filename="inboard_tf_coils.stp",
stl_filename="inboard_tf_coils.stl",
name="inboard_tf_coils",
material_tag="inboard_tf_coils_mat",
)
list_of_components.append(self._inboard_tf_coils)
if self._tf_info_provided:
self._tf_coil = paramak.ToroidalFieldCoilCoatHanger(
with_inner_leg=False,
horizontal_start_point=(
self._inboard_tf_coils_start_radius,
self._blanket_rear_wall_end_height,
),
vertical_mid_point=(self._outboard_tf_coil_start_radius, 0),
thickness=self.outboard_tf_coil_radial_thickness,
number_of_coils=self.number_of_tf_coils,
distance=self.outboard_tf_coil_poloidal_thickness,
stp_filename="outboard_tf_coil.stp",
stl_filename="outboard_tf_coil.stl",
rotation_angle=self.rotation_angle,
horizontal_length=self._outboard_tf_coils_horizontal_length,
vertical_length=self._outboard_tf_coils_vertical_height)
list_of_components.append(self._tf_coil)
if self._pf_info_provided:
self._pf_coil = paramak.PoloidalFieldCoilSet(
heights=self.pf_coil_vertical_thicknesses,
widths=self.pf_coil_radial_thicknesses,
center_points=self._pf_coils_xy_values,
rotation_angle=self.rotation_angle,
stp_filename='pf_coils.stp',
stl_filename='pf_coils.stl',
name="pf_coil",
material_tag="pf_coil_mat",
)
list_of_components.append(self._pf_coil)
self._pf_coils_casing = paramak.PoloidalFieldCoilCaseSetFC(
pf_coils=self._pf_coil,
casing_thicknesses=self.pf_coil_case_thickness,
rotation_angle=self.rotation_angle,
stp_filename='pf_coil_cases.stp',
stl_filename='pf_coil_cases.stl',
name="pf_coil_case",
material_tag="pf_coil_case_mat",
)
list_of_components.append(self._pf_coils_casing)
return list_of_components
def test_PoloidalFieldCoilCaseSetFC_from_list(self):
"""Creates a set of PF coil cases from a list of PF coils"""
pf_coils_1 = paramak.PoloidalFieldCoil(height=10,
width=10,
center_point=(100, 100),
rotation_angle=180)
pf_coils_2 = paramak.PoloidalFieldCoil(height=10,
width=10,
center_point=(100, 150),
rotation_angle=180)
pf_coils_3 = paramak.PoloidalFieldCoil(height=20,
width=20,
center_point=(50, 200),
rotation_angle=180)
pf_coils_4 = paramak.PoloidalFieldCoil(height=20,
width=40,
center_point=(50, 50),
rotation_angle=180)
test_shape = paramak.PoloidalFieldCoilCaseSetFC(
pf_coils=[pf_coils_1, pf_coils_2, pf_coils_3, pf_coils_4],
casing_thicknesses=[5, 5, 10, 10],
rotation_angle=180)
assert test_shape.solid is not None
assert len(test_shape.solid.Solids()) == 4
def test_PoloidalFieldCoilSet_incorrect_lengths():
"""Checks PoloidalFieldCoilSet with the wrong number of casing thicknesses using a list"""
pf_coils_1 = paramak.PoloidalFieldCoil(height=10,
width=10,
center_point=(100, 100),
rotation_angle=180)
paramak.PoloidalFieldCoilCaseSetFC(
pf_coils=[pf_coils_1],
casing_thicknesses=[5, 5, 10, 10],
rotation_angle=180).solid
def test_PoloidalFieldCoilSet_incorrect_lengths_FC():
"""Checks PoloidalFieldCoilSet with the wrong number of casing thicknesses using a coil set object"""
pf_coils_set = paramak.PoloidalFieldCoilSet(
heights=[10, 10, 20, 20],
widths=[10, 10, 20, 40],
center_points=[(100, 100), (100, 150), (50, 200), (50, 50)],
rotation_angle=180)
paramak.PoloidalFieldCoilCaseSetFC(pf_coils=pf_coils_set,
casing_thicknesses=[5, 5, 10],
rotation_angle=180).solid
def setUp(self):
self.pf_coils_set = paramak.PoloidalFieldCoilSet(
heights=[10, 10, 20, 20],
widths=[10, 10, 20, 40],
center_points=[(100, 100), (100, 150), (50, 200), (50, 50)],
)
self.test_shape = paramak.PoloidalFieldCoilCaseSetFC(
pf_coils=self.pf_coils_set,
casing_thicknesses=[5, 10, 5, 10],
)
def main():
number_of_toroidal_field_coils = 8
angle_offset = (360 / number_of_toroidal_field_coils) / 2.
tf_coil_thickness = 50
tf_coil_distance = 130
tf_coil_casing_distance = tf_coil_distance + 40
tf_coil_casing_thickness = 20
inner_tf_case = paramak.InnerTfCoilsFlat(
height=1800,
inner_radius=330,
outer_radius=430,
number_of_coils=number_of_toroidal_field_coils,
gap_size=5,
rotation_angle=180,
azimuth_start_angle=angle_offset)
tf_coils = paramak.ToroidalFieldCoilPrincetonD(
R1=400,
R2=1500, # height
thickness=tf_coil_thickness,
distance=tf_coil_distance,
number_of_coils=number_of_toroidal_field_coils,
rotation_angle=180,
)
tf_coil_casing = paramak.TFCoilCasing(
magnet=tf_coils,
distance=tf_coil_casing_distance,
inner_offset=tf_coil_casing_thickness,
outer_offset=tf_coil_casing_thickness,
vertical_section_offset=tf_coil_casing_thickness,
# rotation_angle=180, # producing occational errors with this arg
)
pf_coils = paramak.PoloidalFieldCoilSet(
heights=[100, 120, 80, 80, 120, 180],
widths=[100, 120, 80, 80, 120, 180],
center_points=[(530, 1030), (1370, 790), (1740, 250), (1750, -250),
(1360, -780), (680, -1050)],
rotation_angle=180)
pf_coils_casing = paramak.PoloidalFieldCoilCaseSetFC(
pf_coils=pf_coils, casing_thicknesses=[10] * 6, rotation_angle=180)
pf_coils.export_stp('pf_coils.stp')
pf_coils_casing.export_stp('pf_coils_case.stp')
tf_coils.export_stp('tf_coil.stp')
tf_coil_casing.export_stp('tf_coil_casing.stp')
inner_tf_case.export_stp('inner_tf_case.stp')
def test_PoloidalFieldCoilCaseSetFC_from_pf_coil_set(self):
"""Creates a set of PF coil cases from a PF coils object"""
pf_coils_set = paramak.PoloidalFieldCoilSet(heights=[10, 10, 20, 20],
widths=[10, 10, 20, 40],
center_points=[(100, 100),
(100, 150),
(50, 200),
(50, 50)],
rotation_angle=180)
test_shape = paramak.PoloidalFieldCoilCaseSetFC(
pf_coils=pf_coils_set,
casing_thicknesses=[5, 5, 10, 10],
rotation_angle=180)
assert test_shape.solid is not None
assert len(test_shape.solid.Solids()) == 4
assert len(pf_coils_set.solid.Solids()) == 4
def main():
rot_angle = 180
all_components = []
plasma = paramak.Plasma(
# default parameters
rotation_angle=rot_angle,
stp_filename="plasma_shape.stp",
)
all_components.append(plasma)
component = paramak.BlanketFP(
plasma=plasma,
thickness=100,
stop_angle=90,
start_angle=-90,
offset_from_plasma=30,
rotation_angle=rot_angle,
stp_filename="blanket_constant_thickness_outboard_plasma.stp",
)
all_components.append(component)
component = paramak.BlanketCutterStar(height=2000,
width=2000,
distance=100)
all_components.append(component)
component = paramak.BlanketFP(
plasma=plasma,
thickness=100,
stop_angle=90,
start_angle=250,
offset_from_plasma=30,
rotation_angle=rot_angle,
stp_filename="blanket_constant_thickness_inboard_plasma.stp",
)
all_components.append(component)
component = paramak.BlanketFP(
plasma=plasma,
thickness=100,
stop_angle=250,
start_angle=-90,
offset_from_plasma=30,
rotation_angle=rot_angle,
stp_filename="blanket_constant_thickness_plasma.stp",
)
all_components.append(component)
CenterColumnShieldCylinder = paramak.CenterColumnShieldCylinder(
inner_radius=80,
outer_radius=100,
height=300,
rotation_angle=rot_angle,
stp_filename="center_column_shield_cylinder.stp",
)
all_components.append(CenterColumnShieldCylinder)
component = paramak.InboardFirstwallFCCS(
central_column_shield=CenterColumnShieldCylinder,
thickness=50,
rotation_angle=rot_angle,
stp_filename="firstwall_from_center_column_shield_cylinder.stp",
)
all_components.append(component)
CenterColumnShieldHyperbola = paramak.CenterColumnShieldHyperbola(
inner_radius=50,
mid_radius=75,
outer_radius=100,
height=300,
rotation_angle=rot_angle,
stp_filename="center_column_shield_hyperbola.stp",
)
all_components.append(CenterColumnShieldHyperbola)
component = paramak.InboardFirstwallFCCS(
central_column_shield=CenterColumnShieldHyperbola,
thickness=50,
rotation_angle=rot_angle,
stp_filename="firstwall_from_center_column_shield_hyperbola.stp",
)
all_components.append(component)
CenterColumnShieldCircular = paramak.CenterColumnShieldCircular(
inner_radius=50,
mid_radius=75,
outer_radius=100,
height=300,
rotation_angle=rot_angle,
stp_filename="center_column_shield_circular.stp",
)
all_components.append(CenterColumnShieldCircular)
component = paramak.InboardFirstwallFCCS(
central_column_shield=CenterColumnShieldCircular,
thickness=50,
rotation_angle=rot_angle,
stp_filename="firstwall_from_center_column_shield_circular.stp",
)
all_components.append(component)
CenterColumnShieldFlatTopHyperbola = paramak.CenterColumnShieldFlatTopHyperbola(
inner_radius=50,
mid_radius=75,
outer_radius=100,
arc_height=220,
height=300,
rotation_angle=rot_angle,
stp_filename="center_column_shield_flat_top_hyperbola.stp",
)
all_components.append(CenterColumnShieldFlatTopHyperbola)
component = paramak.InboardFirstwallFCCS(
central_column_shield=CenterColumnShieldFlatTopHyperbola,
thickness=50,
rotation_angle=rot_angle,
stp_filename=
"firstwall_from_center_column_shield_flat_top_hyperbola.stp",
)
all_components.append(component)
CenterColumnShieldFlatTopCircular = paramak.CenterColumnShieldFlatTopCircular(
inner_radius=50,
mid_radius=75,
outer_radius=100,
arc_height=220,
height=300,
rotation_angle=rot_angle,
stp_filename="center_column_shield_flat_top_Circular.stp",
)
all_components.append(CenterColumnShieldFlatTopCircular)
component = paramak.InboardFirstwallFCCS(
central_column_shield=CenterColumnShieldFlatTopCircular,
thickness=50,
rotation_angle=rot_angle,
stp_filename=
"firstwall_from_center_column_shield_flat_top_Circular.stp",
)
all_components.append(component)
CenterColumnShieldPlasmaHyperbola = paramak.CenterColumnShieldPlasmaHyperbola(
inner_radius=150,
mid_offset=50,
edge_offset=40,
height=800,
rotation_angle=rot_angle,
stp_filename="center_column_shield_plasma_hyperbola.stp",
)
all_components.append(CenterColumnShieldPlasmaHyperbola)
component = paramak.InboardFirstwallFCCS(
central_column_shield=CenterColumnShieldPlasmaHyperbola,
thickness=50,
rotation_angle=rot_angle,
stp_filename="firstwall_from_center_column_shield_plasma_hyperbola.stp",
)
all_components.append(component)
component = paramak.InnerTfCoilsCircular(
inner_radius=25,
outer_radius=100,
number_of_coils=10,
gap_size=5,
height=300,
stp_filename="inner_tf_coils_circular.stp",
)
all_components.append(component)
component = paramak.InnerTfCoilsFlat(
inner_radius=25,
outer_radius=100,
number_of_coils=10,
gap_size=5,
height=300,
stp_filename="inner_tf_coils_flat.stp",
)
all_components.append(component)
# this makes 4 pf coil cases
pf_coil_set = paramak.PoloidalFieldCoilCaseSet(
heights=[10, 10, 20, 20],
widths=[10, 10, 20, 40],
casing_thicknesses=[5, 5, 10, 10],
center_points=[(100, 100), (100, 150), (50, 200), (50, 50)],
rotation_angle=rot_angle,
stp_filename="pf_coil_case_set.stp")
all_components.append(pf_coil_set)
# this makes 4 pf coils
pf_coil_set = paramak.PoloidalFieldCoilSet(heights=[10, 10, 20, 20],
widths=[10, 10, 20, 40],
center_points=[(100, 100),
(100, 150),
(50, 200),
(50, 50)],
rotation_angle=rot_angle,
stp_filename="pf_coil_set.stp")
all_components.append(pf_coil_set)
# this makes 4 pf coil cases for the 4 pf coils made above
component = paramak.PoloidalFieldCoilCaseSetFC(
pf_coils=pf_coil_set,
casing_thicknesses=[5, 5, 10, 10],
rotation_angle=rot_angle,
stp_filename="pf_coil_cases_set.stp")
all_components.append(component)
# this makes 1 pf coils
pf_coil = paramak.PoloidalFieldCoil(center_point=(100, 100),
height=20,
width=20,
rotation_angle=rot_angle,
stp_filename="poloidal_field_coil.stp")
all_components.append(pf_coil)
# this makes one PF coil case for the provided pf coil
component = paramak.PoloidalFieldCoilCaseSetFC(
pf_coils=[pf_coil],
casing_thicknesses=[10],
rotation_angle=rot_angle,
stp_filename="pf_coil_cases_set_fc.stp")
all_components.append(component)
component = paramak.PoloidalFieldCoilCaseFC(
pf_coil=pf_coil,
casing_thickness=10,
rotation_angle=rot_angle,
stp_filename="poloidal_field_coil_case_fc.stp",
)
all_components.append(component)
component = paramak.PoloidalFieldCoilCase(
center_point=(100, 100),
coil_height=20,
coil_width=20,
casing_thickness=10,
rotation_angle=rot_angle,
stp_filename="poloidal_field_coil_case.stp",
)
all_components.append(component)
component = paramak.BlanketConstantThicknessArcV(
inner_lower_point=(300, -200),
inner_mid_point=(500, 0),
inner_upper_point=(300, 200),
thickness=100,
rotation_angle=rot_angle,
stp_filename="blanket_arc_v.stp",
)
all_components.append(component)
component = paramak.BlanketConstantThicknessArcH(
inner_lower_point=(300, -200),
inner_mid_point=(400, 0),
inner_upper_point=(300, 200),
thickness=100,
rotation_angle=rot_angle,
stp_filename="blanket_arc_h.stp",
)
all_components.append(component)
component = paramak.ToroidalFieldCoilRectangle(
horizontal_start_point=(100, 700),
vertical_mid_point=(800, 0),
thickness=150,
distance=60,
stp_filename="tf_coil_rectangle.stp",
number_of_coils=1,
)
all_components.append(component)
component = paramak.ToroidalFieldCoilCoatHanger(
horizontal_start_point=(200, 500),
horizontal_length=400,
vertical_mid_point=(700, 50),
vertical_length=500,
thickness=50,
distance=50,
stp_filename="toroidal_field_coil_coat_hanger.stp",
number_of_coils=1,
)
all_components.append(component)
component = paramak.ToroidalFieldCoilTripleArc(
R1=80,
h=200,
radii=(70, 100),
coverages=(60, 60),
thickness=30,
distance=30,
number_of_coils=1,
stp_filename="toroidal_field_coil_triple_arc.stp")
all_components.append(component)
magnet = paramak.ToroidalFieldCoilPrincetonD(
R1=80,
R2=300,
thickness=30,
distance=30,
number_of_coils=1,
stp_filename="toroidal_field_coil_princeton_d.stp")
all_components.append(magnet)
component = paramak.ITERtypeDivertor(
# default parameters
rotation_angle=rot_angle,
stp_filename="ITER_type_divertor.stp",
)
all_components.append(component)
component = paramak.PortCutterRotated(center_point=(450, 0),
polar_coverage_angle=20,
rotation_angle=10,
polar_placement_angle=45,
azimuth_placement_angle=0)
all_components.append(component)
component = paramak.PortCutterRectangular(
distance=3,
z_pos=0,
height=0.2,
width=0.4,
fillet_radius=0.02,
azimuth_placement_angle=[0, 45, 90, 180])
all_components.append(component)
component = paramak.PortCutterCircular(
distance=3,
z_pos=0.25,
radius=0.1,
# azimuth_placement_angle=[0, 45, 90, 180], # TODO: fix issue #548
azimuth_placement_angle=[0, 45, 90],
)
all_components.append(component)
component = paramak.VacuumVessel(height=2,
inner_radius=1,
thickness=0.2,
rotation_angle=270)
all_components.append(component)
component = paramak.CoolantChannelRingStraight(
height=200,
channel_radius=10,
ring_radius=70,
number_of_coolant_channels=8,
workplane="XY",
rotation_axis="Z",
stp_filename="coolant_channel_ring_straight.stp",
)
all_components.append(component)
component = paramak.CoolantChannelRingCurved(
height=200,
channel_radius=10,
ring_radius=70,
mid_offset=-20,
number_of_coolant_channels=8,
workplane="XY",
path_workplane="XZ",
stp_filename="coolant_channel_ring_curved.stp",
force_cross_section=True)
all_components.append(component)
component = paramak.RotatedIsoscelesTriangle(
height=20,
base_length=15,
pivot_angle=0,
pivot_point=(100, 50),
rotation_angle=70,
workplane='XY',
stp_filename='rotated_isosceles_triangle.stp')
all_components.append(component)
component = paramak.RotatedTrapezoid(length_1=10,
length_2=20,
length_3=30,
pivot_angle=0,
pivot_point=(100, 50),
rotation_angle=45,
stp_filename='rotated_trapezoid.stp')
all_components.append(component)
component = paramak.PoloidalSegments(number_of_segments=5,
center_point=(400, 50))
all_components.append(component)
component = paramak.TFCoilCasing(magnet=magnet,
inner_offset=10,
outer_offset=15,
vertical_section_offset=20,
distance=40)
all_components.append(component)
component = paramak.HexagonPin(length_of_side=5,
distance=10,
center_point=(10, 20))
all_components.append(component)
return all_components
def test_PoloidalFieldCoilSet_incorrect_pf_coil():
"""Checks PoloidalFieldCoilSet with the pf_coils as an incorrect entry"""
paramak.PoloidalFieldCoilCaseSetFC(
pf_coils=20,
casing_thicknesses=[5, 5, 10, 10],
rotation_angle=180).solid