TestNumericalFluxGasMixHLLC Class Reference

Tests NumericalFluxGasMixHLLC. More...

#include <TestNumericalFluxGasMixHLLC.h>

Inheritance diagram for TestNumericalFluxGasMixHLLC:

Public Member Functions
	TestNumericalFluxGasMixHLLC ()

Protected Member Functions
virtual const NumericalFlux1D &	createFluxObject () override
	Creates the flux object to be tested. More...
virtual std::vector< std::pair< std::vector< ADReal >, std::vector< ADReal > > >	getPrimitiveSolutionsSymmetryTest () const override
	Gets a vector of pairs of primitive solution vectors to use for symmetry test. More...
virtual std::vector< std::vector< ADReal > >	getPrimitiveSolutionsConsistencyTest () const override
	Gets a vector of primitive solution vectors to use for consistency test. More...
virtual std::vector< ADReal >	computeConservativeSolution (const std::vector< ADReal > &W, const ADReal &A) const override
	Computes the conservative solution from the primitive solution. More...
virtual std::vector< ADReal >	computeFluxFromPrimitive (const std::vector< ADReal > &W, const ADReal &A) const override
	Computes the 1D flux vector from the primitive solution. More...
void	addFluidProperties ()
	Adds fluid properties objects needed for testing. More...
void	testConsistency ()
	Runs the consistency test(s) More...
void	testSymmetry ()
	Runs the symmetry test(s) More...
void	buildObjects ()
T &	addObject (const std::string &type, const std::string &name, InputParameters &params)

Protected Attributes
const UserObjectName	_fp_mix_name
	Mixture fluid properties name. More...
const IdealGasMixtureFluidProperties *	_fp_mix
	Fluid properties user object. More...
std::unique_ptr< MooseMesh >	_mesh
std::shared_ptr< MooseApp >	_app
Factory &	_factory
std::shared_ptr< FEProblem >	_fe_problem

Detailed Description

Tests NumericalFluxGasMixHLLC.

Definition at line 17 of file TestNumericalFluxGasMixHLLC.h.

Constructor & Destructor Documentation

TestNumericalFluxGasMixHLLC()

TestNumericalFluxGasMixHLLC::TestNumericalFluxGasMixHLLC

(

)

Definition at line 24 of file TestNumericalFluxGasMixHLLC.C.

                                                         : TestNumericalFluxGasMixBase()
{
}

Member Function Documentation

addFluidProperties()

void TestNumericalFluxGasMixBase::addFluidProperties ( )

protectedinherited

Adds fluid properties objects needed for testing.

Definition at line 36 of file TestNumericalFluxGasMixBase.C.

Referenced by TestNumericalFluxGasMixBase::TestNumericalFluxGasMixBase().

{
  const std::string fp_steam_name = "fp_steam";
  const std::string fp_nitrogen_name = "fp_nitrogen";

  // steam fluid properties; parameters correspond to T in range 298 K to 473 K
  {
    const std::string class_name = "IdealGasFluidProperties";
    InputParameters params = _factory.getValidParams(class_name);
    params.set<Real>("gamma") = 1.43;
    params.set<Real>("molar_mass") = 0.01801488;
    params.set<Real>("mu") = 0.000013277592; // at 400 K and 1.e5 Pa
    params.set<Real>("k") = 0.026824977826;  // at 400 K and 1.e5 Pa
    _fe_problem->addUserObject(class_name, fp_steam_name, params);
  }

  // nitrogen fluid properties
  {
    const std::string class_name = "IdealGasFluidProperties";
    InputParameters params = _factory.getValidParams(class_name);
    params.set<Real>("gamma") = 1.4;
    params.set<Real>("molar_mass") = 0.028012734746133888;
    params.set<Real>("mu") = 0.0000222084; // at 400 K and 1.e5 Pa
    params.set<Real>("k") = 0.032806168;   // at 400 K and 1.e5 Pa
    _fe_problem->addUserObject(class_name, fp_nitrogen_name, params);
  }

  // mixture fluid properties
  {
    const std::string class_name = "IdealGasMixtureFluidProperties";
    InputParameters params = _factory.getValidParams(class_name);
    params.set<std::vector<UserObjectName>>("component_fluid_properties") = {fp_steam_name, fp_nitrogen_name};
    _fe_problem->addUserObject(class_name, _fp_mix_name, params);
    _fp_mix = &_fe_problem->getUserObject<IdealGasMixtureFluidProperties>(_fp_mix_name);
  }
}

computeConservativeSolution()

std::vector< ADReal > TestNumericalFluxGasMixBase::computeConservativeSolution ( const std::vector< ADReal > &  W, const ADReal &  A  ) const

