NASA7Polynomials.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 "MooseTypes.h"
#include "metaphysicl/raw_type.h"
 
#include <array>
#include <cmath>
#include <string>
 
namespace NASA7
{
using Coefficients = std::array<Real, 7>;
 
struct SpeciesData
{
  std::string name;
  Real molar_mass;
  Real T_min;
  Real T_mid;
  Real T_max;
  Coefficients low;
  Coefficients high;
};
 
constexpr Real reference_pressure = 101325.0;
 
const SpeciesData & speciesData(const std::string & name);
 
std::string speciesNames();
 
template <typename CppType>
const Coefficients &
coefficients(const SpeciesData & data, const CppType & T)
{
  return MetaPhysicL::raw_value(T) <= data.T_mid ? data.low : data.high;
}
 
template <typename CppType>
CppType
cpOverR(const SpeciesData & data, const CppType & T)
{
  const auto & a = coefficients(data, T);
  return a[0] + T * (a[1] + T * (a[2] + T * (a[3] + T * a[4])));
}
 
template <typename CppType>
CppType
hOverRT(const SpeciesData & data, const CppType & T)
{
  const auto & a = coefficients(data, T);
  return a[0] + T * (a[1] / 2 + T * (a[2] / 3 + T * (a[3] / 4 + T * a[4] / 5))) + a[5] / T;
}
 
template <typename CppType>
CppType
sOverR(const SpeciesData & data, const CppType & T)
{
  using std::log;
  const auto & a = coefficients(data, T);
  return a[0] * log(T) + T * (a[1] + T * (a[2] / 2 + T * (a[3] / 3 + T * a[4] / 4))) + a[6];
}
}