Skip to content

Updates to reflect compound material layers #133

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
wants to merge 2 commits into
base: master
Choose a base branch
from
Open

Updates to reflect compound material layers #133

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
29d85c0
Update to reflect compound material layers
adrianaghiozzi Apr 3, 2024
67af9fd
Initiate coil_tech as compound materials and add expressions for frac…
adrianaghiozzi Apr 19, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
54 changes: 54 additions & 0 deletions src/expressions/dynamic.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -463,6 +463,60 @@ dyexp["build.tf.ripple"] =
dyexp["build.tf.wedge_thickness"] =
(; build, _...) -> 2π * get_build_layer(build.layer; type=_tf_, fs=_hfs_).end_radius / build.tf.coils_n

dyexp["build.tf.technology.fraction_steel"] =
(; build, _...) -> begin
material = IMAS.get_build_layer(build.layer; type=_tf_, fs=_hfs_).material
return sum([mat.name == "steel" ? mat.composition : 0.0 for mat in material])
end

dyexp["build.tf.technology.ratio_SC_to_copper"] =
(; build, _...) -> begin
material = IMAS.get_build_layer(build.layer; type=_tf_, fs=_hfs_).material
return sum([contains(mat.name, "rebco") || contains(mat.name, "nb3sn") || contains(mat.name, "nbti") ? mat.composition : 0 for mat in material]) / sum([mat.name == "copper" ? mat.composition : 0 for mat in material])
end

dyexp["build.tf.technology.fraction_void"] =
(; build, _...) -> begin
material = IMAS.get_build_layer(build.layer; type=_tf_, fs=_hfs_).material
return sum([mat.name == "vacuum" ? mat.composition : 0.0 for mat in material])
end

dyexp["build.oh.technology.fraction_steel"] =
(; build, _...) -> begin
material = IMAS.get_build_layer(build.layer; type=_oh_).material
return sum([mat.name == "steel" ? mat.composition : 0.0 for mat in material])
end

dyexp["build.oh.technology.ratio_SC_to_copper"] =
(; build, _...) -> begin
material = IMAS.get_build_layer(build.layer; type=_oh_).material
return sum([contains(mat.name, "rebco") || contains(mat.name, "nb3sn") || contains(mat.name, "nbti") ? mat.composition : 0 for mat in material]) / sum([mat.name == "copper" ? mat.composition : 0 for mat in material])
end

dyexp["build.oh.technology.fraction_void"] =
(; build, _...) -> begin
material = IMAS.get_build_layer(build.layer; type=_oh_).material
return sum([mat.name == "vacuum" ? mat.composition : 0.0 for mat in material])
end

dyexp["build.pf_active.technology.fraction_steel"] =
(; build, _...) -> begin
material = build.pf_active.technology.material
return sum([mat.name == "steel" ? mat.composition : 0.0 for mat in material])
end

dyexp["build.pf_active.technology.ratio_SC_to_copper"] =
(; build, _...) -> begin
material = build.pf_active.technology.material
return sum([contains(mat.name, "rebco") || contains(mat.name, "nb3sn") || contains(mat.name, "nbti") ? mat.composition : 0 for mat in material]) / sum([mat.name == "copper" ? mat.composition : 0 for mat in material])
end

dyexp["build.pf_active.technology.fraction_void"] =
(; build, _...) -> begin
material = build.pf_active.technology.material
return sum([mat.name == "vacuum" ? mat.composition : 0.0 for mat in material])
end

