ife_variables

Module containing global variables relating to the inertial fusion energy model

Default IFE builds and material volumes are those for the SOMBRERO device. The 2-dimensional arrays have indices (region, material), where 'region' is the region and maxmat is the 'material':

• 'region' = 1 radially outside chamber
• 'region' = 2 above chamber
• 'region' = 3 below chamber

MAXMAT = 8 module-attribute

Total number of materials in IFE device. Material numbers are as follows: - =0 void - =1 steel - =2 carbon cloth - =3 FLiBe - =4 lithium oxide Li2O - =5 concrete - =6 helium - =7 xenon - =8 lithium

bldr = None module-attribute

radial thickness of IFE blanket (m; calculated if ifetyp=4)

bldrc = None module-attribute

radial thickness of IFE curtain (m; ifetyp=4)

bldzl = None module-attribute

vertical thickness of IFE blanket below chamber (m)

bldzu = None module-attribute

vertical thickness of IFE blanket above chamber (m)

blmatf = None module-attribute

IFE blanket material fractions

blmatm = None module-attribute

IFE blanket material masses (kg)

blmatv = None module-attribute

IFE blanket material volumes (m3)

blvol = None module-attribute

IFE blanket volume (m3)

cdriv0 = None module-attribute

IFE generic/laser driver cost at edrive=0 (M$)

cdriv1 = None module-attribute

IFE low energy heavy ion beam driver cost extrapolated to edrive=0 (M$)

cdriv2 = None module-attribute

IFE high energy heavy ion beam driver cost extrapolated to edrive=0 (M$)

cdriv3 = None module-attribute

IFE driver cost ($/J wall plug) (ifedrv==3)

chdzl = None module-attribute

vertical thickness of IFE chamber below centre (m)

chdzu = None module-attribute

vertical thickness of IFE chamber above centre (m)

chmatf = None module-attribute

IFE chamber material fractions

chmatm = None module-attribute

IFE chamber material masses (kg)

chmatv = None module-attribute

IFE chamber material volumes (m3)

chrad = None module-attribute

radius of IFE chamber (m) (iteration variable 84)

chvol = None module-attribute

IFE chamber volume (m3)

dcdrv0 = None module-attribute

IFE generic/laser driver cost gradient (M$/MJ)

dcdrv1 = None module-attribute

HIB driver cost gradient at low energy (M$/MJ)

dcdrv2 = None module-attribute

HIB driver cost gradient at high energy (M$/MJ)

drveff = None module-attribute

IFE driver wall plug to target efficiency (ifedrv=0,3) (iteration variable 82)

edrive = None module-attribute

IFE driver energy (J) (iteration variable 81)

etadrv = None module-attribute

IFE driver wall plug to target efficiency

etali = None module-attribute

IFE lithium pump wall plug efficiency (ifetyp=4)

etave = None module-attribute

IFE driver efficiency vs driver energy (ifedrv=-1)

fauxbop = None module-attribute

fraction of gross electric power to balance-of-plant (IFE)

fbreed = None module-attribute

fraction of breeder external to device core

fburn = None module-attribute

IFE burn fraction (fraction of tritium fused/target)

flirad = None module-attribute

radius of FLiBe/lithium inlet (m) (ifetyp=3,4)

fwdr = None module-attribute

radial thickness of IFE first wall (m)

fwdzl = None module-attribute

vertical thickness of IFE first wall below chamber (m)

fwdzu = None module-attribute

vertical thickness of IFE first wall above chamber (m)

fwmatf = None module-attribute

IFE first wall material fractions

fwmatm = None module-attribute

IFE first wall material masses (kg)

fwmatv = None module-attribute

IFE first wall material volumes (kg)

fwvol = None module-attribute

IFE first wall volume (m3)

gain = None module-attribute

IFE target gain

gainve = None module-attribute

IFE target gain vs driver energy (ifedrv=-1)

htpmw_ife = None module-attribute

IFE heat transport system electrical pump power (MW)

ife = None module-attribute

Switch for IFE option: - =0 use tokamak, RFP or stellarator model - =1 use IFE model

ifedrv = None module-attribute

Switch for type of IFE driver: - =-1 use gainve, etave for gain and driver efficiency - =0 use tgain, drveff for gain and driver efficiency - =1 use laser driver based on SOMBRERO design - =2 use heavy ion beam driver based on OSIRIS - =3 Input pfusife, rrin and drveff

