Test dpa user objects (idaholab#410)

tests/userobjects/dpa/HG_CC.csv

@@ -0,0 +1,55 @@
+112.4652547,2.918285883
+122.5917992,3.484523459
+136.4022678,3.977275199
+161.6585176,4.534702104
+208.5226252,5.398668283
+247.3852619,6.730959266
+287.0859117,7.677145867
+386.623276,9.98723922
+439.3262076,11.39535849
+478.5170656,12.72398509
+630.6854106,15.83488531
+701.7347069,18.07412043
+780.5886265,20.1740281
+905.1645745,21.51766057
+1073.038882,25.08793828
+1218.688618,27.37372582
+1676.130442,35.5975299
+1946.111373,42.45768326
+2210.268258,46.32604591
+2564.974457,52.83820595
+2976.604284,60.26579547
+3381.499348,67.24293071
+4007.617911,76.66722672
+4650.764693,87.44451712
+5395.746221,97.53233334
+6260.06242,108.783905
+7111.592596,121.3781139
+8081.016324,138.491438
+9574.852126,154.411359
+11108.59731,172.2246361
+12619.65346,192.1635511
+14952.49,214.2532096
+19707.37927,266.6361974
+26526.69668,331.7043507
+35696.57584,403.530559
+48048.61301,502.0055168
+64691.32733,638.6271097
+75034.726,696.5568719
+100973.0996,847.3870883
+135808.5951,985.8106468
+182568.9871,1096.709403
+245492.1813,1247.660509
+323063.253,1357.836259
+443533.3228,1510.030623
+583383.7555,1571.531636
+800518.1909,1671.274866
+1075595.894,1778.00241
+1194933.813,1774.736155
+1639267.446,1845.660216
+2201997.569,1920.124613
+3427334.218,1992.457702
+4602698.934,2027.029058
+6052429.461,2062.95877
+8303012.666,2145.401117
+9620736.086,2139.885493

tests/userobjects/dpa/HG_CSi.csv

@@ -0,0 +1,58 @@
+247.4105339,1.076365284
+319.378375,1.370307271
+420.8331581,1.667642417
+458.4913855,1.904166376
+555.224911,2.269456705
+701.6271903,2.826396534
+986.2639098,3.928979217
+1073.970992,4.290106676
+1246.322999,4.893176951
+1415.855117,5.459673371
+1541.76515,5.961492766
+1751.93238,6.80201462
+2033.084351,7.758189638
+2410.145757,9.045452786
+2796.214326,10.08895864
+3110.424892,11.26112531
+3532.620743,12.28713788
+4099.539475,14.01436944
+4758.653031,16.34568749
+5761.172693,19.05080808
+7276.567166,22.68877208
+9386.01263,27.01151066
+10662.75277,30.1387068
+11610.9766,32.90887024
+13187.0012,35.90723085
+14976.94868,39.17877513
+17009.85602,42.74839313
+22418.98732,53.19999178
+25462.04535,58.04709708
+30176.59875,66.18256975
+34272.64197,72.21253092
+40608.19404,80.51353355
+46120.17765,87.84920339
+54631.88096,95.78276035
+62047.38023,104.509625
+71986.42362,116.5661143
+81736.67826,124.3753914
+112330.9197,151.251612
+151046.259,172.0698873
+190679.4117,195.9696719
+250994.7339,218.0954664
+278985.1864,227.6468728
+330388.8757,242.7193555
+391064.0282,247.4762533
+463118.5127,263.8616378
+548309.1644,275.1136733
+691473.748,286.5292536
+909734.5401,304.9391191
+1147267.988,317.5922779
+1477590.23,330.6488506
+2113620.201,336.0150675
+2722175.28,349.8290216
+3290618.1,356.5539466
+4062354.522,363.2745309
+5455492.122,369.577748
+6734948.395,376.5438143
+8858544.073,391.8798119
+10050693.01,391.0160938

tests/userobjects/dpa/HG_SiC.csv

@@ -0,0 +1,60 @@
+105.5578735,2.856933711
+133.3574219,3.479401561
+172.1929716,4.529702019
+204.1283529,5.28128438
+247.2589403,6.436701407
+287.1592366,7.850667644
+319.2642081,8.379699454
+362.8776814,9.777273905
+412.2384186,10.90921553
+448.7835854,11.64863621
+520.8048523,13.28611209
+591.7988771,15.15934553
+686.77139,17.29032999
+764.1390017,19.73538175
+905.1645745,21.51766057
+1029.341285,26.25415175
+1194.531069,29.94475894
+1357.364745,34.1667257
+1868.290271,47.51269996
+2168.669241,55.41653151
+2516.056957,61.80951993
+3247.938435,78.68893189
+3769.171047,89.75041809
+4374.051624,102.3668431
+5412.307079,127.5352418
+6279.276075,142.2480234
+6986.66226,162.3635318
+8103.749141,177.0915017
+9406.649367,206.551023
+10916.23846,235.5863808
+12401.12883,262.8608576
+14697.32238,299.7015857
+18967.69821,373.1131203
+22963.67991,434.8612715
+26093.99648,496.1731637
+34383.10107,603.8349438
+43438.17817,735.398318
+47288.997,785.2432148
+62342.80894,999.3158196
+70786.87703,1066.264385
+78741.19745,1190.146306
+85699.76005,1242.725225
+112952.2691,1546.560394
+125612.6312,1688.089756
+152114.6616,2011.928508
+204750.6893,2502.906379
+237730.6892,2985.252251
+269930.2991,3185.24744
+363241.1344,3874.97082
+439541.6738,4318.81944
+499075.7224,4608.156165
+671255.74,5360.914123
+847170.7114,5970.573403
+1213074.605,6634.909073
+1530201.462,7066.405703
+2056016.25,7517.66608
+2593508.312,8006.572189
+3557898.94,8326.540096
+4486873.236,8672.041959
+7285715.646,9195.353861

tests/userobjects/dpa/HG_SiSi.csv

@@ -0,0 +1,57 @@
+200.4092889,1.05645251
+227.7282779,1.205403697
+264.3418527,1.405924985
+341.1476157,1.750304579
+413.2292098,2.133228801
+533.2948916,2.655760571
+632.2012265,3.096412691
+718.3803575,3.532981624
+851.6131511,4.119185012
+967.4545353,4.596074277
+1099.33405,5.244083261
+1222.866203,5.853357207
+1389.56241,6.678632829
+1578.578788,7.451836337
+1718.959611,8.136763026
+1994.310708,9.075440175
+2266.16713,10.35500321
+2687.142859,12.34601947
+2988.333813,13.47583261
+3394.824395,15.03596645
+4024.437848,17.53078114
+4382.325528,19.14210207
+4873.522637,20.89384064
+5419.776039,22.80588489
+5900.242222,24.35165533
+7780.491165,31.69084782
+9029.112658,36.14570402
+13240.99427,51.3441192
+16723.84934,61.14885067
+20682.9724,72.85269112
+25040.31604,84.90940727
+33722.2017,110.4590921
+51578.5743,156.7891732
+59840.67476,174.8767608
+66530.98073,186.6611559
+78849.85995,212.82227
+87665.43048,227.1636937
+99590.19327,253.4630531
+115543.0224,282.703179
+131226.341,308.4605529
+172912.1313,375.391646
+227781.7886,446.7481624
+275698.8425,509.1739605
+333695.9122,580.3227499
+379087.1751,647.5082948
+555213.5677,822.5178999
+657680.9345,896.7985122
+795626.7517,977.4278077
+1139850.186,1135.840169
+1500405.27,1264.081268
+2103174.496,1374.190175
+2597087.342,1431.737308
+3488619.278,1489.501708
+4399504.52,1551.307164
+5911282.28,1650.373576
+7940515.895,1716.959003
+9398773.265,1750.608559

tests/userobjects/dpa/gold/huang_ghoniem_92_out.csv

@@ -0,0 +1,3 @@
+time,dpa_from_function,dpa_from_parkin_coulter
+0,0,0
+1,788.16396790273,273.08075993042

tests/userobjects/dpa/gold/verification_1_out.csv

@@ -0,0 +1,3 @@
+time,dpa
+0,0
+1,1183432202.0959

tests/userobjects/dpa/gold/verification_2_out.csv

@@ -0,0 +1,3 @@
+time,dpa
+0,0
+1,8.1747346200188

tests/userobjects/dpa/huang_ghoniem_92.i

@@ -0,0 +1,153 @@
+#
+# This test follows the Huang & Ghoniem "Neutron displacement damage cross sections for SiC"
+# J. of Nucl. Mat. 199, (1993), 221-230
+#
+[Mesh]
+ type = GeneratedMesh
+ dim = 1
+ xmin = 0
+ xmax = 1
+ nx = 2
+[]
+
+[Problem]
+  type = FEProblem
+  kernel_coverage_check = false
+[]
+
+[Variables]
+  [./test_var]
+  [../]
+[]
+
+[Functions]
+  [./HG_SiSi]
+    type = PiecewiseLinear
+    data_file = 'HG_SiSi.csv'
+    format = columns
+  [../]
+
+  [./HG_SiC]
+    type = PiecewiseLinear
+    data_file = 'HG_SiC.csv'
+    format = columns
+  [../]
+
+  [./HG_CSi]
+    type = PiecewiseLinear
+    data_file = 'HG_CSi.csv'
+    format = columns
+  [../]
+
+  [./HG_CC]
+    type = PiecewiseLinear
+    data_file = 'HG_CC.csv'
+    format = columns
+  [../]
+[]
+
+[VectorPostprocessors]
+  [./atomic_numbers]
+    type = ConstantVectorPostprocessor
+    value = '6 14'
+  [../]
+
+  [./mass_numbers]
+    type = ConstantVectorPostprocessor
+    value = '12 28'
+  [../]
+
+  # these are number fractions, but that does not
+  # matter for dpa
+  [./number_densities]
+    type = ConstantVectorPostprocessor
+    value = '0.5 0.5'
+  [../]
+
+  # for this problem scalar fluxes are normalized to sum to 1
+  [./scalar_fluxes]
+    type = ConstantVectorPostprocessor
+    value = '0.048138819 0.16809268 0.142519665 0.122144812 0.104682829
+             0.08971709 0.04904444 0.04213553 0.233524136'
+  [../]
+
+  [./energy_group_boundaries]
+    type = ConstantVectorPostprocessor
+    value = '3.455e6 1e6 1e5 1e4 1e3 1e2 10.0 2.4 0.625 1e-5'
+  [../]
+
+  [./cross_sections]
+    type = ConstantVectorPostprocessor
+    vector_names = 'carbon_elastic si_elastic carbon_inelastic si_inelastic'
+    value = '2.136 3.7 4.622 4.737 4.748 4.75 4.75 4.757 13.198;
+             2.938 4.205 1.980 1.918 1.952 1.956 1.957 1.959 4.034;
+             0     0     0     0     0     0     0     0     0;
+             0.10848 0   0     0     0     0     0     0     0'
+  [../]
+[]
+
+[UserObjects]
+  [./parkin_coulter]
+     type = ParkinCoulterDPAUserObject
+     damage_reaction_types = 'elastic inelastic'
+
+     irradiation_time = 1
+     Z = 'atomic_numbers'
+     Z_column = 'value'
+     A = 'mass_numbers'
+     A_column = 'value'
+     number_densities = 'number_densities'
+     number_densities_column = 'value'
+
+     scalar_flux = 'scalar_fluxes'
+     scalar_flux_column = 'value'
+     energy_group_boundaries = 'energy_group_boundaries'
+     energy_group_boundaries_column = 'value'
+     cross_section = 'cross_sections'
+     cross_section_column = 'carbon_elastic si_elastic; carbon_inelastic si_inelastic'
+     Q = '0 0; 0 -1.779e+6'
+     displacement_thresholds = '16.3 92.6'
+     logarithmic_energy_spacing = 1.25
+  [../]
+
+  [./huang_ghoniem_CSV]
+     type = FunctionDPAUserObject
+     damage_reaction_types = 'elastic inelastic'
+     damage_functions = 'HG_CC HG_CSi; HG_SiC HG_SiSi'
+     irradiation_time = 1
+     Z = 'atomic_numbers'
+     Z_column = 'value'
+     A = 'mass_numbers'
+     A_column = 'value'
+     number_densities = 'number_densities'
+     number_densities_column = 'value'
+
+     scalar_flux = 'scalar_fluxes'
+     scalar_flux_column = 'value'
+     energy_group_boundaries = 'energy_group_boundaries'
+     energy_group_boundaries_column = 'value'
+     cross_section = 'cross_sections'
+     cross_section_column = 'carbon_elastic si_elastic; carbon_inelastic si_inelastic'
+     Q = '0 0; 0 -1.779e+6'
+  [../]
+[]
+
+[Postprocessors]
+  [./dpa_from_function]
+    type = DPAPostprocessor
+    dpa_object = huang_ghoniem_CSV
+  [../]
+
+  [./dpa_from_parkin_coulter]
+    type = DPAPostprocessor
+    dpa_object = parkin_coulter
+  [../]
+[]
+
+[Executioner]
+  type = Steady
+[]
+
+[Outputs]
+  csv = true
+[]