#= ======= =#
# costing #
#= ======= =#
Expand Down
4 changes: 2 additions & 2 deletions src/physics/build.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -154,7 +154,7 @@ function structures_mask(bd::IMAS.build; ngrid::Int=257, border_fraction::Real=0
# start from the first vacuum that goes to zero outside of the TF
start_from = -1
for k in get_build_indexes(bd.layer; fs=_out_)
if bd.layer[k].material == "vacuum" && minimum(bd.layer[k].outline.r) < bd.layer[1].end_radius
if "vacuum" ∈ [bd.layer[k].material[i].name for i in length(bd.layer[k].material)] && minimum(bd.layer[k].outline.r) < bd.layer[1].end_radius
start_from = k
break
end
Expand All @@ -168,7 +168,7 @@ function structures_mask(bd::IMAS.build; ngrid::Int=257, border_fraction::Real=0
end
if valid && !ismissing(layer.outline, :r)
outline = collect(zip(layer.outline.r, layer.outline.z))
if (layer.material == "vacuum") && (layer.side != Int(_in_))
if ("vacuum" ∈ [layer.material[i].name for i in length(layer.material)]) && (layer.side != Int(_in_))
for (kr, rr) in enumerate(rmask)
for (kz, zz) in enumerate(zmask)
if PolygonOps.inpolygon((rr, zz), outline) != 0
Expand Down
125 changes: 97 additions & 28 deletions src/physics/technology.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -41,53 +41,131 @@ function coil_technology(coil_tech::Union{IMAS.build__pf_active__technology,IMAS
end

if technology == :copper
coil_tech.material = "copper"
resize!(coil_tech.material, 2)
coil_tech.temperature = 293.0
coil_tech.fraction_steel = 0.0
coil_tech.ratio_SC_to_copper = 0.0
coil_tech.fraction_void = 0.2
coil_tech.material[1].name = "copper"
coil_tech.material[1].composition = 0.8

coil_tech.material[2].name = "vacuum"
coil_tech.material[2].composition = 0.2

elseif technology ∈ (:nb3sn, :nbti, :nb3sn_iter, :nb3sn_kdemo, :rebco)
resize!(coil_tech.material, 4)
if technology == :nb3sn
coil_tech.temperature = 4.2
coil_tech.material = "nb3sn"
coil_tech.fraction_void = 0.1
coil_tech.material[1].name = "nb3sn"
coil_tech.material[1].composition = 0.2

coil_tech.material[2].name = "vacuum"
coil_tech.material[2].composition = 0.1

coil_tech.material[3].name = "steel"
coil_tech.material[3].composition = 0.5

coil_tech.material[4].name = "copper"
coil_tech.material[4].composition = 0.2
elseif technology == :nbti
coil_tech.temperature = 4.2
coil_tech.material = "nbti"
coil_tech.fraction_void = 0.2 # from Supercond. Sci. Technol. 36 (2023) 075009
coil_tech.material[1].name = "nbti"
coil_tech.material[1].composition = 0.15

coil_tech.material[2].name = "vacuum"
coil_tech.material[2].composition = 0.2 # from Supercond. Sci. Technol. 36 (2023) 075009

coil_tech.material[3].name = "steel"
coil_tech.material[3].composition = 0.5

coil_tech.material[4].name = "copper"
coil_tech.material[4].composition = 0.15
elseif technology == :nb3sn_iter
coil_tech.temperature = 4.2
coil_tech.material = "nb3sn_iter"
coil_tech.fraction_void = 0.1
coil_tech.material[1].name = "nb3sn"
coil_tech.material[1].composition = 0.2

coil_tech.material[2].name = "vacuum"
coil_tech.material[2].composition = 0.1

coil_tech.material[3].name = "steel"
coil_tech.material[3].composition = 0.5

coil_tech.material[4].name = "copper"
coil_tech.material[4].composition = 0.2
elseif technology == :nb3sn_kdemo
coil_tech.temperature = 4.2
coil_tech.material = "nb3sn_kdemo"
coil_tech.material[1].name = "nb3sn_kdemo"
coil_tech.material[1].composition = 0.2

coil_tech.material[2].name = "vacuum"
coil_tech.material[2].composition = 0.1

coil_tech.material[3].name = "steel"
coil_tech.material[3].composition = 0.5

coil_tech.material[4].name = "copper"
coil_tech.material[4].composition = 0.2

if coil_type == :tf
coil_tech.fraction_void = 0.26 # from NF 55 (2015) 053027, Table 2
coil_tech.material[1].name = "nb3sn_kdemo"
coil_tech.material[1].composition = 0.12

coil_tech.material[2].name = "vacuum"
coil_tech.material[2].composition = 0.26 # from NF 55 (2015) 053027, Table 2

coil_tech.material[3].name = "steel"
coil_tech.material[3].composition = 0.5

coil_tech.material[4].name = "copper"
coil_tech.material[4].composition = 0.12
end
else
coil_tech.temperature = 4.2
coil_tech.material = "rebco"
coil_tech.material[1].name = "rebco"
coil_tech.material[1].composition = 0.2

coil_tech.material[2].name = "vacuum"
coil_tech.material[2].composition = 0.1

coil_tech.material[3].name = "steel"
coil_tech.material[3].composition = 0.5

coil_tech.material[4].name = "copper"
coil_tech.material[4].composition = 0.2

end
coil_tech.fraction_steel = 0.5
coil_tech.ratio_SC_to_copper = 1.0
coil_tech.fraction_void = 0.1
end

if technology == :nb3sn_iter
resize!(coil_tech.material, 4)
coil_tech.material[1].name = "nb3sn_iter"

coil_tech.material[2].name = "vacuum"
coil_tech.material[2].composition = 0.1

coil_tech.material[3].name = "steel"

coil_tech.material[4].name = "copper"

if coil_type == :oh
coil_tech.thermal_strain = -0.64
coil_tech.JxB_strain = -0.05
coil_tech.fraction_steel = 0.46

coil_tech.material[1].composition = 0.22
coil_tech.material[3].composition = 0.46
coil_tech.material[4].composition = 0.22
elseif coil_type == :tf
coil_tech.thermal_strain = -0.69
coil_tech.JxB_strain = -0.13
coil_tech.fraction_steel = 0.55

coil_tech.material[1].composition = 0.175
coil_tech.material[3].composition = 0.55
coil_tech.material[4].composition = 0.175
elseif coil_type == :pf_active
coil_tech.thermal_strain = -0.64
coil_tech.JxB_strain = -0.05
coil_tech.fraction_steel = 0.46

coil_tech.material[1].composition = 0.22
coil_tech.material[3].composition = 0.46
coil_tech.material[4].composition = 0.22
end
end

Expand All @@ -97,15 +175,6 @@ function coil_technology(coil_tech::Union{IMAS.build__pf_active__technology,IMAS
return coil_tech
end

function fraction_conductor(coil_tech::Union{IMAS.build__pf_active__technology,IMAS.build__oh__technology,IMAS.build__tf__technology})
frac = 1.0 - coil_tech.fraction_steel - coil_tech.fraction_void # fraction of coil that is a conductor
@assert frac > 0.0 "coil technology has no room for conductor"
if coil_tech.material == "copper"
return frac
else
return frac * coil_tech.ratio_SC_to_copper / (1.0 + coil_tech.ratio_SC_to_copper) # fraction of coil that is Nb3Sn superconductor
end
end

function GAMBL_blanket(bm::IMAS.blanket__module)
layers = resize!(bm.layer, 3)
Expand Down