FlowChannelGasMix

This component implements the gas mixture flow model on a flow channel, as FlowChannel1Phase does for the single-phase flow model.

Usage

The usage for this component is the same as for FlowChannel1Phase but with the addition of the initial condition parameter "initial_mass_fraction" for the mass fraction of the secondary gas.

Input Parameters

AArea of the flow channel, can be a constant or a function
C++ Type:FunctionName
Unit:(no unit assumed)
Controllable:Yes
Description:Area of the flow channel, can be a constant or a function
fpFluid properties user object
C++ Type:UserObjectName
Controllable:No
Description:Fluid properties user object
lengthLength of each axial section [m]
C++ Type:std::vector<double>
Unit:(no unit assumed)
Controllable:No
Description:Length of each axial section [m]
n_elemsNumber of elements in each axial section
C++ Type:std::vector<unsigned int>
Controllable:No
Description:Number of elements in each axial section
orientationDirection of flow channel from start position to end position (no need to normalize). For curved flow channels, it is the (tangent) direction at the start position.
C++ Type:libMesh::VectorValue<double>
Unit:(no unit assumed)
Controllable:No
Description:Direction of flow channel from start position to end position (no need to normalize). For curved flow channels, it is the (tangent) direction at the start position.
positionStart position of axis in 3-D space [m]
C++ Type:libMesh::Point
Controllable:No
Description:Start position of axis in 3-D space [m]

Required Parameters

D_hHydraulic diameter [m]
C++ Type:FunctionName
Unit:(no unit assumed)
Controllable:Yes
Description:Hydraulic diameter [m]
PoD1Pitch-to-diameter ratio for parallel bundle heat transfer [-]
Default:1
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Pitch-to-diameter ratio for parallel bundle heat transfer [-]
axial_region_namesNames to assign to axial regions
C++ Type:std::vector<std::string>
Controllable:No
Description:Names to assign to axial regions
closuresClosures object(s). This is optional since closure relations can be supplied directly by Materials as well.
C++ Type:std::vector<std::string>
Controllable:No
Description:Closures object(s). This is optional since closure relations can be supplied directly by Materials as well.
fWall friction factor [-]
C++ Type:FunctionName
Unit:(no unit assumed)
Controllable:Yes
Description:Wall friction factor [-]
gravity_vector0 0 -9.81Gravitational acceleration vector [m/s^2]
Default:0 0 -9.81
C++ Type:libMesh::VectorValue<double>
Unit:(no unit assumed)
Controllable:No
Description:Gravitational acceleration vector [m/s^2]
heat_transfer_geomPIPEConvective heat transfer geometry
Default:PIPE
C++ Type:MooseEnum
Options:HEX_ROD_BUNDLE, PIPE, ROD_BUNDLE
Controllable:No
Description:Convective heat transfer geometry
name_multiple_ht_by_indexTrueIf true, when there are multiple heat transfer components connected to this flow channel, use their index for naming related quantities; otherwise, use the name of the heat transfer component.
Default:True
C++ Type:bool
Controllable:No
Description:If true, when there are multiple heat transfer components connected to this flow channel, use their index for naming related quantities; otherwise, use the name of the heat transfer component.
pipe_locationINTERIORPipe location within the bundle
Default:INTERIOR
C++ Type:MooseEnum
Options:CORNER, EDGE, INTERIOR
Controllable:No
Description:Pipe location within the bundle
pipe_pars_transferredFalseSet to true if Dh, P_hf and A are going to be transferred in from an external source
Default:False
C++ Type:bool
Controllable:No
Description:Set to true if Dh, P_hf and A are going to be transferred in from an external source
rotation0Angle of rotation about the x-axis [degrees]
Default:0
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Angle of rotation about the x-axis [degrees]
roughness0Roughness [m]
Default:0
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Roughness [m]

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.

Advanced Parameters

initial_TInitial temperature in the flow channel [K]
C++ Type:FunctionName
Unit:(no unit assumed)
Controllable:Yes
Description:Initial temperature in the flow channel [K]
initial_mass_fractionInitial mass fraction of the secondary gas
C++ Type:FunctionName
Unit:(no unit assumed)
Controllable:No
Description:Initial mass fraction of the secondary gas
initial_pInitial pressure in the flow channel [Pa]
C++ Type:FunctionName
Unit:(no unit assumed)
Controllable:Yes
Description:Initial pressure in the flow channel [Pa]
initial_velInitial velocity in the flow channel [m/s]
C++ Type:FunctionName
Unit:(no unit assumed)
Controllable:Yes
Description:Initial velocity in the flow channel [m/s]

