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

TEST_F(FlibeFluidPropertiesTest, fluidName) { EXPECT_EQ(_fp->fluidName(), "flibe"); }

TEST_F(FlibeFluidPropertiesTest, molarMass)
{
  ABS_TEST(_fp->molarMass(), 99.037703E-3, REL_TOL_SAVED_VALUE);
}

TEST_F(FlibeFluidPropertiesTest, thermalConductivity)
{
  const Real T = 800.0;
  const Real p = 2.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 = 1.030000000000000027e+00;

  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(FlibeFluidPropertiesTest, viscosity)
{
  const Real T = 800.0;
  const Real p = 2.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 = 1.267440888533459542e-02;

  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(FlibeFluidPropertiesTest, isobaricSpecificHeat)
{
  const Real T = 800.0;
  const Real p = 2.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 = 2416.;

  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(FlibeFluidPropertiesTest, isochoricSpecificHeat)
{
  const Real T = 800.0;
  const Real p = 2.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), 1.0219369260016974e+03, REL_TOL_SAVED_VALUE);
  DERIV_TEST(_fp->cv_from_v_e, v, e, REL_TOL_DERIVATIVE);
}

TEST_F(FlibeFluidPropertiesTest, density)
{
  const Real T = 800.0;
  const Real p = 2.0 * 101325;

  ABS_TEST(_fp->rho_from_p_T(p, T), 2.0222975535430000e+03, REL_TOL_SAVED_VALUE);
  DERIV_TEST(_fp->rho_from_p_T, p, T, REL_TOL_DERIVATIVE);
}

TEST_F(FlibeFluidPropertiesTest, specificInternalEnergy)
{
  const Real T = 800.0;
  const Real p = 2.0 * 101325;
  const Real e = 1.9326997921944547e+06;

  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(FlibeFluidPropertiesTest, pressure)
{
  const Real T = 800.0;
  const Real p = 2.0 * 101325;
  const Real e = _fp->e_from_p_T(p, T);
  const Real v = 1. / _fp->rho_from_p_T(p, T);

  // A slightly larger tolerance is required for the 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, 2e-11);
  DERIV_TEST(_fp->p_from_v_e, v, e, REL_TOL_DERIVATIVE);
}

TEST_F(FlibeFluidPropertiesTest, temperature)
{
  const Real T = 800.0;
  const Real p = 2.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(FlibeFluidPropertiesTest, specificEnthalpy)
{
  const Real T = 800.0;
  const Real p = 2.0 * 101325;

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