TemperaturePressureFunctionFluidProperties

Single-phase fluid properties that allows to provide thermal conductivity, density, and viscosity as functions of temperature and pressure.

The temperature and pressure will be passed as respectively the x and y spatial arguments of the functions. The following relations hold true:

$var element = document.getElementById("moose-equation-b3d4c187-8323-46e6-bc73-0bcc9e761dbb");katex.render("\\begin{aligned} x = T \\\\ y = P \\\\ \\rho = \\rho^{user}(time=0, (x=T, y=P, z=0)) \\\\ \\mu = \\mu^{user}(time=0, (x=T, y=P, z=0)) \\\\ k = k^{user}(time=0, (x=T, y=P, z=0)) \\\\ c_v = c_v^{user} \\\\ e = c_v * T \\end{aligned}", element, {displayMode:true,throwOnError:false});$

with $var element = document.getElementById("moose-equation-1d4ddf3c-41a4-4a69-9cba-33c7113cc311");katex.render("T", element, {displayMode:false,throwOnError:false});$ the temperature, $var element = document.getElementById("moose-equation-cbc58d5c-ed28-468e-aa87-e29345b085a2");katex.render("P", element, {displayMode:false,throwOnError:false});$ the pressure, $var element = document.getElementById("moose-equation-2f7198b9-c0fd-45c8-8a70-fb3099ce4e28");katex.render("\\rho", element, {displayMode:false,throwOnError:false});$ the density, $var element = document.getElementById("moose-equation-ae0bbe0f-84bc-415e-a140-c9ec61e3350b");katex.render("\\mu", element, {displayMode:false,throwOnError:false});$ the dynamic viscosity, $var element = document.getElementById("moose-equation-4c70d93f-ffa3-4cae-a160-573daca35d0d");katex.render("k", element, {displayMode:false,throwOnError:false});$ the thermal conductivity, $var element = document.getElementById("moose-equation-c6a6f55f-504e-40fe-9999-0130ac2d5074");katex.render("c_v", element, {displayMode:false,throwOnError:false});$ the specific isochoric heat capacity, $var element = document.getElementById("moose-equation-6e68da58-78e7-440d-bcd2-7143c5ec4161");katex.render("e", element, {displayMode:false,throwOnError:false});$ the specific energy and the $var element = document.getElementById("moose-equation-5b27b3b0-769a-4bb4-b342-dd322e7d4626");katex.render("user", element, {displayMode:false,throwOnError:false});$ exponent indicating a user-passed parameter.

The derivatives of the fluid properties are obtained using the Function(s) gradient components and the appropriate derivative chaining for derived properties.

warning

The range of validity of the property is based off of the validity of the functions that are input, and is not checked by this FluidProperties object.

note

Support for the conservative (specific volume, internal energy) variable set is only partial. Notable missing implementations are routines for entropy, the speed of sound, and some conversions between specific enthalpy and specific energy.

warning

Due to the approximations made when computing the isobaric heat capacity from the constant isochoric heat capacity, this material should only be used for nearly-incompressible fluids.

Example Input File Syntax

In this example, temperature and pressure dependent density, dynamic viscosity and thermal conductivity of a fluid are being set using three ParsedFunctions. The functions are specified with the x and y coordinates, and are evaluated with respectively the temperature and pressure variables.

[Functions]
  # This demonstrates how to define fluid properties that are functions
  # of the LOCAL value of the (p,T) variables
  # x for temperature
  # y for pressure
  [k]
    type = ParsedFunction
    expression = '14 + 1e-2 * x + 1e-5 * y'
  []
  [rho]
    type = ParsedFunction
    expression = '1.5e3 + 0.13 * x - 1.5e-4 * y'
  []
  [mu]
    type = ParsedFunction
    expression = '1e-3 + 2e-6 * x - 3e-9 * y'
  []
[]

[FluidProperties]
  [fp]
    type = TemperaturePressureFunctionFluidProperties
    cv = ${cv}
    k = k
    rho = rho
    mu = mu
  []
[]

Input Parameters

