SCMDetailedQuadAssemblyMeshGenerator

Creates a detailed mesh of subchannels and fuel pins in a square lattice arrangement

Overview

This mesh generator creates one detailed 3D visualization mesh containing both the subchannel volume and the fuel pin volume for a square lattice arrangement. The generated subdomains are named subchannel and fuel_pins.

This object replaces the old chained SCMDetailedQuadSubChannelMeshGenerator and SCMDetailedQuadPinMeshGenerator input pattern. See the SubChannel mesh generator migration page for old-to-new input examples.

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

num_sectors16Number of azimuthal sectors used to discretize each circular pin cross section.
Default:16
C++ Type:unsigned int
Range:num_sectors>=4
Controllable:No
Description:Number of azimuthal sectors used to discretize each circular pin cross section.
pin_block_id1Fuel pin block id.
Default:1
C++ Type:unsigned int
Controllable:No
Description:Fuel pin block id.
subchannel_block_id0Subchannel block id.
Default:0
C++ Type:unsigned int
Controllable:No
Description:Subchannel block id.
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/verification/friction_model_verification/3d.i)
(modules/subchannel/validation/psbt/psbt_null_transient/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/verification/enthalpy_mixing_verification/3d.i)
(modules/subchannel/examples/heating_test/3d.i)
(modules/subchannel/validation/Blockage/PNNL_7x7/detailedMesh.i)
(modules/subchannel/validation/PNNL_12_pin/steady_state/3d.i)
(modules/subchannel/test/tests/problems/coupling/3d.i)
(modules/subchannel/validation/PNNL_12_pin/transient/3d.i)
(modules/subchannel/test/tests/mesh/detailed_quad_assembly/coords.i)
(modules/subchannel/validation/psbt/psbt_ss/3d.i)
(modules/subchannel/examples/mesh_generator/3d.i)

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

[Mesh]
    [assembly]
      type = SCMDetailedQuadAssemblyMeshGenerator
      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]
  [assembly]
    type = SCMDetailedQuadAssemblyMeshGenerator
    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/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]
    [assembly]
        type = SCMDetailedQuadAssemblyMeshGenerator
    []
[]

[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]
  [assembly]
    type = SCMDetailedQuadAssemblyMeshGenerator
    side_gap = 0.1
  []
[]

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

[Problem]
  type = NoSolveProblem
[]

[Executioner]
  type = Transient
[]

[Outputs]
  exodus = true
[]

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

[Mesh]
    [assembly]
      type = SCMDetailedQuadAssemblyMeshGenerator
      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/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]
  [assembly]
    type = SCMDetailedQuadAssemblyMeshGenerator
    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]
    block = fuel_pins
  []
[]

[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
[]

(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]
    [assembly]
      type = SCMDetailedQuadAssemblyMeshGenerator
    []
  []

  [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/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]
  [assembly]
    type = SCMDetailedQuadAssemblyMeshGenerator
  []
[]

[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/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]
    [assembly]
        type = SCMDetailedQuadAssemblyMeshGenerator
    []
[]

[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/transient/3d.i)

[Mesh]
  [assembly]
    type = SCMDetailedQuadAssemblyMeshGenerator
    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]
    block = fuel_pins
  []
[]

[Problem]
  type = NoSolveProblem
[]

[Outputs]
  exodus = true
[]

[Executioner]
  type = Transient
[]

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

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

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

[Mesh]
  [assembly]
    type = SCMDetailedQuadAssemblyMeshGenerator
    nx = 6
    ny = 6
    n_cells = 20
    pitch = 0.0126
    pin_diameter = 0.00950
    side_gap = 0.00095
    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/examples/mesh_generator/3d.i)

[Mesh]
  [assembly]
    type = SCMDetailedQuadAssemblyMeshGenerator
    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
[]

Overview
Input Parameters
Input Files