SCMTriAssemblyMeshGenerator

Creates a mesh of 1D subchannels and 1D pins in a triangular lattice arrangement

Overview

This mesh generator creates the 1D subchannel and pin meshes for a triangular lattice arrangement. The generated subdomains are named subchannel and fuel_pins. The center of the mesh is the origin.

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

Example Input File Syntax

[TriSubChannelMesh<<<{"href": "../../syntax/TriSubChannelMesh/index.html"}>>>]
  [subchannel]
    type = SCMTriAssemblyMeshGenerator<<<{"description": "Creates a mesh of 1D subchannels and 1D pins in a triangular lattice arrangement", "href": "SCMTriAssemblyMeshGenerator.html"}>>>
    pin_diameter<<<{"description": "Rod diameter [m]"}>>> = 0.01
    dwire<<<{"description": "Wire diameter [m]"}>>> = 0.002
    hwire<<<{"description": "Wire lead length [m]"}>>> = 0.0833
    spacer_z<<<{"description": "Axial location of spacers/vanes/mixing vanes [m]"}>>> = '0'
    spacer_k<<<{"description": "K-loss coefficient of spacers/vanes/mixing vanes [-]"}>>> = '5.0'
  []
  [duct]
    type = SCMTriDuctMeshGenerator<<<{"description": "Creates a mesh of 2D duct cells around a triangular lattice subassembly", "href": "SCMTriDuctMeshGenerator.html"}>>>
    input<<<{"description": "The corresponding subchannel mesh"}>>> = subchannel
  []
[]

Input Parameters

dwireWire diameter [m]
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Wire diameter [m]
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]
hwireWire lead length [m]
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Wire lead 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

