HydrostaticPressureFunction

This function computes the hydrostatic pressure distribution in a continuous fluid, given a reference point and pressure value:

var element = document.getElementById("moose-equation-1d841dd0-b8fc-4f35-80e6-eb1c05896bf6");katex.render("p(\\mathbf{r}) = p_\\text{ref} - \\rho_\\text{ref} \\mathbf{g} (\\mathbf{r}_\\text{ref} - \\mathbf{r})", 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-758ceeb1-7599-431a-86f9-41acdaa4419a");katex.render("\\mathbf{r}_\\text{ref}", 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 reference point,
$var element = document.getElementById("moose-equation-0ce83279-e325-47e2-b4fa-b7553611b1fe");katex.render("p_\\text{ref}", 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 pressure at that point,
$var element = document.getElementById("moose-equation-c81c647b-a0f7-48e9-b4ad-72c9a7302511");katex.render("T_\\text{ref}", 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 temperature at that point,
$var element = document.getElementById("moose-equation-381b4a81-6d9b-4d1a-b1c1-5ed32697a15f");katex.render("\\rho_\\text{ref} = \\rho(p_\\text{ref}, T_\\text{ref})", 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 reference density, and
$var element = document.getElementById("moose-equation-0f3cd9f3-8f28-4da6-b000-deafb516940d");katex.render("\\mathbf{g}", 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 gravity acceleration vector.

warning

Note that this formulation is only exact if the density does not change between $var element = document.getElementById("moose-equation-15e6da15-619e-4207-b9cf-2f1a1636f44a");katex.render("\\mathbf{r}_\\text{ref}", 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-69753c89-524d-4392-b1ab-b23852862d14");katex.render("\\mathbf{r}", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ .

Input Parameters

fluid_propertiesName of the SinglePhaseFluidProperties user object
C++ Type:UserObjectName
Controllable:No
Description:Name of the SinglePhaseFluidProperties user object
reference_pointReference point where pressure is specified [m]
C++ Type:libMesh::Point
Controllable:No
Description:Reference point where pressure is specified [m]
reference_pressurePressure at the reference point [Pa]
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Pressure at the reference point [Pa]
reference_temperatureReference temperature for density evaluation [K]
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Reference temperature for density evaluation [K]

Required Parameters

gravity_directionDirection of gravitational acceleration. This will be normalized and multiplied by 'gravity_magnitude'. This parameter is mutually exclusive with 'gravity_vector'.
C++ Type:libMesh::VectorValue<double>
Unit:(no unit assumed)
Controllable:No
Description:Direction of gravitational acceleration. This will be normalized and multiplied by 'gravity_magnitude'. This parameter is mutually exclusive with 'gravity_vector'.
gravity_magnitude9.80665Magnitude of the gravitational acceleration [m/s^2]
Default:9.80665
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Magnitude of the gravitational acceleration [m/s^2]
gravity_vectorGravitational acceleration vector [m/s^2]. This parameter is mutually exclusive with 'gravity_direction'.
C++ Type:libMesh::VectorValue<double>
Unit:(no unit assumed)
Controllable:No
Description:Gravitational acceleration vector [m/s^2]. This parameter is mutually exclusive with 'gravity_direction'.

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/fluid_properties/test/tests/functions/hydrostatic_pressure_function/hydrostatic_pressure_function.i)

(modules/fluid_properties/test/tests/functions/hydrostatic_pressure_function/hydrostatic_pressure_function.i)

# rho_ref @ (p_ref, T_ref) = 1.163115975 kg/m^3
# g = 9.80665 m/s^2
#
# Expected values in gold file:
# p(0 m)    = 1e5 Pa
# p(500 m)  = 1e5 + rho_ref * g * 500  = 105703.13563811687 Pa
# p(1000 m) = 1e5 + rho_ref * g * 1000 = 111406.27127623375 Pa

[Mesh]
  type = GeneratedMesh
  dim = 1
  nx = 2
  xmax = 1e3
  allow_renumbering = false # for consistent node ids in VPP
[]

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

[Functions]
  [test_fn]
    type = HydrostaticPressureFunction
    reference_point = '0 0 0'
    reference_pressure = 1e5
    reference_temperature = 300
    fluid_properties = fp
    gravity_direction = '1 0 0'
  []
[]

[AuxVariables]
  [p]
  []
[]

[AuxKernels]
  [p_kernel]
    type = FunctionAux
    variable = p
    function = test_fn
    execute_on = 'INITIAL'
  []
[]

[VectorPostprocessors]
  [p_vpp]
    type = NodalValueSampler
    variable = p
    sort_by = x
    execute_on = 'INITIAL'
  []
[]

[Problem]
  solve = false
[]

[Executioner]
  type = Steady
[]

[Outputs]
  csv = true
  execute_on = 'INITIAL'
[]

Input Parameters
Input Files