SCMDetailedTriAssemblyMeshGenerator

Creates a detailed mesh of subchannels and pins in a triangular lattice arrangement

Overview

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

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

Input Parameters

flat_to_flatFlat to flat distance for the hexagonal assembly [m]
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Flat to flat distance for the hexagonal assembly [m]
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
nringsNumber of fuel-pin rings per assembly, counting the center pin as the first ring [-]
C++ Type:unsigned int
Controllable:No
Description:Number of fuel-pin rings per assembly, counting the center pin as the first ring [-]
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]

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]
verbose_flagFalseFlag to print out the mesh coordinates.
Default:False
C++ Type:bool
Controllable:No
Description:Flag to print out the mesh coordinates.

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/validation/areva_FCTF/3D.i)
(modules/combined/test/tests/subchannel_thm_coupling/3D.i)
(modules/subchannel/validation/EBR-II/3d_SCM_SS.i)
(modules/subchannel/validation/Blockage/THORS/FFM-5B_viz.i)
(modules/subchannel/examples/MultiApp/3d.i)
(modules/subchannel/test/tests/problems/Lead-LBE-19pin/3D_LBE-19pin.i)
(modules/subchannel/validation/Toshiba_37_pin/toshiba_37_pin_viz.i)
(modules/subchannel/validation/ORNL_19_pin/3d_ORNL_19.i)
(modules/subchannel/validation/Blockage/THORS/FFM-3A_viz.i)
(modules/subchannel/validation/EBR-II/3d_SCM_TR.i)
(modules/subchannel/examples/duct/3d.i)
(modules/subchannel/test/tests/mesh/detailed_tri_assembly/coords.i)

(modules/subchannel/validation/areva_FCTF/3D.i)

###################################################
# Geometric parameters
#####################################
fuel_pin_pitch = 0.01125 #m
fuel_pin_diameter = 0.0095 #m
inner_duct_in = 0.09164 #m
n_rings = 5
unheated_length_entry = 1.14 #m
heated_length = 1.71 #m
unheated_length_exit = 0.855 #m
###################################################

[TriSubChannelMesh]
  [subchannel]
    type = SCMDetailedTriAssemblyMeshGenerator
    nrings = '${n_rings}'
    n_cells = 65
    flat_to_flat = '${inner_duct_in}'
    unheated_length_entry = ${unheated_length_entry}
    heated_length = ${heated_length}
    unheated_length_exit = ${unheated_length_exit}
    pin_diameter = '${fuel_pin_diameter}'
    pitch = '${fuel_pin_pitch}'
  []

[]

[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
  []
  [displacement]
    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/combined/test/tests/subchannel_thm_coupling/3D.i)

###################################################
# Geometric parameters
###################################################
n_cells = 25
n_rings = 3
fuel_pin_pitch = 7.26e-3
fuel_pin_diameter =  5.84e-3
inner_duct_in = 3.41e-2
heated_length = 1.0
###################################################
[TriSubChannelMesh]
  [subchannel]
    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}
  []
[]

[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
  []
  [Sij]
    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/EBR-II/3d_SCM_SS.i)

# Following Advanced Burner Test Reactor Preconceptual Design Report
# Vailable at: https://www.ne.anl.gov/eda/ABTR_1cv2_ws.pdf
# This file creates the 3D mesh where subchannel projects to
###################################################
# Geometric parameters
##################################################
scale_factor = 0.01
fuel_pin_pitch = '${fparse 0.5664*scale_factor}'
fuel_pin_diameter = '${fparse 0.4419*scale_factor}'
inner_duct_in = '${fparse 4.64*scale_factor}'
n_rings = 5
heated_length = '${fparse 34.3*scale_factor}'
unheated_length_exit = '${fparse 26.9*scale_factor}'
###################################################

[TriSubChannelMesh]
  [subchannel]
    type = SCMDetailedTriAssemblyMeshGenerator
    nrings = '${n_rings}'
    n_cells = 50
    flat_to_flat = '${inner_duct_in}'
    unheated_length_exit = '${unheated_length_exit}'
    heated_length = '${heated_length}'
    pin_diameter = '${fuel_pin_diameter}'
    pitch = '${fuel_pin_pitch}'
  []

[]

[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
  []
  [displacement]
    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/THORS/FFM-5B_viz.i)

