# Problem when adding temperature to interFoam

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

 November 2, 2012, 01:58 Problem when adding temperature to interFoam #1 Member   HouKen Join Date: Jul 2011 Posts: 66 Rep Power: 7 Hi foamers! I tried to apply temperature calculation to interFoam according to the hot thread: Diverging result for Temperature field in interFoam And what I did was exactly same with the author "ovie" wrote in the #14 reply to the above thread. Everything compiled well, however when I run it on a two-phase flow case, temperature results near interface was higher than heat source temperature which is surely not reliable result. Also I test the modified solver to solve for single phase heat transfer calculation, although I didnot do precise validation, temperature field seems to be reasonable. Thus I think the main problem here lies in the interface. On the other hand, from a theoretical view, no matter interface properties(cp, Pr...) is linear interpolated or harmonic interpolated, interface temperature should be no higher than phase temperature surrounding it. I cant figure out what is the reason for this problem, plz help! (Also I am trying to debug it by output cp and other properties on the fly, but I am not sure how to do it elegantly, directly to the console or paraView, plz help ...)

 November 5, 2012, 00:10 #2 Member   HouKen Join Date: Jul 2011 Posts: 66 Rep Power: 7 I read the thread post by "ovie" again, and I do find that move the Cp to RHS (i.e., now the TEqn becomes: fvScalarMatrix TEqn ( fvm::ddt(rho, T) + fvm::div(rhoPhi, T) - fvm::laplacian((kappa/Cp), T) ); )help the convergence. But from a theoretical view, I believe it's mathematically equal. So I got even more confused, I think there must be sth. I misunderstood/didn't understood in the whole process. plz help!

 November 5, 2012, 01:29 #3 Senior Member   Nima Samkhaniani Join Date: Sep 2009 Location: Tehran, Iran Posts: 1,195 Blog Entries: 1 Rep Power: 16 , i guess it returns to the sensitivity of numerical procedure to calculation of convection terms (flux), those are mathematical the same on the whole domain except at interface but at interface they are a little different, and it can be accepted a new model __________________ Telegram channel (https://telegram.me/openfoam4Iranian) My Weblog (http://openfoam.blogfa.com/) Training Course on OpenFOAM at (http://www.isme.ir/)

November 5, 2012, 02:06
#4
Member

HouKen
Join Date: Jul 2011
Posts: 66
Rep Power: 7
Quote:
 Originally Posted by nimasam , i guess it returns to the sensitivity of numerical procedure to calculation of convection terms (flux), those are mathematical the same on the whole domain except at interface but at interface they are a little different, and it can be accepted a new model
Could you be a little bit more explicit on the difference of the interface.(or where can I find related information)
Basically I have only read Patankar and Versteeg's book on cfd and both these books didn't focus much on detailed numerical scheme.
(Another newbie question: I want to output Cp in order to see the space distribution of it(just like T and alpha), how can I do it? plz help!)

 November 5, 2012, 03:35 #5 Senior Member   Nima Samkhaniani Join Date: Sep 2009 Location: Tehran, Iran Posts: 1,195 Blog Entries: 1 Rep Power: 16 if you want to see Cp just define a volScalarField like T in "creatFields.H" __________________ Telegram channel (https://telegram.me/openfoam4Iranian) My Weblog (http://openfoam.blogfa.com/) Training Course on OpenFOAM at (http://www.isme.ir/)

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post David_010 OpenFOAM Programming & Development 7 February 18, 2013 09:42 pbryant OpenFOAM Running, Solving & CFD 22 October 29, 2012 04:43 houkensjtu OpenFOAM 5 September 26, 2012 22:20 cfd_user2011 OpenFOAM 1 July 7, 2011 12:04 mactech001 CFX 1 June 10, 2010 08:08

All times are GMT -4. The time now is 13:02.