Costing Accounts Module Equations
Revision: 2/2012 by L. Carlson
Code Version: ASC v.35
Code Modules: CostingAccount.cpp, CostingAccount.h
History: Z. Dragojlovic received from L. Waganer 12/2009. Since then, L. Waganer and L. El-Guebaly have made a few modifications, mainly related to replaceable and permanent blanket costs.
Module note: This module applies detailed and intricate costing algorithms to the power plant. The CostingAccounts.cpp module creates an output file and lists every equation that it is using to calculate the costs. There are no numbers in this file, only the equations used within. Viewing this output file is the best way to stay up to date on the current algorithms in use. The costing accounts vary slightly depending on the blanket in use.
References: Complete costing account documentation by L. Waganer 10/2010 “Cost Basis Documentation (11-25-10).doc”
MOST RECENT costing accounts and algorithms in spreadsheet from L. Waganer 12/2010 “Algorithm Comparison (RevD 120910).xls”
Costing Algorithms
The following table is a current list of the costing accounts and associated algorithms used in CostingAccount.cpp. These are printed out for display in the CostingAccounts.data and CostingAlgorithms,data output file. Both SCLL and DCLL algorithms are included and indicated where they differ.
Account Name / Title / Algorithm*20.0 / Land and Land Rights / Sum 20.1 – 20.2
20.1 / Land and Privilege Acquisition / c_20[1] = 20.0e6*pow(Pe/1000.0, 0.3)
20.2 / Relocation of Buildings, Utilities, Highways and Other Services / c_20[2] = (20.01e6)*0.2
Additional inputs required for algorithms / Pe = (*Machine).net_electrical_power()
21.0 / Structures and Site Facilities / Sum 21.1 – 21.17
21.1 / Site Improvements and Facilities / c_21[1] = 27.0e6*pow(Asite/1000.0, 0.2)
21.2 / Power Core Building / c_21[2] = 111.846e6*pow(Vrb/80.0e3, 0.62)
21.3 / Turbine - Generator Building / c_21[3] = 79.03e6*pow(Pet/1246.0, 0.5)
21.4 / Heat Rejection Structures and Facilities / c_21[4] = 16.832e6*pow((Pth-Pet)/1860.0, 0.5)
21.5 / Electrical Equipment and Power Supply Building / c_21[5] = 22.96e6*pow(Pe/1000.0, 0.3)
21.6 / Plant Auxiliary Systems Building (including Switchgear) / c_21[6] = 22.00e6*pow(Pet/1246.0, 0.3)
21.7 / Hot Cell Building (incl. Maintenance and Radioactive Material Storage and Reprocessing) / c_21[7] = 0.34*c_21[2]
21.8 / Power Core Service Building (Non-Radioactive Service?) / c_21[8] = 4.316e6*pow(Pe/1000.0, 0.3)
21.9 / Service Water Building / c_21[9] = 1.515e6*pow(Pe/1000.0, 0.3)
21.10 / Fuel Handling and Storage Building / c_21[10] = 25.0e6*pow(Pfusion/1759.0, 0.3)
21.11 / Control Room Building / c_21[11] = 7.12e6
21.12 / On-Site A/C Power Supply Building / c_21[12] = 4.71e6*pow(Pe/1000.0, 0.3)
21.13 / Administration Building / c_21[13] = 2.00e6
21.14 / Site Service Building / c_21[14] = 2.00e6
21.15 / Cryogenics and Inert Gas Storage Building / c_21[15] = 2.09e6
21.16 / Security Building / c_21[16] = 0.71e6
21.17 / Ventilation Stack / c_21[17] = 4.16e6
Additional inputs required for algorithms / Asite = 1000.0 [acres]
Pet = (*Machine).gross_electrical_power()
Pth = (*Machine).thermal_power()
Pfusion = (*Machine).fusion_power()
Vrb = (*Machine).volume() = 140120
vol_reactor_bldg = (pow((ksihi + 9.0), 2.0)*6.0*etahi + 1.55e5)*cATfit
ksihi = (*Core).bounding_box(2, 'outer')
etahi = (*Core).bounding_box(3, 'outer')
cATfit = 138698.0/181655.0
22.0 / Power Core Equipment / Sum 22.1.0 – 22.11.0
22.1.0 / Fusion Energy Capture and Conversion / Sum 22.1.1 – 22.1.7
22.1.1 / First Wall and Blanket, replaceable, IB & OB / c_22_1[1] (SCLL) = ( (*FWInboard).cost() + (*FWOutboard).cost() + (*BlanketInboard).cost() + (*BlanketOutboard).cost() )*16.0
c_22_1[1] (DCLL) = ( (*FWInboard).cost() + (*FWOutboard).cost() + (*BlanketInboard).cost() + (*BlanketOutboard).cost() + (*BackWallInboard).cost() + (*BackWallOutboard).cost() )*16.0
22.1.2 / Second Blanket, life of plant / c_22_1[2] (SCLL)= (*BlanketII).cost()*16.0
c_22_1[2] (DCLL) = 0.0*16.0 (No second blanket in DCLL model)
22.1.3 / Divertor Assembly, upper and lower, replaceable / c_22_1[3] = DivCost*16.0
22.1.4 / High Temperature Shield, outside of blanket and/or divertor, replaceable / c_22_1[4] (SCLL) = ( (*RepHTShield)[0].cost() + (*RepHTShield)[1].cost() )*16.0
c_22_1[4] (DCLL) = ( (*RepHTShield)[0].cost() + (*RepHTShield)[1].cost() + (*SkeletonRing).cost() )*16.0
22.1.5 / High Temperature Shield, outside of replaceable HT shield, life of plant / c_22_1[5] (SCLL) = ( (*HTShield)[0].cost() + (*HTShield)[1].cost() )*16.0
c_22_1[5] (DCLL) = ( (*HTShieldDCLL).cost() )*16.0
22.1.6 / Penetration Shielding (low temperature, active and passive) / c_22_1[6] = pen_shielding*c_22_1[5]
22.1.7 / Shield cooling, low temperature / c_22_1[7] = 0.0
Additional inputs required for algorithms / DivCost = for(int i=0; i<6; i++) { DivCost = DivCost + (*Div)[i].cost() }
22.2.0 / Plasma Confinement / Sum 22.2.1 – 22.2.6
22.2.1 / Toroidal Field Coils (conductor, case, bucking cyl, cryostat, struts, structure) / c_22_2[1] = ((*TFCoil).cost() + (*BuckingCylinder).cost() + (*UpperCap).cost() + (*LowerCap).cost() + (*CryoDome).cost())*16.0
22.2.1.1 / Conductor and Cable / c_22_2_1[1] = (*TFCoil).sub_cost(0)*16.0
22.2.1.2 / Structural Coil Case / c_22_2_1[2] = (*TFCoil).sub_cost(1)*16.0
22.2.1.3 / Bucking Cylinder / c_22_2_1[3] = (*BuckingCylinder).cost()*16.0
22.2.1.4 / Cryostat and Thermal Shield / c_22_2_1[4] = (*CryoDome).cost()*16.0
22.2.1.5 / Thermal Isolation Struts / c_22_2_1[5] = 0.0*16.0
22.2.1.6 / Anti-torque structure (upper + lower cap) / c_22_2_1[6] = ((*UpperCap).cost() + (*LowerCap).cost())*16.0
Additional inputs required for algorithms / (*TFCoil).sub_cost(0) = 1.95875e+06
(*TFCoil).sub_cost(1) = 6.31544e+06
22.2.2 / Poloidal Field Coils (central solenoid and outer poloidal), lower spares incl / c_22_2[2] = PFCoilCost*16.0
22.2.3 / Feedback Control Coils (normal conducting) / c_22_2[3] = 0.0
22.2.4 / Cryogenicics for Plasma Confinement (incl dewar vacuum pumping) / c_22_2[4] = 0.0
22.2.5 / Power Supplies for Plasma Confinement / c_22_2[3] = 0.0
Additional inputs required for algorithms / PFCoilCost (main): for(int i=0; i<36; i++) { PFCoilCost = PFCoilCost + (*PoloFCoil)[i].cost() }
PFCoilCost (spares): for(int i=4; i<18; i++) { PFCoilCost = PFCoilCost + (*PoloFCoil)[i].cost() }
22.3.0 / Plasma Formation and Sustainment / Sum 22.3.1 – 22.3.4
22.3.1 / Heating and Current Drive with power supplies (Primarily Steady State) / c_22_3[1] = (*Machine).cd_cost()
22.3.1.1 / Ion Cyclotron Resonance Frequency (ICRF) Fast Wave Heating and CD / c_22_3_1[1] = (*Machine).cd_fast_wave_cost()
22.3.1.2 / Lower Hybrid (LH) Wave Plasma Heating and Current Drive / c_22_3_1[2] = (*Machine).cd_lower_hybrid_cost()
22.3.1.3 / Electron Cyclotron Resonance Frequency (ECRF) Plasma Heating and CD / c_22_3_1[3] = 0.0
22.3.1.4 / Neutral Particle Beam Plasma Heating, Current Drive and Rotation / c_22_3_1[4] = (*Machine).cd_neutral_beam_cost()
22.3.2 / Startup subsystems with power supplies / c_22_3[2] = 7.0e6
22.3.2.1 / Electron Cyclotron Resonance Freq. (ECRF) Wave Plasma Breakdown / c_22_3_2[1] = 0.0
22.3.2.2 / Ion Cyclotron Resonance Freq. (ICRF) Wave Current Initiation, Ramp-Up / c_22_3_2[2] = 0.0
22.3.2.3 / Lower Hybrid (LH) Wave Current Initiation and Ramp-Up / c_22_3_2[3] = 0.0
22.3.3 / Stability Control subsystem with power supplies (Primarily Transient) / c_22_3[3] = 0.0
22.3.3.1 / Electron Cyclotron Resonance Frequency (ECRF) Wave Plasma Control / c_22_3_3[1] = 0.0
22.3.3.2 / Neutral Particle Beam Plasma Control / c_22_3_3[2] = 0.0
22.3.4 / Plasma Fueling and Constituent Control / c_22_3_4[0] = 16.0e6
22.3.4.1 / Pellet Injection Fueling and Constituent Control / c_22_3_4[1] = 0.0
22.3.4.2 / Neutral (Particle) Beam Injection (NBI) Fueling and Constituent Control / c_22_3_4[2] = 0.0
22.4.0 / Vacuum, Power Core / Sum 22.4.1 – 22.4.6
22.4.1 / Vacuum Vessel / c_22_4[1] = (cVV + cMP)*16.0
22.4.2 / Helium Liquefier-Refrigerators / c_22_4[2] = 15.0e6
22.4.3 / Primary Vacuum Pumps (cryocondensation, cryosorption, cryodiffusion, and/or turbo-molecular) / c_22_4[3] = 0.0
22.4.4 / Roughing or Backing Pumps / c_22_4[4] = 0.0
22.4.5 / Vacuum Pumping Ducts (vac duct shielding in 22.01.06) / c_22_4[5] = cDuct*16.0
22.4.6 / Plumbing, Cryogenic / c_22_4[6] = 0.0
Additional inputs required for algorithms / cVV (SCLL) = (*VacuumVessel)[0].cost() + (*VacuumVessel)[1].cost()
cVV (DCLL) = (*VacuumVesselDCLL)[0].cost() + (*VacuumVesselDCLL)[1].cost() + (*VacuumVesselDCLL)[2].cost() + (*VacuumVesselDCLL)[3].cost()
cMP = (*MPort).cost()
cDuct = (*Duct)[0].cost() + (*Duct)[1].cost() + 2.0* s(*Duct)[3].cost()
22.5.0 / Primary Structure and Support, Power Core / c_22_5[0] = support_structure_volume*0.20e6
22.5.1 / Carry-Through Structure / c_22_5[1] = 0.0
22.5.2 / Structural Pedestal / c_22_5[2] = 0.0
22.5.3 / Equipment Support Structure / c_22_5[3] = 0.0
Additional inputs required for algorithms / support_structure_volume = 0.20*fpc_vol
fpc_vol (SCLL) = ((*FWInboard).volume() + (*FWOutboard).volume() + (*BlanketInboard).volume() + (*BlanketOutboard).volume() + (*BlanketII).volume() + (*RepHTShield)[0].volume() + (*RepHTShield)[1].volume() + 1.1*((*HTShield)[0].volume() + (*HTShield)[1].volume()) + (*VacuumVessel)[0].volume() + (*VacuumVessel)[1].volume() + DivVol + (*TFCoil).volume() + PoloFCVolume + PoloFCVolumeSpares) *16.0
fpc_vol (DCLL) = ((*FWInboard).volume() + (*FWOutboard).volume() + (*BlanketInboard).volume() + (*BlanketOutboard).volume() + (*BlanketII).volume() + (*RepHTShield)[0].volume() + (*RepHTShield)[1].volume() + 1.1*((*HTShield)[0].volume() + (*HTShield)[1].volume()) + (*VacuumVessel)[0].volume() + (*VacuumVessel)[1].volume() + DivVol + (*TFCoil).volume() + PoloFCVolume + PoloFCVolumeSpares) *16.0
DivVol: for(int i=0; i<6; i++) { DivVol = DivVol + (*Div)[i].volume(); }
PoloFCVolume: for(int i=0; i<36; i++) { PoloFCVolume = PoloFCVolume + (*PoloFCoil)[i].volume() }
PoloFCVolumeSpares: for(int i=4; i<18; i++) { PoloFCVolumeSpares = PoloFCVolumeSpares + (*PoloFCoil)[i].volume() }
22.6.0 / Main Heat Transfer and Transport (MHTT) / Sum 22.6.1 – 22.6.4
22.6.1 / Primary (PbLi) Heat Transfer Loop: servicing FW, blanket & HT shield / c_22_6[1] = 125.0e6*pow(Pth*f_PbLi/2000, 0.55) + 0.01e6*Pth*f_PbLi
22.6.2 / Primary (Helium) Heat Transfer Loop: servicing FW, blanket, HT shield, divertor & Nb IHX) / c_22_6[2] = 110.0e6*pow(Pth*f_He/2000, 0.55)
22.6.3 / Limiter or Divertor Primary Heat Transfer Loop (may tie into same IHX, SG) / c_22_6[3] = 0.0
22.6.4 / Intermediate Heat Transfer Loop (Sodium or Helium) / c_22_6[4] = 0.0
Additional inputs required for algorithms / f_PbLi (SCLL) = 0.8
f_He (SCLL) = 0.2
f_PbLi (DCLL) = 0.5
f_He (DCLL) = 0.5
22.7.0 / Radioactive Materials Treatment and Management (off-line) / c_22_7[0] = 15.0e6*pow(Pfusion/1758.0, 0.8)
22.7.1 / Liquid Materials Processing Equipment (including final detritiation) / c_22_7[1] = 0.0
22.7.2 / Gaseous Materials Processing Equipment (including final detritiation) / c_22_7[2] = 0.0
22.7.3 / Solid Materials Processing Equipment (including detritiation) / c_22_7[3] = 0.0
22.8.0 / Fuel Handling and Storage (on-line) / c_22_8[0] = 70.0e6*pow(Pfusion/1758.0, 0.8)
22.8.1 / Chamber Exhaust Gas Handling and Processing Equipment / c_22_8[1] = 0.0
22.8.2 / Purge and Cover Gas Handling Processing Equipment / c_22_8[2] = 0.0
22.8.3 / Primary Coolant Stream Handling and Processing Equip. (incl. on-line adj. of Li enrichment) / c_22_8[3] = 0.0
22.8.4 / Other Liquid and Gaseous Coolant Stream Handling and Processing Equipment / c_22_8[4] = 0.0
22.8.5 / Purification and Isotope Separation Equipment / c_22_8[5] = 0.0
22.8.6 / Tritium, Deuterium, and DT Storage Equipment / c_22_8[6] = 0.0
22.8.7 / Atmospheric Tritium Recovery Equipment (Power Core, Hot Cell, Fuel Handling Bldgs.) / c_22_8[7] = 0.0
22.9.0 / Maintenance Equipment / c_22_9[0] = 100.0e6*pow(Pe/1000.0, 0.8)
22.9.1 / Power Core Maintenance Equipment / c_22_9[1] = 0.0
22.9.2 / Hot Cell Maintenance Equipment / c_22_9[2] = 0.0
22.9.3 / Fuel Handling Maintenance Equipment / c_22_9[3] = 0.0
22.9.4 / Other Plant Equipment Maintenance Equipment / c_22_9[4] = 0.0
22.10.0 / Instrumentation and Control / c_22_10[0] = 60.0e6
22.10.1 / Power Core Instrumentation & Control Equipment / c_22_10[1] = 0.0
22.10.2 / Radiation Monitoring Equipment / c_22_10[2] = 0.0
22.10.3 / Isolated Indicating & Recording Equipment / c_22_10[3] = 0.0
22.10.4 / Data Acquisition and Recording Equipment / c_22_10[4] = 0.0
22.10.5 / Communications Equipment / c_22_10[5] = 0.0
22.11.0 / Other Power Core Equipment / c_22_11[0] = 8.0e6*pow(Pe/1000.0, 0.8)
22.11.1 / Special Heating Equipment / c_22_11[1] = 0.0
22.11.2 / Special Cooling Equipment (low temp shield, vacuum vessel, other structures) / c_22_11[2] = 0.0
22.11.3 / Coolant Receiving, Storage and Makeup Equipment / c_22_11[3] = 0.0
22.11.4 / Gas Systems Equipment / c_22_11[4] = 0.0
22.11.5 / Inert Atmosphere Equipment (Power Core, Hot Cell, and Fuel Handling Buildings) / c_22_11[5] = 0.0
22.98.0 / Spare Parts Allowance (include in individual costs) / c_22_11[98] = 0.0
22.99.0 / Contingency Allowance (include in Acct 96) / c_22_11[99] = 0.0
23.0 / Turbine - Generator Equipment / c_23[0] = 360.0e6*pow(P_thermal/2000,0.8)*(eta_brayton/0.6)
23.1 / Turbine - Generators / c_23[1] = 0.0
23.2 / Main Steam or other Main Heat Transfer Fluid System / c_23[2] = 0.0
23.3 / Condensing or Heat Sink Heat Exchanger Systems / c_23[3] = 0.0
23.4 / Feedwater Heating or Heat Recovery System / c_23[4] = 0.0
23.5 / Other Turbine Plant Equipment / c_23[5] = 0.0
23.6 / Turbine Plant Instrumentation & Control / c_23[6] = 0.0
23.98 / Spare Parts Allowance (include in individual costs) / c_23[98] = 0.0
23.99 / Contingency Allowance (include in Acct 96) / c_23[99] = 0.0
Additional inputs required for algorithms / eta_brayton = (*Machine).gross_electrical_power()/Pth
24.0 / Electric Plant Equipment / c_24[0] = 187.0e6*pow(Pet/1200.0, 0.5)
24.1 / Switchgear / c_24[1] = 0.0
24.2 / Station Service Equipment / c_24[2] = 0.0
24.3 / Switchboards / c_24[3] = 0.0
24.4 / Protective Equipment / c_24[4] = 0.0
24.5 / Electrical Structures and Wiring Containers / c_24[5] = 0.0
24.6 / Power and Control Wiring / c_24[6] = 0.0
24.7 / Electrical Lighting / c_24[7] = 0.0
24.98 / Spare Parts Allowance (include in individual costs) / c_24[98] = 0.0
24.99 / Contingency Allowance (include in Acct 96) / c_24[99] = 0.0
25.0 / Heat Rejection Equipment / c_25[0] = 89.17e6*((Pth-Pet)/2300.0)
25.1 / Water Intake Common Facilities / c_25[1] = 0.0
25.2 / Circulating Water Systems / c_25[2] = 0.0
25.3 / Cooling Towers / c_25[3] = 0.0
25.4 / Other Heat Rejection Systems / c_25[4] = 0.0
25.98 / Spare Parts Allowance (incl. in individual costs) / c_25[98] = 0.0
25.99 / Contingency Allowance (incl. in Acct 96) / c_25[99] = 0.0
26.0 / Miscellaneous Plant Equipment / c_26[0] = 85.0e6*pow(Pet/1200.0, 0.6)
26.1 / Transportation and Lifting Equipment / c_26[1] = 0.0
26.2 / Air and Water Service Systems / c_26[2] = 0.0
26.3 / Communications Equipment / c_26[3] = 0.0
26.4 / Furnishings and Fixtures / c_26[4] = 0.0
26.98 / Spare Parts Allowance (incl. in individual costs) / c_26[98] = 0.0
26.99 / Contingency Allowance (incl. in Acct 96) / c_26[99] = 0.0
27.0 / Special Materials
27.1 / Primary coolant (PbLi) w/ 2.5x BOP factor / c_27[1] = (PbLi_cost/1000.0)*m_LM_total*1.0e6
27.2 / Divertor coolant (Helium) incl. storage & pressurization subsystem / c_27[2] = (he_vol_in_div + he_vol_in_mhtt + he_vol_in_storage)*he_cost*rho_ave
27.3 / Intermediate Loop Coolant / c_27[3] = 0.0
27.4 / Turbine cycle Working Fluids / c_27[4] = 0.0
27.5 / Other Special Materials / c_27[5] = 0.6e6*pow(Pe/1000.0, 0.8)
27.6 / Argon Power Core Cover Gas / c_27[6] = 0.4e6*pow(Pe/1000.0, 0.8)
Additional inputs required for algorithms / rho_ave = 4.215 = ave He density between STP and 700C, 15 MPa
he_cost = (*Div)[0].cost_unit_mass('He')
m_LM_total = BOP*(m_LM + m_LM_div + m_LM_rhts + m_LM_hts)*16.0/1000.0
m_LM (SCLL) = (*FWInboard).liquid_metal_mass() + (*FWOutboard).liquid_metal_mass() + (*BlanketInboard).liquid_metal_mass() + (*BlanketOutboard).liquid_metal_mass() + (*BlanketII).liquid_metal_mass()
m_LM (DCLL) = (*BlanketInboard).liquid_metal_mass() + (*BlanketOutboard).liquid_metal_mass()
m_LM_div: for(int i=0; i<6; i++) { m_LM_div = m_LM_div + (*Div)[i].liquid_metal_mass() }
m_LM_rhts: for(int i=0; i<2; i++) { m_LM_rhts = m_LM_rhts + (*RepHTShield)[i].liquid_metal_mass() }
m_LM_hts: for(int i=0; i<2; i++) { m_LM_hts = m_LM_hts + (*HTShield)[i].liquid_metal_mass() }
BOP factor = 2.5
he_vol_in_div = div_vol*16.0*BOP*he_frac_div
he_frac_div = 0.526
div_vol: for(int i=0; i<6; i++) { div_vol = div_vol + (*Div)[i].volume(); }
he_vol_in_storage = he_vol_in_div*0.10
mhtt_vol (SCLL) = ((*FWInboard).volume() + (*FWOutboard).volume() + (*BlanketInboard).volume() + (*BlanketOutboard).volume() + (*BlanketII).volume() + (*RepHTShield)[0].volume() + (*RepHTShield)[1].volume() + pen_shielding*((*HTShield)[0].volume() + (*HTShield)[1].volume()) )*16.0*BOP
mhtt_vol (DCLL) = ((*FWInboard).volume() + (*FWOutboard).volume() + (*BlanketInboard).volume() + (*BlanketOutboard).volume() + (*BackWallInboard).volume() + (*BackWallOutboard).volume() + (*SkeletonRing).volume() + (*RepHTShield)[1].volume() + pen_shielding*((*HTShieldDCLL).volume()) )*16.0*BOP
90.0 / Total Direct Cost / c_90 = c_20[0] + c_21[0] + c_22[0] + c_23[0] + c_24[0] + c_25[0] + c_26[0] + c_27[0]
91.0 / Construction Facilities, Equipment and Services / c_91[0] = 0.1130*c_90
91.1 / Temporary Facilities / c_91[1] = 0.0
91.2 / Construction Equipment / c_91[2] = 0.0
91.3 / Construction Services / c_91[3] = 0.0
92.0 / Home Office Engineering and Services / c_92[0] = 0.052*c_90
92.1 / Systems Engineering / c_92[1] = 0.0
92.2 / Management Services / c_92[2] = 0.0
92.3 / Quality Assurance / c_92[3] = 0.0
92.4 / Safety and Environmental Engineering / c_92[4] = 0.0
93.0 / Field Office Engineering and Services / c_93[0] = 0.052*c_90
93.1 / Construction Management / c_93[1] = 0.0
93.2 / Inspection / c_93[2] = 0.0
93.3 / Pre-Operational Training / c_93[3] = 0.0
94.0 / Owner's Cost / c_94[0] = 0.15*(c_90 + c_91[0] + c_92[0] + c_93[0])
94.1 / Project Administration / c_94[1] = 0.0
94.2 / Staff Training and Plant Startup / c_94[2] = 0.0
94.3 / Inventories and Spares Administration / c_94[3] = 0.0
95.0 / Process (design) Contingency / c_95 = 0.0*c_90
96.0 / Project Contingency / c_96 = 0.1465*(c_90 + c_91[0] + c_92[0] + c_93[0] + c_94[0] + c_95)
97.0 / Interest During Construction (IDC) / c_97 = fidc*(c_90 + c_91[0] + c_92[0] + c_93[0] + c_94[0] + c_95 + c_96)
fidc = 0.1303
98.0 / Escalation During Construction (EDC) / c_98 = fedc*c_90
fedc = 0.0
99.0 / Total Project Capital Cost / c_99 = c_90 + c_91[0] + c_92[0] + c_93[0] + c_94[0] + c_95 + c_96 + c_97 + c_98
Cost of Electricity (COE) [mills/kWh]: / cost_Of_Electricity = (c_AC + (c_OM + c_SCR + c_F) * pow(1.0 + y,Y)) / (8760.0*Pe*pf) + c_DD
** See input file costing.data for input parameters: Xrdr, Lecon, neutron_fluence_life, y, Y, pf, c_F, c_DD.
Additional inputs required for algorithms / fcr = Xrdr/( 1.0 - pow(1.0+Xrdr, -Lecon))
c_AC = fcr*c_99
Pe = (*Machine).net_electrical_power()
c_OM = 80.0e6*pow(Pe/1200.0, 0.5)
FWB_life = neutron_fluence_life/nwlmax
nwlmax = (*Machine).neutron_wall_load('maximum')
c_SCR = (c_22_1[1] + c_22_1[3] + c_22_1[4] + c_22_3_1[0] + c_22_3_4[0] + c_22_4[5]) / (FWB_life/pf)
c_AC = fcr*c_99
Subcosts of COE normalized to mills/kWh:
... capital return / coe_subcost[0] = c_AC / (8760.0*Pe*pf)
... operation and maintenance / coe_subcost[1] = c_OM*pow(1.0+y, Y) /(8760.0*Pe*pf)
... component replacement / coe_subcost[2] = c_SCR*pow(1.0+y, Y) /(8760.0*Pe*pf)
... fuel cost / coe_subcost[3] = c_F*pow(1.0+y, Y) /(8760.0*Pe*pf)
... decontamination and decommissioning / coe_subcost[4] = c_DD
* Every account is multiplied by an inflation factor depending on the cost base year in use.