Sodium Fluid Properties

Fluid properties for sodium

SodiumProperties inherits from the base class FluidProperties and provides the needed functionality for determining sodium fluid properties that may be used for coolant channel simulations.

The following methods are provided in this class:

Thermal Conductivity: k(temperature)
Enthalpy: h(temperature)
Heat capacity: heatCapacity(temperature)
Temperature: temperature(h)
Density: rho(temperature)
Derivative of Density w.r.t temperature: drho_dT(temperature)
Derivative of Density w.r.t enthalpy: drho_dh(h)

Properties of liquid sodium are obtained from Fink and Leibowitz (1995).

[FluidProperties<<<{"href": "../../syntax/FluidProperties/index.html"}>>>]
  [sodium]
    type = SodiumProperties<<<{"description": "Fluid properties for sodium", "href": "SodiumProperties.html"}>>>
  []
[]

Input Parameters

specific_heatOptional value for specific heat that overrides interal calculations
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Optional value for specific heat that overrides interal calculations
thermal_conductivityOptional value for thermal conductivity that overrides interal calculations
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Optional value for thermal conductivity that overrides interal calculations

Optional 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_typesodium-specific-fpType of the fluid property object
Default:sodium-specific-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

References

J K Fink and L Leibowitz. Thermodynamic and transport properties of sodium liquid and vapor. Technical Report, Argonne National Laboratory (ANL), Argonne, IL, January 1995.

@techreport{Fink:1995bf,
    author = "Fink, J K and Leibowitz, L",
    title = "{Thermodynamic and transport properties of sodium liquid and vapor}",
    institution = "Argonne National Laboratory (ANL), Argonne, IL",
    year = "1995",
    month = "January"
}

(moose/modules/fluid_properties/test/tests/sodium/exact.i)

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

[Mesh]
  type = GeneratedMesh
  dim = 1
[]

[Problem]
  solve = false
[]

[AuxVariables]
  [./temperature]
  [../]
[]

[AuxKernels]
  [./temperature_aux]
    type = FunctionAux
    variable = temperature
    function = '400 + 200 * t'
  [../]
[]

[Functions]
  [./k]
    type = ParsedFunction
    symbol_names = 'T'
    symbol_values = 'temperature'
    expression = '124.67 - 0.11381 * T + 5.5226e-5 * T^2 - 1.1842e-8 * T^3'
  [../]
  [./h]
    type = ParsedFunction
    symbol_names = 'T'
    symbol_values = 'temperature'
    expression = '1.0e3 * (-365.77 + 1.6582 * T - 4.2395e-4 * T^2 + 1.4847e-7 * T^3 + 2992.6 / T)'
  [../]
  [./cp]
    type = ParsedFunction
    symbol_names = 'T'
    symbol_values = 'temperature'
    expression = '1.0e3 * (1.6582 - 8.4790e-4 * T + 4.4541e-7 * T^2 - 2992.6 / T^2)'
  [../]
  [./rho]
    type = ParsedFunction
    symbol_names = 'T'
    symbol_values = 'temperature'
    expression = '219.0 + 275.32 * (1.0 - T / 2503.7) + 511.58 * (1.0 - T / 2503.7)^(0.5)'
  [../]
  [./drho_dT]
    type = ParsedFunction
    symbol_names = 'T'
    symbol_values = 'temperature'
    expression = '-(2.0 * 275.32 + 511.58 / (1.0 - T / 2503.7)^(0.5)) / 2.0 / 2503.7'
  [../]
  [./drho_dh]
    type = ParsedFunction
    symbol_names = 'drho_dT_exact cp_exact'
    symbol_values = 'drho_dT_exact cp_exact'
    expression = 'drho_dT_exact/cp_exact'
  [../]
[]

[FluidProperties/sodium]
  type = SodiumProperties
[]

[Materials]
  [./fp_mat]
    type = SodiumPropertiesMaterial
    temperature = temperature
    outputs = all
  [../]
[]


[Executioner]
  type = Transient

  num_steps = 10
[]

[Postprocessors]
  [./temperature]
    type = ElementAverageValue
    variable = temperature
    outputs = none
  [../]
  [./k_exact]
    type = FunctionValuePostprocessor
    function = k
    outputs = none
  [../]
  [./h_exact]
    type = FunctionValuePostprocessor
    function = h
    outputs = none
  [../]
  [./cp_exact]
    type = FunctionValuePostprocessor
    function = cp
    outputs = none
  [../]
  [./rho_exact]
    type = FunctionValuePostprocessor
    function = rho
    outputs = none
  [../]
  [./drho_dT_exact]
    type = FunctionValuePostprocessor
    function = drho_dT
    outputs = none
  [../]
  [./drho_dh_exact]
    type = FunctionValuePostprocessor
    function = drho_dh
    outputs = none
  [../]
  [./k_avg]
    type = ElementAverageValue
    variable = k
    outputs = none
  [../]
  [./h_avg]
    type = ElementAverageValue
    variable = h
    outputs = none
  [../]
  [./cp_avg]
    type = ElementAverageValue
    variable = cp
    outputs = none
  [../]
  [./t_from_h_avg]
    type = ElementAverageValue
    variable = temperature
    outputs = none
  [../]
  [./rho_avg]
    type = ElementAverageValue
    variable = rho
    outputs = none
  [../]
  [./drho_dT_avg]
    type = ElementAverageValue
    variable = drho_dT
    outputs = none
  [../]
  [./drho_dh_avg]
    type = ElementAverageValue
    variable = drho_dh
    outputs = none
  [../]
  [./k_diff]
    type = DifferencePostprocessor
    value1 = k_exact
    value2 = k_avg
  [../]
  [./h_diff]
    type = DifferencePostprocessor
    value1 = h_exact
    value2 = h_avg
  [../]
  [./cp_diff]
    type = DifferencePostprocessor
    value1 = cp_exact
    value2 = cp_avg
  [../]
  [./t_from_h_diff]
    type = DifferencePostprocessor
    value1 = temperature
    value2 = t_from_h_avg
  [../]
  [./rho_avg_diff]
    type = DifferencePostprocessor
    value1 = rho_exact
    value2 = rho_avg
  [../]
  [./drho_dT_avg_diff]
    type = DifferencePostprocessor
    value1 = drho_dT_exact
    value2 = drho_dT_avg
  [../]
  [./drho_dh_avg_diff]
    type = DifferencePostprocessor
    value1 = drho_dh_exact
    value2 = drho_dh_avg
  [../]
[]

[Outputs]
  csv = true
[]

Input Parameters
References