PorousFlowWaterVaporTest.C File Reference

Go to the source code of this file.

Functions
	TEST_F (PorousFlowWaterVaporTest, name)
	Verify that the correct name is supplied. More...
	TEST_F (PorousFlowWaterVaporTest, properties)

Function Documentation

TEST_F() [1/2]

TEST_F	(	PorousFlowWaterVaporTest	,
		name
	)

Verify that the correct name is supplied.

Definition at line 16 of file PorousFlowWaterVaporTest.C.

{ EXPECT_EQ("water-vapor", _fp->fluidStateName()); }

TEST_F() [2/2]

TEST_F	(	PorousFlowWaterVaporTest	,
		properties
	)

Definition at line 18 of file PorousFlowWaterVaporTest.C.

{
  FluidStatePhaseEnum phase_state;
  std::vector<FluidStateProperties> fsp(2, FluidStateProperties(1));

  // Single phase liquid region
  Real p = 1.0e6;
  Real h = 1.0e5;

  const Real tol = REL_TOL_CONSISTENCY;

  _fp->thermophysicalProperties(p, h, 0, phase_state, fsp);
  EXPECT_EQ(phase_state, FluidStatePhaseEnum::LIQUID);

  Real T = _water_fp->T_from_p_h(p, h);

  ABS_TEST(fsp[0].density.value(), _water_fp->rho_from_p_T(p, T), tol);
  ABS_TEST(fsp[0].viscosity.value(), _water_fp->mu_from_p_T(p, T), tol);
  ABS_TEST(fsp[0].enthalpy.value(), h, tol);
  ABS_TEST(fsp[0].internal_energy.value(), _water_fp->e_from_p_T(p, T), tol);
  ABS_TEST(fsp[0].saturation.value(), 1.0, tol);
  ABS_TEST(fsp[0].pressure.value(), p, tol);

  // Single phase gas region
  h = 3.0e6;
  Real pg = p + _pc->capillaryPressure(0.0);

  _fp->thermophysicalProperties(p, h, 0, phase_state, fsp);
  EXPECT_EQ(phase_state, FluidStatePhaseEnum::GAS);

  T = _water_fp->T_from_p_h(pg, h);
  ABS_TEST(fsp[1].density.value(), _water_fp->rho_from_p_T(pg, T), tol);
  ABS_TEST(fsp[1].viscosity.value(), _water_fp->mu_from_p_T(pg, T), tol);
  ABS_TEST(fsp[1].enthalpy.value(), h, tol);
  ABS_TEST(fsp[1].internal_energy.value(), _water_fp->e_from_p_T(pg, T), tol);
  ABS_TEST(fsp[1].saturation.value(), 1.0, tol);
  ABS_TEST(fsp[1].pressure.value(), pg, tol);

  // Two-phase region
  h = 8.0e5;

  _fp->thermophysicalProperties(p, h, 0, phase_state, fsp);
  EXPECT_EQ(phase_state, FluidStatePhaseEnum::TWOPHASE);

  const Real Tsat = _water_fp->vaporTemperature(p);
  const Real dT = 1.0e-6;

  // Liquid phase
  ABS_TEST(fsp[0].density.value(), _water_fp->rho_from_p_T(p, Tsat - dT), tol);
  ABS_TEST(fsp[0].viscosity.value(), _water_fp->mu_from_p_T(p, Tsat - dT), tol);
  ABS_TEST(fsp[0].enthalpy.value(), _water_fp->h_from_p_T(p, Tsat - dT), tol);
  ABS_TEST(fsp[0].internal_energy.value(), _water_fp->e_from_p_T(p, Tsat - dT), tol);
  ABS_TEST(fsp[0].pressure.value(), p, tol);

  // Gas phase
  pg = p + _pc->capillaryPressure(fsp[0].saturation.value());
  ABS_TEST(fsp[1].density.value(), _water_fp->rho_from_p_T(p, Tsat + dT), tol);
  ABS_TEST(fsp[1].viscosity.value(), _water_fp->mu_from_p_T(p, Tsat + dT), tol);
  ABS_TEST(fsp[1].enthalpy.value(), _water_fp->h_from_p_T(p, Tsat + dT), tol);
  ABS_TEST(fsp[1].internal_energy.value(), _water_fp->e_from_p_T(p, Tsat + dT), tol);
  ABS_TEST(fsp[1].pressure.value(), pg, tol);
}