CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

Problem when adding temperature to interFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 2, 2012, 01:58
Default Problem when adding temperature to interFoam
  #1
Member
 
HouKen
Join Date: Jul 2011
Posts: 66
Rep Power: 6
houkensjtu is on a distinguished road
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 ...)
houkensjtu is offline   Reply With Quote

Old   November 5, 2012, 00:10
Default
  #2
Member
 
HouKen
Join Date: Jul 2011
Posts: 66
Rep Power: 6
houkensjtu is on a distinguished road
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!
houkensjtu is offline   Reply With Quote

Old   November 5, 2012, 01:29
Default
  #3
Senior Member
 
Nima Sam
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,123
Blog Entries: 1
Rep Power: 14
nimasam is on a distinguished road
Send a message via Yahoo to 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
__________________
Training Course on OpenFOAM at (http://www.isme.ir/)
My Weblog (http://openfoam.blogfa.com/)
nimasam is online now   Reply With Quote

Old   November 5, 2012, 02:06
Default
  #4
Member
 
HouKen
Join Date: Jul 2011
Posts: 66
Rep Power: 6
houkensjtu is on a distinguished road
Quote:
Originally Posted by nimasam View Post
, 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!)
houkensjtu is offline   Reply With Quote

Old   November 5, 2012, 03:35
Default
  #5
Senior Member
 
Nima Sam
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,123
Blog Entries: 1
Rep Power: 14
nimasam is on a distinguished road
Send a message via Yahoo to nimasam
if you want to see Cp
just define a volScalarField like T in "creatFields.H"
__________________
Training Course on OpenFOAM at (http://www.isme.ir/)
My Weblog (http://openfoam.blogfa.com/)
nimasam is online now   Reply With Quote

Reply

Thread Tools
Display Modes

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 Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Adding temperature equation to porousSimpleFoam David_010 OpenFOAM Programming & Development 7 February 18, 2013 09:42
InterFoam viscosity problem / possible bug pbryant OpenFOAM Running, Solving & CFD 22 October 29, 2012 04:43
adding temperature field to icofoam houkensjtu OpenFOAM 5 September 26, 2012 22:20
Interface location problem in interFoam cfd_user2011 OpenFOAM 1 July 7, 2011 12:04
[CFXv11]fluid temperature distr problem mactech001 CFX 1 June 10, 2010 08:08


All times are GMT -4. The time now is 15:30.