diff --git a/process/geometry/__init__.py b/process/geometry/__init__.py new file mode 100644 index 0000000000..e69de29bb2 diff --git a/process/geometry/blanket_geometry.py b/process/geometry/blanket_geometry.py new file mode 100644 index 0000000000..9256ccfec2 --- /dev/null +++ b/process/geometry/blanket_geometry.py @@ -0,0 +1,216 @@ +""" +Calculate radial and vertical coordinates for the geometry of the blanket +""" +from typing import Tuple +import numpy as np +from process.geometry.utils import dh_vertices, dhgap_vertices +from process.geometry.geometry_parameterisations import ArbitraryGeometry + + +def blanket_geometry_single_null( + radx_outer: float, + rminx_outer: float, + radx_inner: float, + rminx_inner: float, + cumulative_upper: dict, + triang: float, + cumulative_lower: dict, + blnktth: float, + c_shldith: float, + c_blnkoth: float, + blnkith: float, + blnkoth: float, +) -> ArbitraryGeometry: + """Calculates radial and vertical distances for the geometry of blanket in a single null configuration + + :param radx_outer: outboard radius of outer surface of blanket + :type radx_outer: float + :param rminx_outer: inboard radius of outer surface of blanket + :type rminx_outer: float + :param radx_inner: outboard radius of inner surface of blanket + :type radx_inner: float + :param rminx_inner: inboard radius of inner surface of blanket + :type rminx_inner: float + :param cumulative_upper: cumulative vertical thicknesses of components above the midplane + :type cumulative_upper: dict + :param triang: plasma triangularity + :type triang: float + :param cumulative_lower: cumulative vertical thicknesses of components below the midplane + :type cumulative_lower: dict + :param blnktth: top blanket vertical thickness + :type blnktth: float + :param c_shldith: inboard shield thickness + :type c_shldith: float + :param c_blnkoth: outboard blanket radial thickness + :type c_blnkoth: float + :param blnkith: inboard blanket radial thickness + :type blnkith: float + :param blnkoth: outboard blanket radial thickness + :type blnkoth: float + :return: dataclass returning radial and vertical coordinates + :rtype: ArbitraryGeometry + """ + # Upper blanket outer surface + kapx = cumulative_upper["blnktth"] / rminx_outer + rs_upper_outboard, zs_upper_outboard = dh_vertices( + radx_outer, rminx_outer, triang, kapx + ) + + # Upper blanket inner surface + kapx = cumulative_upper["fwtth"] / rminx_inner + rs_upper_inboard, zs_upper_inboard = dh_vertices( + radx_inner, rminx_inner, triang, kapx + ) + + # Lower blanket + divgap = cumulative_lower["divfix"] + + ( + rs_lower_outboard, + zs_lower_outboard, + rs_lower_inboard, + zs_lower_inboard, + ) = blanket_geometry_lower( + triang=triang, + blnktth=blnktth, + c_shldith=c_shldith, + c_blnkoth=c_blnkoth, + blnkith=blnkith, + blnkoth=blnkoth, + divgap=divgap, + ) + + rs = np.concatenate( + [ + rs_upper_outboard, + rs_lower_inboard, + rs_upper_inboard[::-1], + rs_lower_outboard[::-1], + ] + ) + zs = np.concatenate( + [ + zs_upper_outboard, + zs_lower_inboard, + zs_upper_inboard[::-1], + zs_lower_outboard[::-1], + ] + ) + + return ArbitraryGeometry(rs=rs, zs=zs) + + +def blanket_geometry_lower( + triang: float, + blnktth: float, + c_shldith: float, + c_blnkoth: float, + blnkith: float, + blnkoth: float, + divgap: float, +) -> Tuple[np.ndarray, np.ndarray, np.ndarray, np.ndarray]: + """Calculates radial and vertical distances for the geometry of section of blanket below the midplane + + :param triang: plasma triangularity + :type triang: float + :param blnktth: top blanket vertical thickness + :type blnktth: float + :param c_shldith: inboard shield thickness + :type c_shldith: float + :param c_blnkoth: outboard blanket radial thickness + :type c_blnkoth: float + :param blnkith: inboard blanket radial thickness + :type blnkith: float + :param blnkoth: outboard blanket radial thickness + :type blnkoth: float + :param divgap: divertor structure vertical thickness + :type divgap: float + :return: tuple containing the R coordinates for the outboard, Z coordinates for the outboard, R coordinates for the inboard, Z coordinates for the inboard of the blanket geometry below the midplane + :rtype: Tuple[np.ndarray, np.ndarray, np.ndarray, np.ndarray] + """ + # Lower blanket + rs1, rs2, rs3, rs4, zs1, zs2, zs3, zs4 = dhgap_vertices( + c_shldith, c_blnkoth, blnkith, blnkoth, divgap, -blnktth, triang + ) + # outboard radial and vertical coordinates + rs_lower_outboard = np.concatenate([rs1, rs2[::-1]]) + zs_lower_outboard = np.concatenate([zs1, zs2[::-1]]) + # inboard radial and vertical coordinates + rs_lower_inboard = np.concatenate([rs3, rs4[::-1]]) + zs_lower_inboard = -np.concatenate([zs3, zs4[::-1]]) + + return rs_lower_outboard, zs_lower_outboard, rs_lower_inboard, zs_lower_inboard + + +def blanket_geometry_double_null( + cumulative_lower: dict, + triang: float, + blnktth: float, + c_shldith: float, + c_blnkoth: float, + blnkith: float, + blnkoth: float, +) -> ArbitraryGeometry: + """Calculates radial and vertical distances for the geometry of blanket in a double null configuration + In a double null configuration, the geometry of the lower blanket is reflected across the midplane to create the section of blanket above the midplane + + :param cumulative_lower: cumulative vertical thicknesses of components below the midplane + :type cumulative_lower: dict + :param triang: plasma triangularity + :type triang: float + :param blnktth: top blanket vertical thickness + :type blnktth: float + :param c_shldith: inboard shield thickness + :type c_shldith: float + :param c_blnkoth: outboard blanket radial thickness + :type c_blnkoth: float + :param blnkith: inboard blanket radial thickness + :type blnkith: float + :param blnkoth: outboard blanket radial thickness + :type blnkoth: float + :return: dataclass returning radial and vertical coordinates + :rtype: ArbitraryGeometry + """ + # Lower blanket + divgap = cumulative_lower["divfix"] + + ( + rs_lower_outboard, + zs_lower_outboard, + rs_lower_inboard, + zs_lower_inboard, + ) = blanket_geometry_lower( + triang=triang, + blnktth=blnktth, + c_shldith=c_shldith, + c_blnkoth=c_blnkoth, + blnkith=blnkith, + blnkoth=blnkoth, + divgap=divgap, + ) + + # Upper blanket + divgap = -1 * divgap + + ( + rs_upper_outboard, + zs_upper_outboard, + rs_upper_inboard, + zs_upper_inboard, + ) = blanket_geometry_lower( + triang=triang, + blnktth=blnktth, + c_shldith=c_shldith, + c_blnkoth=c_blnkoth, + blnkith=blnkith, + blnkoth=blnkoth, + divgap=divgap, + ) + rs_1 = np.concatenate([rs_upper_outboard, rs_lower_outboard[::-1]]) + rs_2 = np.concatenate([rs_upper_inboard, rs_lower_inboard[::-1]]) + zs_1 = np.concatenate([zs_upper_outboard, zs_lower_outboard[::-1]]) + zs_2 = np.concatenate([zs_upper_inboard, zs_lower_inboard[::-1]]) + rs = [rs_1, rs_2] + zs = [zs_1, zs_2] + + return ArbitraryGeometry(rs=rs, zs=zs) diff --git a/process/geometry/cryostat_geometry.py b/process/geometry/cryostat_geometry.py new file mode 100644 index 0000000000..35535998ed --- /dev/null +++ b/process/geometry/cryostat_geometry.py @@ -0,0 +1,41 @@ +""" +Calculate cryostat geometries +""" +from typing import List + +from process.geometry.geometry_parameterisations import RectangleGeometry + + +def cryostat_geometry( + rdewex: float, ddwex: float, zdewex: float +) -> List[RectangleGeometry]: + """Calculates rectangular geometries of the cryostat + + :param rdewex: cryostat internal radius + :type rdewex: float + :param ddwex: external cryostat thickness + :type ddwex: float + :param zdewex: cryostat internal half-height + :type zdewex: float + :return: list of RectangleGeometry - dataclass returning rectangular geometry parameters + :rtype: List[RectangleGeometry] + """ + # The cryostat is represented by 4 rectangular sections as follows: + + # rectangle representing vertical part of cryostat above the midplane + rect1 = RectangleGeometry( + anchor_x=rdewex, anchor_z=0, width=ddwex, height=(zdewex + ddwex) + ) + + # rectangle representing vertical part of cryostat below the midplane + rect2 = RectangleGeometry( + anchor_x=rdewex, anchor_z=0, width=ddwex, height=-(zdewex + ddwex) + ) + + # rectangle representing horizontal part of cryostat above the midplane + rect3 = RectangleGeometry(anchor_x=0, anchor_z=zdewex, width=rdewex, height=ddwex) + + # rectangle representing horizontal part of cryostat below the midplane + rect4 = RectangleGeometry(anchor_x=0, anchor_z=-zdewex, width=rdewex, height=-ddwex) + return_rects = [rect1, rect2, rect3, rect4] + return return_rects diff --git a/process/geometry/firstwall_geometry.py b/process/geometry/firstwall_geometry.py new file mode 100644 index 0000000000..b47c54b1c1 --- /dev/null +++ b/process/geometry/firstwall_geometry.py @@ -0,0 +1,227 @@ +""" +Calculate radial and vertical coordinates for the geometry of the first wall +""" +from typing import Tuple +import numpy as np +from process.geometry.utils import dh_vertices, dhgap_vertices +from process.geometry.geometry_parameterisations import ArbitraryGeometry + + +def first_wall_geometry_single_null( + radx_outer: float, + rminx_outer: float, + radx_inner: float, + rminx_inner: float, + cumulative_upper: dict, + triang: float, + cumulative_lower: dict, + blnktth: float, + c_blnkith: float, + c_fwoth: float, + fwith: float, + fwoth: float, + tfwvt: float, +) -> ArbitraryGeometry: + """Calculates radial and vertical distances for the geometry of first wall in a single null configuration + + :param radx_outer: outboard radius of outer surface of first wall + :type radx_outer: float + :param rminx_outer: inboard radius of outer surface of first wall + :type rminx_outer: float + :param radx_inner: outboard radius of inner surface of first wall + :type radx_inner: float + :param rminx_inner: inboard radius of inner surface of first wall + :type rminx_inner: float + :param cumulative_upper: cumulative vertical thicknesses of components above the midplane + :type cumulative_upper: dict + :param triang: plasma triangularity + :type triang: float + :param cumulative_lower: cumulative vertical thicknesses of components below the midplane + :type cumulative_lower: dict + :param blnktth: top blanket vertical thickness + :type blnktth: float + :param c_blnkith: inboard blanket vertical thickness + :type c_blnkith: float + :param c_fwoth: outboard first wall vertical thickness + :type c_fwoth: float + :param fwith: inboard first wall radial thickness + :type fwith: float + :param fwoth: outboard first wall radial thickness + :type fwoth: float + :param tfwvt: top first wall vertical thickness + :type tfwvt: float + :return: dataclass returning radial and vertical coordinates + :rtype: ArbitraryGeometry + """ + # Upper first wall: outer surface + kapx = cumulative_upper["fwtth"] / rminx_outer + rs_upper_outboard, zs_upper_outboard = dh_vertices( + radx_outer, rminx_outer, triang, kapx + ) + + # Upper first wall: inner surface + rs_upper_inboard, zs_upper_inboard = dh_vertices( + radx_inner, rminx_inner, triang, kapx + ) + + # Lower first wall + divgap = cumulative_lower["divfix"] + top_point = divgap + blnktth + ( + rs_lower_outboard, + zs_lower_outboard, + rs_lower_inboard, + zs_lower_inboard, + ) = first_wall_geometry_lower( + triang=triang, + c_blnkith=c_blnkith, + c_fwoth=c_fwoth, + fwith=fwith, + fwoth=fwoth, + tfwvt=tfwvt, + top_point=top_point, + ) + + rs = np.concatenate( + [ + rs_upper_outboard, + rs_lower_inboard, + rs_upper_inboard[::-1], + rs_lower_outboard[::-1], + ] + ) + zs = np.concatenate( + [ + zs_upper_outboard, + zs_lower_inboard, + zs_upper_inboard[::-1], + zs_lower_outboard[::-1], + ] + ) + + return ArbitraryGeometry(rs=rs, zs=zs) + + +def first_wall_geometry_lower( + triang: float, + c_blnkith: float, + c_fwoth: float, + fwith: float, + fwoth: float, + tfwvt: float, + top_point: float, +) -> Tuple[np.ndarray, np.ndarray, np.ndarray, np.ndarray]: + """Calculates radial and vertical distances for the geometry of section of first wall below the midplane + + :param triang: plasma triangularity + :type triang: float + :param c_blnkith: inboard blanket vertical thickness + :type c_blnkith: float + :param c_fwoth: outboard first wall vertical thickness + :type c_fwoth: float + :param fwith: inboard first wall radial thickness + :type fwith: float + :param fwoth: outboard first wall radial thickness + :type fwoth: float + :param tfwvt: top first wall vertical thickness + :type tfwvt: float + :param top_point: top point for plotdhgap, equal to + :type top_point: float + :return: tuple containing the R coordinates for the outboard, Z coordinates for the outboard, R coordinates for the inboard, Z coordinates for the inboard of the first wall geometry below the midplane + :rtype: Tuple[np.ndarray, np.ndarray, np.ndarray, np.ndarray] + """ + # Lower first wall + rs1, rs2, rs3, rs4, zs1, zs2, zs3, zs4 = dhgap_vertices( + c_blnkith, + c_fwoth, + fwith, + fwoth, + top_point, + -tfwvt, + triang, + ) + # outboard radial and vertical coordinates + rs_lower_outboard = np.concatenate([rs1, rs2[::-1]]) + zs_lower_outboard = np.concatenate([zs1, zs2[::-1]]) + # inboard radial and vertical coordinates + rs_lower_inboard = np.concatenate([rs3, rs4[::-1]]) + zs_lower_inboard = -np.concatenate([zs3, zs4[::-1]]) + + return rs_lower_outboard, zs_lower_outboard, rs_lower_inboard, zs_lower_inboard + + +def first_wall_geometry_double_null( + cumulative_lower: dict, + triang: float, + blnktth: float, + c_blnkith: float, + c_fwoth: float, + fwith: float, + fwoth: float, + tfwvt: float, +) -> ArbitraryGeometry: + """Calculates radial and vertical distances for the geometry of first wall in a double null configuration + In a double null configuration, the geometry of the lower first wall is reflected across the midplane to create the section of first wall above the midplane + + :param cumulative_lower: cumulative vertical thicknesses of components below the midplane + :type cumulative_lower: dict + :param triang: plasma triangularity + :type triang: float + :param blnktth: top blanket vertical thickness + :type blnktth: float + :param c_blnkith: inboard blanket vertical thickness + :type c_blnkith: float + :param c_fwoth: outboard first wall vertical thickness + :type c_fwoth: float + :param fwith: inboard first wall radial thickness + :type fwith: float + :param fwoth: outboard first wall radial thickness + :type fwoth: float + :param tfwvt: top first wall vertical thickness + :type tfwvt: float + :return: dataclass returning radial and vertical coordinates + :rtype: ArbitraryGeometry + """ + # Lower first wall + divgap = cumulative_lower["divfix"] + top_point = divgap + blnktth + ( + rs_lower_outboard, + zs_lower_outboard, + rs_lower_inboard, + zs_lower_inboard, + ) = first_wall_geometry_lower( + triang=triang, + c_blnkith=c_blnkith, + c_fwoth=c_fwoth, + fwith=fwith, + fwoth=fwoth, + tfwvt=tfwvt, + top_point=top_point, + ) + + # Upper first wall + top_point = -1 * top_point + ( + rs_upper_outboard, + zs_upper_outboard, + rs_upper_inboard, + zs_upper_inboard, + ) = first_wall_geometry_lower( + triang=triang, + c_blnkith=c_blnkith, + c_fwoth=c_fwoth, + fwith=fwith, + fwoth=fwoth, + tfwvt=tfwvt, + top_point=top_point, + ) + + rs_1 = np.concatenate([rs_upper_outboard, rs_lower_outboard[::-1]]) + rs_2 = np.concatenate([rs_upper_inboard, rs_lower_inboard[::-1]]) + zs_1 = np.concatenate([zs_upper_outboard, zs_lower_outboard[::-1]]) + zs_2 = np.concatenate([zs_upper_inboard, zs_lower_inboard[::-1]]) + rs = [rs_1, rs_2] + zs = [zs_1, zs_2] + + return ArbitraryGeometry(rs=rs, zs=zs) diff --git a/process/geometry/geometry_parameterisations.py b/process/geometry/geometry_parameterisations.py new file mode 100644 index 0000000000..2e312a181c --- /dev/null +++ b/process/geometry/geometry_parameterisations.py @@ -0,0 +1,31 @@ +""" +Module to hold geometry parameterisation dataclasses common to multiple reactor components +""" +from dataclasses import dataclass +import numpy as np + + +@dataclass +class RectangleGeometry: + """Holds data for rectangular geometries""" + + anchor_x: float + """rectangle x coordinate anchor point""" + anchor_z: float + """rectangle z coordinate anchor point""" + width: float + """rectangle width""" + height: float + """rectangle height""" + + +@dataclass +class ArbitraryGeometry: + """Holds radial and vertical coordinates for arbitrary reactor component shapes + + Example: a triangular shaped component with vertices [(0,0), (1,0), (0,1)] would be represented by ArbitraryGeometry(rs=[0,1,0], zs=[0,0,1])""" + + rs: np.ndarray + """outboard and inboard radial coordinates""" + zs: np.ndarray + """outboard and inboard vertical coordinates""" diff --git a/process/geometry/pfcoil_geometry.py b/process/geometry/pfcoil_geometry.py new file mode 100644 index 0000000000..7d3e1974ec --- /dev/null +++ b/process/geometry/pfcoil_geometry.py @@ -0,0 +1,51 @@ +""" +Calculate radial and vertical coordinates for the geometry of the pf coils and central coil +""" +from typing import List, Tuple +import numpy as np +from process.geometry.geometry_parameterisations import RectangleGeometry + + +def pfcoil_geometry( + coils_r: List[float], + coils_z: List[float], + coils_dr: List[float], + coils_dz: List[float], + bore: float, + ohcth: float, + ohdz: float, +) -> Tuple[np.ndarray, np.ndarray, RectangleGeometry]: + """Calculates radial and vertical distances for the geometry of the pf coils and central coil + + :param coils_r: list of pf coil radii + :type coils_r: List[float] + :param coils_z: list of pf coil vertical positions + :type coils_z: List[float] + :param coils_dr: list of pf coil radial thicknesses + :type coils_dr: List[float] + :param coils_dz: list of pf coil vertical thicknesses + :type coils_dz: List[float] + :param bore: central solenoid inboard radius + :type bore: float + :param ohcth: central solenoid thickness + :type ohcth: float + :param ohdz: central solenoid vertical thickness + :type ohdz: float + :return: tuple containing radial and vertical coordinates for pf coils, and dataclass returning coordinates representing a rectangular geometry used to plot the central coil + :rtype: Tuple[np.ndarray, np.ndarray, RectangleGeometry] + """ + r_points = [] + z_points = [] + for i in range(len(coils_r)): + r_1 = float(coils_r[i]) - 0.5 * float(coils_dr[i]) + z_1 = float(coils_z[i]) - 0.5 * float(coils_dz[i]) + z_2 = float(coils_z[i]) + 0.5 * float(coils_dz[i]) + r_2 = float(coils_r[i]) + 0.5 * float(coils_dr[i]) + r_points.append([r_1, r_1, r_2, r_2, r_1]) + z_points.append([z_1, z_2, z_2, z_1, z_1]) + + central_coil = RectangleGeometry( + anchor_x=bore, anchor_z=(-ohdz / 2), width=ohcth, height=ohdz + ) + + return r_points, z_points, central_coil diff --git a/process/geometry/plasma_geometry.py b/process/geometry/plasma_geometry.py new file mode 100644 index 0000000000..070db744af --- /dev/null +++ b/process/geometry/plasma_geometry.py @@ -0,0 +1,77 @@ +""" +Calculate plasma elongation and radial and vertical coordinates for the geometry of the plasma +""" +from dataclasses import dataclass +import math +import numpy as np + + +@dataclass +class PlasmaGeometry: + """Holds radial and vertical coordinates for the geometry of a plasma""" + + rs: np.ndarray + """radial plasma coordinates""" + zs: np.ndarray + """vertical plasma coordinates""" + kappa: float + """plasma separatrix elongation""" + + +def plasma_geometry( + r_0: float, a: float, triang_95: float, kappa_95: float, i_single_null: bool +) -> PlasmaGeometry: + """Calculates radial and vertical distances and plasma elongation for the geometry of the plasma + + :param r_0: plasma major radius + :type r_0: float + :param a: plasma minor radius + :type a: float + :param triang_95: plasma triangularity at 95% surface + :type triang_95: float + :param kappa_95: plasma elongation at 95% surface + :type kappa_95: float + :param i_single_null: switch for single null / double null plasma + :type i_single_null: bool + :return: dataclass returning plasma elongation and radial and vertical coordinates of the plasma + :rtype: PlasmaGeometry + """ + delta = 1.5 * triang_95 + kappa = (1.1 * kappa_95) + 0.04 + + x1 = (2.0 * r_0 * (1.0 + delta) - a * (delta**2 + kappa**2 - 1.0)) / ( + 2.0 * (1.0 + delta) + ) + x2 = (2.0 * r_0 * (delta - 1.0) - a * (delta**2 + kappa**2 - 1.0)) / ( + 2.0 * (delta - 1.0) + ) + r1 = 0.5 * math.sqrt( + (a**2 * ((delta + 1.0) ** 2 + kappa**2) ** 2) / ((delta + 1.0) ** 2) + ) + r2 = 0.5 * math.sqrt( + (a**2 * ((delta - 1.0) ** 2 + kappa**2) ** 2) / ((delta - 1.0) ** 2) + ) + theta1 = np.arcsin((kappa * a) / r1) + theta2 = np.arcsin((kappa * a) / r2) + inang = 1.0 / r1 + outang = 1.5 / r2 + if i_single_null == 0: + angs1 = np.linspace( + -(inang + theta1) + np.pi, (inang + theta1) + np.pi, 256, endpoint=True + ) + angs2 = np.linspace(-(outang + theta2), (outang + theta2), 256, endpoint=True) + else: + angs1 = np.linspace( + -theta1 + np.pi, (inang + theta1) + np.pi, 256, endpoint=True + ) + angs2 = np.linspace(-(outang + theta2), theta2, 256, endpoint=True) + + xs1 = -(r1 * np.cos(angs1) - x1) + ys1 = r1 * np.sin(angs1) + xs2 = -(r2 * np.cos(angs2) - x2) + ys2 = r2 * np.sin(angs2) + + rs = [xs1, xs2] + zs = [ys1, ys2] + + return PlasmaGeometry(rs=rs, zs=zs, kappa=kappa) diff --git a/process/geometry/shield_geometry.py b/process/geometry/shield_geometry.py new file mode 100644 index 0000000000..f31b93e43d --- /dev/null +++ b/process/geometry/shield_geometry.py @@ -0,0 +1,179 @@ +""" +Calculate radial and vertical coordinates for the geometry of the shield +""" +from typing import Tuple +import numpy as np +from process.geometry.utils import dh_vertices +from process.geometry.geometry_parameterisations import ArbitraryGeometry + + +def shield_geometry_single_null( + cumulative_upper: dict, + radx_far: float, + rminx_far: float, + radx_near: float, + rminx_near: float, + triang: float, + cumulative_lower: dict, +) -> ArbitraryGeometry: + """Calculates radial and vertical distances for the geometry of shield in a single null configuration + + :param cumulative_upper: cumulative vertical thicknesses of components above the midplane + :type cumulative_upper: dict + :param radx_far: outboard radius of outer surface of shield + :type radx_far: float + :param rminx_far: inboard radius of outer surface of shield + :type rminx_far: float + :param radx_near: outboard radius of inner surface of shield + :type radx_near: float + :param rminx_near: inboard radius of inner surface of shield + :type rminx_near: float + :param triang: plasma triangularity + :type triang: float + :param cumulative_lower: cumulative vertical thicknesses of components below the midplane + :type cumulative_lower: dict + :return: dataclass returning radial and vertical coordinates + :rtype: ArbitraryGeometry + """ + # Upper shield + # Side furthest from plasma + kapx = cumulative_upper["shldtth"] / rminx_far + rs_upper_outboard, zs_upper_outboard = dh_vertices( + radx_far, rminx_far, triang, kapx + ) + + # Side nearest to plasma + kapx = (cumulative_upper["vvblgap"]) / rminx_near + rs_upper_inboard, zs_upper_inboard = dh_vertices( + radx_near, rminx_near, triang, kapx + ) + + # Lower shield + ( + rs_lower_outboard, + zs_lower_outboard, + rs_lower_inboard, + zs_lower_inboard, + ) = shield_geometry_lower( + cumulative_lower=cumulative_lower, + radx_far=radx_far, + radx_near=radx_near, + rminx_far=rminx_far, + rminx_near=rminx_near, + triang=triang, + ) + + rs = np.concatenate( + [ + rs_lower_inboard, + rs_lower_outboard[::-1], + rs_upper_outboard, + rs_upper_inboard[::-1], + ] + ) + zs = np.concatenate( + [ + zs_lower_inboard, + zs_lower_outboard[::-1], + zs_upper_outboard, + zs_upper_inboard[::-1], + ] + ) + return ArbitraryGeometry( + rs=rs, + zs=zs, + ) + + +def shield_geometry_lower( + cumulative_lower: dict, + radx_far: float, + rminx_far: float, + radx_near: float, + rminx_near: float, + triang: float, +) -> Tuple[np.ndarray, np.ndarray, np.ndarray, np.ndarray]: + """Calculates radial and vertical distances for the geometry of section of shield below the midplane + + :param cumulative_lower: cumulative vertical thicknesses of components below the midplane + :type cumulative_lower: dict + :param radx_far: outboard radius of outer surface of shield + :type radx_far: float + :param rminx_far: inboard radius of outer surface of shield + :type rminx_far: float + :param radx_near: outboard radius of inner surface of shield + :type radx_near: float + :param rminx_near: inboard radius of inner surface of shield + :type rminx_near: float + :param triang: plasma triangularity + :type triang: float + :return: tuple containing the R coordinates for the outboard, Z coordinates for the outboard, R coordinates for the inboard, Z coordinates for the inboard of the shield geometry below the midplane + :rtype: Tuple[np.ndarray, np.ndarray, np.ndarray, np.ndarray] + """ + # Side furthest from plasma + kapx = cumulative_lower["shldlth"] / rminx_far + rs_lower_outboard, zs_lower_outboard = dh_vertices( + radx_far, rminx_far, triang, kapx + ) + + # Side nearest to plasma + kapx = (cumulative_lower["divfix"]) / rminx_near + rs_lower_inboard, zs_lower_inboard = dh_vertices( + radx_near, rminx_near, triang, kapx + ) + + return rs_lower_outboard, zs_lower_outboard, rs_lower_inboard, zs_lower_inboard + + +def shield_geometry_double_null( + cumulative_lower: dict, + radx_far: float, + rminx_far: float, + radx_near: float, + rminx_near: float, + triang: float, +) -> ArbitraryGeometry: + """Calculates radial and vertical distances for the geometry of shield in a double null configuration + In a double null configuration, the geometry of the lower shield is reflected across the midplane to create the section of shield above the midplane + + :param cumulative_lower: cumulative vertical thicknesses of components below the midplane + :type cumulative_lower: dict + :param radx_far: outboard radius of outer surface of shield + :type radx_far: float + :param rminx_far: inboard radius of outer surface of shield + :type rminx_far: float + :param radx_near: outboard radius of inner surface of shield + :type radx_near: float + :param rminx_near: inboard radius of inner surface of shield + :type rminx_near: float + :param triang: plasma triangularity + :type triang: float + :return: dataclass returning radial and vertical coordinates + :rtype: ArbitraryGeometry + """ + # Lower shield + ( + rs_lower_outboard, + zs_lower_outboard, + rs_lower_inboard, + zs_lower_inboard, + ) = shield_geometry_lower( + cumulative_lower=cumulative_lower, + radx_far=radx_far, + radx_near=radx_near, + rminx_far=rminx_far, + rminx_near=rminx_near, + triang=triang, + ) + + rs_lower = np.concatenate([rs_lower_outboard, rs_lower_inboard[::-1]]) + zs_lower = np.concatenate([zs_lower_outboard, zs_lower_inboard[::-1]]) + + # Upper shield + rs_upper = rs_lower + zs_upper = -1 * zs_lower + + rs = np.concatenate([rs_lower[::-1], rs_upper]) + zs = np.concatenate([zs_lower[::-1], zs_upper]) + + return ArbitraryGeometry(rs=rs, zs=zs) diff --git a/process/geometry/tfcoil_geometry.py b/process/geometry/tfcoil_geometry.py new file mode 100644 index 0000000000..23a9f0cb8d --- /dev/null +++ b/process/geometry/tfcoil_geometry.py @@ -0,0 +1,197 @@ +""" +Calculate radial and vertical coordinates for the geometry of the tf coils +""" +from typing import List, Tuple +from process.geometry.geometry_parameterisations import RectangleGeometry +from process.geometry.utils import ellips_fill_vertices + + +def tfcoil_geometry_rectangular_shape( + x1: float, + x2: float, + x4: float, + x5: float, + y1: float, + y2: float, + y4: float, + y5: float, + tfcth: float, +) -> List[RectangleGeometry]: + """Calculates rectangular geometries for tf coils in a picture frame/rectangular shape parametrization + + :param x1: radial location of arc point 1 + :type x1: float + :param x2: radial location of arc point 2 + :type x2: float + :param x4: radial location of arc point 4 + :type x4: float + :param x5: radial location of arc point 5 + :type x5: float + :param y1: vertical location of arc point 1 + :type y1: float + :param y2: vertical location of arc point 2 + :type y2: float + :param y4: vertical location of arc point 4 + :type y4: float + :param y5: vertical location of arc point 5 + :type y5: float + :param tfcth: inboard tf coil thickness + :type tfcth: float + :return: list of RectangleGeometry - dataclass returning rectangular geometry parameters + :rtype: List[RectangleGeometry] + """ + return_rects = [] + # In this geometry, the tf coil is represented by 4 rectangular sections as follows: + + # rectangle representing the inboard part of the tf coil + return_rects.append( + RectangleGeometry( + anchor_x=x5 - tfcth, + anchor_z=y5 - tfcth, + width=tfcth, + height=(y1 - y5 + 2.0 * tfcth), + ) + ) + # rectangle representing the outboard part of the tf coil + return_rects.append( + RectangleGeometry( + anchor_x=x4, + anchor_z=y4 - tfcth, + width=tfcth, + height=(y2 - y4 + 2.0 * tfcth), + ) + ) + # rectangle representing the lower horizontal part of the tf coil + return_rects.append( + RectangleGeometry(anchor_x=x5, anchor_z=y5 - tfcth, width=x4 - x5, height=tfcth) + ) + # rectangle representing the upper horizontal part of the tf coil + return_rects.append( + RectangleGeometry(anchor_x=x1, anchor_z=y1, width=(x2 - x1), height=tfcth) + ) + + return return_rects + + +def tfcoil_geometry_d_shape( + x1: float, + x2: float, + x3: float, + x4: float, + x5: float, + y1: float, + y2: float, + y4: float, + y5: float, + tfcth: float, + rtangle: float, + rtangle2: float, +) -> Tuple[List[RectangleGeometry], List[List[Tuple[float, float]]]]: + """Calculates radial and vertical distances for the geometry of the tf coils in a D-shape parametrization + + :param x1: radial location of arc point 1 + :type x1: float + :param x2: radial location of arc point 2 + :type x2: float + :param x3: radial location of arc point 3 + :type x3: float + :param x4: radial location of arc point 4 + :type x4: float + :param x5: radial location of arc point 5 + :type x5: float + :param y1: vertical location of arc point 1 + :type y1: float + :param y2: vertical location of arc point 2 + :type y2: float + :param y4: vertical location of arc point 4 + :type y4: float + :param y5: vertical location of arc point 5 + :type y5: float + :param tfcth: inboard tf coil thickness + :type tfcth: float + :param rtangle: angle used in tf coil parametrization + :type rtangle: float + :param rtangle2: angle used in tf coil parametrization + :type rtangle2: float + :return: radial and vertical coordinates for tf coils + :rtype: Tuple[List[RectangleGeometry], List[List[Tuple[float, float]]]] + """ + return_rects = [] + return_verts = [] + # Inboard upper arc + x0 = x2 + y0 = y1 + a1 = x2 - x1 + b1 = y2 - y1 + a2 = a1 + tfcth + b2 = b1 + tfcth + verts1 = ellips_fill_vertices( + a1=a1, + a2=a2, + b1=b1, + b2=b2, + x0=x0, + y0=y0, + ang1=rtangle, + ang2=rtangle2, + ) + # Outboard upper arc + x0 = x2 + y0 = 0 + a1 = x3 - x2 + b1 = y2 + a2 = a1 + tfcth + b2 = b1 + tfcth + verts2 = ellips_fill_vertices( + a1=a1, + a2=a2, + b1=b1, + b2=b2, + x0=x0, + y0=y0, + ang1=0, + ang2=rtangle, + ) + # Inboard lower arc + x0 = x4 + y0 = y5 + a1 = x4 - x5 + b1 = y5 - y4 + a2 = a1 + tfcth + b2 = b1 + tfcth + verts3 = ellips_fill_vertices( + a1=a1, + a2=a2, + b1=b1, + b2=b2, + x0=x0, + y0=y0, + ang1=-rtangle, + ang2=-1 * rtangle2, + ) + # Outboard lower arc + x0 = x4 + y0 = 0 + a1 = x3 - x2 + b1 = -y4 + a2 = a1 + tfcth + b2 = b1 + tfcth + verts4 = ellips_fill_vertices( + a1=a1, + a2=a2, + b1=b1, + b2=b2, + x0=x0, + y0=y0, + ang1=0, + ang2=-rtangle, + ) + # Vertical leg + return_rects.append( + RectangleGeometry( + anchor_x=x5 - tfcth, anchor_z=y5, width=tfcth, height=(y1 - y5) + ) + ) + return_verts = [verts1, verts2, verts3, verts4] + + return return_rects, return_verts diff --git a/process/geometry/utils.py b/process/geometry/utils.py new file mode 100644 index 0000000000..506a948a0f --- /dev/null +++ b/process/geometry/utils.py @@ -0,0 +1,132 @@ +""" +Module to hold plotting functions, used in plot_proc.py, which are common to multiple reactor components +""" +from typing import List, Tuple +import numpy as np + + +def dh_vertices( + r0: float, a: float, triang: float, kap: float +) -> Tuple[np.ndarray, np.ndarray]: + """Returns the radial and vertical coordinates which, when plotted, plots half a thin D-section, centred on z = 0 + + :param r0: major radius of centre + :type r0: float + :param a: horizontal radius + :type a: float + :param triang: plasma triangularity + :type triang: float + :param kap: plasma elongation + :type kap: float + :return: tuple containing radial and vertical coordinates which, when plotted, plots a half thin D-section with a gap + :rtype: Tuple[np.ndarray, np.ndarray] + """ + angs = np.linspace(0, np.pi, 50, endpoint=True) + rs = r0 + a * np.cos(angs + triang * np.sin(1.0 * angs)) + zs = kap * a * np.sin(angs) + return rs, zs + + +def dhgap_vertices( + inpt: float, + outpt: float, + inthk: float, + outthk: float, + toppt: float, + topthk: float, + triang: float, +) -> Tuple[ + np.ndarray, + np.ndarray, + np.ndarray, + np.ndarray, + np.ndarray, + np.ndarray, + np.ndarray, + np.ndarray, +]: + """Returns the radial and vertical coordinates which, when plotted, plots a half thick D-section with a gap + + :param inpt: inner point + :type inpt: float + :param outpt: outer point + :type outpt: float + :param inthk: inner thickness + :type inthk: float + :param outthk: outer thickness + :type outthk: float + :param toppt: top point + :type toppt: float + :param topthk: top thickness + :type topthk: float + :param triang: plasma triangularity + :type triang: float + :return: tuple containing radial and vertical coordinates which, when plotted, plots a half thick D-section with a gap + :rtype: Tuple[np.ndarray, np.ndarray, np.ndarray, np.ndarray, np.ndarray, np.ndarray, np.ndarray, np.ndarray,] + """ + arc = np.pi / 4.0 + r01 = (inpt + outpt) / 2.0 + r02 = (inpt + inthk + outpt - outthk) / 2.0 + a1 = r01 - inpt + a2 = r02 - inpt - inthk + kap1 = toppt / a1 + kap2 = (toppt - topthk) / a2 + angs = np.linspace(0.0, (np.pi / 2.0) - arc / 2.0, 50, endpoint=True) + rs1 = r01 + a1 * np.cos(angs + triang * np.sin(angs)) + zs1 = kap1 * a1 * np.sin(angs) + rs2 = r02 + a2 * np.cos(angs + triang * np.sin(angs)) + zs2 = kap2 * a2 * np.sin(angs) + angs = np.linspace(np.pi, np.pi + ((np.pi / 2.0) - arc), 50, endpoint=True) + rs3 = r01 + a1 * np.cos(angs + triang * np.sin(angs)) + zs3 = kap1 * a1 * np.sin(angs) + rs4 = r02 + a2 * np.cos(angs + triang * np.sin(angs)) + zs4 = kap2 * a2 * np.sin(angs) + + return rs1, rs2, rs3, rs4, zs1, zs2, zs3, zs4 + + +def ellips_fill_vertices( + a1: float = 0, + a2: float = 0, + b1: float = 0, + b2: float = 0, + x0: float = 0, + y0: float = 0, + ang1: float = 0, + ang2: float = np.pi / 2, +) -> List[Tuple[float, float]]: + """Returns the vertices of a shape which, when filled, fills the space between two concentric ellipse sectors + + :param a1: horizontal radius to be filled, defaults to 0 + :type a1: float, optional + :param a2: horizontal radius to be filled, defaults to 0 + :type a2: float, optional + :param b1: vertical radius to be filled, defaults to 0 + :type b1: float, optional + :param b2: vertical radius to be filled, defaults to 0 + :type b2: float, optional + :param x0: x coordinate of centre of ellipses, defaults to 0 + :type x0: float, optional + :param y0: y coordinate of centre of ellipses, defaults to 0 + :type y0: float, optional + :param ang1: polar angle at start, defaults to 0 + :type ang1: float, optional + :param ang2: polar angle at end, defaults to np.pi/2 + :type ang2: float, optional + :return: list containing (R,Z) coordinates which, when plotted, fill space between ellipses + :rtype: List[Tuple[float, float]] + """ + angs = np.linspace(ang1, ang2, endpoint=True) + r1 = ((np.cos(angs) / a1) ** 2 + (np.sin(angs) / b1) ** 2) ** (-0.5) + xs1 = r1 * np.cos(angs) + x0 + ys1 = r1 * np.sin(angs) + y0 + angs = np.linspace(ang2, ang1, endpoint=True) + r2 = ((np.cos(angs) / a2) ** 2 + (np.sin(angs) / b2) ** 2) ** (-0.5) + xs2 = r2 * np.cos(angs) + x0 + ys2 = r2 * np.sin(angs) + y0 + verts = list(zip(xs1, ys1)) + verts.extend(list(zip(xs2, ys2))) + endpoint = verts[-1:] + verts.extend(endpoint) + + return verts diff --git a/process/geometry/vacuum_vessel_geometry.py b/process/geometry/vacuum_vessel_geometry.py new file mode 100644 index 0000000000..693699f0f3 --- /dev/null +++ b/process/geometry/vacuum_vessel_geometry.py @@ -0,0 +1,192 @@ +""" +Calculate radial and vertical coordinates for the geometry of the vacuum vessel +""" +from typing import Tuple +import numpy as np +from process.geometry.utils import dh_vertices +from process.geometry.geometry_parameterisations import ArbitraryGeometry + + +def vacuum_vessel_geometry_single_null( + cumulative_upper: dict, + upper: dict, + triang: float, + radx_outer: float, + rminx_outer: float, + radx_inner: float, + rminx_inner: float, + cumulative_lower: dict, + lower: dict, +) -> ArbitraryGeometry: + """Calculates radial and vertical distances for the geometry of vacuum vessel in a single null configuration + + :param cumulative_upper: cumulative vertical thicknesses of components above the midplane + :type cumulative_upper: dict + :param upper: vertical thicknesses of components above the midplane + :type upper: dict + :param triang: plasma triangularity + :type triang: float + :param radx_outer: outboard radius of outer surface of vacuum vessel + :type radx_outer: float + :param rminx_outer: inboard radius of outer surface of vacuum vessel + :type rminx_outer: float + :param radx_inner: outboard radius of inner surface of vacuum vessel + :type radx_inner: float + :param rminx_inner: inboard radius of inner surface of vacuum vessel + :type rminx_inner: float + :param cumulative_lower: cumulative vertical thicknesses of components below the midplane + :type cumulative_lower: dict + :param lower: vertical thicknesses of components below the midplane + :type lower: dict + :return: dataclass returning radial and vertical coordinates + :rtype: ArbitraryGeometry + """ + # Upper vacuum vessel + kapx = cumulative_upper["d_vv_top"] / rminx_outer + rs_upper_outboard, zs_upper_outboard = dh_vertices( + radx_outer, rminx_outer, triang, kapx + ) + + kapx = ( + float(cumulative_upper["d_vv_top"]) - float(upper["d_vv_top"]) + ) / rminx_inner + rs_upper_inboard, zs_upper_inboard = dh_vertices( + radx_inner, rminx_inner, triang, kapx + ) + + # Lower vacuum vessel + ( + rs_lower_outboard, + zs_lower_outboard, + rs_lower_inboard, + zs_lower_inboard, + ) = vacuum_vessel_geometry_lower( + cumulative_lower=cumulative_lower, + lower=lower, + radx_inner=radx_inner, + radx_outer=radx_outer, + rminx_inner=rminx_inner, + rminx_outer=rminx_outer, + triang=triang, + ) + + rs = np.concatenate( + [ + rs_lower_inboard, + rs_lower_outboard[::-1], + rs_upper_outboard, + rs_upper_inboard[::-1], + ] + ) + zs = np.concatenate( + [ + zs_lower_inboard, + zs_lower_outboard[::-1], + zs_upper_outboard, + zs_upper_inboard[::-1], + ] + ) + return ArbitraryGeometry( + rs=rs, + zs=zs, + ) + + +def vacuum_vessel_geometry_lower( + cumulative_lower: dict, + lower: dict, + triang: float, + radx_outer: float, + rminx_outer: float, + radx_inner: float, + rminx_inner: float, +) -> Tuple[np.ndarray, np.ndarray, np.ndarray, np.ndarray]: + """Calculates radial and vertical distances for the geometry of section of vacuum vessel below the midplane + + :param cumulative_lower: cumulative vertical thicknesses of components below the midplane + :type cumulative_lower: dict + :param lower: vertical thicknesses of components below the midplane + :type lower: dict + :param triang: plasma triangularity + :type triang: float + :param radx_outer: outboard radius of outer surface of vacuum vessel + :type radx_outer: float + :param rminx_outer: inboard radius of outer surface of vacuum vessel + :type rminx_outer: float + :param radx_inner: outboard radius of inner surface of vacuum vessel + :type radx_inner: float + :param rminx_inner: inboard radius of inner surface of vacuum vessel + :type rminx_inner: float + :return: tuple containing the R coordinates for the outboard, Z coordinates for the outboard, R coordinates for the inboard, Z coordinates for the inboard of the vacuum vessel geometry below the midplane + :rtype: Tuple[np.ndarray, np.ndarray, np.ndarray, np.ndarray] + """ + kapx = cumulative_lower["d_vv_bot"] / rminx_outer + rs_lower_outboard, zs_lower_outboard = dh_vertices( + radx_outer, rminx_outer, triang, kapx + ) + + kapx = ( + float(cumulative_lower["d_vv_bot"]) + float(lower["d_vv_bot"]) + ) / rminx_inner + rs_lower_inboard, zs_lower_inboard = dh_vertices( + radx_inner, rminx_inner, triang, kapx + ) + + return rs_lower_outboard, zs_lower_outboard, rs_lower_inboard, zs_lower_inboard + + +def vacuum_vessel_geometry_double_null( + cumulative_lower: dict, + lower: dict, + triang: float, + radx_outer: float, + rminx_outer: float, + radx_inner: float, + rminx_inner: float, +) -> ArbitraryGeometry: + """Calculates radial and vertical distances for the geometry of vacuum vessel in a double null configuration + In a double null configuration, the geometry of the lower vacuum vessel is reflected across the midplane to create the section of vacuum vessel above the midplane + + :param cumulative_lower: cumulative vertical thicknesses of components below the midplane + :type cumulative_lower: dict + :param lower: vertical thicknesses of components below the midplane + :type lower: dict + :param triang: plasma triangularity + :type triang: float + :param radx_outer: outboard radius of outer surface of vacuum vessel + :type radx_outer: float + :param rminx_outer: inboard radius of outer surface of vacuum vessel + :type rminx_outer: float + :param radx_inner: outboard radius of inner surface of vacuum vessel + :type radx_inner: float + :param rminx_inner: inboard radius of inner surface of vacuum vessel + :type rminx_inner: float + :return: dataclass returning radial and vertical coordinates + :rtype: ArbitraryGeometry + """ + # Lower vacuum vessel + ( + rs_lower_outboard, + zs_lower_outboard, + rs_lower_inboard, + zs_lower_inboard, + ) = vacuum_vessel_geometry_lower( + cumulative_lower=cumulative_lower, + lower=lower, + radx_inner=radx_inner, + radx_outer=radx_outer, + rminx_inner=rminx_inner, + rminx_outer=rminx_outer, + triang=triang, + ) + rs_lower = np.concatenate([rs_lower_outboard, rs_lower_inboard[::-1]]) + zs_lower = np.concatenate([zs_lower_outboard, zs_lower_inboard[::-1]]) + + # Upper vacuum vessel + rs_upper = rs_lower + zs_upper = -1 * zs_lower + + rs = np.concatenate([rs_lower[::-1], rs_upper]) + zs = np.concatenate([zs_lower[::-1], zs_upper]) + + return ArbitraryGeometry(rs=rs, zs=zs) diff --git a/process/io/plot_proc.py b/process/io/plot_proc.py index 06d07af343..e2daf60bfb 100755 --- a/process/io/plot_proc.py +++ b/process/io/plot_proc.py @@ -20,12 +20,35 @@ import matplotlib import matplotlib.pyplot as plt import matplotlib.backends.backend_pdf as bpdf -import math from matplotlib.path import Path import matplotlib.patches as patches import numpy as np import process.io.mfile as mf + +from process.geometry.shield_geometry import ( + shield_geometry_single_null, + shield_geometry_double_null, +) +from process.geometry.plasma_geometry import plasma_geometry +from process.geometry.vacuum_vessel_geometry import ( + vacuum_vessel_geometry_single_null, + vacuum_vessel_geometry_double_null, +) +from process.geometry.blanket_geometry import ( + blanket_geometry_single_null, + blanket_geometry_double_null, +) +from process.geometry.cryostat_geometry import cryostat_geometry +from process.geometry.tfcoil_geometry import ( + tfcoil_geometry_rectangular_shape, + tfcoil_geometry_d_shape, +) +from process.geometry.pfcoil_geometry import pfcoil_geometry +from process.geometry.firstwall_geometry import ( + first_wall_geometry_single_null, + first_wall_geometry_double_null, +) from process.impurity_radiation import read_impurity_file from process.io.python_fortran_dicts import get_dicts @@ -52,7 +75,7 @@ winding = "blue" nbshield_colour = "gray" -thin = 0 +thin = 0.0 RADIAL_BUILD = [ "bore", @@ -116,82 +139,7 @@ ] rtangle = np.pi / 2 - - -def plotdh(axis, r0, a, delta, kap): - """Plots half a thin D-section, centred on z = 0. - - Arguments: - axis --> axis object to plot to - r0 --> major radius of centre - a --> horizontal radius - delta --> triangularity - kap --> elongation - - Returns: - rs --> radial coordinates of D-section - zs --> vertical coordinates of D-section - """ - angs = np.linspace(0, np.pi, 50, endpoint=True) - rs = r0 + a * np.cos(angs + delta * np.sin(1.0 * angs)) - zs = kap * a * np.sin(angs) - axis.plot(rs, zs, color="black", lw=thin) - return rs, zs - - -def plotdhgap(axis, inpt, outpt, inthk, outthk, toppt, topthk, delta, col): - """Plots half a thick D-section with a gap. - - Arguments: - axis --> axis object to plot to - inpt --> inner points - outpt --> outer points - inthk --> inner thickness - outthk --> outer thickness - toppt --> top points - topthk --> top thickness - delta --> triangularity - col --> color for fill - - """ - arc = np.pi / 4.0 - r01 = (inpt + outpt) / 2.0 - r02 = (inpt + inthk + outpt - outthk) / 2.0 - a1 = r01 - inpt - a2 = r02 - inpt - inthk - kap1 = toppt / a1 - kap2 = (toppt - topthk) / a2 - # angs = ((np.pi/2.) - arc/2.) * findgen(50)/49. - angs = np.linspace(0.0, (np.pi / 2.0) - arc / 2.0, 50, endpoint=True) - rs1 = r01 + a1 * np.cos(angs + delta * np.sin(angs)) - zs1 = kap1 * a1 * np.sin(angs) - rs2 = r02 + a2 * np.cos(angs + delta * np.sin(angs)) - zs2 = kap2 * a2 * np.sin(angs) - # angs = !pi + ((!pi/2.) - arc) * findgen(50)/49. - angs = np.linspace(np.pi, np.pi + ((np.pi / 2.0) - arc), 50, endpoint=True) - rs3 = r01 + a1 * np.cos(angs + delta * np.sin(angs)) - zs3 = kap1 * a1 * np.sin(angs) - rs4 = r02 + a2 * np.cos(angs + delta * np.sin(angs)) - zs4 = kap2 * a2 * np.sin(angs) - - axis.plot( - np.concatenate([rs1, rs2[::-1]]), - np.concatenate([zs1, zs2[::-1]]), - color="black", - lw=thin, - ) - axis.plot( - np.concatenate([rs3, rs4[::-1]]), - -np.concatenate([zs3, zs4[::-1]]), - color="black", - lw=thin, - ) - axis.fill( - np.concatenate([rs1, rs2[::-1]]), np.concatenate([zs1, zs2[::-1]]), color=col - ) - axis.fill( - np.concatenate([rs3, rs4[::-1]]), -np.concatenate([zs3, zs4[::-1]]), color=col - ) +rtangle2 = 2 * rtangle def plot_plasma(axis, mfile_data, scan): @@ -204,58 +152,34 @@ def plot_plasma(axis, mfile_data, scan): """ - r0 = mfile_data.data["rmajor"].get_scan(scan) + r_0 = mfile_data.data["rmajor"].get_scan(scan) a = mfile_data.data["rminor"].get_scan(scan) - delta = 1.5 * mfile_data.data["triang95"].get_scan(scan) - kappa = (1.1 * mfile_data.data["kappa95"].get_scan(scan)) + 0.04 + triang_95 = mfile_data.data["triang95"].get_scan(scan) + kappa_95 = mfile_data.data["kappa95"].get_scan(scan) i_single_null = mfile_data.data["i_single_null"].get_scan(scan) - x1 = (2.0 * r0 * (1.0 + delta) - a * (delta**2 + kappa**2 - 1.0)) / ( - 2.0 * (1.0 + delta) - ) - x2 = (2.0 * r0 * (delta - 1.0) - a * (delta**2 + kappa**2 - 1.0)) / ( - 2.0 * (delta - 1.0) + pg = plasma_geometry( + r_0=r_0, + a=a, + triang_95=triang_95, + kappa_95=kappa_95, + i_single_null=i_single_null, ) - r1 = 0.5 * math.sqrt( - (a**2 * ((delta + 1.0) ** 2 + kappa**2) ** 2) / ((delta + 1.0) ** 2) - ) - r2 = 0.5 * math.sqrt( - (a**2 * ((delta - 1.0) ** 2 + kappa**2) ** 2) / ((delta - 1.0) ** 2) - ) - theta1 = np.arcsin((kappa * a) / r1) - theta2 = np.arcsin((kappa * a) / r2) - inang = 1.0 / r1 - outang = 1.5 / r2 - if i_single_null == 0: - angs1 = np.linspace( - -(inang + theta1) + np.pi, (inang + theta1) + np.pi, 256, endpoint=True - ) - angs2 = np.linspace(-(outang + theta2), (outang + theta2), 256, endpoint=True) - elif i_single_null < 0: - angs1 = np.linspace( - -(inang + theta1) + np.pi, theta1 + np.pi, 256, endpoint=True - ) - angs2 = np.linspace(-theta2, (outang + theta2), 256, endpoint=True) - else: - angs1 = np.linspace( - -theta1 + np.pi, (inang + theta1) + np.pi, 256, endpoint=True - ) - angs2 = np.linspace(-(outang + theta2), theta2, 256, endpoint=True) - - xs1 = -(r1 * np.cos(angs1) - x1) - ys1 = r1 * np.sin(angs1) - xs2 = -(r2 * np.cos(angs2) - x2) - ys2 = r2 * np.sin(angs2) - axis.plot(xs1, ys1, color="black") - axis.plot(xs2, ys2, color="black") + + axis.plot(pg.rs[0], pg.zs[0], color="black") + axis.plot(pg.rs[1], pg.zs[1], color="black") axis.fill_betweenx( - ys1, - xs1, - xs2, - where=(xs2 < xs1) & (ys1 > (-a * kappa)) & (ys1 < (a * kappa)), + pg.zs[0], + pg.rs[0], + pg.rs[1], + where=(pg.rs[1] < pg.rs[0]) + & (pg.zs[0] > (-a * pg.kappa)) + & (pg.zs[0] < (a * pg.kappa)), color=plasma, ) - axis.fill_betweenx(ys1, xs1, xs2, where=(xs2 > xs1), color="none") + axis.fill_betweenx( + pg.zs[0], pg.rs[0], pg.rs[1], where=(pg.rs[1] > pg.rs[0]), color="none" + ) def plot_centre_cross(axis, mfile_data, scan): @@ -385,19 +309,18 @@ def poloidal_cross_section(axis, mfile_data, scan, demo_ranges): def plot_cryostat(axis, mfile_data, scan): """Function to plot cryostat in poloidal cross-section""" - rect = patches.Rectangle( - [rdewex, 0], ddwex, zdewex + ddwex, lw=0, facecolor=cryostat - ) - axis.add_patch(rect) - rect = patches.Rectangle( - [rdewex, 0], ddwex, -(zdewex + ddwex), lw=0, facecolor=cryostat - ) - axis.add_patch(rect) - rect = patches.Rectangle([0, zdewex], rdewex, ddwex, lw=0, facecolor=cryostat) - axis.add_patch(rect) - rect = patches.Rectangle([0, -zdewex], rdewex, -ddwex, lw=0, facecolor=cryostat) - axis.add_patch(rect) + rects = cryostat_geometry(rdewex=rdewex, ddwex=ddwex, zdewex=zdewex) + + for rec in rects: + axis.add_patch( + patches.Rectangle( + xy=(rec.anchor_x, rec.anchor_z), + width=rec.width, + height=rec.height, + facecolor=cryostat, + ) + ) def color_key(axis): @@ -728,36 +651,6 @@ def arc_fill(axis, r1, r2, color="pink"): axis.add_patch(patch) -def ellips_fill( - axis, a1=0, a2=0, b1=0, b2=0, x0=0, y0=0, ang1=0, ang2=rtangle, color="pink" -): - """Fills the space between two concentric ellipse sectors. - - Arguments - - axis: plot object - a1, a2, b1, b2 horizontal and vertical radii to be filled - x0, y0 coordinates of centre of the ellipses - ang1, ang2 are the polar angles of the start and end - - """ - angs = np.linspace(ang1, ang2, endpoint=True) - r1 = ((np.cos(angs) / a1) ** 2 + (np.sin(angs) / b1) ** 2) ** (-0.5) - xs1 = r1 * np.cos(angs) + x0 - ys1 = r1 * np.sin(angs) + y0 - angs = np.linspace(ang2, ang1, endpoint=True) - r2 = ((np.cos(angs) / a2) ** 2 + (np.sin(angs) / b2) ** 2) ** (-0.5) - xs2 = r2 * np.cos(angs) + x0 - ys2 = r2 * np.sin(angs) + y0 - verts = list(zip(xs1, ys1)) - verts.extend(list(zip(xs2, ys2))) - endpoint = verts[-1:] - verts.extend(endpoint) - path = Path(verts, closed=True) - patch = patches.PathPatch(path, facecolor=color, lw=0) - axis.add_patch(patch) - - def plot_nprofile(prof, demo_ranges): """Function to plot density profile Arguments: @@ -1227,62 +1120,71 @@ def plot_vacuum_vessel(axis, mfile_data, scan): scan --> scan number to use """ i_single_null = mfile_data.data["i_single_null"].get_scan(scan) - triang = mfile_data.data["triang95"].get_scan(scan) - temp_array_1 = () - temp_array_2 = () + triang_95 = mfile_data.data["triang95"].get_scan(scan) # Outer side (furthest from plasma) - radx = ( + radx_outer = ( cumulative_radial_build("d_vv_out", mfile_data, scan) + cumulative_radial_build("gapds", mfile_data, scan) ) / 2.0 - rminx = ( + rminx_outer = ( cumulative_radial_build("d_vv_out", mfile_data, scan) - cumulative_radial_build("gapds", mfile_data, scan) ) / 2.0 - kapx = cumulative_upper["d_vv_top"] / rminx - - if i_single_null == 1: - (rs, zs) = plotdh(axis, radx, rminx, triang, kapx) - temp_array_1 = temp_array_1 + ((rs, zs)) - - kapx = cumulative_lower["d_vv_bot"] / rminx - (rs, zs) = plotdh(axis, radx, rminx, triang, kapx) - temp_array_2 = temp_array_2 + ((rs, zs)) - # Inner side (nearest to the plasma) - radx = ( + radx_inner = ( cumulative_radial_build("shldoth", mfile_data, scan) + cumulative_radial_build("d_vv_in", mfile_data, scan) ) / 2.0 - rminx = ( + rminx_inner = ( cumulative_radial_build("shldoth", mfile_data, scan) - cumulative_radial_build("d_vv_in", mfile_data, scan) ) / 2.0 if i_single_null == 1: - kapx = (cumulative_upper["d_vv_top"] - upper["d_vv_top"]) / rminx - (rs, zs) = plotdh(axis, radx, rminx, triang, kapx) - temp_array_1 = temp_array_1 + ((rs, zs)) + vvg_single_null = vacuum_vessel_geometry_single_null( + cumulative_upper=cumulative_upper, + upper=upper, + triang=triang_95, + radx_outer=radx_outer, + rminx_outer=rminx_outer, + radx_inner=radx_inner, + rminx_inner=rminx_inner, + cumulative_lower=cumulative_lower, + lower=lower, + ) - kapx = (cumulative_lower["d_vv_bot"] + lower["d_vv_bot"]) / rminx - (rs, zs) = plotdh(axis, radx, rminx, triang, kapx) - temp_array_2 = temp_array_2 + ((rs, zs)) + axis.plot( + vvg_single_null.rs, + vvg_single_null.zs, + color="black", + lw=thin, + ) + + axis.fill( + vvg_single_null.rs, + vvg_single_null.zs, + color=vessel, + ) - # Single null: Draw top half from output - # Double null: Reflect bottom half to top - if i_single_null == 1: - rs = np.concatenate([temp_array_1[0], temp_array_1[2][::-1]]) - zs = np.concatenate([temp_array_1[1], temp_array_1[3][::-1]]) - axis.fill(rs, zs, color=vessel) - - rs = np.concatenate([temp_array_2[0], temp_array_2[2][::-1]]) - zs = np.concatenate([temp_array_2[1], temp_array_2[3][::-1]]) - axis.fill(rs, zs, color=vessel) - # For double null, reflect shape of lower half to top instead if i_single_null == 0: - axis.fill(rs, -zs, color=vessel) + vvg_double_null = vacuum_vessel_geometry_double_null( + cumulative_lower=cumulative_lower, + lower=lower, + radx_inner=radx_inner, + radx_outer=radx_outer, + rminx_inner=rminx_inner, + rminx_outer=rminx_outer, + triang=triang_95, + ) + axis.plot(vvg_double_null.rs, vvg_double_null.zs, color="black", lw=thin) + + axis.fill( + vvg_double_null.rs, + vvg_double_null.zs, + color=vessel, + ) def plot_shield(axis, mfile_data, scan): @@ -1294,63 +1196,55 @@ def plot_shield(axis, mfile_data, scan): scan --> scan number to use """ i_single_null = mfile_data.data["i_single_null"].get_scan(scan) - triang = mfile_data.data["triang95"].get_scan(scan) - temp_array_1 = () - temp_array_2 = () + triang_95 = mfile_data.data["triang95"].get_scan(scan) # Side furthest from plasma - radx = ( + radx_far = ( cumulative_radial_build("shldoth", mfile_data, scan) + cumulative_radial_build("d_vv_in", mfile_data, scan) ) / 2.0 - rminx = ( + rminx_far = ( cumulative_radial_build("shldoth", mfile_data, scan) - cumulative_radial_build("d_vv_in", mfile_data, scan) ) / 2.0 - if i_single_null == 1: - kapx = cumulative_upper["shldtth"] / rminx - (rs, zs) = plotdh(axis, radx, rminx, triang, kapx) - temp_array_1 = temp_array_1 + ((rs, zs)) - - kapx = cumulative_lower["shldlth"] / rminx - (rs, zs) = plotdh(axis, radx, rminx, triang, kapx) - temp_array_2 = temp_array_2 + ((rs, zs)) - # Side nearest to the plasma - radx = ( + radx_near = ( cumulative_radial_build("vvblgapo", mfile_data, scan) + cumulative_radial_build("shldith", mfile_data, scan) ) / 2.0 - rminx = ( + rminx_near = ( cumulative_radial_build("vvblgapo", mfile_data, scan) - cumulative_radial_build("shldith", mfile_data, scan) ) / 2.0 if i_single_null == 1: - kapx = (cumulative_upper["vvblgap"]) / rminx - (rs, zs) = plotdh(axis, radx, rminx, triang, kapx) - temp_array_1 = temp_array_1 + ((rs, zs)) - - kapx = (cumulative_lower["divfix"]) / rminx - (rs, zs) = plotdh(axis, radx, rminx, triang, kapx) - temp_array_2 = temp_array_2 + ((rs, zs)) - - # Single null: Draw top half from output - # Double null: Reflect bottom half to top - if i_single_null == 1: - rs = np.concatenate([temp_array_1[0], temp_array_1[2][::-1]]) - zs = np.concatenate([temp_array_1[1], temp_array_1[3][::-1]]) - axis.fill(rs, zs, color=shield) + sg_single_null = shield_geometry_single_null( + cumulative_upper=cumulative_upper, + radx_far=radx_far, + rminx_far=rminx_far, + radx_near=radx_near, + rminx_near=rminx_near, + triang=triang_95, + cumulative_lower=cumulative_lower, + ) + axis.plot(sg_single_null.rs, sg_single_null.zs, color="black", lw=thin) + axis.fill(sg_single_null.rs, sg_single_null.zs, color=shield) - rs = np.concatenate([temp_array_2[0], temp_array_2[2][::-1]]) - zs = np.concatenate([temp_array_2[1], temp_array_2[3][::-1]]) - axis.fill(rs, zs, color=shield) if i_single_null == 0: - axis.fill(rs, -zs, color=shield) + sg_double_null = shield_geometry_double_null( + cumulative_lower=cumulative_lower, + radx_far=radx_far, + radx_near=radx_near, + rminx_far=rminx_far, + rminx_near=rminx_near, + triang=triang_95, + ) + axis.plot(sg_double_null.rs, sg_double_null.zs, color="black", lw=thin) + axis.fill(sg_double_null.rs, sg_double_null.zs, color=shield) -def plot_blanket(axis, mfile_data, scan): +def plot_blanket(axis, mfile_data, scan) -> None: """Function to plot blanket Arguments: @@ -1359,65 +1253,74 @@ def plot_blanket(axis, mfile_data, scan): scan --> scan number to use """ - point_array = () - # Single null: Draw top half from output # Double null: Reflect bottom half to top i_single_null = mfile_data.data["i_single_null"].get_scan(scan) + triang_95 = mfile_data.data["triang95"].get_scan(scan) + blnktth = mfile_data.data["blnktth"].get_scan(scan) + c_shldith = cumulative_radial_build("shldith", mfile_data, scan) + c_blnkoth = cumulative_radial_build("blnkoth", mfile_data, scan) + if i_single_null == 1: # Upper blanket: outer surface - radx = ( + radx_outer = ( cumulative_radial_build("blnkoth", mfile_data, scan) + cumulative_radial_build("vvblgapi", mfile_data, scan) ) / 2.0 - rminx = ( + rminx_outer = ( cumulative_radial_build("blnkoth", mfile_data, scan) - cumulative_radial_build("vvblgapi", mfile_data, scan) ) / 2.0 - kapx = cumulative_upper["blnktth"] / rminx - (rs, zs) = plotdh(axis, radx, rminx, triang, kapx) - point_array = point_array + ((rs, zs)) - # Upper blanket: inner surface - radx = ( + radx_inner = ( cumulative_radial_build("fwoth", mfile_data, scan) + cumulative_radial_build("blnkith", mfile_data, scan) ) / 2.0 - rminx = ( + rminx_inner = ( cumulative_radial_build("fwoth", mfile_data, scan) - cumulative_radial_build("blnkith", mfile_data, scan) ) / 2.0 + bg_single_null = blanket_geometry_single_null( + radx_outer=radx_outer, + rminx_outer=rminx_outer, + radx_inner=radx_inner, + rminx_inner=rminx_inner, + cumulative_upper=cumulative_upper, + triang=triang_95, + cumulative_lower=cumulative_lower, + blnktth=blnktth, + c_shldith=c_shldith, + c_blnkoth=c_blnkoth, + blnkith=blnkith, + blnkoth=blnkoth, + ) - kapx = cumulative_upper["fwtth"] / rminx - (rs, zs) = plotdh(axis, radx, rminx, triang, kapx) - point_array = point_array + ((rs, zs)) + # Plot blanket + axis.plot( + bg_single_null.rs, + bg_single_null.zs, + color="black", + lw=thin, + ) - # Plot upper blanket - rs = np.concatenate([point_array[0], point_array[2][::-1]]) - zs = np.concatenate([point_array[1], point_array[3][::-1]]) - axis.fill(rs, zs, color=blanket) + axis.fill(bg_single_null.rs, bg_single_null.zs, color=blanket) - # Lower blanket - blnktth = mfile_data.data["blnktth"].get_scan(scan) - c_shldith = cumulative_radial_build("shldith", mfile_data, scan) - c_blnkoth = cumulative_radial_build("blnkoth", mfile_data, scan) - divgap = cumulative_lower["divfix"] - plotdhgap( - axis, c_shldith, c_blnkoth, blnkith, blnkoth, divgap, -blnktth, triang, blanket - ) if i_single_null == 0: - plotdhgap( - axis, - c_shldith, - c_blnkoth, - blnkith, - blnkoth, - -divgap, - -blnktth, - triang, - blanket, + bg_double_null = blanket_geometry_double_null( + cumulative_lower=cumulative_lower, + triang=triang_95, + blnktth=blnktth, + c_shldith=c_shldith, + c_blnkoth=c_blnkoth, + blnkith=blnkith, + blnkoth=blnkoth, ) + # Plot blanket + axis.plot(bg_double_null.rs[0], bg_double_null.zs[0], color="black", lw=thin) + axis.plot(bg_double_null.rs[1], bg_double_null.zs[1], color="black", lw=thin) + axis.fill(bg_double_null.rs[0], bg_double_null.zs[0], color=blanket) + axis.fill(bg_double_null.rs[1], bg_double_null.zs[1], color=blanket) def plot_firstwall(axis, mfile_data, scan): @@ -1429,75 +1332,70 @@ def plot_firstwall(axis, mfile_data, scan): scan --> scan number to use """ + i_single_null = mfile_data.data["i_single_null"].get_scan(scan) + triang_95 = mfile_data.data["triang95"].get_scan(scan) blnktth = mfile_data.data["blnktth"].get_scan(scan) tfwvt = mfile_data.data["fwtth"].get_scan(scan) - i_single_null = mfile_data.data["i_single_null"].get_scan(scan) - point_array = () + c_blnkith = cumulative_radial_build("blnkith", mfile_data, scan) + c_fwoth = cumulative_radial_build("fwoth", mfile_data, scan) - # Single null: Draw top half from output - # Double null: Reflect bottom half to top if i_single_null == 1: # Upper first wall: outer surface - radx = ( + radx_outer = ( cumulative_radial_build("fwoth", mfile_data, scan) + cumulative_radial_build("blnkith", mfile_data, scan) ) / 2.0 - rminx = ( + rminx_outer = ( cumulative_radial_build("fwoth", mfile_data, scan) - cumulative_radial_build("blnkith", mfile_data, scan) ) / 2.0 - kapx = cumulative_upper["fwtth"] / rminx - (rs, zs) = plotdh(axis, radx, rminx, triang, kapx) - point_array = point_array + ((rs, zs)) - # Upper first wall: inner surface - radx = ( + radx_inner = ( cumulative_radial_build("scraplo", mfile_data, scan) + cumulative_radial_build("fwith", mfile_data, scan) ) / 2.0 - rminx = ( + rminx_inner = ( cumulative_radial_build("scraplo", mfile_data, scan) - cumulative_radial_build("fwith", mfile_data, scan) ) / 2.0 - (rs, zs) = plotdh(axis, radx, rminx, triang, kapx) - point_array = point_array + ((rs, zs)) - - # Plot upper first wall - rs = np.concatenate([point_array[0], point_array[2][::-1]]) - zs = np.concatenate([point_array[1], point_array[3][::-1]]) - axis.fill(rs, zs, color=firstwall) + fwg_single_null = first_wall_geometry_single_null( + radx_outer=radx_outer, + rminx_outer=rminx_outer, + radx_inner=radx_inner, + rminx_inner=rminx_inner, + cumulative_upper=cumulative_upper, + triang=triang_95, + cumulative_lower=cumulative_lower, + blnktth=blnktth, + c_blnkith=c_blnkith, + c_fwoth=c_fwoth, + fwith=fwith, + fwoth=fwoth, + tfwvt=tfwvt, + ) - # Lower first wall - c_blnkith = cumulative_radial_build("blnkith", mfile_data, scan) - c_fwoth = cumulative_radial_build("fwoth", mfile_data, scan) - divgap = cumulative_lower["divfix"] - - plotdhgap( - axis, - c_blnkith, - c_fwoth, - fwith, - fwoth, - divgap + blnktth, - -tfwvt, - triang, - firstwall, - ) + # Plot first wall + axis.plot(fwg_single_null.rs, fwg_single_null.zs, color="black", lw=thin) + axis.fill(fwg_single_null.rs, fwg_single_null.zs, color=firstwall) if i_single_null == 0: - plotdhgap( - axis, - c_blnkith, - c_fwoth, - fwith, - fwoth, - -(divgap + blnktth), - -tfwvt, - triang, - firstwall, + fwg_double_null = first_wall_geometry_double_null( + cumulative_lower=cumulative_lower, + triang=triang_95, + blnktth=blnktth, + c_blnkith=c_blnkith, + c_fwoth=c_fwoth, + fwith=fwith, + fwoth=fwoth, + tfwvt=tfwvt, ) + # Plot blanket + axis.plot(fwg_double_null.rs[0], fwg_double_null.zs[0], color="black", lw=thin) + axis.plot(fwg_double_null.rs[1], fwg_double_null.zs[1], color="black", lw=thin) + axis.fill(fwg_double_null.rs[0], fwg_double_null.zs[0], color=blanket) + axis.fill(fwg_double_null.rs[1], fwg_double_null.zs[1], color=blanket) def angle_check(angle1, angle2): @@ -1543,116 +1441,59 @@ def plot_tf_coils(axis, mfile_data, scan): else: i_tf_shape = int(1) - # D-shaped (i_tf_shape=1), Picture frame (i_tf_shape=2) if i_tf_shape == 2: - # Inboard leg - rect1 = patches.Rectangle( - [x5 - tfcth, y5 - tfcth], - tfcth, - (y1 - y5 + 2.0 * tfcth), - lw=0, - facecolor="cyan", - ) - # Outboard leg vertical - rect2 = patches.Rectangle( - [x4, y4 - tfcth], tfcth, (y2 - y4 + 2.0 * tfcth), lw=0, facecolor="cyan" + rects = tfcoil_geometry_rectangular_shape( + x1=x1, + x2=x2, + x4=x4, + x5=x5, + y1=y1, + y2=y2, + y4=y4, + y5=y5, + tfcth=tfcth, ) - # Outboard leg horizontal bottom - rect3 = patches.Rectangle( - [x5, y5 - tfcth], x4 - x5, tfcth, lw=0, facecolor="cyan" - ) - # Outboard leg horizontal top - rect4 = patches.Rectangle([x1, y1], x2 - x1, tfcth, lw=0, facecolor="cyan") - # Plot it all - axis.add_patch(rect1) - axis.add_patch(rect2) - axis.add_patch(rect3) - axis.add_patch(rect4) + for rec in rects: + axis.add_patch( + patches.Rectangle( + xy=(rec.anchor_x, rec.anchor_z), + width=rec.width, + height=rec.height, + facecolor=tfc, + ) + ) else: - # Inboard upper arc - x0 = x2 - y0 = y1 - a1 = x2 - x1 - b1 = y2 - y1 - a2 = a1 + tfcth - b2 = b1 + tfcth - ellips_fill( - axis, - a1=a1, - a2=a2, - b1=b1, - b2=b2, - x0=x0, - y0=y0, - ang1=rtangle, - ang2=2 * rtangle, - color="cyan", - ) - # Outboard upper arc - x0 = x2 - y0 = 0 - a1 = x3 - x2 - b1 = y2 - a2 = a1 + tfcth - b2 = b1 + tfcth - ellips_fill( - axis, - a1=a1, - a2=a2, - b1=b1, - b2=b2, - x0=x0, - y0=y0, - ang1=0, - ang2=rtangle, - color="cyan", - ) - # Inboard lower arc - x0 = x4 - y0 = y5 - a1 = x4 - x5 - b1 = y5 - y4 - a2 = a1 + tfcth - b2 = b1 + tfcth - ellips_fill( - axis, - a1=a1, - a2=a2, - b1=b1, - b2=b2, - x0=x0, - y0=y0, - ang1=-rtangle, - ang2=-2 * rtangle, - color="cyan", - ) - # Outboard lower arc - x0 = x4 - y0 = 0 - a1 = x3 - x2 - b1 = -y4 - a2 = a1 + tfcth - b2 = b1 + tfcth - ellips_fill( - axis, - a1=a1, - a2=a2, - b1=b1, - b2=b2, - x0=x0, - y0=y0, - ang1=0, - ang2=-rtangle, - color="cyan", + rects, verts = tfcoil_geometry_d_shape( + x1=x1, + x2=x2, + x3=x3, + x4=x4, + x5=x5, + y1=y1, + y2=y2, + y4=y4, + y5=y5, + tfcth=tfcth, + rtangle=rtangle, + rtangle2=rtangle2, ) - # Vertical leg - # Bottom left corner - rect = patches.Rectangle( - [x5 - tfcth, y5], tfcth, (y1 - y5), lw=0, facecolor="cyan" - ) - axis.add_patch(rect) + + for vert in verts: + path = Path(vert, closed=True) + patch = patches.PathPatch(path, facecolor=tfc, lw=0) + axis.add_patch(patch) + + for rec in rects: + axis.add_patch( + patches.Rectangle( + xy=(rec.anchor_x, rec.anchor_z), + width=rec.width, + height=rec.height, + facecolor=tfc, + ) + ) def plot_pf_coils(axis, mfile_data, scan): @@ -1670,16 +1511,16 @@ def plot_pf_coils(axis, mfile_data, scan): coils_dz = [] coil_text = [] + bore = mfile_data.data["bore"].get_scan(scan) + ohcth = mfile_data.data["ohcth"].get_scan(scan) + ohdz = mfile_data.data["ohdz"].get_scan(scan) + # Number of coils (1 is OH coil) number_of_coils = 0 for item in mfile_data.data.keys(): if "rpf[" in item: number_of_coils += 1 - bore = mfile_data.data["bore"].get_scan(scan) - ohcth = mfile_data.data["ohcth"].get_scan(scan) - ohdz = mfile_data.data["ohdz"].get_scan(scan) - # Check for Central Solenoid if "iohcl" in mfile_data.data.keys(): iohcl = mfile_data.data["iohcl"].get_scan(scan) @@ -1700,23 +1541,18 @@ def plot_pf_coils(axis, mfile_data, scan): coils_dz.append(mfile_data.data["pfdz({:01})".format(coil)].get_scan(scan)) coil_text.append(str(coil + 1)) + r_points, z_points, central_coil = pfcoil_geometry( + coils_r=coils_r, + coils_z=coils_z, + coils_dr=coils_dr, + coils_dz=coils_dz, + bore=bore, + ohcth=ohcth, + ohdz=ohdz, + ) + for i in range(len(coils_r)): - r_1 = coils_r[i] - 0.5 * coils_dr[i] - z_1 = coils_z[i] - 0.5 * coils_dz[i] - r_2 = coils_r[i] - 0.5 * coils_dr[i] - z_2 = coils_z[i] + 0.5 * coils_dz[i] - r_3 = coils_r[i] + 0.5 * coils_dr[i] - z_3 = coils_z[i] + 0.5 * coils_dz[i] - r_4 = coils_r[i] + 0.5 * coils_dr[i] - z_4 = coils_z[i] - 0.5 * coils_dz[i] - r_5 = coils_r[i] - 0.5 * coils_dr[i] - z_5 = coils_z[i] - 0.5 * coils_dz[i] - - r_points = [r_1, r_2, r_3, r_4, r_5] - z_points = [z_1, z_2, z_3, z_4, z_5] - axis.plot(r_points, z_points, color="black") - rect = patches.Rectangle([bore, -ohdz / 2], ohcth, ohdz, lw=0, facecolor="pink") - axis.add_patch(rect) + axis.plot(r_points[i], z_points[i], color="black") axis.text( coils_r[i], coils_z[i], @@ -1725,6 +1561,14 @@ def plot_pf_coils(axis, mfile_data, scan): va="center", fontsize=5 * abs((coils_dr[i] * coils_dz[i])), ) + axis.add_patch( + patches.Rectangle( + xy=(central_coil.anchor_x, central_coil.anchor_z), + width=central_coil.width, + height=central_coil.height, + facecolor="pink", + ) + ) def plot_info(axis, data, mfile_data, scan): diff --git a/tests/integration/data/large_tokamak_IN.DAT b/tests/integration/data/large_tokamak_IN.DAT index 65094ff6ec..4acbc76b48 100644 --- a/tests/integration/data/large_tokamak_IN.DAT +++ b/tests/integration/data/large_tokamak_IN.DAT @@ -310,7 +310,7 @@ d_vv_top = 0.3 * Underside vacuum vessel radial thickness [m] d_vv_bot = 0.3 -* Inboard vacuum vessel thickness [m] +* Inboard radiation shield radial thickness [m] shldith = 0.3 * Gap between vacuum vessel and blanket [m] diff --git a/tests/integration/data/large_tokamak_MFILE.DAT b/tests/integration/data/large_tokamak_MFILE.DAT index afdf9c3afa..0644c3c468 100644 --- a/tests/integration/data/large_tokamak_MFILE.DAT +++ b/tests/integration/data/large_tokamak_MFILE.DAT @@ -2,15 +2,15 @@ # Power Reactor Optimisation Code # # PROCESS # # Power Reactor Optimisation Code # - PROCESS_version_number__________________________________________________ (procver)_____________________ "2.4.0 R" - Date_of_run_____________________________________________________________ (date)________________________ "30/06/2023" - Time_of_run_____________________________________________________________ (time)________________________ "11:21" - User____________________________________________________________________ (username)____________________ "timothyn" + PROCESS_version_number__________________________________________________ (procver)_____________________ "3.0.1 R" + Date_of_run_____________________________________________________________ (date)________________________ "11/12/2023" + Time_of_run_____________________________________________________________ (time)________________________ "13:46" + User____________________________________________________________________ (username)____________________ "clair" PROCESS_run_title_______________________________________________________ (runtitle)____________________ "Generic large tokamak - PROCESS_tag_number______________________________________________________ (tagno)_______________________ "v2.4.0-1153-gdf5180f97" - PROCESS_git_branch_name_________________________________________________ (branch_name)_________________ "vv_stress" - PROCESS_last_commit_message_____________________________________________ (commsg)______________________ "Small changes to the VV stress model" - Input_filename__________________________________________________________ (fileprefix)__________________ "/Users/timothynunn/process/tests/regression/scenarios/large-tokamak/IN.DAT " + PROCESS_tag_number______________________________________________________ (tagno)_______________________ "v3.0.1-19-ge6b4c55" + PROCESS_git_branch_name_________________________________________________ (branch_name)_________________ "cmould/updating-plot-proc" + PROCESS_last_commit_message_____________________________________________ (commsg)______________________ "removed todo comments" + Input_filename__________________________________________________________ (fileprefix)__________________ "/home/clair/development/PROCESS/tests/integration/data/large_tokamak_IN.DAT " # Numerics # VMCON_error_flag________________________________________________________ (ifail)_______________________ 1 # PROCESS found a feasible solution # @@ -18,192 +18,192 @@ Number_of_constraints_(total)___________________________________________ (neqns+nineqns)_______________ 26 Optimisation_switch_____________________________________________________ (ioptimz)_____________________ 1 Figure_of_merit_switch__________________________________________________ (minmax)______________________ 1 - Square_root_of_the_sum_of_squares_of_the_constraint_residuals___________ (sqsumsq)_____________________ 3.6555E-08 OP - VMCON_convergence_parameter_____________________________________________ (convergence_parameter)_______ 3.7075E-17 OP + Square_root_of_the_sum_of_squares_of_the_constraint_residuals___________ (sqsumsq)_____________________ 5.3805E-05 OP + VMCON_convergence_parameter_____________________________________________ (convergence_parameter)_______ 0.0000E+00 OP Normalised_objective_function___________________________________________ (norm_objf)___________________ 1.6000E+00 OP - Number_of_VMCON_iterations______________________________________________ (nviter)______________________ 11 OP - beta____________________________________________________________________ (itvar001)____________________ 3.6055E-02 - beta_(final_value/initial_value)________________________________________ (xcm001)______________________ 1.2018E+00 - beta_(range_normalised)_________________________________________________ (nitvar001)___________________ 3.5090E-02 - dene____________________________________________________________________ (itvar002)____________________ 8.0995E+19 - dene_(final_value/initial_value)________________________________________ (xcm002)______________________ 1.0799E+00 - dene_(range_normalised)_________________________________________________ (nitvar002)___________________ 7.1712E-02 - fwalld__________________________________________________________________ (itvar003)____________________ 4.9896E-01 - fwalld_(final_value/initial_value)______________________________________ (xcm003)______________________ 4.9896E-01 - fwalld_(range_normalised)_______________________________________________ (nitvar003)___________________ 4.9846E-01 - ffuspow_________________________________________________________________ (itvar004)____________________ 5.3088E-01 - ffuspow_(final_value/initial_value)_____________________________________ (xcm004)______________________ 5.3088E-01 - ffuspow_(range_normalised)______________________________________________ (nitvar004)___________________ 5.3041E-01 - ftburn__________________________________________________________________ (itvar005)____________________ 1.0000E+00 - ftburn_(final_value/initial_value)______________________________________ (xcm005)______________________ 1.0000E+00 - ftburn_(range_normalised)_______________________________________________ (nitvar005)___________________ 1.0000E+00 - flhthresh_______________________________________________________________ (itvar006)____________________ 1.6436E+00 - flhthresh_(final_value/initial_value)___________________________________ (xcm006)______________________ 1.6436E+00 - flhthresh_(range_normalised)____________________________________________ (nitvar006)___________________ 7.1515E-02 - fpinj___________________________________________________________________ (itvar007)____________________ 3.7606E-01 - fpinj_(final_value/initial_value)_______________________________________ (xcm007)______________________ 3.7606E-01 - fpinj_(range_normalised)________________________________________________ (nitvar007)___________________ 3.7544E-01 + Number_of_VMCON_iterations______________________________________________ (nviter)______________________ 0 OP + beta____________________________________________________________________ (itvar001)____________________ 3.4120E-02 + beta_(final_value/initial_value)________________________________________ (xcm001)______________________ 1.1373E+00 + beta_(range_normalised)_________________________________________________ (nitvar001)___________________ 3.3154E-02 + dene____________________________________________________________________ (itvar002)____________________ 8.0594E+19 + dene_(final_value/initial_value)________________________________________ (xcm002)______________________ 1.0746E+00 + dene_(range_normalised)_________________________________________________ (nitvar002)___________________ 7.1307E-02 + fwalld__________________________________________________________________ (itvar003)____________________ 5.0034E-01 + fwalld_(final_value/initial_value)______________________________________ (xcm003)______________________ 5.0034E-01 + fwalld_(range_normalised)_______________________________________________ (nitvar003)___________________ 4.9984E-01 + ffuspow_________________________________________________________________ (itvar004)____________________ 5.3236E-01 + ffuspow_(final_value/initial_value)_____________________________________ (xcm004)______________________ 5.3236E-01 + ffuspow_(range_normalised)______________________________________________ (nitvar004)___________________ 5.3189E-01 + ftburn__________________________________________________________________ (itvar005)____________________ 9.9849E-01 + ftburn_(final_value/initial_value)______________________________________ (xcm005)______________________ 9.9849E-01 + ftburn_(range_normalised)_______________________________________________ (nitvar005)___________________ 9.9849E-01 + flhthresh_______________________________________________________________ (itvar006)____________________ 1.6284E+00 + flhthresh_(final_value/initial_value)___________________________________ (xcm006)______________________ 1.6284E+00 + flhthresh_(range_normalised)____________________________________________ (nitvar006)___________________ 6.9821E-02 + fpinj___________________________________________________________________ (itvar007)____________________ 3.7699E-01 + fpinj_(final_value/initial_value)_______________________________________ (xcm007)______________________ 3.7699E-01 + fpinj_(range_normalised)________________________________________________ (nitvar007)___________________ 3.7637E-01 fpnetel_________________________________________________________________ (itvar008)____________________ 1.0000E+00 fpnetel_(final_value/initial_value)_____________________________________ (xcm008)______________________ 1.0000E+00 fpnetel_(range_normalised)______________________________________________ (nitvar008)___________________ 1.0000E+00 - fbetatry________________________________________________________________ (itvar009)____________________ 5.3636E-01 - fbetatry_(final_value/initial_value)____________________________________ (xcm009)______________________ 1.0727E+00 - fbetatry_(range_normalised)_____________________________________________ (nitvar009)___________________ 5.3589E-01 - fpeakb__________________________________________________________________ (itvar010)____________________ 8.3188E-01 - fpeakb_(final_value/initial_value)______________________________________ (xcm010)______________________ 8.3188E-01 - fpeakb_(range_normalised)_______________________________________________ (nitvar010)___________________ 8.3171E-01 - coheof__________________________________________________________________ (itvar011)____________________ 2.0066E+07 - coheof_(final_value/initial_value)______________________________________ (xcm011)______________________ 1.3377E+00 - coheof_(range_normalised)_______________________________________________ (nitvar011)___________________ 1.9986E-01 - fjohc___________________________________________________________________ (itvar012)____________________ 4.7704E-01 - fjohc_(final_value/initial_value)_______________________________________ (xcm012)______________________ 7.9506E-01 - fjohc_(range_normalised)________________________________________________ (nitvar012)___________________ 4.7176E-01 - fjohc0__________________________________________________________________ (itvar013)____________________ 5.6968E-01 - fjohc0_(final_value/initial_value)______________________________________ (xcm013)______________________ 9.4947E-01 - fjohc0_(range_normalised)_______________________________________________ (nitvar013)___________________ 5.6925E-01 - fcohbop_________________________________________________________________ (itvar014)____________________ 9.8528E-01 - fcohbop_(final_value/initial_value)_____________________________________ (xcm014)______________________ 1.0948E+00 - fcohbop_(range_normalised)______________________________________________ (nitvar014)___________________ 9.8527E-01 - fiooic__________________________________________________________________ (itvar015)____________________ 7.3084E-01 - fiooic_(final_value/initial_value)______________________________________ (xcm015)______________________ 1.1244E+00 - fiooic_(range_normalised)_______________________________________________ (nitvar015)___________________ 7.3057E-01 - fvdump__________________________________________________________________ (itvar016)____________________ 9.3438E-01 - fvdump_(final_value/initial_value)______________________________________ (xcm016)______________________ 9.3438E-01 - fvdump_(range_normalised)_______________________________________________ (nitvar016)___________________ 9.3431E-01 - vdalw___________________________________________________________________ (itvar017)____________________ 9.3438E+00 - vdalw_(final_value/initial_value)_______________________________________ (xcm017)______________________ 9.3438E-01 - vdalw_(range_normalised)________________________________________________ (nitvar017)___________________ 9.3437E-01 - fjprot__________________________________________________________________ (itvar018)____________________ 1.0000E+00 - fjprot_(final_value/initial_value)______________________________________ (xcm018)______________________ 1.0000E+00 - fjprot_(range_normalised)_______________________________________________ (nitvar018)___________________ 1.0000E+00 - ftmargtf________________________________________________________________ (itvar019)____________________ 9.9341E-01 - ftmargtf_(final_value/initial_value)____________________________________ (xcm019)______________________ 9.9341E-01 - ftmargtf_(range_normalised)_____________________________________________ (nitvar019)___________________ 9.9340E-01 - ftmargoh________________________________________________________________ (itvar020)____________________ 1.0000E+00 - ftmargoh_(final_value/initial_value)____________________________________ (xcm020)______________________ 1.0000E+00 - ftmargoh_(range_normalised)_____________________________________________ (nitvar020)___________________ 1.0000E+00 - ftaulimit_______________________________________________________________ (itvar021)____________________ 7.2549E-01 - ftaulimit_(final_value/initial_value)___________________________________ (xcm021)______________________ 7.2549E-01 - ftaulimit_(range_normalised)____________________________________________ (nitvar021)___________________ 7.2521E-01 - ftaucq__________________________________________________________________ (itvar022)____________________ 9.1920E-01 - ftaucq_(final_value/initial_value)______________________________________ (xcm022)______________________ 1.0213E+00 - ftaucq_(range_normalised)_______________________________________________ (nitvar022)___________________ 9.1912E-01 - foh_stress______________________________________________________________ (itvar023)____________________ 9.1989E-01 - foh_stress_(final_value/initial_value)__________________________________ (xcm023)______________________ 9.1989E-01 - foh_stress_(range_normalised)___________________________________________ (nitvar023)___________________ 9.1981E-01 + fbetatry________________________________________________________________ (itvar009)____________________ 5.1462E-01 + fbetatry_(final_value/initial_value)____________________________________ (xcm009)______________________ 1.0292E+00 + fbetatry_(range_normalised)_____________________________________________ (nitvar009)___________________ 5.1414E-01 + fpeakb__________________________________________________________________ (itvar010)____________________ 8.6061E-01 + fpeakb_(final_value/initial_value)______________________________________ (xcm010)______________________ 8.6061E-01 + fpeakb_(range_normalised)_______________________________________________ (nitvar010)___________________ 8.6047E-01 + coheof__________________________________________________________________ (itvar011)____________________ 2.0403E+07 + coheof_(final_value/initial_value)______________________________________ (xcm011)______________________ 1.3602E+00 + coheof_(range_normalised)_______________________________________________ (nitvar011)___________________ 2.0323E-01 + fjohc___________________________________________________________________ (itvar012)____________________ 5.0665E-01 + fjohc_(final_value/initial_value)_______________________________________ (xcm012)______________________ 8.4442E-01 + fjohc_(range_normalised)________________________________________________ (nitvar012)___________________ 5.0167E-01 + fjohc0__________________________________________________________________ (itvar013)____________________ 5.5864E-01 + fjohc0_(final_value/initial_value)______________________________________ (xcm013)______________________ 9.3106E-01 + fjohc0_(range_normalised)_______________________________________________ (nitvar013)___________________ 5.5820E-01 + fcohbop_________________________________________________________________ (itvar014)____________________ 9.6866E-01 + fcohbop_(final_value/initial_value)_____________________________________ (xcm014)______________________ 1.0763E+00 + fcohbop_(range_normalised)______________________________________________ (nitvar014)___________________ 9.6863E-01 + fiooic__________________________________________________________________ (itvar015)____________________ 7.1738E-01 + fiooic_(final_value/initial_value)______________________________________ (xcm015)______________________ 1.1037E+00 + fiooic_(range_normalised)_______________________________________________ (nitvar015)___________________ 7.1709E-01 + fvdump__________________________________________________________________ (itvar016)____________________ 9.9851E-01 + fvdump_(final_value/initial_value)______________________________________ (xcm016)______________________ 9.9851E-01 + fvdump_(range_normalised)_______________________________________________ (nitvar016)___________________ 9.9851E-01 + vdalw___________________________________________________________________ (itvar017)____________________ 9.9851E+00 + vdalw_(final_value/initial_value)_______________________________________ (xcm017)______________________ 9.9851E-01 + vdalw_(range_normalised)________________________________________________ (nitvar017)___________________ 9.9851E-01 + fjprot__________________________________________________________________ (itvar018)____________________ 9.9818E-01 + fjprot_(final_value/initial_value)______________________________________ (xcm018)______________________ 9.9818E-01 + fjprot_(range_normalised)_______________________________________________ (nitvar018)___________________ 9.9818E-01 + ftmargtf________________________________________________________________ (itvar019)____________________ 9.8389E-01 + ftmargtf_(final_value/initial_value)____________________________________ (xcm019)______________________ 9.8389E-01 + ftmargtf_(range_normalised)_____________________________________________ (nitvar019)___________________ 9.8388E-01 + ftmargoh________________________________________________________________ (itvar020)____________________ 9.9869E-01 + ftmargoh_(final_value/initial_value)____________________________________ (xcm020)______________________ 9.9869E-01 + ftmargoh_(range_normalised)_____________________________________________ (nitvar020)___________________ 9.9868E-01 + ftaulimit_______________________________________________________________ (itvar021)____________________ 7.1296E-01 + ftaulimit_(final_value/initial_value)___________________________________ (xcm021)______________________ 7.1296E-01 + ftaulimit_(range_normalised)____________________________________________ (nitvar021)___________________ 7.1268E-01 + fmaxvvstress____________________________________________________________ (itvar022)____________________ 4.3859E-01 + fmaxvvstress_(final_value/initial_value)________________________________ (xcm022)______________________ 4.3859E-01 + fmaxvvstress_(range_normalised)_________________________________________ (nitvar022)___________________ 4.3803E-01 + foh_stress______________________________________________________________ (itvar023)____________________ 9.2954E-01 + foh_stress_(final_value/initial_value)__________________________________ (xcm023)______________________ 9.2954E-01 + foh_stress_(range_normalised)___________________________________________ (nitvar023)___________________ 9.2947E-01 fne0____________________________________________________________________ (itvar024)____________________ 5.9638E-01 fne0_(final_value/initial_value)________________________________________ (xcm024)______________________ 5.9638E-01 fne0_(range_normalised)_________________________________________________ (nitvar024)___________________ 5.9598E-01 - fpsepbqar_______________________________________________________________ (itvar025)____________________ 9.5570E-01 - fpsepbqar_(final_value/initial_value)___________________________________ (xcm025)______________________ 9.5570E-01 - fpsepbqar_(range_normalised)____________________________________________ (nitvar025)___________________ 9.5566E-01 - fstrcase________________________________________________________________ (itvar026)____________________ 1.0000E+00 - fstrcase_(final_value/initial_value)____________________________________ (xcm026)______________________ 1.0000E+00 - fstrcase_(range_normalised)_____________________________________________ (nitvar026)___________________ 1.0000E+00 - fstrcond________________________________________________________________ (itvar027)____________________ 8.4877E-01 - fstrcond_(final_value/initial_value)____________________________________ (xcm027)______________________ 8.4877E-01 - fstrcond_(range_normalised)_____________________________________________ (nitvar027)___________________ 8.4862E-01 - bt______________________________________________________________________ (itvar028)____________________ 5.0627E+00 - bt_(final_value/initial_value)__________________________________________ (xcm028)______________________ 8.8820E-01 - bt_(range_normalised)___________________________________________________ (nitvar028)___________________ 1.6848E-01 + fpsepbqar_______________________________________________________________ (itvar025)____________________ 9.6274E-01 + fpsepbqar_(final_value/initial_value)___________________________________ (xcm025)______________________ 9.6274E-01 + fpsepbqar_(range_normalised)____________________________________________ (nitvar025)___________________ 9.6271E-01 + fstrcase________________________________________________________________ (itvar026)____________________ 9.9838E-01 + fstrcase_(final_value/initial_value)____________________________________ (xcm026)______________________ 9.9838E-01 + fstrcase_(range_normalised)_____________________________________________ (nitvar026)___________________ 9.9838E-01 + fstrcond________________________________________________________________ (itvar027)____________________ 8.4907E-01 + fstrcond_(final_value/initial_value)____________________________________ (xcm027)______________________ 8.4907E-01 + fstrcond_(range_normalised)_____________________________________________ (nitvar027)___________________ 8.4892E-01 + bt______________________________________________________________________ (itvar028)____________________ 5.2376E+00 + bt_(final_value/initial_value)__________________________________________ (xcm028)______________________ 9.1887E-01 + bt_(range_normalised)___________________________________________________ (nitvar028)___________________ 1.7431E-01 rmajor__________________________________________________________________ (itvar029)____________________ 8.0000E+00 rmajor_(final_value/initial_value)______________________________________ (xcm029)______________________ 1.0000E+00 rmajor_(range_normalised)_______________________________________________ (nitvar029)___________________ 0.0000E+00 - te______________________________________________________________________ (itvar030)____________________ 1.2238E+01 - te_(final_value/initial_value)__________________________________________ (xcm030)______________________ 1.0198E+00 - te_(range_normalised)___________________________________________________ (nitvar030)___________________ 7.1249E-02 - hfact___________________________________________________________________ (itvar031)____________________ 1.1808E+00 - hfact_(final_value/initial_value)_______________________________________ (xcm031)______________________ 1.0734E+00 - hfact_(range_normalised)________________________________________________ (nitvar031)___________________ 9.8250E-01 + te______________________________________________________________________ (itvar030)____________________ 1.2437E+01 + te_(final_value/initial_value)__________________________________________ (xcm030)______________________ 1.0364E+00 + te_(range_normalised)___________________________________________________ (nitvar030)___________________ 7.3355E-02 + hfact___________________________________________________________________ (itvar031)____________________ 1.1947E+00 + hfact_(final_value/initial_value)_______________________________________ (xcm031)______________________ 1.0861E+00 + hfact_(range_normalised)________________________________________________ (nitvar031)___________________ 9.9517E-01 tfcth___________________________________________________________________ (itvar032)____________________ 1.2000E+00 tfcth_(final_value/initial_value)_______________________________________ (xcm032)______________________ 1.0000E+00 tfcth_(range_normalised)________________________________________________ (nitvar032)___________________ 1.1628E-01 - ohcth___________________________________________________________________ (itvar033)____________________ 5.5551E-01 - ohcth_(final_value/initial_value)_______________________________________ (xcm033)______________________ 1.1110E+00 - ohcth_(range_normalised)________________________________________________ (nitvar033)___________________ 2.6342E-02 - q_______________________________________________________________________ (itvar034)____________________ 3.4251E+00 - q_(final_value/initial_value)___________________________________________ (xcm034)______________________ 9.7860E-01 - q_(range_normalised)____________________________________________________ (nitvar034)___________________ 9.0450E-03 - bore____________________________________________________________________ (itvar035)____________________ 2.0588E+00 - bore_(final_value/initial_value)________________________________________ (xcm035)______________________ 1.0294E+00 - bore_(range_normalised)_________________________________________________ (nitvar035)___________________ 1.9786E-01 - fvsbrnni________________________________________________________________ (itvar036)____________________ 4.1495E-01 - fvsbrnni_(final_value/initial_value)____________________________________ (xcm036)______________________ 1.0374E+00 - fvsbrnni_(range_normalised)_____________________________________________ (nitvar036)___________________ 4.1437E-01 - tdmptf__________________________________________________________________ (itvar037)____________________ 1.9607E+01 - tdmptf_(final_value/initial_value)______________________________________ (xcm037)______________________ 7.8427E-01 - tdmptf_(range_normalised)_______________________________________________ (nitvar037)___________________ 1.9526E-01 - thkcas__________________________________________________________________ (itvar038)____________________ 2.2965E-01 - thkcas_(final_value/initial_value)______________________________________ (xcm038)______________________ 4.5931E-01 - thkcas_(range_normalised)_______________________________________________ (nitvar038)___________________ 1.8911E-01 - thwcndut________________________________________________________________ (itvar039)____________________ 8.0000E-03 - thwcndut_(final_value/initial_value)____________________________________ (xcm039)______________________ 1.0000E+00 - thwcndut_(range_normalised)_____________________________________________ (nitvar039)___________________ 0.0000E+00 - fcutfsu_________________________________________________________________ (itvar040)____________________ 8.6211E-01 - fcutfsu_(final_value/initial_value)_____________________________________ (xcm040)______________________ 1.0776E+00 - fcutfsu_(range_normalised)______________________________________________ (nitvar040)___________________ 8.2299E-01 - cpttf___________________________________________________________________ (itvar041)____________________ 8.0573E+04 - cpttf_(final_value/initial_value)_______________________________________ (xcm041)______________________ 1.2396E+00 - cpttf_(range_normalised)________________________________________________ (nitvar041)___________________ 6.2291E-01 - ralpne__________________________________________________________________ (itvar042)____________________ 8.2363E-02 - ralpne_(final_value/initial_value)______________________________________ (xcm042)______________________ 8.2363E-01 - ralpne_(range_normalised)_______________________________________________ (nitvar042)___________________ 6.4725E-01 - oh_steel_frac___________________________________________________________ (itvar043)____________________ 5.3310E-01 - oh_steel_frac_(final_value/initial_value)_______________________________ (xcm043)______________________ 6.6637E-01 - oh_steel_frac_(range_normalised)________________________________________ (nitvar043)___________________ 5.6069E-01 - fimp(13)________________________________________________________________ (itvar044)____________________ 5.6323E-04 - fimp(13)_(final_value/initial_value)____________________________________ (xcm044)______________________ 1.4822E+00 - fimp(13)_(range_normalised)_____________________________________________ (nitvar044)___________________ 5.6322E-02 - dr_tf_wp________________________________________________________________ (itvar045)____________________ 5.0550E-01 - dr_tf_wp_(final_value/initial_value)____________________________________ (xcm045)______________________ 1.0110E+00 - dr_tf_wp_(range_normalised)_____________________________________________ (nitvar045)___________________ 6.5940E-02 - Beta_consistency__________________normalised_residue____________________ (normres001)__________________ -3.7854E-12 - Global_power_balance_consistency__normalised_residue____________________ (normres002)__________________ 1.6275E-09 - Radial_build_consistency__________normalised_residue____________________ (normres003)__________________ -6.6613E-16 - Density_upper_limit_______________normalised_residue____________________ (normres004)__________________ 1.8163E-13 - Neutron_wall_load_upper_limit_____normalised_residue____________________ (normres005)__________________ 3.2079E-11 - Fusion_power_upper_limit__________normalised_residue____________________ (normres006)__________________ 3.2079E-11 - Burn_time_lower_limit_____________normalised_residue____________________ (normres007)__________________ -1.1697E-05 - L-H_power_threshold_limit_________normalised_residue____________________ (normres008)__________________ -1.2914E-08 - Injection_power_upper_limit_______normalised_residue____________________ (normres009)__________________ -2.6298E-08 - Net_electric_power_lower_limit____normalised_residue____________________ (normres010)__________________ 1.4277E-07 - Beta_upper_limit__________________normalised_residue____________________ (normres011)__________________ -4.7070E-12 - Peak_toroidal_field_upper_limit___normalised_residue____________________ (normres012)__________________ 8.3267E-14 - CS_coil_EOF_current_density_limit_normalised_residue____________________ (normres013)__________________ -9.6190E-13 - CS_coil_BOP_current_density_limit_normalised_residue____________________ (normres014)__________________ -1.1113E-12 - I_op_/_I_critical_(TF_coil)_______normalised_residue____________________ (normres015)__________________ 6.9611E-13 - Dump_voltage_upper_limit__________normalised_residue____________________ (normres016)__________________ -3.0620E-13 - J_winding_pack/J_protection_limit_normalised_residue____________________ (normres017)__________________ 2.0761E-13 - TF_coil_temp._margin_lower_limit__normalised_residue____________________ (normres018)__________________ -1.2975E-12 - CS_temperature_margin_lower_limit_normalised_residue____________________ (normres019)__________________ -0.0000E+00 - taup/taueff_______________________normalised_residue____________________ (normres020)__________________ -4.7055E-09 - Dump_time_set_by_VV_stress________normalised_residue____________________ (normres021)__________________ 3.0487E-13 - CS_Tresca_yield_criterion_________normalised_residue____________________ (normres022)__________________ -4.9472E-13 - ne0_>_neped_______________________normalised_residue____________________ (normres023)__________________ 1.2901E-13 - Upper_Lim._on_Psep_*_Bt_/_q_A_R___normalised_residue____________________ (normres024)__________________ 1.2885E-08 - TF_coil_case_stress_upper_limit___normalised_residue____________________ (normres025)__________________ 2.7156E-13 - TF_coil_conduit_stress_upper_lim__normalised_residue____________________ (normres026)__________________ 2.6334E-13 + ohcth___________________________________________________________________ (itvar033)____________________ 5.5625E-01 + ohcth_(final_value/initial_value)_______________________________________ (xcm033)______________________ 1.1125E+00 + ohcth_(range_normalised)________________________________________________ (nitvar033)___________________ 2.6417E-02 + q_______________________________________________________________________ (itvar034)____________________ 3.5610E+00 + q_(final_value/initial_value)___________________________________________ (xcm034)______________________ 1.0174E+00 + q_(range_normalised)____________________________________________________ (nitvar034)___________________ 1.1936E-02 + bore____________________________________________________________________ (itvar035)____________________ 2.0251E+00 + bore_(final_value/initial_value)________________________________________ (xcm035)______________________ 1.0126E+00 + bore_(range_normalised)_________________________________________________ (nitvar035)___________________ 1.9446E-01 + fvsbrnni________________________________________________________________ (itvar036)____________________ 4.1997E-01 + fvsbrnni_(final_value/initial_value)____________________________________ (xcm036)______________________ 1.0499E+00 + fvsbrnni_(range_normalised)_____________________________________________ (nitvar036)___________________ 4.1938E-01 + tdmptf__________________________________________________________________ (itvar037)____________________ 1.7534E+01 + tdmptf_(final_value/initial_value)______________________________________ (xcm037)______________________ 7.0135E-01 + tdmptf_(range_normalised)_______________________________________________ (nitvar037)___________________ 1.7451E-01 + thkcas__________________________________________________________________ (itvar038)____________________ 2.6741E-01 + thkcas_(final_value/initial_value)______________________________________ (xcm038)______________________ 5.3481E-01 + thkcas_(range_normalised)_______________________________________________ (nitvar038)___________________ 2.2885E-01 + thwcndut________________________________________________________________ (itvar039)____________________ 8.0091E-03 + thwcndut_(final_value/initial_value)____________________________________ (xcm039)______________________ 1.0011E+00 + thwcndut_(range_normalised)_____________________________________________ (nitvar039)___________________ 9.9068E-05 + fcutfsu_________________________________________________________________ (itvar040)____________________ 8.4176E-01 + fcutfsu_(final_value/initial_value)_____________________________________ (xcm040)______________________ 1.0522E+00 + fcutfsu_(range_normalised)______________________________________________ (nitvar040)___________________ 7.7673E-01 + cpttf___________________________________________________________________ (itvar041)____________________ 8.4815E+04 + cpttf_(final_value/initial_value)_______________________________________ (xcm041)______________________ 1.3048E+00 + cpttf_(range_normalised)________________________________________________ (nitvar041)___________________ 7.9260E-01 + ralpne__________________________________________________________________ (itvar042)____________________ 8.5281E-02 + ralpne_(final_value/initial_value)______________________________________ (xcm042)______________________ 8.5281E-01 + ralpne_(range_normalised)_______________________________________________ (nitvar042)___________________ 7.0562E-01 + oh_steel_frac___________________________________________________________ (itvar043)____________________ 5.2166E-01 + oh_steel_frac_(final_value/initial_value)_______________________________ (xcm043)______________________ 6.5207E-01 + oh_steel_frac_(range_normalised)________________________________________ (nitvar043)___________________ 5.4864E-01 + fimp(13)________________________________________________________________ (itvar044)____________________ 5.6341E-04 + fimp(13)_(final_value/initial_value)____________________________________ (xcm044)______________________ 1.4826E+00 + fimp(13)_(range_normalised)_____________________________________________ (nitvar044)___________________ 5.6340E-02 + dr_tf_wp________________________________________________________________ (itvar045)____________________ 4.9973E-01 + dr_tf_wp_(final_value/initial_value)____________________________________ (xcm045)______________________ 9.9945E-01 + dr_tf_wp_(range_normalised)_____________________________________________ (nitvar045)___________________ 6.2329E-02 + Beta_consistency__________________normalised_residue____________________ (normres001)__________________ 1.5917E-08 + Global_power_balance_consistency__normalised_residue____________________ (normres002)__________________ 1.3192E-07 + Radial_build_consistency__________normalised_residue____________________ (normres003)__________________ 4.1783E-10 + Density_upper_limit_______________normalised_residue____________________ (normres004)__________________ 3.3764E-08 + Neutron_wall_load_upper_limit_____normalised_residue____________________ (normres005)__________________ 3.0867E-07 + Fusion_power_upper_limit__________normalised_residue____________________ (normres006)__________________ 3.2827E-07 + Burn_time_lower_limit_____________normalised_residue____________________ (normres007)__________________ -6.4872E-05 + L-H_power_threshold_limit_________normalised_residue____________________ (normres008)__________________ -7.3515E-07 + Injection_power_upper_limit_______normalised_residue____________________ (normres009)__________________ 8.9159E-08 + Net_electric_power_lower_limit____normalised_residue____________________ (normres010)__________________ 8.7522E-06 + Beta_upper_limit__________________normalised_residue____________________ (normres011)__________________ -2.5620E-08 + Peak_toroidal_field_upper_limit___normalised_residue____________________ (normres012)__________________ -1.7419E-08 + CS_coil_EOF_current_density_limit_normalised_residue____________________ (normres013)__________________ -3.4768E-08 + CS_coil_BOP_current_density_limit_normalised_residue____________________ (normres014)__________________ -3.1882E-08 + I_op_/_I_critical_(TF_coil)_______normalised_residue____________________ (normres015)__________________ -4.4579E-07 + Dump_voltage_upper_limit__________normalised_residue____________________ (normres016)__________________ -5.2052E-08 + J_winding_pack/J_protection_limit_normalised_residue____________________ (normres017)__________________ -6.3741E-08 + TF_coil_temp._margin_lower_limit__normalised_residue____________________ (normres018)__________________ -1.3486E-06 + CS_temperature_margin_lower_limit_normalised_residue____________________ (normres019)__________________ -2.6193E-08 + taup/taueff_______________________normalised_residue____________________ (normres020)__________________ 1.5765E-07 + Dump_time_set_by_VV_stress________normalised_residue____________________ (normres021)__________________ -1.3625E-07 + CS_Tresca_yield_criterion_________normalised_residue____________________ (normres022)__________________ -4.4229E-08 + ne0_>_neped_______________________normalised_residue____________________ (normres023)__________________ -4.3315E-08 + Upper_Lim._on_Psep_*_Bt_/_q_A_R___normalised_residue____________________ (normres024)__________________ 7.3919E-07 + TF_coil_case_stress_upper_limit___normalised_residue____________________ (normres025)__________________ -1.2835E-07 + TF_coil_conduit_stress_upper_lim__normalised_residue____________________ (normres026)__________________ -1.1371E-07 # Final Feasible Point # # Power Reactor Costs (1990 US$) # - First_wall_/_blanket_life_(years)_______________________________________ (fwbllife)____________________ 6.2630E+00 - Divertor_life_(years)___________________________________________________ (divlife.)____________________ 4.7579E+00 - Cost_of_electricity_(m$/kWh)____________________________________________ (coe)_________________________ 5.0303E+02 + First_wall_/_blanket_life_(years)_______________________________________ (fwbllife)____________________ 6.2457E+00 + Divertor_life_(years)___________________________________________________ (divlife.)____________________ 4.4952E+00 + Cost_of_electricity_(m$/kWh)____________________________________________ (coe)_________________________ 5.0802E+02 # Detailed Costings (1990 US$) # Acc.22_multiplier_for_Nth_of_a_kind_____________________________________ (fkind)_______________________ 1.0000E+00 Level_of_Safety_Assurance_______________________________________________ (lsa)_________________________ 4 # Structures and Site Facilities # Site_improvements,_facilities,_land_(M$)________________________________ (c211)________________________ 3.5200E+01 - Reactor_building_cost_(M$)______________________________________________ (c212)________________________ 4.4905E+02 + Reactor_building_cost_(M$)______________________________________________ (c212)________________________ 4.5113E+02 Turbine_building_cost_(M$)______________________________________________ (c213)________________________ 3.8000E+01 Reactor_maintenance_building_cost_(M$)__________________________________ (c2141)_______________________ 1.0307E+02 Warm_shop_cost_(M$)_____________________________________________________ (c2142)_______________________ 5.7641E+01 Tritium_building_cost_(M$)______________________________________________ (c215)________________________ 1.4800E+01 - Electrical_equipment_building_cost_(M$)_________________________________ (c216)________________________ 1.9338E+01 + Electrical_equipment_building_cost_(M$)_________________________________ (c216)________________________ 1.9722E+01 Additional_buildings_cost_(M$)__________________________________________ (c2171)_______________________ 1.8000E+01 Control_room_buildings_cost_(M$)________________________________________ (c2172)_______________________ 2.1000E+01 Shop_and_warehouses_cost_(M$)___________________________________________ (c2173)_______________________ 1.1500E+01 - Cryogenic_building_cost_(M$)____________________________________________ (c2174)_______________________ 6.7114E+00 - Total_account_21_cost_(M$)______________________________________________ (c21)_________________________ 7.7431E+02 + Cryogenic_building_cost_(M$)____________________________________________ (c2174)_______________________ 6.7990E+00 + Total_account_21_cost_(M$)______________________________________________ (c21)_________________________ 7.7686E+02 # Reactor Systems # First_wall_cost_(M$)____________________________________________________ (c2211)_______________________ 1.9130E+02 Blanket_beryllium_cost_(M$)_____________________________________________ (c22121)______________________ 2.6109E+02 @@ -214,109 +214,109 @@ Bulk_shield_cost_(M$)___________________________________________________ (c22131)______________________ 7.8394E+01 Penetration_shielding_cost_(M$)_________________________________________ (c22132)______________________ 7.8394E+01 Total_shield_cost_(M$)__________________________________________________ (c2213)_______________________ 1.5679E+02 - Total_support_structure_cost_(M$)_______________________________________ (c2214)_______________________ 3.9599E+01 + Total_support_structure_cost_(M$)_______________________________________ (c2214)_______________________ 4.0485E+01 Divertor_cost_(M$)______________________________________________________ (c2215)_______________________ 4.1660E+01 - Total_account_221_cost_(M$)_____________________________________________ (c221)________________________ 1.4110E+03 + Total_account_221_cost_(M$)_____________________________________________ (c221)________________________ 1.4119E+03 # Magnets # - TF_coil_conductor_cost_(M$)_____________________________________________ (c22211)______________________ 1.1908E+02 - TF_coil_winding_cost_(M$)_______________________________________________ (c22212)______________________ 5.7685E+01 - TF_coil_case_cost_(M$)__________________________________________________ (c22213)______________________ 3.9778E+02 - TF_intercoil_structure_cost_(M$)________________________________________ (c22214)______________________ 1.4196E+02 - TF_coil_gravity_support_structure_(M$)__________________________________ (c22215)______________________ 4.3172E+01 - TF_magnet_assemblies_cost_(M$)__________________________________________ (c2221)_______________________ 7.5968E+02 - PF_coil_conductor_cost_(M$)_____________________________________________ (c22221)______________________ 4.6338E+02 - PF_coil_winding_cost_(M$)_______________________________________________ (c22222)______________________ 7.6727E+01 - PF_coil_case_cost_(M$)__________________________________________________ (c22223)______________________ 9.6937E+01 - PF_coil_support_structure_cost_(M$)_____________________________________ (c22224)______________________ 8.0889E+00 - PF_magnet_assemblies_cost_(M$)__________________________________________ (c2222)_______________________ 6.4514E+02 + TF_coil_conductor_cost_(M$)_____________________________________________ (c22211)______________________ 1.2631E+02 + TF_coil_winding_cost_(M$)_______________________________________________ (c22212)______________________ 5.6782E+01 + TF_coil_case_cost_(M$)__________________________________________________ (c22213)______________________ 4.2976E+02 + TF_intercoil_structure_cost_(M$)________________________________________ (c22214)______________________ 1.4787E+02 + TF_coil_gravity_support_structure_(M$)__________________________________ (c22215)______________________ 4.4157E+01 + TF_magnet_assemblies_cost_(M$)__________________________________________ (c2221)_______________________ 8.0487E+02 + PF_coil_conductor_cost_(M$)_____________________________________________ (c22221)______________________ 4.6310E+02 + PF_coil_winding_cost_(M$)_______________________________________________ (c22222)______________________ 7.6566E+01 + PF_coil_case_cost_(M$)__________________________________________________ (c22223)______________________ 9.5238E+01 + PF_coil_support_structure_cost_(M$)_____________________________________ (c22224)______________________ 8.0091E+00 + PF_magnet_assemblies_cost_(M$)__________________________________________ (c2222)_______________________ 6.4292E+02 Vacuum_vessel_assembly_cost_(M$)________________________________________ (c2223)_______________________ 2.5404E+02 - Total_account_222_cost_(M$)_____________________________________________ (c222)________________________ 1.6589E+03 + Total_account_222_cost_(M$)_____________________________________________ (c222)________________________ 1.7018E+03 # Power Injection # - ECH_system_cost_(M$)____________________________________________________ (c2231)_______________________ 2.2564E+02 + ECH_system_cost_(M$)____________________________________________________ (c2231)_______________________ 2.2619E+02 Lower_hybrid_system_cost_(M$)___________________________________________ (c2232)_______________________ 0.0000E+00 Neutral_beam_system_cost_(M$)___________________________________________ (c2233)_______________________ 0.0000E+00 - Total_account_223_cost_(M$)_____________________________________________ (c223)________________________ 2.2564E+02 + Total_account_223_cost_(M$)_____________________________________________ (c223)________________________ 2.2619E+02 # Vacuum Systems # High_vacuum_pumps_cost_(M$)_____________________________________________ (c2241)_______________________ 1.2480E+01 Backing_pumps_cost_(M$)_________________________________________________ (c2242)_______________________ 4.6800E+00 - Vacuum_duct_cost_(M$)___________________________________________________ (c2243)_______________________ 2.5133E+00 - Valves_cost_(M$)________________________________________________________ (c2244)_______________________ 5.7355E+00 + Vacuum_duct_cost_(M$)___________________________________________________ (c2243)_______________________ 2.4485E+00 + Valves_cost_(M$)________________________________________________________ (c2244)_______________________ 5.5707E+00 Duct_shielding_cost_(M$)________________________________________________ (c2245)_______________________ 0.0000E+00 Instrumentation_cost_(M$)_______________________________________________ (c2246)_______________________ 1.3000E+00 - Total_account_224_cost_(M$)_____________________________________________ (c224)________________________ 2.6709E+01 + Total_account_224_cost_(M$)_____________________________________________ (c224)________________________ 2.6479E+01 # Power Conditioning # - TF_coil_power_supplies_cost_(M$)________________________________________ (c22511)______________________ 3.7398E+00 - TF_coil_breakers_cost_(M$)______________________________________________ (c22512)______________________ 3.0527E+01 - TF_coil_dump_resistors_cost_(M$)________________________________________ (c22513)______________________ 1.9349E+01 + TF_coil_power_supplies_cost_(M$)________________________________________ (c22511)______________________ 3.9137E+00 + TF_coil_breakers_cost_(M$)______________________________________________ (c22512)______________________ 3.4719E+01 + TF_coil_dump_resistors_cost_(M$)________________________________________ (c22513)______________________ 2.0798E+01 TF_coil_instrumentation_and_control_(M$)________________________________ (c22514)______________________ 4.8000E+00 - TF_coil_bussing_cost_(M$)_______________________________________________ (c22515)______________________ 3.1645E+01 - Total,_TF_coil_power_costs_(M$)_________________________________________ (c2251)_______________________ 9.0061E+01 - PF_coil_power_supplies_cost_(M$)________________________________________ (c22521)______________________ 2.6080E+00 + TF_coil_bussing_cost_(M$)_______________________________________________ (c22515)______________________ 3.3311E+01 + Total,_TF_coil_power_costs_(M$)_________________________________________ (c2251)_______________________ 9.7542E+01 + PF_coil_power_supplies_cost_(M$)________________________________________ (c22521)______________________ 2.7897E+00 PF_coil_instrumentation_and_control_(M$)________________________________ (c22522)______________________ 3.6000E+00 PF_coil_bussing_cost_(M$)_______________________________________________ (c22523)______________________ 1.1995E+01 - PF_coil_burn_power_supplies_cost_(M$)___________________________________ (c22524)______________________ 1.2675E+00 + PF_coil_burn_power_supplies_cost_(M$)___________________________________ (c22524)______________________ 1.2676E+00 PF_coil_breakers_cost_(M$)______________________________________________ (c22525)______________________ 1.4709E+01 - PF_coil_dump_resistors_cost_(M$)________________________________________ (c22526)______________________ 4.2001E+00 + PF_coil_dump_resistors_cost_(M$)________________________________________ (c22526)______________________ 4.2348E+00 PF_coil_ac_breakers_cost_(M$)___________________________________________ (c22527)______________________ 9.0000E-01 - Total,_PF_coil_power_costs_(M$)_________________________________________ (c2252)_______________________ 3.9280E+01 + Total,_PF_coil_power_costs_(M$)_________________________________________ (c2252)_______________________ 3.9497E+01 Total,_energy_storage_cost_(M$)_________________________________________ (c2253)_______________________ 1.6864E+01 - Total_account_225_cost_(M$)_____________________________________________ (c225)________________________ 1.4621E+02 + Total_account_225_cost_(M$)_____________________________________________ (c225)________________________ 1.5390E+02 # Heat Transport System # - Pumps_and_piping_system_cost_(M$)_______________________________________ (cpp)_________________________ 3.6356E+01 - Primary_heat_exchanger_cost_(M$)________________________________________ (chx)_________________________ 7.0758E+01 - Total,_reactor_cooling_system_cost_(M$)_________________________________ (c2261)_______________________ 1.0711E+02 - Pumps,_piping_cost_(M$)_________________________________________________ (cppa)________________________ 2.7713E+01 - Total,_auxiliary_cooling_system_cost_(M$)_______________________________ (c2262)_______________________ 2.7713E+01 - Total,_cryogenic_system_cost_(M$)_______________________________________ (c2263)_______________________ 1.6452E+02 - Total_account_226_cost_(M$)_____________________________________________ (c226)________________________ 2.9935E+02 + Pumps_and_piping_system_cost_(M$)_______________________________________ (cpp)_________________________ 3.6427E+01 + Primary_heat_exchanger_cost_(M$)________________________________________ (chx)_________________________ 7.0891E+01 + Total,_reactor_cooling_system_cost_(M$)_________________________________ (c2261)_______________________ 1.0732E+02 + Pumps,_piping_cost_(M$)_________________________________________________ (cppa)________________________ 2.7855E+01 + Total,_auxiliary_cooling_system_cost_(M$)_______________________________ (c2262)_______________________ 2.7855E+01 + Total,_cryogenic_system_cost_(M$)_______________________________________ (c2263)_______________________ 1.6741E+02 + Total_account_226_cost_(M$)_____________________________________________ (c226)________________________ 3.0258E+02 # Fuel Handling System # Fuelling_system_cost_(M$)_______________________________________________ (c2271)_______________________ 2.2300E+01 - Fuel_processing_and_purification_cost_(M$)______________________________ (c2272)_______________________ 1.0137E+02 - Atmospheric_recovery_systems_cost_(M$)__________________________________ (c2273)_______________________ 5.7691E+01 - Nuclear_building_ventilation_cost_(M$)__________________________________ (c2274)_______________________ 6.8823E+01 - Total_account_227_cost_(M$)_____________________________________________ (c227)________________________ 2.5019E+02 + Fuel_processing_and_purification_cost_(M$)______________________________ (c2272)_______________________ 1.0152E+02 + Atmospheric_recovery_systems_cost_(M$)__________________________________ (c2273)_______________________ 5.7939E+01 + Nuclear_building_ventilation_cost_(M$)__________________________________ (c2274)_______________________ 6.9059E+01 + Total_account_227_cost_(M$)_____________________________________________ (c227)________________________ 2.5082E+02 # Instrumentation and Control # Instrumentation_and_control_cost_(M$)___________________________________ (c228)________________________ 1.5000E+02 # Maintenance Equipment # Maintenance_equipment_cost_(M$)_________________________________________ (c229)________________________ 1.2500E+02 # Total Account 22 Cost # - Total_account_22_cost_(M$)______________________________________________ (c22)_________________________ 4.2930E+03 + Total_account_22_cost_(M$)______________________________________________ (c22)_________________________ 4.3487E+03 # Turbine Plant Equipment # - Turbine_plant_equipment_cost_(M$)_______________________________________ (c23)_________________________ 1.7341E+02 + Turbine_plant_equipment_cost_(M$)_______________________________________ (c23)_________________________ 1.7380E+02 # Electric Plant Equipment # Switchyard_equipment_cost_(M$)__________________________________________ (c241)________________________ 1.8400E+01 - Transformers_cost_(M$)__________________________________________________ (c242)________________________ 7.7268E+00 - Low_voltage_equipment_cost_(M$)_________________________________________ (c243)________________________ 6.2652E+00 + Transformers_cost_(M$)__________________________________________________ (c242)________________________ 7.8162E+00 + Low_voltage_equipment_cost_(M$)_________________________________________ (c243)________________________ 6.3369E+00 Diesel_backup_equipment_cost_(M$)_______________________________________ (c244)________________________ 6.8000E+00 Auxiliary_facilities_cost_(M$)__________________________________________ (c245)________________________ 1.5000E+00 - Total_account_24_cost_(M$)______________________________________________ (c24)_________________________ 4.0692E+01 + Total_account_24_cost_(M$)______________________________________________ (c24)_________________________ 4.0853E+01 # Miscellaneous Plant Equipment # Miscellaneous_plant_equipment_cost_(M$)_________________________________ (c25)_________________________ 2.5000E+01 # Heat Rejection System # - Heat_rejection_system_cost_(M$)_________________________________________ (c26)_________________________ 4.8951E+01 + Heat_rejection_system_cost_(M$)_________________________________________ (c26)_________________________ 4.9083E+01 # Plant Direct Cost # - Plant_direct_cost_(M$)__________________________________________________ (cdirt)_______________________ 5.3553E+03 + Plant_direct_cost_(M$)__________________________________________________ (cdirt)_______________________ 5.4143E+03 # Reactor Core Cost # - Reactor_core_cost_(M$)__________________________________________________ (crctcore)____________________ 3.2955E+03 + Reactor_core_cost_(M$)__________________________________________________ (crctcore)____________________ 3.3399E+03 # Indirect Cost # - Indirect_cost_(M$)______________________________________________________ (c9)__________________________ 1.7860E+03 + Indirect_cost_(M$)______________________________________________________ (c9)__________________________ 1.8057E+03 # Total Contingency # - Total_contingency_(M$)__________________________________________________ (ccont)_______________________ 1.3926E+03 + Total_contingency_(M$)__________________________________________________ (ccont)_______________________ 1.4079E+03 # Constructed Cost # - Constructed_cost_(M$)___________________________________________________ (concost)_____________________ 8.5339E+03 + Constructed_cost_(M$)___________________________________________________ (concost)_____________________ 8.6279E+03 # Interest during Construction # - Interest_during_construction_(M$)_______________________________________ (moneyint)____________________ 1.4081E+03 + Interest_during_construction_(M$)_______________________________________ (moneyint)____________________ 1.4236E+03 # Total Capital Investment # - Total_capital_investment_(M$)___________________________________________ (capcost)_____________________ 9.9420E+03 + Total_capital_investment_(M$)___________________________________________ (capcost)_____________________ 1.0051E+04 # Plant Availability # Allowable_blanket_neutron_fluence_(MW-yr/m2)____________________________ (abktflnc)____________________ 5.0000E+00 Allowable_divertor_heat_fluence_(MW-yr/m2)______________________________ (adivflnc)____________________ 7.0000E+00 - First_wall_/_blanket_lifetime_(years)___________________________________ (bktlife)_____________________ 6.2630E+00 OP - Divertor_lifetime_(years)_______________________________________________ (divlife)_____________________ 4.7579E+00 OP - Heating/CD_system_lifetime_(years)______________________________________ (cdrlife)_____________________ 6.2630E+00 OP + First_wall_/_blanket_lifetime_(years)___________________________________ (bktlife)_____________________ 6.2457E+00 OP + Divertor_lifetime_(years)_______________________________________________ (divlife)_____________________ 4.4952E+00 OP + Heating/CD_system_lifetime_(years)______________________________________ (cdrlife)_____________________ 6.2457E+00 OP Total_plant_lifetime_(years)____________________________________________ (tlife)_______________________ 3.0000E+01 Total_plant_availability_fraction_______________________________________ (cfactr)______________________ 8.0000E-01 - Number_of_fusion_cycles_to_reach_allowable_fw/blanket_DPA_______________ (bktcycles)___________________ 1.6063E+04 + Number_of_fusion_cycles_to_reach_allowable_fw/blanket_DPA_______________ (bktcycles)___________________ 1.6004E+04 # Plasma # Tokamak_aspect_ratio_=_Conventional,_itart_=_0__________________________ (itart)_______________________ 0.0000E+00 Major_radius_(m)________________________________________________________ (rmajor)______________________ 8.0000E+00 ITV @@ -332,54 +332,54 @@ Plasma_surface_area_(m2)________________________________________________ (sarea)_______________________ 1.1738E+03 OP Plasma_volume_(m3)______________________________________________________ (vol)_________________________ 1.8882E+03 OP Plasma_current_scaling_law_used_________________________________________ (icurr)_______________________ 4 - Plasma_current_(MA)_____________________________________________________ (plascur/1D6)_________________ 1.6699E+01 - Current_density_profile_factor__________________________________________ (alphaj)______________________ 1.8387E+00 - Plasma_internal_inductance,_li__________________________________________ (rli)_________________________ 1.1898E+00 - Vertical_field_at_plasma_(T)____________________________________________ (bvert)_______________________ -7.3601E-01 - Vacuum_toroidal_field_at_R_(T)__________________________________________ (bt)__________________________ 5.0627E+00 ITV - Average_poloidal_field_(T)______________________________________________ (bp)__________________________ 8.7139E-01 - Total_field_(sqrt(bp^2_+_bt^2))_(T)_____________________________________ (btot)________________________ 5.1372E+00 + Plasma_current_(MA)_____________________________________________________ (plascur/1D6)_________________ 1.6616E+01 + Current_density_profile_factor__________________________________________ (alphaj)______________________ 1.9514E+00 + Plasma_internal_inductance,_li__________________________________________ (rli)_________________________ 1.2199E+00 + Vertical_field_at_plasma_(T)____________________________________________ (bvert)_______________________ -7.4089E-01 + Vacuum_toroidal_field_at_R_(T)__________________________________________ (bt)__________________________ 5.2376E+00 ITV + Average_poloidal_field_(T)______________________________________________ (bp)__________________________ 8.6708E-01 + Total_field_(sqrt(bp^2_+_bt^2))_(T)_____________________________________ (btot)________________________ 5.3089E+00 Safety_factor_on_axis___________________________________________________ (q0)__________________________ 1.0000E+00 - Safety_factor_at_95%_flux_surface_______________________________________ (q95)_________________________ 3.4251E+00 ITV - Cylindrical_safety_factor_(qcyl)________________________________________ (qstar)_______________________ 2.8387E+00 - Total_plasma_beta_______________________________________________________ (beta)________________________ 3.6055E-02 ITV - Total_poloidal_beta_____________________________________________________ (betap)_______________________ 1.2531E+00 OP - Total_toroidal_beta_____________________________________________________ ______________________________ 3.7123E-02 OP - Fast_alpha_beta_________________________________________________________ (betaft)______________________ 4.4814E-03 OP + Safety_factor_at_95%_flux_surface_______________________________________ (q95)_________________________ 3.5610E+00 ITV + Cylindrical_safety_factor_(qcyl)________________________________________ (qstar)_______________________ 2.9514E+00 + Total_plasma_beta_______________________________________________________ (beta)________________________ 3.4120E-02 ITV + Total_poloidal_beta_____________________________________________________ (betap)_______________________ 1.2791E+00 OP + Total_toroidal_beta_____________________________________________________ ______________________________ 3.5056E-02 OP + Fast_alpha_beta_________________________________________________________ (betaft)______________________ 4.2465E-03 OP Beam_ion_beta___________________________________________________________ (betanb)______________________ 0.0000E+00 OP - (Fast_alpha_+_beam_beta)/(thermal_beta)_________________________________ (gammaft)_____________________ 1.4194E-01 OP - Thermal_beta____________________________________________________________ ______________________________ 3.1574E-02 OP - Thermal_poloidal_beta___________________________________________________ ______________________________ 1.0974E+00 OP - Thermal_toroidal_beta_(=_beta-exp)______________________________________ ______________________________ 3.2509E-02 OP - 2nd_stability_beta_:_beta_p_/_(R/a)_____________________________________ (eps*betap)___________________ 4.1770E-01 + (Fast_alpha_+_beam_beta)/(thermal_beta)_________________________________ (gammaft)_____________________ 1.4215E-01 OP + Thermal_beta____________________________________________________________ ______________________________ 2.9874E-02 OP + Thermal_poloidal_beta___________________________________________________ ______________________________ 1.1199E+00 OP + Thermal_toroidal_beta_(=_beta-exp)______________________________________ ______________________________ 3.0693E-02 OP + 2nd_stability_beta_:_beta_p_/_(R/a)_____________________________________ (eps*betap)___________________ 4.2636E-01 2nd_stability_beta_upper_limit__________________________________________ (epbetmax)____________________ 1.3800E+00 - Beta_g_coefficient______________________________________________________ (dnbeta)______________________ 4.7593E+00 - Normalised_thermal_beta_________________________________________________ ______________________________ 2.5527E+00 - Normalised_total_beta___________________________________________________ ______________________________ 2.9150E+00 - Normalised_toroidal_beta________________________________________________ (normalised_toroidal_beta)____ 3.0013E+00 - Limit_on_thermal_beta___________________________________________________ (betalim)_____________________ 5.8867E-02 - Plasma_thermal_energy_(J)_______________________________________________ ______________________________ 9.3900E+08 OP - Total_plasma_internal_energy_(J)________________________________________ (total_plasma_internal_energy) 1.0723E+09 OP - Electron_temperature_(keV)______________________________________________ (te)__________________________ 1.2238E+01 ITV - Electron_temperature_on_axis_(keV)______________________________________ (te0)_________________________ 2.5063E+01 - Ion_temperature_(keV)___________________________________________________ (ti)__________________________ 1.2238E+01 - Ion_temperature_on_axis_(keV)___________________________________________ (ti0)_________________________ 2.5063E+01 - Electron_temp.,_density_weighted_(keV)__________________________________ (ten)_________________________ 1.3514E+01 - Electron_density_(/m3)__________________________________________________ (dene)________________________ 8.0995E+19 ITV - Electron_density_on_axis_(/m3)__________________________________________ (ne0)_________________________ 1.0654E+20 OP - Line-averaged_electron_density_(/m3)____________________________________ (dnla)________________________ 8.9697E+19 OP + Beta_g_coefficient______________________________________________________ (dnbeta)______________________ 4.8794E+00 + Normalised_thermal_beta_________________________________________________ ______________________________ 2.5111E+00 + Normalised_total_beta___________________________________________________ ______________________________ 2.8680E+00 + Normalised_toroidal_beta________________________________________________ (normalised_toroidal_beta)____ 2.9466E+00 + Limit_on_thermal_beta___________________________________________________ (betalim)_____________________ 5.8050E-02 + Plasma_thermal_energy_(J)_______________________________________________ ______________________________ 9.4883E+08 OP + Total_plasma_internal_energy_(J)________________________________________ (total_plasma_internal_energy) 1.0837E+09 OP + Electron_temperature_(keV)______________________________________________ (te)__________________________ 1.2437E+01 ITV + Electron_temperature_on_axis_(keV)______________________________________ (te0)_________________________ 2.5614E+01 + Ion_temperature_(keV)___________________________________________________ (ti)__________________________ 1.2437E+01 + Ion_temperature_on_axis_(keV)___________________________________________ (ti0)_________________________ 2.5614E+01 + Electron_temp.,_density_weighted_(keV)__________________________________ (ten)_________________________ 1.3745E+01 + Electron_density_(/m3)__________________________________________________ (dene)________________________ 8.0594E+19 ITV + Electron_density_on_axis_(/m3)__________________________________________ (ne0)_________________________ 1.0601E+20 OP + Line-averaged_electron_density_(/m3)____________________________________ (dnla)________________________ 8.9254E+19 OP Line-averaged_electron_density_/_Greenwald_density______________________ (dnla_gw)_____________________ 1.2000E+00 OP - Ion_density_(/m3)_______________________________________________________ (dnitot)______________________ 7.2127E+19 OP - Fuel_density_(/m3)______________________________________________________ (deni)________________________ 6.5392E+19 OP - Total_impurity_density_with_Z_>_2_(no_He)_(/m3)_________________________ (dnz)_________________________ 4.6024E+16 OP - Helium_ion_density_(thermalised_ions_only)_(/m3)________________________ (dnalp)_______________________ 6.6709E+18 OP - Proton_density_(/m3)____________________________________________________ (dnprot)______________________ 1.8502E+16 OP + Ion_density_(/m3)_______________________________________________________ (dnitot)______________________ 7.1530E+19 OP + Fuel_density_(/m3)______________________________________________________ (deni)________________________ 6.4591E+19 OP + Total_impurity_density_with_Z_>_2_(no_He)_(/m3)_________________________ (dnz)_________________________ 4.5810E+16 OP + Helium_ion_density_(thermalised_ions_only)_(/m3)________________________ (dnalp)_______________________ 6.8731E+18 OP + Proton_density_(/m3)____________________________________________________ (dnprot)______________________ 1.9207E+16 OP Hot_beam_density_(/m3)__________________________________________________ (dnbeam)______________________ 0.0000E+00 OP - Density_limit_from_scaling_(/m3)________________________________________ (dnelimt)_____________________ 7.4748E+19 OP - Density_limit_(enforced)_(/m3)__________________________________________ (boundu(9)*dnelimt)___________ 7.4748E+19 OP - Helium_ion_density_(thermalised_ions_only)_/_electron_density___________ (ralpne)______________________ 8.2363E-02 ITV - H__concentration________________________________________________________ (fimp(01))____________________ 8.0759E-01 OP - He_concentration________________________________________________________ (fimp(02))____________________ 8.2363E-02 + Density_limit_from_scaling_(/m3)________________________________________ (dnelimt)_____________________ 7.4378E+19 OP + Density_limit_(enforced)_(/m3)__________________________________________ (boundu(9)*dnelimt)___________ 7.4378E+19 OP + Helium_ion_density_(thermalised_ions_only)_/_electron_density___________ (ralpne)______________________ 8.5281E-02 ITV + H__concentration________________________________________________________ (fimp(01))____________________ 8.0168E-01 OP + He_concentration________________________________________________________ (fimp(02))____________________ 8.5281E-02 Be_concentration________________________________________________________ (fimp(03))____________________ 0.0000E+00 C__concentration________________________________________________________ (fimp(04))____________________ 0.0000E+00 N__concentration________________________________________________________ (fimp(05))____________________ 0.0000E+00 @@ -390,190 +390,187 @@ Fe_concentration________________________________________________________ (fimp(10))____________________ 0.0000E+00 Ni_concentration________________________________________________________ (fimp(11))____________________ 0.0000E+00 Kr_concentration________________________________________________________ (fimp(12))____________________ 0.0000E+00 - Xe_concentration________________________________________________________ (fimp(13))____________________ 5.6323E-04 + Xe_concentration________________________________________________________ (fimp(13))____________________ 5.6341E-04 W__concentration________________________________________________________ (fimp(14))____________________ 5.0000E-06 - Average_mass_of_all_ions_(amu)__________________________________________ (aion)________________________ 2.7208E+00 OP - Effective_charge________________________________________________________ (zeff)________________________ 2.4867E+00 + Average_mass_of_all_ions_(amu)__________________________________________ (aion)________________________ 2.7265E+00 OP + Effective_charge________________________________________________________ (zeff)________________________ 2.4987E+00 Density_profile_factor__________________________________________________ (alphan)______________________ 1.0000E+00 Plasma_profile_model____________________________________________________ (ipedestal)___________________ 1 Density_pedestal_r/a_location___________________________________________ (rhopedn)_____________________ 9.4000E-01 - Electron_density_pedestal_height_(/m3)__________________________________ (neped)_______________________ 6.3536E+19 OP + Electron_density_pedestal_height_(/m3)__________________________________ (neped)_______________________ 6.3221E+19 OP Electron_density_at_pedestal_/_nGW______________________________________ (fgwped_out)__________________ 8.5000E-01 Temperature_pedestal_r/a_location_______________________________________ (rhopedt)_____________________ 9.4000E-01 Pedestal_scaling_switch_________________________________________________ (ieped)_______________________ 0 Electron_temp._pedestal_height_(keV)____________________________________ (teped)_______________________ 5.5000E+00 Electron_temp._at_separatrix_(keV)______________________________________ (tesep)_______________________ 1.0000E-01 - Electron_density_at_separatrix_(/m3)____________________________________ (nesep)_______________________ 3.7374E+19 + Electron_density_at_separatrix_(/m3)____________________________________ (nesep)_______________________ 3.7189E+19 Electron_density_at_separatrix_/_nGW____________________________________ (fgwsep_out)__________________ 5.0000E-01 Temperature_profile_index_______________________________________________ (alphat)______________________ 1.4500E+00 Temperature_profile_index_beta__________________________________________ (tbeta)_______________________ 2.0000E+00 - Old_ASDEX_model_________________________________________________________ (dlimit(1))___________________ 5.2107E+19 OP - Borrass_ITER_model_I____________________________________________________ (dlimit(2))___________________ 1.0653E+20 OP - Borrass_ITER_model_II___________________________________________________ (dlimit(3))___________________ 4.1967E+19 OP - JET_edge_radiation_model________________________________________________ (dlimit(4))___________________ 3.0404E+21 OP - JET_simplified_model____________________________________________________ (dlimit(5))___________________ 4.0490E+20 OP - Hugill-Murakami_Mq_model________________________________________________ (dlimit(6))___________________ 6.6879E+19 OP - Greenwald_model_________________________________________________________ (dlimit(7))___________________ 7.4748E+19 OP + Old_ASDEX_model_________________________________________________________ (dlimit(1))___________________ 5.2019E+19 OP + Borrass_ITER_model_I____________________________________________________ (dlimit(2))___________________ 1.0743E+20 OP + Borrass_ITER_model_II___________________________________________________ (dlimit(3))___________________ 4.2560E+19 OP + JET_edge_radiation_model________________________________________________ (dlimit(4))___________________ 2.9966E+21 OP + JET_simplified_model____________________________________________________ (dlimit(5))___________________ 4.2146E+20 OP + Hugill-Murakami_Mq_model________________________________________________ (dlimit(6))___________________ 6.6549E+19 OP + Greenwald_model_________________________________________________________ (dlimit(7))___________________ 7.4378E+19 OP Deuterium_fuel_fraction_________________________________________________ (fdeut)_______________________ 5.0000E-01 Tritium_fuel_fraction___________________________________________________ (ftrit)_______________________ 5.0000E-01 - Total_fusion_power_(MW)_________________________________________________ (powfmw)______________________ 1.5926E+03 OP - _=____D-T_fusion_power_(MW)_____________________________________________ (pdt)_________________________ 1.5907E+03 OP - __+___D-D_fusion_power_(MW)_____________________________________________ (pdd)_________________________ 1.8966E+00 OP + Total_fusion_power_(MW)_________________________________________________ (powfmw)______________________ 1.5971E+03 OP + _=____D-T_fusion_power_(MW)_____________________________________________ (pdt)_________________________ 1.5952E+03 OP + __+___D-D_fusion_power_(MW)_____________________________________________ (pdd)_________________________ 1.9174E+00 OP __+_D-He3_fusion_power_(MW)_____________________________________________ (pdhe3)_______________________ 0.0000E+00 OP - Alpha_power:_total_(MW)_________________________________________________ (palpmw)______________________ 3.1815E+02 OP + Alpha_power:_total_(MW)_________________________________________________ (palpmw)______________________ 3.1903E+02 OP Alpha_power:_beam-plasma_(MW)___________________________________________ (palpnb)______________________ 0.0000E+00 OP - Neutron_power_(MW)______________________________________________________ (pneutmw)_____________________ 1.2733E+03 OP - Charged_particle_power_(excluding_alphas)_(MW)__________________________ (pchargemw)___________________ 1.2322E+00 OP - Total_power_deposited_in_plasma_(MW)____________________________________ (tot_power_plasma)____________ 3.7935E+02 OP - Bremsstrahlung_radiation_power_(MW)_____________________________________ (pbrempv*vol)_________________ 6.4356E+01 OP - Line_radiation_power_(MW)_______________________________________________ (plinepv*vol)_________________ 1.4715E+02 OP - Synchrotron_radiation_power_(MW)________________________________________ (psyncpv*vol)_________________ 1.2665E+01 OP + Neutron_power_(MW)______________________________________________________ (pneutmw)_____________________ 1.2768E+03 OP + Charged_particle_power_(excluding_alphas)_(MW)__________________________ (pchargemw)___________________ 1.2453E+00 OP + Total_power_deposited_in_plasma_(MW)____________________________________ (tot_power_plasma)____________ 3.8036E+02 OP + Bremsstrahlung_radiation_power_(MW)_____________________________________ (pbrempv*vol)_________________ 6.4495E+01 OP + Line_radiation_power_(MW)_______________________________________________ (plinepv*vol)_________________ 1.4420E+02 OP + Synchrotron_radiation_power_(MW)________________________________________ (psyncpv*vol)_________________ 1.4563E+01 OP Synchrotron_wall_reflectivity_factor____________________________________ (ssync)_______________________ 6.0000E-01 Normalised_minor_radius_defining_'core'_________________________________ (coreradius)__________________ 7.5000E-01 Fraction_of_core_radiation_subtracted_from_P_L__________________________ (coreradiationfraction)_______ 6.0000E-01 - Radiation_power_from_inner_zone_(MW)____________________________________ (pinnerzoneradmw)_____________ 8.9376E+01 OP - Radiation_power_from_outer_zone_(MW)____________________________________ (pouterzoneradmw)_____________ 1.3480E+02 OP - SOL_radiation_power_as_imposed_by_f_rad_(MW)____________________________ (psolradmw)___________________ 0.0000E+00 OP - Total_radiation_power_from_inside_LCFS_(MW)_____________________________ (pradmw)______________________ 2.2417E+02 OP - LCFS_radiation_fraction_=_total_radiation_in_LCFS_/_total_power_deposite (rad_fraction_LCFS)___________ 5.9094E-01 OP - Nominal_mean_radiation_load_on_inside_surface_of_reactor_(MW/m2)________ (photon_wall)_________________ 1.7570E-01 OP + Radiation_power_from_inner_zone_(MW)____________________________________ (pinnerzoneradmw)_____________ 9.0174E+01 OP + Radiation_power_from_outer_zone_(MW)____________________________________ (pouterzoneradmw)_____________ 1.3309E+02 OP + Total_radiation_power_from_inside_LCFS_(MW)_____________________________ (pradmw)______________________ 2.2326E+02 OP + LCFS_radiation_fraction_=_total_radiation_in_LCFS_/_total_power_deposite (rad_fraction_LCFS)___________ 5.8698E-01 OP + Nominal_mean_radiation_load_on_inside_surface_of_reactor_(MW/m2)________ (photon_wall)_________________ 1.7498E-01 OP Peaking_factor_for_radiation_wall_load__________________________________ (peakfactrad)_________________ 3.3300E+00 IP Maximum_permitted_radiation_wall_load_(MW/m^2)__________________________ (maxradwallload)______________ 1.0000E+00 IP - Peak_radiation_wall_load_(MW/m^2)_______________________________________ (peakradwallload)_____________ 5.8507E-01 OP - Fast_alpha_particle_power_incident_on_the_first_wall_(MW)_______________ (palpfwmw)____________________ 1.5907E+01 OP - Nominal_mean_neutron_load_on_inside_surface_of_reactor_(MW/m2)__________ (wallmw)______________________ 9.9792E-01 OP - Radiation_power_from_SoL_(MW)___________________________________________ (pradsolmw)___________________ 1.2414E+02 OP - SoL_radiation_fraction_=_total_radiation_in_SoL_/_total_power_accross_se (rad_fraction_sol)____________ 8.0000E-01 IP - Radiation_fraction_total_=_SoL_+_LCFS_radiation_/_total_power_deposited_ (rad_fraction_total)__________ 9.1819E-01 OP - Power_incident_on_the_divertor_targets_(MW)_____________________________ (ptarmw)______________________ 3.1035E+01 OP + Peak_radiation_wall_load_(MW/m^2)_______________________________________ (peakradwallload)_____________ 5.8269E-01 OP + Fast_alpha_particle_power_incident_on_the_first_wall_(MW)_______________ (palpfwmw)____________________ 1.5952E+01 OP + Nominal_mean_neutron_load_on_inside_surface_of_reactor_(MW/m2)__________ (wallmw)______________________ 1.0007E+00 OP + Power_incident_on_the_divertor_targets_(MW)_____________________________ (ptarmw)______________________ 3.1419E+01 OP Fraction_of_power_to_the_lower_divertor_________________________________ (ftar)________________________ 1.0000E+00 IP Outboard_side_heat_flux_decay_length_(m)________________________________ (lambdaio)____________________ 1.5700E-03 OP Fraction_of_power_on_the_inner_targets__________________________________ (fio)_________________________ 4.1000E-01 OP Fraction_of_power_incident_on_the_lower_inner_target____________________ (fLI)_________________________ 4.1000E-01 OP Fraction_of_power_incident_on_the_lower_outer_target____________________ (fLO)_________________________ 5.9000E-01 OP - Power_incident_on_the_lower_inner_target_(MW)___________________________ (pLImw)_______________________ 1.2724E+01 OP - Power_incident_on_the_lower_outer_target_(MW)___________________________ (pLOmw)_______________________ 1.8311E+01 OP - Ohmic_heating_power_(MW)________________________________________________ (pohmmw)______________________ 6.6172E-01 OP + Power_incident_on_the_lower_inner_target_(MW)___________________________ (pLImw)_______________________ 1.2882E+01 OP + Power_incident_on_the_lower_outer_target_(MW)___________________________ (pLOmw)_______________________ 1.8537E+01 OP + Ohmic_heating_power_(MW)________________________________________________ (pohmmw)______________________ 6.3634E-01 OP Fraction_of_alpha_power_deposited_in_plasma_____________________________ (falpha)______________________ 9.5000E-01 - Fraction_of_alpha_power_to_electrons____________________________________ (falpe)_______________________ 7.2139E-01 - Fraction_of_alpha_power_to_ions_________________________________________ (falpi)_______________________ 2.7861E-01 - Ion_transport_(MW)______________________________________________________ (ptrimw)______________________ 1.3659E+02 OP - Electron_transport_(MW)_________________________________________________ (ptremw)______________________ 1.5338E+02 OP + Fraction_of_alpha_power_to_electrons____________________________________ (falpe)_______________________ 7.1893E-01 + Fraction_of_alpha_power_to_ions_________________________________________ (falpi)_______________________ 2.8107E-01 + Ion_transport_(MW)______________________________________________________ (ptrimw)______________________ 1.3645E+02 OP + Electron_transport_(MW)_________________________________________________ (ptremw)______________________ 1.5374E+02 OP Injection_power_to_ions_(MW)____________________________________________ (pinjimw)_____________________ 0.0000E+00 OP - Injection_power_to_electrons_(MW)_______________________________________ (pinjemw)_____________________ 7.5213E+01 OP + Injection_power_to_electrons_(MW)_______________________________________ (pinjemw)_____________________ 7.5398E+01 OP Ignited_plasma_switch_(0=not_ignited,_1=ignited)________________________ (ignite)______________________ 0 - Power_into_divertor_zone_via_charged_particles_(MW)_____________________ (pdivt)_______________________ 1.5518E+02 OP - Psep_/_R_ratio_(MW/m)___________________________________________________ (pdivt/rmajor)________________ 1.9397E+01 OP - Psep_Bt_/_qAR_ratio_(MWT/m)_____________________________________________ (pdivtbt/qar)_________________ 9.5570E+00 OP - ITER_1996_scaling:_nominal_(MW)_________________________________________ (pthrmw(1))___________________ 1.2449E+02 OP - ITER_1996_scaling:_upper_bound_(MW)_____________________________________ (pthrmw(2))___________________ 2.7464E+02 OP - ITER_1996_scaling:_lower_bound_(MW)_____________________________________ (pthrmw(3))___________________ 5.5673E+01 OP - ITER_1997_scaling_(1)_(MW)______________________________________________ (pthrmw(4))___________________ 2.0721E+02 OP - ITER_1997_scaling_(2)_(MW)______________________________________________ (pthrmw(5))___________________ 1.6581E+02 OP - Martin_2008_scaling:_nominal_(MW)_______________________________________ (pthrmw(6))___________________ 9.4410E+01 OP - Martin_2008_scaling:_95%_upper_bound_(MW)_______________________________ (pthrmw(7))___________________ 1.2366E+02 OP - Martin_2008_scaling:_95%_lower_bound_(MW)_______________________________ (pthrmw(8))___________________ 6.5163E+01 OP - Snipes_2000_scaling:_nominal_(MW)_______________________________________ (pthrmw(9))___________________ 6.5607E+01 OP - Snipes_2000_scaling:_upper_bound_(MW)___________________________________ (pthrmw(10))__________________ 9.6122E+01 OP - Snipes_2000_scaling:_lower_bound_(MW)___________________________________ (pthrmw(11))__________________ 4.4421E+01 OP - Snipes_2000_scaling_(closed_divertor):_nominal_(MW)_____________________ (pthrmw(12))__________________ 3.0395E+01 OP - Snipes_2000_scaling_(closed_divertor):_upper_bound_(MW)_________________ (pthrmw(13))__________________ 4.2193E+01 OP - Snipes_2000_scaling_(closed_divertor):_lower_bound_(MW)_________________ (pthrmw(14))__________________ 2.1742E+01 OP - Hubbard_2012_L-I_threshold_-_nominal_(MW)_______________________________ (pthrmw(15))__________________ 2.7728E+01 OP - Hubbard_2012_L-I_threshold_-_lower_bound_(MW)___________________________ (pthrmw(16))__________________ 1.4387E+01 OP - Hubbard_2012_L-I_threshold_-_upper_bound_(MW)___________________________ (pthrmw(17))__________________ 5.3439E+01 OP - Hubbard_2017_L-I_threshold______________________________________________ (pthrmw(18))__________________ 2.6004E+02 OP - Martin_2008_aspect_ratio_corrected_scaling:_nominal_(MW)________________ (pthrmw(19))__________________ 9.4410E+01 OP - Martin_2008_aspect_ratio_corrected_scaling:_95%_upper_bound_(MW)________ (pthrmw(20))__________________ 1.2366E+02 OP - Martin_2008_aspect_ratio_corrected_scaling:_95%_lower_bound_(MW)________ (pthrmw(21))__________________ 6.5163E+01 OP - L-H_threshold_power_(enforced)_(MW)_____________________________________ (boundl(103)*plhthresh)_______ 9.4410E+01 OP - L-H_threshold_power_(MW)________________________________________________ (plhthresh)___________________ 9.4410E+01 OP + Power_into_divertor_zone_via_charged_particles_(MW)_____________________ (pdivt)_______________________ 1.5710E+02 OP + Psep_/_R_ratio_(MW/m)___________________________________________________ (pdivt/rmajor)________________ 1.9637E+01 OP + Psep_Bt_/_qAR_ratio_(MWT/m)_____________________________________________ (pdivtbt/qar)_________________ 9.6274E+00 OP + ITER_1996_scaling:_nominal_(MW)_________________________________________ (pthrmw(1))___________________ 1.2831E+02 OP + ITER_1996_scaling:_upper_bound_(MW)_____________________________________ (pthrmw(2))___________________ 2.8272E+02 OP + ITER_1996_scaling:_lower_bound_(MW)_____________________________________ (pthrmw(3))___________________ 5.7453E+01 OP + ITER_1997_scaling_(1)_(MW)______________________________________________ (pthrmw(4))___________________ 2.1237E+02 OP + ITER_1997_scaling_(2)_(MW)______________________________________________ (pthrmw(5))___________________ 1.7028E+02 OP + Martin_2008_scaling:_nominal_(MW)_______________________________________ (pthrmw(6))___________________ 9.6473E+01 OP + Martin_2008_scaling:_95%_upper_bound_(MW)_______________________________ (pthrmw(7))___________________ 1.2643E+02 OP + Martin_2008_scaling:_95%_lower_bound_(MW)_______________________________ (pthrmw(8))___________________ 6.6515E+01 OP + Snipes_2000_scaling:_nominal_(MW)_______________________________________ (pthrmw(9))___________________ 6.7125E+01 OP + Snipes_2000_scaling:_upper_bound_(MW)___________________________________ (pthrmw(10))__________________ 9.8433E+01 OP + Snipes_2000_scaling:_lower_bound_(MW)___________________________________ (pthrmw(11))__________________ 4.5409E+01 OP + Snipes_2000_scaling_(closed_divertor):_nominal_(MW)_____________________ (pthrmw(12))__________________ 3.0806E+01 OP + Snipes_2000_scaling_(closed_divertor):_upper_bound_(MW)_________________ (pthrmw(13))__________________ 4.2835E+01 OP + Snipes_2000_scaling_(closed_divertor):_lower_bound_(MW)_________________ (pthrmw(14))__________________ 2.1999E+01 OP + Hubbard_2012_L-I_threshold_-_nominal_(MW)_______________________________ (pthrmw(15))__________________ 2.7510E+01 OP + Hubbard_2012_L-I_threshold_-_lower_bound_(MW)___________________________ (pthrmw(16))__________________ 1.4304E+01 OP + Hubbard_2012_L-I_threshold_-_upper_bound_(MW)___________________________ (pthrmw(17))__________________ 5.2909E+01 OP + Hubbard_2017_L-I_threshold______________________________________________ (pthrmw(18))__________________ 2.6105E+02 OP + Martin_2008_aspect_ratio_corrected_scaling:_nominal_(MW)________________ (pthrmw(19))__________________ 9.6473E+01 OP + Martin_2008_aspect_ratio_corrected_scaling:_95%_upper_bound_(MW)________ (pthrmw(20))__________________ 1.2643E+02 OP + Martin_2008_aspect_ratio_corrected_scaling:_95%_lower_bound_(MW)________ (pthrmw(21))__________________ 6.6515E+01 OP + L-H_threshold_power_(enforced)_(MW)_____________________________________ (boundl(103)*plhthresh)_______ 9.6473E+01 OP + L-H_threshold_power_(MW)________________________________________________ (plhthresh)___________________ 9.6473E+01 OP Confinement_scaling_law_________________________________________________ (tauelaw)_____________________ "IPB98(y,2)" - Confinement_H_factor____________________________________________________ (hfact)_______________________ 1.1808E+00 ITV - Global_thermal_energy_confinement_time_(s)______________________________ (taueff)______________________ 3.2379E+00 - Ion_energy_confinement_time_(s)_________________________________________ (tauei)_______________________ 3.2379E+00 - Electron_energy_confinement_time_(s)____________________________________ (tauee)_______________________ 3.2379E+00 - n.tau_=_Volume-average_electron_density_x_Energy_confinement_time_(s/m3) (dntau)_______________________ 2.6225E+20 OP - Triple_product__(keV_s/m3)______________________________________________ (dntau*te)____________________ 3.2095E+21 OP - Transport_loss_power_assumed_in_scaling_law_(MW)________________________ (powerht)_____________________ 2.8997E+02 OP + Confinement_H_factor____________________________________________________ (hfact)_______________________ 1.1947E+00 ITV + Global_thermal_energy_confinement_time_(s)______________________________ (taueff)______________________ 3.2694E+00 + Ion_energy_confinement_time_(s)_________________________________________ (tauei)_______________________ 3.2694E+00 + Electron_energy_confinement_time_(s)____________________________________ (tauee)_______________________ 3.2694E+00 + n.tau_=_Volume-average_electron_density_x_Energy_confinement_time_(s/m3) (dntau)_______________________ 2.6349E+20 OP + Triple_product__(keV_s/m3)______________________________________________ (dntau*te)____________________ 3.2771E+21 OP + Transport_loss_power_assumed_in_scaling_law_(MW)________________________ (powerht)_____________________ 2.9018E+02 OP Switch_for_radiation_loss_term_usage_in_power_balance___________________ (iradloss)____________________ 1 - Radiation_power_subtracted_from_plasma_power_balance_(MW)_______________ ______________________________ 8.9376E+01 OP - Alpha_particle_confinement_time_(s)_____________________________________ (taup)________________________ 2.2315E+01 - Alpha_particle/energy_confinement_time_ratio____________________________ (taup/taueff)_________________ 6.8919E+00 + Radiation_power_subtracted_from_plasma_power_balance_(MW)_______________ ______________________________ 9.0174E+01 OP + Alpha_particle_confinement_time_(s)_____________________________________ (taup)________________________ 2.2928E+01 + Alpha_particle/energy_confinement_time_ratio____________________________ (taup/taueff)_________________ 7.0130E+00 Lower_limit_on_taup/taueff______________________________________________ (taulimit)____________________ 5.0000E+00 - Total_energy_confinement_time_including_radiation_loss_(s)______________ (total_energy_conf_time)______ 2.8266E+00 - Normalized_plasma_pressure_beta_as_defined_by_McDonald_et_al____________ (beta_mcdonald)_______________ 3.7123E-02 OP - Normalized_ion_Larmor_radius____________________________________________ (rho_star)____________________ 2.0262E-03 OP - Normalized_collisionality_______________________________________________ (nu_star)_____________________ 3.7167E-03 OP + Total_energy_confinement_time_including_radiation_loss_(s)______________ (total_energy_conf_time)______ 2.8492E+00 + Normalized_plasma_pressure_beta_as_defined_by_McDonald_et_al____________ (beta_mcdonald)_______________ 3.5056E-02 OP + Normalized_ion_Larmor_radius____________________________________________ (rho_star)____________________ 1.9759E-03 OP + Normalized_collisionality_______________________________________________ (nu_star)_____________________ 3.7312E-03 OP Volume_measure_of_elongation____________________________________________ (kappaa_IPB)__________________ 1.6815E+00 OP - Total_volt-second_requirement_(Wb)______________________________________ (vsstt)_______________________ 5.6973E+02 OP - Inductive_volt-seconds_(Wb)_____________________________________________ (vsind)_______________________ 2.3365E+02 OP + Total_volt-second_requirement_(Wb)______________________________________ (vsstt)_______________________ 5.6164E+02 OP + Inductive_volt-seconds_(Wb)_____________________________________________ (vsind)_______________________ 2.3500E+02 OP Ejima_coefficient_______________________________________________________ (gamma)_______________________ 3.0000E-01 - Start-up_resistive_(Wb)_________________________________________________ (vsres)_______________________ 5.0362E+01 OP - Flat-top_resistive_(Wb)_________________________________________________ (vsbrn)_______________________ 2.8572E+02 OP + Start-up_resistive_(Wb)_________________________________________________ (vsres)_______________________ 5.0113E+01 OP + Flat-top_resistive_(Wb)_________________________________________________ (vsbrn)_______________________ 2.7652E+02 OP bootstrap_current_fraction_multiplier___________________________________ (cboot)_______________________ 1.0000E+00 - Bootstrap_fraction_(ITER_1989)__________________________________________ (bscf_iter89)_________________ 3.6783E-01 - Bootstrap_fraction_(Sauter_et_al)_______________________________________ (bscf_sauter)_________________ 4.1436E-01 - Bootstrap_fraction_(Nevins_et_al)_______________________________________ (bscf_nevins)_________________ 3.4713E-01 - Bootstrap_fraction_(Wilson)_____________________________________________ (bscf_wilson)_________________ 4.3036E-01 - Diamagnetic_fraction_(Hender)___________________________________________ (diacf_hender)________________ 1.2877E-02 - Diamagnetic_fraction_(SCENE)____________________________________________ (diacf_scene)_________________ 1.1680E-02 - Pfirsch-Schlueter_fraction_(SCENE)______________________________________ (pscf_scene)__________________ -3.2449E-03 - Bootstrap_fraction_(enforced)___________________________________________ (bootipf.)____________________ 4.1436E-01 + Bootstrap_fraction_(ITER_1989)__________________________________________ (bscf_iter89)_________________ 3.7036E-01 + Bootstrap_fraction_(Sauter_et_al)_______________________________________ (bscf_sauter)_________________ 4.1885E-01 + Bootstrap_fraction_(Nevins_et_al)_______________________________________ (bscf_nevins)_________________ 3.4880E-01 + Bootstrap_fraction_(Wilson)_____________________________________________ (bscf_wilson)_________________ 4.3324E-01 + Diamagnetic_fraction_(Hender)___________________________________________ (diacf_hender)________________ 1.2186E-02 + Diamagnetic_fraction_(SCENE)____________________________________________ (diacf_scene)_________________ 1.1246E-02 + Pfirsch-Schlueter_fraction_(SCENE)______________________________________ (pscf_scene)__________________ -3.0708E-03 + Bootstrap_fraction_(enforced)___________________________________________ (bootipf.)____________________ 4.1885E-01 Diamagnetic_fraction_(enforced)_________________________________________ (diaipf.)_____________________ 0.0000E+00 Pfirsch-Schlueter_fraction_(enforced)___________________________________ (psipf.)______________________ 0.0000E+00 - Loop_voltage_during_burn_(V)____________________________________________ (vburn)_______________________ 3.9627E-02 OP - Plasma_resistance_(ohm)_________________________________________________ (rplas)_______________________ 4.0562E-09 OP - Resistive_diffusion_time_(s)____________________________________________ (res_time)____________________ 3.0010E+03 OP - Plasma_inductance_(H)___________________________________________________ (rlp)_________________________ 1.3992E-05 OP + Loop_voltage_during_burn_(V)____________________________________________ (vburn)_______________________ 3.8296E-02 OP + Plasma_resistance_(ohm)_________________________________________________ (rplas)_______________________ 3.9735E-09 OP + Resistive_diffusion_time_(s)____________________________________________ (res_time)____________________ 3.0634E+03 OP + Plasma_inductance_(H)___________________________________________________ (rlp)_________________________ 1.4143E-05 OP Coefficient_for_sawtooth_effects_on_burn_V-s_requirement________________ (csawth)______________________ 1.0000E+00 Ratio_of_He_and_pellet_particle_confinement_times_______________________ (tauratio)____________________ 1.0000E+00 - Fuelling_rate_(nucleus-pairs/s)_________________________________________ (qfuel)_______________________ 3.3502E+21 OP - Fuel_burn-up_rate_(reactions/s)_________________________________________ (rndfuel)_____________________ 5.6770E+20 OP - Burn-up_fraction________________________________________________________ (burnup)______________________ 1.6946E-01 + Fuelling_rate_(nucleus-pairs/s)_________________________________________ (qfuel)_______________________ 3.2444E+21 OP + Fuel_burn-up_rate_(reactions/s)_________________________________________ (rndfuel)_____________________ 5.6930E+20 OP + Burn-up_fraction________________________________________________________ (burnup)______________________ 1.7547E-01 # Energy confinement times, and required H-factors : # # Current Drive System # Current_drive_efficiency_model__________________________________________ (iefrf)_______________________ 10 + Secondary_current_drive_efficiency_model________________________________ (iefrffix)____________________ 0 + Ratio_of_power_for_flat-top_to_start-up_(MW)____________________________ (startupratio)________________ 1.0000E+00 Auxiliary_power_used_for_plasma_heating_only_(MW)_______________________ (pheat)_______________________ 7.5000E+01 - Power_injected_for_current_drive_(MW)___________________________________ (pcurrentdrivemw)_____________ 2.1293E-01 + Power_injected_for_current_drive_(MW)___________________________________ (pcurrentdrivemw)_____________ 3.9779E-01 Maximum_Allowed_Bootstrap_current_fraction______________________________ (bscfmax)_____________________ 9.5000E-01 - Fusion_gain_factor_Q____________________________________________________ (bigq)________________________ 2.0990E+01 OP - Auxiliary_current_drive_(A)_____________________________________________ (auxiliary_cd)________________ 9.8586E+03 OP - Current_drive_efficiency_(A/W)__________________________________________ (effcd)_______________________ 4.6299E-02 OP + Fusion_gain_factor_Q____________________________________________________ (bigq)________________________ 2.1005E+01 OP + Auxiliary_current_drive_(A)_____________________________________________ (auxiliary_cd)________________ 1.8509E+04 OP + Current_drive_efficiency_(A/W)__________________________________________ (effcd)_______________________ 4.6530E-02 OP Normalised_current_drive_efficiency,_gamma_(10^20_A/W-m2)_______________ (gamcd)_______________________ 3.0000E-01 OP Wall_plug_to_injector_efficiency________________________________________ (etacd)_______________________ 5.0000E-01 ECRH_plasma_heating_efficiency__________________________________________ (gamma_ecrh)__________________ 3.0000E-01 - Bootstrap_fraction______________________________________________________ (bootipf)_____________________ 4.1436E-01 + Bootstrap_fraction______________________________________________________ (bootipf)_____________________ 4.1885E-01 Diamagnetic_fraction____________________________________________________ (diaipf)______________________ 0.0000E+00 Pfirsch-Schlueter_fraction______________________________________________ (psipf)_______________________ 0.0000E+00 - Auxiliary_current_drive_fraction________________________________________ (faccd)_______________________ 5.9038E-04 - Inductive_fraction______________________________________________________ (facoh)_______________________ 5.8505E-01 - Total___________________________________________________________________ (plasipf+faccd+facoh)_________ 1.0000E+00 - Fraction_of_the_plasma_current_produced_by_non-inductive_means__________ (fvsbrnni)____________________ 4.1495E-01 ITV - Electron_cyclotron_injected_power_(MW)__________________________________ (echpwr)______________________ 7.5213E+01 OP + Auxiliary_current_drive_fraction________________________________________ (faccd)_______________________ 1.1139E-03 + Inductive_fraction______________________________________________________ (facoh)_______________________ 5.8003E-01 + Total___________________________________________________________________ (current_drive_variables.plasi 1.0000E+00 + Fraction_of_the_plasma_current_produced_by_non-inductive_means__________ (fvsbrnni)____________________ 4.1997E-01 ITV + Electron_cyclotron_injected_power_(MW)__________________________________ (echpwr)______________________ 7.5398E+01 OP Maximum_allowable_ECRH_power_(MW)_______________________________________ (pinjalw)_____________________ 2.0000E+02 ECH_wall_plug_efficiency________________________________________________ (etaech)______________________ 5.0000E-01 - ECH_wall_plug_power_(MW)________________________________________________ (echwpow)_____________________ 1.5043E+02 OP - Total_V-s_capability_of_Central_Solenoid/PF_coils_(Wb)__________________ (abs(vstot))__________________ 5.8939E+02 - Required_volt-seconds_during_start-up_(Wb)______________________________ (vssoft)______________________ 2.8401E+02 - Available_volt-seconds_during_burn_(Wb)_________________________________ (vsmax)_______________________ 2.8571E+02 + ECH_wall_plug_power_(MW)________________________________________________ (echwpow)_____________________ 1.5080E+02 OP + Total_V-s_capability_of_Central_Solenoid/PF_coils_(Wb)__________________ (abs(vstot))__________________ 5.8162E+02 + Required_volt-seconds_during_start-up_(Wb)______________________________ (vssoft)______________________ 2.8512E+02 + Available_volt-seconds_during_burn_(Wb)_________________________________ (vsmax)_______________________ 2.7652E+02 # Times # Initial_charge_time_for_CS_from_zero_current_(s)________________________ (tramp)_______________________ 5.0000E+02 - Plasma_current_ramp-up_time_(s)_________________________________________ (tohs)________________________ 1.6699E+02 + Plasma_current_ramp-up_time_(s)_________________________________________ (tohs)________________________ 1.6616E+02 Heating_time_(s)________________________________________________________ (theat)_______________________ 1.0000E+01 - Burn_time_(s)___________________________________________________________ (tburn)_______________________ 7.1999E+03 OP - Reset_time_to_zero_current_for_CS_(s)___________________________________ (tqnch)_______________________ 1.6699E+02 + Burn_time_(s)___________________________________________________________ (tburn)_______________________ 7.2104E+03 OP + Reset_time_to_zero_current_for_CS_(s)___________________________________ (tqnch)_______________________ 1.6616E+02 Time_between_pulses_(s)_________________________________________________ (tdwell)______________________ 1.8000E+03 - Total_plant_cycle_time_(s)______________________________________________ (tcycle)______________________ 9.8441E+03 OP - kallenbach_switch_______________________________________________________ (kallenbach_switch)___________ 0 + Total_plant_cycle_time_(s)______________________________________________ (tcycle)______________________ 9.8529E+03 OP # Radial Build # TF_coil_radial_placement_switch_________________________________________ (tf_in_cs)____________________ 0 TF_coil_radial_placement_switch_________________________________________ (tf_in_cs)____________________ 0 - Machine_bore_(m)________________________________________________________ (bore)________________________ 2.0588E+00 ITV - CS_radial_thickness_(m)_________________________________________________ (ohcth)_______________________ 5.5551E-01 ITV - CS_precompression_(m)___________________________________________________ (precomp)_____________________ 6.6224E-02 + Machine_bore_(m)________________________________________________________ (bore)________________________ 2.0251E+00 ITV + CS_radial_thickness_(m)_________________________________________________ (ohcth)_______________________ 5.5625E-01 ITV + CS_precompression_(m)___________________________________________________ (precomp)_____________________ 6.7181E-02 CS_precompresion_to_TF_coil_radial_gap_(m)______________________________ (gapoh)_______________________ 8.0000E-02 - TF_coil_inboard_leg_(m)_________________________________________________ (tfcth)_______________________ 8.6482E-01 ITV + TF_coil_inboard_leg_(m)_________________________________________________ (tfcth)_______________________ 8.9679E-01 ITV TF_coil_inboard_leg_insulation_gap_(m)__________________________________ (tftsgap)_____________________ 5.0000E-02 Thermal_shield,_inboard_(m)_____________________________________________ (thshield_ib)_________________ 5.0000E-02 thermal_shield_to_vessel_radial_gap_(m)_________________________________ (gapds)_______________________ 2.0000E-02 @@ -588,10 +585,10 @@ Outboard_blanket_radial_thickness_(m)___________________________________ (blnkoth)_____________________ 1.0000E+00 Outer_radiation_shield_radial_thickness_(m)_____________________________ (shldoth)_____________________ 8.0000E-01 Outboard_vacuum_vessel_radial_thickness_(m)_____________________________ (d_vv_out)____________________ 3.0000E-01 - Vessel_to_TF_radial_gap_(m)_____________________________________________ (gapsto)______________________ 1.3976E+00 + Vessel_to_TF_radial_gap_(m)_____________________________________________ (gapsto)______________________ 1.3815E+00 Thermal_shield,_outboard_(m)____________________________________________ (thshield_ob)_________________ 5.0000E-02 Gap_(m)_________________________________________________________________ (tftsgap)_____________________ 5.0000E-02 - TF_coil_outboard_leg_radial_thickness_(m)_______________________________ (tfthko)______________________ 8.6482E-01 + TF_coil_outboard_leg_radial_thickness_(m)_______________________________ (tfthko)______________________ 8.9679E-01 # Vertical Build # Divertor_null_switch____________________________________________________ (i_single_null)_______________ 1 Vessel_-_TF_coil_vertical_gap_(m)_______________________________________ (vgap2)_______________________ 1.6300E-01 @@ -644,316 +641,314 @@ # TF coils # Allowable_maximum_shear_stress_in_TF_coil_case_(Tresca_criterion)_(Pa)__ (sig_tf_case_max)_____________ 7.5000E+08 Allowable_maximum_shear_stress_in_TF_coil_conduit_(Tresca_criterion)_(Pa (sig_tf_wp_max)_______________ 7.5000E+08 - WP_transverse_modulus_(GPa)_____________________________________________ (eyoung_wp_trans*1.0d-9)______ 4.4610E+01 OP - WP_vertical_modulus_(GPa)_______________________________________________ (eyoung_wp_axial*1.0d-9)______ 1.1695E+02 OP - WP_transverse_Poissons_ratio____________________________________________ (poisson_wp_trans)____________ 3.0406E-01 OP - WP_vertical-transverse_Pois._rat._______________________________________ (poisson_wp_axial)____________ 3.1562E-01 OP + WP_transverse_modulus_(GPa)_____________________________________________ (eyoung_wp_trans*1.0d-9)______ 4.4524E+01 OP + WP_vertical_modulus_(GPa)_______________________________________________ (eyoung_wp_axial*1.0d-9)______ 1.1613E+02 OP + WP_transverse_Poissons_ratio____________________________________________ (poisson_wp_trans)____________ 3.0405E-01 OP + WP_vertical-transverse_Pois._rat._______________________________________ (poisson_wp_axial)____________ 3.1537E-01 OP Radial____stress_at_maximum_shear_of_layer_1_(Pa)_______________________ (sig_tf_r_max(1))_____________ 4.8849E-08 - toroidal__stress_at_maximum_shear_of_layer_1_(Pa)_______________________ (sig_tf_t_max(1))_____________ -4.9590E+08 - Vertical__stress_at_maximum_shear_of_layer_1_(Pa)_______________________ (sig_tf_z_max(1))_____________ 2.5410E+08 - Von-Mises_stress_at_maximum_shear_of_layer_1_(Pa)_______________________ (sig_tf_vmises_max(1))________ 6.6068E+08 - Maximum_shear_stress_for_the_Tresca_yield_criterion_1_(Pa)______________ (sig_tf_tresca_max(1))________ 7.5000E+08 - Radial____stress_at_maximum_shear_of_layer_2_(Pa)_______________________ (sig_tf_r_max(2))_____________ -1.2706E+08 - toroidal__stress_at_maximum_shear_of_layer_2_(Pa)_______________________ (sig_tf_t_max(2))_____________ -3.8248E+08 - Vertical__stress_at_maximum_shear_of_layer_2_(Pa)_______________________ (sig_tf_z_max(2))_____________ 2.5410E+08 - Von-Mises_stress_at_maximum_shear_of_layer_2_(Pa)_______________________ (sig_tf_vmises_max(2))________ 5.5502E+08 - Maximum_shear_stress_for_the_Tresca_yield_criterion_2_(Pa)______________ (sig_tf_tresca_max(2))________ 6.3658E+08 - Radial____stress_at_maximum_shear_of_layer_3_(Pa)_______________________ (sig_tf_r_max(3))_____________ 6.3351E+06 - toroidal__stress_at_maximum_shear_of_layer_3_(Pa)_______________________ (sig_tf_t_max(3))_____________ -3.7469E+08 - Vertical__stress_at_maximum_shear_of_layer_3_(Pa)_______________________ (sig_tf_z_max(3))_____________ 1.4165E+08 - Von-Mises_stress_at_maximum_shear_of_layer_3_(Pa)_______________________ (sig_tf_vmises_max(3))________ 4.6210E+08 - Maximum_shear_stress_for_the_Tresca_yield_criterion_3_(Pa)______________ (sig_tf_tresca_max(3))________ 5.1634E+08 - Maximum_radial_deflection_at_midplane_(m)_______________________________ (deflect)_____________________ -6.5297E-03 OP - Vertical_strain_on_casing_______________________________________________ (casestr)_____________________ 1.2395E-03 OP - Radial_strain_on_insulator______________________________________________ (insstrain)___________________ -6.3417E-03 OP + toroidal__stress_at_maximum_shear_of_layer_1_(Pa)_______________________ (sig_tf_t_max(1))_____________ -4.8992E+08 + Vertical__stress_at_maximum_shear_of_layer_1_(Pa)_______________________ (sig_tf_z_max(1))_____________ 2.5886E+08 + Von-Mises_stress_at_maximum_shear_of_layer_1_(Pa)_______________________ (sig_tf_vmises_max(1))________ 6.5868E+08 + Maximum_shear_stress_for_the_Tresca_yield_criterion_1_(Pa)______________ (sig_tf_tresca_max(1))________ 7.4879E+08 + Radial____stress_at_maximum_shear_of_layer_2_(Pa)_______________________ (sig_tf_r_max(2))_____________ -1.4373E+08 + toroidal__stress_at_maximum_shear_of_layer_2_(Pa)_______________________ (sig_tf_t_max(2))_____________ -3.7794E+08 + Vertical__stress_at_maximum_shear_of_layer_2_(Pa)_______________________ (sig_tf_z_max(2))_____________ 2.5886E+08 + Von-Mises_stress_at_maximum_shear_of_layer_2_(Pa)_______________________ (sig_tf_vmises_max(2))________ 5.5469E+08 + Maximum_shear_stress_for_the_Tresca_yield_criterion_2_(Pa)______________ (sig_tf_tresca_max(2))________ 6.3680E+08 + Radial____stress_at_maximum_shear_of_layer_3_(Pa)_______________________ (sig_tf_r_max(3))_____________ 6.2042E+06 + toroidal__stress_at_maximum_shear_of_layer_3_(Pa)_______________________ (sig_tf_t_max(3))_____________ -3.6694E+08 + Vertical__stress_at_maximum_shear_of_layer_3_(Pa)_______________________ (sig_tf_z_max(3))_____________ 1.4431E+08 + Von-Mises_stress_at_maximum_shear_of_layer_3_(Pa)_______________________ (sig_tf_vmises_max(3))________ 4.5654E+08 + Maximum_shear_stress_for_the_Tresca_yield_criterion_3_(Pa)______________ (sig_tf_tresca_max(3))________ 5.1125E+08 + Maximum_radial_deflection_at_midplane_(m)_______________________________ (deflect)_____________________ -6.3578E-03 OP + Vertical_strain_on_casing_______________________________________________ (casestr)_____________________ 1.2627E-03 OP + Radial_strain_on_insulator______________________________________________ (insstrain)___________________ -7.1743E-03 OP Conductor_technology____________________________________________________ (i_tf_sup)____________________ 1 Superconductor_material_________________________________________________ (i_tf_sc_mat)_________________ 1 Presence_of_TF_demountable_joints_______________________________________ (itart)_______________________ 0 TF_inboard_leg_support_strategy_________________________________________ (i_tf_bucking)________________ 1 Number_of_TF_coils______________________________________________________ (n_tf)________________________ 16 - Inboard_leg_centre_radius_(m)___________________________________________ (r_tf_inboard_mid)____________ 3.1929E+00 OP + Inboard_leg_centre_radius_(m)___________________________________________ (r_tf_inboard_mid)____________ 3.1769E+00 OP Outboard_leg_centre_radius_(m)__________________________________________ (r_tf_outboard_mid)___________ 1.4985E+01 OP - Total_inboard_leg_radial_thickness_(m)__________________________________ (tfcth)_______________________ 8.6482E-01 ITV - Total_outboard_leg_radial_thickness_(m)_________________________________ (tfthko)______________________ 8.6482E-01 + Total_inboard_leg_radial_thickness_(m)__________________________________ (tfcth)_______________________ 8.9679E-01 ITV + Total_outboard_leg_radial_thickness_(m)_________________________________ (tfthko)______________________ 8.9679E-01 Outboard_leg_toroidal_thickness_(m)_____________________________________ (tftort)______________________ 1.4145E+00 OP Maximum_inboard_edge_height_(m)_________________________________________ (hmax)________________________ 8.8182E+00 OP - Mean_coil_circumference_(including_inboard_leg_length)_(m)______________ (tfleng)______________________ 4.7815E+01 OP + Mean_coil_circumference_(including_inboard_leg_length)_(m)______________ (tfleng)______________________ 4.7891E+01 OP Vertical_TF_shape_______________________________________________________ (i_tf_shape)__________________ 1 TF_coil_arc_point_0_R_(m)_______________________________________________ (xarc(1))_____________________ 3.6253E+00 TF_coil_arc_point_0_Z_(m)_______________________________________________ (yarc(1))_____________________ 4.5506E+00 TF_coil_arc_point_1_R_(m)_______________________________________________ (xarc(2))_____________________ 7.4667E+00 TF_coil_arc_point_1_Z_(m)_______________________________________________ (yarc(2))_____________________ 7.5843E+00 - TF_coil_arc_point_2_R_(m)_______________________________________________ (xarc(3))_____________________ 1.4552E+01 + TF_coil_arc_point_2_R_(m)_______________________________________________ (xarc(3))_____________________ 1.4536E+01 TF_coil_arc_point_2_Z_(m)_______________________________________________ (yarc(3))_____________________ 0.0000E+00 TF_coil_arc_point_3_R_(m)_______________________________________________ (xarc(4))_____________________ 7.4667E+00 TF_coil_arc_point_3_Z_(m)_______________________________________________ (yarc(4))_____________________ -8.8182E+00 TF_coil_arc_point_4_R_(m)_______________________________________________ (xarc(5))_____________________ 3.6253E+00 TF_coil_arc_point_4_Z_(m)_______________________________________________ (yarc(5))_____________________ -5.2909E+00 - TF_cross-section_(total)_(m2)___________________________________________ (tfareain)____________________ 1.7350E+01 - Total_steel_cross-section_(m2)__________________________________________ (a_tf_steel*n_tf)_____________ 1.2056E+01 - Total_steel_TF_fraction_________________________________________________ (f_tf_steel)__________________ 6.9487E-01 - Total_Insulation_cross-section_(total)_(m2)_____________________________ (a_tf_ins*n_tf)_______________ 8.6240E-01 - Total_Insulation_fraction_______________________________________________ (f_tf_ins)____________________ 4.9707E-02 - Casing_cross_section_area_(per_leg)_(m2)________________________________ (acasetf)_____________________ 5.0816E-01 + TF_cross-section_(total)_(m2)___________________________________________ (tfareain)____________________ 1.7901E+01 + Total_steel_cross-section_(m2)__________________________________________ (a_tf_steel*n_tf)_____________ 1.2650E+01 + Total_steel_TF_fraction_________________________________________________ (f_tf_steel)__________________ 7.0666E-01 + Total_Insulation_cross-section_(total)_(m2)_____________________________ (a_tf_ins*n_tf)_______________ 8.5501E-01 + Total_Insulation_fraction_______________________________________________ (f_tf_ins)____________________ 4.7763E-02 + Casing_cross_section_area_(per_leg)_(m2)________________________________ (acasetf)_____________________ 5.4835E-01 Inboard_leg_case_plasma_side_wall_thickness_(m)_________________________ (casthi)______________________ 6.0000E-02 - Inboard_leg_case_inboard_"nose"_thickness_(m)___________________________ (thkcas)______________________ 2.2965E-01 ITV + Inboard_leg_case_inboard_"nose"_thickness_(m)___________________________ (thkcas)______________________ 2.6741E-01 ITV Inboard_leg_case_sidewall_thickness_at_its_narrowest_point_(m)__________ (casths)______________________ 5.0000E-02 - External_case_mass_per_coil_(kg)________________________________________ (whtcas)______________________ 4.9722E+05 OP - WP_cross_section_area_with_insulation_and_insertion_(per_coil)_(m2)_____ (awpc)________________________ 5.7619E-01 - WP_cross_section_area_(per_coil)_(m2)___________________________________ (aswp)________________________ 5.1826E-01 - Winding_pack_radial_thickness_(m)_______________________________________ (dr_tf_wp)____________________ 5.0550E-01 ITV - Winding_pack_toroidal_width_1_(m)_______________________________________ (wwp1)________________________ 1.1901E+00 OP - Winding_pack_toroidal_width_2_(m)_______________________________________ (wwp2)________________________ 1.0896E+00 OP + External_case_mass_per_coil_(kg)________________________________________ (whtcas)______________________ 5.3720E+05 OP + WP_cross_section_area_with_insulation_and_insertion_(per_coil)_(m2)_____ (awpc)________________________ 5.7047E-01 + WP_cross_section_area_(per_coil)_(m2)___________________________________ (aswp)________________________ 5.1268E-01 + Winding_pack_radial_thickness_(m)_______________________________________ (dr_tf_wp)____________________ 4.9973E-01 ITV + Winding_pack_toroidal_width_1_(m)_______________________________________ (wwp1)________________________ 1.1913E+00 OP + Winding_pack_toroidal_width_2_(m)_______________________________________ (wwp2)________________________ 1.0919E+00 OP Ground_wall_insulation_thickness_(m)____________________________________ (tinstf)______________________ 8.0000E-03 Winding_pack_insertion_gap_(m)__________________________________________ (tfinsgap)____________________ 1.0000E-02 - Steel_WP_cross-section_(total)_(m2)_____________________________________ (aswp*n_tf)___________________ 3.9252E+00 - Steel_WP_fraction_______________________________________________________ (aswp/awpc)___________________ 4.2577E-01 - Insulation_WP_fraction__________________________________________________ (aiwp/awpc)___________________ 4.9412E-02 - Cable_WP_fraction_______________________________________________________ ((awpc-aswp-aiwp)/awpc)_______ 5.2482E-01 + Steel_WP_cross-section_(total)_(m2)_____________________________________ (aswp*n_tf)___________________ 3.8763E+00 + Steel_WP_fraction_______________________________________________________ (aswp/awpc)___________________ 4.2469E-01 + Insulation_WP_fraction__________________________________________________ (aiwp/awpc)___________________ 4.9212E-02 + Cable_WP_fraction_______________________________________________________ ((awpc-aswp-aiwp)/awpc)_______ 5.2610E-01 Turn_parametrisation____________________________________________________ (i_tf_turns_integer)__________ 0 - Number_of_turns_per_TF_coil_____________________________________________ (n_tf_turn)___________________ 1.5709E+02 OP - Width_of_turn_including_inter-turn_insulation_(m)_______________________ (t_turn_tf)___________________ 5.7439E-02 OP - Width_of_conductor_(square)_(m)_________________________________________ (t_conductor)_________________ 5.5839E-02 OP - Width_of_space_inside_conductor_(m)_____________________________________ (t_cable)_____________________ 3.9839E-02 OP - Steel_conduit_thickness_(m)_____________________________________________ (thwcndut)____________________ 8.0000E-03 ITV + Number_of_turns_per_TF_coil_____________________________________________ (n_tf_turn)___________________ 1.5438E+02 OP + Width_of_turn_including_inter-turn_insulation_(m)_______________________ (t_turn_tf)___________________ 5.7627E-02 OP + Width_of_conductor_(square)_(m)_________________________________________ (t_conductor)_________________ 5.6027E-02 OP + Width_of_space_inside_conductor_(m)_____________________________________ (t_cable)_____________________ 4.0008E-02 OP + Steel_conduit_thickness_(m)_____________________________________________ (thwcndut)____________________ 8.0091E-03 ITV Inter-turn_insulation_thickness_(m)_____________________________________ (thicndut)____________________ 8.0000E-04 Diameter_of_central_helium_channel_in_cable_____________________________ (dhecoil)_____________________ 1.0000E-02 - internal_area_of_the_cable_space________________________________________ (acstf)_______________________ 1.5562E-03 + internal_area_of_the_cable_space________________________________________ (acstf)_______________________ 1.5697E-03 Coolant_fraction_in_conductor_excluding_central_channel_________________ (vftf)________________________ 3.0000E-01 - Copper_fraction_of_conductor____________________________________________ (fcutfsu)_____________________ 8.6211E-01 ITV - Superconductor_fraction_of_conductor____________________________________ (1-fcutfsu)___________________ 1.3789E-01 + Copper_fraction_of_conductor____________________________________________ (fcutfsu)_____________________ 8.4176E-01 ITV + Superconductor_fraction_of_conductor____________________________________ (1-fcutfsu)___________________ 1.5824E-01 Check_total_area_fractions_in_winding_pack_=_1__________________________ ______________________________ 1.0000E+00 minimum_TF_conductor_temperature_margin__(K)____________________________ (tmargmin_tf)_________________ 1.5000E+00 - TF_conductor_temperature_margin_(K)_____________________________________ (tmargtf)_____________________ 1.5100E+00 + TF_conductor_temperature_margin_(K)_____________________________________ (tmargtf)_____________________ 1.5246E+00 Elastic_properties_behavior_____________________________________________ (i_tf_cond_eyoung_axial)______ 0 Conductor_axial_Youngs_modulus__________________________________________ (eyoung_cond_axial)___________ 0.0000E+00 Conductor_transverse_Youngs_modulus_____________________________________ (eyoung_cond_trans)___________ 0.0000E+00 - Superconductor_mass_per_coil_(kg)_______________________________________ (whtconsc)____________________ 3.2728E+03 OP - Copper_mass_per_coil_(kg)_______________________________________________ (whtconcu)____________________ 5.7530E+04 OP - Steel_conduit_mass_per_coil_(kg)________________________________________ (whtconsh)____________________ 9.1497E+04 OP - Conduit_insulation_mass_per_coil_(kg)___________________________________ (whtconin)____________________ 2.4504E+03 OP - Total_conduit_mass_per_coil_(kg)________________________________________ (whtcon)______________________ 1.5475E+05 OP - Mass_of_each_TF_coil_(kg)_______________________________________________ (whttf/n_tf)__________________ 6.5696E+05 OP - Total_TF_coil_mass_(kg)_________________________________________________ (whttf)_______________________ 1.0511E+07 OP - Nominal_peak_field_assuming_toroidal_symmetry_(T)_______________________ (bmaxtf)______________________ 1.1646E+01 OP - Total_current_in_all_TF_coils_(MA)______________________________________ (ritfc/1.D6)__________________ 2.0251E+02 OP - TF_coil_current_(summed_over_all_coils)_(A)_____________________________ (ritfc)_______________________ 2.0251E+08 - Actual_peak_field_at_discrete_conductor_(T)_____________________________ (bmaxtfrp)____________________ 1.2085E+01 OP - Winding_pack_current_density_(A/m2)_____________________________________ (jwptf)_______________________ 2.4422E+07 OP - Inboard_leg_mid-plane_conductor_current_density_(A/m2)__________________ (oacdcp)______________________ 1.1672E+07 - Total_stored_energy_in_TF_coils_(GJ)____________________________________ (estotftgj)___________________ 1.1034E+02 OP - Inboard_vertical_tension_per_coil_(N)___________________________________ (vforce)______________________ 1.9146E+08 OP - Outboard_vertical_tension_per_coil_(N)__________________________________ (vforce_outboard)_____________ 1.9146E+08 OP + Superconductor_mass_per_coil_(kg)_______________________________________ (whtconsc)____________________ 4.2854E+03 OP + Copper_mass_per_coil_(kg)_______________________________________________ (whtconcu)____________________ 5.5694E+04 OP + Steel_conduit_mass_per_coil_(kg)________________________________________ (whtconsh)____________________ 9.0500E+04 OP + Conduit_insulation_mass_per_coil_(kg)___________________________________ (whtconin)____________________ 2.4201E+03 OP + Total_conduit_mass_per_coil_(kg)________________________________________ (whtcon)______________________ 1.5290E+05 OP + Mass_of_each_TF_coil_(kg)_______________________________________________ (whttf/n_tf)__________________ 6.9508E+05 OP + Total_TF_coil_mass_(kg)_________________________________________________ (whttf)_______________________ 1.1121E+07 OP + Nominal_peak_field_assuming_toroidal_symmetry_(T)_______________________ (bmaxtf)______________________ 1.2048E+01 OP + Total_current_in_all_TF_coils_(MA)______________________________________ (ritfc/1.D6)__________________ 2.0950E+02 OP + TF_coil_current_(summed_over_all_coils)_(A)_____________________________ (ritfc)_______________________ 2.0950E+08 + Actual_peak_field_at_discrete_conductor_(T)_____________________________ (bmaxtfrp)____________________ 1.2501E+01 OP + Winding_pack_current_density_(A/m2)_____________________________________ (jwptf)_______________________ 2.5540E+07 OP + Inboard_leg_mid-plane_conductor_current_density_(A/m2)__________________ (oacdcp)______________________ 1.1703E+07 + Total_stored_energy_in_TF_coils_(GJ)____________________________________ (estotftgj)___________________ 1.1862E+02 OP + Inboard_vertical_tension_per_coil_(N)___________________________________ (vforce)______________________ 2.0466E+08 OP + Outboard_vertical_tension_per_coil_(N)__________________________________ (vforce_outboard)_____________ 2.0466E+08 OP inboard_vertical_tension_fraction_(-)___________________________________ (f_vforce_inboard)____________ 5.0000E-01 OP - Centring_force_per_coil_(N/m)___________________________________________ (cforce)______________________ 7.3703E+07 OP + Centring_force_per_coil_(N/m)___________________________________________ (cforce)______________________ 7.8881E+07 OP Max_allowed_tfcoil_variables.ripple_amplitude_at_plasma_outboard_midplan (ripmax)______________________ 6.0000E-01 Ripple_amplitude_at_plasma_outboard_midplane_(%)________________________ (ripple)______________________ 6.0000E-01 OP - Allowable_stress_in_vacuum_vessel_(VV)_due_to_quench_(Pa)_______________ (sigvvall)____________________ 1.0000E+08 - Minimum_allowed_quench_time_due_to_stress_in_VV_(s)_____________________ (taucq)_______________________ 1.8023E+01 OP - Actual_quench_time_(or_time_constant)_(s)_______________________________ (tdmptf)______________________ 1.9607E+01 ITV - Vacuum_Vassel_stress_on_quench_(Pa)_____________________________________ (vv_stress_quench)____________ 5.5495E+07 OP - Maximum_allowed_voltage_during_quench_due_to_insulation_(kV)____________ (vdalw)_______________________ 9.3438E+00 ITV - Actual_quench_voltage_(kV)______________________________________________ (vtfskv)______________________ 8.7306E+00 OP + Actual_quench_time_(or_time_constant)_(s)_______________________________ (tdmptf)______________________ 1.7534E+01 ITV + Vacuum_Vessel_stress_on_quench_(Pa)_____________________________________ (vv_stress_quench)____________ 6.2719E+07 OP + Maximum_allowed_voltage_during_quench_due_to_insulation_(kV)____________ (vdalw)_______________________ 9.9851E+00 ITV + Actual_quench_voltage_(kV)______________________________________________ (vtfskv)______________________ 9.9703E+00 OP Maximum_allowed_temp_rise_during_a_quench_(K)___________________________ (tmaxpro)_____________________ 1.5000E+02 # Superconducting TF Coils # Superconductor_switch___________________________________________________ (isumat)______________________ 1 Critical_field_at_zero_temperature_and_strain_(T)_______________________ (bc20m)_______________________ 3.2970E+01 Critical_temperature_at_zero_field_and_strain_(K)_______________________ (tc0m)________________________ 1.6060E+01 Helium_temperature_at_peak_field_(=_superconductor_temperature)_(K)_____ (thelium)_____________________ 4.7500E+00 - Total_helium_fraction_inside_cable_space________________________________ (fhetot)______________________ 3.5047E-01 OP - Copper_fraction_of_conductor____________________________________________ (fcutfsu)_____________________ 8.6211E-01 ITV + Total_helium_fraction_inside_cable_space________________________________ (fhetot)______________________ 3.5003E-01 OP + Copper_fraction_of_conductor____________________________________________ (fcutfsu)_____________________ 8.4176E-01 ITV Residual_manufacturing_strain_on_superconductor_________________________ (str_tf_con_res)______________ -5.0000E-03 - Self-consistent_strain_on_superconductor________________________________ (str_wp)______________________ 1.9975E-03 - Critical_current_density_in_superconductor_(A/m2)_______________________ (jcritsc)_____________________ 7.9100E+08 OP - Critical_current_density_in_strand_(A/m2)_______________________________ (jcritstr)____________________ 1.0907E+08 OP - Critical_current_density_in_winding_pack_(A/m2)_________________________ (jwdgcrt)_____________________ 3.3416E+07 OP - Actual_current_density_in_winding_pack_(A/m2)___________________________ (jwdgop)______________________ 2.4422E+07 OP + Self-consistent_strain_on_superconductor________________________________ (str_wp)______________________ 2.0475E-03 + Critical_current_density_in_superconductor_(A/m2)_______________________ (jcritsc)_____________________ 7.3233E+08 OP + Critical_current_density_in_strand_(A/m2)_______________________________ (jcritstr)____________________ 1.1588E+08 OP + Critical_current_density_in_winding_pack_(A/m2)_________________________ (jwdgcrt)_____________________ 3.5602E+07 OP + Actual_current_density_in_winding_pack_(A/m2)___________________________ (jwdgop)______________________ 2.5540E+07 OP Minimum_allowed_temperature_margin_in_superconductor_(K)________________ (tmargmin_tf)_________________ 1.5000E+00 - Actual_temperature_margin_in_superconductor_(K)_________________________ (tmarg)_______________________ 1.5100E+00 OP - Critical_current_(A)____________________________________________________ (icrit)_______________________ 1.1025E+05 OP - Actual_current_(A)______________________________________________________ (cpttf)_______________________ 8.0573E+04 ITV - Actual_current_/_critical_current_______________________________________ (iooic)_______________________ 7.3084E-01 OP + Actual_temperature_margin_in_superconductor_(K)_________________________ (tmarg)_______________________ 1.5246E+00 OP + Critical_current_(A)____________________________________________________ (icrit)_______________________ 1.1823E+05 OP + Actual_current_(A)______________________________________________________ (cpttf)_______________________ 8.4815E+04 ITV + Actual_current_/_critical_current_______________________________________ (iooic)_______________________ 7.1738E-01 OP # Central Solenoid and PF Coils # Central_solenoid_superconductor_material________________________________ (isumatoh)____________________ 1 - Maximum_field_at_Beginning_Of_Pulse_(T)_________________________________ (bmaxoh0)_____________________ 1.3997E+01 OP - Critical_superconductor_current_density_at_BOP_(A/m2)___________________ (jscoh_bop)___________________ 3.5395E+08 OP - Critical_strand_current_density_at_BOP_(A/m2)___________________________ (jstrandoh_bop)_______________ 1.0618E+08 OP - Allowable_overall_current_density_at_BOP_(A/m2)_________________________ (rjohc0)______________________ 3.4704E+07 OP - Actual_overall_current_density_at_BOP_(A/m2)____________________________ (cohbop)______________________ 1.9770E+07 OP - Maximum_field_at_End_Of_Flattop_(T)_____________________________________ (bmaxoh)______________________ 1.3221E+01 OP - Critical_superconductor_current_density_at_EOF_(A/m2)___________________ (jscoh_eof)___________________ 4.2900E+08 OP - Critical_strand_current_density_at_EOF_(A/m2)___________________________ (jstrandoh_eof)_______________ 1.2870E+08 OP - Allowable_overall_current_density_at_EOF_(A/m2)_________________________ (rjohc)_______________________ 4.2063E+07 OP - Actual_overall_current_density_at_EOF_(A/m2)____________________________ (coheof)______________________ 2.0066E+07 ITV - CS_inside_radius_(m)____________________________________________________ (bore)________________________ 2.0588E+00 ITV - CS_thickness_(m)________________________________________________________ (ohcth)_______________________ 5.5551E-01 ITV + Maximum_field_at_Beginning_Of_Pulse_(T)_________________________________ (bmaxoh0)_____________________ 1.4017E+01 OP + Critical_superconductor_current_density_at_BOP_(A/m2)___________________ (jscoh_bop)___________________ 3.5219E+08 OP + Critical_strand_current_density_at_BOP_(A/m2)___________________________ (jstrandoh_bop)_______________ 1.0566E+08 OP + Allowable_overall_current_density_at_BOP_(A/m2)_________________________ (rjohc0)______________________ 3.5378E+07 OP + Actual_overall_current_density_at_BOP_(A/m2)____________________________ (cohbop)______________________ 1.9764E+07 OP + Maximum_field_at_End_Of_Flattop_(T)_____________________________________ (bmaxoh)______________________ 1.3500E+01 OP + Critical_superconductor_current_density_at_EOF_(A/m2)___________________ (jscoh_eof)___________________ 4.0090E+08 OP + Critical_strand_current_density_at_EOF_(A/m2)___________________________ (jstrandoh_eof)_______________ 1.2027E+08 OP + Allowable_overall_current_density_at_EOF_(A/m2)_________________________ (rjohc)_______________________ 4.0271E+07 OP + Actual_overall_current_density_at_EOF_(A/m2)____________________________ (coheof)______________________ 2.0403E+07 ITV + CS_inside_radius_(m)____________________________________________________ (bore)________________________ 2.0251E+00 ITV + CS_thickness_(m)________________________________________________________ (ohcth)_______________________ 5.5625E-01 ITV Gap_between_central_solenoid_and_TF_coil_(m)____________________________ (gapoh)_______________________ 8.0000E-02 - CS_overall_cross-sectional_area_(m2)____________________________________ (areaoh)______________________ 8.8176E+00 OP - CS_conductor+void_cross-sectional_area_(m2)_____________________________ (awpoh)_______________________ 4.1169E+00 OP - ___CS_conductor_cross-sectional_area_(m2)_______________________________ (awpoh*(1-vfohc))_____________ 2.8819E+00 OP - ___CS_void_cross-sectional_area_(m2)____________________________________ (awpoh*vfohc)_________________ 1.2351E+00 OP - CS_steel_cross-sectional_area_(m2)______________________________________ (areaoh-awpoh)________________ 4.7006E+00 OP - CS_steel_area_fraction__________________________________________________ (oh_steel_frac)_______________ 5.3310E-01 ITV + CS_overall_cross-sectional_area_(m2)____________________________________ (areaoh)______________________ 8.8292E+00 OP + CS_conductor+void_cross-sectional_area_(m2)_____________________________ (awpoh)_______________________ 4.2234E+00 OP + ___CS_conductor_cross-sectional_area_(m2)_______________________________ (awpoh*(1-vfohc))_____________ 2.9563E+00 OP + ___CS_void_cross-sectional_area_(m2)____________________________________ (awpoh*vfohc)_________________ 1.2670E+00 OP + CS_steel_cross-sectional_area_(m2)______________________________________ (areaoh-awpoh)________________ 4.6058E+00 OP + CS_steel_area_fraction__________________________________________________ (oh_steel_frac)_______________ 5.2166E-01 ITV Switch_for_CS_stress_calculation________________________________________ (i_cs_stress)_________________ 0 Allowable_stress_in_CS_steel_(Pa)_______________________________________ (alstroh)_____________________ 7.5000E+08 - Hoop_stress_in_CS_steel_(Pa)____________________________________________ (sig_hoop)____________________ 6.8991E+08 OP - Axial_stress_in_CS_steel_(Pa)___________________________________________ (sig_axial)___________________ -7.0866E+08 OP - Maximum_shear_stress_in_CS_steel_for_the_Tresca_criterion_(Pa)__________ (s_tresca_oh)_________________ 6.8991E+08 OP - Axial_force_in_CS_(N)___________________________________________________ (axial_force)_________________ -1.5405E+09 OP + Hoop_stress_in_CS_steel_(Pa)____________________________________________ (sig_hoop)____________________ 6.9716E+08 OP + Axial_stress_in_CS_steel_(Pa)___________________________________________ (sig_axial)___________________ -7.4296E+08 OP + Maximum_shear_stress_in_CS_steel_for_the_Tresca_criterion_(Pa)__________ (s_tresca_oh)_________________ 6.9716E+08 OP + Axial_force_in_CS_(N)___________________________________________________ (axial_force)_________________ -1.5599E+09 OP Residual_manufacturing_strain_in_CS_superconductor_material_____________ (tfcoil_variables.str_cs_con_r -5.0000E-03 Copper_fraction_in_strand_______________________________________________ (fcuohsu)_____________________ 7.0000E-01 Void_(coolant)_fraction_in_conductor____________________________________ (vfohc)_______________________ 3.0000E-01 Helium_coolant_temperature_(K)__________________________________________ (tftmp)_______________________ 4.7500E+00 - CS_temperature_margin_(K)_______________________________________________ (tmargoh)_____________________ 1.5000E+00 OP + CS_temperature_margin_(K)_______________________________________________ (tmargoh)_____________________ 1.5020E+00 OP Minimum_permitted_temperature_margin_(K)________________________________ (tmargmin_cs)_________________ 1.5000E+00 Residual_hoop_stress_in_CS_Steel_(Pa)___________________________________ (residual_sig_hoop)___________ 2.4000E+08 Minimum_burn_time_(s)___________________________________________________ (tbrnmn)______________________ 7.2000E+03 Initial_vertical_crack_size_(m)_________________________________________ (t_crack_vertical)____________ 8.9000E-04 Initial_radial_crack_size_(m)___________________________________________ (t_crack_radial)______________ 2.6700E-03 - CS_turn_area_(m)________________________________________________________ (a_oh_turn)___________________ 1.9935E-03 - CS_turn_length_(m)______________________________________________________ (l_cond_cst)__________________ 7.7333E-02 - CS_turn_internal_cable_space_radius_(m)_________________________________ (r_in_cst)____________________ 7.3195E-03 - CS_turn_width_(m)_______________________________________________________ (d_cond_cst)__________________ 2.5778E-02 - CS_structural_vertical_thickness_(m)____________________________________ (t_structural_vertical)_______ 5.5693E-03 - CS_structural_radial_thickness_(m)______________________________________ (t_structural_radial)_________ 5.5693E-03 - Allowable_number_of_cycles_till_CS_fracture_____________________________ (n_cycle)_____________________ 4.2524E+02 OP + CS_turn_area_(m)________________________________________________________ (a_oh_turn)___________________ 1.9605E-03 + CS_turn_length_(m)______________________________________________________ (l_cond_cst)__________________ 7.6691E-02 + CS_turn_internal_cable_space_radius_(m)_________________________________ (r_in_cst)____________________ 7.4090E-03 + CS_turn_width_(m)_______________________________________________________ (d_cond_cst)__________________ 2.5564E-02 + CS_structural_vertical_thickness_(m)____________________________________ (t_structural_vertical)_______ 5.3728E-03 + CS_structural_radial_thickness_(m)______________________________________ (t_structural_radial)_________ 5.3728E-03 + Allowable_number_of_cycles_till_CS_fracture_____________________________ (n_cycle)_____________________ 1.2206E+02 OP Minimum_number_of_cycles_required_till_CS_fracture______________________ (n_cycle_min)_________________ 2.0000E+04 OP PF_coil_superconductor_material_________________________________________ (isumatpf)____________________ 3 Copper_fraction_in_conductor____________________________________________ (fcupfsu)_____________________ 6.9000E-01 Maximum_permissible_tensile_stress_(MPa)________________________________ (sigpfcalw)___________________ 5.0000E+02 JxB_hoop_force_fraction_supported_by_case_______________________________ (sigpfcf)_____________________ 6.6600E-01 PF_coil_0_radius_(m)____________________________________________________ (rpf[0]_______________________ 5.5667E+00 - PF_coil_0_vertical_position_(m)_________________________________________ (zpf[0])______________________ 9.3091E+00 - PF_coil_0_radial_thickness_(m)__________________________________________ (pfdr(0))_____________________ 1.2994E+00 - PF_coil_0_vertical_thickness_(m)________________________________________ (pfdz(0))_____________________ 1.2994E+00 - PF_coil_0_turns_________________________________________________________ (turns[0])____________________ 4.6436E+02 - PF_coil_0_current_(MA)__________________________________________________ (ric[0])______________________ 1.8574E+01 - PF_coil_0_field_(T)_____________________________________________________ (bpf[0])______________________ 6.6228E+00 + PF_coil_0_vertical_position_(m)_________________________________________ (zpf[0])______________________ 9.3411E+00 + PF_coil_0_radial_thickness_(m)__________________________________________ (pfdr(0))_____________________ 1.2855E+00 + PF_coil_0_vertical_thickness_(m)________________________________________ (pfdz(0))_____________________ 1.2855E+00 + PF_coil_0_turns_________________________________________________________ (turns[0])____________________ 4.5442E+02 + PF_coil_0_current_(MA)__________________________________________________ (ric[0])______________________ 1.8177E+01 + PF_coil_0_field_(T)_____________________________________________________ (bpf[0])______________________ 6.4291E+00 PF_coil_1_radius_(m)____________________________________________________ (rpf[1]_______________________ 5.5667E+00 - PF_coil_1_vertical_position_(m)_________________________________________ (zpf[1])______________________ -1.0543E+01 - PF_coil_1_radial_thickness_(m)__________________________________________ (pfdr(1))_____________________ 1.3991E+00 - PF_coil_1_vertical_thickness_(m)________________________________________ (pfdz(1))_____________________ -1.3991E+00 - PF_coil_1_turns_________________________________________________________ (turns[1])____________________ 5.3830E+02 - PF_coil_1_current_(MA)__________________________________________________ (ric[1])______________________ 2.1532E+01 - PF_coil_1_field_(T)_____________________________________________________ (bpf[1])______________________ 7.1259E+00 - PF_coil_2_radius_(m)____________________________________________________ (rpf[2]_______________________ 1.6706E+01 + PF_coil_1_vertical_position_(m)_________________________________________ (zpf[1])______________________ -1.0575E+01 + PF_coil_1_radial_thickness_(m)__________________________________________ (pfdr(1))_____________________ 1.3815E+00 + PF_coil_1_vertical_thickness_(m)________________________________________ (pfdz(1))_____________________ -1.3815E+00 + PF_coil_1_turns_________________________________________________________ (turns[1])____________________ 5.2482E+02 + PF_coil_1_current_(MA)__________________________________________________ (ric[1])______________________ 2.0993E+01 + PF_coil_1_field_(T)_____________________________________________________ (bpf[1])______________________ 6.9155E+00 + PF_coil_2_radius_(m)____________________________________________________ (rpf[2]_______________________ 1.6722E+01 PF_coil_2_vertical_position_(m)_________________________________________ (zpf[2])______________________ 2.6667E+00 - PF_coil_2_radial_thickness_(m)__________________________________________ (pfdr(2))_____________________ 1.1273E+00 - PF_coil_2_vertical_thickness_(m)________________________________________ (pfdz(2))_____________________ 1.1273E+00 - PF_coil_2_turns_________________________________________________________ (turns[2])____________________ 1.9061E+02 - PF_coil_2_current_(MA)__________________________________________________ (ric[2])______________________ -7.6244E+00 - PF_coil_2_field_(T)_____________________________________________________ (bpf[2])______________________ 2.6246E+00 - PF_coil_3_radius_(m)____________________________________________________ (rpf[3]_______________________ 1.6706E+01 + PF_coil_2_radial_thickness_(m)__________________________________________ (pfdr(2))_____________________ 1.1305E+00 + PF_coil_2_vertical_thickness_(m)________________________________________ (pfdz(2))_____________________ 1.1305E+00 + PF_coil_2_turns_________________________________________________________ (turns[2])____________________ 1.9172E+02 + PF_coil_2_current_(MA)__________________________________________________ (ric[2])______________________ -7.6688E+00 + PF_coil_2_field_(T)_____________________________________________________ (bpf[2])______________________ 2.6325E+00 + PF_coil_3_radius_(m)____________________________________________________ (rpf[3]_______________________ 1.6722E+01 PF_coil_3_vertical_position_(m)_________________________________________ (zpf[3])______________________ -2.6667E+00 - PF_coil_3_radial_thickness_(m)__________________________________________ (pfdr(3))_____________________ 1.1273E+00 - PF_coil_3_vertical_thickness_(m)________________________________________ (pfdz(3))_____________________ -1.1273E+00 - PF_coil_3_turns_________________________________________________________ (turns[3])____________________ 1.9061E+02 - PF_coil_3_current_(MA)__________________________________________________ (ric[3])______________________ -7.6244E+00 - PF_coil_3_field_(T)_____________________________________________________ (bpf[3])______________________ 2.6246E+00 - PF_coil_4_radius_(m)____________________________________________________ (rpf[4]_______________________ 1.5180E+01 + PF_coil_3_radial_thickness_(m)__________________________________________ (pfdr(3))_____________________ 1.1305E+00 + PF_coil_3_vertical_thickness_(m)________________________________________ (pfdz(3))_____________________ -1.1305E+00 + PF_coil_3_turns_________________________________________________________ (turns[3])____________________ 1.9172E+02 + PF_coil_3_current_(MA)__________________________________________________ (ric[3])______________________ -7.6688E+00 + PF_coil_3_field_(T)_____________________________________________________ (bpf[3])______________________ 2.6325E+00 + PF_coil_4_radius_(m)____________________________________________________ (rpf[4]_______________________ 1.5198E+01 PF_coil_4_vertical_position_(m)_________________________________________ (zpf[4])______________________ 7.4667E+00 - PF_coil_4_radial_thickness_(m)__________________________________________ (pfdr(4))_____________________ 8.0922E-01 - PF_coil_4_vertical_thickness_(m)________________________________________ (pfdz(4))_____________________ 8.0922E-01 - PF_coil_4_turns_________________________________________________________ (turns[4])____________________ 1.3097E+02 - PF_coil_4_current_(MA)__________________________________________________ (ric[4])______________________ -5.2388E+00 - PF_coil_4_field_(T)_____________________________________________________ (bpf[4])______________________ 2.5758E+00 - PF_coil_5_radius_(m)____________________________________________________ (rpf[5]_______________________ 1.5180E+01 + PF_coil_4_radial_thickness_(m)__________________________________________ (pfdr(4))_____________________ 8.0816E-01 + PF_coil_4_vertical_thickness_(m)________________________________________ (pfdz(4))_____________________ 8.0816E-01 + PF_coil_4_turns_________________________________________________________ (turns[4])____________________ 1.3063E+02 + PF_coil_4_current_(MA)__________________________________________________ (ric[4])______________________ -5.2250E+00 + PF_coil_4_field_(T)_____________________________________________________ (bpf[4])______________________ 2.5744E+00 + PF_coil_5_radius_(m)____________________________________________________ (rpf[5]_______________________ 1.5198E+01 PF_coil_5_vertical_position_(m)_________________________________________ (zpf[5])______________________ -7.4667E+00 - PF_coil_5_radial_thickness_(m)__________________________________________ (pfdr(5))_____________________ 8.0922E-01 - PF_coil_5_vertical_thickness_(m)________________________________________ (pfdz(5))_____________________ -8.0922E-01 - PF_coil_5_turns_________________________________________________________ (turns[5])____________________ 1.3097E+02 - PF_coil_5_current_(MA)__________________________________________________ (ric[5])______________________ -5.2388E+00 - PF_coil_5_field_(T)_____________________________________________________ (bpf[5])______________________ 2.5758E+00 - Central_solenoid_radius_(m)_____________________________________________ (rpf[nohc-1])_________________ 2.3365E+00 + PF_coil_5_radial_thickness_(m)__________________________________________ (pfdr(5))_____________________ 8.0816E-01 + PF_coil_5_vertical_thickness_(m)________________________________________ (pfdz(5))_____________________ -8.0816E-01 + PF_coil_5_turns_________________________________________________________ (turns[5])____________________ 1.3063E+02 + PF_coil_5_current_(MA)__________________________________________________ (ric[5])______________________ -5.2250E+00 + PF_coil_5_field_(T)_____________________________________________________ (bpf[5])______________________ 2.5744E+00 + Central_solenoid_radius_(m)_____________________________________________ (rpf[nohc-1])_________________ 2.3032E+00 Central_solenoid_vertical_position_(m)__________________________________ (zpf[nohc-1])_________________ 0.0000E+00 - Central_solenoid_radial_thickness_(m)___________________________________ (ohdr)________________________ 5.5551E-01 + Central_solenoid_radial_thickness_(m)___________________________________ (ohdr)________________________ 5.5625E-01 Central_solenoid_vertical_thickness_(m)_________________________________ (ohdz)________________________ 1.5873E+01 - Central_solenoid_turns__________________________________________________ (turns[nohc-1])_______________ 4.4232E+03 - Central_solenoid_current_(MA)___________________________________________ (ric[nohc-1])_________________ -1.7693E+02 - Central_solenoid_field_(T)______________________________________________ (bpf[nohc-1])_________________ 1.3997E+01 - Sum_of_squares_of_residuals_____________________________________________ (ssq0)________________________ 3.8915E-04 OP + Central_solenoid_turns__________________________________________________ (turns[nohc-1])_______________ 4.5036E+03 + Central_solenoid_current_(MA)___________________________________________ (ric[nohc-1])_________________ -1.8014E+02 + Central_solenoid_field_(T)______________________________________________ (bpf[nohc-1])_________________ 1.4017E+01 + Sum_of_squares_of_residuals_____________________________________________ (ssq0)________________________ 3.6878E-04 OP Smoothing_parameter_____________________________________________________ (alfapf)______________________ 5.0000E-10 # Volt Second Consumption # - Total_volt-second_consumption_by_coils_(Wb)_____________________________ (vstot)_______________________ -5.9000E+02 OP + Total_volt-second_consumption_by_coils_(Wb)_____________________________ (vstot)_______________________ -5.8000E+02 OP # Waveforms # - Ratio_of_central_solenoid_current_at_beginning_of_Pulse_/_end_of_flat-to (fcohbop)_____________________ 9.8528E-01 ITV - Ratio_of_central_solenoid_current_at_beginning_of_Flat-top_/_end_of_flat (fcohbof)_____________________ -1.9610E-01 OP + Ratio_of_central_solenoid_current_at_beginning_of_Pulse_/_end_of_flat-to (fcohbop)_____________________ 9.6866E-01 ITV + Ratio_of_central_solenoid_current_at_beginning_of_Flat-top_/_end_of_flat (fcohbof)_____________________ -1.7009E-01 OP # PF Circuit Waveform Data # Number_of_PF_circuits_including_CS_and_plasma___________________________ (ncirt)_______________________ 8 PF_Circuit_0_Time_point_0_(A)___________________________________________ (pfc0t0)______________________ 0.0000E+00 - PF_Circuit_0_Time_point_1_(A)___________________________________________ (pfc0t1)______________________ 1.4447E+07 - PF_Circuit_0_Time_point_2_(A)___________________________________________ (pfc0t2)______________________ 1.8574E+07 - PF_Circuit_0_Time_point_3_(A)___________________________________________ (pfc0t3)______________________ 1.8574E+07 - PF_Circuit_0_Time_point_4_(A)___________________________________________ (pfc0t4)______________________ 1.0356E+06 + PF_Circuit_0_Time_point_1_(A)___________________________________________ (pfc0t1)______________________ 1.4212E+07 + PF_Circuit_0_Time_point_2_(A)___________________________________________ (pfc0t2)______________________ 1.8177E+07 + PF_Circuit_0_Time_point_3_(A)___________________________________________ (pfc0t3)______________________ 1.8177E+07 + PF_Circuit_0_Time_point_4_(A)___________________________________________ (pfc0t4)______________________ 1.0098E+06 PF_Circuit_0_Time_point_5_(A)___________________________________________ (pfc0t5)______________________ 0.0000E+00 PF_Circuit_1_Time_point_0_(A)___________________________________________ (pfc1t0)______________________ 0.0000E+00 - PF_Circuit_1_Time_point_1_(A)___________________________________________ (pfc1t1)______________________ 1.8901E+07 - PF_Circuit_1_Time_point_2_(A)___________________________________________ (pfc1t2)______________________ 2.1532E+07 - PF_Circuit_1_Time_point_3_(A)___________________________________________ (pfc1t3)______________________ 2.1532E+07 - PF_Circuit_1_Time_point_4_(A)___________________________________________ (pfc1t4)______________________ -1.4128E+06 + PF_Circuit_1_Time_point_1_(A)___________________________________________ (pfc1t1)______________________ 1.8587E+07 + PF_Circuit_1_Time_point_2_(A)___________________________________________ (pfc1t2)______________________ 2.0993E+07 + PF_Circuit_1_Time_point_3_(A)___________________________________________ (pfc1t3)______________________ 2.0993E+07 + PF_Circuit_1_Time_point_4_(A)___________________________________________ (pfc1t4)______________________ -1.4588E+06 PF_Circuit_1_Time_point_5_(A)___________________________________________ (pfc1t5)______________________ 0.0000E+00 PF_Circuit_2_Time_point_0_(A)___________________________________________ (pfc2t0)______________________ -0.0000E+00 - PF_Circuit_2_Time_point_1_(A)___________________________________________ (pfc2t1)______________________ 4.9967E+05 - PF_Circuit_2_Time_point_2_(A)___________________________________________ (pfc2t2)______________________ -7.0179E+06 - PF_Circuit_2_Time_point_3_(A)___________________________________________ (pfc2t3)______________________ -7.0179E+06 - PF_Circuit_2_Time_point_4_(A)___________________________________________ (pfc2t4)______________________ -7.6244E+06 + PF_Circuit_2_Time_point_1_(A)___________________________________________ (pfc2t1)______________________ 4.9864E+05 + PF_Circuit_2_Time_point_2_(A)___________________________________________ (pfc2t2)______________________ -7.0665E+06 + PF_Circuit_2_Time_point_3_(A)___________________________________________ (pfc2t3)______________________ -7.0665E+06 + PF_Circuit_2_Time_point_4_(A)___________________________________________ (pfc2t4)______________________ -7.6688E+06 PF_Circuit_2_Time_point_5_(A)___________________________________________ (pfc2t5)______________________ -0.0000E+00 PF_Circuit_3_Time_point_0_(A)___________________________________________ (pfc3t0)______________________ -0.0000E+00 - PF_Circuit_3_Time_point_1_(A)___________________________________________ (pfc3t1)______________________ 4.9967E+05 - PF_Circuit_3_Time_point_2_(A)___________________________________________ (pfc3t2)______________________ -7.0179E+06 - PF_Circuit_3_Time_point_3_(A)___________________________________________ (pfc3t3)______________________ -7.0179E+06 - PF_Circuit_3_Time_point_4_(A)___________________________________________ (pfc3t4)______________________ -7.6244E+06 + PF_Circuit_3_Time_point_1_(A)___________________________________________ (pfc3t1)______________________ 4.9864E+05 + PF_Circuit_3_Time_point_2_(A)___________________________________________ (pfc3t2)______________________ -7.0665E+06 + PF_Circuit_3_Time_point_3_(A)___________________________________________ (pfc3t3)______________________ -7.0665E+06 + PF_Circuit_3_Time_point_4_(A)___________________________________________ (pfc3t4)______________________ -7.6688E+06 PF_Circuit_3_Time_point_5_(A)___________________________________________ (pfc3t5)______________________ -0.0000E+00 PF_Circuit_4_Time_point_0_(A)___________________________________________ (pfc4t0)______________________ -0.0000E+00 - PF_Circuit_4_Time_point_1_(A)___________________________________________ (pfc4t1)______________________ 4.1114E+05 - PF_Circuit_4_Time_point_2_(A)___________________________________________ (pfc4t2)______________________ -4.7396E+06 - PF_Circuit_4_Time_point_3_(A)___________________________________________ (pfc4t3)______________________ -4.7396E+06 - PF_Circuit_4_Time_point_4_(A)___________________________________________ (pfc4t4)______________________ -5.2388E+06 + PF_Circuit_4_Time_point_1_(A)___________________________________________ (pfc4t1)______________________ 3.6106E+05 + PF_Circuit_4_Time_point_2_(A)___________________________________________ (pfc4t2)______________________ -4.7889E+06 + PF_Circuit_4_Time_point_3_(A)___________________________________________ (pfc4t3)______________________ -4.7889E+06 + PF_Circuit_4_Time_point_4_(A)___________________________________________ (pfc4t4)______________________ -5.2250E+06 PF_Circuit_4_Time_point_5_(A)___________________________________________ (pfc4t5)______________________ -0.0000E+00 PF_Circuit_5_Time_point_0_(A)___________________________________________ (pfc5t0)______________________ -0.0000E+00 - PF_Circuit_5_Time_point_1_(A)___________________________________________ (pfc5t1)______________________ 4.1114E+05 - PF_Circuit_5_Time_point_2_(A)___________________________________________ (pfc5t2)______________________ -4.7396E+06 - PF_Circuit_5_Time_point_3_(A)___________________________________________ (pfc5t3)______________________ -4.7396E+06 - PF_Circuit_5_Time_point_4_(A)___________________________________________ (pfc5t4)______________________ -5.2388E+06 + PF_Circuit_5_Time_point_1_(A)___________________________________________ (pfc5t1)______________________ 3.6106E+05 + PF_Circuit_5_Time_point_2_(A)___________________________________________ (pfc5t2)______________________ -4.7889E+06 + PF_Circuit_5_Time_point_3_(A)___________________________________________ (pfc5t3)______________________ -4.7889E+06 + PF_Circuit_5_Time_point_4_(A)___________________________________________ (pfc5t4)______________________ -5.2250E+06 PF_Circuit_5_Time_point_5_(A)___________________________________________ (pfc5t5)______________________ -0.0000E+00 CS_Circuit_Time_point_0_(A)_____________________________________________ (cs_t0)_______________________ -0.0000E+00 - CS_Circuit_Time_point_1_(A)_____________________________________________ (cs_t1)_______________________ 1.7433E+08 - CS_Circuit_Time_point_2_(A)_____________________________________________ (cs_t2)_______________________ 3.4696E+07 - CS_Circuit_Time_point_3_(A)_____________________________________________ (cs_t3)_______________________ 3.4696E+07 - CS_Circuit_Time_point_4_(A)_____________________________________________ (cs_t4)_______________________ -1.7693E+08 + CS_Circuit_Time_point_1_(A)_____________________________________________ (cs_t1)_______________________ 1.7450E+08 + CS_Circuit_Time_point_2_(A)_____________________________________________ (cs_t2)_______________________ 3.0640E+07 + CS_Circuit_Time_point_3_(A)_____________________________________________ (cs_t3)_______________________ 3.0640E+07 + CS_Circuit_Time_point_4_(A)_____________________________________________ (cs_t4)_______________________ -1.8014E+08 CS_Circuit_Time_point_5_(A)_____________________________________________ (cs_t5)_______________________ -0.0000E+00 Plasma_Time_point_0_(A)_________________________________________________ (plasmat0)____________________ 0.0000E+00 Plasma_Time_point_1_(A)_________________________________________________ (plasmat1)____________________ 0.0000E+00 - Plasma_Time_point_2_(A)_________________________________________________ (plasmat2)____________________ 1.6699E+07 - Plasma_Time_point_3_(A)_________________________________________________ (plasmat3)____________________ 1.6699E+07 - Plasma_Time_point_4_(A)_________________________________________________ (plasmat4)____________________ 1.6699E+07 + Plasma_Time_point_2_(A)_________________________________________________ (plasmat2)____________________ 1.6616E+07 + Plasma_Time_point_3_(A)_________________________________________________ (plasmat3)____________________ 1.6616E+07 + Plasma_Time_point_4_(A)_________________________________________________ (plasmat4)____________________ 1.6616E+07 Plasma_Time_point_5_(A)_________________________________________________ (plasmat5)____________________ 0.0000E+00 # Support Structure # - Outer_PF_coil_fence_mass_(kg)___________________________________________ (fncmass)_____________________ 2.3111E+05 OP - Intercoil_support_structure_mass_(kg)___________________________________ (aintmass)____________________ 4.0561E+06 OP - Mass_of_cooled_components_(kg)__________________________________________ (coldmass)____________________ 3.3225E+07 OP - Gravity_support_structure_mass_(kg)_____________________________________ (clgsmass)____________________ 1.2335E+06 OP + Outer_PF_coil_fence_mass_(kg)___________________________________________ (fncmass)_____________________ 2.2883E+05 OP + Intercoil_support_structure_mass_(kg)___________________________________ (aintmass)____________________ 4.2248E+06 OP + Mass_of_cooled_components_(kg)__________________________________________ (coldmass)____________________ 3.3983E+07 OP + Gravity_support_structure_mass_(kg)_____________________________________ (clgsmass)____________________ 1.2616E+06 OP Torus_leg_support_mass_(kg)_____________________________________________ (gsm1)________________________ 8.5820E+04 OP - Ring_beam_mass_(kg)_____________________________________________________ (gsm2)________________________ 3.9215E+05 OP - Ring_legs_mass_(kg)_____________________________________________________ (gsm3)________________________ 6.5343E+05 OP + Ring_beam_mass_(kg)_____________________________________________________ (gsm2)________________________ 3.9796E+05 OP + Ring_legs_mass_(kg)_____________________________________________________ (gsm3)________________________ 6.7293E+05 OP # PF Coil Inductances # # Nuclear Heating Magnets Before Renormalisation # Shield_line_density_(tonne/m2)__________________________________________ (x_shield)____________________ 4.0560E+00 Blanket_line_density_(tonne/m2)_________________________________________ (x_blanket)___________________ 2.2911E+00 - Unit_nuclear_heating_in_TF_coil_(W/GW)__________________________________ (tfc_nuc_heating)_____________ 1.3678E+04 - Total_nuclear_heating_in_TF_coil_(MW)___________________________________ (ptfnuc.)_____________________ 2.1784E-02 - powfmw__________________________________________________________________ (powfmw.)_____________________ 1.5926E+03 - total_mass_of_the_TF_coils_(kg)_________________________________________ (whttf)_______________________ 1.0511E+07 + Unit_nuclear_heating_in_TF_coil_(W/GW)__________________________________ (tfc_nuc_heating)_____________ 1.4472E+04 + Total_nuclear_heating_in_TF_coil_(MW)___________________________________ (ptfnuc.)_____________________ 2.3112E-02 + powfmw__________________________________________________________________ (powfmw.)_____________________ 1.5971E+03 + total_mass_of_the_TF_coils_(kg)_________________________________________ (whttf)_______________________ 1.1121E+07 # Pumping for primary coolant (helium) # Pressure_drop_in_FW_and_blanket_coolant_incl._hx_and_pipes_(Pa)_________ (dp_he)_______________________ 5.5000E+05 Fraction_of_FW_and_blanket_thermal_power_required_for_pumping___________ (fpump)_______________________ 8.9463E-02 OP - Total_power_absorbed_by_FW_&_blanket_(MW)_______________________________ (p_plasma)____________________ 1.6428E+03 OP + Total_power_absorbed_by_FW_&_blanket_(MW)_______________________________ (p_plasma)____________________ 1.6460E+03 OP Inlet_temperature_of_FW_&_blanket_coolant_pump_(K)______________________ (t_in_compressor)_____________ 5.5703E+02 OP Coolant_pump_outlet/Inlet_temperature_of_FW_&_blanket_(K)_______________ (t_in_bb)_____________________ 5.7313E+02 Outlet_temperature_of_FW_&_blanket_(K)__________________________________ (t_out_bb)____________________ 7.7313E+02 - Mechanical_pumping_power_for_FW_and_blanket_cooling_loop_including_heat_ (htpmw_fw_blkt)_______________ 1.6141E+02 OP - Mechanical_pumping_power_for_divertor_(MW)______________________________ (htpmw_div)___________________ 1.6374E+00 OP - Mechanical_pumping_power_for_shield_and_vacuum_vessel_(MW)______________ (htpmw_shld)__________________ 6.9487E-03 OP + Mechanical_pumping_power_for_FW_and_blanket_cooling_loop_including_heat_ (htpmw_fw_blkt)_______________ 1.6173E+02 OP + Mechanical_pumping_power_for_divertor_(MW)______________________________ (htpmw_div)___________________ 1.6485E+00 OP + Mechanical_pumping_power_for_shield_and_vacuum_vessel_(MW)______________ (htpmw_shld)__________________ 6.9679E-03 OP # First wall and blanket : CCFE HCPB model # Titanium_beryllide_fraction_____________________________________________ (fbltibe12)___________________ 3.7500E-01 Lithium_orthosilicate_fraction__________________________________________ (fblli2sio4)__________________ 3.7500E-01 @@ -974,11 +969,11 @@ Total_mass_of_armour,_first_wall_and_blanket_(kg)_______________________ (armour_fw_bl_mass)___________ 3.2350E+06 OP Shield_Mass_(kg)________________________________________________________ (whtshld)_____________________ 2.4498E+06 OP Vacuum_vessel_mass_(kg)_________________________________________________ (vvmass)______________________ 7.9388E+06 OP - Total_nuclear_heating_in_TF+PF_coils_(CS_is_negligible)_(MW)____________ (ptfnuc)______________________ 2.2718E-02 OP - Total_nuclear_heating_in_FW_(MW)________________________________________ (pnucfw)______________________ 1.6038E+02 OP - Total_nuclear_heating_in_the_blanket_(including_emult)_(MW)_____________ (pnucblkt)____________________ 1.2682E+03 OP - Total_nuclear_heating_in_the_shield_(MW)________________________________ (pnucshld)____________________ 1.3897E+00 OP - Total_nuclear_heating_in_the_divertor_(MW)______________________________ (pnucdiv)_____________________ 1.4652E+02 OP + Total_nuclear_heating_in_TF+PF_coils_(CS_is_negligible)_(MW)____________ (ptfnuc)______________________ 2.4103E-02 OP + Total_nuclear_heating_in_FW_(MW)________________________________________ (pnucfw)______________________ 1.6083E+02 OP + Total_nuclear_heating_in_the_blanket_(including_emult)_(MW)_____________ (pnucblkt)____________________ 1.2717E+03 OP + Total_nuclear_heating_in_the_shield_(MW)________________________________ (pnucshld)____________________ 1.3936E+00 OP + Total_nuclear_heating_in_the_divertor_(MW)______________________________ (pnucdiv)_____________________ 1.4693E+02 OP Blanket_exponential_factor______________________________________________ (exp_blanket)_________________ 9.9936E-01 OP Shield:_first_exponential_______________________________________________ (exp_shield1)_________________ 1.9523E-03 OP Shield:_second_exponential______________________________________________ (exp_shield2)_________________ 2.5427E-01 OP @@ -993,119 +988,119 @@ No_of_outboard_blanket_modules_toroidally_______________________________ (nblktmodto)__________________ 48 Isentropic_efficiency_of_first_wall_/_blanket_coolant_pumps_____________ (etaiso)______________________ 9.0000E-01 First_wall_area_(m2)____________________________________________________ (fwarea)______________________ 1.6044E+03 OP - Cryostat_internal_radius_(m)____________________________________________ (rdewex)______________________ 1.7769E+01 OP - Cryostat_internal_half-height_(m)_______________________________________ (zdewex)______________________ 1.5292E+01 OP - Vertical_clearance_from_TF_coil_to_cryostat_(m)_________________________ (clh1)________________________ 5.6091E+00 OP + Cryostat_internal_radius_(m)____________________________________________ (rdewex)______________________ 1.7787E+01 OP + Cryostat_internal_half-height_(m)_______________________________________ (zdewex)______________________ 1.5322E+01 OP + Vertical_clearance_from_TF_coil_to_cryostat_(m)_________________________ (clh1)________________________ 5.6074E+00 OP Divertor_area_(m2)______________________________________________________ (divsur)______________________ 1.4879E+02 OP Divertor_mass_(kg)______________________________________________________ (divmas)______________________ 3.6453E+04 OP # Superconducting TF Coil Power Conversion # - TF_coil_current_(kA)____________________________________________________ (itfka)_______________________ 8.0573E+01 OP + TF_coil_current_(kA)____________________________________________________ (itfka)_______________________ 8.4815E+01 OP Number_of_TF_coils______________________________________________________ (ntfc)________________________ 1.6000E+01 - Voltage_across_a_TF_coil_during_quench_(kV)_____________________________ (vtfskv)______________________ 8.7306E+00 OP + Voltage_across_a_TF_coil_during_quench_(kV)_____________________________ (vtfskv)______________________ 9.9703E+00 OP TF_coil_charge_time_(hours)_____________________________________________ (tchghr)______________________ 4.0000E+00 - Total_inductance_of_TF_coils_(H)________________________________________ (ltfth)_______________________ 3.3992E+01 OP + Total_inductance_of_TF_coils_(H)________________________________________ (ltfth)_______________________ 3.2978E+01 OP Total_resistance_of_TF_coils_(ohm)______________________________________ (rcoils)______________________ 0.0000E+00 OP - TF_coil_charging_voltage_(V)____________________________________________ (tfcv)________________________ 2.9477E+02 + TF_coil_charging_voltage_(V)____________________________________________ (tfcv)________________________ 2.9881E+02 Number_of_DC_circuit_breakers___________________________________________ (ntfbkr)______________________ 1.6000E+01 Number_of_dump_resistors________________________________________________ (ndumpr)______________________ 6.4000E+01 - Resistance_per_dump_resistor_(ohm)______________________________________ (r1dump)______________________ 1.0836E-01 OP - Dump_resistor_peak_power_(MW)___________________________________________ (r1ppmw)______________________ 1.7586E+02 OP - Energy_supplied_per_dump_resistor_(MJ)__________________________________ (r1emj)_______________________ 1.7240E+03 OP - TF_coil_L/R_time_constant_(s)___________________________________________ (ttfsec)______________________ 1.9607E+01 OP - Power_supply_voltage_(V)________________________________________________ (tfpsv)_______________________ 3.0951E+02 OP - Power_supply_current_(kA)_______________________________________________ (tfpska)______________________ 8.4601E+01 OP - DC_power_supply_rating_(kW)_____________________________________________ (tfckw)_______________________ 2.6185E+04 OP - AC_power_for_charging_(kW)______________________________________________ (tfackw)______________________ 2.9094E+04 OP - TF_coil_resistive_power_(MW)____________________________________________ (rpower)______________________ 8.4257E+00 OP - TF_coil_inductive_power_(MVA)___________________________________________ (xpower)______________________ 1.5325E+01 OP + Resistance_per_dump_resistor_(ohm)______________________________________ (r1dump)______________________ 1.1755E-01 OP + Dump_resistor_peak_power_(MW)___________________________________________ (r1ppmw)______________________ 2.1141E+02 OP + Energy_supplied_per_dump_resistor_(MJ)__________________________________ (r1emj)_______________________ 1.8534E+03 OP + TF_coil_L/R_time_constant_(s)___________________________________________ (ttfsec)______________________ 1.7534E+01 OP + Power_supply_voltage_(V)________________________________________________ (tfpsv)_______________________ 3.1375E+02 OP + Power_supply_current_(kA)_______________________________________________ (tfpska)______________________ 8.9056E+01 OP + DC_power_supply_rating_(kW)_____________________________________________ (tfckw)_______________________ 2.7942E+04 OP + AC_power_for_charging_(kW)______________________________________________ (tfackw)______________________ 3.1046E+04 OP + TF_coil_resistive_power_(MW)____________________________________________ (rpower)______________________ 8.8694E+00 OP + TF_coil_inductive_power_(MVA)___________________________________________ (xpower)______________________ 1.6475E+01 OP Aluminium_bus_current_density_(kA/cm2)__________________________________ (djmka)_______________________ 1.2500E-01 - Aluminium_bus_cross-sectional_area_(cm2)________________________________ (albusa)______________________ 6.4458E+02 OP + Aluminium_bus_cross-sectional_area_(cm2)________________________________ (albusa)______________________ 6.7852E+02 OP Total_length_of_TF_coil_bussing_(m)_____________________________________ (tfbusl)______________________ 3.1931E+03 OP - Aluminium_bus_weight_(tonnes)___________________________________________ (albuswt)_____________________ 5.5571E+02 OP - Total_TF_coil_bus_resistance_(ohm)______________________________________ (rtfbus)______________________ 1.2979E-03 OP + Aluminium_bus_weight_(tonnes)___________________________________________ (albuswt)_____________________ 5.8497E+02 OP + Total_TF_coil_bus_resistance_(ohm)______________________________________ (rtfbus)______________________ 1.2330E-03 OP TF_coil_bus_voltage_drop_(V)____________________________________________ (vtfbus)______________________ 1.0457E+02 OP - Dump_resistor_floor_area_(m2)___________________________________________ (drarea)______________________ 4.6142E+03 OP - TF_coil_power_conversion_floor_space_(m2)_______________________________ (tfcfsp)______________________ 1.6483E+03 OP - TF_coil_power_conv._building_volume_(m3)________________________________ (tfcbv)_______________________ 9.8896E+03 OP - TF_coil_AC_inductive_power_demand_(MW)__________________________________ (xpwrmw)______________________ 1.7027E+01 OP - Total_steady_state_AC_power_demand_(MW)_________________________________ (tfacpd)______________________ 9.3619E+00 OP + Dump_resistor_floor_area_(m2)___________________________________________ (drarea)______________________ 4.8422E+03 OP + TF_coil_power_conversion_floor_space_(m2)_______________________________ (tfcfsp)______________________ 1.8165E+03 OP + TF_coil_power_conv._building_volume_(m3)________________________________ (tfcbv)_______________________ 1.0899E+04 OP + TF_coil_AC_inductive_power_demand_(MW)__________________________________ (xpwrmw)______________________ 1.8305E+01 OP + Total_steady_state_AC_power_demand_(MW)_________________________________ (tfacpd)______________________ 9.8548E+00 OP # PF Coils and Central Solenoid: Power and Energy # Number_of_PF_coil_circuits______________________________________________ (pfckts)______________________ 1.2000E+01 - Sum_of_PF_power_supply_ratings_(MVA)____________________________________ (spsmva)______________________ 2.8939E+02 OP + Sum_of_PF_power_supply_ratings_(MVA)____________________________________ (spsmva)______________________ 2.9630E+02 OP Total_PF_coil_circuit_bus_length_(m)____________________________________ (spfbusl)_____________________ 2.4480E+03 OP - Total_PF_coil_bus_resistive_power_(kW)__________________________________ (pfbuspwr)____________________ 9.6444E+02 OP - Total_PF_coil_resistive_power_(kW)______________________________________ (srcktpm)_____________________ 9.6444E+02 OP + Total_PF_coil_bus_resistive_power_(kW)__________________________________ (pfbuspwr)____________________ 9.6460E+02 OP + Total_PF_coil_resistive_power_(kW)______________________________________ (srcktpm)_____________________ 9.6460E+02 OP Maximum_PF_coil_voltage_(kV)____________________________________________ (vpfskv)______________________ 2.0000E+01 Efficiency_of_transfer_of_PF_stored_energy_into_or_out_of_storage_______ (etapsu)______________________ 9.0000E-01 - Maximum_stored_energy_in_poloidal_field_(MJ)____________________________ (ensxpfm)_____________________ 2.8001E+04 OP - Peak_absolute_rate_of_change_of_stored_energy_in_poloidal_field_(MW)____ (peakpoloidalpower)___________ 1.6768E+02 OP + Maximum_stored_energy_in_poloidal_field_(MJ)____________________________ (ensxpfm)_____________________ 2.8232E+04 OP + Peak_absolute_rate_of_change_of_stored_energy_in_poloidal_field_(MW)____ (peakpoloidalpower)___________ 1.6991E+02 OP # Vacuum System # First_wall_outgassing_rate_(Pa_m/s)_____________________________________ (rat)_________________________ 1.3000E-08 Total_outgassing_load_(Pa_m3/s)_________________________________________ (ogas)________________________ 1.6691E-04 OP Base_pressure_required_(Pa)_____________________________________________ (pbase)_______________________ 5.0000E-04 Required_N2_pump_speed_(m3/s)___________________________________________ (s(1))________________________ 3.3381E-01 OP - N2_pump_speed_provided_(m3/s)___________________________________________ (snet(1))_____________________ 2.3496E+01 OP + N2_pump_speed_provided_(m3/s)___________________________________________ (snet(1))_____________________ 2.2321E+01 OP Plasma_chamber_volume_(m3)______________________________________________ (volume)______________________ 2.2588E+03 OP - Chamber_pressure_after_burn_(Pa)________________________________________ (pend)________________________ 1.6766E-01 OP - Chamber_pressure_before_burn_(Pa)_______________________________________ (pstart)______________________ 1.6766E-03 + Chamber_pressure_after_burn_(Pa)________________________________________ (pend)________________________ 1.6683E-01 OP + Chamber_pressure_before_burn_(Pa)_______________________________________ (pstart)______________________ 1.6683E-03 Allowable_pumping_time_switch___________________________________________ (dwell_pump)__________________ 0 Dwell_time_between_burns_(s)____________________________________________ (tdwell.)_____________________ 1.8000E+03 CS_ramp-up_time_burns_(s)_______________________________________________ (tramp.)______________________ 5.0000E+02 Allowable_pumping_time_between_burns_(s)________________________________ (tpump)_______________________ 1.8000E+03 Required_D-T_pump_speed_(m3/s)__________________________________________ (s(2))________________________ 5.7790E+00 OP - D-T_pump_speed_provided_(m3/s)__________________________________________ (snet(2))_____________________ 5.6929E+01 OP + D-T_pump_speed_provided_(m3/s)__________________________________________ (snet(2))_____________________ 5.4016E+01 OP Divertor_chamber_gas_pressure_(Pa)______________________________________ (prdiv)_______________________ 3.6000E-01 - Helium_gas_fraction_in_divertor_chamber_________________________________ (fhe)_________________________ 1.6884E-01 OP - Required_helium_pump_speed_(m3/s)_______________________________________ (s(3))________________________ 3.8516E+01 OP - Helium_pump_speed_provided_(m3/s)_______________________________________ (snet(3))_____________________ 3.8516E+01 OP - D-T_fuelling_rate_(kg/s)________________________________________________ (frate)_______________________ 2.7815E-05 OP - Required_D-T_pump_speed_(m3/s)__________________________________________ (s(4))________________________ 3.8516E+01 OP - D-T_pump_speed_provided_(m3/s)__________________________________________ (snet(4))_____________________ 5.6929E+01 OP + Helium_gas_fraction_in_divertor_chamber_________________________________ (fhe)_________________________ 1.7483E-01 OP + Required_helium_pump_speed_(m3/s)_______________________________________ (s(3))________________________ 3.7300E+01 OP + Helium_pump_speed_provided_(m3/s)_______________________________________ (snet(3))_____________________ 3.7300E+01 OP + D-T_fuelling_rate_(kg/s)________________________________________________ (frate)_______________________ 2.6937E-05 OP + Required_D-T_pump_speed_(m3/s)__________________________________________ (s(4))________________________ 3.7300E+01 OP + D-T_pump_speed_provided_(m3/s)__________________________________________ (snet(4))_____________________ 5.4016E+01 OP Number_of_large_pump_ducts______________________________________________ (nduct)_______________________ 16 - Passage_diameter,_divertor_to_ducts_(m)_________________________________ (d(imax))_____________________ 4.7824E-01 OP - Passage_length_(m)______________________________________________________ (l1)__________________________ 1.6648E+00 OP - Diameter_of_ducts_(m)___________________________________________________ (dout)________________________ 5.7389E-01 OP + Passage_diameter,_divertor_to_ducts_(m)_________________________________ (d(imax))_____________________ 4.6838E-01 OP + Passage_length_(m)______________________________________________________ (l1)__________________________ 1.6968E+00 OP + Diameter_of_ducts_(m)___________________________________________________ (dout)________________________ 5.6206E-01 OP Duct_length,_divertor_to_elbow_(m)______________________________________ (l2)__________________________ 4.8000E+00 OP Duct_length,_elbow_to_pumps_(m)_________________________________________ (l3)__________________________ 2.0000E+00 Number_of_pumps_________________________________________________________ (pumpn)_______________________ 3.2000E+01 OP # Plant Buildings System # - Internal_volume_of_reactor_building_(m3)________________________________ (vrci)________________________ 9.8942E+05 - Dist_from_centre_of_torus_to_bldg_wall_(m)______________________________ (wrbi)________________________ 3.9426E+01 - Effective_floor_area_(m2)_______________________________________________ (efloor)______________________ 3.3497E+05 - Reactor_building_volume_(m3)____________________________________________ (rbv)_________________________ 1.1226E+06 + Internal_volume_of_reactor_building_(m3)________________________________ (vrci)________________________ 9.9419E+05 + Dist_from_centre_of_torus_to_bldg_wall_(m)______________________________ (wrbi)________________________ 3.9491E+01 + Effective_floor_area_(m2)_______________________________________________ (efloor)______________________ 3.3604E+05 + Reactor_building_volume_(m3)____________________________________________ (rbv)_________________________ 1.1278E+06 Reactor_maintenance_building_volume_(m3)________________________________ (rmbv)________________________ 3.9641E+05 Warmshop_volume_(m3)____________________________________________________ (wsv)_________________________ 1.2531E+05 Tritium_building_volume_(m3)____________________________________________ (triv)________________________ 4.0000E+04 - Electrical_building_volume_(m3)_________________________________________ (elev)________________________ 5.0890E+04 + Electrical_building_volume_(m3)_________________________________________ (elev)________________________ 5.1899E+04 Control_building_volume_(m3)____________________________________________ (conv)________________________ 6.0000E+04 - Cryogenics_building_volume_(m3)_________________________________________ (cryv)________________________ 1.4590E+04 + Cryogenics_building_volume_(m3)_________________________________________ (cryv)________________________ 1.4780E+04 Administration_building_volume_(m3)_____________________________________ (admv)________________________ 1.0000E+05 Shops_volume_(m3)_______________________________________________________ (shov)________________________ 1.0000E+05 - Total_volume_of_nuclear_buildings_(m3)__________________________________ (volnucb)_____________________ 1.5657E+06 + Total_volume_of_nuclear_buildings_(m3)__________________________________ (volnucb)_____________________ 1.5707E+06 # Electric Power Requirements # Facility_base_load_(MW)_________________________________________________ (basemw)______________________ 5.0000E+00 Divertor_coil_power_supplies_(MW)_______________________________________ (bdvmw)_______________________ 0.0000E+00 - Cryoplant_electric_power_(MW)___________________________________________ (crymw)_______________________ 3.4704E+01 OP - Primary_coolant_pumps_(MW)______________________________________________ (htpmw..)_____________________ 1.8742E+02 OP - PF_coil_power_supplies_(MW)_____________________________________________ (ppfmw)_______________________ 7.5479E+01 OP - TF_coil_power_supplies_(MW)_____________________________________________ (ptfmw)_______________________ 9.3619E+00 OP - Plasma_heating_supplies_(MW)____________________________________________ (pheatingmw)__________________ 1.5043E+02 OP + Cryoplant_electric_power_(MW)___________________________________________ (crymw)_______________________ 3.5616E+01 OP + Primary_coolant_pumps_(MW)______________________________________________ (htpmw..)_____________________ 1.8780E+02 OP + PF_coil_power_supplies_(MW)_____________________________________________ (ppfmw)_______________________ 8.0671E+01 OP + TF_coil_power_supplies_(MW)_____________________________________________ (ptfmw)_______________________ 9.8548E+00 OP + Plasma_heating_supplies_(MW)____________________________________________ (pheatingmw)__________________ 1.5080E+02 OP Tritium_processing_(MW)_________________________________________________ (trithtmw..)__________________ 1.5000E+01 Vacuum_pumps__(MW)______________________________________________________ (vachtmw..)___________________ 5.0000E-01 - Total_pulsed_power_(MW)_________________________________________________ (pacpmw)______________________ 5.2814E+02 OP - Total_base_power_required_at_all_times_(MW)_____________________________ (fcsht)_______________________ 5.5246E+01 OP + Total_pulsed_power_(MW)_________________________________________________ (pacpmw)______________________ 5.3564E+02 OP + Total_base_power_required_at_all_times_(MW)_____________________________ (fcsht)_______________________ 5.5405E+01 OP # Cryogenics # - Conduction_and_radiation_heat_loads_on_cryogenic_components_(MW)________ (qss/1.0d6)___________________ 1.4287E-02 OP + Conduction_and_radiation_heat_loads_on_cryogenic_components_(MW)________ (qss/1.0d6)___________________ 1.4613E-02 OP Nuclear_heating_of_cryogenic_components_(MW)____________________________ (qnuc/1.0d6)__________________ 1.3000E-02 OP - AC_losses_in_cryogenic_components_(MW)__________________________________ (qac/1.0d6)___________________ 3.7116E-03 OP - Resistive_losses_in_current_leads_(MW)__________________________________ (qcl/1.0d6)___________________ 1.7533E-02 OP - 45%_allowance_for_heat_loads_in_transfer_lines,_storage_tanks_etc_(MW)__ (qmisc/1.0d6)_________________ 2.1839E-02 OP - Sum_=_Total_heat_removal_at_cryogenic_temperatures_(tfcoil_variables.tmp (helpow_+_helpow_cryal/1.0d6)_ 7.0370E-02 OP + AC_losses_in_cryogenic_components_(MW)__________________________________ (qac/1.0d6)___________________ 3.7379E-03 OP + Resistive_losses_in_current_leads_(MW)__________________________________ (qcl/1.0d6)___________________ 1.8456E-02 OP + 45%_allowance_for_heat_loads_in_transfer_lines,_storage_tanks_etc_(MW)__ (qmisc/1.0d6)_________________ 2.2413E-02 OP + Sum_=_Total_heat_removal_at_cryogenic_temperatures_(tfcoil_variables.tmp (helpow_+_helpow_cryal/1.0d6)_ 7.2219E-02 OP Temperature_of_cryogenic_superconducting_components_(K)_________________ (tmpcry)______________________ 4.5000E+00 Temperature_of_cryogenic_aluminium_components_(K)_______________________ (tcoolin)_____________________ 3.1315E+02 Efficiency_(figure_of_merit)_of_cryogenic_plant_is_13%_of_ideal_Carnot_v ______________________________ 2.0277E-03 OP Efficiency_(figure_of_merit)_of_cryogenic_aluminium_plant_is_40%_of_idea ______________________________ -2.0203E+00 OP - Electric_power_for_cryogenic_plant_(MW)_________________________________ (crypmw)______________________ 3.4704E+01 OP + Electric_power_for_cryogenic_plant_(MW)_________________________________ (crypmw)______________________ 3.5616E+01 OP # Plant Power / Heat Transport Balance # Neutron_power_multiplication_in_blanket_________________________________ (emult)_______________________ 1.2690E+00 Divertor_area_fraction_of_whole_toroid_surface__________________________ (fdiv)________________________ 1.1500E-01 @@ -1114,77 +1109,77 @@ Switch_for_pumping_of_primary_coolant___________________________________ (primary_pumping)_____________ 3 Mechanical_pumping_power_for_FW_cooling_loop_including_heat_exchanger_(M (htpmw_fw)____________________ 0.0000E+00 OP Mechanical_pumping_power_for_blanket_cooling_loop_including_heat_exchang (htpmw_blkt)__________________ 0.0000E+00 OP - Mechanical_pumping_power_for_FW_and_blanket_cooling_loop_including_heat_ (htpmw_fw_blkt)_______________ 1.6141E+02 OP - Mechanical_pumping_power_for_divertor_(MW)______________________________ (htpmw_div)___________________ 1.6374E+00 OP - Mechanical_pumping_power_for_shield_and_vacuum_vessel_(MW)______________ (htpmw_shld)__________________ 6.9487E-03 OP - Electrical_pumping_power_for_FW_and_blanket_(MW)________________________ (htpmwe_fw_blkt)______________ 1.8553E+02 OP - Electrical_pumping_power_for_shield_(MW)________________________________ (htpmwe_shld)_________________ 7.9870E-03 OP - Electrical_pumping_power_for_divertor_(MW)______________________________ (htpmwe_div)__________________ 1.8821E+00 OP - Total_electrical_pumping_power_for_primary_coolant_(MW)_________________ (htpmw)_______________________ 1.8742E+02 OP + Mechanical_pumping_power_for_FW_and_blanket_cooling_loop_including_heat_ (htpmw_fw_blkt)_______________ 1.6173E+02 OP + Mechanical_pumping_power_for_divertor_(MW)______________________________ (htpmw_div)___________________ 1.6485E+00 OP + Mechanical_pumping_power_for_shield_and_vacuum_vessel_(MW)______________ (htpmw_shld)__________________ 6.9679E-03 OP + Electrical_pumping_power_for_FW_and_blanket_(MW)________________________ (htpmwe_fw_blkt)______________ 1.8590E+02 OP + Electrical_pumping_power_for_shield_(MW)________________________________ (htpmwe_shld)_________________ 8.0091E-03 OP + Electrical_pumping_power_for_divertor_(MW)______________________________ (htpmwe_div)__________________ 1.8948E+00 OP + Total_electrical_pumping_power_for_primary_coolant_(MW)_________________ (htpmw)_______________________ 1.8780E+02 OP Coolant_pump_power_/_non-pumping_thermal_power_in_shield________________ (fpumpshld)___________________ 5.0000E-03 Coolant_pump_power_/_non-pumping_thermal_power_in_divertor______________ (fpumpdiv)____________________ 5.0000E-03 Electrical_efficiency_of_heat_transport_coolant_pumps___________________ (etahtp)______________________ 8.7000E-01 Thermal_to_electric_conversion_efficiency_of_the_power_conversion_cycle_ (etath)_______________________ 4.0000E-01 - Fraction_of_total_high-grade_thermal_power_to_divertor__________________ (pdivfraction)________________ 1.5417E-01 - Total_power_leaving_reactor_(across_vacuum_vessel_boundary)_(MW)________ ______________________________ 1.9734E+03 - Heat_removal_from_cryogenic_plant_(MW)__________________________________ (crypmw)______________________ 3.4704E+01 - Heat_removal_from_facilities_(MW)_______________________________________ (fachtmw)_____________________ 5.5246E+01 - Coolant_pumping_efficiency_losses_(MW)__________________________________ (htpsecmw)____________________ 2.4365E+01 - Heat_removal_from_injection_power_(MW)__________________________________ (pinjht)______________________ 7.5213E+01 + Fraction_of_total_high-grade_thermal_power_to_divertor__________________ (pdivfraction)________________ 1.5480E-01 + Total_power_leaving_reactor_(across_vacuum_vessel_boundary)_(MW)________ ______________________________ 1.9788E+03 + Heat_removal_from_cryogenic_plant_(MW)__________________________________ (crypmw)______________________ 3.5616E+01 + Heat_removal_from_facilities_(MW)_______________________________________ (fachtmw)_____________________ 5.5405E+01 + Coolant_pumping_efficiency_losses_(MW)__________________________________ (htpsecmw)____________________ 2.4414E+01 + Heat_removal_from_injection_power_(MW)__________________________________ (pinjht)______________________ 7.5398E+01 Heat_removal_from_tritium_plant_(MW)____________________________________ (trithtmw)____________________ 1.5000E+01 Heat_removal_from_vacuum_pumps_(MW)_____________________________________ (vachtmw)_____________________ 5.0000E-01 TF_coil_resistive_power_(MW)____________________________________________ (tfcmw)_______________________ 0.0000E+00 - Total_low-grade_thermal_power_(MW)______________________________________ (psechtmw)____________________ 2.1566E+02 - Total_High-grade_thermal_power_(MW)_____________________________________ (pthermmw)____________________ 2.1348E+03 + Total_low-grade_thermal_power_(MW)______________________________________ (psechtmw)____________________ 2.1745E+02 + Total_High-grade_thermal_power_(MW)_____________________________________ (pthermmw)____________________ 2.1405E+03 Number_of_primary_heat_exchangers_______________________________________ (nphx)________________________ 3 OP - Transport_power_from_scaling_law_(MW)___________________________________ (pscalingmw)__________________ 2.8997E+02 - Radiation_power_from_inside_"coreradius"_(MW)___________________________ (pcoreradmw.)_________________ 8.9376E+01 - Total_(MW)______________________________________________________________ ______________________________ 3.7935E+02 - Alpha_power_deposited_in_plasma_(MW)____________________________________ (falpha*palpmw)_______________ 3.0224E+02 - Power_from_charged_products_of_DD_and/or_D-He3_fusion_(MW)______________ (pchargemw.)__________________ 1.2322E+00 - Ohmic_heating_(MW)______________________________________________________ (pohmmw.)_____________________ 6.6172E-01 - Injected_power_deposited_in_plasma_(MW)_________________________________ (pinjmw)______________________ 7.5213E+01 - Total_(MW)______________________________________________________________ ______________________________ 3.7935E+02 - Fusion_power_(MW)_______________________________________________________ (powfmw)______________________ 1.5926E+03 - Power_from_energy_multiplication_in_blanket_and_shield_(MW)_____________ (emultmw)_____________________ 3.0312E+02 - Injected_power_(MW)_____________________________________________________ (pinjmw.)_____________________ 7.5213E+01 - Ohmic_power_(MW)________________________________________________________ (pohmmw.)_____________________ 6.6172E-01 - Power_deposited_in_primary_coolant_by_pump_(MW)_________________________ (htpmw_mech)__________________ 1.6306E+02 - Total_(MW)______________________________________________________________ ______________________________ 2.1347E+03 - Heat_extracted_from_first_wall_and_blanket_(MW)_________________________ (pthermfw_blkt)_______________ 1.8043E+03 - Heat_extracted_from_shield__(MW)________________________________________ (pthermshld)__________________ 1.3967E+00 - Heat_extracted_from_divertor_(MW)_______________________________________ (pthermdiv)___________________ 3.2912E+02 + Transport_power_from_scaling_law_(MW)___________________________________ (pscalingmw)__________________ 2.9018E+02 + Radiation_power_from_inside_"coreradius"_(MW)___________________________ (pcoreradmw.)_________________ 9.0174E+01 + Total_(MW)______________________________________________________________ ______________________________ 3.8036E+02 + Alpha_power_deposited_in_plasma_(MW)____________________________________ (falpha*palpmw)_______________ 3.0308E+02 + Power_from_charged_products_of_DD_and/or_D-He3_fusion_(MW)______________ (pchargemw.)__________________ 1.2453E+00 + Ohmic_heating_(MW)______________________________________________________ (pohmmw.)_____________________ 6.3634E-01 + Injected_power_deposited_in_plasma_(MW)_________________________________ (pinjmw)______________________ 7.5398E+01 + Total_(MW)______________________________________________________________ ______________________________ 3.8036E+02 + Fusion_power_(MW)_______________________________________________________ (powfmw)______________________ 1.5971E+03 + Power_from_energy_multiplication_in_blanket_and_shield_(MW)_____________ (emultmw)_____________________ 3.0396E+02 + Injected_power_(MW)_____________________________________________________ (pinjmw.)_____________________ 7.5398E+01 + Ohmic_power_(MW)________________________________________________________ (pohmmw.)_____________________ 6.3634E-01 + Power_deposited_in_primary_coolant_by_pump_(MW)_________________________ (htpmw_mech)__________________ 1.6338E+02 + Total_(MW)______________________________________________________________ ______________________________ 2.1404E+03 + Heat_extracted_from_first_wall_and_blanket_(MW)_________________________ (pthermfw_blkt)_______________ 1.8078E+03 + Heat_extracted_from_shield__(MW)________________________________________ (pthermshld)__________________ 1.4006E+00 + Heat_extracted_from_divertor_(MW)_______________________________________ (pthermdiv)___________________ 3.3135E+02 Nuclear_and_photon_power_lost_to_H/CD_system_(MW)_______________________ (psechcd)_____________________ 0.0000E+00 - Nuclear_power_lost_to_TF_(MW)___________________________________________ (ptfnuc)______________________ 2.2718E-02 - Total_(MW)______________________________________________________________ ______________________________ 2.1348E+03 + Nuclear_power_lost_to_TF_(MW)___________________________________________ (ptfnuc)______________________ 2.4103E-02 + Total_(MW)______________________________________________________________ ______________________________ 2.1405E+03 Net_electric_power_output(MW)___________________________________________ (pnetelmw.)___________________ 4.0000E+02 Required_Net_electric_power_output(MW)__________________________________ (pnetelin)____________________ 4.0000E+02 - Electric_power_for_heating_and_current_drive_(MW)_______________________ (pinjwp)______________________ 1.5043E+02 - Electric_power_for_primary_coolant_pumps_(MW)___________________________ (htpmw)_______________________ 1.8742E+02 + Electric_power_for_heating_and_current_drive_(MW)_______________________ (pinjwp)______________________ 1.5080E+02 + Electric_power_for_primary_coolant_pumps_(MW)___________________________ (htpmw)_______________________ 1.8780E+02 Electric_power_for_vacuum_pumps_(MW)____________________________________ (vachtmw)_____________________ 5.0000E-01 Electric_power_for_tritium_plant_(MW)___________________________________ (trithtmw)____________________ 1.5000E+01 - Electric_power_for_cryoplant_(MW)_______________________________________ (crypmw)______________________ 3.4704E+01 - Electric_power_for_TF_coils_(MW)________________________________________ (tfacpd)______________________ 9.3619E+00 - Electric_power_for_PF_coils_(MW)________________________________________ (pfwpmw)______________________ 1.2444E+00 - All_other_internal_electric_power_requirements_(MW)_____________________ (fachtmw)_____________________ 5.5246E+01 - Total_(MW)______________________________________________________________ (tot_plant_power)_____________ 8.5391E+02 - Total_(MW)______________________________________________________________ ______________________________ 8.5391E+02 - Gross_electrical_output*_(MW)___________________________________________ (pgrossmw)____________________ 8.5391E+02 - Fusion_power_(MW)_______________________________________________________ (powfmw)______________________ 1.5926E+03 - Power_from_energy_multiplication_in_blanket_and_shield_(MW)_____________ (emultmw)_____________________ 3.0312E+02 - Total_(MW)______________________________________________________________ ______________________________ 1.8958E+03 + Electric_power_for_cryoplant_(MW)_______________________________________ (crypmw)______________________ 3.5616E+01 + Electric_power_for_TF_coils_(MW)________________________________________ (tfacpd)______________________ 9.8548E+00 + Electric_power_for_PF_coils_(MW)________________________________________ (pfwpmw)______________________ 1.2352E+00 + All_other_internal_electric_power_requirements_(MW)_____________________ (fachtmw)_____________________ 5.5405E+01 + Total_(MW)______________________________________________________________ (tot_plant_power)_____________ 8.5621E+02 + Total_(MW)______________________________________________________________ ______________________________ 8.5621E+02 + Gross_electrical_output*_(MW)___________________________________________ (pgrossmw)____________________ 8.5621E+02 + Fusion_power_(MW)_______________________________________________________ (powfmw)______________________ 1.5971E+03 + Power_from_energy_multiplication_in_blanket_and_shield_(MW)_____________ (emultmw)_____________________ 3.0396E+02 + Total_(MW)______________________________________________________________ ______________________________ 1.9010E+03 Net_electrical_output_(MW) _____________________________________________ (pnetelmw)____________________ 4.0000E+02 - Heat_rejected_by_main_power_conversion_circuit_(MW)_____________________ (rejected_main)_______________ 1.2809E+03 - Heat_rejected_by_other_cooling_circuits_(MW)____________________________ (psechtmw)____________________ 2.1566E+02 - Total_(MW)______________________________________________________________ ______________________________ 1.8965E+03 - Net_electric_power_/_total_nuclear_power_(%)____________________________ (pnetelmw/(powfmw+emultmw)____ 2.1100E+01 - Net_electric_power_/_total_fusion_power_(%)_____________________________ (pnetelmw/powfmw)_____________ 2.5116E+01 + Heat_rejected_by_main_power_conversion_circuit_(MW)_____________________ (rejected_main)_______________ 1.2843E+03 + Heat_rejected_by_other_cooling_circuits_(MW)____________________________ (psechtmw)____________________ 2.1745E+02 + Total_(MW)______________________________________________________________ ______________________________ 1.9018E+03 + Net_electric_power_/_total_nuclear_power_(%)____________________________ (pnetelmw/(powfmw+emultmw)____ 2.1041E+01 + Net_electric_power_/_total_fusion_power_(%)_____________________________ (pnetelmw/powfmw)_____________ 2.5046E+01 Gross_electric_power*_/_high_grade_heat_(%)_____________________________ (etath)_______________________ 4.0000E+01 - Recirculating_power_fraction____________________________________________ (cirpowfr)____________________ 5.3156E-01 + Recirculating_power_fraction____________________________________________ (cirpowfr)____________________ 5.3282E-01 # Water usage during plant operation (secondary cooling) # - Volume_used_in_cooling_tower_(m3/day)___________________________________ (waterusetower)_______________ 5.4455E+04 OP - Volume_used_in_recirculating_water_system_(m3/day)______________________ (wateruserecirc)______________ 1.8271E+04 OP - Volume_used_in_once-through_water_system_(m3/day)_______________________ (wateruseonethru)_____________ 1.7906E+06 OP + Volume_used_in_cooling_tower_(m3/day)___________________________________ (waterusetower)_______________ 5.4602E+04 OP + Volume_used_in_recirculating_water_system_(m3/day)______________________ (wateruserecirc)______________ 1.8320E+04 OP + Volume_used_in_once-through_water_system_(m3/day)_______________________ (wateruseonethru)_____________ 1.7954E+06 OP # Errors and Warnings # # Errors and Warnings # PROCESS_error_status_flag_______________________________________________ (error_status)________________ 2 @@ -1356,9 +1351,8 @@ taulimit = 5.0 * dump time constraint for VV stresses * *--------------------------------------* icc = 65 -ixc = 113 * ftaucq -ftaucq = 0.9 -sigvvall = 1.0e8 +ixc = 113 * fmaxvvstress +fmaxvvstress = 1.0 * CS stress limit * *-----------------* @@ -1505,7 +1499,7 @@ d_vv_top = 0.3 * Underside vacuum vessel radial thickness [m] d_vv_bot = 0.3 -* Inboard vacuum vessel thickness [m] +* Inboard radiation shield radial thickness [m] shldith = 0.3 * Gap between vacuum vessel and blanket [m]