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

When to stop the iteration?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 13, 2013, 03:22
Red face When to stop the iteration?
  #1
New Member
 
House xu
Join Date: Aug 2010
Posts: 10
Rep Power: 5
oldisbest is on a distinguished road
Hello,
I have one similified model for micro-channel heat exchanger. Please find the attachment of detailed discription.

only about 0.8 million meshes but quality is bad (several sharp angles near tangency area). Although I'm proficient in Gambit and know the skills to improve the volume mesh quality, I consider my problem as a common case and can be calculated without diverging issue.

So for saving time, I gave up to improve my meshes in all 6 options(1 baseline and 5 modified).

I choose SST-KW with low-RE correlation, and reduce relaxing factors. and what's important, I deactive energy equation for first-round iteration, and all 6 option got perfact results for flow and turbulence equations, converged at about 500 steps.
Then in second-round iteration, I deactive flow and turbulence while active energy equation. As usual, I reduce the convergence residual of energy equation to 1e-8, and add a monitor for mass-averaged temperature at outlet.

The monitored temperature seems change quickly in first 5-10 steps while became almost flat when residual is about 2e-7, then, every 5-6 step, the temperature will still increased by 0.01 degree C.

So I try to resuming the iteration, the residual curve of energy start rising. and then, temperature is out of limitation, that's a signal of divergence.

I found this problem happened in all the 6 options, since the final difference of monitored temperature of each option is smaller then 0.3-0.5 degree C. So in my simulation the results will influenced by whole iteration steps, I am afraid it can't reflect the actual phenomenon.
Attached Images
File Type: jpg my_model.JPG (44.0 KB, 18 views)
oldisbest is offline   Reply With Quote

Old   May 13, 2013, 03:49
Default
  #2
Senior Member
 
Join Date: Aug 2011
Posts: 313
Rep Power: 10
blackmask will become famous soon enough
There is no option other than improving your mesh. What's the temperate difference between the inlet and outlet? What is the Reynolds number based on channel width? What are the discretization schemes? You should try structured mesh for simple geometry.
blackmask is offline   Reply With Quote

Old   May 13, 2013, 04:08
Default Re
  #3
New Member
 
House xu
Join Date: Aug 2010
Posts: 10
Rep Power: 5
oldisbest is on a distinguished road
Reynolds number is about 250.
Temperature of air inlet is 35 deg C.
Temperature of heating wall is 50 deg C.

Air outlet temperature based on current CFD results is about 47 deg C.

SIMPLE method and second-order up-stream discretization is adopted except pressure (Standard)

Quote:
Originally Posted by blackmask View Post
There is no option other than improving your mesh. What's the temperate difference between the inlet and outlet? What is the Reynolds number based on channel width? What are the discretization schemes? You should try structured mesh for simple geometry.
oldisbest is offline   Reply With Quote

Old   May 13, 2013, 05:17
Default
  #4
Senior Member
 
Join Date: Aug 2011
Posts: 313
Rep Power: 10
blackmask will become famous soon enough
The flow should be laminar. I did not compute but a 0.3-0.5K temperate different at the different results in large log mean temperature difference, right? The pressure difference is a much more useful indicator for mesh independence. Please find the pressure difference for all six cases.
blackmask is offline   Reply With Quote

Old   May 13, 2013, 07:35
Default
  #5
Senior Member
 
OJ
Join Date: Apr 2012
Location: United Kindom
Posts: 462
Rep Power: 10
oj.bulmer will become famous soon enough
You have a conjugate heat transfer, along with not-so-simple flow passage- why do you say you know the simulation is a simple case? Shouldn't the mesh independence study, along with the quality improvement, be on the top of your to-do list?

OJ
oj.bulmer is offline   Reply With Quote

Old   May 20, 2013, 00:25
Default
  #6
New Member
 