overrideprotectedvirtualinherited

Computes the conservative solution from the primitive solution.

Parameters

[in]	W	Primitive solution vector
[in]	A	Cross-sectional area

Implements TestNumericalFlux1D.

Definition at line 22 of file TestNumericalFluxGasMixBase.C.

{
  return FlowModelGasMixUtils::computeConservativeSolution<true>(W, A, *_fp_mix);
}

computeFluxFromPrimitive()

std::vector< ADReal > TestNumericalFluxGasMixBase::computeFluxFromPrimitive ( const std::vector< ADReal > &  W, const ADReal &  A  ) const

overrideprotectedvirtualinherited

Computes the 1D flux vector from the primitive solution.

Parameters

[in]	W	Primitive solution vector
[in]	A	Cross-sectional area

Implements TestNumericalFlux1D.

Definition at line 29 of file TestNumericalFluxGasMixBase.C.

{
  return FlowModelGasMixUtils::computeFluxFromPrimitive<true>(W, A, *_fp_mix);
}

createFluxObject()

const NumericalFlux1D & TestNumericalFluxGasMixHLLC::createFluxObject ( )

overrideprotectedvirtual

Creates the flux object to be tested.

Implements TestNumericalFlux1D.

Definition at line 29 of file TestNumericalFluxGasMixHLLC.C.

{
  const std::string class_name = "NumericalFluxGasMixHLLC";
  const std::string flux_mix_name = "flux_mix";
  InputParameters params = _factory.getValidParams(class_name);
  params.set<UserObjectName>("fluid_properties") = _fp_mix_name;
  _fe_problem->addUserObject(class_name, flux_mix_name, params);

  return _fe_problem->getUserObject<NumericalFluxGasMixHLLC>(flux_mix_name);
}

getPrimitiveSolutionsConsistencyTest()

std::vector< std::vector< ADReal > > TestNumericalFluxGasMixHLLC::getPrimitiveSolutionsConsistencyTest ( ) const

overrideprotectedvirtual

Gets a vector of primitive solution vectors to use for consistency test.

Implements TestNumericalFlux1D.

Definition at line 77 of file TestNumericalFluxGasMixHLLC.C.

{
  std::vector<ADReal> W1(THMGasMix1D::N_PRIM_VARS);
  W1[THMGasMix1D::MASS_FRACTION] = 0.4;
  W1[THMGasMix1D::PRESSURE] = 1e5;
  W1[THMGasMix1D::TEMPERATURE] = 300;
  W1[THMGasMix1D::VELOCITY] = 1.5;

  std::vector<std::vector<ADReal>> W_list;
  W_list.push_back(W1);

  return W_list;
}

getPrimitiveSolutionsSymmetryTest()

std::vector< std::pair< std::vector< ADReal >, std::vector< ADReal > > > TestNumericalFluxGasMixHLLC::getPrimitiveSolutionsSymmetryTest ( ) const

overrideprotectedvirtual

Gets a vector of pairs of primitive solution vectors to use for symmetry test.

Implements TestNumericalFlux1D.

Definition at line 41 of file TestNumericalFluxGasMixHLLC.C.

{
  // sL < 0 < sM
  std::vector<ADReal> W1(THMGasMix1D::N_PRIM_VARS);
  W1[THMGasMix1D::MASS_FRACTION] = 0.4;
  W1[THMGasMix1D::PRESSURE] = 1e5;
  W1[THMGasMix1D::TEMPERATURE] = 300;
  W1[THMGasMix1D::VELOCITY] = 20;

  std::vector<ADReal> W2(THMGasMix1D::N_PRIM_VARS);
  W2[THMGasMix1D::MASS_FRACTION] = 0.7;
  W2[THMGasMix1D::PRESSURE] = 2e5;
  W2[THMGasMix1D::TEMPERATURE] = 310;
  W2[THMGasMix1D::VELOCITY] = 1.2;

  // sL > 0
  std::vector<ADReal> W3(THMGasMix1D::N_PRIM_VARS);
  W3[THMGasMix1D::MASS_FRACTION] = 0.4;
  W3[THMGasMix1D::PRESSURE] = 1e5;
  W3[THMGasMix1D::TEMPERATURE] = 300;
  W3[THMGasMix1D::VELOCITY] = 500;

  std::vector<ADReal> W4(THMGasMix1D::N_PRIM_VARS);
  W4[THMGasMix1D::MASS_FRACTION] = 0.7;
  W4[THMGasMix1D::PRESSURE] = 2e5;
  W4[THMGasMix1D::TEMPERATURE] = 310;
  W4[THMGasMix1D::VELOCITY] = 550;

  std::vector<std::pair<std::vector<ADReal>, std::vector<ADReal>>> W_pairs;
  W_pairs.push_back(std::pair<std::vector<ADReal>, std::vector<ADReal>>(W1, W2));
  W_pairs.push_back(std::pair<std::vector<ADReal>, std::vector<ADReal>>(W3, W4));

  return W_pairs;
}

