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

unbounded reversed flow on outlets

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 8, 2009, 16:08
Default unbounded reversed flow on outlets
  #1
Senior Member
 
Kent Wardle
Join Date: Mar 2009
Location: Illinois, USA
Posts: 201
Rep Power: 11
kwardle is on a distinguished road
All,
I have had problems with unbounded reversed flow on outlet boundaries in a number of different simulations I have tried using interFoam and interDyMFoam--particularly in rotating domains.

One example. I have been trying to set up a rotating domain problem using GGI. The general geometry is similar to the one described in another post for using MRF (http://www.cfd-online.com/Forums/ope...ost213325.html) except that there are multiple radial outlets from the rotating portion. I have tried a variety of different outlet types for U (pressureInletOutletVelocity, intletOutlet, zeroGradient, fixedValue) and pd (totalPressure, fixededValue ) on these boundaries and get pretty much the same behavior every time. The solution appears to go as normal for a few time steps and then the velocity begins to blow up--but still converge--until the magnitude of the velocity gets too high. Looking at the results in paraview, I see that it is the inward velocity on the outlets that is to blame. I have tried using very small time steps (~1e-7) and a variety of solutions schemes as well as upwinding. I have also explored different meshes, and reproduced the same behavior using sliding mesh instead of GGI. It seems this is a BC issue.

I have tried initializing from a uniformly rotating field and gotten the exact same behavior. It would seem that setting an outlet boundary condition for U of inletOutlet with an inletValue of (0 0 0) should impose a limit on reversed flow. In my experience, it has absolutely no effect which is quite confusing.

I have also worked a lot with Fluent and somewhat with Star-CD and Star-CCM+ for similar simulations and in general, openFoam seems to be much more sensitive to reversed flow and blows up very easily in these conditions. What are the commercial codes doing differently to deal with reversed flow? How can I get a solution using OpenFOAM for such cases?

Any insight would be very much appreciated. Thanks.

Last edited by kwardle; May 8, 2009 at 16:12. Reason: link to wrong post
kwardle is offline   Reply With Quote

Old   May 15, 2009, 10:20
Default
  #2
Senior Member
 
Aram Amouzandeh
Join Date: Mar 2009
Location: Vienna, Vienna, Austria
Posts: 186
Rep Power: 8
mabinty is on a distinguished road
Hi Kent!

I have the same problem as you describe even though I m working on another problem (see can t reach steady state chtMultiregionFaom). Was it already possible for you to find a solution for that??

Aram
mabinty is offline   Reply With Quote

Old   May 15, 2009, 11:31
Default
  #3
Senior Member
 
Kent Wardle
Join Date: Mar 2009
Location: Illinois, USA
Posts: 201
Rep Power: 11
kwardle is on a distinguished road
No I have not, but I really need one. I am pretty certain that this is a common complaint as I have heard it from others--it seems that OF is very sensitive to reversed flow--and I really hope that someone out there can (and is willing to) address this.
kwardle is offline   Reply With Quote

Old   May 16, 2009, 02:12
Default
  #4
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,763
Rep Power: 21
hjasak will become famous soon enough
Why don't you slap an inletOutlet boundary condition on U and kill the back-flow; otherwise, use a total pressure boundary condition at the outlet, which will stabilise things for you.

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   May 16, 2009, 03:16
Default
  #5
Senior Member
 
Aram Amouzandeh
Join Date: Mar 2009
Location: Vienna, Vienna, Austria
Posts: 186
Rep Power: 8
mabinty is on a distinguished road
Hi!

What my case concerns I have an "open" boundary where hot air exits at the top and ambient air at the lower part of the boundary. Hence I prescribed pressureDirectedInletOutletVelocity for U and experimented with the BCs for p (calculated, fixedValue, outletInlet, totalPressure) and pd (fixedFluxBouyandPressure, zeroGradient). But everytime the in-comming velocity blows up; beside this also the values for k and epsilon (standard k-eps) constantly increase.

Thx in advance for any comments!
Aram
mabinty is offline   Reply With Quote

Old   May 18, 2009, 09:38
Default inletOutlet
  #6
Senior Member
 
Kent Wardle
Join Date: Mar 2009
Location: Illinois, USA
Posts: 201
Rep Power: 11
kwardle is on a distinguished road
Hrv,
Thanks for the suggestions. Unfortunately, as I mentioned in the original post I have tried this and for some reason it does not work. If the inletOutlet BC is used with an inletValue of (0 0 0) then there should not be any inflow, right? Any inward directed flow would be stopped. For some strange reason this does not work for one case I am working on. In my simplified version the inletOutlet BC works fine, but in my more complex geometry the inflow still happens and blows things up. I have tried to use a totalPressure BC in combination with a pressureInletOutletVelocity and gotten the same problem. I thought maybe it was a mesh problem, but the mesh checks out with no errors and I have also tried from scratch with other meshes (the original one is converted from a Gambit .msh)

It may also be that this is a problem with using GGI for the rotating region--but as you said after the webcast last Thurs, this shouldn't be the problem. Even so, my geometry does not require GGI and I could get by with MRF but unfortunately I don't have a working MRFinterFoam solver... I get the smae behavior with sliding mesh too.

Note the the outlets (there are two separate patches each with four distinct channels) are all oriented normal to the rotation axis. Not sure if this complicates things.

What really baffles me is that I get the same behavior if I set the rotation speed to 0.1, 1, 10, 100 rpm. (I actually need the solution for > 3000 rpm). This makes me think it is a solution convergence problem, but everything is converging fine--just to a blown up solution. I tried changing a few things in fvSolution, but to no avail.

Any help would be very much appreciated.
THANKS!
kwardle is offline   Reply With Quote

Old   May 18, 2009, 10:47
Default
  #7
Senior Member
 
Aram Amouzandeh
Join Date: Mar 2009
Location: Vienna, Vienna, Austria
Posts: 186
Rep Power: 8
mabinty is on a distinguished road
Hi!

Finally I was able to reach a stable calculation (see can t reach steady state chtMultiregionFaom).

Best regards,
Aram
mabinty is offline   Reply With Quote

Old   May 18, 2009, 15:23
Default
  #8
Senior Member
 
Kent Wardle
Join Date: Mar 2009
Location: Illinois, USA
Posts: 201
Rep Power: 11
kwardle is on a distinguished road
Hrv,
I looked closer at the results and with inletOutlet BCs what I am actually getting is oscillatory behavior the flow is out at first and then in, then out, etc. -- oscillating and growing with each time step after the first 5 or so which seem stable. I also tried to use the fixedMeanValue BC that Aram mentions worked for him--same behavior. Something must be wrong with my case but I can't seem to find it.
kwardle is offline   Reply With Quote

Reply

Tags
interdymfoam, interfoam, pressure bc, rotation

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
Reversed flow at pressure outlet Seeker Phil FLUENT 9 January 2, 2010 06:21
reversed flow Lakshmi Kanth FLUENT 1 September 26, 2005 00:41
Skewness -> Reversed Flow krishna FLUENT 1 May 27, 2004 08:49
reversed flow mukkarum FLUENT 1 April 13, 2004 11:17
Error Message: Reversed flow Mnu & Luc FLUENT 3 February 5, 2003 08:02


All times are GMT -4. The time now is 17:49.