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-cb69e573-35b4-4339-8123-fa541900e0f8");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-fb41c208-8be1-41bf-9169-0030ab925db2");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-312a8c59-51b4-4f60-90d3-79e867c5ac96");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-ee9f9fa7-5697-49aa-a144-d4ce8bc44fc8");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-b2aea6f3-b6bb-4afc-acbb-f7c10846d3e9");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-53273221-7399-4872-98d2-3691c90886bc");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-05dcd0d6-3ee6-419a-a03b-d73f6827fad8");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-33beac7b-81b8-4d51-a3b6-a788d85928b1");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-45363944-95fa-4405-a88b-e60b75c6e001");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-6061e16e-6f60-4ba8-a2a1-aff2a0193236");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-43e33548-d135-4528-96de-b98468301b3b");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-23529cbe-d7a4-4424-8244-3d3e5c1d51be");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-63c1ab3a-738a-4f6f-9d4f-fbce29d04d4b");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 child app
`energy_rate_<suffix>`	$var element = document.getElementById("moose-equation-695d033e-7763-4284-a190-e4a96d409a2b");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 child app
`p_<suffix>`	$var element = document.getElementById("moose-equation-7773e5c2-b3ad-4f2c-b0e9-bb439eadc01c");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 child app
`T_<suffix>`	$var element = document.getElementById("moose-equation-c78e240f-8ecf-460c-a599-2ee9ab62b261");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 child app

MultiAppPostprocessorTransfer objects are used to transfer the post-processor values in both directions.

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]
multi_appMultiApp to transfer to and from
C++ Type:MultiAppName
Controllable:No
Description:MultiApp to transfer to and from
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

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