StiffenedGasFluidProperties

Fluid properties for a stiffened gas

A simple fluid class that implements a stiffened equation of state (Métayer et al., 2004) $var element = document.getElementById("moose-equation-14067132-305a-46fb-a843-a74f11e8ba22");katex.render(" p = (\\gamma - 1) \\rho (e - q) - \\gamma p_{\\infty},", element, {displayMode:true,throwOnError:false,macros:{"\\pf":"\\frac{\\partial #1}{\\partial #2}"}});$ where $var element = document.getElementById("moose-equation-35d1908f-c67f-432f-9b92-8a7f96e29e9c");katex.render("\\gamma = c_p/c_v", element, {displayMode:false,throwOnError:false,macros:{"\\pf":"\\frac{\\partial #1}{\\partial #2}"}});$ is the ratio of specific heat capacities, $var element = document.getElementById("moose-equation-889c4040-730e-415e-b79f-53f33ced9efa");katex.render("q", element, {displayMode:false,throwOnError:false,macros:{"\\pf":"\\frac{\\partial #1}{\\partial #2}"}});$ is a constant that defines the zero reference state for internal energy, and $var element = document.getElementById("moose-equation-74a7f4e6-e425-44d1-b004-ed1a7dda5997");katex.render("p_{\\infty}", element, {displayMode:false,throwOnError:false,macros:{"\\pf":"\\frac{\\partial #1}{\\partial #2}"}});$ is a constant representing the attraction between fluid molecules that makes the fluid stiff in comparison to an ideal gas. This equation of state is typically used to represent water that is under very high pressure.

Input Parameters

cvConstant volume specific heat
C++ Type:double
Options:
Description:Constant volume specific heat
gammaHeat capacity ratio
C++ Type:double
Options:
Description:Heat capacity ratio
p_infStiffness parameter
C++ Type:double
Options:
Description:Stiffness parameter
qParameter defining zero point of internal energy
C++ Type:double
Options:
Description:Parameter defining zero point of internal energy

Required Parameters

M0Molar mass, kg/mol
Default:0
C++ Type:double
Options:
Description:Molar mass, kg/mol
T_c0Critical temperature, K
Default:0
C++ Type:double
Options:
Description:Critical temperature, K
allow_imperfect_jacobiansFalsetrue to allow unimplemented property derivative terms to be set to zero for the AD API
Default:False
C++ Type:bool
Options:
Description:true to allow unimplemented property derivative terms to be set to zero for the AD API
allow_nonphysical_statesTrueAllows for non-physical states, e.g., negative density.
Default:True
C++ Type:bool
Options:
Description:Allows for non-physical states, e.g., negative density.
e_c0Internal energy at the critical point, J/kg
Default:0
C++ Type:double
Options:
Description:Internal energy at the critical point, J/kg
emit_on_nannoneRaise mooseWarning or mooseError?
Default:none
C++ Type:MooseEnum
Options:none warning error
Description:Raise mooseWarning or mooseError?
execute_onTIMESTEP_ENDThe list of flag(s) indicating when this object should be executed, the available options include NONE, INITIAL, LINEAR, NONLINEAR, TIMESTEP_END, TIMESTEP_BEGIN, FINAL, CUSTOM.
Default:TIMESTEP_END
C++ Type:ExecFlagEnum
Options:NONE INITIAL LINEAR NONLINEAR TIMESTEP_END TIMESTEP_BEGIN FINAL CUSTOM
Description:The list of flag(s) indicating when this object should be executed, the available options include NONE, INITIAL, LINEAR, NONLINEAR, TIMESTEP_END, TIMESTEP_BEGIN, FINAL, CUSTOM.
fp_typesingle-phase-fpType of the fluid property object
Default:single-phase-fp
C++ Type:FPType
Options:
Description:Type of the fluid property object
k0.6Thermal conductivity, W/(m-K)
Default:0.6
C++ Type:double
Options:
Description:Thermal conductivity, W/(m-K)
mu0.001Dynamic viscosity, Pa.s
Default:0.001
C++ Type:double
Options:
Description:Dynamic viscosity, Pa.s
q_prime0Parameter
Default:0
C++ Type:double
Options:
Description:Parameter
rho_c0Critical density, kg/m3
Default:0
C++ Type:double
Options:
Description:Critical density, kg/m3

