Flow of Calculations in the Aries Systems Code

Revision: 10/2011 by L. Carlson

Code Version:ASC v.34

Code Modules:sdfgsf

History:sdfgsfd

Module note:sdfgsdf

This is the current understanding of how the systems code works.

See also diagram

An overview of the systems code details is partially explained in “Aries systems code user’s guide” by Zoran Dragojlovic, May 7, 2009. Also included are sample files and the definitions and symbols of the parameters.

Additionally, L. Carlson has a running Excel sheet “ASC Legend #.xls” that contains a host of information relevant to the Aries design study that began ~2009. Available upon request .

(by Z.D.)

  1. Physics module of the systems code scans the parametric space for a large number of viable physical operating points. The outputs from this module are physics parameters and radial inboard builds that define the Tokamak reactor for a given operating point.
  2. Power core elements are generated from plasma contour to the cryogenic dome. In this process, geometrical contours are defined, as well as material composition, structure and cost. The objects are generated in the following order:
  3. Plasma contour. (which module.cpp are each of these done in???)
  4. Tangents to the left and right of the plasma tip are estimated in order to determine the X-point and the locations of the divertor plates.
  5. First Wall.
  6. Blanket.
  7. Blanket II.
  8. Divertor plates.
  9. HT Shield.
  10. Vacuum Vessel.
  11. TF coil with bucking cylinder and toroidal caps.
  12. PF coil.
  13. Cryogenic dome.
  14. Power flow is calculated from the fusion power at the input to the net electric power at the output. Net electric power, gross electric power and thermal power are used as inputs to the costing account module of the systems code.
  15. Costing accounts are calculated with arbitrary cost base and level of safety assurance. (LSA should be removed ~Aug 2009 by Zoran D., right? (LCC))
  16. Cost of electricity (COE) is estimated [mill/kWh].

Systems Code Modules: (by LCC)

There are 2 main steps of the systems code:

  1. Physics – input files with user-specified parameters are called “inpar” (fixed parameters) and “inpar2” (scanned parameters). The plasmas that remain satisfy the power and particle balance and are described by 55 parameters. This output file becomes the input file for the engineering module. “sysplas_PC.F” is the name of the code that computes these parameters and “corona” is a look-up table that gives the calculation of the radiative cooling rates and other relevant parameters of impurities in the plasma.
  2. Engineering – at this step, there are 4 different modules, each with a different blanket design. COE is also calculated here. These modules are very similar with only slight differences between them.
  3. ARIES-AT – benchmark model with new (Aug 2009) costing accounts
  4. DCLL – dual-coolant lead-lithium blanket without manifold
  5. DCLL_M – with manifold (not used often)
  6. SiC – next-step version with a silicon-carbide blanket

2. Engineering:

There are 5 input files required to run the engineering module. They are called by the certain subroutines during the engineering module:

  1. magnets.data – called by geometry.cpp, contains constants and assumptions of the magnets
  2. costing.data – called by costingaccount.cpp, input parameters that are used to define costing accounts and COE
  3. powerflow.data – called by designpoint.cpp, constants provided by Rene Raffray, very similar between SiC and DCLL blankets
  4. sysout.data – called by aries.cpp, contains the output database from the physics module consisting of 55 plasma parameters
  5. builds.data – is unique only to the DCLL case and includes a detailed inboard, outboard and divertor build of the blanket. This is read by the Aries.cpp program

There are 5 subroutines of code in the engineering module under the main program “Aries.cpp”:

Aries.cpp – main program

  1. Cost.cpp – does unit costing analysis with specific material costs
  2. CostingAccount.cpp – goes through detailed costing for every account
  3. DesignPoint.cpp – might also be known as “power flow,” computes powers, heat fluxes, pumping, CD, etc
  4. Geometry.cpp – computes the geometries of the machine, as well as magnets
  5. Part.cpp – makes output parts under “/OutputParts”

And 5 corresponding header files. These call various lines in the appropriate .cpp code:

Cost.h, CostingAccount.h, DesignPoint.h, Geometry.h, Part.h

There are 2 modes in which the engineering module can be run by modifying line 50 in “Aries.cpp”):

  1. “System” – the code prints out the detailed geometry, costing accounts and other detailed data for a single data point for the purpose of debugging, etc. This output is put in the “OutputParts” subdirectory.
  2. “COE” – all the data points in the physics input files have the engineering module applied. All the unnecessary geometrical and detailed design is skipped. The output is generated in “COEMaps/SysoutFinal.data” and contains 84 parameters.

There are 5 output files in “/OutputParts” from the engineering module when “COE” is turned on for detailed analysis:

  1. "Part Name.data” – geometrical contours of the power core components used by Matlab scripts to complete the radial build of the machine
  2. “Part Name.ibl” – geometrical contours of the power core components for exporting to ProEngineer
  3. “System.data” – all system output parameters
  4. CostingAccounts.data – printout of costing accounts, cost of electricity and components of the COE.
  5. Powers.data – a list of the powers, Gross Electrical Power, Net Electrical Power, Thermal, Fusion, Divertor and Blanket Pumping Power.

