temperature_response_functions.h

Go to the documentation of this file.

#ifndef TEMPERATURE_RESPONSE_FUNCTIONS_H
#define TEMPERATURE_RESPONSE_FUNCTIONS_H
 
#include <cmath>                     // for exp, pow
#include "../framework/constants.h"  // for ideal_gas_constant
 
inline double arrhenius_exponential(
    double c,                  // dimensionless
    double activation_energy,  // J / mol
    double temperature_k       // Kelvin
)
{
    using physical_constants::ideal_gas_constant;  // J / k / mol
    return exp(c - activation_energy / (ideal_gas_constant * temperature_k));
}
 
inline double Q10_temperature_response(
    double temperature,  // degrees C
    double Tref          // degrees C
)
{
    constexpr double Q10 = 2.0;
    return pow(Q10, (temperature - Tref) / 10.0);
}
 
inline double johnson_eyring_williams_response(
    double c,             // dimensionless
    double Ha,            // J / mol
    double Hd,            // J / mol
    double S,             // J / K / mol
    double temperature_k  // K
)
{
    using physical_constants::ideal_gas_constant;  // J / k / mol
 
    double const top =
        temperature_k * arrhenius_exponential(c, Ha, temperature_k);
 
    double const bot =
        1.0 + arrhenius_exponential(S / ideal_gas_constant, Hd, temperature_k);
 
    return top / bot;
}
 
inline double polynomial_response(
    double c0,            // output units
    double c1,            // output units * (degrees C)^(-1)
    double c2,            // output units * (degrees C)^(-2)
    double temperature_c  // degrees C
)
{
    return c0 + c1 * temperature_c + c2 * pow(temperature_c, 2);
}
 
#endif