CoupledPressureIncompressibleMomentumSPScalarKernel

Overview

This class derives from IncompressibleMomentumSPBase and implements the steady-state residual of the momentum equation for the Path-integrated incompressible flow model, the right-hand side of Path-Integrated Incompressible Flow Model, Eq. (4). It requires a coupled variable characteristic mass flow rate and $var element = document.getElementById("moose-equation-86648603-a87e-4b54-8e93-a3b892c9938d");katex.render("N", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ coupled fluid temperature variables for each segment, given as coupled ScalarVariables. It operates on the characteristic pressure drop, $var element = document.getElementById("moose-equation-f1ddc7e9-f96c-430a-b97d-a61f08887c72");katex.render("\\Delta P_c", 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-6e0b73d6-b9ce-400f-ab5f-13da285cc31a");katex.render("0 = \\frac{1}{\\sum_{i=1}^{n} \\frac{L_i}{A_i}} \\left [ - \\Delta P_c + \\sum_{i=1}^{N} \\left ( \\frac{f_i L_i}{D_{h,i}} \\frac{\\dot{m}|\\dot{m}|}{2 \\rho_c A_i^2} + K_i \\frac{\\dot{m}|\\dot{m}|}{2 \\rho_c A_i^2} + \\rho_i g L_i \\sin{\\alpha_i} - \\Delta P_p \\right) \\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}"}});

This kernel is intended for the use case of creating a closed loop flow path with one or more complementary instances of IncompressibleMomentumSPScalarKernel and potentially a ParsedODEKernel defining the relationship between mass flow rates in the different flow paths. Thereby creating a closed system of equations for mass flow rates and characteristic pressure drop. Please note, due to the intended use case where this kernel and complementary instances of IncompressibleMomentumSPScalarKernel share mass flow rate information and an equal, but opposite pressure gradient, the sign of the pressure gradient term in the residual is opposite that of the IncompressibleMomentumSPScalarKernel. Furthermore, this kernel is steady-state by default. A transient use case also necessitates a CoupledODETimeDerivative, which includes the time derivative of a coupled variable, $var element = document.getElementById("moose-equation-41976f69-8340-4e43-8197-7176ba9a4db3");katex.render("\\frac{dv}{dt}", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ , with $var element = document.getElementById("moose-equation-5490c4f1-2f1f-4a55-9f15-8be44a41f5ba");katex.render("v", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ being the coupled reference mass flow rate, which adds the time derivative of the mass flow rate to the residual.

As a reminder, the system of variables should be defined with the "family" attribute set to SCALAR for each variable.

Input Parameters

fpThe name of the user object for fluid properties
C++ Type:UserObjectName
Controllable:No
Description:The name of the user object for fluid properties
reference_pressuresystem reference pressure [Pa]. A functor is any of the following: a variable, a functor material property, a function, a postprocessor or a number.
C++ Type:MooseFunctorName
Unit:(no unit assumed)
Controllable:No
Description:system reference pressure [Pa]. A functor is any of the following: a variable, a functor material property, a function, a postprocessor or a number.
variableThe name of the variable that this residual object operates on
C++ Type:NonlinearVariableName
Unit:(no unit assumed)
Controllable:No
Description:The name of the variable that this residual object operates on

Required Parameters

alphasComponent flow angles per segment with respect to horizontal (-pi/2 downward to pi/2 upward). Takes a vector of functors. A functor is any of the following: a variable, a functor material property, a function, a postprocessor or a number.
C++ Type:std::vector<MooseFunctorName>
Unit:(no unit assumed)
Controllable:No
Description:Component flow angles per segment with respect to horizontal (-pi/2 downward to pi/2 upward). Takes a vector of functors. A functor is any of the following: a variable, a functor material property, a function, a postprocessor or a number.
areasComponent flow areas per segment. Takes a vector of functors. A functor is any of the following: a variable, a functor material property, a function, a postprocessor or a number.
C++ Type:std::vector<MooseFunctorName>
Unit:(no unit assumed)
Controllable:No
Description:Component flow areas per segment. Takes a vector of functors. A functor is any of the following: a variable, a functor material property, a function, a postprocessor or a number.
coupled_mass_flow_ratecoupled_mass_flow_rate. Takes a scalar variable name
C++ Type:std::vector<VariableName>
Unit:(no unit assumed)
Controllable:No
Description:coupled_mass_flow_rate. Takes a scalar variable name
forms_lossesForms loss coefficients per segment. Takes a vector of functors. A functor is any of the following: a variable, a functor material property, a function, a postprocessor or a number.
C++ Type:std::vector<MooseFunctorName>
Unit:(no unit assumed)
Controllable:No
Description:Forms loss coefficients per segment. Takes a vector of functors. A functor is any of the following: a variable, a functor material property, a function, a postprocessor or a number.
g9.80665Gravitational acceleration [m/s]. A functor is any of the following: a variable, a functor material property, a function, a postprocessor or a number.
Default:9.80665
C++ Type:MooseFunctorName
Unit:(no unit assumed)
Controllable:No
Description:Gravitational acceleration [m/s]. A functor is any of the following: a variable, a functor material property, a function, a postprocessor or a number.
is_implicitFalseWhether an explicit (previous value calculation) or implicit (current value) is used
Default:False
C++ Type:bool
Controllable:No
Description:Whether an explicit (previous value calculation) or implicit (current value) is used
lengthsComponent flow lengths per segment. Takes a vector of functors. A functor is any of the following: a variable, a functor material property, a function, a postprocessor or a number.
C++ Type:std::vector<MooseFunctorName>
Unit:(no unit assumed)
Controllable:No
Description:Component flow lengths per segment. Takes a vector of functors. A functor is any of the following: a variable, a functor material property, a function, a postprocessor or a number.
perimetersComponent flow perimeters per segment. Takes a vector of functors. A functor is any of the following: a variable, a functor material property, a function, a postprocessor or a number.
C++ Type:std::vector<MooseFunctorName>
Unit:(no unit assumed)
Controllable:No
Description:Component flow perimeters per segment. Takes a vector of functors. A functor is any of the following: a variable, a functor material property, a function, a postprocessor or a number.
pump_pressuresPump pressure gains per segment [Pa]. Takes a vector of functors. A functor is any of the following: a variable, a functor material property, a function, a postprocessor or a number.
C++ Type:std::vector<MooseFunctorName>
Unit:(no unit assumed)
Controllable:No
Description:Pump pressure gains per segment [Pa]. Takes a vector of functors. A functor is any of the following: a variable, a functor material property, a function, a postprocessor or a number.
roughnessesComponent wall roughnesses per segment [m]. Takes a vector of functors. A functor is any of the following: a variable, a functor material property, a function, a postprocessor or a number.
C++ Type:std::vector<MooseFunctorName>
Unit:(no unit assumed)
Controllable:No
Description:Component wall roughnesses per segment [m]. Takes a vector of functors. A functor is any of the following: a variable, a functor material property, a function, a postprocessor or a number.
temperaturesFluid temperature in each segment of this component. Takes a list of scalar variable names
C++ Type:std::vector<VariableName>
Unit:(no unit assumed)
Controllable:No
Description:Fluid temperature in each segment of this component. Takes a list of scalar variable names

Optional Parameters

absolute_value_vector_tagsThe tags for the vectors this residual object should fill with the absolute value of the residual contribution
C++ Type:std::vector<TagName>
Controllable:No
Description:The tags for the vectors this residual object should fill with the absolute value of the residual contribution
extra_matrix_tagsThe extra tags for the matrices this Kernel should fill
C++ Type:std::vector<TagName>
Controllable:No
Description:The extra tags for the matrices this Kernel should fill
extra_vector_tagsThe extra tags for the vectors this Kernel should fill
C++ Type:std::vector<TagName>
Controllable:No
Description:The extra tags for the vectors this Kernel should fill
matrix_onlyFalseWhether this object is only doing assembly to matrices (no vectors)
Default:False
C++ Type:bool
Controllable:No
Description:Whether this object is only doing assembly to matrices (no vectors)
matrix_tagssystemThe tag for the matrices this Kernel should fill
Default:system
C++ Type:MultiMooseEnum
Options:nontime, system
Controllable:No
Description:The tag for the matrices this Kernel should fill
vector_tagsnontimeThe tag for the vectors this Kernel should fill
Default:nontime
C++ Type:MultiMooseEnum
Options:nontime, time
Controllable:No
Description:The tag for the vectors this Kernel should fill

Contribution To Tagged Field Data 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.
implicitTrueDetermines whether this object is calculated using an implicit or explicit form
Default:True
C++ Type:bool
Controllable:No
Description:Determines whether this object is calculated using an implicit or explicit form
seed0The seed for the master random number generator
Default:0
C++ Type:unsigned int
Controllable:No
Description:The seed for the master random number generator
use_displaced_meshFalseWhether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.
Default:False
C++ Type:bool
Controllable:No
Description:Whether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.

Advanced Parameters

Input Files

(modules/thermal_hydraulics/test/tests/scalarkernels/IncompressibleVerification/ADfullCP_exp.i)
(modules/thermal_hydraulics/test/tests/scalarkernels/IncompressibleVerification/ADcoupledpressure_twoseg.i)
(modules/thermal_hydraulics/test/tests/scalarkernels/IncompressibleVerification/ADformsCP.i)
(modules/thermal_hydraulics/test/tests/scalarkernels/IncompressibleVerification/ADfullCP.i)
(modules/thermal_hydraulics/test/tests/scalarkernels/IncompressibleVerification/coupledpressure.i)
(modules/thermal_hydraulics/test/tests/scalarkernels/IncompressibleVerification/coupledpressure_twoseg.i)
(modules/thermal_hydraulics/test/tests/scalarkernels/IncompressibleVerification/fullCP.i)
(modules/thermal_hydraulics/test/tests/scalarkernels/IncompressibleVerification/coupledpressure_reverse.i)
(modules/thermal_hydraulics/test/tests/scalarkernels/IncompressibleVerification/ADcoupledpressure_reverse.i)
(modules/thermal_hydraulics/test/tests/scalarkernels/IncompressibleVerification/ADpumpCP.i)
(modules/thermal_hydraulics/test/tests/scalarkernels/IncompressibleVerification/ADgravCP.i)
(modules/thermal_hydraulics/test/tests/scalarkernels/IncompressibleVerification/ADfrictionCP.i)

(modules/thermal_hydraulics/test/tests/scalarkernels/IncompressibleVerification/ADfullCP_exp.i)



length = 1.0
R = 0.025
min = 1.0
Pin = 101300
dP = 10130
Tin = 293.15
rho = 998.2
mu = 0.001002
alpha = 0.0
epsilon = 0.0001
forms = 0.0
pump = 0.0
gravity = 0.0
area = '${fparse 3.14159* ${R}^2}'

[Mesh]
  type = GeneratedMesh
  dim = 1
  xmin = 0
  xmax = ${length}
  nx = 1
[]

[Variables]
  [m1]
    family = SCALAR
    initial_condition = ${min}
  []
  [dPc]
    family = SCALAR
    initial_condition = ${dP}
  []
  [T1]
    family = SCALAR
    initial_condition = ${Tin}
  []
[]

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

[ScalarKernels]
  [pipe1_mom]
    type = CoupledPressureIncompressibleMomentumSPScalarKernel
    variable = 'dPc'
    coupled_mass_flow_rate = 'm1'
    temperatures = 'T1'
    reference_pressure = ${Pin}
    fp = 'water'
    areas = '${area}'
    perimeters = '${fparse 2*3.14159* ${R}}'
    lengths = '${length}'
    alphas = '${alpha}'
    forms_losses = '${forms}'
    pump_pressures = '${pump}'
    roughnesses = '${epsilon}'
    g = ${gravity}
    is_implicit = True
  []
  [pipe1_coupledDT]
    type = CoupledODETimeDerivative
    variable = 'dPc'
    v = 'm1'
  []
  [temp]
    type = ParsedODEKernel
    expression = 'T1 - ${Tin}'
    variable = T1
  []
  [m1k]
    type = ParsedODEKernel
    expression = 'm1 - ${min}'
    variable = m1
  []
[]

[Postprocessors]
  [mdot]
    type = ScalarVariable
    variable = m1
    execute_on = 'TIMESTEP_END'
  []
  [dP]
    type = ScalarVariable
    variable = dPc
    execute_on = 'TIMESTEP_END'
  []
  [Re]
    type = ParsedPostprocessor
    expression = 'mdot / ${area} * 2 * ${R} / ${mu}'
    pp_names = 'mdot'
    execute_on = 'TIMESTEP_END'
  []
  [f]
    type = ParsedPostprocessor
    expression = '0.25 / ((log10(${epsilon} / 2 / 3.7 / ${R} + 5.74 / Re ^ 0.9 ))^2)'
    pp_names = 'Re'
    execute_on = 'TIMESTEP_END'
  []
  [analytical_dP]
    type = ParsedPostprocessor
    expression = '-${pump} + f * ${length} * mdot^2 / 4 / ${R} / ${rho} / ${area}^2 + ${forms} * mdot^2 / 2 / ${rho} / ${area}^2 + ${rho} * ${gravity} * ${length} * sin(${alpha})'
    pp_names = 'f mdot'
    execute_on = 'TIMESTEP_END'
  []
  [relative_error]
    type = ParsedPostprocessor
    expression = 'abs((analytical_dP - dP)/analytical_dP)'
    pp_names = 'analytical_dP dP'
    execute_on = 'TIMESTEP_END'
  []
[]

[Executioner]
  type = Transient
  start_time = 0
  end_time = 100.0
  [TimeStepper]
    type = IterationAdaptiveDT
    growth_factor = 1.4
    dt = 5
  []
  solve_type = 'PJFNK'
  nl_abs_tol = 1e-09
  l_tol = 1e-07
[]

[Outputs]
  perf_graph = true
  [out]
    type = CSV
    execute_on = 'FINAL'
  []
[]

(modules/thermal_hydraulics/test/tests/scalarkernels/IncompressibleVerification/ADcoupledpressure_twoseg.i)



length1 = 1.0
length2 = 2.0
R1 = 0.025
R2 = 0.05
min = 1.0
Pin = 101300
dP = -10130
Tin = 293.15
rho = 998.2
mu = 0.001002
alpha1 = 0.0
alpha2 = '${fparse 3.14159/2}'
epsilon1 = 0.0001
epsilon2 = 0.0002
forms1 = 0.0
forms2 = 1.0
pump1 = '${fparse - ${dP}}'
pump2 = 0.0
gravity = 9.81
area1 = '${fparse 3.14159* ${R1}^2}'
area2 = '${fparse 3.14159* ${R2}^2}'

[Mesh]
  type = GeneratedMesh
  dim = 1
  xmin = 0
  xmax = ${length1}
  nx = 1
[]

[Variables]
  [m1]
    family = SCALAR
    initial_condition = ${min}
  []
  [dPc]
    family = SCALAR
    initial_condition = ${dP}
  []
  [T1]
    family = SCALAR
    initial_condition = ${Tin}
  []
[]

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

[ScalarKernels]
  [pipe1_mom]
    type = CoupledPressureIncompressibleMomentumSPScalarKernel
    variable = 'dPc'
    coupled_mass_flow_rate = 'm1'
    temperatures = 'T1 T1'
    reference_pressure = ${Pin}
    fp = 'water'
    areas = '${area1} ${area2}'
    perimeters = '${fparse 2*3.14159* ${R1}} ${fparse 2*3.14159* ${R2}}'
    lengths = '${length1} ${length2}'
    alphas = '${alpha1} ${alpha2}'
    forms_losses = '${forms1} ${forms2}'
    pump_pressures = '${pump1} ${pump2}'
    roughnesses = '${epsilon1} ${epsilon2}'
    g = ${gravity}
    is_implicit = True
  []
  [pipe1_coupledDT]
    type = CoupledODETimeDerivative
    variable = 'dPc'
    v = 'm1'
  []
  [temp]
    type = ParsedODEKernel
    expression = 'T1 - ${Tin}'
    variable = T1
  []
  [m1k]
    type = ParsedODEKernel
    expression = 'm1 - ${min}'
    variable = m1
  []
[]

[Postprocessors]
  [mdot]
    type = ScalarVariable
    variable = m1
    execute_on = 'TIMESTEP_END'
  []
  [dP]
    type = ScalarVariable
    variable = dPc
    execute_on = 'TIMESTEP_END'
  []
  [Re1]
    type = ParsedPostprocessor
    expression = 'mdot / ${area1} * 2 * ${R1} / ${mu}'
    pp_names = 'mdot'
    execute_on = 'TIMESTEP_END'
  []
  [Re2]
    type = ParsedPostprocessor
    expression = 'mdot / ${area2} * 2 * ${R2} / ${mu}'
    pp_names = 'mdot'
    execute_on = 'TIMESTEP_END'
  []
  [f1]
    type = ParsedPostprocessor
    expression = '0.25 / ((log10(${epsilon1} / 2 / 3.7 / ${R1} + 5.74 / Re1 ^ 0.9 ))^2)'
    pp_names = 'Re1'
    execute_on = 'TIMESTEP_END'
  []
  [f2]
    type = ParsedPostprocessor
    expression = '0.25 / ((log10(${epsilon2} / 2 / 3.7 / ${R2} + 5.74 / Re2 ^ 0.9 ))^2)'
    pp_names = 'Re2'
    execute_on = 'TIMESTEP_END'
  []
  [analytical_dP]
    type = ParsedPostprocessor
    expression = '- ${pump1} - ${pump2}
                  + f1 * ${length1} * mdot^2 / 4 / ${R1} / ${rho} / ${area1}^2
                  + f2 * ${length2} * mdot^2 / 4 / ${R2} / ${rho} / ${area2}^2
                  + ${forms1} * mdot^2 / 2 / ${rho} / ${area1}^2
                  + ${forms2} * mdot^2 / 2 / ${rho} / ${area2}^2
                  + ${rho} * ${gravity} * ${length1} * sin(${alpha1})
                  + ${rho} * ${gravity} * ${length2} * sin(${alpha2})'
    pp_names = 'f1 f2 mdot'
    execute_on = 'TIMESTEP_END'
  []
  [relative_error]
    type = ParsedPostprocessor
    expression = 'abs((analytical_dP - dP)/analytical_dP)'
    pp_names = 'analytical_dP dP'
    execute_on = 'TIMESTEP_END'
  []
[]

[Executioner]
  type = Transient
  start_time = 0
  end_time = 100.0
  [TimeStepper]
    type = IterationAdaptiveDT
    growth_factor = 1.4
    dt = 5
  []
  solve_type = 'PJFNK'
  nl_abs_tol = 1e-12
  l_tol = 1e-07
[]

[Outputs]
  perf_graph = true
  [out]
    type = CSV
    execute_on = 'FINAL'
  []
[]

(modules/thermal_hydraulics/test/tests/scalarkernels/IncompressibleVerification/ADformsCP.i)



length = 1.0
R = 0.025
min = 1.0
Pin = 101300
dP = 10130
Tin = 293.15
rho = 998.2
mu = 0.001002
alpha = 0.0
epsilon = 0.0001
forms = 0.0
pump = 0.0
gravity = 0.0
area = '${fparse 3.14159* ${R}^2}'

[Mesh]
  type = GeneratedMesh
  dim = 1
  xmin = 0
  xmax = ${length}
  nx = 1
[]

[Variables]
  [m1]
    family = SCALAR
    initial_condition = ${min}
  []
  [dPc]
    family = SCALAR
    initial_condition = ${dP}
  []
  [T1]
    family = SCALAR
    initial_condition = ${Tin}
  []
[]

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

[ScalarKernels]
  [pipe1_mom]
    type = CoupledPressureIncompressibleMomentumSPScalarKernel
    variable = 'dPc'
    coupled_mass_flow_rate = 'm1'
    temperatures = 'T1'
    reference_pressure = ${Pin}
    fp = 'water'
    areas = '${area}'
    perimeters = '${fparse 2*3.14159* ${R}}'
    lengths = '${length}'
    alphas = '${alpha}'
    forms_losses = '${forms}'
    pump_pressures = '${pump}'
    roughnesses = '${epsilon}'
    g = ${gravity}
    is_implicit = True
  []
  [pipe1_coupledDT]
    type = CoupledODETimeDerivative
    variable = 'dPc'
    v = 'm1'
  []
  [temp]
    type = ParsedODEKernel
    expression = 'T1 - ${Tin}'
    variable = T1
  []
  [m1k]
    type = ParsedODEKernel
    expression = 'm1 - ${min}'
    variable = m1
  []
[]

[Postprocessors]
  [mdot]
    type = ScalarVariable
    variable = m1
    execute_on = 'TIMESTEP_END'
  []
  [dP]
    type = ScalarVariable
    variable = dPc
    execute_on = 'TIMESTEP_END'
  []
  [Re]
    type = ParsedPostprocessor
    expression = 'mdot / ${area} * 2 * ${R} / ${mu}'
    pp_names = 'mdot'
    execute_on = 'TIMESTEP_END'
  []
  [f]
    type = ParsedPostprocessor
    expression = '0.25 / ((log10(${epsilon} / 2 / 3.7 / ${R} + 5.74 / Re ^ 0.9 ))^2)'
    pp_names = 'Re'
    execute_on = 'TIMESTEP_END'
  []
  [analytical_dP]
    type = ParsedPostprocessor
    expression = '-${pump} + f * ${length} * mdot^2 / 4 / ${R} / ${rho} / ${area}^2 + ${forms} * mdot^2 / 2 / ${rho} / ${area}^2 + ${rho} * ${gravity} * ${length} * sin(${alpha})'
    pp_names = 'f mdot'
    execute_on = 'TIMESTEP_END'
  []
  [relative_error]
    type = ParsedPostprocessor
    expression = 'abs((analytical_dP - dP)/analytical_dP)'
    pp_names = 'analytical_dP dP'
    execute_on = 'TIMESTEP_END'
  []
[]

[Executioner]
  type = Transient
  start_time = 0
  end_time = 100.0
  [TimeStepper]
    type = IterationAdaptiveDT
    growth_factor = 1.4
    dt = 5
  []
  solve_type = 'PJFNK'
  nl_abs_tol = 1e-09
  l_tol = 1e-07
[]

[Outputs]
  perf_graph = true
  [out]
    type = CSV
    execute_on = 'FINAL'
  []
[]

(modules/thermal_hydraulics/test/tests/scalarkernels/IncompressibleVerification/ADfullCP.i)



length = 1.0
R = 0.025
min = 1.0
Pin = 101300
dP = 10130
Tin = 293.15
rho = 998.2
mu = 0.001002
alpha = 0.0
epsilon = 0.0001
forms = 0.0
pump = 0.0
gravity = 0.0
area = '${fparse 3.14159* ${R}^2}'

[Mesh]
  type = GeneratedMesh
  dim = 1
  xmin = 0
  xmax = ${length}
  nx = 1
[]

[Variables]
  [m1]
    family = SCALAR
    initial_condition = ${min}
  []
  [dPc]
    family = SCALAR
    initial_condition = ${dP}
  []
  [T1]
    family = SCALAR
    initial_condition = ${Tin}
  []
[]

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

[ScalarKernels]
  [pipe1_mom]
    type = CoupledPressureIncompressibleMomentumSPScalarKernel
    variable = 'dPc'
    coupled_mass_flow_rate = 'm1'
    temperatures = 'T1'
    reference_pressure = ${Pin}
    fp = 'water'
    areas = '${area}'
    perimeters = '${fparse 2*3.14159* ${R}}'
    lengths = '${length}'
    alphas = '${alpha}'
    forms_losses = '${forms}'
    pump_pressures = '${pump}'
    roughnesses = '${epsilon}'
    g = ${gravity}
    is_implicit = True
  []
  [pipe1_coupledDT]
    type = CoupledODETimeDerivative
    variable = 'dPc'
    v = 'm1'
  []
  [temp]
    type = ParsedODEKernel
    expression = 'T1 - ${Tin}'
    variable = T1
  []
  [m1k]
    type = ParsedODEKernel
    expression = 'm1 - ${min}'
    variable = m1
  []
[]

[Postprocessors]
  [mdot]
    type = ScalarVariable
    variable = m1
    execute_on = 'TIMESTEP_END'
  []
  [dP]
    type = ScalarVariable
    variable = dPc
    execute_on = 'TIMESTEP_END'
  []
  [Re]
    type = ParsedPostprocessor
    expression = 'mdot / ${area} * 2 * ${R} / ${mu}'
    pp_names = 'mdot'
    execute_on = 'TIMESTEP_END'
  []
  [f]
    type = ParsedPostprocessor
    expression = '0.25 / ((log10(${epsilon} / 2 / 3.7 / ${R} + 5.74 / Re ^ 0.9 ))^2)'
    pp_names = 'Re'
    execute_on = 'TIMESTEP_END'
  []
  [analytical_dP]
    type = ParsedPostprocessor
    expression = '-${pump} + f * ${length} * mdot^2 / 4 / ${R} / ${rho} / ${area}^2 + ${forms} * mdot^2 / 2 / ${rho} / ${area}^2 + ${rho} * ${gravity} * ${length} * sin(${alpha})'
    pp_names = 'f mdot'
    execute_on = 'TIMESTEP_END'
  []
  [relative_error]
    type = ParsedPostprocessor
    expression = 'abs((analytical_dP - dP)/analytical_dP)'
    pp_names = 'analytical_dP dP'
    execute_on = 'TIMESTEP_END'
  []
[]

[Executioner]
  type = Transient
  start_time = 0
  end_time = 100.0
  [TimeStepper]
    type = IterationAdaptiveDT
    growth_factor = 1.4
    dt = 5
  []
  solve_type = 'PJFNK'
  nl_abs_tol = 1e-09
  l_tol = 1e-07
[]

[Outputs]
  perf_graph = true
  [out]
    type = CSV
    execute_on = 'FINAL'
  []
[]

(modules/thermal_hydraulics/test/tests/scalarkernels/IncompressibleVerification/coupledpressure.i)



length = 1.0
R = 0.025
min = 1.0
Pin = 101300
dP = 10130
Tin = 293.15
rho = 998.2
mu = 0.001002
alpha = 0.0
epsilon = 0.0001
forms = 0.0
pump = 0.0
gravity = 0.0
area = '${fparse 3.14159* ${R}^2}'

[Mesh]
  type = GeneratedMesh
  dim = 1
  xmin = 0
  xmax = ${length}
  nx = 1
[]

[Variables]
  [m1]
    family = SCALAR
    initial_condition = ${min}
  []
  [dPc]
    family = SCALAR
    initial_condition = ${dP}
  []
  [T1]
    family = SCALAR
    initial_condition = ${Tin}
  []
[]

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

[ScalarKernels]
  [pipe1_mom]
    type = CoupledPressureIncompressibleMomentumSPScalarKernel
    variable = 'dPc'
    coupled_mass_flow_rate = 'm1'
    temperatures = 'T1'
    reference_pressure = ${Pin}
    fp = 'water'
    areas = '${area}'
    perimeters = '${fparse 2*3.14159* ${R}}'
    lengths = '${length}'
    alphas = '${alpha}'
    forms_losses = '${forms}'
    pump_pressures = '${pump}'
    roughnesses = '${epsilon}'
    g = ${gravity}
    is_implicit = True
  []
  [pipe1_coupledDT]
    type = CoupledODETimeDerivative
    variable = 'dPc'
    v = 'm1'
  []
  [temp]
    type = ParsedODEKernel
    expression = 'T1 - ${Tin}'
    variable = T1
  []
  [m1k]
    type = ParsedODEKernel
    expression = 'm1 - ${min}'
    variable = m1
  []
[]

[Postprocessors]
  [mdot]
    type = ScalarVariable
    variable = m1
    execute_on = 'TIMESTEP_END'
  []
  [dP]
    type = ScalarVariable
    variable = dPc
    execute_on = 'TIMESTEP_END'
  []
  [Re]
    type = ParsedPostprocessor
    expression = 'mdot / ${area} * 2 * ${R} / ${mu}'
    pp_names = 'mdot'
    execute_on = 'TIMESTEP_END'
  []
  [f]
    type = ParsedPostprocessor
    expression = '0.25 / ((log10(${epsilon} / 2 / 3.7 / ${R} + 5.74 / Re ^ 0.9 ))^2)'
    pp_names = 'Re'
    execute_on = 'TIMESTEP_END'
  []
  [analytical_dP]
    type = ParsedPostprocessor
    expression = '-${pump} + f * ${length} * mdot^2 / 4 / ${R} / ${rho} / ${area}^2 + ${forms} * mdot^2 / 2 / ${rho} / ${area}^2 + ${rho} * ${gravity} * ${length} * sin(${alpha})'
    pp_names = 'f mdot'
    execute_on = 'TIMESTEP_END'
  []
  [relative_error]
    type = ParsedPostprocessor
    expression = 'abs((analytical_dP - dP)/analytical_dP)'
    pp_names = 'analytical_dP dP'
    execute_on = 'TIMESTEP_END'
  []
[]

[Executioner]
  type = Transient
  start_time = 0
  end_time = 100.0
  [TimeStepper]
    type = IterationAdaptiveDT
    growth_factor = 1.4
    dt = 5
  []
  solve_type = 'PJFNK'
  nl_abs_tol = 1e-09
  l_tol = 1e-07
[]

[Outputs]
  perf_graph = true
  [out]
    type = CSV
    execute_on = 'FINAL'
  []
[]

(modules/thermal_hydraulics/test/tests/scalarkernels/IncompressibleVerification/coupledpressure_twoseg.i)



length1 = 1.0
length2 = 2.0
R1 = 0.025
R2 = 0.05
min = 1.0
Pin = 101300
dP = -10130
Tin = 293.15
rho = 998.2
mu = 0.001002
alpha1 = 0.0
alpha2 = '${fparse 3.14159/2}'
epsilon1 = 0.0001
epsilon2 = 0.0002
forms1 = 0.0
forms2 = 1.0
pump1 = '${fparse - ${dP}}'
pump2 = 0.0
gravity = 9.81
area1 = '${fparse 3.14159* ${R1}^2}'
area2 = '${fparse 3.14159* ${R2}^2}'

[Mesh]
  type = GeneratedMesh
  dim = 1
  xmin = 0
  xmax = ${length1}
  nx = 1
[]

[Variables]
  [m1]
    family = SCALAR
    initial_condition = ${min}
  []
  [dPc]
    family = SCALAR
    initial_condition = ${dP}
  []
  [T1]
    family = SCALAR
    initial_condition = ${Tin}
  []
[]

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

[ScalarKernels]
  [pipe1_mom]
    type = CoupledPressureIncompressibleMomentumSPScalarKernel
    variable = 'dPc'
    coupled_mass_flow_rate = 'm1'
    temperatures = 'T1 T1'
    reference_pressure = ${Pin}
    fp = 'water'
    areas = '${area1} ${area2}'
    perimeters = '${fparse 2*3.14159* ${R1}} ${fparse 2*3.14159* ${R2}}'
    lengths = '${length1} ${length2}'
    alphas = '${alpha1} ${alpha2}'
    forms_losses = '${forms1} ${forms2}'
    pump_pressures = '${pump1} ${pump2}'
    roughnesses = '${epsilon1} ${epsilon2}'
    g = ${gravity}
    is_implicit = True
  []
  [pipe1_coupledDT]
    type = CoupledODETimeDerivative
    variable = 'dPc'
    v = 'm1'
  []
  [temp]
    type = ParsedODEKernel
    expression = 'T1 - ${Tin}'
    variable = T1
  []
  [m1k]
    type = ParsedODEKernel
    expression = 'm1 - ${min}'
    variable = m1
  []
[]

[Postprocessors]
  [mdot]
    type = ScalarVariable
    variable = m1
    execute_on = 'TIMESTEP_END'
  []
  [dP]
    type = ScalarVariable
    variable = dPc
    execute_on = 'TIMESTEP_END'
  []
  [Re1]
    type = ParsedPostprocessor
    expression = 'mdot / ${area1} * 2 * ${R1} / ${mu}'
    pp_names = 'mdot'
    execute_on = 'TIMESTEP_END'
  []
  [Re2]
    type = ParsedPostprocessor
    expression = 'mdot / ${area2} * 2 * ${R2} / ${mu}'
    pp_names = 'mdot'
    execute_on = 'TIMESTEP_END'
  []
  [f1]
    type = ParsedPostprocessor
    expression = '0.25 / ((log10(${epsilon1} / 2 / 3.7 / ${R1} + 5.74 / Re1 ^ 0.9 ))^2)'
    pp_names = 'Re1'
    execute_on = 'TIMESTEP_END'
  []
  [f2]
    type = ParsedPostprocessor
    expression = '0.25 / ((log10(${epsilon2} / 2 / 3.7 / ${R2} + 5.74 / Re2 ^ 0.9 ))^2)'
    pp_names = 'Re2'
    execute_on = 'TIMESTEP_END'
  []
  [analytical_dP]
    type = ParsedPostprocessor
    expression = '- ${pump1} - ${pump2}
                  + f1 * ${length1} * mdot^2 / 4 / ${R1} / ${rho} / ${area1}^2
                  + f2 * ${length2} * mdot^2 / 4 / ${R2} / ${rho} / ${area2}^2
                  + ${forms1} * mdot^2 / 2 / ${rho} / ${area1}^2
                  + ${forms2} * mdot^2 / 2 / ${rho} / ${area2}^2
                  + ${rho} * ${gravity} * ${length1} * sin(${alpha1})
                  + ${rho} * ${gravity} * ${length2} * sin(${alpha2})'
    pp_names = 'f1 f2 mdot'
    execute_on = 'TIMESTEP_END'
  []
  [relative_error]
    type = ParsedPostprocessor
    expression = 'abs((analytical_dP - dP)/analytical_dP)'
    pp_names = 'analytical_dP dP'
    execute_on = 'TIMESTEP_END'
  []
[]

[Executioner]
  type = Transient
  start_time = 0
  end_time = 100.0
  [TimeStepper]
    type = IterationAdaptiveDT
    growth_factor = 1.4
    dt = 5
  []
  solve_type = 'PJFNK'
  nl_abs_tol = 1e-09
  l_tol = 1e-07
[]

[Outputs]
  perf_graph = true
  [out]
    type = CSV
    execute_on = 'FINAL'
  []
[]

(modules/thermal_hydraulics/test/tests/scalarkernels/IncompressibleVerification/fullCP.i)



length = 1.0
R = 0.025
min = 1.0
Pin = 101300
dP = 10130
Tin = 293.15
rho = 998.2
mu = 0.001002
alpha = 0.0
epsilon = 0.0001
forms = 0.0
pump = 0.0
gravity = 0.0
area = '${fparse 3.14159* ${R}^2}'

[Mesh]
  type = GeneratedMesh
  dim = 1
  xmin = 0
  xmax = ${length}
  nx = 1
[]

[Variables]
  [m1]
    family = SCALAR
    initial_condition = ${min}
  []
  [dPc]
    family = SCALAR
    initial_condition = ${dP}
  []
  [T1]
    family = SCALAR
    initial_condition = ${Tin}
  []
[]

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

[ScalarKernels]
  [pipe1_mom]
    type = CoupledPressureIncompressibleMomentumSPScalarKernel
    variable = 'dPc'
    coupled_mass_flow_rate = 'm1'
    temperatures = 'T1'
    reference_pressure = ${Pin}
    fp = 'water'
    areas = '${area}'
    perimeters = '${fparse 2*3.14159* ${R}}'
    lengths = '${length}'
    alphas = '${alpha}'
    forms_losses = '${forms}'
    pump_pressures = '${pump}'
    roughnesses = '${epsilon}'
    g = ${gravity}
    is_implicit = True
  []
  [pipe1_coupledDT]
    type = CoupledODETimeDerivative
    variable = 'dPc'
    v = 'm1'
  []
  [temp]
    type = ParsedODEKernel
    expression = 'T1 - ${Tin}'
    variable = T1
  []
  [m1k]
    type = ParsedODEKernel
    expression = 'm1 - ${min}'
    variable = m1
  []
[]

[Postprocessors]
  [mdot]
    type = ScalarVariable
    variable = m1
    execute_on = 'TIMESTEP_END'
  []
  [dP]
    type = ScalarVariable
    variable = dPc
    execute_on = 'TIMESTEP_END'
  []
  [Re]
    type = ParsedPostprocessor
    expression = 'mdot / ${area} * 2 * ${R} / ${mu}'
    pp_names = 'mdot'
    execute_on = 'TIMESTEP_END'
  []
  [f]
    type = ParsedPostprocessor
    expression = '0.25 / ((log10(${epsilon} / 2 / 3.7 / ${R} + 5.74 / Re ^ 0.9 ))^2)'
    pp_names = 'Re'
    execute_on = 'TIMESTEP_END'
  []
  [analytical_dP]
    type = ParsedPostprocessor
    expression = '-${pump} + f * ${length} * mdot^2 / 4 / ${R} / ${rho} / ${area}^2 + ${forms} * mdot^2 / 2 / ${rho} / ${area}^2 + ${rho} * ${gravity} * ${length} * sin(${alpha})'
    pp_names = 'f mdot'
    execute_on = 'TIMESTEP_END'
  []
  [relative_error]
    type = ParsedPostprocessor
    expression = 'abs((analytical_dP - dP)/analytical_dP)'
    pp_names = 'analytical_dP dP'
    execute_on = 'TIMESTEP_END'
  []
[]

[Executioner]
  type = Transient
  start_time = 0
  end_time = 100.0
  [TimeStepper]
    type = IterationAdaptiveDT
    growth_factor = 1.4
    dt = 5
  []
  solve_type = 'PJFNK'
  nl_abs_tol = 1e-09
  l_tol = 1e-07
[]

[Outputs]
  perf_graph = true
  [out]
    type = CSV
    execute_on = 'FINAL'
  []
[]

(modules/thermal_hydraulics/test/tests/scalarkernels/IncompressibleVerification/coupledpressure_reverse.i)



length = 1.0
R = 0.025
min = 1.0
Pin = 101300
dP = 10130
Tin = 293.15
rho = 998.2
mu = 0.001002
alpha = 0.0
epsilon = 0.0001
forms = 0.0
pump = 0.0
gravity = 0.0
area = '${fparse 3.14159* ${R}^2}'

[Mesh]
  type = GeneratedMesh
  dim = 1
  xmin = 0
  xmax = ${length}
  nx = 1
[]

[Variables]
  [m1]
    family = SCALAR
    initial_condition = ${min}
  []
  [dPc]
    family = SCALAR
    initial_condition = ${dP}
  []
  [T1]
    family = SCALAR
    initial_condition = ${Tin}
  []
[]

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

[ScalarKernels]
  [pipe1_mom]
    type = CoupledPressureIncompressibleMomentumSPScalarKernel
    variable = 'dPc'
    coupled_mass_flow_rate = 'm1'
    temperatures = 'T1'
    reference_pressure = ${Pin}
    fp = 'water'
    areas = '${area}'
    perimeters = '${fparse 2*3.14159* ${R}}'
    lengths = '${length}'
    alphas = '${alpha}'
    forms_losses = '${forms}'
    pump_pressures = '${pump}'
    roughnesses = '${epsilon}'
    g = ${gravity}
    is_implicit = True
  []
  [pipe1_coupledDT]
    type = CoupledODETimeDerivative
    variable = 'dPc'
    v = 'm1'
  []
  [temp]
    type = ParsedODEKernel
    expression = 'T1 - ${Tin}'
    variable = T1
  []
  [m1k]
    type = ParsedODEKernel
    expression = 'm1 - ${min}'
    variable = m1
  []
[]

[Postprocessors]
  [mdot]
    type = ScalarVariable
    variable = m1
    execute_on = 'TIMESTEP_END'
  []
  [dP]
    type = ScalarVariable
    variable = dPc
    execute_on = 'TIMESTEP_END'
  []
  [Re]
    type = ParsedPostprocessor
    expression = 'abs(mdot) / ${area} * 2 * ${R} / ${mu}'
    pp_names = 'mdot'
    execute_on = 'TIMESTEP_END'
  []
  [f]
    type = ParsedPostprocessor
    expression = '0.25 / ((log10(${epsilon} / 2 / 3.7 / ${R} + 5.74 / Re ^ 0.9 ))^2)'
    pp_names = 'Re'
    execute_on = 'TIMESTEP_END'
  []
  [analytical_dP]
    type = ParsedPostprocessor
    expression = '-${pump} + f * ${length} * mdot * abs(mdot) / 4 / ${R} / ${rho} / ${area}^2 + ${forms} * mdot * abs(mdot) / 2 / ${rho} / ${area}^2 + ${rho} * ${gravity} * ${length} * sin(${alpha})'
    pp_names = 'f mdot'
    execute_on = 'TIMESTEP_END'
  []
  [relative_error]
    type = ParsedPostprocessor
    expression = 'abs((analytical_dP - dP)/analytical_dP)'
    pp_names = 'analytical_dP dP'
    execute_on = 'TIMESTEP_END'
  []
[]

[Executioner]
  type = Transient
  start_time = 0
  end_time = 100.0
  [TimeStepper]
    type = IterationAdaptiveDT
    growth_factor = 1.4
    dt = 5
  []
  solve_type = 'PJFNK'
  nl_abs_tol = 1e-09
  l_tol = 1e-07
[]

[Outputs]
  perf_graph = true
  [out]
    type = CSV
    execute_on = 'FINAL'
  []
[]

(modules/thermal_hydraulics/test/tests/scalarkernels/IncompressibleVerification/ADcoupledpressure_reverse.i)



length = 1.0
R = 0.025
min = 1.0
Pin = 101300
dP = 10130
Tin = 293.15
rho = 998.2
mu = 0.001002
alpha = 0.0
epsilon = 0.0001
forms = 0.0
pump = 0.0
gravity = 0.0
area = '${fparse 3.14159* ${R}^2}'

[Mesh]
  type = GeneratedMesh
  dim = 1
  xmin = 0
  xmax = ${length}
  nx = 1
[]

[Variables]
  [m1]
    family = SCALAR
    initial_condition = ${min}
  []
  [dPc]
    family = SCALAR
    initial_condition = ${dP}
  []
  [T1]
    family = SCALAR
    initial_condition = ${Tin}
  []
[]

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

[ScalarKernels]
  [pipe1_mom]
    type = CoupledPressureIncompressibleMomentumSPScalarKernel
    variable = 'dPc'
    coupled_mass_flow_rate = 'm1'
    temperatures = 'T1'
    reference_pressure = ${Pin}
    fp = 'water'
    areas = '${area}'
    perimeters = '${fparse 2*3.14159* ${R}}'
    lengths = '${length}'
    alphas = '${alpha}'
    forms_losses = '${forms}'
    pump_pressures = '${pump}'
    roughnesses = '${epsilon}'
    g = ${gravity}
    is_implicit = True
  []
  [pipe1_coupledDT]
    type = CoupledODETimeDerivative
    variable = 'dPc'
    v = 'm1'
  []
  [temp]
    type = ParsedODEKernel
    expression = 'T1 - ${Tin}'
    variable = T1
  []
  [m1k]
    type = ParsedODEKernel
    expression = 'm1 - ${min}'
    variable = m1
  []
[]

[Postprocessors]
  [mdot]
    type = ScalarVariable
    variable = m1
    execute_on = 'TIMESTEP_END'
  []
  [dP]
    type = ScalarVariable
    variable = dPc
    execute_on = 'TIMESTEP_END'
  []
  [Re]
    type = ParsedPostprocessor
    expression = 'abs(mdot) / ${area} * 2 * ${R} / ${mu}'
    pp_names = 'mdot'
    execute_on = 'TIMESTEP_END'
  []
  [f]
    type = ParsedPostprocessor
    expression = '0.25 / ((log10(${epsilon} / 2 / 3.7 / ${R} + 5.74 / Re ^ 0.9 ))^2)'
    pp_names = 'Re'
    execute_on = 'TIMESTEP_END'
  []
  [analytical_dP]
    type = ParsedPostprocessor
    expression = '-${pump} + f * ${length} * mdot * abs(mdot) / 4 / ${R} / ${rho} / ${area}^2 + ${forms} * mdot * abs(mdot) / 2 / ${rho} / ${area}^2 + ${rho} * ${gravity} * ${length} * sin(${alpha})'
    pp_names = 'f mdot'
    execute_on = 'TIMESTEP_END'
  []
  [relative_error]
    type = ParsedPostprocessor
    expression = 'abs((analytical_dP - dP)/analytical_dP)'
    pp_names = 'analytical_dP dP'
    execute_on = 'TIMESTEP_END'
  []
[]

[Executioner]
  type = Transient
  start_time = 0
  end_time = 100.0
  [TimeStepper]
    type = IterationAdaptiveDT
    growth_factor = 1.4
    dt = 5
  []
  solve_type = 'PJFNK'
  nl_abs_tol = 1e-12
  l_tol = 1e-07
[]

[Outputs]
  perf_graph = true
  [out]
    type = CSV
    execute_on = 'FINAL'
  []
[]

(modules/thermal_hydraulics/test/tests/scalarkernels/IncompressibleVerification/ADpumpCP.i)



length = 1.0
R = 0.025
min = 1.0
Pin = 101300
dP = 10130
Tin = 293.15
rho = 998.2
mu = 0.001002
alpha = 0.0
epsilon = 0.0001
forms = 0.0
pump = 0.0
gravity = 0.0
area = '${fparse 3.14159* ${R}^2}'

[Mesh]
  type = GeneratedMesh
  dim = 1
  xmin = 0
  xmax = ${length}
  nx = 1
[]

[Variables]
  [m1]
    family = SCALAR
    initial_condition = ${min}
  []
  [dPc]
    family = SCALAR
    initial_condition = ${dP}
  []
  [T1]
    family = SCALAR
    initial_condition = ${Tin}
  []
[]

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

[ScalarKernels]
  [pipe1_mom]
    type = CoupledPressureIncompressibleMomentumSPScalarKernel
    variable = 'dPc'
    coupled_mass_flow_rate = 'm1'
    temperatures = 'T1'
    reference_pressure = ${Pin}
    fp = 'water'
    areas = '${area}'
    perimeters = '${fparse 2*3.14159* ${R}}'
    lengths = '${length}'
    alphas = '${alpha}'
    forms_losses = '${forms}'
    pump_pressures = '${pump}'
    roughnesses = '${epsilon}'
    g = ${gravity}
    is_implicit = True
  []
  [pipe1_coupledDT]
    type = CoupledODETimeDerivative
    variable = 'dPc'
    v = 'm1'
  []
  [temp]
    type = ParsedODEKernel
    expression = 'T1 - ${Tin}'
    variable = T1
  []
  [m1k]
    type = ParsedODEKernel
    expression = 'm1 - ${min}'
    variable = m1
  []
[]

[Postprocessors]
  [mdot]
    type = ScalarVariable
    variable = m1
    execute_on = 'TIMESTEP_END'
  []
  [dP]
    type = ScalarVariable
    variable = dPc
    execute_on = 'TIMESTEP_END'
  []
  [Re]
    type = ParsedPostprocessor
    expression = 'mdot / ${area} * 2 * ${R} / ${mu}'
    pp_names = 'mdot'
    execute_on = 'TIMESTEP_END'
  []
  [f]
    type = ParsedPostprocessor
    expression = '0.25 / ((log10(${epsilon} / 2 / 3.7 / ${R} + 5.74 / Re ^ 0.9 ))^2)'
    pp_names = 'Re'
    execute_on = 'TIMESTEP_END'
  []
  [analytical_dP]
    type = ParsedPostprocessor
    expression = '-${pump} + f * ${length} * mdot^2 / 4 / ${R} / ${rho} / ${area}^2 + ${forms} * mdot^2 / 2 / ${rho} / ${area}^2 + ${rho} * ${gravity} * ${length} * sin(${alpha})'
    pp_names = 'f mdot'
    execute_on = 'TIMESTEP_END'
  []
  [relative_error]
    type = ParsedPostprocessor
    expression = 'abs((analytical_dP - dP)/analytical_dP)'
    pp_names = 'analytical_dP dP'
    execute_on = 'TIMESTEP_END'
  []
[]

[Executioner]
  type = Transient
  start_time = 0
  end_time = 100.0
  [TimeStepper]
    type = IterationAdaptiveDT
    growth_factor = 1.4
    dt = 5
  []
  solve_type = 'PJFNK'
  nl_abs_tol = 1e-09
  l_tol = 1e-07
[]

[Outputs]
  perf_graph = true
  [out]
    type = CSV
    execute_on = 'FINAL'
  []
[]

(modules/thermal_hydraulics/test/tests/scalarkernels/IncompressibleVerification/ADgravCP.i)



length = 1.0
R = 0.025
min = 1.0
Pin = 101300
dP = 10130
Tin = 293.15
rho = 998.2
mu = 0.001002
alpha = 0.0
epsilon = 0.0001
forms = 0.0
pump = 0.0
gravity = 0.0
area = '${fparse 3.14159* ${R}^2}'

[Mesh]
  type = GeneratedMesh
  dim = 1
  xmin = 0
  xmax = ${length}
  nx = 1
[]

[Variables]
  [m1]
    family = SCALAR
    initial_condition = ${min}
  []
  [dPc]
    family = SCALAR
    initial_condition = ${dP}
  []
  [T1]
    family = SCALAR
    initial_condition = ${Tin}
  []
[]

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

[ScalarKernels]
  [pipe1_mom]
    type = CoupledPressureIncompressibleMomentumSPScalarKernel
    variable = 'dPc'
    coupled_mass_flow_rate = 'm1'
    temperatures = 'T1'
    reference_pressure = ${Pin}
    fp = 'water'
    areas = '${area}'
    perimeters = '${fparse 2*3.14159* ${R}}'
    lengths = '${length}'
    alphas = '${alpha}'
    forms_losses = '${forms}'
    pump_pressures = '${pump}'
    roughnesses = '${epsilon}'
    g = ${gravity}
    is_implicit = True
  []
  [pipe1_coupledDT]
    type = CoupledODETimeDerivative
    variable = 'dPc'
    v = 'm1'
  []
  [temp]
    type = ParsedODEKernel
    expression = 'T1 - ${Tin}'
    variable = T1
  []
  [m1k]
    type = ParsedODEKernel
    expression = 'm1 - ${min}'
    variable = m1
  []
[]

[Postprocessors]
  [mdot]
    type = ScalarVariable
    variable = m1
    execute_on = 'TIMESTEP_END'
  []
  [dP]
    type = ScalarVariable
    variable = dPc
    execute_on = 'TIMESTEP_END'
  []
  [Re]
    type = ParsedPostprocessor
    expression = 'mdot / ${area} * 2 * ${R} / ${mu}'
    pp_names = 'mdot'
    execute_on = 'TIMESTEP_END'
  []
  [f]
    type = ParsedPostprocessor
    expression = '0.25 / ((log10(${epsilon} / 2 / 3.7 / ${R} + 5.74 / Re ^ 0.9 ))^2)'
    pp_names = 'Re'
    execute_on = 'TIMESTEP_END'
  []
  [analytical_dP]
    type = ParsedPostprocessor
    expression = '-${pump} + f * ${length} * mdot^2 / 4 / ${R} / ${rho} / ${area}^2 + ${forms} * mdot^2 / 2 / ${rho} / ${area}^2 + ${rho} * ${gravity} * ${length} * sin(${alpha})'
    pp_names = 'f mdot'
    execute_on = 'TIMESTEP_END'
  []
  [relative_error]
    type = ParsedPostprocessor
    expression = 'abs((analytical_dP - dP)/analytical_dP)'
    pp_names = 'analytical_dP dP'
    execute_on = 'TIMESTEP_END'
  []
[]

[Executioner]
  type = Transient
  start_time = 0
  end_time = 100.0
  [TimeStepper]
    type = IterationAdaptiveDT
    growth_factor = 1.4
    dt = 5
  []
  solve_type = 'PJFNK'
  nl_abs_tol = 1e-09
  l_tol = 1e-07
[]

[Outputs]
  perf_graph = true
  [out]
    type = CSV
    execute_on = 'FINAL'
  []
[]

(modules/thermal_hydraulics/test/tests/scalarkernels/IncompressibleVerification/ADfrictionCP.i)



length = 1.0
R = 0.025
min = 1.0
Pin = 101300
dP = 10130
Tin = 293.15
rho = 998.2
mu = 0.001002
alpha = 0.0
epsilon = 0.0001
forms = 0.0
pump = 0.0
gravity = 0.0
area = '${fparse 3.14159* ${R}^2}'

[Mesh]
  type = GeneratedMesh
  dim = 1
  xmin = 0
  xmax = ${length}
  nx = 1
[]

[Variables]
  [m1]
    family = SCALAR
    initial_condition = ${min}
  []
  [dPc]
    family = SCALAR
    initial_condition = ${dP}
  []
  [T1]
    family = SCALAR
    initial_condition = ${Tin}
  []
[]

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

[ScalarKernels]
  [pipe1_mom]
    type = CoupledPressureIncompressibleMomentumSPScalarKernel
    variable = 'dPc'
    coupled_mass_flow_rate = 'm1'
    temperatures = 'T1'
    reference_pressure = ${Pin}
    fp = 'water'
    areas = '${area}'
    perimeters = '${fparse 2*3.14159* ${R}}'
    lengths = '${length}'
    alphas = '${alpha}'
    forms_losses = '${forms}'
    pump_pressures = '${pump}'
    roughnesses = '${epsilon}'
    g = ${gravity}
    is_implicit = True
  []
  [pipe1_coupledDT]
    type = CoupledODETimeDerivative
    variable = 'dPc'
    v = 'm1'
  []
  [temp]
    type = ParsedODEKernel
    expression = 'T1 - ${Tin}'
    variable = T1
  []
  [m1k]
    type = ParsedODEKernel
    expression = 'm1 - ${min}'
    variable = m1
  []
[]

[Postprocessors]
  [mdot]
    type = ScalarVariable
    variable = m1
    execute_on = 'TIMESTEP_END'
  []
  [dP]
    type = ScalarVariable
    variable = dPc
    execute_on = 'TIMESTEP_END'
  []
  [Re]
    type = ParsedPostprocessor
    expression = 'mdot / ${area} * 2 * ${R} / ${mu}'
    pp_names = 'mdot'
    execute_on = 'TIMESTEP_END'
  []
  [f]
    type = ParsedPostprocessor
    expression = '0.25 / ((log10(${epsilon} / 2 / 3.7 / ${R} + 5.74 / Re ^ 0.9 ))^2)'
    pp_names = 'Re'
    execute_on = 'TIMESTEP_END'
  []
  [analytical_dP]
    type = ParsedPostprocessor
    expression = '-${pump} + f * ${length} * mdot^2 / 4 / ${R} / ${rho} / ${area}^2 + ${forms} * mdot^2 / 2 / ${rho} / ${area}^2 + ${rho} * ${gravity} * ${length} * sin(${alpha})'
    pp_names = 'f mdot'
    execute_on = 'TIMESTEP_END'
  []
  [relative_error]
    type = ParsedPostprocessor
    expression = 'abs((analytical_dP - dP)/analytical_dP)'
    pp_names = 'analytical_dP dP'
    execute_on = 'TIMESTEP_END'
  []
[]

[Executioner]
  type = Transient
  start_time = 0
  end_time = 100.0
  [TimeStepper]
    type = IterationAdaptiveDT
    growth_factor = 1.4
    dt = 5
  []
  solve_type = 'PJFNK'
  nl_abs_tol = 1e-09
  l_tol = 1e-07
[]

[Outputs]
  perf_graph = true
  [out]
    type = CSV
    execute_on = 'FINAL'
  []
[]

Overview
Input Parameters
Input Files