Refactor and regold ICs tests Refs #32847

Author: kyriv1980

modules/combined/test/tests/subchannel_thm_coupling/3D.i

@@ -10,24 +10,14 @@ heated_length = 1.0
 ###################################################
 [TriSubChannelMesh]
   [subchannel]
-    type = SCMDetailedTriSubChannelMeshGenerator
+    type = SCMDetailedTriAssemblyMeshGenerator
     nrings = ${n_rings}
     n_cells = ${n_cells}
     flat_to_flat = ${inner_duct_in}
     heated_length = ${heated_length}
     pin_diameter = ${fuel_pin_diameter}
     pitch = ${fuel_pin_pitch}
   []
-
-  [fuel_pins]
-    type = SCMDetailedTriPinMeshGenerator
-    input = subchannel
-    nrings = ${n_rings}
-    n_cells = ${n_cells}
-    pin_diameter = ${fuel_pin_diameter}
-    heated_length = ${heated_length}
-    pitch = ${fuel_pin_pitch}
-  []
 []
 
 [AuxVariables]

modules/combined/test/tests/subchannel_thm_coupling/subchannel.i

@@ -17,7 +17,7 @@ heated_length = 1.0
 ###################################################
 [TriSubChannelMesh]
   [subchannel]
-    type = SCMTriSubChannelMeshGenerator
+    type = SCMTriAssemblyMeshGenerator
     nrings = ${n_rings}
     n_cells = ${n_cells}
     flat_to_flat = ${inner_duct_in}
@@ -29,15 +29,6 @@ heated_length = 1.0
     spacer_z = '0.0'
     spacer_k = '0.0'
   []
-
-  [fuel_pins]
-    type = SCMTriPinMeshGenerator
-    input = subchannel
-    nrings = ${n_rings}
-    n_cells = ${n_cells}
-    heated_length = ${heated_length}
-    pitch = ${fuel_pin_pitch}
-  []
 []
 
 [AuxVariables]
@@ -273,12 +264,14 @@ heated_length = 1.0
 [Transfers]
   [subchannel_transfer]
     type = SCMSolutionTransfer
+    transfer_type = subchannel
     to_multi_app = viz
     variable = 'mdot SumWij P DP h T rho mu S'
   []
   [pin_transfer]
-    type = SCMPinSolutionTransfer
+    type = SCMSolutionTransfer
+    transfer_type = pin
     to_multi_app = viz
-    variable = 'Dpin Tpin q_prime'
+    variable = 'Tpin q_prime Dpin'
   []
 []

modules/subchannel/doc/content/source/meshgenerators/SCMQuadAssemblyMeshGenerator.md

@@ -0,0 +1,21 @@
+# SCMQuadAssemblyMeshGenerator
+
+!syntax description /Mesh/SCMQuadAssemblyMeshGenerator
+
+## Overview
+
+!! Intentional comment to provide extra spacing
+
+This mesh generator creates the 1D subchannel and pin meshes for a square lattice arrangement.
+The generated subdomains are named `subchannel` and `fuel_pins`.
+The center of the mesh is the origin.
+
+## Example Input File Syntax
+
+!listing /test/tests/problems/psbt/psbt_explicit.i block=QuadSubChannelMesh language=moose
+
+!syntax parameters /Mesh/SCMQuadAssemblyMeshGenerator
+
+!syntax inputs /Mesh/SCMQuadAssemblyMeshGenerator
+
+!syntax children /Mesh/SCMQuadAssemblyMeshGenerator

modules/subchannel/doc/content/source/postprocessors/SCMPinPowerPostprocessor.md