Optional 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
Options:
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).
control_tagsAdds user-defined labels for accessing object parameters via control logic.
C++ Type:std::vector
Options:
Description:Adds user-defined labels for accessing object parameters via control logic.
enableTrueSet the enabled status of the MooseObject.
Default:True
C++ Type:bool
Options:
Description:Set the enabled status of the MooseObject.
force_preauxFalseForces the GeneralUserObject to be executed in PREAUX
Default:False
C++ Type:bool
Options:
Description:Forces the GeneralUserObject to be executed in PREAUX
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
Options:
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/fluid_properties/test/tests/auxkernels/stagnation_temperature_aux.i
modules/fluid_properties/test/tests/auxkernels/stagnation_pressure_aux.i
modules/fluid_properties/test/tests/auxkernels/specific_enthalpy_aux.i
modules/fluid_properties/test/tests/stiffened_gas/test.i

modules/fluid_properties/test/tests/auxkernels/stagnation_temperature_aux.i

[Mesh]
  type = GeneratedMesh
  dim = 1
  nx = 2
[]

[Variables]
  [./u]
    order = FIRST
    family = LAGRANGE
  [../]
[]

[AuxVariables]
  [./specific_internal_energy]
  [../]
  [./specific_volume]
  [../]
  [./velocity]
  [../]
  [./stagnation_temperature]
  [../]
[]

[Kernels]
  [./diff_u]
    type = Diffusion
    variable = u
  [../]
[]

[AuxKernels]
  [./specific_internal_energy_ak]
    type = ConstantAux
    variable = specific_internal_energy
    value = 1026.2e3
  [../]
  [./specific_volume_ak]
    type = ConstantAux
    variable = specific_volume
    value = 0.0012192
  [../]
  [./velocity_ak]
    type = ConstantAux
    variable = velocity
    value = 10.0
  [../]
  [./stagnation_temperature_ak]
    type = StagnationTemperatureAux
    variable = stagnation_temperature
    e = specific_internal_energy
    v = specific_volume
    vel = velocity
    fp = eos
  [../]
[]

[Modules]
  [./FluidProperties]
    [./eos]
      type = StiffenedGasFluidProperties
      gamma = 2.35
      q = -1167e3
      q_prime = 0.0
      p_inf = 1e9
      cv = 1816.0
    [../]
  []
[]

[BCs]
  [./left_u]
    type = DirichletBC
    variable = u
    boundary = 0
    value = 1
  [../]

  [./right_u]
    type = DirichletBC
    variable = u
    boundary = 1
    value = 2
  [../]
[]

[Executioner]
  type = Steady

  solve_type = 'PJFNK'
[]

[Outputs]
  exodus = true
[]

modules/fluid_properties/test/tests/auxkernels/stagnation_pressure_aux.i

[Mesh]
  type = GeneratedMesh
  dim = 1
  nx = 2
[]

[Variables]
  [./u]
    order = FIRST
    family = LAGRANGE
  [../]
[]

[AuxVariables]
  [./specific_internal_energy]
  [../]
  [./specific_volume]
  [../]
  [./velocity]
  [../]
  [./stagnation_pressure]
  [../]
[]

[Kernels]
  [./diff_u]
    type = Diffusion
    variable = u
  [../]
[]

[AuxKernels]
  [./specific_internal_energy_ak]
    type = ConstantAux
    variable = specific_internal_energy
    value = 1026.2e3
  [../]
  [./specific_volume_ak]
    type = ConstantAux
    variable = specific_volume
    value = 0.0012192
  [../]
  [./velocity_ak]
    type = ConstantAux
    variable = velocity
    value = 10.0
  [../]
  [./stagnation_pressure_ak]
    type = StagnationPressureAux
    variable = stagnation_pressure
    e = specific_internal_energy
    v = specific_volume
    vel = velocity
    fp = eos
  [../]
[]

[Modules]
  [./FluidProperties]
    [./eos]
      type = StiffenedGasFluidProperties
      gamma = 2.35
      q = -1167e3
      q_prime = 0.0
      p_inf = 1e9
      cv = 1816.0
    [../]
  []
[]

[BCs]
  [./left_u]
    type = DirichletBC
    variable = u
    boundary = 0
    value = 1
  [../]

  [./right_u]
    type = DirichletBC
    variable = u
    boundary = 1
    value = 2
  [../]
[]

[Executioner]
  type = Steady

  solve_type = 'PJFNK'
[]

[Outputs]
  exodus = true
[]

modules/fluid_properties/test/tests/auxkernels/specific_enthalpy_aux.i

[Mesh]
  type = GeneratedMesh
  dim = 1
  nx = 2
[]

[Variables]
  [./u]
    order = FIRST
    family = LAGRANGE
  [../]
[]

[AuxVariables]
  [./pressure]
  [../]
  [./temperature]
  [../]
  [./specific_enthalpy]
  [../]
[]

[Kernels]
  [./diff_u]
    type = Diffusion
    variable = u
  [../]
[]