################################################################################
## THORS bundle 5B partial edge blockage benchmark                            ##
## SCM simulation                                                             ##
## POC : Vasileios Kyriakopoulos, vasileios.kyriakopoulos@inl.gov             ##
################################################################################
# Details on the experimental facility modeled can be found at:
# Han, J. T. "Blockages in LMFBR fuel assemblies: A review of experimental and theoretical studies." (1977).
# This input file creates a 3D subchannel mesh that gets the subchannel information projected on,
# for visualization purposes

[Mesh]
  [subchannel]
    type = SCMDetailedTriAssemblyMeshGenerator
    nrings = 3
    n_cells = 50
    flat_to_flat = 3.41e-2
    heated_length = 0.4572
    unheated_length_entry = 0.4064
    unheated_length_exit = 0.1524
    pin_diameter = 5.84e-3
    pitch = 7.26e-3
  []
[]

[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
  []
  [displacement]
    block = subchannel
  []
  [q_prime]
    block = fuel_pins
  []
[]

[Problem]
  type = NoSolveProblem
[]

[Outputs]
  exodus = true
[]

[Executioner]
  type = Steady
[]

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

# Following Advanced Burner Test Reactor Preconceptual Design Report
# Available at: https://www.ne.anl.gov/eda/ABTR_1cv2_ws.pdf

###################################################
# Geometric parameters
###################################################

# units are cm - do not forget to convert to meter
scale_factor = 0.01
fuel_element_pitch = ${fparse 14.598*scale_factor}
inter_assembly_gap = ${fparse 0.4*scale_factor}
duct_thickness = ${fparse 0.3*scale_factor}
fuel_pin_pitch = ${fparse 0.904*scale_factor}
fuel_pin_diameter= ${fparse 0.8*scale_factor}
# wire_z_spacing = ${fparse 20.32*scale_factor}
# wire_diameter = ${fparse 0.103*scale_factor}

# Reduced height for convenience
length_entry_fuel = '${fparse 20*scale_factor}'
length_heated_fuel = '${fparse 40*scale_factor}'
length_outlet_fuel = '${fparse 20*scale_factor}'
height = ${fparse length_entry_fuel+length_heated_fuel+length_outlet_fuel}
# orifice_plate_height = ${fparse 5*scale_factor}
duct_outside = ${fparse fuel_element_pitch - inter_assembly_gap}
duct_inside = ${fparse duct_outside - 2 * duct_thickness}
n_rings = 9

n_cells = 100

[Mesh]
  [subchannel]
    type = SCMDetailedTriAssemblyMeshGenerator
    nrings = '${fparse n_rings}'
    n_cells = ${n_cells}
    flat_to_flat = '${fparse duct_inside}'
    unheated_length_entry = '${fparse length_entry_fuel}'
    heated_length = '${fparse length_heated_fuel}'
    unheated_length_exit = '${fparse length_outlet_fuel}'
    pin_diameter = '${fparse fuel_pin_diameter}'
    pitch = '${fparse fuel_pin_pitch}'
  []

[]

[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
  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_center]
    type = LineValueSampler
    variable = 'mdot SumWij P DP h T rho mu S w_perim'
    sort_by = 'z'
    # n_nodes is n_elems + 1
    num_points = ${fparse n_cells + 1}
    start_point = '0 ${fuel_pin_pitch} 0'
    end_point = '0 ${fuel_pin_pitch} ${height}'
  []
  [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}
    # the offset places the sample inside the channel
    start_point = '${fparse fuel_pin_pitch / 2} ${fparse (fuel_element_pitch - inter_assembly_gap - 2 * duct_thickness) / 2 - 1e-3} 0'
    end_point = '${fparse fuel_pin_pitch / 2} ${fparse (fuel_element_pitch - inter_assembly_gap - 2 * duct_thickness) / 2 - 1e-3} ${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 (n_rings - 0.5) * fuel_pin_pitch} 0 0'
    end_point = '${fparse (n_rings - 0.5) * fuel_pin_pitch} 0 ${height}'
  []
  [sample_pin_center]
    type = LineValueSampler
    variable = 'q_prime Tpin Dpin'
    sort_by = 'z'
    num_points = ${fparse n_cells + 1}
    start_point = '0 0 0'
    end_point = '0 0 ${height}'
  []
[]

[Executioner]
  type = Transient
[]

(modules/subchannel/test/tests/problems/Lead-LBE-19pin/3D_LBE-19pin.i)

[Mesh]
  [subchannel]
    type = SCMDetailedTriAssemblyMeshGenerator
    nrings = 3
    n_cells = 50
    flat_to_flat = 0.05319936
    heated_length = 0.87
    unheated_length_entry = 0.0
    unheated_length_exit = 0.402
    pin_diameter = 8.2e-3
    pitch = 0.01148
  []
[]

[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/Toshiba_37_pin/toshiba_37_pin_viz.i)

[Mesh]
  [subchannel]
    type = SCMDetailedTriAssemblyMeshGenerator
    nrings = 4
    n_cells = 20
    flat_to_flat = 0.085
    heated_length = 1.0
    pin_diameter = 0.01
    pitch = 0.012
  []
[]

[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/ORNL_19_pin/3d_ORNL_19.i)

[Mesh]
  [subchannel]
    type = SCMDetailedTriAssemblyMeshGenerator
    nrings = 3
    n_cells = 40
    flat_to_flat = 3.41e-2
    heated_length = 0.5334
    unheated_length_entry = 0.4064
    unheated_length_exit = 0.0762
    pin_diameter = 5.84e-3
    pitch = 7.26e-3
  []
[]

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

(modules/subchannel/validation/Blockage/THORS/FFM-3A_viz.i)

################################################################################
## THORS bundle 3A central blockage benchmark                                 ##
## SCM simulation                                                             ##
## POC : Vasileios Kyriakopoulos, vasileios.kyriakopoulos@inl.gov             ##
################################################################################
# Details on the experimental facility modeled can be found at:
# Han, J. T. "Blockages in LMFBR fuel assemblies: A review of experimental and theoretical studies." (1977).
# This input file creates a 3D subchannel mesh that gets the subchannel information projected on,
# for visualization purposes

[Mesh]
  [subchannel]
    type = SCMDetailedTriAssemblyMeshGenerator
    nrings = 3
    n_cells = 36
    flat_to_flat = 3.41e-2
    heated_length = 0.5334
    unheated_length_entry = 0.3048
    unheated_length_exit = 0.0762
    pin_diameter = 5.84e-3
    pitch = 7.26e-3
  []
[]

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

(modules/subchannel/validation/EBR-II/3d_SCM_TR.i)

# Following Advanced Burner Test Reactor Preconceptual Design Report
# Vailable at: https://www.ne.anl.gov/eda/ABTR_1cv2_ws.
# This file creates the 3D mesh where subchannel projects to
###################################################
# Geometric parameters
###################################################
# units are cm - do not forget to convert to meter
scale_factor = 0.01
fuel_pin_pitch = '${fparse 0.5664*scale_factor}'
fuel_pin_diameter = '${fparse 0.4419*scale_factor}'
inner_duct_in = '${fparse 4.64*scale_factor}'
n_rings = 5
heated_length = '${fparse 34.3*scale_factor}'
unheated_length_exit = '${fparse 26.9*scale_factor}'
###################################################

[TriSubChannelMesh]
  [subchannel]
    type = SCMDetailedTriAssemblyMeshGenerator
    nrings = '${n_rings}'
    n_cells = 50
    flat_to_flat = '${inner_duct_in}'
    unheated_length_exit = '${unheated_length_exit}'
    heated_length = '${heated_length}'
    pin_diameter = '${fuel_pin_diameter}'
    pitch = '${fuel_pin_pitch}'
  []

[]

[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 = Transient
  start_time = -1.0
  dt = 100.0
  end_time = 900.0
[]

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

[Mesh]
  [subchannel]
    type = SCMDetailedTriAssemblyMeshGenerator
    nrings = 4
    n_cells = 100
    flat_to_flat = 0.085
    heated_length = 1.0
    pin_diameter = 0.01
    pitch = 0.012
  []
[]

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

[Problem]
  type = NoSolveProblem
[]

[Outputs]
  exodus = true
[]

[Executioner]
  type = Steady
[]

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

[GlobalParams]
  nrings = 3
  n_cells = 20
  flat_to_flat = 0.056
  heated_length = 0.2
  pitch = 0.012
[]

[TriSubChannelMesh]
  [subchannel]
    type = SCMDetailedTriAssemblyMeshGenerator
    pin_diameter = 0.01
    verbose_flag = true
  []
[]

Overview
Input Parameters
Input Files