# How thermal conductivity is calculated in Sutherland transport type in reactingFoam

 Register Blogs Members List Search Today's Posts Mark Forums Read

 September 26, 2016, 17:20 How thermal conductivity is calculated in Sutherland transport type in reactingFoam #1 Member   Sheikh Ahmed Join Date: Dec 2015 Location: South Carolina, USA Posts: 88 Rep Power: 10 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. anthony761 likes this.

 September 27, 2016, 09:57 #2 Member   Sheikh Ahmed Join Date: Dec 2015 Location: South Carolina, USA Posts: 88 Rep Power: 10 Could anyone PLEASE give me any idea regarding the thermal conductivity calculation?

 October 16, 2016, 18:04 #3 Member   Chris Cloney Join Date: Jun 2016 Location: Halifax, Canada Posts: 62 Rep Power: 10 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) Code: ```template inline Foam::scalar Foam::sutherlandTransport::mu ( const scalar p, const scalar T ) const { return As_*::sqrt(T)/(1.0 + Ts_/T); } template inline Foam::scalar Foam::sutherlandTransport::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 inline Foam::scalar Foam::sutherlandTransport::alphah ( const scalar p, const scalar T ) const { return kappa(p, T)/this->Cpv(p, T); }``` hofbeer, anthony761, Tobermory and 1 others like this.

 Tags reactingfoam, sutherland, thermal conductivity, viscosity