SCMTriInterWrapperMeshGenerator

Creates a mesh for the inter-wrapper around triangular subassemblies

Overview

This mesh generator creates the mesh were the SCM solution variables associated with the inter-wrapper live. The inter-wrapper is the flow area that wraps arround the subchannel sub-assemblies. The user needs to provide the "nrings" parameter which defines the number of sub-assemblies the inter-wrapper wraps around. For triangular sub-assemblies one nring would define one center sub-assembly and six neighboring ones. The number of cells in the -z direction is given by "n_cells". The distance of the sub-assemblies from eachother is given by the "assembly_pitch" parameter and the total length of the inter-wrapper in the -z direction is defined by the parameters: "heated_length","unheated_length_entry","unheated_length_entry". The size of each subchannel sub-assembly including the duct is defined by "flat_to_flat". Last, "side_bypass" is the extra size of the bypass for the side sub-assemblies. The center of the mesh is the origin.

Example Input File Syntax

[TriInterWrapperMesh<<<{"href": "../../syntax/TriInterWrapperMesh/index.html"}>>>]
  [sub_channel]
    type = SCMTriInterWrapperMeshGenerator<<<{"description": "Creates a mesh for the inter-wrapper around triangular subassemblies", "href": "SCMTriInterWrapperMeshGenerator.html"}>>>
    nrings<<<{"description": "Number of fuel Pin rings per assembly [-]"}>>> = 3
    n_cells<<<{"description": "The number of cells in the axial direction"}>>> = 10
    flat_to_flat<<<{"description": "Flat to flat distance for the hexagonal assembly [m]"}>>> = 0.16
    heated_length<<<{"description": "Heated length [m]"}>>> = 1.0
    assembly_pitch<<<{"description": "Pitch [m]"}>>> = 0.17
    side_bypass<<<{"description": "Extra size of the bypass for the side assemblies [m]"}>>> = 0.01
  []
[]

Input Parameters

assembly_pitchPitch [m]
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Pitch [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]
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 [-]
C++ Type:unsigned int
Controllable:No
Description:Number of fuel Pin rings per assembly [-]
side_bypassExtra size of the bypass for the side assemblies [m]
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Extra size of the bypass for the side assemblies [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 [-]
tight_side_bypassFalseWhether the side bypass shape follows the hexagonal shape of the assemblies
Default:False
C++ Type:bool
Controllable:No
Description:Whether the side bypass shape follows the hexagonal shape of the assemblies
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

(moose/modules/subchannel/test/tests/problems/interwrapper/tri_interwrapper_monolithic.i)

T_in = 360.0
P_out = 4.923e6 # Pa

[TriInterWrapperMesh]
  [sub_channel]
    type = SCMTriInterWrapperMeshGenerator
    nrings = 3
    n_cells = 10
    flat_to_flat = 0.16
    heated_length = 1.0
    assembly_pitch = 0.17
    side_bypass = 0.01
  []
[]

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

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

[ICs]
  [S_IC]
    type = TriInterWrapperFlowAreaIC
    variable = S
  []

  [w_perim_IC]
    type = TriInterWrapperWettedPerimIC
    variable = w_perim
  []

  [q_prime_IC]
    type = TriInterWrapperPowerIC
    variable = q_prime
  []

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

  [P_ic]
    type = ConstantIC
    variable = P
    value = 0.0
  []

  [DP_ic]
    type = ConstantIC
    variable = DP
    value = 0.0
  []

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

[SubChannel]
  type = TriInterWrapper1PhaseProblem
  fp = sodium
  n_blocks = 1
  beta = 0.1
  P_out = 2.0e5
  CT = 1.0
  compute_density = false
  compute_viscosity = false
  compute_power = false
  P_tol = 1.0e-6
  T_tol = 1.0e-6
  implicit = true
  segregated = false
  staggered_pressure = false
  monolithic_thermal = false
  interpolation_scheme = 'central_difference'
[]

[AuxKernels]
[]

[Outputs]
  exodus = true
  checkpoint = false
[]

[Executioner]
  type = Steady
[]

Overview
Example Input File Syntax
Input Parameters