There is one output files in “/COEMaps” from the engineering module when “system” is turned on for detailed analysis:

  1. "SysoutFinal.data” – contains 84 power plant design parameters per line and each line is a different design.

INPUT FILES into Engineering Module (referenced above)

(OLD now)

Magnets.data (SiC & DCLL)

Maximum stress in protective sheath around strands [Pa]*

8.000E+08

Maximum allowed stress in the coil casing [Pa] *

4.000E+08

Current density in copper [A/m^2] *

2.240E+08

Volume fraction of He coolant [-] *

2.500E-01

Volume fraction of the insulator [-] *

1.000E-01

Nominal current in TF coil [A] *

4.000E+04

Side thickness of the coil casing X 2 [m] *

2.000E-02

Clearance between V Vessel and TF coil outmost edge [m]*

1.000E-01

Yield stress for the PF coil case [Pa] *

1.000E+09

Upper limiter of the current (qmax) *

4.500E+00

Lower limiter of the current (qmin) *

3.150E+00

Maximum allowed magnetic field on PF coils [T] *

1.800E+01

Quench protection parameter J2tau *

5.000E+16

Quench time [s] *

2.000E+00

Costing.data (SiC & DCLL)

Real Discount Rate *

5.000E-02

Economic Life of Plant (Years) *

4.000E+01

Neutron Fluence Life ( MW*year/(m^2) ) *

1.500E+07

Annual Escalation Rate *

5.000E-02

Construction Period (Years) *

6.000E+00

Plant Capacity Factor *

8.500E-01

Annual Fuel Cost (mill/kWhe) *

1.000E+06

Decontamination & Decommissioning Allowance (mill/kWhe)*

3.600E-01

Powerflow.data

SiC / DCLL
// Provided by Rene Raffray
Fraction of loss of alpha particles [-] *
0.000E+00
Fraction of alpha loss to divertor [-] *
0.500E+00
Fraction of edge radiation (typical range 0.75-0.9) [-]*
0.900E+00
Fraction of edge radiation to divertor [-] *
0.900E+00
Neutron energy multiplication factor FW/blanket inboard*
1.100E+00
Neutron energy multiplication factor FW/blanket outb. *
1.100E+00
Neutron energy multiplication factor shield inboard *
1.100E+00
Neutron energy multiplication factor shield outboard *
1.100E+00
Neutron energy multiplication factor divertor plates *
1.100E+00
Neutron energy multiplication factor divertor blanket *
1.100E+00
Neutron energy multiplication factor divertor shield *
1.100E+00
Efficiency of the pump [-] *
0.900E+00
Generic current drive source efficiency [-] *
0.500E+00
Current drive source efficiency (NB) [-] *
0.680E+00
Current drive source efficiency (LH) [-] *
0.680E+00
Current drive source efficiency (ICRF) [-] *
0.840E+00
Current drive source efficiency (EC) [-] *
0.430E+00
Power to auxiliary functions [W] *
5.000E+07
Cryogenic power [W] *
2.000E+06
Efficiency of plasma heating [-] *
0.430E+00
Outboard fraction of incoming power [-] *
0.600E+00
Inboard fraction of incoming power [-] *
0.400E+00
Fraction of neutron flux absorbed by FW/blanket [-] *
0.900E+00
Fraction of neutron flux absorbed by the shield [-] *
0.100E+00
Fraction of neutron pow. to divertor in div. plates [-]*
0.500E+00
Fraction of neutron pow. to divertor in div. blnkt. [-]*
0.400E+00
Fraction of neutron pow. to divertor in div. shield [-]*
0.100E+00
Up-down symmetry factor. 1 for symmetric, 0 otherwise *
1.000E+00
Fraction of conducted power to divertor => outboard [-]*
0.800E+00
Fraction of conducted power to divertor => inboard [-] *
0.200E+00
Fraction of radiated power to divertor => outboard [-] *
0.800E+00
Fraction of radiated power to divertor => inboard [-] *
0.200E+00
Poloidal magnetic flux expansion on outboard [-] *
1.000E+01
Poloidal magnetic flux expansion on inboard [-] *
1.000E+01 / // Provided by Rene Raffray
Fraction of loss of alpha particles [-] *
0.000E+00
Fraction of alpha loss to divertor [-] *
0.500E+00
Fraction of edge radiation (typical range 0.75-0.9) [-]*
0.950E+00
Fraction of edge radiation to divertor [-] *
0.900E+00
Neutron energy multiplication factor FW/blanket inboard*
1.150E+00
Neutron energy multiplication factor FW/blanket outb. *
1.150E+00
Neutron energy multiplication factor shield inboard *
1.150E+00
Neutron energy multiplication factor shield outboard *
1.150E+00
Neutron energy multiplication factor divertor plates *
1.150E+00
Neutron energy multiplication factor divertor blanket *
1.150E+00
Neutron energy multiplication factor divertor shield *
1.150E+00
Efficiency of the pump [-] *
0.900E+00
Generic current drive source efficiency [-] *
0.500E+00
Current drive source efficiency (NB) [-] *
0.680E+00
Current drive source efficiency (LH) [-] *
0.680E+00
Current drive source efficiency (ICRF) [-] *
0.840E+00
Current drive source efficiency (EC) [-] *
0.430E+00
Power to auxiliary functions [W] *
5.000E+07
Cryogenic power [W] *
2.000E+06
Efficiency of plasma heating [-] *
0.430E+00
Outboard fraction of incoming power [-] *
0.600E+00
Inboard fraction of incoming power [-] *
0.400E+00
Fraction of neutron flux absorbed by FW/blanket [-] *
0.900E+00
Fraction of neutron flux absorbed by the shield [-] *
0.100E+00
Fraction of neutron pow. to divertor in div. plates [-]*
0.500E+00
Fraction of neutron pow. to divertor in div. blnkt. [-]*
0.400E+00
Fraction of neutron pow. to divertor in div. shield [-]*
0.100E+00
Up-down symmetry factor. 1 for symmetric, 0 otherwise *
1.000E+00
Fraction of conducted power to divertor => outboard [-]*
0.800E+00
Fraction of conducted power to divertor => inboard [-] *
0.200E+00
Fraction of radiated power to divertor => outboard [-] *
0.800E+00
Fraction of radiated power to divertor => inboard [-] *
0.200E+00
Poloidal magnetic flux expansion on outboard [-] *
1.000E+01
Poloidal magnetic flux expansion on inboard [-] *
1.000E+01