ifetyp = None module-attribute

Switch for type of IFE device build: - =0 generic (cylindrical) build - =1 OSIRIS-like build - =2 SOMBRERO-like build - =3 HYLIFE-II-like build - =4 2019 build

lipmw = None module-attribute

IFE lithium pump power (MW; ifetyp=4)

mcdriv = None module-attribute

IFE driver cost multiplier

mflibe = None module-attribute

total mass of FLiBe (kg)

pdrive = None module-attribute

IFE driver power reaching target (W) (iteration variable 85)

pfusife = None module-attribute

IFE input fusion power (MW) (ifedrv=3 only; itv 155)

pifecr = None module-attribute

IFE cryogenic power requirements (MW)

ptargf = None module-attribute

IFE target factory power at 6 Hz repetition rate (MW)

r1 = None module-attribute

IFE device radial build (m)

r2 = None module-attribute

IFE device radial build (m)

r3 = None module-attribute

IFE device radial build (m)

r4 = None module-attribute

IFE device radial build (m)

r5 = None module-attribute

IFE device radial build (m)

r6 = None module-attribute

IFE device radial build (m)

r7 = None module-attribute

IFE device radial build (m)

reprat = None module-attribute

IFE driver repetition rate (Hz)

rrin = None module-attribute

Input IFE repetition rate (Hz) (ifedrv=3 only; itv 156)

rrmax = None module-attribute

maximum IFE repetition rate (Hz)

shdr = None module-attribute

radial thickness of IFE shield (m)

shdzl = None module-attribute

vertical thickness of IFE shield below chamber (m)

shdzu = None module-attribute

vertical thickness of IFE shield above chamber (m)

shmatf = None module-attribute

IFE shield material fractions

shmatm = None module-attribute

IFE shield material masses (kg)

shmatv = None module-attribute

IFE shield material volumes (kg)

shvol = None module-attribute

IFE shield volume (m3)

sombdr = None module-attribute

radius of cylindrical blanket section below chamber (ifetyp=2)

somtdr = None module-attribute

radius of cylindrical blanket section above chamber (ifetyp=2)

taufall = None module-attribute

Lithium Fall Time (s)

tdspmw = None module-attribute

IFE target delivery system power (MW)

tfacmw = None module-attribute

IFE target factory power (MW)

tgain = None module-attribute

IFE target gain (if ifedrv = 0) (iteration variable 83)

uccarb = None module-attribute

cost of carbon cloth ($/kg)

ucconc = None module-attribute

cost of concrete ($/kg)

ucflib = None module-attribute

cost of FLiBe ($/kg)

uctarg = None module-attribute

cost of IFE target ($/target)

v1dr = None module-attribute

radial thickness of IFE void between first wall and blanket (m)

v1dzl = None module-attribute

vertical thickness of IFE void 1 below chamber (m)

v1dzu = None module-attribute

vertical thickness of IFE void 1 above chamber (m)

v1matf = None module-attribute

IFE void 1 material fractions

v1matm = None module-attribute

IFE void 1 material masses (kg)

v1matv = None module-attribute

IFE void 1 material volumes (kg)

v1vol = None module-attribute

IFE void 1 volume (m3)

v2dr = None module-attribute

radial thickness of IFE void between blanket and shield (m)

v2dzl = None module-attribute

vertical thickness of IFE void 2 below chamber (m)

v2dzu = None module-attribute

vertical thickness of IFE void 2 above chamber (m)

v2matf = None module-attribute

IFE void 2 material fractions

v2matm = None module-attribute

IFE void 2 material masses (kg)

v2matv = None module-attribute

IFE void 2 material volumes (kg)

v2vol = None module-attribute

IFE void 2 volume (m3)

v3dr = None module-attribute

radial thickness of IFE void outside shield (m)

v3dzl = None module-attribute

vertical thickness of IFE void 3 below chamber (m)

v3dzu = None module-attribute

vertical thickness of IFE void 3 above chamber (m)

v3matf = None module-attribute

IFE void 3 material fractions

v3matm = None module-attribute

IFE void 3 material masses (kg)

v3matv = None module-attribute

IFE void 3 material volumes (kg)

