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:

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
virtual void	threadJoin (const UserObject &) final
virtual void	subdomainSetup () final

Static Public Attributes
static const Real	_R = 8.3144598
	Universal gas constant (J/mol/K) More...

Protected Attributes
const Real	_T_c2k
	Conversion of temperature from Celsius to Kelvin. More...
const bool	_allow_imperfect_jacobians
	Flag to set unimplemented Jacobian entries to zero. More...

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 24 of file SodiumProperties.h.

Constructor & Destructor Documentation

SodiumProperties()

SodiumProperties::SodiumProperties

(

const InputParameters &

parameters

)

Definition at line 24 of file SodiumProperties.C.

                                                                     : FluidProperties(parameters)
{
}

Member Function Documentation

drho_dh()

Real SodiumProperties::drho_dh

(

Real

enthalpy

)

const

Derivative of density w.r.t enthalpy.

Parameters

h

enthalpy

Definition at line 99 of file SodiumProperties.C.

{
  Real T = 0.0;
  T = temperature(h);
  return drho_dT(T) / heatCapacity(T);
}

drho_dT()

Real SodiumProperties::drho_dT

(

Real

temperature

)

const

Derivative of density w.r.t temperature.

Parameters

T	temperature

Definition at line 88 of file SodiumProperties.C.

{
  const Real f = 275.32;
  const Real g = 511.58;
  const Real Tc = 2503.7; // critical temperature, K
  mooseAssert(T < Tc, "Temperature is greater than critical temperature 2503.7 K ");
 
  return -(f + g * (0.5) / std::sqrt(1 - T / Tc)) / Tc;
}

Referenced by drho_dh().

execute()

virtual void FluidProperties::execute ( )

inlinefinalvirtualinherited

Definition at line 34 of file FluidProperties.h.

{}

finalize()

virtual void FluidProperties::finalize ( )

inlinefinalvirtualinherited

Definition at line 36 of file FluidProperties.h.

{}

h()

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

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

Definition at line 37 of file SodiumProperties.C.

{
  const Real T2 = T * T;
  const Real T3 = T2 * T;
 
  // Converted from kJ/kg to J/kg.
  return -365.77e3 + 1.6582e3 * T - 4.2395e-1 * T2 + 1.4847e-4 * T3 + 2992.6e3 / T;
}

Referenced by drho_dh(), and temperature().

heatCapacity()

Real SodiumProperties::heatCapacity

(

Real

temperature

)

const

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

Parameters

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

Definition at line 47 of file SodiumProperties.C.

{
  const Real T2 = T * T;
  // Converted from kJ/kg-K to J/kg-K.
  return 1.6582e3 - 8.4790e-1 * T + 4.4541e-4 * T2 - 2992.6e3 / T2;
}

Referenced by drho_dh(), and temperature().

initialize()

virtual void FluidProperties::initialize ( )

inlinefinalvirtualinherited

Definition at line 35 of file FluidProperties.h.

{}

k()

Real SodiumProperties::k

(

Real

temperature

)

const

Thermal conductivity as a function of temperature.

Parameters

T	temperature From page 181. Valid from 371 to 1500 K.

Definition at line 29 of file SodiumProperties.C.

{
  const Real T2 = T * T;
  const Real T3 = T2 * T;
  return 124.67 - 0.11381 * T + 5.5226e-5 * T2 - 1.1842e-8 * T3;
}

rho()

Real SodiumProperties::rho

(

Real

temperature

)

const

Density as a function of temperature.

Parameters

T	temperature

Definition at line 76 of file SodiumProperties.C.

{
  const Real rhoc = 219.0; // kg/m^3
  const Real f = 275.32;
  const Real g = 511.58;
  const Real Tc = 2503.7; // critical temperature, K
  mooseAssert(T < Tc, "Temperature is greater than critical temperature 2503.7 K ");
 
  return rhoc + f * (1 - T / Tc) + g * std::sqrt(1 - T / Tc);
}

subdomainSetup()

virtual void FluidProperties::subdomainSetup ( )

inlinefinalvirtualinherited

Definition at line 39 of file FluidProperties.h.

{}

temperature()

Real SodiumProperties::temperature

(

Real

enthalpy

)

const

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

Parameters

h

enthalpy

Definition at line 55 of file SodiumProperties.C.

{
  // Estimate initial guess from linear part of enthalpy.
  Real T = (H + 365.77e3) / 1.6582e3;
 
  // Newton-Raphson for this equation: enthalpy(T) - H = 0 = residual. This is easy because
  // dResidual/dT is just dH/dT, which is heat capacity.
  for (unsigned iteration = 0; iteration < 10; ++iteration)
  {
    Real residual = h(T) - H;
    T -= residual / heatCapacity(T);
    if (std::abs(residual / H) < 1e-6)
      break;
  }
  // If we get here, enthalpy is probably out of bounds. However, due to the nature of the JFNK
  // calculation, we probably just want to ignore the error and spit out a bogus T so that the
  // solver keeps rolling.
  return T;
}

Referenced by drho_dh().

threadJoin()

virtual void FluidProperties::threadJoin ( const UserObject &  )

inlinefinalvirtualinherited

Definition at line 38 of file FluidProperties.h.

{}

Member Data Documentation

_allow_imperfect_jacobians

const bool FluidProperties::_allow_imperfect_jacobians

protectedinherited

Flag to set unimplemented Jacobian entries to zero.

Definition at line 48 of file FluidProperties.h.

Referenced by SinglePhaseFluidProperties::fluidPropError().

_R

const Real FluidProperties::_R = 8.3144598

staticinherited

Universal gas constant (J/mol/K)

Definition at line 42 of file FluidProperties.h.

Referenced by HelmholtzFluidProperties::c_from_p_T(), IdealGasFluidProperties::c_from_p_T(), HelmholtzFluidProperties::cp_from_p_T(), HelmholtzFluidProperties::cv_from_p_T(), HelmholtzFluidProperties::e_from_p_T(), HelmholtzFluidProperties::h_from_p_T(), HelmholtzFluidProperties::p_from_rho_T(), HelmholtzFluidProperties::rho_from_p_T(), IdealGasFluidProperties::rho_from_p_T(), and HelmholtzFluidProperties::s_from_p_T().

_T_c2k

const Real FluidProperties::_T_c2k

protectedinherited

Conversion of temperature from Celsius to Kelvin.

Definition at line 46 of file FluidProperties.h.

Referenced by NaClFluidProperties::cp_from_p_T(), BrineFluidProperties::cp_from_p_T_X(), NaClFluidProperties::h_from_p_T(), BrineFluidProperties::h_from_p_T_X(), BrineFluidProperties::haliteSolubility(), NaClFluidProperties::k_from_p_T(), BrineFluidProperties::k_from_p_T_X(), BrineFluidProperties::mu_from_p_T_X(), CO2FluidProperties::partialDensity(), NaClFluidProperties::rho_from_p_T(), and BrineFluidProperties::rho_from_p_T_X().

---

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

• fluid_properties/include/userobjects/SodiumProperties.h
• fluid_properties/src/userobjects/SodiumProperties.C