www.mooseframework.org
Public Member Functions | List of all members
SodiumProperties Class Reference

Properties of liquid sodium from ANL/RE-95/2 report "Thermodynamic and Transport Properties of Sodium Liquid and Vapor" from ANL Reactor Engineering Division. More...

#include <SodiumProperties.h>

Inheritance diagram for SodiumProperties:
[legend]

Public Member Functions

 SodiumProperties (const InputParameters &parameters)
 
Real k (Real temperature) const
 Thermal conductivity as a function of temperature. More...
 
Real h (Real temperature) const
 Enthalpy of liquid Na (relative to solid Na at STP) in J/kg as a function of temperature. More...
 
Real heatCapacity (Real temperature) const
 Heat capacity of liquid Na in J/kg-K as a function of temperature. More...
 
Real temperature (Real enthalpy) const
 Inverse solve for temperature [K] from enthalpy [J/kg]. More...
 
Real rho (Real temperature) const
 Density as a function of temperature. More...
 
Real drho_dT (Real temperature) const
 Derivative of density w.r.t temperature. More...
 
Real drho_dh (Real enthalpy) const
 Derivative of density w.r.t enthalpy. More...
 
virtual void execute () final
 
virtual void initialize () final
 
virtual void finalize () final
 

Detailed Description

Properties of liquid sodium from ANL/RE-95/2 report "Thermodynamic and Transport Properties of Sodium Liquid and Vapor" from ANL Reactor Engineering Division.

Definition at line 23 of file SodiumProperties.h.

Constructor & Destructor Documentation

SodiumProperties::SodiumProperties ( const InputParameters &  parameters)

Definition at line 20 of file SodiumProperties.C.

20  : FluidProperties(parameters)
21 {
22 }
FluidProperties(const InputParameters &parameters)

Member Function Documentation

Real SodiumProperties::drho_dh ( Real  enthalpy) const

Derivative of density w.r.t enthalpy.

Parameters
henthalpy

Definition at line 95 of file SodiumProperties.C.

96 {
97  Real T = 0.0;
98  T = temperature(h);
99  return drho_dT(T) / heatCapacity(T);
100 }
Real heatCapacity(Real temperature) const
Heat capacity of liquid Na in J/kg-K as a function of temperature.
Real drho_dT(Real temperature) const
Derivative of density w.r.t temperature.
Real temperature(Real enthalpy) const
Inverse solve for temperature [K] from enthalpy [J/kg].
Real h(Real temperature) const
Enthalpy of liquid Na (relative to solid Na at STP) in J/kg as a function of temperature.
Real SodiumProperties::drho_dT ( Real  temperature) const

Derivative of density w.r.t temperature.

Parameters
Ttemperature

Definition at line 84 of file SodiumProperties.C.

Referenced by drho_dh().

85 {
86  const Real f = 275.32;
87  const Real g = 511.58;
88  const Real Tc = 2503.7; // critical temperature, K
89  mooseAssert(T < Tc, "Temperature is greater than critical temperature 2503.7 K ");
90 
91  return -(f + g * (0.5) / std::sqrt(1 - T / Tc)) / Tc;
92 }
virtual void FluidProperties::execute ( )
inlinefinalvirtualinherited

Definition at line 25 of file FluidProperties.h.

25 {}
virtual void FluidProperties::finalize ( )
inlinefinalvirtualinherited

Definition at line 27 of file FluidProperties.h.

27 {}
Real SodiumProperties::h ( Real  temperature) const

Enthalpy of liquid Na (relative to solid Na at STP) in J/kg as a function of temperature.

Parameters
Ttemperature From page 4. Valid from 371 to 2000 K, and relative to enthalpy of solid Na at 298.15.

Definition at line 33 of file SodiumProperties.C.

Referenced by temperature().

34 {
35  const Real T2 = T * T;
36  const Real T3 = T2 * T;
37 
38  // Converted from kJ/kg to J/kg.
39  return -365.77e3 + 1.6582e3 * T - 4.2395e-1 * T2 + 1.4847e-4 * T3 + 2992.6e3 / T;
40 }
Real SodiumProperties::heatCapacity ( Real  temperature) const

Heat capacity of liquid Na in J/kg-K as a function of temperature.

Parameters
Ttemperature From page 29 (or by differentiating enthalpy). Valid from 371 to 2000 K.

Definition at line 43 of file SodiumProperties.C.

Referenced by drho_dh(), and temperature().

44 {
45  const Real T2 = T * T;
46  // Converted from kJ/kg-K to J/kg-K.
47  return 1.6582e3 - 8.4790e-1 * T + 4.4541e-4 * T2 - 2992.6e3 / T2;
48 }
virtual void FluidProperties::initialize ( )
inlinefinalvirtualinherited

Definition at line 26 of file FluidProperties.h.

26 {}
Real SodiumProperties::k ( Real  temperature) const

Thermal conductivity as a function of temperature.

Parameters
Ttemperature From page 181. Valid from 371 to 1500 K.

Definition at line 25 of file SodiumProperties.C.

26 {
27  const Real T2 = T * T;
28  const Real T3 = T2 * T;
29  return 124.67 - 0.11381 * T + 5.5226e-5 * T2 - 1.1842e-8 * T3;
30 }
Real SodiumProperties::rho ( Real  temperature) const

Density as a function of temperature.

Parameters
Ttemperature

Definition at line 72 of file SodiumProperties.C.

73 {
74  const Real rhoc = 219.0; // kg/m^3
75  const Real f = 275.32;
76  const Real g = 511.58;
77  const Real Tc = 2503.7; // critical temperature, K
78  mooseAssert(T < Tc, "Temperature is greater than critical temperature 2503.7 K ");
79 
80  return rhoc + f * (1 - T / Tc) + g * std::sqrt(1 - T / Tc);
81 }
Real SodiumProperties::temperature ( Real  enthalpy) const

Inverse solve for temperature [K] from enthalpy [J/kg].

Parameters
henthalpy

Definition at line 51 of file SodiumProperties.C.

Referenced by drho_dh().

52 {
53  // Estimate initial guess from linear part of enthalpy.
54  Real T = (H + 365.77e3) / 1.6582e3;
55 
56  // Newton-Raphson for this equation: enthalpy(T) - H = 0 = residual. This is easy because
57  // dResidual/dT is just dH/dT, which is heat capacity.
58  for (unsigned iteration = 0; iteration < 10; ++iteration)
59  {
60  Real residual = h(T) - H;
61  T -= residual / heatCapacity(T);
62  if (std::abs(residual / H) < 1e-6)
63  break;
64  }
65  // If we get here, enthalpy is probably out of bounds. However, due to the nature of the JFNK
66  // calculation, we probably just want to ignore the error and spit out a bogus T so that the
67  // solver keeps rolling.
68  return T;
69 }
Real heatCapacity(Real temperature) const
Heat capacity of liquid Na in J/kg-K as a function of temperature.
Real h(Real temperature) const
Enthalpy of liquid Na (relative to solid Na at STP) in J/kg as a function of temperature.

The documentation for this class was generated from the following files: