LinearFluidProperties

Fluid properties for a fluid with density linearly dependent on temperature and pressure

All other fluid properties are either constant, such as:

dynamic viscosity
thermal conductivity
specific heat
speed of sound
Prandtl number

or dependent on pressure and temperature through density's dependence:

specific energy
specific enthalpy

Input Parameters

T_0Reference internal energy
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Reference internal energy
a2dp/d(rho)
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:dp/d(rho)
betaCoefficient of thermal expansion
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Coefficient of thermal expansion
cvSpecific heat
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Specific heat
e_0Reference internal energy
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Reference internal energy
kThermal conductivity, W/(m-K)
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Thermal conductivity, W/(m-K)
muDynamic viscosity, Pa.s
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Dynamic viscosity, Pa.s
p_0Reference pressure
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Reference pressure
rho_0Reference density
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Reference density

Required Parameters

T_initial_guess400Temperature initial guess for Newton Method variable set conversion
Default:400
C++ Type:double
Unit:(no unit assumed)
Range:T_initial_guess > 0
Controllable:No
Description:Temperature initial guess for Newton Method variable set conversion
max_newton_its100Maximum number of Newton iterations for variable set conversions
Default:100
C++ Type:unsigned int
Controllable:No
Description:Maximum number of Newton iterations for variable set conversions
p_initial_guess200000Pressure initial guess for Newton Method variable set conversion
Default:200000
C++ Type:double
Unit:(no unit assumed)
Range:p_initial_guess > 0
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
Unit:(no unit assumed)
Range:tolerance > 0
Controllable:No
Description:Tolerance for 2D Newton variable set conversion
verbose_newtonFalseWhether to output Newton inversion iterations to console
Default:False
C++ Type:bool
Controllable:No
Description:Whether to output Newton inversion iterations to console

Variable Set Conversions Newton Solve Parameters

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.
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

Advanced Parameters

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
Unit:(no unit assumed)
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.

Material Property Retrieval Parameters

Input Files

(modules/thermal_hydraulics/test/tests/components/flow_channel_1phase/err.wrong_end.i)

(modules/thermal_hydraulics/test/tests/components/flow_channel_1phase/err.wrong_end.i)

[GlobalParams]
  gravity_vector = '0 0 0'
  scaling_factor_1phase = '1. 1e-4'

  closures = simple_closures
[]

[FluidProperties]
  [barotropic]
    type = LinearFluidProperties
    p_0 = 1.e5     # Pa
    rho_0 = 1.e3   # kg/m^3
    a2 = 1.e7      # m^2/s^2
    beta = .46e-3 # K^{-1}
    cv = 4.18e3    # J/kg-K, could be a global parameter?
    e_0 = 1.254e6  # J/kg
    T_0 = 300      # K
  []
[]

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

[Components]
  [pipe]
    type = FlowChannel1Phase
    fp = barotropic
    # geometry
    position = '0 0 0'
    orientation = '1 0 0'
    A = 1.
    D_h = 1.12837916709551
    f = 0.01
    length = 1
    n_elems = 100
  []

  [inlet]
    type = InletStagnationPressureTemperature1Phase
    input = 'pipe:asdf'      # this is an error we are checking for
    p0 = 1e5
    T0 = 300
  []

  [outlet]
    type = Outlet1Phase
    input = 'pipe:out'
    p = 9.5e4
  []
[]

[Preconditioning]
  [FDP_PJFNK]
    type = FDP
    full = true

    petsc_options_iname = '-mat_fd_coloring_err'
    petsc_options_value = '1.e-10'
  []
[]


[Executioner]
  type = Transient
  scheme = 'bdf2'

  dt = 1e-4
  dtmin = 1.e-7

  solve_type = 'PJFNK'
  nl_rel_tol = 1e-9
  nl_abs_tol = 1e-8
  nl_max_its = 10

  l_tol = 1e-8
  l_max_its = 100

  start_time = 0.0
  num_steps = 10
[]


[Outputs]
  [out]
    type = Exodus
  []
[]

Input Parameters
Input Files