Skip to content

fwbs_variables

Module containing global variables relating to the first wall, blanket and shield components

References

life_blkt_fpy = None module-attribute

Full power blanket lifetime (years)

life_blkt = None module-attribute

Calendar year blanket lifetime (years)

m_fw_blkt_div_coolant_total = None module-attribute

mass of water coolant (in shield, blanket, first wall, divertor) [kg]

m_vv = None module-attribute

vacuum vessel mass [kg]

den_steel = None module-attribute

density of steel [kg m^-3]

denwc = None module-attribute

density of tungsten carbide [kg m^-3]

dewmkg = None module-attribute

total mass of vacuum vessel + cryostat [kg] (calculated if blktmodel>0)

f_p_blkt_multiplication = None module-attribute

energy multiplication in blanket and shield

p_blkt_multiplication_mw = None module-attribute

power due to energy multiplication in blanket and shield [MW]

fblss = None module-attribute

KIT blanket model: steel fraction of breeding zone

f_ster_div_single = None module-attribute

Solid angle fraction taken by one divertor

f_a_fw_outboard_hcd = None module-attribute

area fraction of first wall covered by heating/current drive apparatus plus diagnostics

fhole = None module-attribute

area fraction taken up by other holes (IFE)

i_fw_blkt_vv_shape = None module-attribute

switch for first wall, blanket, shield and vacuum vessel shape: - =1 D-shaped (cylinder inboard + ellipse outboard) - =2 defined by two ellipses

life_fw_fpy = None module-attribute

first wall full-power year lifetime (y)

m_fw_total = None module-attribute

first wall mass [kg]

fw_armour_mass = None module-attribute

first wall armour mass [kg]

fw_armour_thickness = None module-attribute

first wall armour thickness [m]

fw_armour_vol = None module-attribute

first wall armour volume [m^3]

i_blanket_type = None module-attribute

switch for blanket model: - =1 CCFE HCPB model - =2 KIT HCPB model # REMOVED, no longer usable - =3 CCFE HCPB model with Tritium Breeding Ratio calculation - =4 KIT HCLL model # REMOVED, no longer usable - =5 DCLL model - no nutronics model included (in development) please check/choose values for 'dual-coolant blanket' fractions (provided in this file). - please use i_p_coolant_pumping = 0 or 1.

i_blkt_inboard = None module-attribute

switch for inboard blanket: - =0 No inboard blanket (dr_blkt_inboard=0.0) - =1 Inboard blanket present

inuclear = None module-attribute

switch for nuclear heating in the coils: - =0 Frances Fox model (default) - =1 Fixed by user (qnuc)

qnuc = None module-attribute

nuclear heating in the coils (W) (inuclear=1)

f_blkt_li6_enrichment = None module-attribute

lithium-6 enrichment of breeding material (%)

p_blkt_nuclear_heat_total_mw = None module-attribute

nuclear heating in the blanket [MW]

pnuc_cp = None module-attribute

Total nuclear heating in the ST centrepost [MW]

p_cp_shield_nuclear_heat_mw = None module-attribute

Neutronic shield nuclear heating in the ST centrepost [MW]

pnuc_cp_tf = None module-attribute

TF neutronic nuclear heating in the ST centrepost [MW]

p_div_nuclear_heat_total_mw = None module-attribute

nuclear heating in the divertor [MW]

p_fw_nuclear_heat_total_mw = None module-attribute

nuclear heating in the first wall [MW]

p_fw_hcd_nuclear_heat_mw = None module-attribute

Nuclear heating in the HCD apparatus and diagnostics on the first wall [MW]

pnucloss = None module-attribute

nuclear heating lost via holes [MW]

pnucvvplus = None module-attribute

nuclear heating to vacuum vessel and beyond [MW]

p_shld_nuclear_heat_mw = None module-attribute

nuclear heating in the shield [MW]

m_blkt_total = None module-attribute

mass of blanket [kg]

m_blkt_steel_total = None module-attribute

mass of blanket - steel part [kg]

armour_fw_bl_mass = None module-attribute

Total mass of armour, first wall and blanket [kg]

breeder_f = None module-attribute

Volume ratio: Li4SiO4/(Be12Ti+Li4SiO4) (iteration variable 108)

breeder_multiplier = None module-attribute

combined breeder/multipler fraction of blanket by volume

vfcblkt = None module-attribute

He coolant fraction of blanket by volume (i_blanket_type= 1,3 (CCFE HCPB))

vfpblkt = None module-attribute

He purge gas fraction of blanket by volume (i_blanket_type= 1,3 (CCFE HCPB))

m_blkt_li4sio4 = None module-attribute

mass of lithium orthosilicate in blanket [kg] (i_blanket_type=1,3 (CCFE HCPB))

m_blkt_tibe12 = None module-attribute

mass of titanium beryllide in blanket [kg] (i_blanket_type=1,3 (CCFE HCPB))

neut_flux_cp = None module-attribute

Centrepost TF fast neutron flux (E > 0.1 MeV) [m(-2).(-1)] This variable is only calculated for superconducting (i_tf_sup = 1 ) spherical tokamal magnet designs (itart = 0)

f_neut_shield = None module-attribute

Fraction of nuclear power shielded before the CP magnet (ST) ( neut_absorb = -1 → a fit on simplified MCNP neutronic calculation is used assuming water cooled (13%) tungesten carbyde )

