SimpleFluidProperties

Fluid properties for a simple fluid with a constant bulk density

This is a computationally simple fluid based on a constant bulk modulus density fluid, with density given by ρ=ρ0exp(P/KfαfT), \rho = \rho_{0}\exp(P/K_{f} - \alpha_{f} T), where KfK_{f} (the bulk modulus) and αf\alpha_{f} (the thermal expansion coefficient) are constants.

In this formulation, viscosity and thermal conductivity are constant (with values specified in the input file), while internal energy and enthalpy are given by e=cvT e = c_v T and h=e+γpρ h = e + \gamma\frac{p}{\rho} respectively. Here γ\gamma is called the porepressure_coefficient: usually it should be set to 11 but many analytical studies assume γ=0\gamma=0.

Input Parameters

  • bulk_modulus2e+09Constant bulk modulus (Pa)

    Default:2e+09

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Constant bulk modulus (Pa)

  • cp4194Constant specific heat capacity at constant pressure (J/kg/K)

    Default:4194

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Constant specific heat capacity at constant pressure (J/kg/K)

  • cv4186Constant specific heat capacity at constant volume (J/kg/K)

    Default:4186

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Constant specific heat capacity at constant volume (J/kg/K)

  • density01000Density at zero pressure and zero temperature

    Default:1000

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Density at zero pressure and zero temperature

  • molar_mass0.018Constant molar mass of the fluid (kg/mol)

    Default:0.018

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Constant molar mass of the fluid (kg/mol)

  • porepressure_coefficient1The enthalpy is internal_energy + P / density * porepressure_coefficient. Physically this should be 1.0, but analytic solutions are simplified when it is zero

    Default:1

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:The enthalpy is internal_energy + P / density * porepressure_coefficient. Physically this should be 1.0, but analytic solutions are simplified when it is zero

  • specific_entropy300Constant specific entropy (J/kg/K)

    Default:300

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Constant specific entropy (J/kg/K)

  • thermal_conductivity0.6Constant thermal conductivity (W/m/K)

    Default:0.6

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Constant thermal conductivity (W/m/K)

  • thermal_expansion0.000214Constant coefficient of thermal expansion (1/K)

    Default:0.000214

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Constant coefficient of thermal expansion (1/K)

  • viscosity0.001Constant dynamic viscosity (Pa.s)

    Default:0.001

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Constant dynamic viscosity (Pa.s)

Optional Parameters

  • T_initial_guess400Temperature initial guess for Newton Method variable set conversion

    Default:400

    C++ Type:double

    Unit:(no unit assumed)

    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)

    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)

    Controllable:No

    Description:Tolerance for 2D Newton variable set conversion

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