Test DPAUO (idaholab#410)

tests/userobjects/dpa/huang_ghoniem_92.i

@@ -0,0 +1,83 @@
+#
+# This test follows the Huang & Ghoniem "Neutron displacement damage cross sections for SiC"
+# J. of Nucl. Mat. 199, (1993), 221-230
+# Results (with extended energy range) from this test case are included in the writeup
+#
+[Mesh]
+ type = GeneratedMesh
+ dim = 1
+ xmin = 0
+ xmax = 1
+ nx = 2
+[]
+
+[Problem]
+  type = FEProblem
+  kernel_coverage_check = false
+[]
+
+[Variables]
+  [./test_var]
+  [../]
+[]
+
+[VectorPostprocessors]
+  [./atomic_numbers]
+    type = ConstantVectorPostprocessor
+    value = '6 14'
+  [../]
+
+  [./mass_numbers]
+    type = ConstantVectorPostprocessor
+    value = '12 28'
+  [../]
+
+  [./number_densities]
+    type = ConstantVectorPostprocessor
+    value = '0.66 0.66'
+  [../]
+
+  [./scalar_fluxes]
+    type = ConstantVectorPostprocessor
+    value = '1 2 3'
+  [../]
+
+  [./energy_group_boundaries]
+    type = ConstantVectorPostprocessor
+    value = '2e6 20 0.625 0'
+  [../]
+
+  [./cross_sections]
+    type = ConstantVectorPostprocessor
+    vector_names = 'carbon_elastic si_elastic'
+    value = '0.5 0.8 1 ; 0.2 0.4 1.2'
+  [../]
+[]
+
+[UserObjects]
+  [./parkin_coulter]
+     type = DPAUserObject
+     damage_reaction_types = 'elastic'
+
+     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'
+
+     displacement_thresholds = '16.3 92.6'
+     logarithmic_energy_spacing = 1.1
+  [../]
+[]
+
+[Executioner]
+  type = Steady
+[]