v3vol = None module-attribute

IFE void 3 volume (m3)

zl1 = None module-attribute

IFE vertical build below centre (m)

zl2 = None module-attribute

IFE vertical build below centre (m)

zl3 = None module-attribute

IFE vertical build below centre (m)

zl4 = None module-attribute

IFE vertical build below centre (m)

zl5 = None module-attribute

IFE vertical build below centre (m)

zl6 = None module-attribute

IFE vertical build below centre (m)

zl7 = None module-attribute

IFE vertical build below centre (m)

zu1 = None module-attribute

IFE vertical build above centre (m)

zu2 = None module-attribute

IFE vertical build above centre (m)

zu3 = None module-attribute

IFE vertical build above centre (m)

zu4 = None module-attribute

IFE vertical build above centre (m)

zu5 = None module-attribute

IFE vertical build above centre (m)

zu6 = None module-attribute

IFE vertical build above centre (m)

zu7 = None module-attribute

IFE vertical build above centre (m)

init_ife_variables()

Source code in process/data_structure/ife_variables.py

def init_ife_variables():
    global \
        bldr, \
        bldrc, \
        bldzl, \
        bldzu, \
        blmatf, \
        blmatm, \
        blmatv, \
        blvol, \
        cdriv0, \
        cdriv1, \
        cdriv2, \
        cdriv3, \
        chdzl, \
        chdzu, \
        chmatf, \
        chmatm, \
        chmatv, \
        chrad, \
        chvol, \
        dcdrv0, \
        dcdrv1, \
        dcdrv2, \
        drveff, \
        edrive, \
        etadrv, \
        etali, \
        etave, \
        fauxbop, \
        fbreed, \
        fburn, \
        flirad, \
        fwdr, \
        fwdzl, \
        fwdzu, \
        fwmatf, \
        fwmatm, \
        fwmatv, \
        fwvol, \
        gain, \
        gainve, \
        htpmw_ife, \
        ife, \
        ifedrv, \
        ifetyp, \
        lipmw, \
        mcdriv, \
        mflibe, \
        pdrive, \
        pfusife, \
        pifecr, \
        ptargf, \
        r1, \
        r2, \
        r3, \
        r4, \
        r5, \
        r6, \
        r7, \
        reprat, \
        rrin, \
        rrmax, \
        shdr, \
        shdzl, \
        shdzu, \
        shmatf, \
        shmatm, \
        shmatv, \
        shvol, \
        sombdr, \
        somtdr, \
        taufall, \
        tdspmw, \
        tfacmw, \
        tgain, \
        uccarb, \
        ucconc, \
        ucflib, \
        uctarg, \
        v1dr, \
        v1dzl, \
        v1dzu, \
        v1matf, \
        v1matm, \
        v1matv, \
        v1vol, \
        v2dr, \
        v2dzl, \
        v2dzu, \
        v2matf, \
        v2matm, \
        v2matv, \
        v2vol, \
        v3dr, \
        v3dzl, \
        v3dzu, \
        v3matf, \
        v3matm, \
        v3matv, \
        v3vol, \
        zl1, \
        zl2, \
        zl3, \
        zl4, \
        zl5, \
        zl6, \
        zl7, \
        zu1, \
        zu2, \
        zu3, \
        zu4, \
        zu5, \
        zu6, \
        zu7
    """Initialise IFE variables"""
    bldr = 1.0
    bldrc = 1.0
    bldzl = 4.0
    bldzu = 4.0
    blmatf = np.reshape(
        [
            0.05,
            0.05,
            0.05,
            0.0,
            0.0,
            0.0,
            0.45,
            0.45,
            0.45,
            0.0,
            0.0,
            0.0,
            0.20,
            0.20,
            0.20,
            0.0,
            0.0,
            0.0,
            0.30,
            0.30,
            0.30,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
        ],
        (3, MAXMAT + 1),
    )
    blmatm = np.zeros((3, MAXMAT + 1))
    blmatv = np.zeros((3, MAXMAT + 1))
    blvol = np.zeros(3)
    cdriv0 = 154.3
    cdriv1 = 163.2
    cdriv2 = 244.9
    cdriv3 = 1.463
    chdzl = 9.0
    chdzu = 9.0
    chmatf = np.array([1.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0])
    chmatm = np.zeros(MAXMAT + 1)
    chmatv = np.zeros(MAXMAT + 1)
    chrad = 6.5
    chvol = 0.0
    dcdrv0 = 111.4
    dcdrv1 = 78.0
    dcdrv2 = 59.9
    drveff = 0.28
    edrive = 5.0e6
    etadrv = 0.0
    etali = 0.4
    etave = np.array([
        0.082,
        0.079,
        0.076,
        0.073,
        0.069,
        0.066,
        0.062,
        0.059,
        0.055,
        0.051,
    ])
    fauxbop = 0.06
    fbreed = 0.51
    fburn = 0.3333
    flirad = 0.78
    fwdr = 0.01
    fwdzl = 0.01
    fwdzu = 0.01
    fwmatf = np.reshape(
        [
            0.05,
            0.05,
            0.05,
            0.0,
            0.0,
            0.0,
            0.95,
            0.95,
            0.95,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
        ],
        (3, MAXMAT + 1),
    )
    fwmatm = np.zeros((3, MAXMAT + 1))
    fwmatv = np.zeros((3, MAXMAT + 1))
    fwvol = np.zeros(3)
    gain = 0.0
    gainve = np.array([
        60.0,
        95.0,
        115.0,
        125.0,
        133.0,
        141.0,
        152.0,
        160.0,
        165.0,
        170.0,
    ])
    htpmw_ife = 0.0
    ife = 0
    ifedrv = 2
    ifetyp = 0
    lipmw = 0.0
    mcdriv = 1.0
    mflibe = 0.0
    pdrive = 23.0e6
    pfusife = 1000.0
    pifecr = 10.0
    ptargf = 2.0
    r1 = 0.0
    r2 = 0.0
    r3 = 0.0
    r4 = 0.0
    r5 = 0.0
    r6 = 0.0
    r7 = 0.0
    reprat = 0.0
    rrin = 6.0
    rrmax = 20.0
    shdr = 1.7
    shdzl = 5.0
    shdzu = 5.0
    shmatf = np.reshape(
        [
            0.05,
            0.05,
            0.05,
            0.19,
            0.19,
            0.19,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.665,
            0.665,
            0.665,
            0.095,
            0.095,
            0.095,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
        ],
        (3, MAXMAT + 1),
    )
    shmatm = np.zeros((3, MAXMAT + 1))
    shmatv = np.zeros((3, MAXMAT + 1))
    shvol = np.zeros(3)
    sombdr = 2.7
    somtdr = 2.7
    taufall = 0.0
    tdspmw = 0.01
    tfacmw = 0.0
    tgain = 85.0
    uccarb = 50.0
    ucconc = 0.1
    ucflib = 84.0
    uctarg = 0.3
    v1dr = 0.0
    v1dzl = 0.0
    v1dzu = 0.0
    v1matf = np.reshape(
        [
            1.0,
            1.0,
            1.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
        ],
        (3, MAXMAT + 1),
    )
    v1matm = np.zeros((3, MAXMAT + 1))
    v1matv = np.zeros((3, MAXMAT + 1))
    v1vol = np.zeros(3)
    v2dr = 2.0
    v2dzl = 7.0
    v2dzu = 7.0
    v2matf = np.reshape(
        [
            1.0,
            1.0,
            1.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
        ],
        (3, MAXMAT + 1),
    )
    v2matm = np.zeros((3, MAXMAT + 1))
    v2matv = np.zeros((3, MAXMAT + 1))
    v2vol = np.zeros(3)
    v3dr = 43.3
    v3dzl = 30.0
    v3dzu = 20.0
    v3matf = np.reshape(
        [
            1.0,
            1.0,
            1.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
            0.0,
        ],
        (3, MAXMAT + 1),
    )
    v3matm = np.zeros((3, MAXMAT + 1))
    v3matv = np.zeros((3, MAXMAT + 1))
    v3vol = np.zeros(3)
    zl1 = 0.0
    zl2 = 0.0
    zl3 = 0.0
    zl4 = 0.0
    zl5 = 0.0
    zl6 = 0.0
    zl7 = 0.0
    zu1 = 0.0
    zu2 = 0.0
    zu3 = 0.0
    zu4 = 0.0
    zu5 = 0.0
    zu6 = 0.0
    zu7 = 0.0