TestNumericalFlux3EqnCentered.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 "TestNumericalFlux3EqnCentered.h"
#include "ADNumericalFlux3EqnCentered.h"
#include "THMIndicesVACE.h"

TEST_F(TestNumericalFlux3EqnCentered, testSymmetry) { testSymmetry(); }
TEST_F(TestNumericalFlux3EqnCentered, testConsistency) { testConsistency(); }

TestNumericalFlux3EqnCentered::TestNumericalFlux3EqnCentered() : TestNumericalFlux3EqnBase()
{
}

const NumericalFlux1D &
TestNumericalFlux3EqnCentered::createFluxObject()
{
  const std::string class_name = "ADNumericalFlux3EqnCentered";
  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);
}

std::vector<std::pair<std::vector<ADReal>, std::vector<ADReal>>>
TestNumericalFlux3EqnCentered::getPrimitiveSolutionsSymmetryTest() const
{
  std::vector<ADReal> W1(THMVACE1D::N_PRIM_VARS);
  W1[THMVACE1D::PRESSURE] = 1e5;
  W1[THMVACE1D::TEMPERATURE] = 300;
  W1[THMVACE1D::VELOCITY] = 1.5;

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

  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));

  return W_pairs;
}

std::vector<std::vector<ADReal>>
TestNumericalFlux3EqnCentered::getPrimitiveSolutionsConsistencyTest() const
{
  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;
}