Variable Initialization Parameters

rdg_slope_reconstructionNONESlope reconstruction type for rDG spatial discretization
Default:NONE
C++ Type:MooseEnum
Options:FULL, MC, MINMOD, NONE, SUPERBEE
Controllable:No
Description:Slope reconstruction type for rDG spatial discretization
scaling_factor_rhoA1Scaling factor for the mixture mass equation
Default:1
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Scaling factor for the mixture mass equation
scaling_factor_rhoEA1Scaling factor for the energy equation
Default:1
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Scaling factor for the energy equation
scaling_factor_rhouA1Scaling factor for the momentum equation
Default:1
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Scaling factor for the momentum equation
scaling_factor_xirhoA1Scaling factor for the secondary component mass equation
Default:1
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Scaling factor for the secondary component mass equation

Numerical Scheme Parameters

Mesh

Mesh-related parameters of this component match FlowChannel1Phase.

Variables

Compared to FlowChannel1Phase, this component adds the following solution variable:

Variable	Symbol	Description
`xirhoA`	$var element = document.getElementById("moose-equation-eabe7b9c-f274-4364-a680-2fe2b97f261c");katex.render("\\xi \\rho A", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$	Secondary gas mass per unit length [kg/m]

and the following auxiliary variable:

Variable	Symbol	Description
`mass_fraction`	$var element = document.getElementById("moose-equation-cf0efcc7-ee9a-40b9-a622-0791603c7b7b");katex.render("\\xi", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$	Mass fraction of the secondary gas

Material Properties

Compared to FlowChannel1Phase, this component adds the material property:

Material Property	Symbol	Description
`mass_fraction`	$var element = document.getElementById("moose-equation-b45f3240-cfc4-407a-bfb6-1f8d0707d727");katex.render("\\xi", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$	Mass fraction of the secondary gas

Input Files

(modules/thermal_hydraulics/test/tests/components/flow_channel_gasmix/flow_channel_gasmix.i)

References

No citations exist within this document.

(modules/thermal_hydraulics/test/tests/components/flow_channel_gasmix/flow_channel_gasmix.i)

initial_p = 1e5
initial_T = 500

[FluidProperties]
  [fp1]
    type = IdealGasFluidProperties
    gamma = 1.4
    molar_mass = 0.029
  []
  [fp2]
    type = IdealGasFluidProperties
    gamma = 1.5
    molar_mass = 0.04
  []
  [fp_mix]
    type = IdealGasMixtureFluidProperties
    component_fluid_properties = 'fp1 fp2'
  []
[]

[Closures]
  [closures]
    type = FunctorClosures
    properties = 'f_D mass_diffusion_coefficient'
    functors = '0 0.26e-4'
  []
[]

[Functions]
  [initial_mass_fraction_fn]
    type = PiecewiseConstant
    axis = x
    x = '0 5.0'
    y = '0.2 0.4'
  []
[]

[Components]
  [pipe]
    type = FlowChannelGasMix
    position = '0 0 0'
    orientation = '1 0 0'
    length = 10.0
    n_elems = 50
    A = 0.2

    initial_mass_fraction = initial_mass_fraction_fn
    initial_p = ${initial_p}
    initial_T = ${initial_T}
    initial_vel = 0

    fp = fp_mix
    closures = 'closures'

    scaling_factor_rhoEA = 1e-5
  []
  [inlet]
    type = SolidWallGasMix
    input = 'pipe:in'
  []
  [outlet]
    type = SolidWallGasMix
    input = 'pipe:out'
  []
[]

[Preconditioning]
  [pc]
    type = SMP
    full = true
  []
[]

[Executioner]
  type = Transient
  scheme = bdf2

  dt = 1000
  num_steps = 5

  solve_type = NEWTON
  petsc_options_iname = '-pc_type'
  petsc_options_value = 'lu'
  nl_rel_tol = 1e-8
  nl_abs_tol = 1e-8
  nl_max_its = 15

  l_tol = 1e-3
  l_max_its = 10
[]

[Outputs]
  [exodus]
    type = Exodus
    file_base = flow_channel_gasmix
    show = 'mass_fraction p T'
  []
[]

Usage
Input Parameters
Mesh
Variables
Material Properties
Input Files
References