cvConstant isochoric specific heat [J/(kg-K)]
C++ Type:double
Controllable:No
Description:Constant isochoric specific heat [J/(kg-K)]
kThermal conductivity function of temperature and pressure [W/(m-K)]
C++ Type:FunctionName
Controllable:No
Description:Thermal conductivity function of temperature and pressure [W/(m-K)]
muDynamic viscosity function of temperature and pressure [Pa-s]
C++ Type:FunctionName
Controllable:No
Description:Dynamic viscosity function of temperature and pressure [Pa-s]
rhoDensity function of temperature and pressure [kg/m^3]
C++ Type:FunctionName
Controllable:No
Description:Density function of temperature and pressure [kg/m^3]

Required Parameters

execute_onTIMESTEP_ENDThe list of flag(s) indicating when this object should be executed, the available options include NONE, INITIAL, LINEAR, NONLINEAR, POSTCHECK, TIMESTEP_END, TIMESTEP_BEGIN, MULTIAPP_FIXED_POINT_END, MULTIAPP_FIXED_POINT_BEGIN, FINAL, CUSTOM.
Default:TIMESTEP_END
C++ Type:ExecFlagEnum
Options:NONE, INITIAL, LINEAR, NONLINEAR, POSTCHECK, TIMESTEP_END, TIMESTEP_BEGIN, MULTIAPP_FIXED_POINT_END, MULTIAPP_FIXED_POINT_BEGIN, FINAL, CUSTOM
Controllable:No
Description:The list of flag(s) indicating when this object should be executed, the available options include NONE, INITIAL, LINEAR, NONLINEAR, POSTCHECK, TIMESTEP_END, TIMESTEP_BEGIN, MULTIAPP_FIXED_POINT_END, MULTIAPP_FIXED_POINT_BEGIN, FINAL, CUSTOM.
prop_getter_suffixAn optional suffix parameter that can be appended to any attempt to retrieve/get material properties. The suffix will be prepended with a '_' character.
C++ Type:MaterialPropertyName
Controllable:No
Description:An optional suffix parameter that can be appended to any attempt to retrieve/get material properties. The suffix will be prepended with a '_' character.
use_interpolated_stateFalseFor the old and older state use projected material properties interpolated at the quadrature points. To set up projection use the ProjectedStatefulMaterialStorageAction.
Default:False
C++ Type:bool
Controllable:No
Description:For the old and older state use projected material properties interpolated at the quadrature points. To set up projection use the ProjectedStatefulMaterialStorageAction.

Optional Parameters

T_initial_guess400Temperature initial guess for Newton Method variable set conversion
Default:400
C++ Type:double
Controllable:No
Description:Temperature initial guess for Newton Method variable set conversion
p_initial_guess200000Pressure initial guess for Newton Method variable set conversion
Default:200000
C++ Type:double
Controllable:No
Description:Pressure initial guess for Newton Method variable set conversion
tolerance1e-08Tolerance for 2D Newton variable set conversion
Default:1e-08
C++ Type:double
Controllable:No
Description:Tolerance for 2D Newton variable set conversion

Variable Set Conversions Newton Solve Parameters

