VolumeJunctionCoupledFlux1Phase

This component couples a VolumeJunction1Phase to another application via a MultiApp by applying a numerical flux computed using the junction state and the external state.

The junction state $var element = document.getElementById("moose-equation-a7bda03b-98b9-4a61-a70d-0df5c7d991f5");katex.render("\\mathbf{U}_\\text{J}", 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

var element = document.getElementById("moose-equation-6d52e488-884a-48cd-8af5-5b558c852a71");katex.render("\\mathbf{U} = \\left[\\begin{array}{c} \\rho_\\text{J} A_\\text{coupled}\\\\ \\rho_\\text{J} u_\\text{J} A_\\text{coupled}\\\\ \\rho_\\text{J} v_\\text{J} A_\\text{coupled}\\\\ \\rho_\\text{J} w_\\text{J} A_\\text{coupled}\\\\ \\rho_\\text{J} E_\\text{J} A_\\text{coupled} \\end{array}\\right]", 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-bd396d80-6103-47e0-aeef-1989ec0b8e07");katex.render("A_\\text{coupled}", 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 coupling area; see VolumeJunction1Phase for the definition of junction quantities.

The external state $var element = document.getElementById("moose-equation-99395015-827f-475c-9e47-8a3a65a5ceb8");katex.render("\\mathbf{U}_\\text{ext}", 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 defined by a pressure $var element = document.getElementById("moose-equation-6e98bf59-6bfb-4f31-af63-986a0bd68241");katex.render("p_\\text{ext}", 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 a temperature $var element = document.getElementById("moose-equation-28d48e2b-8344-408e-85cf-bbc37acebf56");katex.render("T_\\text{ext}", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ , while assuming the velocity is zero:

var element = document.getElementById("moose-equation-b03711f3-e0ad-415c-98ef-dd60b3aa5266");katex.render("\\mathbf{U}_\\text{ext} = \\left[\\begin{array}{c} \\rho_\\text{ext} A_\\text{coupled}\\\\ 0\\\\ 0\\\\ 0\\\\ \\rho_\\text{ext} e_\\text{ext} A_\\text{coupled} \\end{array}\\right]", 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-1cec67ec-fbe8-4097-b151-9f088f7fe95f");katex.render("\\rho_\\text{ext} = \\rho(p_\\text{ext}, T_\\text{ext})", 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-938e40e7-08a8-4ce8-9351-576f685a6872");katex.render("e_\\text{ext} = e(p_\\text{ext}, T_\\text{ext})", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ .

The numerical flux vector is computed as

var element = document.getElementById("moose-equation-57e69f0b-f260-4d13-888f-72b6f4ae6e61");katex.render("\\mathbf{F} = \\mathcal{F}(\\mathbf{U}_\\text{J}, \\mathbf{U}_\\text{ext}, \\mathbf{n}_\\text{ext,J})", 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-ccb58e0f-e711-47a9-aa6b-5fb39b9d659a");katex.render("\\mathcal{F}", 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 3-D numerical flux function and $var element = document.getElementById("moose-equation-3fb049b7-7b9f-4b43-b916-eb9e9ee52689");katex.render("\\mathbf{n}_\\text{ext,J}", 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 normal vector from the exterior application to the junction.

The following table gives the post-processors created by this component, where <suffix> is the suffix provided by "pp_suffix":

Name	Description	Type	Transfer Direction
`mass_rate_<suffix>`	$var element = document.getElementById("moose-equation-70724990-229b-4511-8eb2-55146e9caf4f");katex.render("F_\\text{mass}", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$	VolumeJunctionCoupledFlux1PhasePostprocessor	To other app
`energy_rate_<suffix>`	$var element = document.getElementById("moose-equation-a6f280ff-12b3-4044-9be0-56594ee8a236");katex.render("F_\\text{energy}", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$	VolumeJunctionCoupledFlux1PhasePostprocessor	To other app
`p_<suffix>`	$var element = document.getElementById("moose-equation-ecfe9f19-7f82-4bc3-bb80-e3add6519686");katex.render("p_\\text{ext}", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$	Receiver	From other app
`T_<suffix>`	$var element = document.getElementById("moose-equation-446ee7a2-7ebf-47b4-afdd-8675dc2952fa");katex.render("T_\\text{ext}", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$	Receiver	From other app

If this application is the main application, then the "multi_app" parameter should be provided, which creates MultiAppPostprocessorTransfer objects; otherwise the other application is responsible for creating the transfers.

Input Parameters

A_coupledArea of the coupled surface [m^2]
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Area of the coupled surface [m^2]
normal_from_junctionDirection vector pointing from the junction to the coupled volume. This vector will be normalized.
C++ Type:libMesh::VectorValue<double>
Unit:(no unit assumed)
Controllable:No
Description:Direction vector pointing from the junction to the coupled volume. This vector will be normalized.
pp_suffixSuffix to append to post-processor names
C++ Type:std::string
Controllable:No
Description:Suffix to append to post-processor names
volume_junctionName of volume junction component on which to apply flux
C++ Type:std::string
Controllable:No
Description:Name of volume junction component on which to apply flux

Required Parameters

multi_appIf provided, transfers occur with this MultiApp. The following would be transferred: mass flow rate, energy flow rate, temperature, and pressure at the junction.
C++ Type:MultiAppName
Controllable:No
Description:If provided, transfers occur with this MultiApp. The following would be transferred: mass flow rate, energy flow rate, temperature, and pressure at the junction.

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

Input Files

(modules/thermal_hydraulics/test/tests/components/volume_junction_coupled_flux_1phase/main.i)

References

No citations exist within this document.

(modules/thermal_hydraulics/test/tests/components/volume_junction_coupled_flux_1phase/main.i)

outlet_p = 1e5
initial_T = 300
R_pipe = 0.1
A_pipe = ${fparse pi * R_pipe^2}
A_coupled = ${A_pipe}
V_junction = ${fparse 4/3 * pi * R_pipe^3}

[GlobalParams]
  gravity_vector = '0 0 0'
  scaling_factor_1phase = '1 1 1e-5'

  initial_T = ${initial_T}
  initial_p = ${outlet_p}
  initial_vel = 0
  initial_vel_x = 0
  initial_vel_y = 0
  initial_vel_z = 0
[]

[FluidProperties]
  [fp]
    type = IdealGasFluidProperties
  []
[]

[Closures]
  [simple_closures]
    type = Closures1PhaseSimple
  []
[]

[Components]
  [outlet]
    type = Outlet1Phase
    input = 'pipe:in'
    p = ${outlet_p}
  []
  [pipe]
    type = FlowChannel1Phase
    position = '0 0 0'
    orientation = '1 0 0'
    length = 1.0
    n_elems = 5
    A = ${A_pipe}
    fp = fp
    closures = simple_closures
    f = 0
  []
  [junction]
    type = VolumeJunction1Phase
    connections = 'pipe:out'
    position = '1.0 0 0'
    volume = ${V_junction}
  []
  [junction_flux]
    type = VolumeJunctionCoupledFlux1Phase
    volume_junction = junction
    A_coupled = ${A_coupled}
    pp_suffix = test
    multi_app = child
    normal_from_junction = '1 0 0'
  []
[]

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

[MultiApps]
  [child]
    type = TransientMultiApp
    app_type = ThermalHydraulicsApp
    input_files = child.i
    execute_on = 'INITIAL TIMESTEP_END'
  []
[]

[Executioner]
  type = Transient
  scheme = bdf2
  end_time = 0.5
  dt = 0.1

  solve_type = NEWTON
  nl_rel_tol = 1e-8
  nl_abs_tol = 1e-8
  nl_max_its = 10

  l_tol = 1e-3
  l_max_its = 10

  petsc_options_iname = '-pc_type'
  petsc_options_value = 'lu'
[]

[Outputs]
  exodus = true
  csv = true
[]

Input Parameters
Input Files
References