TemperaturePressureFunctionFluidPropertiesTest.h

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
 
#pragma once
 
#include "MooseObjectUnitTest.h"
#include "TemperaturePressureFunctionFluidProperties.h"
#include "MooseParsedFunction.h"
 
class TemperaturePressureFunctionFluidPropertiesTest : public MooseObjectUnitTest
{
public:
  TemperaturePressureFunctionFluidPropertiesTest() : MooseObjectUnitTest("FluidPropertiesApp")
  {
    InputParameters params = _factory.getValidParams("ParsedFunction");
    params.set<std::string>(MooseBase::name_param) = "k_function";
    params.set<FEProblem *>("_fe_problem") = _fe_problem.get();
    params.set<FEProblemBase *>("_fe_problem_base") = _fe_problem.get();
    params.set<SubProblem *>("_subproblem") = _fe_problem.get();
    params.set<std::string>("expression") = "14 + 2e-2 * x + 3e-5 * y";
    params.set<std::string>(MooseBase::type_param) = "MooseParsedFunction";
    _fe_problem->addFunction("ParsedFunction", "k_function", params);
 
    params.set<std::string>("expression") = "1400 + 2.5 * x + 32e-5 * y";
    params.set<std::string>(MooseBase::name_param) = "rho_function";
    _fe_problem->addFunction("ParsedFunction", "rho_function", params);
 
    params.set<std::string>("expression") = "1e-3 + 1e-5 * x - 3e-9 * y";
    params.set<std::string>(MooseBase::name_param) = "mu_function";
    _fe_problem->addFunction("ParsedFunction", "mu_function", params);
 
    params.set<std::string>("expression") = "3000. + 3*x + 5e-4*y";
    params.set<std::string>(MooseBase::name_param) = "cp_function";
    _fe_problem->addFunction("ParsedFunction", "cp_function", params);
 
    _fe_problem->getFunction("k_function").initialSetup();
    _fe_problem->getFunction("rho_function").initialSetup();
    _fe_problem->getFunction("mu_function").initialSetup();
    _fe_problem->getFunction("cp_function").initialSetup();
 
    buildObjects();
  }
 
protected:
  void buildObjects()
  {
    InputParameters uo_params =
        _factory.getValidParams("TemperaturePressureFunctionFluidProperties");
    // Set the three functions and the specific isobaric heat capacity
    uo_params.set<FunctionName>("k") = "k_function";
    uo_params.set<FunctionName>("rho") = "rho_function";
    uo_params.set<FunctionName>("mu") = "mu_function";
    uo_params.set<Real>("cv") = 4186;
 
    _fe_problem->addUserObject("TemperaturePressureFunctionFluidProperties", "fp", uo_params);
    _fp = &_fe_problem->getUserObject<TemperaturePressureFunctionFluidProperties>("fp");
 
    InputParameters uo_params2 =
        _factory.getValidParams("TemperaturePressureFunctionFluidProperties");
    // Set the three functions and the specific heat capacity
    uo_params2.set<FunctionName>("k") = "k_function";
    uo_params2.set<FunctionName>("rho") = "rho_function";
    uo_params2.set<FunctionName>("mu") = "mu_function";
    uo_params2.set<FunctionName>("cp") = "cp_function";
    uo_params2.set<Real>("T_initial_guess") = 250;
    uo_params2.set<Real>("p_initial_guess") = 1e7;
    uo_params2.set<Real>("dT_integration_intervals") = 0.01;
 
    _fe_problem->addUserObject("TemperaturePressureFunctionFluidProperties", "fp_cp", uo_params2);
    _fp_cp = &_fe_problem->getUserObject<TemperaturePressureFunctionFluidProperties>("fp_cp");
 
    InputParameters uo_params3 =
        _factory.getValidParams("TemperaturePressureFunctionFluidProperties");
    // Set the three functions and the specific heat capacity
    uo_params3.set<FunctionName>("k") = "k_function";
    uo_params3.set<FunctionName>("rho") = "2000";
    uo_params3.set<FunctionName>("mu") = "mu_function";
    uo_params3.set<FunctionName>("cp") = "cp_function";
    uo_params3.set<Real>("T_initial_guess") = 250;
    uo_params3.set<Real>("p_initial_guess") = 1e7;
    uo_params3.set<Real>("dT_integration_intervals") = 0.001;
 
    _fe_problem->addUserObject(
        "TemperaturePressureFunctionFluidProperties", "fp_incompressible", uo_params3);
    _fp_incompressible = &_fe_problem->getUserObject<TemperaturePressureFunctionFluidProperties>(
        "fp_incompressible");
  }
 
  TemperaturePressureFunctionFluidProperties * _fp;
  TemperaturePressureFunctionFluidProperties * _fp_cp;
  TemperaturePressureFunctionFluidProperties * _fp_incompressible;
};