SodiumSaturationFluidPropertiesTest.C

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

#include "SinglePhaseFluidPropertiesTestUtils.h"
#include "SodiumSaturationFluidPropertiesTest.h"

TEST_F(SodiumSaturationFluidPropertiesTest, fluidName) { EXPECT_EQ(_fp->fluidName(), "sodium_sat"); }

TEST_F(SodiumSaturationFluidPropertiesTest, molarMass)
{
  ABS_TEST(_fp->molarMass(), 22.989769E-3, REL_TOL_SAVED_VALUE);
}

TEST_F(SodiumSaturationFluidPropertiesTest, thermalConductivity)
{
  const Real T = 800.0;
  const Real p = 101325;

  // Fink and Leibowitz (1979) list conductivity as 66.8 at 800 K, so the fit agrees well
  REL_TEST(_fp->k_from_p_T(p, T), 66.9752096, REL_TOL_SAVED_VALUE);
  DERIV_TEST(_fp->k_from_p_T, p, T, REL_TOL_DERIVATIVE);
}

TEST_F(SodiumSaturationFluidPropertiesTest, viscosity)
{
  const Real T = 800.0;
  const Real p = 101325;

  // Fink and Leibowitz (1979) list viscosity as 0.000229 at 800 K, so the fit agrees well
  REL_TEST(_fp->mu_from_p_T(p, T), 0.00022907910937500003, REL_TOL_SAVED_VALUE);
  DERIV_TEST(_fp->mu_from_p_T, p, T, REL_TOL_DERIVATIVE);
}

TEST_F(SodiumSaturationFluidPropertiesTest, isobaricSpecificHeat)
{
  const Real T1 = 400.0;
  const Real T2 = 800.0;
  const Real p = 101325;

  // Fink and Leibowitz (1979) show Cp = 1373.6 at 400 K, so the fit agrees well
  REL_TEST(_fp->cp_from_p_T(p, T1), 1383.985792, REL_TOL_SAVED_VALUE);
  DERIV_TEST(_fp->cp_from_p_T, p, T1, REL_TOL_DERIVATIVE);

  // Fink and Leibowitz (1979) show Cp = 1263.6 at 800 K, so the fit agrees well
  REL_TEST(_fp->cp_from_p_T(p, T2), 1260.719872, REL_TOL_SAVED_VALUE);
  DERIV_TEST(_fp->cp_from_p_T, p, T2, REL_TOL_DERIVATIVE);
}

TEST_F(SodiumSaturationFluidPropertiesTest, isochoricSpecificHeat)
{
  const Real T1 = 400.0;
  const Real T2 = 800.0;
  const Real p = 101325;

  // Fink and Leibowitz (1979) show Cv = 1230.1 at 400 K, so the fit agrees well
  REL_TEST(_fp->cv_from_p_T(p, T1), 1222.9660159999999, REL_TOL_SAVED_VALUE);
  DERIV_TEST(_fp->cv_from_p_T, p, T1, REL_TOL_DERIVATIVE);

  // Fink and Leibowitz (1979) show Cv = 982.611 at 800 K, so the fit agrees well
  REL_TEST(_fp->cv_from_p_T(p, T2), 986.2385279999999, REL_TOL_SAVED_VALUE);
  DERIV_TEST(_fp->cv_from_p_T, p, T2, REL_TOL_DERIVATIVE);
}

TEST_F(SodiumSaturationFluidPropertiesTest, density)
{
  const Real Tm = 370.98;
  const Real T = 800.0;
  const Real Tb = 1156.5;
  const Real p = 101325;

  // Fink and Leibowitz (1979) list density at 370.98 K is 927.3, so the fit agrees well
  ABS_TEST(_fp->rho_from_p_T(p, Tm), 923.3571594322216, REL_TOL_SAVED_VALUE);
  DERIV_TEST(_fp->rho_from_p_T, p, Tm, REL_TOL_DERIVATIVE);
  ABS_TEST(_fp->v_from_p_T(p, Tm), 0.001083004544649772, REL_TOL_SAVED_VALUE);
  DERIV_TEST(_fp->v_from_p_T, p, Tm, REL_TOL_DERIVATIVE);

  // Fink and Leibowitz (1979) list density at 800 K is 825.6, so the fit agrees well
  ABS_TEST(_fp->rho_from_p_T(p, T), 826.04056, REL_TOL_SAVED_VALUE);
  DERIV_TEST(_fp->rho_from_p_T, p, T, REL_TOL_DERIVATIVE);
  ABS_TEST(_fp->v_from_p_T(p, T), 0.0012105943078630425, REL_TOL_SAVED_VALUE);
  DERIV_TEST(_fp->v_from_p_T, p, T, REL_TOL_DERIVATIVE);

  // Fink and Leibowitz (1979) list density at 1156.5 K is 739.4, so the fit agrees well
  ABS_TEST(_fp->rho_from_p_T(p, Tb), 742.0807106065, REL_TOL_SAVED_VALUE);
  DERIV_TEST(_fp->rho_from_p_T, p, Tb, REL_TOL_DERIVATIVE);
  ABS_TEST(_fp->v_from_p_T(p, Tb), 0.0013475623145933863, REL_TOL_SAVED_VALUE);
  DERIV_TEST(_fp->v_from_p_T, p, Tb, REL_TOL_DERIVATIVE);
}

TEST_F(SodiumSaturationFluidPropertiesTest, specificInternalEnergy)
{
  const Real T = 800.0;
  const Real p = 101325;

  ABS_TEST(_fp->e_from_p_T(p, T),
           _fp->h_from_p_T(p, T) - p * _fp->v_from_p_T(p, T),
           REL_TOL_SAVED_VALUE);
  DERIV_TEST(_fp->e_from_p_T, p, T, REL_TOL_DERIVATIVE);
}

TEST_F(SodiumSaturationFluidPropertiesTest, specificEnthalpy)
{
  const Real T = 800.0;
  const Real p = 101325;

  // Fink and Leibowitz (1979) compute enthalpy at 800 K as 768602, so the fit agrees well
  ABS_TEST(_fp->h_from_p_T(p, T), 767034.6715200001, REL_TOL_SAVED_VALUE);
  DERIV_TEST(_fp->h_from_p_T, p, T, REL_TOL_DERIVATIVE);

  // Fink and Leibowitz (1979) compute enthalpy at 1156.5 K as 1218803, so the fit agrees well
  const Real T2 = 1156.5;
  ABS_TEST(_fp->h_from_p_T(p, T2), 1218897.6776327428, REL_TOL_SAVED_VALUE);
  DERIV_TEST(_fp->h_from_p_T, p, T2, REL_TOL_DERIVATIVE);
}

TEST_F(SodiumSaturationFluidPropertiesTest, PTFromVE)
{
  // Starting pressure and temperature
  const Real T = 800.0;
  const Real p = 101325;

  // Obtain specific volume and specific energy
  const Real v = _fp->v_from_p_T(p, T);
  const Real e = _fp->e_from_p_T(p, T);

  // Verify that the conversion back succeeds
  ABS_TEST(_fp->T_from_v_e(v, e), T, TOLERANCE);
  ABS_TEST(_fp->p_from_v_e(v, e), p, 1e-6);
}