TriInterWrapperPowerIC

Computes linear power rate (W/m) that goes into interwrapper cells in a triangular subchannel lattice

Overview

This IC creates an initial field for the linear heat rate ("variable" = q_prime), that gets added or removed to the inter-wrapper flow area, that wraps around sub-assemblies with a triangular subchannel/pin arrangement.

Input Parameters

variableThe variable this initial condition is supposed to provide values for.
C++ Type:VariableName
Unit:(no unit assumed)
Controllable:No
Description:The variable this initial condition is supposed to provide values for.

Required Parameters

axial_heat_rate1.0user provided normalized function of axial heat rate [Unitless]. The integral over pin length should equal the heated length
Default:1.0
C++ Type:FunctionName
Unit:(no unit assumed)
Controllable:No
Description:user provided normalized function of axial heat rate [Unitless]. The integral over pin length should equal the heated length
blockThe list of blocks (ids or names) that this object will be applied
C++ Type:std::vector<SubdomainName>
Controllable:No
Description:The list of blocks (ids or names) that this object will be applied
boundaryThe list of boundaries (ids or names) from the mesh where this object applies
C++ Type:std::vector<BoundaryName>
Controllable:No
Description:The list of boundaries (ids or names) from the mesh where this object applies
filenamefile_was_not_foundname of power profile .txt file (should be a single column). It's a Radial Power Profile. [UnitLess]
Default:file_was_not_found
C++ Type:std::string
Controllable:No
Description:name of power profile .txt file (should be a single column). It's a Radial Power Profile. [UnitLess]
power0The total heating power [W]
Default:0
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:The total heating power [W]
stateCURRENTThis parameter is used to set old state solutions at the start of simulation. If specifying multiple states at the start of simulation, use one IC object for each state being specified. The states are CURRENT=0 OLD=1 OLDER=2. States older than 2 are not currently supported. When the user only specifies current state, the solution is copied to the old and older states, as expected. This functionality is mainly used for dynamic simulations with explicit time integration schemes, where old solution states are used in the velocity and acceleration approximations.
Default:CURRENT
C++ Type:MooseEnum
Options:CURRENT, OLD, OLDER
Controllable:No
Description:This parameter is used to set old state solutions at the start of simulation. If specifying multiple states at the start of simulation, use one IC object for each state being specified. The states are CURRENT=0 OLD=1 OLDER=2. States older than 2 are not currently supported. When the user only specifies current state, the solution is copied to the old and older states, as expected. This functionality is mainly used for dynamic simulations with explicit time integration schemes, where old solution states are used in the velocity and acceleration approximations.

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:No
Description:Set the enabled status of the MooseObject.
ignore_uo_dependencyFalseWhen set to true, a UserObject retrieved by this IC will not be executed before the this IC
Default:False
C++ Type:bool
Controllable:No
Description:When set to true, a UserObject retrieved by this IC will not be executed before the this IC

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/interwrapper/tri_interwrapper_implicit.i)
(modules/subchannel/test/tests/problems/interwrapper/tri_interwrapper_monolithic_staggered.i)
(modules/subchannel/examples/mesh_generator/inter_wrapper_tri_mesh_generator_test.i)
(modules/subchannel/test/tests/problems/interwrapper/tri_interwrapper.i)
(modules/subchannel/test/tests/problems/interwrapper/tri_interwrapper_monolithic.i)

(modules/subchannel/test/tests/problems/interwrapper/tri_interwrapper_implicit.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 = 2
  beta = 0.1
  P_out = 2.0e5
  CT = 1.0
  compute_density = true
  compute_viscosity = true
  compute_power = false
  P_tol = 1.0e-6
  T_tol = 1.0e-6
  implicit = true
  segregated = true
  staggered_pressure = false
  monolithic_thermal = false
  interpolation_scheme = 'upwind'
[]

[AuxKernels]
[]

[Outputs]
  exodus = true
  checkpoint = false
[]

[Executioner]
  type = Steady
[]

(modules/subchannel/test/tests/problems/interwrapper/tri_interwrapper_monolithic_staggered.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 = true
  monolithic_thermal = false
  interpolation_scheme = 'central_difference'
[]

[AuxKernels]
[]

[Outputs]
  exodus = true
  checkpoint = false
[]

[Executioner]
  type = Steady
[]

(modules/subchannel/examples/mesh_generator/inter_wrapper_tri_mesh_generator_test.i)

T_in = 360.0
P_out = 4.923e6 # Pa

[TriInterWrapperMesh]
  [sub_channel]
    type = SCMTriInterWrapperMeshGenerator
    nrings = 4
    n_cells = 50
    flat_to_flat = 0.16
    heated_length = 3.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 = true
  compute_viscosity = true
  compute_power = false
  #  P_tol = 1.0e-5
  #  T_tol = 1.0e-5
  implicit = false
  segregated = true
  staggered_pressure = false
  monolithic_thermal = false
[]

[AuxKernels]
[]

[Outputs]
  exodus = true
  checkpoint = false
[]

[Executioner]
  type = Steady
[]

(modules/subchannel/test/tests/problems/interwrapper/tri_interwrapper.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 = true
  compute_viscosity = true
  compute_power = false
  P_tol = 1.0e-6
  T_tol = 1.0e-6
  implicit = false
  segregated = true
  staggered_pressure = false
  monolithic_thermal = false
  interpolation_scheme = 'upwind'
[]

[AuxKernels]
[]

[Outputs]
  exodus = true
  checkpoint = false
[]

[Executioner]
  type = Steady
[]

(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
Input Parameters
Input Files