[sahmed]

Dear Foamers
I have a very basic question which might be helpful for the community. I was using reactingFoam with Sutherland transport type, which calculates dynamic viscosity as a function of temperature. But in that case, could anyone tell me how does it calculate the thermal conductivity value k and where can I find that part of the code?
Instead of Sutherland, if I take type const, it assumes a constant mu and Pr and cp is calculated from JANAF and then the conductivity can easily be calculated. But how it is calculated for Sutherland type of transport?

Please help me with your comments or suggestions.

[sahmed]

Could anyone PLEASE give me any idea regarding the thermal conductivity calculation?

[DustExplosion]

Hi Sahmed,

Thermal conductivity (kappa) is calculated from viscosity in sutherlandTransportI.H. It is calculated using the Eucken approximation (See Thermal Conductivty Chapter of The Properties of Gases and Liquids, Poling, Prausnitz, and Oconnell, Page 10.2 in the Fith Edition). You may be able to get more information by searching Eucken in the forum search box (e.g., see http://www.cfd-online.com/Forums/ope...ctingfoam.html, post by carl_xw4)

template<class Thermo>
inline Foam::scalar Foam::sutherlandTransport<Thermo>::mu
(
    const scalar p,
    const scalar T
) const
{
    return As_*::sqrt(T)/(1.0 + Ts_/T);
}


template<class Thermo>
inline Foam::scalar Foam::sutherlandTransport<Thermo>::kappa
(
    const scalar p, const scalar T
) const
{
    scalar Cv_ = this->Cv(p, T);
    return mu(p, T)*Cv_*(1.32 + 1.77*this->R()/Cv_);
}


template<class Thermo>
inline Foam::scalar Foam::sutherlandTransport<Thermo>::alphah
(
    const scalar p,
    const scalar T
) const
{

    return kappa(p, T)/this->Cpv(p, T);
}