[AuxKernels]
  [./pressure_ak]
    type = ConstantAux
    variable = pressure
    value = 10e6
  [../]
  [./temperature_ak]
    type = ConstantAux
    variable = temperature
    value = 400.0
  [../]
  [./specific_enthalpy_ak]
    type = SpecificEnthalpyAux
    variable = specific_enthalpy
    fp = eos
    p = pressure
    T = temperature
  [../]
[]

[Modules]
  [./FluidProperties]
    [./eos]
      type = StiffenedGasFluidProperties
      gamma = 2.35
      q = -1167e3
      q_prime = 0.0
      p_inf = 1e9
      cv = 1816.0
    [../]
  []
[]

[BCs]
  [./left_u]
    type = DirichletBC
    variable = u
    boundary = 0
    value = 1
  [../]

  [./right_u]
    type = DirichletBC
    variable = u
    boundary = 1
    value = 2
  [../]
[]

[Executioner]
  type = Steady

  solve_type = 'PJFNK'
[]

[Outputs]
  exodus = true
[]

modules/fluid_properties/test/tests/stiffened_gas/test.i

[Mesh]
  type = GeneratedMesh
  dim = 2
  nx = 2
  ny = 2
  elem_type = QUAD4
[]

[Functions]
  [./f_fn]
    type = ParsedFunction
    value = -4
  [../]
  [./bc_fn]
    type = ParsedFunction
    value = 'x*x+y*y'
  [../]
[]

[Variables]
  [./u]
  [../]
[]

[AuxVariables]
  [./e]
    initial_condition = 113206.45935406466
  [../]
  [./v]
    initial_condition = 0.0007354064593540647
  [../]

  [./p]
    family = MONOMIAL
    order = CONSTANT
  [../]
  [./T]
    family = MONOMIAL
    order = CONSTANT
  [../]
  [./cp]
    family = MONOMIAL
    order = CONSTANT
  [../]
  [./cv]
    family = MONOMIAL
    order = CONSTANT
  [../]
  [./c]
    family = MONOMIAL
    order = CONSTANT
  [../]
  [./mu]
    family = MONOMIAL
    order = CONSTANT
  [../]
  [./k]
    family = MONOMIAL
    order = CONSTANT
  [../]
  [./g]
    family = MONOMIAL
    order = CONSTANT
  [../]
[]

[AuxKernels]
  [./p]
    type = MaterialRealAux
     variable = p
     property = pressure
  [../]
  [./T]
    type = MaterialRealAux
     variable = T
     property = temperature
  [../]
  [./cp]
    type = MaterialRealAux
     variable = cp
     property = cp
  [../]
  [./cv]
    type = MaterialRealAux
     variable = cv
     property = cv
  [../]
  [./c]
    type = MaterialRealAux
     variable = c
     property = c
  [../]
  [./mu]
    type = MaterialRealAux
     variable = mu
     property = mu
  [../]
  [./k]
    type = MaterialRealAux
     variable = k
     property = k
  [../]
  [./g]
    type = MaterialRealAux
     variable = g
     property = g
  [../]
[]

[Modules]
  [./FluidProperties]
    [./sg]
      type = StiffenedGasFluidProperties
      gamma = 2.35
      q = -1167e3
      q_prime = 0
      p_inf = 1.e9
      cv = 1816

      mu = 0.9
      k = 0.6
    [../]
  []
[]

[Materials]
  [./fp_mat]
    type = FluidPropertiesMaterial
    e = e
    v = v
    fp = sg
  [../]
[]

[Kernels]
  [./diff]
    type = Diffusion
    variable = u
  [../]
  [./ffn]
    type = BodyForce
    variable = u
    function = f_fn
  [../]
[]

[BCs]
  [./all]
    type = FunctionDirichletBC
    variable = u
    boundary = 'left right top bottom'
    function = bc_fn
  [../]
[]

[Executioner]
  type = Steady
  solve_type = NEWTON
[]

[Outputs]
  exodus = true
[]

References

O. L. Métayer, J. Massoni, and R. Saurel. Elaborating equations of state of a liquid and its vapor for two-phase flow models. Int. J. Therm. Sci., 43:265–276, 2004.

@article{metayer2004,
    author = "M{\'e}tayer, O. L. and Massoni, J. and Saurel, R.",
    title = "Elaborating equations of state of a liquid and its vapor for two-phase flow models",
    journal = "Int. J. Therm. Sci.",
    volume = "43",
    pages = "265--276",
    year = "2004"
}

Input Parameters
Input Files
References

Application Usage

Physics and Syntax

Application Development

Framework Development

MOOSEDocs

Infrastructure

Questions

Information and Tools