TestNumericalFlux1D Class Reference

Base class for testing NumericalFlux1D objects. More...

#include <TestNumericalFlux1D.h>

Inheritance diagram for TestNumericalFlux1D:

Public Member Functions
	TestNumericalFlux1D ()

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

Protected Attributes
std::unique_ptr< MooseMesh >	_mesh
std::shared_ptr< MooseApp >	_app
Factory &	_factory
std::shared_ptr< FEProblem >	_fe_problem

Detailed Description

Base class for testing NumericalFlux1D objects.

Definition at line 18 of file TestNumericalFlux1D.h.

Constructor & Destructor Documentation

TestNumericalFlux1D()

TestNumericalFlux1D::TestNumericalFlux1D

(

)

Definition at line 13 of file TestNumericalFlux1D.C.

: MooseObjectUnitTest("ThermalHydraulicsApp") {}

Member Function Documentation

computeConservativeSolution()

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

protectedpure virtual

Computes the conservative solution from the primitive solution.

Parameters

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

Implemented in TestNumericalFlux3EqnBase, and TestNumericalFluxGasMixBase.

Referenced by testConsistency(), and testSymmetry().

computeFluxFromPrimitive()

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

protectedpure virtual

Computes the 1D flux vector from the primitive solution.

Parameters

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

Implemented in TestNumericalFlux3EqnBase, and TestNumericalFluxGasMixBase.

Referenced by testConsistency().

createFluxObject()

virtual const NumericalFlux1D & TestNumericalFlux1D::createFluxObject ( )

protectedpure virtual

Creates the flux object to be tested.

Implemented in TestNumericalFlux3EqnCentered, TestNumericalFlux3EqnHLLC, and TestNumericalFluxGasMixHLLC.

Referenced by testConsistency(), and testSymmetry().

getPrimitiveSolutionsConsistencyTest()

virtual std::vector< std::vector< ADReal > > TestNumericalFlux1D::getPrimitiveSolutionsConsistencyTest ( ) const

protectedpure virtual

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

Implemented in TestNumericalFlux3EqnCentered, TestNumericalFlux3EqnHLLC, and TestNumericalFluxGasMixHLLC.

Referenced by testConsistency().

getPrimitiveSolutionsSymmetryTest()

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

protectedpure virtual

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

Implemented in TestNumericalFlux3EqnCentered, TestNumericalFlux3EqnHLLC, and TestNumericalFluxGasMixHLLC.

Referenced by testSymmetry().

testConsistency()

void TestNumericalFlux1D::testConsistency ( )

protected

Runs the consistency test(s)

Definition at line 57 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 ( )

protected

Runs the symmetry test(s)

Definition at line 16 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());
}

---

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

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