Sysout.data (sample physics output file of 55 plasma parameters)

8.000E+00 4.000E+00 4.500E+00 1.486E+01 4.000E-02 3.200E+00 2.210E+00 4.700E-01 3.000E-01 8.000E-01 7.000E-01 1.500E+01 2.200E+00 1.581E+00 9.329E+02 2.381E+00 5.028E+00 7.457E+19 2.358E+01 1.000E+30 4.902E+01 1.492E+00 3.509E+01 5.158E-01 3.300E-01 1.301E+02 1.340E+02 3.279E+01 1.280E-01 9.123E+00 1.552E+00 1.624E+00 0.000E+00 0.000E+00 5.885E-44 -4.319E-26 0.000E+00 0.000E+00 1.500E-03 0.000E+00 0.000E+00 3.808E+01 4.842E-01 1.000E+00 8.258E-02 8.078E-01 1.109E+03 3.000E+02 0.000E+00 1.295E+00 1.796E+00 1.768E+00 6.605E+00 4.025E+02 9.501E+00

Builds.data (DCLL only, SiC is included in Aries.cpp)

0 1 2 3 4 5 6 7 8 9 10 11

IB SOL First Wall Gap Blanket Gap Back Wall Gap HT Shield Skel. Ring Gap Vac. Vessel Gap

5.000E-02 3.800E-02 0.000E-02 3.620E-01 0.000E-02 5.000E-02 0.000E-02 1.500E-01 2.000E-01 2.000E-02 3.100E-01 2.000E-02

OB SOL First Wall Gap Blanket Gap Back Wall Gap HT Shield Skel. Ring Gap Vac. Vessel Gap

5.000E-02 3.800E-02 0.000E-02 7.120E-01 0.000E-02 5.000E-02 0.000E-02 0.000E-01 2.000E-01 2.000E-02 3.000E-01 1.000E-00

VB Div Plate Rep. HTS Gap HT Shield Skel. Ring Gap Vac. Vessel

1.000E-01 2.000E-01 0.000E-02 2.000E-01 4.500E-01 2.000E-02 2.500E-01

Details of the main program and 5 modules:

GENERAL:

  • Not sure what data_point_filter does. It adds zeros to the database output file from [65] to [83].

1.Aries.cpp

  • CD parameters currently in use – Pow_NB = 0.0, Pow_LH = 34.149e6, Pow_FW = 3.284e6, freq_FW = 95.0e6
  • For SiC, the radial build is set here. For DCLL, the radial build is an input file builds.data

2. Cost.cpp

  • This is where the material components are specified.

3. CostingAccount.cpp

  • Outputs detailed CostingAccounts.data

4. DesignPoint.cpp

  • Outputs detailed system design parameters in System.data

5. Geometry.cpp

  • TF coil made here

6. Part.cpp

  • Generalized formatting for finding masses, volumes, surface areas, etc for all the parts.

Costs and Volumes Calculations:

  • Unit costing is done in Cost.cpp with specific material costs for different parts
  • Costing is accounted for in CostingAccount.cpp x16, written to CA.data
  • In DesignPoint.cpp, costs are recalculated by cost.cpp x16, written to System.data
  • In Geometry.cpp, contours are used to calculate volumes, then /16 for unit volume.
  • In DP.cpp, volumes are recalculated by Geometry.cpp, then x16