SCMTriDuctMeshGenerator

Creates a mesh of 2D duct cells around a triangular lattice subassembly

Overview

This mesh generator creates the mesh (with Quad4 elements) where the variables associated with the duct live.

The subdomain name used for the duct will match the name of this mesh generator.

Example Input File Syntax

[TriSubChannelMesh<<<{"href": "../../syntax/TriSubChannelMesh/index.html"}>>>]
  [subchannel]
    type = SCMTriSubChannelMeshGenerator<<<{"description": "Creates a mesh of 1D subchannels in a triangular lattice arrangement", "href": "SCMTriSubChannelMeshGenerator.html"}>>>
    nrings<<<{"description": "Number of fuel Pin rings per assembly [-]"}>>> = 4
    n_cells<<<{"description": "The number of cells in the axial direction"}>>> = 100
    flat_to_flat<<<{"description": "Flat to flat distance for the hexagonal assembly [m]"}>>> = 0.085
    heated_length<<<{"description": "Heated length [m]"}>>> = 1.0
    pin_diameter<<<{"description": "Rod diameter [m]"}>>> = 0.01
    pitch<<<{"description": "Pitch [m]"}>>> = 0.012
    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 0.2 0.4 0.6 0.8'
    spacer_k<<<{"description": "K-loss coefficient of spacers/vanes/mixing_vanes [-]"}>>> = '0.1 0.1 0.1 0.1 0.10'
  []

  [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
    nrings<<<{"description": "Number of fuel Pin rings in the assembly [-]"}>>> = 4
    n_cells<<<{"description": "The number of cells in the axial direction"}>>> = 100
    flat_to_flat<<<{"description": "Flat to flat distance for the hexagonal assembly [m]"}>>> = 0.085
    heated_length<<<{"description": "Heated length [m]"}>>> = 1.0
    pitch<<<{"description": "Pitch is the distance between adjacent pins [m]"}>>> = 0.012
  []
[]

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]
inputThe corresponding subchannel mesh
C++ Type:MeshGeneratorName
Controllable:No
Description:The corresponding subchannel mesh
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 in the assembly [-]
C++ Type:unsigned int
Controllable:No
Description:Number of fuel Pin rings in the assembly [-]
pitchPitch is the distance between adjacent pins [m]
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Pitch is the distance between adjacent pins [m]

Required Parameters

block_id2Domain Index
Default:2
C++ Type:unsigned int
Controllable:No
Description:Domain Index
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

(contrib/moose/modules/subchannel/examples/duct/test.i)

T_in = 660
mass_flux_in = '${fparse 1e+6 * 37.00 / 36000.*0.5}'
P_out = 2.0e5 # Pa
[TriSubChannelMesh]
  [subchannel]
    type = SCMTriSubChannelMeshGenerator
    nrings = 4
    n_cells = 100
    flat_to_flat = 0.085
    heated_length = 1.0
    pin_diameter = 0.01
    pitch = 0.012
    dwire = 0.002
    hwire = 0.0833
    spacer_z = '0 0.2 0.4 0.6 0.8'
    spacer_k = '0.1 0.1 0.1 0.1 0.10'
  []

  [duct]
    type = SCMTriDuctMeshGenerator
    input = subchannel
    nrings = 4
    n_cells = 100
    flat_to_flat = 0.085
    heated_length = 1.0
    pitch = 0.012
  []
[]

[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
  []
  [q_prime]
    block = subchannel
  []
  [mu]
    block = subchannel
  []
  [displacement]
    block = subchannel
  []
  [duct_heat_flux]
    block = duct
  []
  [Tduct]
    block = duct
  []
[]

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

[Problem]
  type = TriSubChannel1PhaseProblem
  fp = sodium
  n_blocks = 1
  P_out = 2.0e5
  CT = 1.0
  compute_density = false
  compute_viscosity = false
  compute_power = true
  P_tol = 1.0e-5
  T_tol = 1.0e-5
  implicit = true
  segregated = false
  staggered_pressure = false
  verbose_multiapps = true
  verbose_subchannel = false
[]

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

  [w_perim_IC]
    type = SCMTriWettedPerimIC
    variable = w_perim
  []

  [q_prime_IC]
    type = SCMTriPowerIC
    variable = q_prime
    power = 1e5 #1.000e5 # W
    filename = "pin_power_profile37.txt"
  []

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

  [T_duct_ic]
    type = ConstantIC
    variable = Tduct
    value = ${T_in}
  []
[]

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

[UserObjects]
  [Tduct_avg_uo]
    type = NearestPointLayeredAverage
    direction = z
    num_layers = 1000
    variable = Tduct
    block = duct
    points = '0 0 0'
    execute_on = 'TIMESTEP_END'
  []
[]

[Outputs]
  exodus = true
[]

[Executioner]
  type = Steady
  petsc_options_iname = '-pc_type -pc_hypre_type'
  petsc_options_value = 'hypre boomeramg'
  fixed_point_max_its = 2
  fixed_point_min_its = 2
  fixed_point_rel_tol = 1e-6
[]

################################################################################
# A multiapp that projects data to a detailed mesh
################################################################################

[MultiApps]
  # Multiapp to duct heat conduction module
  [duct_map]
    type = FullSolveMultiApp
    input_files = wrapper.i # seperate file for multiapps due to radial power profile
    execute_on = 'timestep_end'
    positions = '0   0   0' #center of assembly
    bounding_box_padding = '10.0 10.0 10.0'
  []

  # Multiapp to detailed mesh for vizualization
  [viz]
    type = FullSolveMultiApp
    input_files = "3d.i"
    execute_on = 'timestep_end'
  []
[]

[Transfers]
  [duct_temperature_transfer] # Send duct temperature to heat conduction
    type = MultiAppInterpolationTransfer
    to_multi_app = duct_map
    source_variable = Tduct
    variable = duct_surface_temperature
  []

  [displacement_transfer]
    type = MultiAppGeneralFieldNearestNodeTransfer
    from_multi_app = duct_map
    source_variable = disp_magnitude
    variable = displacement
  []

  [q_prime] # Recover q_prime from heat conduction solve
    type = MultiAppInterpolationTransfer
    from_multi_app = duct_map
    source_variable = q_prime
    variable = duct_heat_flux
  []

  [xfer]
    type = SCMSolutionTransfer
    to_multi_app = viz
    variable = 'mdot SumWij P DP h T rho mu q_prime S displacement'
  []
[]

Overview
Example Input File Syntax
Input Parameters