progress of aries systems code development zoran dragojlovic a. rené raffray farrokh najmabadi...
TRANSCRIPT
Progress of ARIES Systems Code Development
Zoran Dragojlovic A. René Raffray
Farrokh Najmabadi
ARIES-“TNS” Project MeetingJune 14 and 15, 2007
General Atomics, La Jolla CA
Progress Overview• Previous Work
• Realistic 3-D geometry of ARIES-AT Tokamak generated, while using plasma parameters and inboard radial builds as input data.
• Power flow incorporated, with the input based on an external physics code provided by Chuck Kessel.
• Tested geometry and power flow scheme by comparison with ARIES-AT data. As a test run, generated geometries and power flows for 120 physics operating points, which had passed the engineering filters.
• New Accomplishments• Adopted a new object-oriented algorithm and replaced the previously written
Fortran code with C++. The new algorithm allows modular programming, where generic building blocks can be re-used many times and in different combinations.
• Costing algorithms included. • Achieved a reasonable agreement with previously published ARIES-AT data.• Included Les Waganer’s updates for accounts 20 and 21.
• Cost of Electricity obtained.• Agrees within 2-3% with the published ARIES-AT value.• Generated contour maps of COE for 1,452 data points provided by Chuck
Kessel. The contour maps are generated in several different 2-D planes in order to demonstrate potential control knobs for the cost of electricity.
• Engineering Algorithms are in Progress• Include the engineering limits in order to calculate radial and axial builds.
Currently inboard radial builds are taken as an input and remaining builds and outboard/axial builds are scaled to match ARIES-AT.
• Incorporate the magnet algorithms in collaboration with Leslie Bromberg.
ARIES Systems Code• Foundation for the algorithm is a
general-purpose systems analysis toolbox.
– Consists of ready-to-use generic objects (classes) that serve as building blocks for different systems algorithms with different objectives.
– Class DesignPoint holds design-specific data that describe the entire machine, such as plasma parameters, builds, power flow, building volumes, etc. These data are accessed, operated on and displayed by special functions that belong to the same class.
– Class Part holds part-specific data such as contours, areas, volumes, etc.
– Class CostingAccount holds the costing account structure for the selected machine design.
– Simple declaration statement such as “Part Blanket_II;” declares all the arrays and single-valued variables needed to define this particular object.
• Systems Code is generated by connecting the elements of the toolbox together.
Class DesignPoint{data;
functions that operate on data;};
Class Part{data;
functions;};
Class CostingAccount
{data;
Functions; };
Systems Analysis Algorithm
• An example configuration that calculates the cost of electricity for ARIES-AT is simple and straightforward.
• Once the objects are well defined, they can be used to build many different codes.
Input Dataphysics parametersradial builds
DesignPoint Aries_AT;Aries_AT.get_physics();Aries_AT.get_builds();
Geometry
Part First_Wall_Inboard;
Part Blanket_II;
Part PFCoil[29];
Display OutputAries_AT.show();
Costing Algorithms
COE
1
2
3
4
Different Arbitrary Configurations are Easy to Assemble
• Configuration 1: different design points can be run together, at the same time, in case they need to share/exchange any information during the run.
• Configuration 2: different parts or costing accounts can be plugged in and out as the code runs.
Input Data
Starfire
ARIES SPPS
ARIES AT
Prometheus-L
Geom.
Costing
Geom.Costing
Geom.Costing
Geom.Costing
Configuration 1
ARIES-TNS
Geometry
Blanket A
Blanket B
Blanket C
Configuration 2
Costing Algorithms
Output
Costing Algorithm
• At recommendation from Ronald Miller, forwarded by Laila El-Guebaly, costing algorithms were based on an unpublished Chapter 2: System Study of the ARIES II-IV report. A copy is available upon request.
• After the algorithms were implemented, the new code was validated by comparing the results to the data published in the ARIES-AT design book. The cost base for these data is from 1992 and the case with LSA=4 was selected.
• New values of Total Direct Cost (TDC) and Total Capital Cost (TCC) differ from published data by 10.7% each . Cost of Electricity (COE) estimated by the new code differs from the previous data by 1.64%.
Side-By-Side Comparison of Costing Accounts for ARIES-AT
Costing Accounts New DataPrevious
DataRelative Error
LSA = 4 [M$, 1992] [M$, 1992] [%]
20. Land and Land Rights 10.44 10.589 1.42
21.0 Structures and Site Facilities 335.458 374.226 10.9321.1 Site Improvements and Facilities 18.4421.2 Reactor Building 127.90221.3 Turbine Building 33.844121.4 Cooling Structures 12.231721.5 Power Supply and Energy Storage Building 14.9821.6 Miscelaneous Buildings 125.121.7 Ventilation Stack 2.96
22.0 Reactor Plant Equipment (RPE) 849.538 940.939 10.21
22.1.0 Reactor Equipment 519.599 549.989 5.6822.1.1 First Wall/Blanket/Reflector 79.3427 71.422 10.5122.1.2 Shield 82.9521 77.121 7.2922.1.3 Magnet Coils 154.39 140.762 9.2322.1.4 Suplemental Heating Systems 54.6586 43.6 22.5122.1.5 Primary Structure and Support 31.6866 31.686 0.0022.1.6 Reactor Vacuum System 98.3727 116.203 16.6222.1.7 Power Supply 9.26 59.702 146.2922.1.8 Impurity Control System 6.3366 4.817 27.2522.1.9 Direct Energy Conversion 0 022.1.10 ECRH Breakdown System 2.6 4.677 57.08
Questions:1. Account 22.1.5
Primary Structure and Support uses Vstr volume which I am not sure where to get from.
2. Account 22.1.7 Power Supply uses a formula in terms of $/kVA. Where does the “kVA” come from?
Side-By-Side Comparison of Costing Accounts for ARIES-AT
Costing Accounts New DataPrevious
DataRelative Error
LSA = 4 [M$, 1992] [M$, 1992] [%]
22.2.0 Main Heat-Transfer System 255.663 209.947 19.6422.2.1 Primary Coolant 171.10122.2.2 Intermediate Coolant System 29.699422.2.3 Secondary Coolant System 54.8634
22.3 Auxilary Cooling Systems 2.1806722.4 Radioactive Waste Treatment 3.88556
22.5.0 Fuel Handling and Storage 26.3922.5.1 Pellet Injectors 12.1422.5.2 Fuel Processing System 022.5.3 Fuel Storage 6.0722.5.4 Atmospheric Tritium Recovery 022.5.5 Water Detritiation System 8.18
22.6 Other Reactor Plant Equipment 3.5287322.7 Instrumentation and Control 38.29
23.0 Turbine Plant Equipment 235.242 243.034 3.26
24.0 Electric Plant Equipment 124.822 131.34 5.09
25.0 Miscellaneous Plant Equipment 59.592 55.71 6.73
26.0 Special Materials 91.4584 83.766 8.78
3. Account 22.2.0 Main Heat Transfer System uses straightforward formulas which are dependent on thermal power PTH. No idea where the discrepancy comes from.
Side-By-Side Comparison of Costing Accounts for ARIES-AT
Costing Accounts New DataPrevious
DataRelative Error
LSA = 4 [M$, 1992] [M$, 1992] [%]
90.0 Total Direct Cost (TDC) 1674.5 1862.92 10.65
91.0 Construction Services and Equipment 252.853 281.301 10.65
92.0 Home Office Engineering and Services 87.0753 96.872 10.65
93.0 Field Office Engineering and Services 145.684 162.074 10.65
94.0 Owner's Cost 324.02 360.475 10.65
95.0 Process Contingency 0 0
96.0 Project Contingency 484.411 523.108 7.68
97.0 Interest During Construction 472.383 542.971 13.90
98.0 Escalation During Construction 0 0
99.0 Total Capital Cost (TCC) 3441 3829.72 10.69
Cost of Electricity [mill/kVeh] 65.628 64.56 1.64
•Accounts 90 through 99 are proportional to the sum of the costs defined by the accounts 20 – 26, therefore they propagate the same cumulative error of 10.69%.
•With the few exceptions, the new code provided a reasonable match with the existing ARIES-AT data.
Comparison Between ARIES-AT and ARIES-TNS for Cost Base of 2007
Costing Accounts ARIES-AT ARIES-TNS Relative Error
LSA = 4 [M$, 2007] [M$, 2007] [%]
20. Land and Land Rights 15.315 11 32.79
21.0 Structures and Site Facilities 445.121 383.992 14.7521.1 Site Improvements and Facilities 27.0506 25 7.8821.2 Reactor Building 140.646 132.684 5.8321.3 Turbine Building 49.6478 76.6075 42.7121.4 Heat Rejection Structures 17.9433 11.1325 46.8521.5 Electrical Equipment and Power Supply Bldg. 21.975 21.23 3.4521.6 Plant Auxiliary Systems Building 20.831
21.7 Hot Cell Building 54.0026
21.8 Reactor Service Building 4.16
21.9 Service Water Building 1.46
21.10 Fuel Handling and Storage Building 15.1035
21.11 Control Room Building 6.8
21.12 On-Site AC Power Supply Building 4.5
21.13 Administrative Building 1.9
21.14 Site Service Building 1.9
21.15 Cryogenic and Inert Gas Storage Building 2
21.16 Security Building 0.68
21.17 Ventilation Stack 4.34218 4 8.2021.x Miscelaneous Buildings 183.516
• Accounts 20 and 21 for ARIES-TNS were updated by Les Waganer, who will provide a more detailed information in his talk.
Comparison Between ARIES-AT and ARIES-TNS for Cost Base of 2007
Costing Accounts ARIES-AT ARIES-TNS Relative Error
LSA = 4 [M$, 2007] [M$, 2007] [%]
22.0 Reactor Plant Equipment (RPE) 1326.71 1326.71 0.00
23.0 Turbine Plant Equipment 345.089 345.089 0.00
24.0 Electric Plant Equipment 183.109 183.109 0.00
25.0 Miscellaneous Plant Equipment 87.4187 87.4187 0.00
26.0 Special Materials 134.165 134.165 0.00
90.0 Total Direct Cost (TDC) 2536.9 2471.49 2.61
91.0 Construction Services and Equipment 383.077 373.195 2.61
92.0 Home Office Engineering and Services 131.92 128.517 2.61
93.0 Field Office Engineering and Services 220.713 215.019 2.61
94.0 Owner's Cost 490.896 478.233 2.61
95.0 Process Contingency 0 0
96.0 Project Contingency 733.89 714.958 2.61
97.0 Interest During Construction 715.669 697.207 2.61
98.0 Escalation During Construction 0 0
99.0 Total Capital Cost (TCC) 5213.1 5078.62 2.61
Cost of Electricity [mill/kVeh] 91.821 89.785 2.24
New accounts 20 and 21 did not significantly impact the cost of electricity.
Contour Maps of Cost of Electricity Versus Several Selected Parameters
• The data for this test was provided by Charles Kessel, PPPL.– Generated a large number of physics operating points close to ARIES-
AT specs.– All the viable points were filtered through engineering limits, which are:
• First wall heat flux. • Divertor peak heat flux. • TF and PF coils
– peak field limit and superconducting current density limit. • Bucking cylinder criteria.
– This resulted in an inboard radial build with assumed thicknesses for the First Wall, Blanket, Shield, and Vacuum Vessel from ARIES-AT neutronics analysis.
• Final outcome: 1,452 surviving data points were provided as an input to the new ARIES Systems Code.
• The new systems code estimated the cost of electricity for all the data points in 2 hours of CPU time. The data was used to generate several contour maps.
• COE was plotted against several pairs of parameters, at Chuck’s suggestion:
– Toroidal field BT and plasma major radius R.– Bootstrap fraction and n.– Power density and plasma major radius R.
Cost of Electricity Versus Toroidal Field and Plasma Major Radius
• For each pair of coordinates (BT, R), there was a small variation (< 8 mill/kWeh) of COE within the same “geometrical” point. The lower limit of that variation was shown on the left, while the upper limit was shown on the right.
• COE increases with toroidal field and plasma major radius, as expected.
• The dependence of COE on BT is higher for larger radius R.
65 70 75 80 85 [mill/kWeh]
Lower Limit of COE Upper Limit of COE
(5.5, 6.0) (5.5, 6.0)
Cost of Electricity Versus Bootstrap Fraction and n
• Data implies that the cost of electricity depends only on the bootstrap fraction.
65 70 75 80 85 [mill/kWeh]
Lower Limit of COE Upper Limit of COE
Cost of Electricity Versus Power Density Pdiv,cond/R and Plasma Major Radius
• Power density is based on the power conducted in the divertor divided by the plasma major radius.• In the space of these two parameters, the COE behaves as expected.
65 70 75 80 85 [mill/kWeh]
Lower Limit of COE Upper Limit of COE
Cost of Electricity Versus Power Density Pdiv,cond/R2 and Plasma Major Radius
• Power density is based on the power conducted in the divertor divided by the plasma major radius squared.• P/R2 tends to over-emphasize the effect of radius, compared to P/R shown in the previous slide.
65 70 75 80 85 [mill/kWeh]
Lower Limit of COE Upper Limit of COE
Local Variation of COE Across the Planes of Different Parameters
• Local variation of COE is obtained by subtracting the lower from the upper limit at each “geometrical” point.
• In the plane defined by power density Pdiv,cond/R and plasma major radius R, variation of COE at each point is almost negligible, which indicates that these two parameters might be the strongest knobs of impact on the COE.
8
7
6
5
4
3
2
1
0
[mill/kWeh]
Action ItemsDeliverables Criteria of Success Duration Date
ExpectedUrgency
Engineering Algorithms
• Implement the engineering criteria in order arrive at realistic radial and axial builds.• Include the new magnet algorithms from Leslie Bromberg.
1 - 1.5 months
July 30, 2007
HIGH
Other Updates
Of Systems Algorithms +
Tests
•Costing•Physics•…etc …
up to 1 month
Aug. 30, 2007
HIGH
Code Documentation
•Provide a web-based manual. 1-2 months
Sept. 15, 2007
MEDIUM
Presentation of Data
•Interactive COE maps. 2-4 weeks
July 30, 2007
MEDIUM
Ready for TNS studies
August 30, 2007
Discussion and Conclusions• ARIES Systems Code development is on the schedule.
– Costing algorithms are fully implemented by June 2007, as promised at the last meeting.
– Cost of electricity is now within few percents from the existing ARIES-AT data. More tests are possible and can be accomplished within a month, if needed.
– Engineering algorithms will be implemented in order to make the code ready for the ARIES-TNS study by August 30, possibly earlier.
• Object-oriented programming style was adopted since the last meeting, in order to make the code modular and more flexible with respect to different objectives of the system study.
• Data visualization is undergoing development, in order to provide an insight into trends of the cost of electricity and help to optimize the power plant.– COE contour maps were plotted against several different parameters,
such as major plasma radius, power density, bootstrap fraction, toroidal field and n.
– Data included in this analysis imply that there could be a pair of strong knobs that control the COE, such as plasma major radius and power density, for example. A better resolution and more accurate calculations will be done to test this finding.
Acknowledgement to Co-Workers – Thanks!Who Help Provided When
Laila El-Guebaly • Shield thickness scaling with the average neutron wall load.• TF Coil thickness and material densities.• Costing documentation.
May 5, 2007
May 11, 2007
May 17, 2007
Les Waganer • Correspondence about costing details, approach and background.• Updated accounts 20 and 21.• Commented on new results, answered questions about costing accounts.
4/3 – 5/22, 2007
May 30, 2007
June 5, 2007
Leslie Bromberg • Costing data sheets for magnets made of different materials.
May 15, 2007
Chuck Kessel • Provided data points for COE maps. May 31, 2007
Xueren Wang • Provided tables of gross efficiency and pumping power as functions of neutron wall load and surface heat flux
May, June 2007
Tak-Kuen Mau • Provided discussion and offered help with divertor heat transfer modeling.• Explained how current drive works.
May 2007