f_a_fw_coolant_inboard = None module-attribute

Inboard FW coolant cross-sectional area void fraction

f_a_fw_coolant_outboard = None module-attribute

Outboard FW coolant cross-sectional area void fraction

psurffwi = None module-attribute

Surface heat flux on first wall [MW] (sum = p_fw_rad_total_mw)

psurffwo = None module-attribute

Surface heat flux on first wall [MW] (sum = p_fw_rad_total_mw)

vol_fw_total = None module-attribute

First wall volume [m3]

f_vol_blkt_steel = None module-attribute

Fractions of blanket by volume: steel

f_vol_blkt_li4sio4 = None module-attribute

Fractions of blanket by volume: lithium orthosilicate

f_vol_blkt_tibe12 = None module-attribute

Fractions of blanket by volume: titanium beryllide

breedmat = None module-attribute

breeder material switch (i_blanket_type=2 (KIT HCPB)): - =1 Lithium orthosilicate - =2 Lithium methatitanate - =3 Lithium zirconate

densbreed = None module-attribute

density of breeder material [kg m^-3] (i_blanket_type=2 (KIT HCPB))

fblbe = None module-attribute

beryllium fraction of blanket by volume (if i_blanket_type=2, is Be fraction of breeding zone)

fblbreed = None module-attribute

breeder fraction of blanket breeding zone by volume (i_blanket_type=2 (KIT HCPB))

fblhebmi = None module-attribute

helium fraction of inboard blanket box manifold by volume (i_blanket_type=2 (KIT HCPB))

fblhebmo = None module-attribute

helium fraction of outboard blanket box manifold by volume (i_blanket_type=2 (KIT HCPB))

fblhebpi = None module-attribute

helium fraction of inboard blanket back plate by volume (i_blanket_type=2 (KIT HCPB))

fblhebpo = None module-attribute

helium fraction of outboard blanket back plate by volume (i_blanket_type=2 (KIT HCPB))

hcdportsize = None module-attribute

switch for size of heating/current drive ports (i_blanket_type=2 (KIT HCPB)): - =1 'small' - =2 'large'

nflutf = None module-attribute

peak fast neutron fluence on TF coil superconductor [n m^-2] (i_blanket_type=2 (KIT HCPB))

npdiv = None module-attribute

number of divertor ports (i_blanket_type=2 (KIT HCPB))

nphcdin = None module-attribute

number of inboard ports for heating/current drive (i_blanket_type=2 (KIT HCPB))

nphcdout = None module-attribute

number of outboard ports for heating/current drive (i_blanket_type=2 (KIT HCPB))

tbr = None module-attribute

tritium breeding ratio (i_blanket_type=2,3 (KIT HCPB/HCLL))

tritprate = None module-attribute

tritium production rate [g day^-1] (i_blanket_type=2 (KIT HCPB))

wallpf = None module-attribute

neutron wall load peaking factor (i_blanket_type=2 (KIT HCPB))

whtblbreed = None module-attribute

mass of blanket - breeder part [kg] (i_blanket_type=2 (KIT HCPB))

m_blkt_beryllium = None module-attribute

mass of blanket - beryllium part [kg]

i_p_coolant_pumping = None module-attribute

Switch for pumping power for primary coolant (mechanical power only and peak first wall temperature is only calculated if i_p_coolant_pumping=2): - =0 User sets pump power directly (p_blkt_coolant_pump_mw, p_fw_coolant_pump_mw, p_div_coolant_pump_mw, p_shld_coolant_pump_mw) - =1 User sets pump power as a fraction of thermal power (f_p_blkt_coolant_pump_total_heat, f_p_fw_coolant_pump_total_heat, f_p_div_coolant_pump_total_heat, f_p_shld_coolant_pump_total_heat) - =2 Mechanical pumping power is calculated - =3 Mechanical pumping power is calculated using specified pressure drop

i_shield_mat = None module-attribute

Switch for shield material - currently only applied in costing routines cost_model = 2 - =0 Tungsten (default) - =1 Tungsten carbide

i_thermal_electric_conversion = None module-attribute

Switch for power conversion cycle: - =0 Set efficiency for chosen blanket, from detailed models (divertor heat not used) - =1 Set efficiency for chosen blanket, from detailed models (divertor heat used) - =2 user input thermal-electric efficiency (eta_turbine) - =3 steam Rankine cycle - =4 supercritical CO2 cycle

secondary_cycle_liq = None module-attribute

Switch for power conversion cycle for the liquid breeder component of the blanket: - =2 user input thermal-electric efficiency (eta_turbine) - =4 supercritical CO2 cycle

i_blkt_coolant_type = None module-attribute

Switch for blanket coolant (set via blkttype): - =1 helium - =2 pressurized water

i_fw_coolant_type = None module-attribute

switch for first wall coolant (can be different from blanket coolant): - 'helium' - 'water'

dr_fw_wall = None module-attribute

wall thickness of first wall coolant channels [m]

radius_fw_channel = None module-attribute

radius of first wall cooling channels [m]

dx_fw_module = None module-attribute

Width of a FW module containing a cooling channel [m]

temp_fw_coolant_in = None module-attribute

inlet temperature of first wall coolant [K]

