|
[Sponsors] |
July 5, 2014, 18:43 |
Temperature inside domain is more than BC!
|
#1 |
Senior Member
Joe
Join Date: Feb 2012
Location: Canada
Posts: 111
Rep Power: 14 |
Hello,
I am simulating natural convection in a closed domain. My BC is 370K and 280 K. I run the case in steady laminar mood. when I check the temperature distribution inside the fluid domain, it is 450K!!! I set the initial temperature at 325K.I am pretty sure with the quality of my mesh...there si no source term Does anybody know how come I have such a this problem? Thank you in advance! |
|
July 11, 2014, 11:21 |
|
#2 |
Member
David
Join Date: Aug 2012
Posts: 48
Rep Power: 13 |
It may happen that temperature goes up the BC value at the early stage of the iteration, try to limit it to the BC value and see if the number of limited cells drops down while the case is converging
|
|
July 11, 2014, 11:31 |
|
#3 |
Senior Member
Joe
Join Date: Feb 2012
Location: Canada
Posts: 111
Rep Power: 14 |
GM_XIII,
I run the case with second order for Momentum and Energy...when I switch to first order, every thing looks alright... In second order simulation, i run the case for 200,000 iterations but there were fluctuations. Hooman |
|
July 11, 2014, 11:39 |
|
#4 |
Member
David
Join Date: Aug 2012
Posts: 48
Rep Power: 13 |
You mean fluctuations in residuals or in max temperature?
|
|
July 11, 2014, 11:40 |
|
#5 |
Senior Member
Joe
Join Date: Feb 2012
Location: Canada
Posts: 111
Rep Power: 14 |
In max temperature
|
|
July 11, 2014, 11:43 |
|
#6 |
Senior Member
Join Date: Feb 2013
Location: Germany
Posts: 200
Rep Power: 24 |
Fluctuations are characteristc for 2nd order discretization methods when gradients aren't resolved good enough by the grid. Maybe you can try to resolve your grid finer especially where large gradients occur in the flow and temperature field.
|
|
July 11, 2014, 11:44 |
|
#7 |
Member
David
Join Date: Aug 2012
Posts: 48
Rep Power: 13 |
Which model are you using for air? Have you tried underrelaxing the solution of the energy equation or density? Also check what kad said, recently i read in a topic here to use adaptative grid adaptation based on temperature gradient, but be careful not to adapt too many cells in one time as it can unstabilize your run.
|
|
July 11, 2014, 11:48 |
|
#8 |
Senior Member
Joe
Join Date: Feb 2012
Location: Canada
Posts: 111
Rep Power: 14 |
David,
I am sure about the quality of grid...I created several cases with different mesh....all resulted in fluctuations... I use steady laminar model with water... I played with URF...it was helpful regards to fluctuations but caused a new problem...the residual of continuity went up to 10e3!! When I run the case with second order transientو every thing is alright Hooman |
|
July 11, 2014, 12:03 |
|
#9 |
Senior Member
Join Date: Feb 2013
Location: Germany
Posts: 200
Rep Power: 24 |
Do you use double precision solver?
|
|
July 11, 2014, 12:05 |
|
#10 |
Senior Member
Joe
Join Date: Feb 2012
Location: Canada
Posts: 111
Rep Power: 14 |
yes I do.
yes I do. yes I do. |
|
July 11, 2014, 18:11 |
|
#11 |
Senior Member
Behrooz Jamshidi
Join Date: Apr 2013
Posts: 110
Rep Power: 13 |
Using second order discretization almost cause convergency problem, which one of its solutions is to switch to second order after complete convergence of first order. But I suppose you have done this.
High order schemes better capture physical changes. So if your physic is inherently time dependent and you are modeling it as steady state, its easier to converge with first order. According to your information that you said with transient run is ok, I think your convergency problem with your second order discretization is because of your physic. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
whats the cause of error? | immortality | OpenFOAM Running, Solving & CFD | 13 | March 24, 2021 07:15 |
Static Temperature / Opening Temperature | JulianP | CFX | 12 | April 10, 2019 18:00 |
UDF for Back-flow Temperature | G340 | Fluent UDF and Scheme Programming | 3 | August 21, 2013 04:56 |
injection problem | Mark New | FLUENT | 0 | August 4, 2013 01:30 |
is internalField(U) equivalent to zeroGradient? | immortality | OpenFOAM Running, Solving & CFD | 7 | March 29, 2013 01:27 |