From 8b938565b92a3a80ab61244298c9ea50884b8d92 Mon Sep 17 00:00:00 2001
From: Sebastian Schunert <snschune@ncsu.edu>
Date: Mon, 13 Jan 2020 14:25:52 -0700
Subject: [PATCH] Test DPAUO (#410)

---
 tests/userobjects/dpa/huang_ghoniem_92.i | 83 ++++++++++++++++++++++++
 1 file changed, 83 insertions(+)
 create mode 100644 tests/userobjects/dpa/huang_ghoniem_92.i

diff --git a/tests/userobjects/dpa/huang_ghoniem_92.i b/tests/userobjects/dpa/huang_ghoniem_92.i
new file mode 100644
index 00000000..74339fb4
--- /dev/null
+++ b/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
+[]