temp_fw_coolant_out = None module-attribute

outlet temperature of first wall coolant [K]

pres_fw_coolant = None module-attribute

first wall coolant pressure [Pa] (i_thermal_electric_conversion>1)

temp_fw_peak = None module-attribute

peak first wall temperature [K]

roughness_fw_channel = None module-attribute

first wall channel roughness epsilon [m]

len_fw_channel = None module-attribute

Length of a single first wall channel (all in parallel) [m] (iteration variable 114, useful for constraint equation 39)

f_fw_peak = None module-attribute

peaking factor for first wall heat loads. (Applied separately to inboard and outboard loads. Applies to both neutron and surface loads. Only used to calculate peak temperature - not the coolant flow rate.)

pres_blkt_coolant = None module-attribute

blanket coolant pressure [Pa] (i_thermal_electric_conversion>1)

temp_blkt_coolant_in = None module-attribute

inlet temperature of blanket coolant [K] (i_thermal_electric_conversion>1)

temp_blkt_coolant_out = None module-attribute

Outlet temperature of blanket coolant [K] (i_thermal_electric_conversion>1) - input if i_blkt_coolant_type=1 (helium) - calculated if i_blkt_coolant_type=2 (water)

coolp = None module-attribute

blanket coolant pressure [Pa] (stellarator only)

n_blkt_outboard_modules_poloidal = None module-attribute

number of outboard blanket modules in poloidal direction (i_thermal_electric_conversion>1)

n_blkt_inboard_modules_poloidal = None module-attribute

number of inboard blanket modules in poloidal direction (i_thermal_electric_conversion>1)

n_blkt_outboard_modules_toroidal = None module-attribute

number of outboard blanket modules in toroidal direction (i_thermal_electric_conversion>1)

n_blkt_inboard_modules_toroidal = None module-attribute

number of inboard blanket modules in toroidal direction (i_thermal_electric_conversion>1)

temp_fw_max = None module-attribute

maximum temperature of first wall material [K] (i_thermal_electric_conversion>1)

fw_th_conductivity = None module-attribute

thermal conductivity of first wall material at 293 K (W/m/K) (Temperature dependence is as for unirradiated Eurofer)

fvoldw = None module-attribute

area coverage factor for vacuum vessel volume

fvolsi = None module-attribute

area coverage factor for inboard shield volume

fvolso = None module-attribute

area coverage factor for outboard shield volume

fwclfr = None module-attribute

first wall coolant fraction (calculated if i_pulsed_plant=1 or ipowerflow=1)

p_div_rad_total_mw = None module-attribute

Total radiation power incident on the divertor(s) (MW)

p_fw_rad_total_mw = None module-attribute

Radiation power incident on the first wall (MW)

p_fw_hcd_rad_total_mw = None module-attribute

Radiation power incident on the heating and current drive systems on the first wall (MW)

pradloss = None module-attribute

Radiation power lost through holes (eventually hits shield) (MW) Only used for stellarator

p_tf_nuclear_heat_mw = None module-attribute

nuclear heating in the TF coil (MW)

ptfnucpm3 = None module-attribute

nuclear heating in the TF coil (MW/m3) (blktmodel>0)

r_cryostat_inboard = None module-attribute

cryostat radius [m]

z_cryostat_half_inside = None module-attribute

cryostat height [m]

dr_pf_cryostat = None module-attribute

Radial distance between outer edge of furthest away PF coil (or stellarator modular coil) and cryostat [m]

vol_cryostat = None module-attribute

Cryostat structure volume [m^3]

vol_cryostat_internal = None module-attribute

Internal volume of the cryostat [m^3]

vol_vv = None module-attribute

vacuum vessel volume [m^3]

vfshld = None module-attribute

coolant void fraction in shield

vol_blkt_total = None module-attribute

volume of blanket [m^3]

vol_blkt_total_full_coverage = None module-attribute

Volume of blanket with no holes or ports (toroidally continuous) [m³]

vol_blkt_inboard = None module-attribute

volume of inboard blanket [m^3]

vol_blkt_inboard_full_coverage = None module-attribute

Volume of inboard blanket with no holes or ports (toroidally continuous) [m³]

vol_blkt_outboard = None module-attribute

volume of outboard blanket [m^3]

vol_blkt_outboard_full_coverage = None module-attribute

Volume of outboard blanket with no holes or ports (toroidally continuous) [m³]

vol_shld_total = None module-attribute

volume of shield [m^3]

whtshld = None module-attribute

mass of shield [kg]

wpenshld = None module-attribute

mass of the penetration shield [kg]

wtshldi = None module-attribute

mass of inboard shield [kg]

wtshldo = None module-attribute

mass of outboard shield [kg]

irefprop = None module-attribute

Switch to use REFPROP routines (stellarator only)

fblli = None module-attribute

lithium fraction of blanket by volume (stellarator only)

fblli2o = None module-attribute

lithium oxide fraction of blanket by volume (stellarator only)

fbllipb = None module-attribute

lithium lead fraction of blanket by volume (stellarator only)

fblvd = None module-attribute

vanadium fraction of blanket by volume (stellarator only)

m_blkt_li2o = None module-attribute

mass of blanket - Li_2O part [kg]

wtbllipb = None module-attribute

mass of blanket - Li-Pb part [kg]

m_blkt_vanadium = None module-attribute