testConsistency()

void TestNumericalFlux1D::testConsistency ( )

protectedinherited

Runs the consistency test(s)

Definition at line 60 of file TestNumericalFlux1D.C.

{
  const auto & flux = createFluxObject();

  unsigned int i_side = 0;
  const auto W_list = getPrimitiveSolutionsConsistencyTest();
  for (const auto & W : W_list)
  {
    const ADReal A = 2.0;

    const auto U = computeConservativeSolution(W, A);
    const auto F_expected = computeFluxFromPrimitive(W, A);

    const auto & FL_computed_pos = flux.getFlux(i_side, 0, true, U, U, 1.0);
    const auto & FR_computed_pos = flux.getFlux(i_side, 0, false, U, U, 1.0);
    i_side++;

    const auto & FL_computed_neg = flux.getFlux(i_side, 0, true, U, U, -1.0);
    const auto & FR_computed_neg = flux.getFlux(i_side, 0, false, U, U, -1.0);
    i_side++;

    for (unsigned int i = 0; i < FL_computed_pos.size(); ++i)
    {
      REL_TEST(FL_computed_pos[i], F_expected[i], REL_TOL_CONSISTENCY);
      REL_TEST(FR_computed_pos[i], F_expected[i], REL_TOL_CONSISTENCY);

      REL_TEST(FL_computed_neg[i], F_expected[i], REL_TOL_CONSISTENCY);
      REL_TEST(FR_computed_neg[i], F_expected[i], REL_TOL_CONSISTENCY);
    }
  }
}

testSymmetry()

void TestNumericalFlux1D::testSymmetry ( )

protectedinherited

Runs the symmetry test(s)

Definition at line 19 of file TestNumericalFlux1D.C.

{
  const auto & flux = createFluxObject();

  std::set<unsigned int> flux_regions;

  unsigned int i_side = 0;
  const auto W_pairs = getPrimitiveSolutionsSymmetryTest();
  for (const auto & W_pair : W_pairs)
  {
    const auto & WL = W_pair.first;
    const auto & WR = W_pair.second;

    const ADReal AL = 1.0;
    const ADReal AR = 1.5;

    const std::vector<ADReal> UL = computeConservativeSolution(WL, AL);
    const std::vector<ADReal> UR = computeConservativeSolution(WR, AR);

    const auto FLR = flux.getFlux(i_side, 0, true, UL, UR, 1.0);
    const auto FRL = flux.getFlux(i_side, 0, false, UL, UR, 1.0);
    i_side++;
    flux_regions.insert(flux.getLastRegionIndex());

    const auto FRL_flipped = flux.getFlux(i_side, 0, true, UR, UL, -1.0);
    const auto FLR_flipped = flux.getFlux(i_side, 0, false, UR, UL, -1.0);
    i_side++;
    flux_regions.insert(flux.getLastRegionIndex());

    for (unsigned int i = 0; i < FLR.size(); ++i)
    {
      REL_TEST(FLR[i], FLR_flipped[i], REL_TOL_CONSISTENCY);
      REL_TEST(FRL[i], FRL_flipped[i], REL_TOL_CONSISTENCY);
    }
  }

  // Check that all of the regions in the flux have been tested.
  EXPECT_TRUE(flux_regions.size() == flux.getNumberOfRegions());
}

Member Data Documentation

_fp_mix

const IdealGasMixtureFluidProperties* TestNumericalFluxGasMixBase::_fp_mix

protectedinherited

Fluid properties user object.

Definition at line 39 of file TestNumericalFluxGasMixBase.h.

Referenced by TestNumericalFluxGasMixBase::addFluidProperties(), TestNumericalFluxGasMixBase::computeConservativeSolution(), and TestNumericalFluxGasMixBase::computeFluxFromPrimitive().

_fp_mix_name

const UserObjectName TestNumericalFluxGasMixBase::_fp_mix_name

protectedinherited

Mixture fluid properties name.

Definition at line 36 of file TestNumericalFluxGasMixBase.h.

Referenced by TestNumericalFluxGasMixBase::addFluidProperties(), and createFluxObject().

---

The documentation for this class was generated from the following files:

• thermal_hydraulics/unit/include/TestNumericalFluxGasMixHLLC.h
• thermal_hydraulics/unit/src/TestNumericalFluxGasMixHLLC.C