Adding synthetic regression test.

Author: wdhawkins

framework/materials/multi_group_xs/cepxs_import.cc

@@ -243,8 +243,8 @@ PieceWiseLinearSlabMapping::MakeSurfaceFiniteElementData(size_t face_index) cons
   for (size_t qp = 0; qp < num_srf_qpoints; ++qp)
   {
     F_JxW.push_back(JxW);
-
-    F_qpoints_xyz.emplace_back(0.0, 0.0, f);
+    const auto face_xyz = (f == 0) ? grid_->vertices[v0i_] : grid_->vertices[v1i_];
+    F_qpoints_xyz.push_back(face_xyz);
   }
 
   F_num_nodes = 1;

framework/math/spatial_discretization/cell_mappings/finite_element/piecewise_linear/piecewise_linear_slab_mapping.cc

@@ -507,6 +507,15 @@ LBSProblem::GetPowerFieldFunction() const
   return field_functions_[power_gen_fieldfunc_local_handle_];
 }
 
+std::shared_ptr<FieldFunctionGridBased>
+LBSProblem::GetEnergyDepositionFieldFunction() const
+{
+  OpenSnLogicalErrorIf(not options_.energy_deposition_field_function_on,
+                       "Called when options_.energy_deposition_field_function_on == false");
+
+  return field_functions_[energy_dep_fieldfunc_local_handle_];
+}
+
 InputParameters
 LBSProblem::GetOptionsBlock()
 {
@@ -556,6 +565,11 @@ LBSProblem::GetOptionsBlock()
                               "Flag to control the creation of the power generation field "
                               "function. If set to `true` then a field function will be created "
                               "with the general name <solver_name>_power_generation`.");
+  params.AddOptionalParameter("energy_deposition_field_function_on",
+                              false,
+                              "Flag to control the creation of an energy-deposition field "
+                              "function. If set to `true` then a field function will be created "
+                              "with the general name <solver_name>_energy_deposition`.");
   params.AddOptionalParameter("power_default_kappa",
                               3.20435e-11,
                               "Default `kappa` value (Energy released per fission) to use for "
@@ -657,6 +671,9 @@ LBSProblem::ParseOptions(const InputParameters& input)
     {"power_field_function_on",
      [this](const ParameterBlock& spec)
      { options_.power_field_function_on = spec.GetValue<bool>(); }},
+    {"energy_deposition_field_function_on",
+     [this](const ParameterBlock& spec)
+     { options_.energy_deposition_field_function_on = spec.GetValue<bool>(); }},
     {"power_default_kappa",
      [this](const ParameterBlock& spec)
      { options_.power_default_kappa = spec.GetValue<double>(); }},
@@ -1217,6 +1234,27 @@ LBSProblem::InitializeFieldFunctions()
 
     power_gen_fieldfunc_local_handle_ = field_functions_.size() - 1;
   }
+
+  if (options_.energy_deposition_field_function_on)
+  {
+    std::string prefix;
+    if (options_.field_function_prefix_option == "prefix")
+    {
+      prefix = options_.field_function_prefix;
+      if (not prefix.empty())
+        prefix += "_";
+    }
+    if (options_.field_function_prefix_option == "solver_name")
+      prefix = GetName() + "_";
+
+    auto edep_ff = std::make_shared<FieldFunctionGridBased>(
+      prefix + "energy_deposition", discretization_, Unknown(UnknownType::SCALAR));
+
+    field_function_stack.push_back(edep_ff);
+    field_functions_.push_back(edep_ff);
+
+    energy_dep_fieldfunc_local_handle_ = field_functions_.size() - 1;
+  }
 }
 
 void
@@ -1310,6 +1348,39 @@ LBSProblem::UpdateFieldFunctions()
     ff_ptr->UpdateFieldVector(data_vector_local);
   }
 