mass of blanket - vanadium part [kg]

m_blkt_lithium = None module-attribute

mass of blanket - lithium part [kg]

f_a_blkt_cooling_channels = None module-attribute

coolant void fraction in blanket.

blktmodel = None module-attribute

switch for blanket/tritium breeding model (see i_blanket_type): - =0 original simple model - =1 KIT model based on a helium-cooled pebble-bed blanket (HCPB) reference design

declblkt = None module-attribute

neutron power deposition decay length of blanket structural material [m] (stellarators only)

declfw = None module-attribute

neutron power deposition decay length of first wall structural material [m] (stellarators only)

declshld = None module-attribute

neutron power deposition decay length of shield structural material [m] (stellarators only)

blkttype = None module-attribute

Switch for blanket type: - =1 WCLL; - =2 HCLL; efficiency taken from M. Kovari 2016 "PROCESS": A systems code for fusion power plants - Part 2: Engineering https://www.sciencedirect.com/science/article/pii/S0920379616300072 Feedheat & reheat cycle assumed - =3 HCPB; efficiency taken from M. Kovari 2016 "PROCESS": A systems code for fusion power plants - Part 2: Engineering https://www.sciencedirect.com/science/article/pii/S0920379616300072 Feedheat & reheat cycle assumed

etaiso = None module-attribute

isentropic efficiency of FW and blanket coolant pumps

eta_coolant_pump_electric = None module-attribute

electrical efficiency of primary coolant pumps

i_fw_blkt_shared_coolant = None module-attribute

Switch for whether the FW and BB are on the same pump system i.e. do they have the same primary coolant or not - =0 FW and BB have the same primary coolant, flow = FWin->FWout->BBin->BBout - =1 FW and BB have the different primary coolant and are on different pump systems

i_blkt_liquid_breeder_type = None module-attribute

Switch for Liquid Metal Breeder Material - =0 PbLi - =1 Li

i_blkt_dual_coolant = None module-attribute

Switch to specify whether breeding blanket is single-cooled or dual-coolant. - =0 Single coolant used for FW and Blanket (H2O or He). Solid Breeder. - =1 Single coolant used for FW and Blanket (H2O or He). Liquid metal breeder circulted for tritium extraction. - =2 Dual coolant: primary coolant (H2O or He) for FW and blanket structure; secondary coolant is self-cooled liquid metal breeder.

i_blkt_liquid_breeder_channel_type = None module-attribute

Switch for Flow Channel Insert (FCI) type if liquid metal breeder blanket. - =0 Thin conducting walls, default electrical conductivity (bz_channel_conduct_liq) is Eurofer - =1 Insulating Material, assumed perfect electrical insulator, default density (den_ceramic) is for SiC - =2 Insulating Material, electrical conductivity (bz_channel_conduct_liq) is input (default Eurofer), default density (den_ceramic) is for SiC

i_blkt_module_segmentation = None module-attribute

Switch for Multi Module Segment (MMS) or Single Modle Segment (SMS) - =0 MMS - =1 SMS

n_liq_recirc = None module-attribute

Number of liquid metal breeder recirculations per day, for use with i_blkt_dual_coolant=1

r_f_liq_ib = None module-attribute

Radial fraction of BZ liquid channels

r_f_liq_ob = None module-attribute

Radial fraction of BZ liquid channels

w_f_liq_ib = None module-attribute

Toroidal fraction of BZ liquid channels

w_f_liq_ob = None module-attribute

Toroidal fraction of BZ liquid channels

den_ceramic = None module-attribute

FCI material density

th_wall_secondary = None module-attribute

Liquid metal coolant/breeder wall thickness thin conductor or FCI [m]

bz_channel_conduct_liq = None module-attribute

Liquid metal coolant/breeder thin conductor or FCI wall conductance [A V^-1 m^-1]

a_bz_liq = None module-attribute

Toroidal width of the rectangular cooling channel [m] for long poloidal sections of blanket breeding zone

b_bz_liq = None module-attribute

Radial width of the rectangular cooling channel [m] for long poloidal sections of blanket breeding zone

nopol = None module-attribute

Number of poloidal sections in a liquid metal breeder/coolant channel for module/segment

nopipes = None module-attribute

Number of Liquid metal breeder/coolant channels per module/segment

den_liq = None module-attribute

Liquid metal breeder/coolant density [kg m^-3]

wht_liq = None module-attribute

Liquid metal

wht_liq_ib = None module-attribute

Liquid metal

wht_liq_ob = None module-attribute

Liquid metal

specific_heat_liq = None module-attribute

Liquid metal breeder/coolant specific heat [J kg^-1 K^-1]

thermal_conductivity_liq = None module-attribute

Liquid metal breeder/coolant thermal conductivity [W m^-1 K^-1]

dynamic_viscosity_liq = None module-attribute

Liquid metal breeder/coolant dynamic viscosity [Pa s]

electrical_conductivity_liq = None module-attribute

Liquid metal breeder/coolant electrical conductivity [Ohm m]

hartmann_liq = None module-attribute

Hartmann number

b_mag_blkt = None module-attribute

Toroidal Magnetic feild strength for IB/OB blanket [T]

etaiso_liq = None module-attribute

Isentropic efficiency of blanket liquid breeder/coolant pumps

