318 StateType Pr = M * k0/kinf;
319 StateType dPr_dT = Pr * (dk0_dT/k0 - dkinf_dT/kinf);
320 StateType log10Pr = Constants::log10_to_log<CoeffType>() * ant_log(Pr);
321 StateType dlog10Pr_dT = Constants::log10_to_log<CoeffType>()*dPr_dT/Pr;
323 for(
unsigned int ip = 0; ip < dlog10Pr_dX.size(); ip++)
325 dlog10Pr_dX[ip] = Constants::log10_to_log<CoeffType>()/M;
331 StateType dlog10Fcent_dT = Constants::log10_to_log<CoeffType>()*dFcent_dT/Fcent;
333 StateType logFcent = ant_log(Fcent);
336 StateType c = - CoeffType(0.4L) -
_c_coeff * logFcent;
337 StateType n = CoeffType(0.75L) -
_n_coeff * logFcent;
342 StateType logF = logFcent/(1 + ant_pow(((log10Pr + c)/(n - d*(log10Pr + c) )),2));
343 StateType dlogF_dT = logF * (dlog10Fcent_dT / Fcent
344 - 2 * ant_pow((log10Pr + c)/(n - d * (log10Pr + c)),2)
345 * ((dlog10Pr_dT + dc_dT)/(log10Pr + c) -
346 (dn_dT - d * (dlog10Pr_dT + dc_dT))/(n - d * (log10Pr + c))
348 / (1 + ant_pow((log10Pr + c)/(n - d * (log10Pr + c)),2))
351 for(
unsigned int ip = 0; ip < dlog10Pr_dX.size(); ip++)
353 dlogF_dX[ip] = - ant_pow(logF,2)/logFcent * dlog10Pr_dX[ip] *(1 - 1/(n - d * (log10Pr + c))) * (log10Pr + c);
357 dF_dT =
F * dlogF_dT;
358 for(
unsigned int ip = 0; ip < dlog10Pr_dX.size(); ip++)
360 dF_dX[ip] =
F * dlogF_dX[ip];
StateType Fcent(const StateType &T) const
#define antioch_assert_equal_to(expr1, expr2)
Scalar F(const Scalar &x)
_Matrix< _Scalar, _Rows, _Cols, _Options, _MaxRows, _MaxCols > constant_clone(const _Matrix< _Scalar, _Rows, _Cols, _Options, _MaxRows, _MaxCols > &ex, const Scalar &value)
CoeffType _n_coeff
Precompute coefficient for log conversion.
CoeffType _c_coeff
Precompute coefficient for log conversion.
const ANTIOCH_AUTO(StateType) KineticsTheoryThermalConductivity< ThermoEvaluator
The parameters are reduced parameters.
_Matrix< _Scalar, _Rows, _Cols, _Options, _MaxRows, _MaxCols > zero_clone(const _Matrix< _Scalar, _Rows, _Cols, _Options, _MaxRows, _MaxCols > &ex)
void Fcent_and_derivatives(const StateType &T, StateType &Fc, StateType &dFc_dT) const