#include <IdealRealGasMixtureFluidPropertiesTest.h>

Inheritance diagram for IdealRealGasMixtureFluidPropertiesTest:

Public Member Functions
	IdealRealGasMixtureFluidPropertiesTest ()

Protected Member Functions
void	buildObjects ()

T &	addObject (const std::string &type, const std::string &name, InputParameters &params)

Protected Attributes
const StiffenedGasFluidProperties *	_fp_steam

const IdealGasFluidProperties *	_fp_nitrogen

const IdealRealGasMixtureFluidProperties *	_fp_mix

std::unique_ptr< MooseMesh >	_mesh

std::shared_ptr< MooseApp >	_app

Factory &	_factory

std::shared_ptr< FEProblem >	_fe_problem

Detailed Description

Definition at line 17 of file IdealRealGasMixtureFluidPropertiesTest.h.

Constructor & Destructor Documentation

◆ IdealRealGasMixtureFluidPropertiesTest()

IdealRealGasMixtureFluidPropertiesTest::IdealRealGasMixtureFluidPropertiesTest ( )

inline

Definition at line 20 of file IdealRealGasMixtureFluidPropertiesTest.h.

                                            : MooseObjectUnitTest("FluidPropertiesApp")
   {
     buildObjects();
   }

Member Function Documentation

◆ buildObjects()

void IdealRealGasMixtureFluidPropertiesTest::buildObjects ( )

inlineprotected

Definition at line 26 of file IdealRealGasMixtureFluidPropertiesTest.h.

Referenced by IdealRealGasMixtureFluidPropertiesTest().

   {
     const std::string fp_steam_name = "fp_steam";
     const std::string fp_nitrogen_name = "fp_nitrogen";
     const std::string fp_mix_name = "fp_mix";
 
     // steam; parameters correspond to T in range 298 K to 473 K
     {
       const std::string class_name = "StiffenedGasFluidProperties";
       InputParameters params = _factory.getValidParams(class_name);
       params.set<Real>("gamma") = 1.43;
       params.set<Real>("q") = 2030e3;
       params.set<Real>("q_prime") = -23e3;
       // This parameter is supposed to be 0 for steam, but then de/dp = 0. For
       // the used mixture model, the derivative dp/de is computed as (de/dp)^-1,
       // so this is artificially changed to a non-zero number here.
       params.set<Real>("p_inf") = 1e4;
       params.set<Real>("cv") = 1040;
       params.set<Real>("M") = 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
       params.set<Real>("T_c") = 647.096;
       params.set<Real>("rho_c") = 322.;
       params.set<Real>("e_c") = 2015734.52419;
       _fe_problem->addUserObject(class_name, fp_steam_name, params);
       _fp_steam = &_fe_problem->getUserObject<StiffenedGasFluidProperties>(fp_steam_name);
     }
 
     // nitrogen
     {
       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);
       _fp_nitrogen = &_fe_problem->getUserObject<IdealGasFluidProperties>(fp_nitrogen_name);
     }
 
     // mixture
     {
       const std::string class_name = "IdealRealGasMixtureFluidProperties";
       InputParameters params = _factory.getValidParams(class_name);
       params.set<UserObjectName>("fp_primary") = fp_steam_name;
       params.set<std::vector<UserObjectName>>("fp_secondary") = {fp_nitrogen_name};
       params.set<bool>("allow_imperfect_jacobians") = true;
       _fe_problem->addUserObject(class_name, fp_mix_name, params);
       _fp_mix = &_fe_problem->getUserObject<IdealRealGasMixtureFluidProperties>(fp_mix_name);
     }
   }