blpressure_liq = None module-attribute

blanket liquid metal breeder/coolant pressure [Pa]

inlet_temp_liq = None module-attribute

Inlet (scan var 68) temperature of the liquid breeder/coolant [K]

outlet_temp_liq = None module-attribute

Outlet (scan var 69) temperature of the liquid breeder/coolant [K]

den_fw_coolant = None module-attribute

Density of the FW primary coolant

visc_fw_coolant = None module-attribute

Viscosity of the FW primary coolant

den_blkt_coolant = None module-attribute

Density of the blanket primary coolant

visc_blkt_coolant = None module-attribute

Viscosity of the blanket primary coolant

cp_fw = None module-attribute

Spesific heat for FW and blanket primary coolant(s)

cv_fw = None module-attribute

Spesific heat for FW and blanket primary coolant(s)

cp_bl = None module-attribute

Spesific heat for FW and blanket primary coolant(s)

cv_bl = None module-attribute

Spesific heat for FW and blanket primary coolant(s)

f_nuc_pow_bz_struct = None module-attribute

For a dual-coolant blanket, fraction of BZ power cooled by primary coolant

f_nuc_pow_bz_liq = None module-attribute

For a dual-coolant blanket, fraction of BZ self-cooled power (secondary coolant)

pnuc_fw_ratio_dcll = None module-attribute

For a dual-coolant blanket, ratio of FW nuclear power as fraction of total

pnuc_blkt_ratio_dcll = None module-attribute

For a dual-coolant blanket, ratio of Blanket nuclear power as fraction of total

n_blkt_inboard_module_coolant_sections_radial = None module-attribute

Number of radial and poloidal sections that make up the total primary coolant flow length in a blanket module (IB and OB)

n_blkt_inboard_module_coolant_sections_poloidal = None module-attribute

Number of radial and poloidal sections that make up the total primary coolant flow length in a blanket module (IB and OB)

n_blkt_outboard_module_coolant_sections_radial = None module-attribute

Number of radial and poloidal sections that make up the total primary coolant flow length in a blanket module (IB and OB)

n_blkt_outboard_module_coolant_sections_poloidal = None module-attribute

Number of radial and poloidal sections that make up the total primary coolant flow length in a blanket module (IB and OB)

bzfllengi_n_rad_liq = None module-attribute

Number of radial and poloidal sections that make up the total secondary coolant/breeder flow length in a blanket module (IB and OB)

bzfllengi_n_pol_liq = None module-attribute

Number of radial and poloidal sections that make up the total secondary coolant/breeder flow length in a blanket module (IB and OB)

bzfllengo_n_rad_liq = None module-attribute

Number of radial and poloidal sections that make up the total secondary coolant/breeder flow length in a blanket module (IB and OB)

bzfllengo_n_pol_liq = None module-attribute

Number of radial and poloidal sections that make up the total secondary coolant/breeder flow length in a blanket module (IB and OB)

radius_blkt_channel = None module-attribute

Radius of blanket cooling channels [m]

radius_blkt_channel_90_bend = None module-attribute

Radius of blanket cooling channel 90° bend [m]

radius_fw_channel_90_bend = None module-attribute

Radius of first wall cooling channel 90° bend [m]

radius_fw_channel_180_bend = None module-attribute

Radius of first wall cooling channel 180° bend [m]

radius_blkt_channel_180_bend = None module-attribute

Radius of blanket cooling channel 180° bend [m]

dz_fw_half = None module-attribute

Half-height of first wall structure [m]

init_fwbs_variables()

Initialise FWBS variables

Source code in process/data_structure/fwbs_variables.py
 914
 915
 916
 917
 918
 919
 920
 921
 922
 923
 924
 925
 926
 927
 928
 929
 930
 931
 932
 933
 934
 935
 936
 937
 938
 939
 940
 941
 942
 943
 944
 945
 946
 947
 948
 949
 950
 951
 952
 953
 954
 955
 956
 957
 958
 959
 960
 961
 962
 963
 964
 965
 966
 967
 968
 969
 970
 971
 972
 973
 974
 975
 976
 977
 978
 979
 980
 981
 982
 983
 984
 985
 986
 987
 988
 989
 990
 991
 992
 993
 994
 995
 996
 997
 998
 999
