FlinakFluidPropertiesTest.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 "FlinakFluidPropertiesTest.h"

TEST_F(FlinakFluidPropertiesTest, fluidName) { EXPECT_EQ(_fp->fluidName(), "flinak"); }

TEST_F(FlinakFluidPropertiesTest, molarMass)
{
  ABS_TEST(_fp->molarMass(), 41.291077435E-3, REL_TOL_SAVED_VALUE);
}

TEST_F(FlinakFluidPropertiesTest, thermalConductivity)
{
  const Real T = 800.0;
  const Real p = 3.0 * 101325;
  const Real e = _fp->e_from_p_T(p, T);
  const Real v = 1. / _fp->rho_from_p_T(p, T);
  const Real k = 0.83;

  REL_TEST(_fp->k_from_v_e(v, e), k, REL_TOL_SAVED_VALUE);
  REL_TEST(_fp->k_from_p_T(p, T), k, REL_TOL_SAVED_VALUE);
  DERIV_TEST(_fp->k_from_p_T, p, T, REL_TOL_DERIVATIVE);
}

TEST_F(FlinakFluidPropertiesTest, viscosity)
{
  const Real T = 800.0;
  const Real p = 3.0 * 101325;
  const Real e = _fp->e_from_p_T(p, T);
  const Real v = 1. / _fp->rho_from_p_T(p, T);
  const Real mu = 0.0073420946394096;

  REL_TEST(_fp->mu_from_v_e(v, e), mu, REL_TOL_SAVED_VALUE);
  REL_TEST(_fp->mu_from_p_T(p, T), mu, REL_TOL_SAVED_VALUE);
  DERIV_TEST(_fp->mu_from_p_T, p, T, REL_TOL_DERIVATIVE);
}

TEST_F(FlinakFluidPropertiesTest, isobaricSpecificHeat)
{
  const Real T = 800.0;
  const Real p = 3.0 * 101325;
  const Real e = _fp->e_from_p_T(p, T);
  const Real v = 1. / _fp->rho_from_p_T(p, T);
  const Real cp = 2010.;

  REL_TEST(_fp->cp_from_v_e(v, e), cp, REL_TOL_SAVED_VALUE);
  REL_TEST(_fp->cp_from_p_T(p, T), cp, REL_TOL_SAVED_VALUE);
  DERIV_TEST(_fp->cp_from_p_T, p, T, REL_TOL_DERIVATIVE);
  DERIV_TEST(_fp->cp_from_v_e, v, e, REL_TOL_DERIVATIVE);
}

TEST_F(FlinakFluidPropertiesTest, isochoricSpecificHeat)
{
  const Real T = 800.0;
  const Real p = 3.0 * 101325;
  const Real e = _fp->e_from_p_T(p, T);
  const Real v = 1. / _fp->rho_from_p_T(p, T);

  REL_TEST(_fp->cv_from_v_e(v, e), 45.55316141330309, REL_TOL_SAVED_VALUE);
  DERIV_TEST(_fp->cv_from_v_e, v, e, REL_TOL_DERIVATIVE);
}

TEST_F(FlinakFluidPropertiesTest, density)
{
  const Real T = 800.0;
  const Real p = 3.0 * 101325;

  ABS_TEST(_fp->rho_from_p_T(p, T), 2145.035107086000153, REL_TOL_SAVED_VALUE);
  DERIV_TEST(_fp->rho_from_p_T, p, T, 4.0e-6);
}

TEST_F(FlinakFluidPropertiesTest, specificInternalEnergy)
{
  const Real T = 800.0;
  const Real p = 3.0 * 101325;
  const Real e = 1607858.289032661588863;

  ABS_TEST(_fp->e_from_p_T(p, T), e, REL_TOL_SAVED_VALUE);
  DERIV_TEST(_fp->e_from_p_T, p, T, REL_TOL_DERIVATIVE);

  ABS_TEST(_fp->e_from_p_rho(p, _fp->rho_from_p_T(p, T)), e, REL_TOL_SAVED_VALUE);
}

TEST_F(FlinakFluidPropertiesTest, pressure)
{
  const Real T = 800.0;
  const Real p = 3.0 * 101325;
  const Real e = _fp->e_from_p_T(p, T);
  const Real v = 1. / _fp->rho_from_p_T(p, T);

  // test pressure and its derivatives with respect to v and e. We need a slightly
  // larger tolerance for pressure because the compressibility is so small we are
  // polluted by a little numerical error.
  REL_TEST(_fp->p_from_v_e(v, e), p, 6e-11);
  DERIV_TEST(_fp->p_from_v_e, v, e, 2e-9);
}

TEST_F(FlinakFluidPropertiesTest, temperature)
{
  const Real T = 800.0;
  const Real p = 3.0 * 101325;
  const Real e = _fp->e_from_p_T(p, T);
  const Real h = _fp->h_from_p_T(p, T);
  const Real v = 1. / _fp->rho_from_p_T(p, T);

  REL_TEST(_fp->T_from_v_e(v, e), T, REL_TOL_SAVED_VALUE);
  DERIV_TEST(_fp->T_from_v_e, v, e, REL_TOL_DERIVATIVE);

  REL_TEST(_fp->T_from_p_rho(p, 1. / v), T, REL_TOL_SAVED_VALUE);
  REL_TEST(_fp->T_from_p_h(p, h), T, REL_TOL_SAVED_VALUE);
}

TEST_F(FlinakFluidPropertiesTest, specificEnthalpy)
{
  const Real T = 800.0;
  const Real p = 3.0 * 101325;

  ABS_TEST(_fp->h_from_p_T(p, T), 1.608e+06, REL_TOL_SAVED_VALUE);
  DERIV_TEST(_fp->h_from_p_T, p, T, REL_TOL_DERIVATIVE);
}