allow_duplicate_execution_on_initialFalseIn the case where this UserObject is depended upon by an initial condition, allow it to be executed twice during the initial setup (once before the IC and again after mesh adaptivity (if applicable).
Default:False
C++ Type:bool
Controllable:No
Description:In the case where this UserObject is depended upon by an initial condition, allow it to be executed twice during the initial setup (once before the IC and again after mesh adaptivity (if applicable).
allow_imperfect_jacobiansFalsetrue to allow unimplemented property derivative terms to be set to zero for the AD API
Default:False
C++ Type:bool
Controllable:No
Description:true to allow unimplemented property derivative terms to be set to zero for the AD API
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.
execution_order_group0Execution order groups are executed in increasing order (e.g., the lowest number is executed first). Note that negative group numbers may be used to execute groups before the default (0) group. Please refer to the user object documentation for ordering of user object execution within a group.
Default:0
C++ Type:int
Controllable:No
Description:Execution order groups are executed in increasing order (e.g., the lowest number is executed first). Note that negative group numbers may be used to execute groups before the default (0) group. Please refer to the user object documentation for ordering of user object execution within a group.
force_postauxFalseForces the UserObject to be executed in POSTAUX
Default:False
C++ Type:bool
Controllable:No
Description:Forces the UserObject to be executed in POSTAUX
force_preauxFalseForces the UserObject to be executed in PREAUX
Default:False
C++ Type:bool
Controllable:No
Description:Forces the UserObject to be executed in PREAUX
force_preicFalseForces the UserObject to be executed in PREIC during initial setup
Default:False
C++ Type:bool
Controllable:No
Description:Forces the UserObject to be executed in PREIC during initial setup
fp_typesingle-phase-fpType of the fluid property object
Default:single-phase-fp
C++ Type:FPType
Controllable:No
Description:Type of the fluid property object
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

(moose/modules/fluid_properties/test/tests/temperature_pressure_function/example.i)

# Test implementation of TemperaturePressureFunctionFluidProperties properties by comparison to analytical functions.

cv = 4000
T_initial = 400
p_initial = 1e5

[Mesh]
  type = GeneratedMesh
  dim = 1
[]

[Problem]
  solve = false
[]

[AuxVariables]
  [temperature]
    initial_condition = ${T_initial}
  []
  [pressure]
    initial_condition = 1e5
  []
[]

[Functions]
  # This demonstrates how to define fluid properties that are functions
  # of the LOCAL value of the (p,T) variables
  # x for temperature
  # y for pressure
  [k]
    type = ParsedFunction
    expression = '14 + 1e-2 * x + 1e-5 * y'
  []
  [rho]
    type = ParsedFunction
    expression = '1.5e3 + 0.13 * x - 1.5e-4 * y'
  []
  [mu]
    type = ParsedFunction
    expression = '1e-3 + 2e-6 * x - 3e-9 * y'
  []
[]

[FluidProperties]
  [fp]
    type = TemperaturePressureFunctionFluidProperties
    cv = ${cv}
    k = k
    rho = rho
    mu = mu
  []
[]

[Materials]
  [to_vars]
    type = FluidPropertiesMaterialPT
    fp = fp
    outputs = 'all'
    output_properties = 'density k cp cv viscosity e h'
    pressure = pressure
    temperature = temperature

    compute_entropy = false
    compute_sound_speed = false
  []
[]

[Executioner]
  type = Steady
[]

[Postprocessors]
  [k_exact]
    type = FunctionValuePostprocessor
    function = k
    outputs = none
    point = '${T_initial} ${p_initial} 0'
  []
  [rho_exact]
    type = FunctionValuePostprocessor
    function = rho
    outputs = none
    point = '${T_initial} ${p_initial} 0'
  []
  [mu_exact]
    type = FunctionValuePostprocessor
    function = mu
    outputs = none
    point = '${T_initial} ${p_initial} 0'
  []
  [e_exact]
    type = Receiver
    default = '${fparse cv * T_initial}'
    outputs = none
  []
  [cv_exact]
    type = Receiver
    default = '${fparse cv}'
    outputs = none
  []

  # Postprocessors to get from the fluid property object
  [k_avg]
    type = ElementAverageValue
    variable = k
    outputs = none
  []
  [rho_avg]
    type = ElementAverageValue
    variable = density
    outputs = none
  []
  [mu_avg]
    type = ElementAverageValue
    variable = viscosity
    outputs = none
  []
  [cv_avg]
    type = ElementAverageValue
    variable = cv
    outputs = none
  []
  [e_avg]
    type = ElementAverageValue
    variable = e
    outputs = none
  []

  # We output these directly, cant compare to anything analytical though
  [cp_avg]
    type = ElementAverageValue
    variable = cp
  []
  [h_avg]
    type = ElementAverageValue
    variable = h
  []

  # Postprocessors to compare the two
  [k_diff]
    type = DifferencePostprocessor
    value1 = k_exact
    value2 = k_avg
  []
  [mu_diff]
    type = DifferencePostprocessor
    value1 = mu_exact
    value2 = mu_avg
  []
  [rho_diff]
    type = DifferencePostprocessor
    value1 = rho_exact
    value2 = rho_avg
  []
  [e_diff]
    type = DifferencePostprocessor
    value1 = e_exact
    value2 = e_avg
  []
  [cv_diff]
    type = DifferencePostprocessor
    value1 = cv_exact
    value2 = cv_avg
  []
[]

[Outputs]
  # Note that diffs wont be settled until timestep 2 because of order of execution
  csv = true
[]

Example Input File Syntax
Input Parameters