1000
1001
1002
1003
1004
1005
1006
1007
1008
1009
1010
1011
1012
1013
1014
1015
1016
1017
1018
1019
1020
1021
1022
1023
1024
1025
1026
1027
1028
1029
1030
1031
1032
1033
1034
1035
1036
1037
1038
1039
1040
1041
1042
1043
1044
1045
1046
1047
1048
1049
1050
1051
1052
1053
1054
1055
1056
1057
1058
1059
1060
1061
1062
1063
1064
1065
1066
1067
1068
1069
1070
1071
1072
1073
1074
1075
1076
1077
1078
1079
1080
1081
1082
1083
1084
1085
1086
1087
1088
1089
1090
1091
1092
1093
1094
1095
1096
1097
1098
1099
1100
1101
1102
1103
1104
1105
1106
1107
1108
1109
1110
1111
1112
1113
1114
1115
1116
1117
1118
1119
1120
1121
1122
1123
1124
1125
1126
1127
1128
1129
1130
1131
1132
1133
1134
1135
1136
1137
1138
1139
1140
1141
1142
1143
1144
1145
1146
1147
1148
1149
1150
1151
1152
1153
1154
1155
1156
1157
1158
1159
1160
1161
1162
1163
1164
1165
1166
1167
1168
1169
1170
1171
1172
1173
1174
1175
1176
1177
1178
1179
1180
1181
1182
1183
1184
1185
1186
1187
1188
1189
1190
1191
1192
1193
1194
1195
1196
1197
1198
1199
1200
1201
1202
1203
1204
1205
1206
1207
1208
1209
1210
1211
1212
1213
1214
1215
1216
1217
1218
1219
1220
1221
1222
1223
1224
1225
1226
1227
1228
1229
1230
1231
1232
1233
1234
1235
1236
1237
1238
1239
1240
1241
1242
1243
1244
1245
1246
1247
1248
1249
1250
1251
1252
1253
1254
1255
1256
1257
1258
1259
1260
1261
1262
1263
1264
1265
1266
1267
1268
1269
1270
1271
1272
1273
1274
1275
1276
1277
1278
1279
1280
1281
1282
1283
1284
1285
1286
1287
1288
1289
1290
1291
1292
1293
1294
1295
1296
1297
1298
1299
1300
1301
1302
1303
1304
1305
1306
1307
1308
1309
1310
1311
1312
1313
1314
1315
def init_fwbs_variables():
    """Initialise FWBS variables"""
    global \
        life_blkt_fpy, \
        life_blkt, \
        m_fw_blkt_div_coolant_total, \
        m_vv, \
        den_steel, \
        denwc, \
        dewmkg, \
        f_p_blkt_multiplication, \
        p_blkt_multiplication_mw, \
        fblss, \
        f_ster_div_single, \
        f_a_fw_outboard_hcd, \
        fhole, \
        i_fw_blkt_vv_shape, \
        life_fw_fpy, \
        m_fw_total, \
        fw_armour_mass, \
        fw_armour_thickness, \
        fw_armour_vol, \
        i_blanket_type, \
        i_blkt_inboard, \
        inuclear, \
        qnuc, \
        f_blkt_li6_enrichment, \
        p_blkt_nuclear_heat_total_mw, \
        pnuc_cp, \
        p_cp_shield_nuclear_heat_mw, \
        pnuc_cp_tf, \
        p_div_nuclear_heat_total_mw, \
        p_fw_nuclear_heat_total_mw, \
        p_fw_hcd_nuclear_heat_mw, \
        pnucloss, \
        pnucvvplus, \
        p_shld_nuclear_heat_mw, \
        m_blkt_total, \
        m_blkt_steel_total, \
        armour_fw_bl_mass, \
        breeder_f, \
        breeder_multiplier, \
        vfcblkt, \
        vfpblkt, \
        m_blkt_li4sio4, \
        m_blkt_tibe12, \
        neut_flux_cp, \
        f_neut_shield, \
        f_a_fw_coolant_inboard, \
        f_a_fw_coolant_outboard, \
        psurffwi, \
        psurffwo, \
        vol_fw_total, \
        f_vol_blkt_steel, \
        f_vol_blkt_li4sio4, \
        f_vol_blkt_tibe12, \
        breedmat, \
        densbreed, \
        fblbe, \
        fblbreed, \
        fblhebmi, \
        fblhebmo, \
        fblhebpi, \
        fblhebpo, \
        hcdportsize, \
        nflutf, \
        npdiv, \
        nphcdin, \
        nphcdout, \
        tbr, \
        tritprate, \
        wallpf, \
        whtblbreed, \
        m_blkt_beryllium, \
        i_p_coolant_pumping, \
        i_shield_mat, \
        i_thermal_electric_conversion, \
        secondary_cycle_liq, \
        i_blkt_coolant_type, \
        i_fw_coolant_type, \
        dr_fw_wall, \
        radius_fw_channel, \
        dx_fw_module, \
        temp_fw_coolant_in, \
        temp_fw_coolant_out, \
        pres_fw_coolant, \
        temp_fw_peak, \
        roughness_fw_channel, \
        len_fw_channel, \
        f_fw_peak, \
        pres_blkt_coolant, \
        temp_blkt_coolant_in, \
        temp_blkt_coolant_out, \
        coolp, \
        n_blkt_outboard_modules_poloidal, \
        n_blkt_inboard_modules_poloidal, \
        n_blkt_outboard_modules_toroidal, \
        n_blkt_inboard_modules_toroidal, \
        temp_fw_max, \
        fw_th_conductivity, \
        fvoldw, \
        fvolsi, \
        fvolso, \
        fwclfr, \
        p_div_rad_total_mw, \
        p_fw_rad_total_mw, \
        p_fw_hcd_rad_total_mw, \
        pradloss, \
        p_tf_nuclear_heat_mw, \
        ptfnucpm3, \
        r_cryostat_inboard, \
        z_cryostat_half_inside, \
        dr_pf_cryostat, \
        vol_cryostat, \
        vol_cryostat_internal, \
        vol_vv, \
        vfshld, \
        vol_blkt_total, \
        vol_blkt_total_full_coverage, \
        vol_blkt_inboard, \
        vol_blkt_inboard_full_coverage, \
        vol_blkt_outboard, \
        vol_blkt_outboard_full_coverage, \
        vol_shld_total, \
        whtshld, \
        wpenshld, \
        wtshldi, \
        wtshldo, \
        irefprop, \
        fblli, \
        fblli2o, \
        fbllipb, \
        fblvd, \
        m_blkt_li2o, \
        wtbllipb, \
        m_blkt_vanadium, \
        m_blkt_lithium, \
        f_a_blkt_cooling_channels, \
        blktmodel, \
        declblkt, \
        declfw, \
        declshld, \
        blkttype, \
        etaiso, \
        eta_coolant_pump_electric, \
        i_fw_blkt_shared_coolant, \
        i_blkt_liquid_breeder_type, \
        i_blkt_dual_coolant, \
        i_blkt_liquid_breeder_channel_type, \
        i_blkt_module_segmentation, \
        n_liq_recirc, \
        r_f_liq_ib, \
        r_f_liq_ob, \
        w_f_liq_ib, \
        w_f_liq_ob, \
        den_ceramic, \
        th_wall_secondary, \
        bz_channel_conduct_liq, \
        a_bz_liq, \
        b_bz_liq, \
        nopol, \
        nopipes, \
        den_liq, \
        wht_liq, \
        wht_liq_ib, \
        wht_liq_ob, \
        specific_heat_liq, \
        thermal_conductivity_liq, \
        dynamic_viscosity_liq, \
        electrical_conductivity_liq, \
        hartmann_liq, \
        b_mag_blkt, \
        etaiso_liq, \
        blpressure_liq, \
        inlet_temp_liq, \
        outlet_temp_liq, \
        den_fw_coolant, \
        visc_fw_coolant, \
        den_blkt_coolant, \
        visc_blkt_coolant, \
        cp_fw, \
        cv_fw, \
        cp_bl, \
        cv_bl, \
        f_nuc_pow_bz_struct, \
        f_nuc_pow_bz_liq, \
        pnuc_fw_ratio_dcll, \
        pnuc_blkt_ratio_dcll, \
        n_blkt_inboard_module_coolant_sections_radial, \
        n_blkt_inboard_module_coolant_sections_poloidal, \
        n_blkt_outboard_module_coolant_sections_radial, \
        n_blkt_outboard_module_coolant_sections_poloidal, \
        bzfllengi_n_rad_liq, \
        bzfllengi_n_pol_liq, \
        bzfllengo_n_rad_liq, \
        bzfllengo_n_pol_liq, \
        radius_blkt_channel, \
        radius_fw_channel_90_bend, \
        radius_fw_channel_180_bend, \
        radius_blkt_channel_90_bend, \
        radius_blkt_channel_180_bend, \
        dz_fw_half

    life_blkt_fpy = 0.0
    life_blkt = 0.0
    m_fw_blkt_div_coolant_total = 0.0
    m_vv = 0.0
    den_steel = 7800.0
    denwc = 15630.0
    dewmkg = 0.0
    f_p_blkt_multiplication = 1.269
    p_blkt_multiplication_mw = 0.0
    fblss = 0.09705
    f_ster_div_single = 0.115
    f_a_fw_outboard_hcd = 0.0
    fhole = 0.0
    i_fw_blkt_vv_shape = 2
    life_fw_fpy = 0.0
    m_fw_total = 0.0
    fw_armour_mass = 0.0
    fw_armour_thickness = 0.005
    fw_armour_vol = 0.0
    i_blanket_type = 1
    i_blkt_inboard = 1
    inuclear = 0
    qnuc = 0.0
    f_blkt_li6_enrichment = 30.0
    p_blkt_nuclear_heat_total_mw = 0.0
    p_div_nuclear_heat_total_mw = 0.0
    p_fw_nuclear_heat_total_mw = 0.0
    p_fw_hcd_nuclear_heat_mw = 0.0
    pnucloss = 0.0
    pnucvvplus = 0.0
    p_shld_nuclear_heat_mw = 0.0
    m_blkt_total = 0.0
    m_blkt_steel_total = 0.0
    armour_fw_bl_mass = 0.0
    breeder_f = 0.5
    breeder_multiplier = 0.75
    vfcblkt = 0.05295
    vfpblkt = 0.1
    m_blkt_li4sio4 = 0.0
    m_blkt_tibe12 = 0.0
    f_neut_shield = -1.0
    f_a_fw_coolant_inboard = 0.0
    f_a_fw_coolant_outboard = 0.0
    psurffwi = 0.0
    psurffwo = 0.0
    vol_fw_total = 0.0
    f_vol_blkt_steel = 0.0
    f_vol_blkt_li4sio4 = 0.0
    f_vol_blkt_tibe12 = 0.0
    breedmat = 1
    densbreed = 0.0
    fblbe = 0.6
    fblbreed = 0.154
    fblhebmi = 0.4
    fblhebmo = 0.4
    fblhebpi = 0.6595
    fblhebpo = 0.6713
    hcdportsize = 1
    nflutf = 0.0
    npdiv = 2
    nphcdin = 2
    nphcdout = 2
    tbr = 0.0
    tritprate = 0.0
    wallpf = 1.21
    whtblbreed = 0.0
    m_blkt_beryllium = 0.0
    i_p_coolant_pumping = 2
    i_shield_mat = 0
    i_thermal_electric_conversion = 0
    secondary_cycle_liq = 4
    i_blkt_coolant_type = 1
    i_fw_coolant_type = "helium"
    dr_fw_wall = 0.003
    radius_fw_channel = 0.006
    dx_fw_module = 0.02
    temp_fw_coolant_in = 573.0
    temp_fw_coolant_out = 823.0
    pres_fw_coolant = 15.5e6
    temp_fw_peak = 873.0
    roughness_fw_channel = 1.0e-6
    len_fw_channel = 4.0
    f_fw_peak = 1.0
    pres_blkt_coolant = 15.50e6
    temp_blkt_coolant_in = 573.0
    temp_blkt_coolant_out = 823.0
    coolp = 15.5e6
    n_blkt_outboard_modules_poloidal = 8
    n_blkt_inboard_modules_poloidal = 7
    n_blkt_outboard_modules_toroidal = 48
    n_blkt_inboard_modules_toroidal = 32
    temp_fw_max = 823.0
    fw_th_conductivity = 28.34
    fvoldw = 1.74
    fvolsi = 1.0
    fvolso = 0.64
    fwclfr = 0.15
    p_div_rad_total_mw = 0.0
    p_fw_rad_total_mw = 0.0
    p_fw_hcd_rad_total_mw = 0.0
    pradloss = 0.0
    p_tf_nuclear_heat_mw = 0.0
    ptfnucpm3 = 0.0
    r_cryostat_inboard = 0.0
    z_cryostat_half_inside = 0.0
    dr_pf_cryostat = 0.5
    vol_cryostat = 0.0
    vol_cryostat_internal = 0.0
    vol_vv = 0.0
    vfshld = 0.25
    vol_blkt_total = 0.0
    vol_blkt_total_full_coverage = 0.0
    vol_blkt_inboard = 0.0
    vol_blkt_inboard_full_coverage = 0.0
    vol_blkt_outboard = 0.0
    vol_blkt_outboard_full_coverage = 0.0
    vol_shld_total = 0.0
    whtshld = 0.0
    wpenshld = 0.0
    wtshldi = 0.0
    wtshldo = 0.0
    irefprop = 1
    fblli = 0.0
    fblli2o = 0.08
    fbllipb = 0.68
    fblvd = 0.0
    m_blkt_li2o = 0.0
    wtbllipb = 0.0
    m_blkt_vanadium = 0.0
    m_blkt_lithium = 0.0
    f_a_blkt_cooling_channels = 0.25
    blktmodel = 0
    declblkt = 0.075
    declfw = 0.075
    declshld = 0.075
    blkttype = 3
    etaiso = 0.85
    eta_coolant_pump_electric = 0.95
    pnuc_cp = 0.0
    p_cp_shield_nuclear_heat_mw = 0.0
    pnuc_cp_tf = 0.0
    neut_flux_cp = 0.0
    i_fw_blkt_shared_coolant = 0
    i_blkt_liquid_breeder_type = 0
    i_blkt_dual_coolant = 0
    i_blkt_liquid_breeder_channel_type = 0
    i_blkt_module_segmentation = 0
    n_liq_recirc = 10
    r_f_liq_ib = 0.5
    r_f_liq_ob = 0.5
    w_f_liq_ib = 0.5
    w_f_liq_ob = 0.5
    den_ceramic = 3.21e3
    th_wall_secondary = 1.25e-2
    bz_channel_conduct_liq = 8.33e5
    a_bz_liq = 0.2
    b_bz_liq = 0.2
    nopol = 2
    nopipes = 4
    den_liq = 9.5e3
    specific_heat_liq = 1.9e2
    thermal_conductivity_liq = 30.0
    wht_liq = 0.0
    wht_liq_ib = 0.0
    wht_liq_ob = 0.0
    dynamic_viscosity_liq = 0.0
    electrical_conductivity_liq = 0.0
    hartmann_liq = [0.0, 0.0]
    b_mag_blkt = [5.0, 5.0]
    etaiso_liq = 0.85
    blpressure_liq = 1.7e6
    inlet_temp_liq = 570.0
    outlet_temp_liq = 720.0
    den_fw_coolant = 0.0
    visc_fw_coolant = 0.0
    den_blkt_coolant = 0.0
    visc_blkt_coolant = 0.0
    cp_fw = 0.0
    cv_fw = 0.0
    cp_bl = 0.0
    cv_bl = 0.0
    f_nuc_pow_bz_struct = 0.34
    f_nuc_pow_bz_liq = 0.66
    pnuc_fw_ratio_dcll = 0.14
    pnuc_blkt_ratio_dcll = 0.86
    n_blkt_inboard_module_coolant_sections_radial = 4
    n_blkt_inboard_module_coolant_sections_poloidal = 2
    n_blkt_outboard_module_coolant_sections_radial = 4
    n_blkt_outboard_module_coolant_sections_poloidal = 2
    bzfllengi_n_rad_liq = 2
    bzfllengi_n_pol_liq = 2
    bzfllengo_n_rad_liq = 2
    bzfllengo_n_pol_liq = 2
    radius_blkt_channel = 0.0
    radius_fw_channel_90_bend = 0.0
    radius_fw_channel_180_bend = 0.0
    radius_blkt_channel_90_bend = 0.0
    radius_blkt_channel_180_bend = 0.0
    dz_fw_half = 0.0