House xu
Join Date: Aug 2010
Posts: 10
Rep Power: 5
oldisbest is on a distinguished road
Quote:
Originally Posted by blackmask View Post
The flow should be laminar. I did not compute but a 0.3-0.5K temperate different at the different results in large log mean temperature difference, right? The pressure difference is a much more useful indicator for mesh independence. Please find the pressure difference for all six cases.
I saw some persons were using standard k-e model in similar phisical model. and they got the results with 10-15% error (compared with test).

Yes, to adopt laminate flow model with viscous heating, the calculation can be easily converged, and I found the results of flow field is very close to the one with previous sst-kw model.

anyway, I will use laminate flow model later in this kind of problem since its calculation effiency.

To stop the calculation when considering heat transfer I now set the energy residue to default 1e-6.
oldisbest is offline   Reply With Quote

Old   May 20, 2013, 00:34
Default
  #7
New Member
 
House xu
Join Date: Aug 2010
Posts: 10
Rep Power: 5
oldisbest is on a distinguished road
Quote:
Originally Posted by oj.bulmer View Post
You have a conjugate heat transfer, along with not-so-simple flow passage- why do you say you know the simulation is a simple case? Shouldn't the mesh independence study, along with the quality improvement, be on the top of your to-do list?

OJ
It's a regular case with less than 1 million meshes, incompressible fluid, constant material properties, only conductive and convetive heat transfer involved. In my work, it's really a simple case.

Mesh sensitivity and quality is really important for a new simulation, thanks for your suggestion.
oldisbest is offline   Reply With Quote

Old   May 20, 2013, 00:50
Default
  #8
New Member
 
House xu
Join Date: Aug 2010
Posts: 10
Rep Power: 5
oldisbest is on a distinguished road
My second-round iteration will do some modifications on the model.
I will consider the tube-inside refrigerant flow which is two-phase.
So their will be two kinds of fluid, one is air, the other is two-phase refrigerant(FREON).

Do you know how to do the definition in FLUENT?
I found if two-phase model is activated, the material option in fluid cell boundary condition would be missing. That means the outside air can not be defined.

Can I define it as three-phase flow with air, liquid refrigerant, gas refrigerant.

Quote:
Originally Posted by oldisbest View Post
Hello,
I have one similified model for micro-channel heat exchanger. Please find the attachment of detailed discription.

only about 0.8 million meshes but quality is bad (several sharp angles near tangency area). Although I'm proficient in Gambit and know the skills to improve the volume mesh quality, I consider my problem as a common case and can be calculated without diverging issue.

So for saving time, I gave up to improve my meshes in all 6 options(1 baseline and 5 modified).

I choose SST-KW with low-RE correlation, and reduce relaxing factors. and what's important, I deactive energy equation for first-round iteration, and all 6 option got perfact results for flow and turbulence equations, converged at about 500 steps.
Then in second-round iteration, I deactive flow and turbulence while active energy equation. As usual, I reduce the convergence residual of energy equation to 1e-8, and add a monitor for mass-averaged temperature at outlet.

The monitored temperature seems change quickly in first 5-10 steps while became almost flat when residual is about 2e-7, then, every 5-6 step, the temperature will still increased by 0.01 degree C.

So I try to resuming the iteration, the residual curve of energy start rising. and then, temperature is out of limitation, that's a signal of divergence.

I found this problem happened in all the 6 options, since the final difference of monitored temperature of each option is smaller then 0.3-0.5 degree C. So in my simulation the results will influenced by whole iteration steps, I am afraid it can't reflect the actual phenomenon.
oldisbest is offline   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
How To Stop The Iteration Process? sefde FLUENT 6 July 25, 2012 02:18
Stop iteration and continuation of the process martapajon OpenFOAM Running, Solving & CFD 3 November 5, 2007 15:06
Parallel runs slower with MTU=9000 than MTU=1500 Javier Larrondo FLUENT 0 October 28, 2007 23:30
How to stop iteration with TUI? sophie FLUENT 4 November 23, 2006 05:18
Heat exchanger problem chiseung FLUENT 16 October 20, 2001 04:36


All times are GMT -4. The time now is 07:20.