ConstCpPoly.cpp

Go to the documentation of this file.

/**
 *  @file ConstCpPoly.cpp
 * Declarations for the \link Cantera::SpeciesThermoInterpType SpeciesThermoInterpType \endlink object that
 * employs a constant heat capacity assumption (see \ref spthermo and
 * \link Cantera::ConstCpPoly ConstCpPoly \endlink).
 */
/*
 * $Author: hkmoffa $
 * $Revision: 279 $
 * $Date: 2009-12-05 14:08:43 -0500 (Sat, 05 Dec 2009) $
 */

// Copyright 2001  California Institute of Technology



#include "ConstCpPoly.h"

#include <cmath>

namespace Cantera {

  ConstCpPoly::ConstCpPoly() 
    :  m_t0(0.0),
       m_cp0_R(0.0),
       m_h0_R(0.0),
       m_s0_R(0.0),
       m_logt0(0.0),
       m_lowT(0.0),
       m_highT(0.0),
       m_Pref(0.0),
       m_index(0) {
  }
    
  ConstCpPoly::ConstCpPoly(int n, doublereal tlow, doublereal thigh, 
                           doublereal pref,
                           const doublereal* coeffs) :
    m_lowT      (tlow),
    m_highT     (thigh),
    m_Pref      (pref),
    m_index     (n) {
    m_t0 = coeffs[0];
    m_h0_R = coeffs[1]  / GasConstant;
    m_s0_R = coeffs[2]  / GasConstant;
    m_cp0_R = coeffs[3] / GasConstant;
    m_logt0 = log(m_t0);
  }

  ConstCpPoly::ConstCpPoly(const ConstCpPoly& b) :
    m_t0        (b.m_t0),
    m_cp0_R     (b.m_cp0_R),
    m_h0_R      (b.m_h0_R),
    m_s0_R      (b.m_s0_R),
    m_logt0     (b.m_logt0),
    m_lowT      (b.m_lowT),
    m_highT     (b.m_highT),
    m_Pref      (b.m_Pref),
    m_index     (b.m_index)
  {
  }    

  ConstCpPoly& ConstCpPoly::operator=(const ConstCpPoly& b) {
    if (&b != this) {
      m_t0    = b.m_t0;
      m_cp0_R = b.m_cp0_R;
      m_h0_R  = b.m_h0_R;
      m_s0_R  = b.m_s0_R;
      m_logt0 = b.m_logt0;
      m_lowT  = b.m_lowT;
      m_highT = b.m_highT;
      m_Pref  = b.m_Pref;
      m_index = b.m_index;
    }
    return *this;
  }
    
  ConstCpPoly::~ConstCpPoly(){}

  SpeciesThermoInterpType *
  ConstCpPoly::duplMyselfAsSpeciesThermoInterpType() const {
    ConstCpPoly* newCCP = new ConstCpPoly(*this);
    return (SpeciesThermoInterpType*) newCCP;
  }

  doublereal ConstCpPoly::minTemp() const {
    return m_lowT;
  }
  doublereal ConstCpPoly::maxTemp() const { 
    return m_highT;
  }
  doublereal ConstCpPoly::refPressure() const { 
    return m_Pref;
  }
      
  void ConstCpPoly::updateProperties(const doublereal* tt, 
                                     doublereal* cp_R,
                                     doublereal* h_RT,
                                     doublereal* s_R) const {
    double t = *tt;
    doublereal logt = log(t);
    doublereal rt = 1.0/t;
    cp_R[m_index] = m_cp0_R;
    h_RT[m_index] = rt*(m_h0_R + (t - m_t0) * m_cp0_R);
    s_R[m_index]  = m_s0_R + m_cp0_R * (logt - m_logt0);
  }

  void ConstCpPoly::updatePropertiesTemp(const doublereal temp, 
                                         doublereal* cp_R,
                                         doublereal* h_RT,
                                         doublereal* s_R) const {
    doublereal logt = log(temp);
    doublereal rt = 1.0/temp;
    cp_R[m_index] = m_cp0_R;
    h_RT[m_index] = rt*(m_h0_R + (temp - m_t0) * m_cp0_R);
    s_R[m_index]  = m_s0_R + m_cp0_R * (logt - m_logt0);
  }

  void ConstCpPoly::reportParameters(int &n, int &type,
                                     doublereal &tlow, doublereal &thigh,
                                     doublereal &pref,
                                     doublereal* const coeffs) const {
    n = m_index;
    type = CONSTANT_CP;
    tlow = m_lowT;
    thigh = m_highT;
    pref = m_Pref;
    coeffs[0] = m_t0;
    coeffs[1] = m_h0_R * GasConstant;
    coeffs[2] = m_s0_R * GasConstant;
    coeffs[3] = m_cp0_R * GasConstant;
  }

  void ConstCpPoly::modifyParameters(doublereal* coeffs) {
    m_t0 = coeffs[0];
    m_h0_R = coeffs[1]  / GasConstant;
    m_s0_R = coeffs[2]  / GasConstant;
    m_cp0_R = coeffs[3] / GasConstant;
    m_logt0 = log(m_t0);
  }

#ifdef H298MODIFY_CAPABILITY

  doublereal ConstCpPoly::reportHf298(doublereal* const h298) const {
    double temp = 298.15;
    doublereal h = GasConstant * (m_h0_R + (temp - m_t0) * m_cp0_R);
    if (h298) {
      h298[m_index] = h;
    }
    return h;
  }
  
  void ConstCpPoly::modifyOneHf298(const int k, const doublereal Hf298New) {
    if (k != m_index) return;
    doublereal hnow = reportHf298();
    doublereal delH = Hf298New - hnow;
    m_h0_R += delH / GasConstant;
  }
  
#endif

}