Kij0.5Lateral form loss coefficient [-]
Default:0.5
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Lateral form loss coefficient [-]
index_blockage0 index of subchannels affected by blockage
Default:0
C++ Type:std::vector<unsigned int>
Controllable:No
Description:index of subchannels affected by blockage
k_blockage0 Form loss coefficient of subchannels affected by blockage
Default:0
C++ Type:std::vector<double>
Unit:(no unit assumed)
Controllable:No
Description:Form loss coefficient of subchannels affected by blockage
pin_block_id1Fuel Pin block id
Default:1
C++ Type:unsigned int
Controllable:No
Description:Fuel Pin block id
reduction_blockage1 Area reduction of subchannels affected by blockage (number to muliply the area)
Default:1
C++ Type:std::vector<double>
Unit:(no unit assumed)
Controllable:No
Description:Area reduction of subchannels affected by blockage (number to muliply the area)
spacer_kK-loss coefficient of spacers/vanes/mixing vanes [-]
C++ Type:std::vector<double>
Unit:(no unit assumed)
Controllable:No
Description:K-loss coefficient of spacers/vanes/mixing vanes [-]
spacer_zAxial location of spacers/vanes/mixing vanes [m]
C++ Type:std::vector<double>
Unit:(no unit assumed)
Controllable:No
Description:Axial location of spacers/vanes/mixing vanes [m]
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]
z_blockage0 0 axial location of blockage (inlet, outlet) [m]
Default:0 0
C++ Type:std::vector<double>
Unit:(no unit assumed)
Controllable:No
Description:axial location of blockage (inlet, outlet) [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

(contrib/moose/modules/subchannel/test/tests/problems/SFR/sodium-19pin/test19_monolithic.i)

T_in = 660
mass_flux_in = '${fparse 1e+6 * 300.00 / 36000.*0.5}'
P_out = 2.0e5 # Pa

[GlobalParams]
  nrings = 3
  n_cells = 5
  flat_to_flat = 0.056
  heated_length = 0.5
  pitch = 0.012
[]

[TriSubChannelMesh]
  [subchannel]
    type = SCMTriAssemblyMeshGenerator
    pin_diameter = 0.01
    dwire = 0.002
    hwire = 0.0833
    spacer_z = '0'
    spacer_k = '5.0'
  []
  [duct]
    type = SCMTriDuctMeshGenerator
    input = subchannel
  []
[]

[FluidProperties]
  [sodium]
    type = PBSodiumFluidProperties
  []
[]

[SubChannel]
  type = TriSubChannel1PhaseProblem
  fp = sodium
  n_blocks = 1
  P_out = 2.0e5
  compute_density = true
  compute_viscosity = true
  compute_power = true
  T_tol = 1.0e-6
  P_tol = 1.0e-6
  implicit = true
  segregated = false
  duct_HTC_closure = 'gnielinski'
  pin_HTC_closure = 'Dittus-Boelter'
  friction_closure = 'cheng'
  mixing_closure = 'cheng_todreas'
[]

[SCMClosures]
  [cheng]
    type = SCMFrictionUpdatedChengTodreas
  []
  [gnielinski]
    type = SCMHTCGnielinski
  []
  [cheng_todreas]
    type = SCMMixingChengTodreas
  []
  [Dittus-Boelter]
    type = SCMHTCDittusBoelter
  []
[]

[ICs]


  [q_prime_IC]
    type = SCMTriPowerIC
    variable = q_prime
    power = 1000.0 # W
    filename = "pin_power_profile19.txt"
  []

  [Dpin_ic]
    type = ConstantIC
    variable = Dpin
    value = 0.01
  []

  [T_ic]
    type = ConstantIC
    variable = T
    value = ${T_in}
  []

  [Viscosity_ic]
    type = ViscosityIC
    variable = mu
    p = ${P_out}
    T = T
    fp = sodium
  []

  [rho_ic]
    type = RhoFromPressureTemperatureIC
    variable = rho
    p = ${P_out}
    T = T
    fp = sodium
  []

  [h_ic]
    type = SpecificEnthalpyFromPressureTemperatureIC
    variable = h
    p = ${P_out}
    T = T
    fp = sodium
  []

  [mdot_ic]
    type = ConstantIC
    variable = mdot
    value = 0.0
  []
[]

[AuxKernels]
  [T_in_bc]
    type = ConstantAux
    variable = T
    boundary = inlet
    value = ${T_in}
    execute_on = 'timestep_begin'
  []
  [mdot_in_bc]
    type = SCMMassFlowRateAux
    variable = mdot
    boundary = inlet
    area = S
    mass_flux = ${mass_flux_in}
    execute_on = 'timestep_begin'
  []
[]

[Outputs]
  exodus = true
  csv = true
[]

[Postprocessors]
  [T1]
    type = SubChannelPointValue
    variable = T
    index = 37
    execute_on = "timestep_end"
    height = 0.5
  []
  [T2]
    type = SubChannelPointValue
    variable = T
    index = 36
    execute_on = "timestep_end"
    height = 0.5
  []
  [T3]
    type = SubChannelPointValue
    variable = T
    index = 20
    execute_on = "timestep_end"
    height = 0.5
  []
  [T4]
    type = SubChannelPointValue
    variable = T
    index = 10
    execute_on = "timestep_end"
    height = 0.5
  []
  [T5]
    type = SubChannelPointValue
    variable = T
    index = 4
    execute_on = "timestep_end"
    height = 0.5
  []
  [T6]
    type = SubChannelPointValue
    variable = T
    index = 1
    execute_on = "timestep_end"
    height = 0.5
  []
  [T7]
    type = SubChannelPointValue
    variable = T
    index = 14
    execute_on = "timestep_end"
    height = 0.5
  []
  [T8]
    type = SubChannelPointValue
    variable = T
    index = 28
    execute_on = "timestep_end"
    height = 0.5
  []
  ####### Assembly pressure drop
  [DP_SubchannelDelta]
    type = SubChannelDelta
    variable = P
    execute_on = 'TIMESTEP_END'
  []
  #####
  [Mean_Temp]
    type = SCMPlanarMean
    variable = T
    height = 2
  []
  [Total_power]
    type = ElementIntegralVariablePostprocessor
    variable = q_prime
    block = fuel_pins
  []
  [mdot-8]
    type = SubChannelPointValue
    variable = mdot
    index = 28
    execute_on = 'TIMESTEP_END'
    height = 0.5
  []
[]

[Executioner]
  type = Steady
[]

Overview
Example Input File Syntax
Input Parameters