calculation of thermal conductivity in reactingFoam
hello,
i am comparing simulation results in fluent to simulation results in openfoam. in my case i am working with two species which stream over a heated wall in reactingFoam. (in OF 2.0.1) now i want to set all bc's and coefficients for calculations to the same value. in fluent the calculations of cp, mu and kappa are done using a polynomial. i managed to fit the values for cp and mu in reactingFoam. this was simply done in my thermo.compressibleGasfile by changing the coeffs and the transportvalues As and Ts. but i can not find the place where i can see and change the calculation of the thermal conductivity kappa(?) in reactingFoam. could anybody help me and explain how and where thermal conductivity is calculated? thank you.... tomek 
Hey,
Now I am working with OpenFoam for almost two months and I think after some time of searching for the calculations used by reactingFoam I understood how it is managed with its subdirectories. Please correct me in my thoughts if i am wrong: When i have choosen hsPsiMixtureThermo<reactingMixture<sutherlandTrans port<specieThermo<janafThermo<perfectGas>>>>> I have to look for the calculations in this Files? So conductivity kappa is calculated in sutherlandTransport by using mu, R and Cv. Also the viscosity mu and the thermal diffusivity alpha is calculated here. The turbulent parts alphat and mut are calculated according to the turbulenceModel which is choosen (here: kOmegaSSt) and calculated in the belonging directory. In the end it is added to alphaEff and muEff. Is this right? Now I went on searching for the mass diffusivity but I could not find any calculation handling this problem. In fluent, diffusivity is calculated using Fick with ChapmanEnskog and LennardJones parameters. Is there something comparable in OpenFoam? The only thing I found was the effective Diffusivity, but this is used for the k and omega terms in the turbulencecalculation.... So please help me if you know how and where mass diffusivity is calculated, or what i have to change in my settings to achieve that. thanks... tomek 
Hi Tomek,
Do you have any update with the calculation including mass diffusivity? I'm looking into it too but to do this, it seems like we need to modify the code. :o Pang. 
hey pang,
there are no separate calculations for mass diffusivity and thermal diffusivity. so if u need them u'll have to construct a new solver which includes the terms.... tomek 
Morning Tomek,
Yes, I realised that too. There are some 'addion' codes which can be integrated into reactingFoam but most of them seem to work for the olderversionOpenFOAM. I'm using 2.0.x and hence thinking how to do it. :rolleyes: Thanks. Best regards, Pang. 
Hi all,
In laminar combustion cases, in order to calculate viscosity (mu) and thermal conductivity (alpha), we define As and Ts in the OpenFOAM thermodynamic properties format (as given in the reactingFoam tutorial). As such, mu is calculated by Sutherland Transport equation while alpha should be calculated using Eucken model based on mu. Kindly correct me if I am wrong. However, I still have a question regarding this mu and alpha calculation. If we provide the chemistry and thermo data in the directory $FOAM_RUN/chemkin, thermo data is in the form of O2 121386O 2 G 200.00 5000.00 1000.000 and seven coefficients for Janaff interpolation are given for first and second temperature ranges, respectively. This is in the form of typical Chemkin format (as given in the dieselFoam tutorial). In this case, how are both mu and alpha calculated? Hope that all the senior Foamers can advise. Thanks in advance. Regards, Carl. 
2 Attachment(s)
These are my files. Can anyone please advise? :o
Carl 
Hello,
In case anyone of you is still interested to know: The As and Ts values are fixed at 1.67212e06 and 170.672, respectively. Check chemkinLexer.L Best regards, Pang 
So what do you want to say? Don't you know why the values are the same?

regarding reactingFoam solver
greetings oll ,
i am using OF 2.2 and trying to solve my case with the reacting Foam solver the doubt in my mind is that do i need to use chemkin folder along with 0, constant and system in OF 2.2 as it already have a reaction and thermo. file in constant folder as given in counterflowflame2D. thanks in advance Regards sonu . 
All times are GMT 4. The time now is 14:56. 