@@ -10,8 +10,8 @@ The user needs to specify a subchannel problem. Either a [QuadSubChannel1PhasePr
 
 ## Example Input File Syntax
 
-!listing /test/tests/SCMQuadPower/test.i block=Postprocessors language=moose
-!listing /test/tests/SCMTriPower/test.i block=Postprocessors language=moose
+!listing /test/tests/SCMQuadPower/test_with_pins.i block=Postprocessors language=moose
+!listing /test/tests/SCMTriPower/test_with_pins.i block=Postprocessors language=moose
 
 !syntax parameters /Postprocessors/SCMPinPowerPostprocessor
 

modules/subchannel/include/mesh/TriSubChannelMesh.h

@@ -41,7 +41,7 @@ class TriSubChannelMesh : public SubChannelMesh
   const Real & getDuctToPinGap() const { return _duct_to_pin_gap; }
 
   /**
-   * Return the number of rings
+   * Return the number of fuel-pin rings, counting the center pin as the first ring
    */
   const unsigned int & getNumOfRings() const { return _n_rings; }
 
@@ -131,7 +131,7 @@ class TriSubChannelMesh : public SubChannelMesh
   unsigned int pinIndex(const Point & p) const override;
 
 protected:
-  /// number of rings of fuel pins
+  /// number of fuel-pin rings, counting the center pin as the first ring
   unsigned int _n_rings;
   /// number of subchannels
   unsigned int _n_channels;

modules/subchannel/src/mesh/TriSubChannelMesh.C

@@ -82,17 +82,15 @@ TriSubChannelMesh::channelIndex(const Point & p) const
   // subchannel
   Real distance_outer_ring = _flat_to_flat / 2 - _duct_to_pin_gap - _pin_diameter / 2;
   Real channel_distance = std::sqrt(std::pow(p(0), 2) + std::pow(p(1), 2));
-  Real angle = std::abs(std::atan(p(1) / p(0)));
+  Real angle = std::abs(std::atan2(p(1), p(0)));
   Real projection_angle =
       angle - libMesh::pi / 6 - std::trunc(angle / (libMesh::pi / 3)) * (libMesh::pi / 3);
   channel_distance = channel_distance * std::cos(projection_angle);
 
   // Projecting point on top edge to determine if the point is a corner or edge subchannel by x
   // coordinate
-  Real loc_angle = std::atan(p(1) / p(0));
-  if (p(0) <= 0)
-    loc_angle += libMesh::pi;
-  else if (p(0) >= 0 && p(1) <= 0)
+  Real loc_angle = std::atan2(p(1), p(0));
+  if (loc_angle < 0.0)
     loc_angle += 2 * libMesh::pi;
   Real rem_ang = std::trunc(loc_angle / (libMesh::pi / 3)) * (libMesh::pi / 3) - libMesh::pi / 3;
   Real x_coord_new = (std::cos(-rem_ang) * p(0) - std::sin(-rem_ang) * p(1));
@@ -101,46 +99,36 @@ TriSubChannelMesh::channelIndex(const Point & p) const
   // looping over all channels
   for (unsigned int i = 0; i < _n_channels; i++)
   {
+    // Distance from the point to subchannel
+    distance1 = std::sqrt(std::pow((p(0) - _subchannel_position[i][0]), 2.0) +
+                          std::pow((p(1) - _subchannel_position[i][1]), 2.0));
 
-    if (_n_rings == 1)
+    // If subchannel belongs to center ring
+    if (channel_distance < distance_outer_ring)
     {
-      Real angle = std::atan(p(1) / p(0));
-      if ((i * libMesh::pi / 6.0 < angle) && (angle <= (i + 1) * libMesh::pi / 6.0))
-        return i;
-    }
+      if ((distance1 < distance0) && (_subch_type[i] == EChannelType::CENTER))
+      {
+        j = i;
+        distance0 = distance1;
+      } // if
+    } // if
+    // If subchannel belongs to outer ring
     else
     {
-      // Distance from the point to subchannel
-      distance1 = std::sqrt(std::pow((p(0) - _subchannel_position[i][0]), 2.0) +
-                            std::pow((p(1) - _subchannel_position[i][1]), 2.0));
-
-      // If subchannel belongs to center ring
-      if (channel_distance < distance_outer_ring)
+      if ((distance1 < distance0) &&
+          (_subch_type[i] == EChannelType::EDGE || _subch_type[i] == EChannelType::CORNER))
       {
-        if ((distance1 < distance0) && (_subch_type[i] == EChannelType::CENTER))
+        if (((x_coord_new > x_lim) || (x_coord_new < -x_lim)) &&
+            _subch_type[i] == EChannelType::CORNER)
         {
           j = i;
           distance0 = distance1;
         } // if
-      }   // if
-      // If subchannel belongs to outer ring
-      else
-      {
-        if ((distance1 < distance0) &&
-            (_subch_type[i] == EChannelType::EDGE || _subch_type[i] == EChannelType::CORNER))
+        else if ((x_coord_new <= x_lim && x_coord_new >= -x_lim) &&
+                 _subch_type[i] == EChannelType::EDGE)
         {
-          if (((x_coord_new > x_lim) || (x_coord_new < -x_lim)) &&
-              _subch_type[i] == EChannelType::CORNER)
-          {
-            j = i;
-            distance0 = distance1;
-          } // if
-          else if (((x_coord_new > x_lim) || (x_coord_new > -x_lim)) &&
-                   _subch_type[i] == EChannelType::EDGE)
-          {
-            j = i;
-            distance0 = distance1;
-          }
+          j = i;
+          distance0 = distance1;
         }
       }
     }
@@ -205,7 +193,7 @@ TriSubChannelMesh::pinPositions(std::vector<Point> & positions,
   Real theta_corrected = 0.0;
   Real pi = libMesh::pi;
   unsigned int k = 0;
-  positions.emplace_back(0.0, 0.0);
+  positions.emplace_back(center(0), center(1));
   for (unsigned int i = 1; i < nrings; i++)
   {
     dtheta = 2.0 * pi / (i * 6);

modules/subchannel/src/meshgenerators/SCMDetailedTriAssemblyMeshGenerator.C

@@ -27,7 +27,9 @@ SCMDetailedTriAssemblyMeshGenerator::validParams()
   params.addParam<Real>("unheated_length_entry", 0.0, "Unheated length at entry [m]");
   params.addRequiredParam<Real>("heated_length", "Heated length [m]");
   params.addParam<Real>("unheated_length_exit", 0.0, "Unheated length at exit [m]");
-  params.addRequiredParam<unsigned int>("nrings", "Number of fuel Pin rings per assembly [-]");
+  params.addRequiredParam<unsigned int>(
+      "nrings",
+      "Number of fuel-pin rings per assembly, counting the center pin as the first ring [-]");
   params.addRequiredParam<Real>("flat_to_flat",
                                 "Flat to flat distance for the hexagonal assembly [m]");
   params.addRequiredParam<unsigned int>("n_cells", "The number of cells in the axial direction");
@@ -63,6 +65,11 @@ SCMDetailedTriAssemblyMeshGenerator::SCMDetailedTriAssemblyMeshGenerator(
     _verbose(getParam<bool>("verbose_flag")),
     _elem_id(0)
 {
+  if (_n_rings < 2)
+    mooseError(name(),
+               ": 'nrings' must be at least 2. In this mesh generator, the center pin counts as "
+               "the first ring, so a 7-pin bundle uses nrings = 2.");
+
   Real L = _unheated_length_entry + _heated_length + _unheated_length_exit;
   Real dz = L / _n_cells;
   for (unsigned int i = 0; i < _n_cells + 1; i++)
@@ -293,18 +300,6 @@ SCMDetailedTriAssemblyMeshGenerator::SCMDetailedTriAssemblyMeshGenerator(
       _subchannel_position[i][0] = (a2 * x1 - a1 * x0) / (a2 - a1);
       _subchannel_position[i][1] = (a2 * y1 - a1 * y0) / (a2 - a1);
     }
-
-    /// Special case _n_rings == 1
-    if (_n_rings == 1)
-    {
-      for (unsigned int i = 0; i < _n_channels; i++)
-      {
-        Real angle = (2 * i + 1) * libMesh::pi / 6.0;
-        _subch_type[i] = EChannelType::CORNER;
-        _subchannel_position[i][0] = std::cos(angle) * _flat_to_flat / 2.0;
-        _subchannel_position[i][1] = std::sin(angle) * _flat_to_flat / 2.0;
-      }
-    }
   }
 }
 
@@ -399,19 +394,9 @@ SCMDetailedTriAssemblyMeshGenerator::generate()
   }
 
   // specify number and type of sub-channel
-  unsigned int n_center, n_side, n_corner;
-  if (_n_rings == 1)
-  {
-    n_corner = 6;
-    n_side = 0;
-    n_center = _n_channels - n_side - n_corner;
-  }
-  else
-  {
-    n_corner = 6;
-    n_side = (_n_rings - 1) * 6;
-    n_center = _n_channels - n_side - n_corner;
-  }
+  unsigned int n_corner = 6;
+  unsigned int n_side = (_n_rings - 1) * 6;
+  unsigned int n_center = _n_channels - n_side - n_corner;
   if (_verbose)
   {
     _console << "Centers: " << n_center << std::endl;

modules/subchannel/src/meshgenerators/SCMTriAssemblyMeshGenerator.C

@@ -27,7 +27,9 @@ SCMTriAssemblyMeshGenerator::validParams()
   params.addParam<Real>("unheated_length_entry", 0.0, "Unheated length at entry [m]");
   params.addRequiredParam<Real>("heated_length", "Heated length [m]");
   params.addParam<Real>("unheated_length_exit", 0.0, "Unheated length at exit [m]");
-  params.addRequiredParam<unsigned int>("nrings", "Number of fuel Pin rings per assembly [-]");
+  params.addRequiredParam<unsigned int>(
+      "nrings",
+      "Number of fuel-pin rings per assembly, counting the center pin as the first ring [-]");
   params.addRequiredParam<Real>("flat_to_flat",
                                 "Flat to flat distance for the hexagonal assembly [m]");
   params.addRequiredParam<Real>("dwire", "Wire diameter [m]");
@@ -86,6 +88,11 @@ SCMTriAssemblyMeshGenerator::SCMTriAssemblyMeshGenerator(const InputParameters &
     _n_gaps(0),
     _elem_id(0)
 {
+  if (_n_rings < 2)
+    mooseError(name(),
+               ": 'nrings' must be at least 2. In this mesh generator, the center pin counts as "
+               "the first ring, so a 7-pin bundle uses nrings = 2.");
+
   if (_spacer_z.size() != _spacer_k.size())
     mooseError(name(), ": Size of vector spacer_z should be equal to size of vector spacer_k");
 
@@ -739,18 +746,6 @@ SCMTriAssemblyMeshGenerator::SCMTriAssemblyMeshGenerator(const InputParameters &
     }
   }
 
-  /// Special case _n_rings == 1
-  if (_n_rings == 1)
-  {
-    for (unsigned int i = 0; i < _n_channels; i++)
-    {
-      Real angle = (2 * i + 1) * libMesh::pi / 6.0;
-      _subch_type[i] = EChannelType::CORNER;
-      _subchannel_position[i][0] = std::cos(angle) * _flat_to_flat / 2.0;
-      _subchannel_position[i][1] = std::sin(angle) * _flat_to_flat / 2.0;
-    }
-  }
-
   // Reduce reserved memory in the channel-to-gap map.
   for (auto & gap : _chan_to_gap_map)
   {

modules/subchannel/test/tests/ics/FCTFdisplacementIC/gold/test_out.csv

@@ -0,0 +1,3 @@
+time,center_S,center_w_perim,corner_S,corner_w_perim,edge_S,edge_w_perim
+0,0,0,0,0,0,0
+1,1.78771169811e-05,0.017645070068785,1.6955350355555e-05,0.014091834788374,4.3125203206028e-05,0.028871590068785