TestNumericalFlux3EqnHLLC Class Reference

Tests NumericalFlux3EqnHLLC. More...

#include <TestNumericalFlux3EqnHLLC.h>

Inheritance diagram for TestNumericalFlux3EqnHLLC:

Public Member Functions
	TestNumericalFlux3EqnHLLC ()

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...
const SinglePhaseFluidProperties &	getFluidPropertiesObject ()
	Builds and gets the fluid properties user object. 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_name
	Fluid properties user object name. More...
const SinglePhaseFluidProperties &	_fp
	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 NumericalFlux3EqnHLLC.

Definition at line 17 of file TestNumericalFlux3EqnHLLC.h.

Constructor & Destructor Documentation

TestNumericalFlux3EqnHLLC()

TestNumericalFlux3EqnHLLC::TestNumericalFlux3EqnHLLC

(

)

Definition at line 17 of file TestNumericalFlux3EqnHLLC.C.

                                                     : TestNumericalFlux3EqnBase()
{
}

Member Function Documentation

computeConservativeSolution()

std::vector< ADReal > TestNumericalFlux3EqnBase::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 TestNumericalFlux3EqnBase.C.

{
  return FlowModel1PhaseUtils::computeConservativeSolutionVector<true>(W, A, _fp);
}

computeFluxFromPrimitive()

std::vector< ADReal > TestNumericalFlux3EqnBase::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 TestNumericalFlux3EqnBase.C.

{
  return FlowModel1PhaseUtils::computeFluxFromPrimitive<true>(W, A, _fp);
}

createFluxObject()

const NumericalFlux1D & TestNumericalFlux3EqnHLLC::createFluxObject ( )

overrideprotectedvirtual

Creates the flux object to be tested.

Implements TestNumericalFlux1D.

Definition at line 22 of file TestNumericalFlux3EqnHLLC.C.

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

  return _fe_problem->getUserObject<NumericalFlux1D>(class_name);
}

getFluidPropertiesObject()

const SinglePhaseFluidProperties & TestNumericalFlux3EqnBase::getFluidPropertiesObject ( )

protectedinherited

Builds and gets the fluid properties user object.

Definition at line 36 of file TestNumericalFlux3EqnBase.C.

{
  const std::string class_name = "IdealGasFluidProperties";
  InputParameters params = _factory.getValidParams(class_name);
  params.set<Real>("gamma") = 1.4;
  params.set<Real>("molar_mass") = 11.640243719999999;
  _fe_problem->addUserObject(class_name, _fp_name, params);

  return _fe_problem->getUserObject<SinglePhaseFluidProperties>(_fp_name);
}

getPrimitiveSolutionsConsistencyTest()

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

overrideprotectedvirtual

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

Implements TestNumericalFlux1D.

Definition at line 65 of file TestNumericalFlux3EqnHLLC.C.

{
  std::vector<ADReal> W1(THMVACE1D::N_PRIM_VARS);
  W1[THMVACE1D::PRESSURE] = 1e5;
  W1[THMVACE1D::TEMPERATURE] = 300;
  W1[THMVACE1D::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 > > > TestNumericalFlux3EqnHLLC::getPrimitiveSolutionsSymmetryTest ( ) const

overrideprotectedvirtual

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

Implements TestNumericalFlux1D.

Definition at line 33 of file TestNumericalFlux3EqnHLLC.C.

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

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

  // sL > 0
  std::vector<ADReal> W3(THMVACE1D::N_PRIM_VARS);
  W3[THMVACE1D::PRESSURE] = 1e5;
  W3[THMVACE1D::TEMPERATURE] = 300;
  W3[THMVACE1D::VELOCITY] = 20;

  std::vector<ADReal> W4(THMVACE1D::N_PRIM_VARS);
  W4[THMVACE1D::PRESSURE] = 2e5;
  W4[THMVACE1D::TEMPERATURE] = 310;
  W4[THMVACE1D::VELOCITY] = 25;

  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

const SinglePhaseFluidProperties& TestNumericalFlux3EqnBase::_fp

protectedinherited

Fluid properties user object.

Definition at line 38 of file TestNumericalFlux3EqnBase.h.

Referenced by TestNumericalFlux3EqnBase::computeConservativeSolution(), and TestNumericalFlux3EqnBase::computeFluxFromPrimitive().

_fp_name

const UserObjectName TestNumericalFlux3EqnBase::_fp_name

protectedinherited

Fluid properties user object name.

Definition at line 36 of file TestNumericalFlux3EqnBase.h.

Referenced by TestNumericalFlux3EqnCentered::createFluxObject(), createFluxObject(), and TestNumericalFlux3EqnBase::getFluidPropertiesObject().

---

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

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