kooij_rate_unit.C File Reference

#include <limits>
#include <vector>
#include "antioch/kooij_rate.h"
#include "antioch/physical_constants.h"
#include "antioch/units.h"

Go to the source code of this file.

Functions
template<typename Scalar >
int	check_rate_and_derivative (const Scalar &rate_exact, const Scalar &derive_exact, const Scalar &rate, const Scalar &derive, const Scalar &T)
template<typename Scalar >
int	test_values (const Scalar &Cf, const Scalar &eta, const Scalar &Ea, const Scalar &Tref, const Scalar &R, const Antioch::KooijRate< Scalar > &kooij_rate)
template<typename Scalar >
int	tester ()
int	main ()

Function Documentation

template<typename Scalar >

int check_rate_and_derivative	(	const Scalar &	rate_exact,
		const Scalar &	derive_exact,
		const Scalar &	rate,
		const Scalar &	derive,
		const Scalar &	T
	)

Definition at line 36 of file kooij_rate_unit.C.

Referenced by test_values().

{
    const Scalar tol = std::numeric_limits<Scalar>::epsilon() * 2;
    int return_flag(0);
    if( abs( (rate - rate_exact)/rate_exact ) > tol )
      {
        std::cout << std::scientific << std::setprecision(16)
                  << "Error: Mismatch in rate values." << std::endl
                  << "T = " << T << " K" << std::endl
                  << "rate(T) = " << rate << std::endl
                  << "rate_exact = " << rate_exact << std::endl
                  << "relative difference = " <<  abs( (rate - rate_exact)/rate_exact ) << std::endl
                  << "tolerance = " <<  tol << std::endl;

        return_flag = 1;
      }
    if( abs( (derive - derive_exact)/derive_exact ) > tol )
      {
        std::cout << std::scientific << std::setprecision(16)
                  << "Error: Mismatch in rate derivative values." << std::endl
                  << "T = " << T << " K" << std::endl
                  << "drate_dT(T) = " << derive << std::endl
                  << "derive_exact = " << derive_exact << std::endl
                  << "relative difference = " <<  abs( (derive - derive_exact)/derive_exact ) << std::endl
                  << "tolerance = " <<  tol << std::endl;

        return_flag = 1;
      }

     return return_flag;
}

int main

(

)

Definition at line 165 of file kooij_rate_unit.C.

{
  return (tester<double>() ||
          tester<long double>() ||
          tester<float>());
}

template<typename Scalar >

int test_values	(	const Scalar &	Cf,
		const Scalar &	eta,
		const Scalar &	Ea,
		const Scalar &	Tref,
		const Scalar &	R,
		const Antioch::KooijRate< Scalar > &	kooij_rate
	)

Definition at line 70 of file kooij_rate_unit.C.

References check_rate_and_derivative(), Antioch::KooijRate< CoeffType >::derivative(), std::pow(), and Antioch::KooijRate< CoeffType >::rate_and_derivative.

Referenced by tester().

{
  using std::abs;
  using std::exp;
  using std::pow;

  int return_flag = 0;

  for(Scalar T = 300.1L; T <= 2500.1L; T += 10.L)
  {

    const Scalar rate_exact = Cf*pow(T/Tref,eta)*exp(-Ea/(R*T));
    const Scalar derive_exact = exp(-Ea/(R*T)) * pow(T/Tref,eta) * Cf * (Ea/(R*T*T) + eta/T );
    Antioch::KineticsConditions<Scalar> cond(T);

// KineticsConditions
    Scalar rate = kooij_rate(cond);
    Scalar deriveRate = kooij_rate.derivative(cond);

    return_flag = check_rate_and_derivative(rate_exact,derive_exact,rate,deriveRate,T) || return_flag;

    kooij_rate.rate_and_derivative(cond,rate,deriveRate);

    return_flag = check_rate_and_derivative(rate_exact,derive_exact,rate,deriveRate,T) || return_flag;

// T
    rate = kooij_rate(T);
    deriveRate = kooij_rate.derivative(T);

    return_flag = check_rate_and_derivative(rate_exact,derive_exact,rate,deriveRate,T) || return_flag;

    kooij_rate.rate_and_derivative(T,rate,deriveRate);

    return_flag = check_rate_and_derivative(rate_exact,derive_exact,rate,deriveRate,T) || return_flag;

  }
  return return_flag;
}

template<typename Scalar >

int tester

(

)

Definition at line 110 of file kooij_rate_unit.C.

References Antioch::Units< T >::get_SI_factor(), Antioch::KooijRate< CoeffType >::reset_coefs(), Antioch::KooijRate< CoeffType >::set_Cf(), Antioch::KooijRate< CoeffType >::set_Ea(), Antioch::KooijRate< CoeffType >::set_eta(), Antioch::KooijRate< CoeffType >::set_rscale(), Antioch::KooijRate< CoeffType >::set_Tref(), test_values(), and Antioch::KineticsModel::Tref().

{

  Scalar Cf = 1.4L;
  Scalar eta = 1.2L;
  Scalar Ea = 298.0L;
  Scalar Tref = 1.L;
  Scalar R = 1.L;

  Antioch::KooijRate<Scalar> kooij_rate(Cf,eta,Ea,Tref,R);

  int return_flag = test_values(Cf,eta,Ea,Tref,R,kooij_rate);

  Cf = 1e-7L;
  eta = 0.7L;
  Ea = 36000.0L;
  Tref = 298.;
  R = Antioch::Constants::R_universal<Scalar>() * Antioch::Units<Scalar>("cal").get_SI_factor();
  kooij_rate.set_Cf(Cf);
  kooij_rate.set_eta(eta);
  kooij_rate.set_Ea(Ea);
  kooij_rate.set_Tref(Tref);
  kooij_rate.set_rscale(R);

  return_flag = test_values(Cf,eta,Ea,Tref,R,kooij_rate) || return_flag;

  Cf = 2.1e-7L;
  eta = 1.1L;
  Ea = 23000.0L;
  std::vector<Scalar> values(3);
  values[0] = Cf;
  values[1] = eta;
  values[2] = Ea;
  kooij_rate.reset_coefs(values);

  return_flag = test_values(Cf,eta,Ea,Tref,R,kooij_rate) || return_flag;

  Cf = 2.1e-11L;
  eta = -0.01L;
  Ea = 100000.0L;
  R = Antioch::Constants::R_universal<Scalar>();
  Tref= 300.L;
  values.resize(5);
  values[0] = Cf;
  values[1] = eta;
  values[2] = Ea;
  values[3] = Tref;
  values[4] = R;
  kooij_rate.reset_coefs(values);

  return_flag = test_values(Cf,eta,Ea,Tref,R,kooij_rate) || return_flag;

  return return_flag;
}