SCMDetailedQuadSubChannelMeshGenerator

Creates a detailed mesh of subchannels in a square lattice arrangement

Overview

This is a mesh generator used for visualization purposes only. It is part of an input file that creates a 3D realistic mesh of the subchannels in a square lattice arrangement. This input file has a Problem block with the NoSolveProblem because it is not used to perform any calculations. It is used to create a detailed 3D mesh that the SCM solution gets projected on. The center of the mesh is the origin. Last the "side_gap" is an added distance between a perimetric pin and the duct, such that: distance(edge/corner pin center, duct wall) = pitch / 2 + side_gap.

Note that:

the vector (in the XY plane) from a corner pin center to the corner subchannel centroid is $var element = document.getElementById("moose-equation-c019d1f2-fb2e-42d3-a063-c1e8de3973b8");katex.render("Pitch/2 \\vec{i} + Pitch/2 \\vec{j}", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ .
the vector (in the XY plane) from a corner pin center to the corner of the duct is $var element = document.getElementById("moose-equation-f6e77181-c3d3-45e1-be98-1576c44e4a19");katex.render("(Pitch/2 + side\\_gap) \\vec{i} + (Pitch/2 + side\\_gap) \\vec{j}", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ .

Figure 1 presents a sketch of the SCM geometry nomenclature near the duct corner.

Figure 1: Geometric features of the quadrilateral SCM mesh

Input Parameters

heated_lengthHeated length [m]
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Heated length [m]
n_cellsThe number of cells in the axial direction
C++ Type:unsigned int
Controllable:No
Description:The number of cells in the axial direction
nxNumber of channels in the x direction [-]
C++ Type:unsigned int
Controllable:No
Description:Number of channels in the x direction [-]
nyNumber of channels in the y direction [-]
C++ Type:unsigned int
Controllable:No
Description:Number of channels in the y direction [-]
pin_diameterRod diameter [m]
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Rod diameter [m]
pitchPitch [m]
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Pitch [m]
side_gapThe side gap, not to be confused with the gap between pins, this refers to the gap next to the duct or else the distance between the subchannel centroid to the duct wall.distance(edge pin center, duct wall) = pitch / 2 + side_gap [m]
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:The side gap, not to be confused with the gap between pins, this refers to the gap next to the duct or else the distance between the subchannel centroid to the duct wall.distance(edge pin center, duct wall) = pitch / 2 + side_gap [m]

Required Parameters

block_id0Block ID used for the mesh subdomain.
Default:0
C++ Type:unsigned int
Controllable:No
Description:Block ID used for the mesh subdomain.
unheated_length_entry0Unheated length at entry [m]
Default:0
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Unheated length at entry [m]
unheated_length_exit0Unheated length at exit [m]
Default:0
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Unheated length at exit [m]

Optional Parameters

enableTrueSet the enabled status of the MooseObject.
Default:True
C++ Type:bool
Controllable:No
Description:Set the enabled status of the MooseObject.
save_with_nameKeep the mesh from this mesh generator in memory with the name specified
C++ Type:std::string
Controllable:No
Description:Keep the mesh from this mesh generator in memory with the name specified

Advanced Parameters

nemesisFalseWhether or not to output the mesh file in the nemesisformat (only if output = true)
Default:False
C++ Type:bool
Controllable:No
Description:Whether or not to output the mesh file in the nemesisformat (only if output = true)
outputFalseWhether or not to output the mesh file after generating the mesh
Default:False
C++ Type:bool
Controllable:No
Description:Whether or not to output the mesh file after generating the mesh
show_infoFalseWhether or not to show mesh info after generating the mesh (bounding box, element types, sidesets, nodesets, subdomains, etc)
Default:False
C++ Type:bool
Controllable:No
Description:Whether or not to show mesh info after generating the mesh (bounding box, element types, sidesets, nodesets, subdomains, etc)

Debugging Parameters

Input Files

(modules/subchannel/examples/mesh_generator/3d.i)
(modules/subchannel/validation/PNNL_12_pin/transient/3d.i)
(modules/subchannel/verification/enthalpy_mixing_verification/3d.i)
(modules/subchannel/test/tests/problems/coupling/3d.i)
(modules/subchannel/validation/PNNL_12_pin/steady_state/3d.i)
(modules/subchannel/examples/coupling/thermo_mech/quad/3d.i)
(modules/subchannel/test/tests/transfers/multiapp_detailed_solution_transfer/quad_viz.i)
(modules/subchannel/test/tests/mesh/detailed_quad_sub_channel_mesh/coords.i)
(modules/subchannel/verification/friction_model_verification/3d.i)
(modules/subchannel/validation/psbt/psbt_null_transient/3d.i)
(modules/subchannel/test/tests/mesh/detailed_quad_pin_mesh/coords.i)
(modules/subchannel/validation/psbt/psbt_ss/3d.i)
(modules/subchannel/validation/Blockage/PNNL_7x7/detailedMesh.i)
(modules/subchannel/examples/heating_test/3d.i)

(modules/subchannel/examples/mesh_generator/3d.i)

[Mesh]
  [subchannel]
    type = SCMDetailedQuadSubChannelMeshGenerator
    nx = 6
    ny = 6
    n_cells = 20
    pitch = 0.0126
    pin_diameter = 0.00950
    side_gap = 0.00095
    heated_length = 1.00
  []
[]

[AuxVariables]
  [mdot]
  []
  [SumWij]
  []
  [P]
  []
  [DP]
  []
  [h]
  []
  [T]
  []
  [rho]
  []
  [mu]
  []
  [S]
  []
  [w_perim]
  []
  [q_prime]
  []
[]

[Problem]
  type = NoSolveProblem
[]

[Outputs]
  exodus = true
[]

[Executioner]
  type = Steady
  nl_rel_tol = 0.9
  l_tol = 0.9
[]

(modules/subchannel/validation/PNNL_12_pin/transient/3d.i)

[Mesh]
  [subchannel]
    type = SCMDetailedQuadSubChannelMeshGenerator
    nx = 7
    ny = 3
    n_cells = 48
    pitch = 0.014605
    pin_diameter = 0.012065
    side_gap = 0.0015875
    heated_length = 1.2192
  []
[]

[AuxVariables]
  [mdot]
  []
  [SumWij]
  []
  [P]
  []
  [DP]
  []
  [h]
  []
  [T]
  []
  [rho]
  []
  [mu]
  []
  [S]
  []
  [w_perim]
  []
  [q_prime]
  []
[]

[Problem]
  type = NoSolveProblem
[]

[Outputs]
  exodus = true
[]

[Executioner]
  type = Transient
[]

(modules/subchannel/verification/enthalpy_mixing_verification/3d.i)

[Mesh]
    [subchannel]
      type = SCMDetailedQuadSubChannelMeshGenerator
      nx = 2
      ny = 1
      n_cells = 100
      pitch = 0.0126
      pin_diameter = 0.00950
      side_gap = 0.00095 #
      heated_length = 10.0
    []
  []

  [AuxVariables]
    [mdot]
    []
    [SumWij]
    []
    [P]
    []
    [DP]
    []
    [h]
    []
    [T]
    []
    [rho]
    []
    [mu]
    []
    [S]
    []
    [w_perim]
    []
    [q_prime]
    []
  []

  [Problem]
    type = NoSolveProblem
  []

  [Outputs]
    exodus = true
  []

  [Executioner]
    type = Steady
  []

(modules/subchannel/test/tests/problems/coupling/3d.i)

[GlobalParams]
    nx = 2
    ny = 2
    n_cells = 10
    pitch = 0.014605
    pin_diameter = 0.012065
    side_gap = 0.0015875
    heated_length = 1.0
[]

[Mesh]
    [subchannel]
        type = SCMDetailedQuadSubChannelMeshGenerator
    []
    [fuel_pins]
        type = SCMDetailedQuadPinMeshGenerator
        input = subchannel
    []
[]

[AuxVariables]
    [mdot]
        block = subchannel
    []
    [SumWij]
        block = subchannel
    []
    [P]
        block = subchannel
    []
    [DP]
        block = subchannel
    []
    [h]
        block = subchannel
    []
    [T]
        block = subchannel
    []
    [rho]
        block = subchannel
    []
    [mu]
        block = subchannel
    []
    [S]
        block = subchannel
    []
    [w_perim]
        block = subchannel
    []
    [q_prime]
        block = fuel_pins
    []
    [Tpin]
        block = fuel_pins
    []
    [Dpin]
        block = fuel_pins
    []
[]

[Problem]
    type = NoSolveProblem
[]

[Outputs]
    csv = true
    exodus = true
[]

[Executioner]
    type = Steady
[]

[VectorPostprocessors]
  # These samplers are used for testing purposes.
  [sample_center_pin]
    type = LineValueSampler
    variable = 'q_prime Tpin Dpin'
    sort_by = 'z'
    num_points = 11
    start_point = '0 0 0'
    end_point = '0 0 1.0'
  []
[]

(modules/subchannel/validation/PNNL_12_pin/steady_state/3d.i)

[GlobalParams]
  ######## Geometry #
  nx = 7
  ny = 3
  n_cells = 48
  pitch = 0.014605
  pin_diameter = 0.012065
  side_gap = 0.0015875
  heated_length = 1.2192
[]

[Mesh]
  [subchannel]
    type = SCMDetailedQuadSubChannelMeshGenerator
  []
  [pins]
    type = SCMDetailedQuadPinMeshGenerator
    input = subchannel
  []
[]

[AuxVariables]
  [mdot]
  []
  [SumWij]
  []
  [P]
  []
  [DP]
  []
  [h]
  []
  [T]
  []
  [Tpin]
  []
  [Dpin]
  []
  [rho]
  []
  [mu]
  []
  [S]
  []
  [w_perim]
  []
  [q_prime]
  []
[]

[Problem]
  type = NoSolveProblem
[]

[Outputs]
  exodus = true
[]

[Executioner]
  type = Steady
[]

(modules/subchannel/examples/coupling/thermo_mech/quad/3d.i)

[GlobalParams]
    nx = 2
    ny = 2
    n_cells = 25
    pitch = 0.0126
    pin_diameter = 0.00950
    side_gap = 0.00095
    heated_length = 1.0
[]

[Mesh]
    [subchannel]
        type = SCMDetailedQuadSubChannelMeshGenerator
    []
    [fuel_pins]
        type = SCMDetailedQuadPinMeshGenerator
        input = subchannel
    []
[]

[AuxVariables]
    [mdot]
        block = subchannel
    []
    [SumWij]
        block = subchannel
    []
    [P]
        block = subchannel
    []
    [DP]
        block = subchannel
    []
    [h]
        block = subchannel
    []
    [T]
        block = subchannel
    []
    [rho]
        block = subchannel
    []
    [mu]
        block = subchannel
    []
    [S]
        block = subchannel
    []
    [w_perim]
        block = subchannel
    []
    [q_prime]
        block = fuel_pins
    []
    [Tpin]
        block = fuel_pins
    []
    [Dpin]
        block = fuel_pins
    []
[]

[Problem]
    type = NoSolveProblem
[]

[Outputs]
    exodus = true
[]

[Executioner]
    type = Steady
[]

(modules/subchannel/test/tests/transfers/multiapp_detailed_solution_transfer/quad_viz.i)

[GlobalParams]
  nx = 3
  ny = 3
  n_cells = 3
  pitch = 1
  heated_length = 0.2
  pin_diameter = 0.5
[]

[Mesh]
  [sub_channel]
    type = SCMDetailedQuadSubChannelMeshGenerator
    side_gap = 0.1
  []

  [fuel_pins]
    type = SCMDetailedQuadPinMeshGenerator
    input = sub_channel
  []
[]

[AuxVariables]
  [P]
    block = sub_channel
  []
  [T]
    block = fuel_pins
  []
[]

[Problem]
  type = NoSolveProblem
[]

[Executioner]
  type = Transient
[]

[Outputs]
  exodus = true
[]

(modules/subchannel/test/tests/mesh/detailed_quad_sub_channel_mesh/coords.i)

[Mesh]
  [subchannel]
    type = SCMDetailedQuadSubChannelMeshGenerator
    nx = 3
    ny = 3
    n_cells = 10
    pitch = 0.25
    pin_diameter = 0.125
    side_gap = 0.1
    heated_length = 1
  []
[]

(modules/subchannel/verification/friction_model_verification/3d.i)

[Mesh]
    [subchannel]
      type = SCMDetailedQuadSubChannelMeshGenerator
      nx = 2
      ny = 1
      n_cells = 100
      pitch = 0.0126
      pin_diameter = 0.00950
      side_gap = 0.00095
      heated_length = 10.0
    []
  []

  [AuxVariables]
    [mdot]
    []
    [SumWij]
    []
    [P]
    []
    [DP]
    []
    [h]
    []
    [T]
    []
    [rho]
    []
    [mu]
    []
    [S]
    []
    [w_perim]
    []
    [q_prime]
    []
  []

  [Problem]
    type = NoSolveProblem
  []

  [Outputs]
    exodus = true
  []

  [Executioner]
    type = Steady
  []

(modules/subchannel/validation/psbt/psbt_null_transient/3d.i)

[Mesh]
  [subchannel]
    type = SCMDetailedQuadSubChannelMeshGenerator
    nx = 6
    ny = 6
    n_cells = 50
    pitch = 0.0126
    pin_diameter = 0.00950
    side_gap = 0.00095
    heated_length = 3.658
  []
[]

[AuxVariables]
  [mdot]
  []
  [SumWij]
  []
  [P]
  []
  [DP]
  []
  [h]
  []
  [T]
  []
  [rho]
  []
  [mu]
  []
  [S]
  []
  [w_perim]
  []
  [q_prime]
  []
[]

[Problem]
  type = NoSolveProblem
[]

[Outputs]
  exodus = true
[]

[Executioner]
  type = Steady
  nl_rel_tol = 0.9
  l_tol = 0.9
[]

(modules/subchannel/test/tests/mesh/detailed_quad_pin_mesh/coords.i)

[GlobalParams]
  nx = 3
  ny = 3
  n_cells = 2
  pitch = 1
  side_gap = 0.1
  heated_length = 1
  pin_diameter = 0.5
[]

[Mesh]
  [subchannel]
    type = SCMDetailedQuadSubChannelMeshGenerator
  []
  [pins]
    type = SCMDetailedQuadPinMeshGenerator
    input = subchannel
  []
[]

(modules/subchannel/validation/psbt/psbt_ss/3d.i)

[Mesh]
  [subchannel]
    type = SCMDetailedQuadSubChannelMeshGenerator
    nx = 6
    ny = 6
    n_cells = 20
    pitch = 0.0126
    pin_diameter = 0.00950
    side_gap = 0.00095
    heated_length = 1.0
  []

  [fuel_pins]
    type = SCMDetailedQuadPinMeshGenerator
    input = subchannel
    nx = 6
    ny = 6
    n_cells = 20
    pitch = 0.0126
    pin_diameter = 0.00950
    heated_length = 1.0
  []
[]

[AuxVariables]
  [mdot]
    block = subchannel
  []
  [SumWij]
    block = subchannel
  []
  [P]
    block = subchannel
  []
  [DP]
    block = subchannel
  []
  [h]
    block = subchannel
  []
  [T]
    block = subchannel
  []
  [rho]
    block = subchannel
  []
  [S]
    block = subchannel
  []
  [w_perim]
    block = subchannel
  []
  [mu]
    block = subchannel
  []
  [q_prime]
    block = fuel_pins
  []
  [Tpin]
    block = fuel_pins
  []
  [Dpin]
    block = fuel_pins
  []
[]

[Problem]
  type = NoSolveProblem
[]

[Outputs]
  exodus = true
[]

[Executioner]
  type = Steady
[]

(modules/subchannel/validation/Blockage/PNNL_7x7/detailedMesh.i)

[GlobalParams]
    nx = 8
    ny = 8
    n_cells = 36
    pitch = 0.0136906
    pin_diameter = 0.0099568
    side_gap = 0.0036957
    heated_length = 1.4224
  []

  [Mesh]
    [subchannel]
      type = SCMDetailedQuadSubChannelMeshGenerator
    []
  []

  [AuxVariables]
    [mdot]
    []
    [SumWij]
    []
    [P]
    []
    [DP]
    []
    [h]
    []
    [T]
    []
    [Tpin]
    []
    [Dpin]
    []
    [rho]
    []
    [mu]
    []
    [S]
    []
    [w_perim]
    []
    [q_prime]
    []
  []

  [Problem]
    type = NoSolveProblem
  []

  [Outputs]
    exodus = true
  []

  [Executioner]
    type = Steady
  []

(modules/subchannel/examples/heating_test/3d.i)

pitch = 0.0126
unheated_length_entry = 2.5
heated_length = 5.0
unheated_length_exit = 2.5
height = ${fparse unheated_length_entry + heated_length + unheated_length_exit}
n_cells = 20

[Mesh]
  [subchannel]
    type = SCMDetailedQuadSubChannelMeshGenerator
    nx = 3
    ny = 3
    n_cells = '${n_cells}'
    pitch = '${pitch}'
    pin_diameter = 0.00950
    side_gap = 0.00095
    unheated_length_entry = '${unheated_length_entry}'
    heated_length = '${heated_length}'
    unheated_length_exit = '${unheated_length_exit}'
  []
[]

[AuxVariables]
  [mdot]
  []
  [SumWij]
  []
  [P]
  []
  [DP]
  []
  [h]
  []
  [T]
  []
  [rho]
  []
  [mu]
  []
  [S]
  []
  [w_perim]
  []
  [q_prime]
  []
[]

[Problem]
  type = NoSolveProblem
[]

[Outputs]
  exodus = true
  csv = true
[]

[VectorPostprocessors]
  # These samplers are used for testing purposes.
  # Sampling over all nodes would creates megabytes of data
  # Ordering the samples by z-coordinate is consistent in parallel
  [sample_channel_edge]
    type = LineValueSampler
    variable = 'mdot SumWij P DP h T rho mu S w_perim'
    sort_by = 'z'
    num_points = ${fparse n_cells + 1}
    start_point = '${fparse pitch} 0 0'
    end_point = '${fparse pitch} 0  ${fparse height}'
  []
  [sample_channel_corner]
    type = LineValueSampler
    variable = 'mdot SumWij P DP h T rho mu S w_perim'
    sort_by = 'z'
    num_points = ${fparse n_cells + 1}
    start_point = '${fparse pitch} ${fparse pitch} 0'
    end_point = '${fparse pitch} ${fparse pitch}  ${fparse height}'
  []
[]

[Executioner]
  type = Steady
[]

Overview
Input Parameters
Input Files