SCMMixingKimAndChung

Class that models the turbulent mixing coefficient using the Kim and Chung correlations.

Overview

This closure class is used to model the turbulent mixing coefficient $var element = document.getElementById("moose-equation-c1dc9452-73ba-4142-8a51-cd471f1bd646");katex.render("\\beta", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ using the Kim and Chung correlations. Specifically this closure model applies to triangular and quadrilateral assemblies with bare pins. The implementation followed:

A scale analysis of the turbulent mixing rate for various Prandtl number flow fields in rod bundles eq 25,Kim and Chung (2001) Kim and Chung (2001).
Modeling of flow blockage in a liquid metal-cooled reactor subassembly with a subchannel analysis code eq 19, Jeong et. al (2005)Jeong et al. (2005).

The implemented mixing Stanton number is

$var element = document.getElementById("moose-equation-c0ad8db4-8882-4407-816d-79a1c9f4135e");katex.render("St_g = \\frac{2}{\\gamma^2}\\sqrt{\\frac{a}{8}}\\frac{D_h}{g} \\left[ \\frac{\\lambda}{Pr_t} + a_x \\frac{z_{FP}}{D} Str \\right] Re^{-b/2},", element, {displayMode:true,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$

where $var element = document.getElementById("moose-equation-3d8d0034-a3b5-48f1-b44a-9ff3e7964338");katex.render("Pr_t = Pr(Re/\\gamma)\\sqrt{f/8}", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ , $var element = document.getElementById("moose-equation-7fc29bfc-2c7c-4858-afb5-e7b33a64b6a5");katex.render("\\lambda = g/L_x", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ , and $var element = document.getElementById("moose-equation-dc924fa0-082c-40bc-a885-143ee9b9dc5d");katex.render("L_x", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ is the lattice-dependent axial length scale. Because this closure is intended for low-Prandtl-number liquid-metal applications, the implementation flags a solution warning when the local average $var element = document.getElementById("moose-equation-4e2883bb-7be9-40a3-8e10-51738089e22f");katex.render("Pr", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ is outside the expected low- $var element = document.getElementById("moose-equation-93639e8e-d2e6-423a-8c1c-187cff9a3eae");katex.render("Pr", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ range.

Information about the use of $var element = document.getElementById("moose-equation-2fa75da3-b099-4f98-88e3-437061daf644");katex.render("\\beta", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ can be found in Turbulent crossflow. Additionally, the user may opt to provide the turbulent momentum mixing parameter CT. Information about the use of this parameter can be found in Turbulent momentum transfer.

Input Parameters

CT1Turbulent momentum, modeling parameter [-].
Default:1
C++ Type:Real
Unit:(no unit assumed)
Controllable:No
Description:Turbulent momentum, modeling parameter [-].

Optional Parameters

control_tagsAdds user-defined labels for accessing object parameters via control logic.
C++ Type:std::vector<std::string>
Controllable:No
Description:Adds user-defined labels for accessing object parameters via control logic.
enableTrueSet the enabled status of the MooseObject.
Default:True
C++ Type:bool
Controllable:Yes
Description:Set the enabled status of the MooseObject.

Advanced Parameters

prop_getter_suffixAn optional suffix parameter that can be appended to any attempt to retrieve/get material properties. The suffix will be prepended with a '_' character.
C++ Type:MaterialPropertyName
Unit:(no unit assumed)
Controllable:No
Description:An optional suffix parameter that can be appended to any attempt to retrieve/get material properties. The suffix will be prepended with a '_' character.
use_interpolated_stateFalseFor the old and older state use projected material properties interpolated at the quadrature points. To set up projection use the ProjectedStatefulMaterialStorageAction.
Default:False
C++ Type:bool
Controllable:No
Description:For the old and older state use projected material properties interpolated at the quadrature points. To set up projection use the ProjectedStatefulMaterialStorageAction.

Material Property Retrieval Parameters

Input Files

(modules/subchannel/test/tests/problems/Lead-LBE-19pin/test_LEAD-19pin.i)
(modules/subchannel/test/tests/problems/SFR/sodium-7pin/case1_scm.i)
(modules/subchannel/test/tests/multiapp/sc_core.i)
(modules/subchannel/test/tests/problems/psbt/psbt_explicit_v2.i)
(modules/subchannel/test/tests/problems/Lead-LBE-19pin/test_LBE-19pin.i)

References

Hae-Yong Jeong, Kwi-Seok Ha, Won-Pyo Chang, Young-Min Kwon, and Yong-Bum Lee. Modeling of flow blockage in a liquid metal-cooled reactor subassembly with a subchannel analysis code. Nuclear technology, 149(1):71–87, 2005.

BibTeX

@article{jeong2005modeling,
    author = "Jeong, Hae-Yong and Ha, Kwi-Seok and Chang, Won-Pyo and Kwon, Young-Min and Lee, Yong-Bum",
    title = "Modeling of flow blockage in a liquid metal-cooled reactor subassembly with a subchannel analysis code",
    journal = "Nuclear technology",
    volume = "149",
    number = "1",
    pages = "71--87",
    year = "2005",
    publisher = "Taylor \\& Francis"
}

RIS

TY  - JOUR
AU  - Jeong, Hae-Yong
AU  - Ha, Kwi-Seok
AU  - Chang, Won-Pyo
AU  - Kwon, Young-Min
AU  - Lee, Yong-Bum
TI  - Modeling of flow blockage in a liquid metal-cooled reactor subassembly with a subchannel analysis code
JO  - Nuclear technology
PY  - 2005
VL  - 149
IS  - 1
PB  - Taylor & Francis
SP  - 71
EP  - 87
ER  -

Plain Text

Hae-Yong Jeong, Kwi-Seok Ha, Won-Pyo Chang, Young-Min Kwon, and Yong-Bum Lee. Modeling of flow blockage in a liquid metal-cooled reactor subassembly with a subchannel analysis code. Nuclear technology, 149(1):71–87, 2005.

Sin Kim and Bum-Jin Chung. A scale analysis of the turbulent mixing rate for various prandtl number flow fields in rod bundles. Nuclear engineering and design, 205(3):281–294, 2001.

BibTeX

@article{kim2001scale,
    author = "Kim, Sin and Chung, Bum-Jin",
    title = "A scale analysis of the turbulent mixing rate for various Prandtl number flow fields in rod bundles",
    journal = "Nuclear engineering and design",
    volume = "205",
    number = "3",
    pages = "281--294",
    year = "2001",
    publisher = "Elsevier"
}

RIS

TY  - JOUR
AU  - Kim, Sin
AU  - Chung, Bum-Jin
TI  - A scale analysis of the turbulent mixing rate for various Prandtl number flow fields in rod bundles
JO  - Nuclear engineering and design
PY  - 2001
VL  - 205
IS  - 3
PB  - Elsevier
SP  - 281
EP  - 294
ER  -

Plain Text

Sin Kim and Bum-Jin Chung. A scale analysis of the turbulent mixing rate for various prandtl number flow fields in rod bundles. Nuclear engineering and design, 205(3):281–294, 2001.

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

T_in = 673.15
flow_area = 0.00128171 #m2
rho_in = 10453.21705
# [10 m^3/hour] turns into kg/m^2-sec
mass_flux_in = '${fparse 10*rho_in/3600/flow_area}'
P_out = 1.0e5 # Pa
[TriSubChannelMesh]
  [subchannel]
    type = SCMTriAssemblyMeshGenerator
    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
    dwire = 0.0
    hwire = 0.0
    spacer_z = '0.177 0.547 0.870'
    spacer_k = '1.1719 1.1719 1.1719'
  []
[]

[FluidProperties]
  [LEAD]
    type = LeadFluidProperties
  []
[]

[SubChannel]
  type = TriSubChannel1PhaseProblem
  fp = LEAD
  n_blocks = 1
  P_out = 1.0e5
  compute_density = true
  compute_viscosity = true
  compute_power = true
  P_tol = 1.0e-4
  T_tol = 1.0e-4
  implicit = true
  segregated = false
  staggered_pressure = false
  verbose_multiapps = true
  verbose_subchannel = true
  interpolation_scheme = upwind
  pin_HTC_closure = 'Dittus-Boelter'

  # friction model
  friction_closure = 'cheng'

  full_output = true
  mixing_closure = 'Kim_and_Chung'
[]

[SCMClosures]
  [cheng]
    type = SCMFrictionUpdatedChengTodreas
  []
  [Kim_and_Chung]
    type = SCMMixingKimAndChung
  []
  [Dittus-Boelter]
    type = SCMHTCDittusBoelter
  []
[]

[ICs]


  [q_prime_IC]
    type = SCMTriPowerIC
    variable = q_prime
    power = '${fparse 250000}'
    filename = "pin_power_profile19.txt"
  []

  [Dpin_ic]
    type = ConstantIC
    variable = Dpin
    value = 8.2e-3
  []

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

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

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

  [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]
  csv = true
[]

[Postprocessors]
  [T1]
    type = SubChannelPointValue
    variable = T
    index = 37
    execute_on = "timestep_end"
    height = 0.87
  []
  [T2]
    type = SubChannelPointValue
    variable = T
    index = 36
    execute_on = "timestep_end"
    height = 0.87
  []
  [T3]
    type = SubChannelPointValue
    variable = T
    index = 20
    execute_on = "timestep_end"
    height = 0.87
  []
  [T4]
    type = SubChannelPointValue
    variable = T
    index = 10
    execute_on = "timestep_end"
    height = 0.87
  []
  [T5]
    type = SubChannelPointValue
    variable = T
    index = 4
    execute_on = "timestep_end"
    height = 0.87
  []
  [T6]
    type = SubChannelPointValue
    variable = T
    index = 1
    execute_on = "timestep_end"
    height = 0.87
  []
  [T7]
    type = SubChannelPointValue
    variable = T
    index = 14
    execute_on = "timestep_end"
    height = 0.87
  []
  [T8]
    type = SubChannelPointValue
    variable = T
    index = 28
    execute_on = "timestep_end"
    height = 0.87
  []
  ####### 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
  []
[]

[Executioner]
  type = Steady
[]

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

# [MultiApps]
#   [viz]
#     type = FullSolveMultiApp
#     input_files = "3d_LBE_19.i"
#     execute_on = "timestep_end"
#     max_procs_per_app = 1
#   []
# []

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

(modules/subchannel/test/tests/problems/SFR/sodium-7pin/case1_scm.i)

T_in = 667.15
outlet_pressure = 1.01325e5
assembly_power = 4000.0
volume_flow_rate = 0.002090242
rho_in = '${fparse 1004.23 - 0.21390 * T_in - 1.1046e-5 * T_in * T_in}'
mass_flow_rate = '${fparse volume_flow_rate * rho_in}'
core_length = 1.145
pin_pitch = 0.027813
pin_diameter = 0.024638
duct_inner_ftf = 0.0769620
n_pins = 7
flow_area = ${fparse sqrt(3)/2 * duct_inner_ftf ^2 - n_pins * pi * pin_diameter^2 / 4}
mass_flux_in = '${fparse mass_flow_rate / flow_area}'
axial_shape_scale = 1.33449883449883

[TriSubChannelMesh]
  [assembly]
    type = SCMTriAssemblyMeshGenerator
    nrings = 2
    n_cells = 20
    heated_length = ${core_length}
    flat_to_flat = ${duct_inner_ftf}
    pin_diameter = ${pin_diameter}
    pitch = ${pin_pitch}
    dwire = 0.0
    hwire = 0.0
    spacer_z = '0.0'
    spacer_k = '0.0'
  []
[]

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

[Functions]
  [axial_heat_rate]
    type = PiecewiseConstant
    axis = z
    x = '0.000 0.118 0.641 0.692 1.027'
    y = '0.0 ${axial_shape_scale} 0.0 ${axial_shape_scale} 0.0'
    direction = left
  []
[]

[SubChannel]
  type = TriSubChannel1PhaseProblem
  fp = sodium
  n_blocks = 1
  P_out = ${outlet_pressure}
  compute_density = true
  compute_viscosity = true
  compute_power = true
  implicit = true
  segregated = false
  interpolation_scheme = upwind
  gravity = counter_flow
  friction_closure = uctd_friction
  mixing_closure = kim_chung_mixing
  pin_HTC_closure = kazimi_carelli_htc
  P_tol = 1.0e-8
  T_tol = 1.0e-8
  full_output = true
  verbose_subchannel = true
[]

[SCMClosures]
  [uctd_friction]
    type = SCMFrictionUpdatedChengTodreas
  []
  [kim_chung_mixing]
    type = SCMMixingKimAndChung
  []
  [kazimi_carelli_htc]
    type = SCMHTCKazimiCarelli
  []
[]

[ICs]
  [q_prime_ic]
    type = SCMTriPowerIC
    variable = q_prime
    power = ${assembly_power}
    filename = "pin_power_profile7.txt"
    axial_heat_rate = axial_heat_rate
  []
  [Dpin_ic]
    type = ConstantIC
    variable = Dpin
    value = ${pin_diameter}
  []
  [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
  []
  [mu_ic]
    type = ViscosityIC
    variable = mu
    p = ${outlet_pressure}
    T = T
    fp = sodium
  []
  [rho_ic]
    type = RhoFromPressureTemperatureIC
    variable = rho
    p = ${outlet_pressure}
    T = T
    fp = sodium
  []
  [h_ic]
    type = SpecificEnthalpyFromPressureTemperatureIC
    variable = h
    p = ${outlet_pressure}
    T = T
    fp = sodium
  []
  [mdot_ic]
    type = ConstantIC
    variable = mdot
    value = 0.0
  []
[]

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

[Postprocessors]
  # For nrings = 2, channels 0-5 are center channels and the perimeter channels then alternate
  # edge/corner. Sample one center, one edge, and one corner channel at the outlet.
  [T_center_out]
    type = SubChannelPointValue
    variable = T
    index = 0
    execute_on = 'timestep_end'
    height = ${core_length}
  []
  [T_edge_out]
    type = SubChannelPointValue
    variable = T
    index = 6
    execute_on = 'timestep_end'
    height = ${core_length}
  []
  [T_corner_out]
    type = SubChannelPointValue
    variable = T
    index = 7
    execute_on = 'timestep_end'
    height = ${core_length}
  []
  [mdot_center_out]
    type = SubChannelPointValue
    variable = mdot
    index = 0
    execute_on = 'timestep_end'
    height = ${core_length}
  []
  [mdot_edge_out]
    type = SubChannelPointValue
    variable = mdot
    index = 6
    execute_on = 'timestep_end'
    height = ${core_length}
  []
  [mdot_corner_out]
    type = SubChannelPointValue
    variable = mdot
    index = 7
    execute_on = 'timestep_end'
    height = ${core_length}
  []
[]

[Executioner]
  type = Steady
[]

[Outputs]
  csv = true
  exodus = false
[]

(modules/subchannel/test/tests/multiapp/sc_core.i)

# Following Advanced Burner Test Reactor Preconceptual Design Report
# Vailable at: https://www.ne.anl.gov/eda/ABTR_1cv2_ws.pdf
###################################################
# Thermal-hydraulics parameters
###################################################
T_in = 866.0
P_out = 253727.1 # Pa
reactor_power = 671337.24 #WTh
mass_flow = '${fparse 6.15}' # kg/(s)

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

# units are cm - do not forget to convert to meter
scale_factor = 0.01
fuel_pin_pitch = '${fparse 1.4478*scale_factor}'
fuel_pin_diameter = '${fparse 1.4268*scale_factor}'
wire_z_spacing = '${fparse 0*scale_factor}'
wire_diameter = '${fparse 0*scale_factor}'
n_rings = 8
length_heated_fuel = '${fparse 35.56*scale_factor}'
entry_length = 0
duct_inside = '${fparse 11.43*2*scale_factor}'

###################################################

[TriSubChannelMesh]
  [subchannel]
    type = SCMTriAssemblyMeshGenerator
    nrings = '${fparse n_rings}'
    n_cells = 10
    flat_to_flat = '${fparse duct_inside}'
    heated_length = '${fparse length_heated_fuel}'
    pin_diameter = '${fparse fuel_pin_diameter}'
    pitch = '${fparse fuel_pin_pitch}'
    dwire = '${fparse wire_diameter}'
    hwire = '${fparse wire_z_spacing}'
    spacer_z = '0'
    spacer_k = '0'
  []


  [duct]
    type = SCMTriDuctMeshGenerator
    input = subchannel
    nrings = '${fparse n_rings}'
    n_cells = 10
    flat_to_flat = '${fparse duct_inside}'
    heated_length = '${fparse length_heated_fuel}'
    pitch = '${fparse fuel_pin_pitch}'
  []
[]

# All needed aux variables are automatically loaded if [SubChannel] block exists
# [AuxVariables]
#   [mdot]
#     block = subchannel
#   []
#   [SumWij]
#     block = subchannel
#   []
#   [P]
#     block = subchannel
#   []
#   [DP]
#     block = subchannel
#   []
#   [h]
#     block = subchannel
#   []
#   [T]
#     block = subchannel
#   []
#   [Tpin]
#     block = fuel_pins
#   []
#   [Dpin]
#     block = fuel_pins
#   []
#   [rho]
#     block = subchannel
#   []
#   [S]
#     block = subchannel
#   []
#   [w_perim]
#     block = subchannel
#   []
#   [q_prime]
#     block = fuel_pins
#   []
#   [mu]
#     block = subchannel
#   []
#   [q_prime_duct]
#     block = duct
#     initial_condition = 0
#   []
#   [Tduct]
#     block = duct
#   []
#   [displacement]
#     block = subchannel
#     initial_condition = 0
#   []
# []

[FluidProperties]
  [sodium]
    type = SimpleFluidProperties
    molar_mass = 0.0355
    cp = 873.0
    cv = 873.0
    specific_entropy = 1055
    viscosity = 0.0001582
    thermal_conductivity = 25.9
    thermal_expansion = 2.77e-4
  []
[]

[SubChannel]
  type = TriSubChannel1PhaseProblem
  fp = sodium
  n_blocks = 1
  P_out = ${P_out}
  P_tol = 1.0e-2
  T_tol = 1.0e-2

  # Solver settings
  implicit = true
  segregated = false

  # Output
  verbose_multiapps = true
  verbose_subchannel = true
  compute_density = false
  compute_viscosity = false
  compute_power = false

  # Heat Transfer Correlations
  pin_HTC_closure = 'gnielinski'
  duct_HTC_closure = 'gnielinski'
  # Friction Correlation
  friction_closure = 'Cheng'
  full_output = true
  # Mixing Correlation
  mixing_closure = 'Kim'

[]

[SCMClosures]
  [Cheng]
    type = SCMFrictionUpdatedChengTodreas
  []
  [gnielinski]
    type = SCMHTCGnielinski
  []
  [Kim]
    type = SCMMixingKimAndChung
  []
[]

[ICs]


  [q_prime_IC]
    type = SCMTriPowerIC
    variable = q_prime
    power = ${reactor_power} # W
    filename = 'pin_p.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
  []

  [T_duct_ic]
    type = ConstantIC
    variable = Tduct
    value = ${T_in}
  []
  [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'
    block = subchannel
  []
  [mdot_in_bc]
    type = SCMFlatMassFlowRateAux
    variable = mdot
    boundary = inlet
    mass_flow = ${mass_flow}
    execute_on = 'timestep_begin'
    block = subchannel
  []
[]

[Executioner]
  type = Steady
[]

[VectorPostprocessors]
  [sub]
    type = LineValueSampler
    start_point = '0 -0.00835888 ${entry_length}'
    end_point = '0 -0.00835888 ${fparse entry_length + length_heated_fuel}'
    num_points = 10
    variable = 'h rho P'
    sort_by = 'z'
    execute_on = 'timestep_end'
  []
[]

[Outputs]
  csv = true
[]

(modules/subchannel/test/tests/problems/psbt/psbt_explicit_v2.i)

T_in = 359.15
# [1e+6 kg/m^2-hour] turns into kg/m^2-sec
mass_flux_in = '${fparse 1e+6 * 17.00 / 3600.}'
P_out = 4.923e6 # Pa
pin_diameter = 0.00950

[QuadSubChannelMesh]
  [sub_channel]
    type = SCMQuadAssemblyMeshGenerator
    nx = 6
    ny = 6
    n_cells = 10
    pitch = 0.0126
    pin_diameter = ${pin_diameter}
    side_gap = 0.00095
    heated_length = 1.0
    spacer_z = '0.0'
    spacer_k = '0.0'
  []
[]

[FluidProperties]
  [water]
    type = Water97FluidProperties
  []
[]

[SubChannel]
  type = QuadSubChannel1PhaseProblem
  fp = water
  n_blocks = 1
  compute_density = true
  compute_viscosity = true
  compute_power = true
  P_out = report_pressure_outlet
  verbose_subchannel = true
  mixing_closure = 'Kim_and_Chung'
  friction_closure = 'Cheng'
  pin_HTC_closure = 'Dittus-Boelter'
  full_output = true
[]

[SCMClosures]
  [Cheng]
    type = SCMFrictionUpdatedChengTodreas
  []
  [Dittus-Boelter]
    type = SCMHTCDittusBoelter
  []
  [Kim_and_Chung]
    type = SCMMixingKimAndChung
    CT = 2.6
  []
[]

[ICs]


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

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

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

  [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'
    block = subchannel
  []
  [mdot_in_bc]
    type = SCMMassFlowRateAux
    variable = mdot
    boundary = inlet
    area = S
    mass_flux = report_mass_flux_inlet
    execute_on = 'timestep_begin'
    block = subchannel
  []
  [q_prime_IC]
    type = SCMQuadPowerAux
    variable = q_prime
    power = 1.0e6 # W
    filename = "power_profile.txt" #type in name of file that describes radial power profile
    execute_on = 'initial timestep_begin'
  []
[]

[Postprocessors]
  [report_mass_flux_inlet]
    type = Receiver
    default = ${mass_flux_in}
  []
  [report_pressure_outlet]
    type = Receiver
    default = ${P_out}
  []
  [total_pressure_drop]
    type = SubChannelDelta
    variable = P
    execute_on = "timestep_end"
  []
  [T1]
    type = SubChannelPointValue
    variable = T
    index = 0
    execute_on = "timestep_end"
    height = 1
  []
  [T2]
    type = SubChannelPointValue
    variable = T
    index = 7
    execute_on = "timestep_end"
    height = 1
  []
  [T3]
    type = SubChannelPointValue
    variable = T
    index = 14
    execute_on = "timestep_end"
    height = 1
  []
  [T4]
    type = SubChannelPointValue
    variable = T
    index = 21
    execute_on = "timestep_end"
    height = 1
  []
  [T5]
    type = SubChannelPointValue
    variable = T
    index = 28
    execute_on = "timestep_end"
    height = 1
  []
  [T6]
    type = SubChannelPointValue
    variable = T
    index = 35
    execute_on = "timestep_end"
    height = 1
  []
[]

[Outputs]
  csv = true
  [Temp_Out_MATRIX]
    type = QuadSubChannelNormalSliceValues
    variable = T
    execute_on = final
    file_base = "Temp_Out_Explicit.txt"
    height = 1.0
  []
[]

[Executioner]
  type = Steady
[]

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

T_in = 673.15
flow_area = 0.00128171 #m2
rho_in = 10453.21705
# [10 m^3/hour] turns into kg/m^2-sec
mass_flux_in = '${fparse 10*rho_in/3600/flow_area}'
P_out = 1.0e5 # Pa
[TriSubChannelMesh]
  [subchannel]
    type = SCMTriAssemblyMeshGenerator
    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
    dwire = 0.0
    hwire = 0.0
    spacer_z = '0.177 0.547 0.870'
    spacer_k = '1.1719 1.1719 1.1719'
  []
[]

[FluidProperties]
  [LBE]
    type = LeadBismuthFluidProperties
  []
[]

[SubChannel]
  type = TriSubChannel1PhaseProblem
  fp = LBE
  n_blocks = 1
  P_out = 1.0e5
  compute_density = true
  compute_viscosity = true
  compute_power = true
  implicit = true
  segregated = false
  verbose_subchannel = true
  interpolation_scheme = upwind
  pin_HTC_closure = 'Dittus-Boelter'

  # friction model
  friction_closure = 'cheng'

  full_output = true
  mixing_closure = 'Kim_and_Chung'
[]

[SCMClosures]
  [cheng]
    type = SCMFrictionUpdatedChengTodreas
  []
  [Kim_and_Chung]
    type = SCMMixingKimAndChung
  []
  [Dittus-Boelter]
    type = SCMHTCDittusBoelter
  []
[]

[ICs]


  [q_prime_IC]
    type = SCMTriPowerIC
    variable = q_prime
    power = '${fparse 250000}'
    filename = "pin_power_profile19.txt"
  []

  [Dpin_ic]
    type = ConstantIC
    variable = Dpin
    value = 8.2e-3
  []

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

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

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

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

[Postprocessors]
  [T1]
    type = SubChannelPointValue
    variable = T
    index = 37
    execute_on = "timestep_end"
    height = 0.87
  []
  [T2]
    type = SubChannelPointValue
    variable = T
    index = 36
    execute_on = "timestep_end"
    height = 0.87
  []
  [T3]
    type = SubChannelPointValue
    variable = T
    index = 20
    execute_on = "timestep_end"
    height = 0.87
  []
  [T4]
    type = SubChannelPointValue
    variable = T
    index = 10
    execute_on = "timestep_end"
    height = 0.87
  []
  [T5]
    type = SubChannelPointValue
    variable = T
    index = 4
    execute_on = "timestep_end"
    height = 0.87
  []
  [T6]
    type = SubChannelPointValue
    variable = T
    index = 1
    execute_on = "timestep_end"
    height = 0.87
  []
  [T7]
    type = SubChannelPointValue
    variable = T
    index = 14
    execute_on = "timestep_end"
    height = 0.87
  []
  [T8]
    type = SubChannelPointValue
    variable = T
    index = 28
    execute_on = "timestep_end"
    height = 0.87
  []
  ####### 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
  []
[]

[Outputs]
  csv = true
[]

[Executioner]
  type = Steady
[]

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

# [MultiApps]
#   [viz]
#     type = FullSolveMultiApp
#     input_files = "3d_LBE_19.i"
#     execute_on = "timestep_end"
#     max_procs_per_app = 1
#   []
# []

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

Overview
Input Parameters
Input Files
References