ideal_gas_micro_thermo.h

Go to the documentation of this file.

//-----------------------------------------------------------------------bl-
//--------------------------------------------------------------------------
//
// Antioch - A Gas Dynamics Thermochemistry Library
//
// Copyright (C) 2014-2016 Paul T. Bauman, Benjamin S. Kirk,
//                         Sylvain Plessis, Roy H. Stonger
//
// Copyright (C) 2013 The PECOS Development Team
//
// This library is free software; you can redistribute it and/or
// modify it under the terms of the Version 2.1 GNU Lesser General
// Public License as published by the Free Software Foundation.
//
// This library is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
// Lesser General Public License for more details.
//
// You should have received a copy of the GNU Lesser General Public
// License along with this library; if not, write to the Free Software
// Foundation, Inc. 51 Franklin Street, Fifth Floor,
// Boston, MA  02110-1301  USA
//
//-----------------------------------------------------------------------el-

#ifndef ANTIOCH_IDEAL_GAS_MICRO_THERMO_H
#define ANTIOCH_IDEAL_GAS_MICRO_THERMO_H

namespace Antioch
{

  template <typename CoeffType>
  class ChemicalMixture;

  template <typename CoeffType>
  class TempCache;


  template <typename MacroThermo, typename CoeffType = double>
  class IdealGasMicroThermo
  {
       public:
         IdealGasMicroThermo(const MacroThermo & ext_thermo, const ChemicalMixture<CoeffType> & chem_mix);
         ~IdealGasMicroThermo();

        template <typename StateType>
        const ANTIOCH_AUTO(StateType)
           cv_vib(unsigned int s, const StateType & T) const
         ANTIOCH_AUTOFUNC(StateType, (_chem_mix.chemical_species()[s]->n_tr_dofs() < CoeffType(2.))?zero_clone(T):
                                                                                                    _ext_therm.cv(TempCache<StateType>(T),s) - this->cv_tr(s))

        template <typename StateType>
        const ANTIOCH_AUTO(StateType)
           cv_vib_over_R(unsigned int s, const StateType & T) const
         ANTIOCH_AUTOFUNC(StateType, (_chem_mix.chemical_species()[s]->n_tr_dofs() < CoeffType(2.))?zero_clone(T):
                                                                                                    _ext_therm.cv_over_R(TempCache<StateType>(T),s) - this->cv_tr_over_R(s))

        const CoeffType cv_rot(unsigned int s) const;

        const CoeffType cv_rot_over_R(unsigned int s) const;

        const CoeffType cv_trans(unsigned int s) const;

        const CoeffType cv_trans_over_R(unsigned int s) const;

        const CoeffType cv_tr(unsigned int s) const;

        const CoeffType cv_tr_over_R(unsigned int s) const;

       private:
        const ChemicalMixture<CoeffType> & _chem_mix;
        const MacroThermo             & _ext_therm;
  };

  template <typename MacroThermo, typename CoeffType>
  inline
  IdealGasMicroThermo<MacroThermo, CoeffType>::IdealGasMicroThermo(const MacroThermo & ext_therm, const ChemicalMixture<CoeffType> & chem_mix):
    _chem_mix(chem_mix),
    _ext_therm(ext_therm)
  {
     return;
  }

  template <typename MacroThermo, typename CoeffType>
  inline
  IdealGasMicroThermo<MacroThermo, CoeffType>::~IdealGasMicroThermo()
  {
     return;
  }

  template <typename MacroThermo, typename CoeffType>
  inline
  const CoeffType IdealGasMicroThermo<MacroThermo, CoeffType>::cv_tr_over_R(unsigned int s) const
  {
     return _chem_mix.chemical_species()[s]->n_tr_dofs();
  }

  template <typename MacroThermo, typename CoeffType>
  inline
  const CoeffType IdealGasMicroThermo<MacroThermo, CoeffType>::cv_tr(unsigned int s) const
  {
    return _chem_mix.R(s) * (_chem_mix.chemical_species()[s])->n_tr_dofs();
  }

  template <typename MacroThermo, typename CoeffType>
  inline
  const CoeffType IdealGasMicroThermo<MacroThermo, CoeffType>::cv_rot(unsigned int s) const
  {
    using std::max;

    return max(this->cv_tr(s) - this->cv_trans(s), CoeffType(0) ); 
  }

  template <typename MacroThermo, typename CoeffType>
  inline
  const CoeffType IdealGasMicroThermo<MacroThermo, CoeffType>::cv_rot_over_R(unsigned int s) const
  {
    return std::max(this->cv_tr_over_R(s) - this->cv_trans_over_R(s), CoeffType(0) ); 
  }

  template <typename MacroThermo, typename CoeffType>
  inline
  const CoeffType IdealGasMicroThermo<MacroThermo, CoeffType>::cv_trans( const unsigned int species ) const
  {
    return CoeffType(1.5)*_chem_mix.R(species);
  }

  template <typename MacroThermo, typename CoeffType>
  inline
  const CoeffType IdealGasMicroThermo<MacroThermo, CoeffType>::cv_trans_over_R( const unsigned int /*species*/ ) const
  {
    return CoeffType(1.5);
  }


}


#endif