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

Grid size, convergence issues

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 17, 2009, 05:19
Default Grid size, convergence issues
  #1
New Member
 
Franz Roman
Join Date: Jun 2009
Posts: 28
Rep Power: 17
franzdrs is on a distinguished road
Hi,
I have ran a simulation for a simple geometry, not much more than a parallelepiped with dimensions 3.6x1.2x1.1 m. The air flows in through a channel in one of the 1.2x1.1m faces and leaves through one of the 3.6x1.2m faces. Theres is a porous media between inlet and outlet. My grid has about 120000 cells. First I ran a simulation using the Laminar model and it wouldnt converge at all. Then I changed to the k-epsilon realizable model and it converges in less than 170 iterations. However, from the beginning I got the warning "turbulent viscosity limited to viscosity ratio of 1e+5" in a ever decreasing number of cells (starting at about 2000 cells) until the warning dissapears after about 80 iterations.
So I have some questions about this. Maybe somebody can help me:
1. What does it mean that the simulation doesnt converge with the laminar model and it does with the turbulent, and so fast? Shouldnt be always easier a convergence in Laminar model? Does this has to do with the grid? Is it too coarse?
2. If this warning "turbulent viscosity limited to viscosity ratio of 1e+5" dissapears before convergence, can I still trust my solution or is it flawed already?
3. I expected a different solution. I am interested in the velocity distribution at the outlet. With this simple configuration I expected a very non-uniform distribution, but it turns to be pretty good (and it does not give much chance for improvement ). So, given the above explained issues, can I trust this solution? Should the grid be much finer?
Can anyone help?

Thanks very much,
Franz

franzdrs is offline   Reply With Quote

Old   June 17, 2009, 08:17
Default
  #2
Senior Member
 
Aroon
Join Date: Apr 2009
Location: Racine WI
Posts: 148
Rep Power: 17
vishyaroon is on a distinguished road
Hello Franz,

How have you defined the porous media? do you have both the inertial and viscous co-efficients? If your flow is laminar then you'll the pressure drop co-relation is linear and so the inertial co-efficient should disappear. However if you still have defined the inertial co-efficient this might cause convergence problems (based on a couple of my experiences).

The porous media imposes a pressure drop across the region and hence the solution in the porous region is not very sensitive to grid refinement. So refining the grid may not give you a better solution than what you have now.

Again to add to this, you may want to remember that the porous media if dense enough causes enough back-pressure that the flow through it is more uniform than what it will be in the absence of porous media.
vishyaroon is offline   Reply With Quote

Old   June 17, 2009, 12:57
Default
  #3
New Member
 
Franz Roman
Join Date: Jun 2009
Posts: 28
Rep Power: 17
franzdrs is on a distinguished road
Hi Aroon,
thanks for your answer. Yes, I have defined both coefficients for the porous media. I calculated them from experimental data on pressure drop versus air speed through the medium. The plot points of deltaP vs air speed adjusted much better to a quadratic function than to a linear one, so I judged both coefficients were needed.
I know for sure that the flow is turbulent at the entrance to the system. I have an air speed of 3.6 m/s in the inlet, which is a rectangular aperture with dimensions 0.4x0.3 m. If the flow turns laminar when it flows through the porous media, I really dont know. the superficial velocity through it is about 0.1m/s.
About your last point, yes, the more the flow resistance the porous media imposes, the more uniform the flow through it will be. I thought, though, that the resistance of my porous media was quite low. At 0.1m/s (superficial velocity) the deltaP is about 6 Pa.

What about the warning "turbulent viscosity limited to viscosity ratio of 1e+5"? if it goes away should there be a problem?

Thanks a lot

Franz
franzdrs is offline   Reply With Quote

Old   June 18, 2009, 07:57
Default
  #4
Senior Member
 
Aroon
Join Date: Apr 2009
Location: Racine WI
Posts: 148
Rep Power: 17
vishyaroon is on a distinguished road
Hello Franz,

Have you checked the difference in the velocity profiles at the outlet with and without the porous media? Do you see a difference? If you see that the flow is more uniform with the porous media then I'd say that atleast the trend is correct.

For the limiter warning, do you know how many cells are being affected? Can you try plotting the turbulent viscosity and see how large this region is? If this is small then you can trust your results.
vishyaroon is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
mesh size and convergence ertan Main CFD Forum 2 April 3, 2009 16:15
Phase locked average in run time panara OpenFOAM 2 February 20, 2008 14:37
Limitation in grid size Omkar Main CFD Forum 6 November 18, 2007 00:16
grid adaptation for better convergence. co2 FLUENT 9 April 2, 2004 14:08
Grid Independent Solution Chuck Leakeas Main CFD Forum 2 May 26, 2000 11:18


All times are GMT -4. The time now is 12:53.