+  // Update energy deposition
+  if (options_.energy_deposition_field_function_on)
+  {
+    std::vector<double> data_vector_edep_local(local_node_count_, 0.0);
+
+    for (const auto& cell : grid_->local_cells)
+    {
+      const auto& cell_mapping = sdm.GetCellMapping(cell);
+      const size_t num_nodes = cell_mapping.GetNumNodes();
+      const auto& xs = block_id_to_xs_map_.at(cell.block_id);
+      const auto& sigma_edep = xs->GetEnergyDeposition();
+
+      if (sigma_edep.empty())
+        continue;
+
+      for (size_t i = 0; i < num_nodes; ++i)
+      {
+        const auto imapA = sdm.MapDOFLocal(cell, i);
+        const auto imapB = sdm.MapDOFLocal(cell, i, phi_uk_man, 0, 0);
+
+        double nodal_edep = 0.0;
+        for (unsigned int g = 0; g < num_groups_; ++g)
+          nodal_edep += sigma_edep[g] * phi_new_local_[imapB + g];
+
+        data_vector_edep_local[imapA] = nodal_edep;
+      } // for node
+    } // for cell
+
+    const size_t ff_index = energy_dep_fieldfunc_local_handle_;
+    auto& ff_ptr = field_functions_.at(ff_index);
+    ff_ptr->UpdateFieldVector(data_vector_edep_local);
+  }
+
   // Update power generation and scalar flux
   if (options_.power_field_function_on)
   {

modules/linear_boltzmann_solvers/lbs_problem/lbs_problem.cc

@@ -256,6 +256,8 @@ class LBSProblem : public Problem
 
   /// Returns the power generation field function, if enabled.
   std::shared_ptr<FieldFunctionGridBased> GetPowerFieldFunction() const;
+  /// Returns the deposited-energy field function, if enabled.
+  std::shared_ptr<FieldFunctionGridBased> GetEnergyDepositionFieldFunction() const;
 
   bool TriggerRestartDump() const
   {
@@ -363,6 +365,7 @@ class LBSProblem : public Problem
 
   std::map<std::pair<unsigned int, unsigned int>, size_t> phi_field_functions_local_map_;
   size_t power_gen_fieldfunc_local_handle_ = 0;
+  size_t energy_dep_fieldfunc_local_handle_ = 0;
 
   /// Data carriers needed to run the sweep on GPU.
   std::shared_ptr<TotalXSCarrier> total_xs_carrier_ = nullptr;

modules/linear_boltzmann_solvers/lbs_problem/lbs_problem.h

@@ -56,6 +56,7 @@ struct LBSOptions
   double ags_tolerance = 1.0e-6;
 
   bool power_field_function_on = false;
+  bool energy_deposition_field_function_on = false;
   double power_default_kappa = 3.20435e-11; // 200MeV to Joule
   double power_normalization = -1.0;
 

modules/linear_boltzmann_solvers/lbs_problem/lbs_structs.h

@@ -345,6 +345,37 @@
       }
     ]
   },
+  {
+    "file": "transport_1d_cepxs_synthetic_3g.py",
+    "comment": "Synthetic CEPXS 3-group reader test with analytic infinite-medium solution and energy deposition",
+    "num_procs": 2,
+    "checks": [
+      {
+        "type": "KeyValuePair",
+        "key": "Phi0=",
+        "goldvalue": 0.950124,
+        "abs_tol": 0.001
+      },
+      {
+        "type": "KeyValuePair",
+        "key": "Phi1=",
+        "goldvalue": 0.933380,
+        "abs_tol": 0.001
+      },
+      {
+        "type": "KeyValuePair",
+        "key": "Phi2=",
+        "goldvalue": 0.562868,
+        "abs_tol": 0.001
+      },
+      {
+        "type": "KeyValuePair",
+        "key": "Edep=",
+        "goldvalue": 0.867650,
+        "abs_tol": 0.001
+      }
+    ]
+  },
   {
     "file": "transport_1d_3a_dsa_ortho.py",
     "comment": "1D LinearBSolver test of a block of graphite with an air cavity. DSA and TG",