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

TEST_F(HeliumFluidPropertiesTest, fluidName) { EXPECT_EQ(_fp->fluidName(), "helium"); }

TEST_F(HeliumFluidPropertiesTest, molarMass)
{
  ABS_TEST(_fp->molarMass(), 4.002602e-3, REL_TOL_SAVED_VALUE);
}

TEST_F(HeliumFluidPropertiesTest, thermalConductivity)
{
  const Real T = 120.0 + 273.15;
  const Real p = 101325.0;
  const Real e = _fp->e_from_p_T(p, T);
  const Real v = 1. / _fp->rho_from_p_T(p, T);

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

TEST_F(HeliumFluidPropertiesTest, viscosity)
{
  const Real T = 120.0 + 273.15;
  const Real p = 101325.0;
  const Real e = _fp->e_from_p_T(p, T);
  const Real v = 1. / _fp->rho_from_p_T(p, T);

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

TEST_F(HeliumFluidPropertiesTest, isobaricSpecificHeat)
{
  const Real T = 120.0 + 273.15;
  const Real p = 101325.0;
  const Real e = _fp->e_from_p_T(p, T);
  const Real v = 1. / _fp->rho_from_p_T(p, T);

  REL_TEST(_fp->cp_from_p_T(p, T), 5195.0, REL_TOL_SAVED_VALUE);
  REL_TEST(_fp->cp_from_v_e(v, e), 5195.0, 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(HeliumFluidPropertiesTest, isochoricSpecificHeat)
{
  const Real T = 120.0 + 273.15;
  const Real p = 101325.0;
  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), 3117.0, REL_TOL_SAVED_VALUE);
  REL_TEST(_fp->cv_from_p_T(p, T), 3117.0, REL_TOL_SAVED_VALUE);
  DERIV_TEST(_fp->cv_from_p_T, p, T, REL_TOL_DERIVATIVE);
  DERIV_TEST(_fp->cv_from_v_e, v, e, REL_TOL_DERIVATIVE);
}

TEST_F(HeliumFluidPropertiesTest, density)
{
  const Real T = 120.0 + 273.15;
  const Real p = 101325.0;

  ABS_TEST(_fp->rho_from_p_T(p, T), 1.2402629878970449e-01, REL_TOL_SAVED_VALUE);
  DERIV_TEST(_fp->rho_from_p_T, p, T, REL_TOL_DERIVATIVE);
}

TEST_F(HeliumFluidPropertiesTest, specificInternalEnergy)
{
  const Real T = 120.0 + 273.15;
  const Real p = 101325.0;
  const Real rho = _fp->rho_from_p_T(p, T);

  ABS_TEST(_fp->e_from_p_T(p, T), 1.2254485499999998e6, REL_TOL_SAVED_VALUE);
  DERIV_TEST(_fp->e_from_p_T, p, T, REL_TOL_DERIVATIVE);
  ABS_TEST(_fp->e_from_p_rho(p, rho), 1.2254485499999998e6, REL_TOL_SAVED_VALUE);

  // Equivalent (v,h), (p,T)
  const Real v = 1. / rho;
  const Real h = _fp->h_from_p_T(p, T);
  ABS_TEST(_fp->e_from_v_h(v, h), 1.2254485499999998e6, REL_TOL_SAVED_VALUE);
  DERIV_TEST(_fp->e_from_v_h, v, h, REL_TOL_DERIVATIVE);
}

TEST_F(HeliumFluidPropertiesTest, pressure)
{
  const Real T = 120.0 + 273.15;
  const Real p = 101325.0;
  const Real e = _fp->e_from_p_T(p, T);
  const Real v = 1. / _fp->rho_from_p_T(p, T);

  REL_TEST(_fp->p_from_v_e(v, e), p, REL_TOL_SAVED_VALUE);
  DERIV_TEST(_fp->p_from_v_e, v, e, REL_TOL_DERIVATIVE);
}

TEST_F(HeliumFluidPropertiesTest, temperature)
{
  const Real T = 120.0 + 273.15;
  const Real p = 101325.0;
  const Real e = _fp->e_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);

  const Real h = _fp->h_from_p_T(p, T);
  REL_TEST(_fp->T_from_p_h(p, h), T, REL_TOL_SAVED_VALUE);
}

TEST_F(HeliumFluidPropertiesTest, specificEnthalpy)
{
  const Real T = 120.0 + 273.15;
  const Real p = 101325.0;

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

  // Derivatives
  const ADReal T_ad = T;
  const ADReal p_ad = p;
  Real h, dhdT, dhdp;
  ADReal h_ad, dhdT_ad, dhdp_ad;
  _fp->h_from_p_T(p, T, h, dhdT, dhdp);
  _fp->h_from_p_T(p_ad, T_ad, h_ad, dhdT_ad, dhdp_ad);
  DERIV_TEST(_fp->h_from_p_T, p, T, REL_TOL_DERIVATIVE);
  REL_TEST(h, h_ad.value(), REL_TOL_CONSISTENCY);
  REL_TEST(dhdT, dhdT_ad.value(), REL_TOL_CONSISTENCY);
  REL_TEST(dhdp, dhdp_ad.value(), REL_TOL_CONSISTENCY);
}