CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Laminar field as initial state for turbulent two phase pipe flow

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

Like Tree1Likes
  • 1 Post By kjetil

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 17, 2009, 05:51
Question Laminar field as initial state for turbulent two phase pipe flow
  #1
Senior Member
 
Join Date: Mar 2009
Location: Norway
Posts: 137
Rep Power: 17
kjetil is on a distinguished road
Hello,
I am setting up a two-phase 40m pipe case, and have so far been successful using the laminar model.

I will also have to do a turbulent case, and this is where it gets interesting - say difficult. Not awfully challenging per se, having a working laminar case. So I have been trying to tune the fvSolutions and PISO settings adding correctors to pressure, and using smaller and smaller delta T's, but still - the simulation keeps exploding at about 1.01 seconds every run.

I tried to use a an exisiting field for alpha and the phase velocities Ua and Ub - using the result from laminar case @ 1.6 seconds. It seems to run nicely, but then suddenly the Courant number increases from 0.37 to 81 in two iterations, and it fails again.

Is this a good way to do it - using a laminar field as initial condition? And should it be chosen closer to 0 or could it be later? Should tune my fvSolutions differently? The laminar case is steady and stabile after about 30 seconds. I attach two photos from the laminar case - the one showing the entire pipe visualize how far the gas has reached by the 1.6 seconds.

Any suggestions to this?
Attached Images
File Type: jpg 40sec_39-5m.jpg (45.5 KB, 37 views)
File Type: jpg 1600msec_allpipe.jpg (21.0 KB, 25 views)
kjetil is offline   Reply With Quote

Old   July 18, 2009, 06:17
Default
  #2
Member
 
Nugroho Adi
Join Date: Mar 2009
Location: norway
Posts: 79
Rep Power: 17
mahaputra is on a distinguished road
Quote:
Originally Posted by kjetil View Post
Hello,
I am setting up a two-phase 40m pipe case, and have so far been successful using the laminar model.

I will also have to do a turbulent case, and this is where it gets interesting - say difficult. Not awfully challenging per se, having a working laminar case. So I have been trying to tune the fvSolutions and PISO settings adding correctors to pressure, and using smaller and smaller delta T's, but still - the simulation keeps exploding at about 1.01 seconds every run.

I tried to use a an exisiting field for alpha and the phase velocities Ua and Ub - using the result from laminar case @ 1.6 seconds. It seems to run nicely, but then suddenly the Courant number increases from 0.37 to 81 in two iterations, and it fails again.

Is this a good way to do it - using a laminar field as initial condition? And should it be chosen closer to 0 or could it be later? Should tune my fvSolutions differently? The laminar case is steady and stabile after about 30 seconds. I attach two photos from the laminar case - the one showing the entire pipe visualize how far the gas has reached by the 1.6 seconds.

Any suggestions to this?
Hei Kjetil




maybethis thesis report useful for you.Bay, M. O. (2008). "Development of Transient One-Dimensional Solver for Severe Slugging Simulation". M.Sc. Thesis, Aalborg University Esbjerg.


Hilsen,



Nugroho Adi
Stavanger
mahaputra is offline   Reply With Quote

Old   July 21, 2009, 06:26
Default
  #3
Senior Member
 
Join Date: Mar 2009
Location: Norway
Posts: 137
Rep Power: 17
kjetil is on a distinguished road
Thanks, Nugroho Adi
that was an interesting paper for sure - and I believe it may become useful to me later. In this case though, I think my problem is quite simpler than the case in that paper - and the author doesn't seem to use an initial laminar field either
kjetil is offline   Reply With Quote

Old   July 21, 2009, 10:15
Lightbulb
  #4
Senior Member
 
Join Date: Mar 2009
Location: Norway
Posts: 137
Rep Power: 17
kjetil is on a distinguished road
I have found a solution to this, and as courtesy to other users that have (or will have) a similar problem, I'll post the quick solution:

When it comes to slowly increasing Courant numbers, members on this forum have been referring to boundary conditions, and to verify that they are set properly. Though, this was not a solution to me - my case (just a pipe w/two phase flow) is also fairly simple.

Until this point I had not been able to use a turbulence model starting from T=0 without crashing in an early stage. As mentioned - the laminar worked well. However, first thing I did was to refine my mesh - so the cells are distributed more evenly.

Then - the major change I did was to alter the Divergence schemes in fvSchemes from limitedLinear or limitedLinearV to upwind. Now I am also being told that using upwind (which is a lower order scheme) initially is common. I changed back to the limitedLinear or limitedLinearV after a while (I had just commented them out when adding the upwind alternative) - the simulation is still running so I can't say 100% if it's successful - but so far it hasn't crashed again.
pconen likes this.
kjetil is offline   Reply With Quote

Reply

Tags
explode, initial, laminar, turbulent, two phase

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
TurbFoam problemlarge Co number sunnysun OpenFOAM Running, Solving & CFD 6 March 10, 2009 09:05
MRFSimpleFoam amp cyclic patches david OpenFOAM Running, Solving & CFD 36 October 21, 2008 22:55
how 2 freeze 1 phase flow field & start lagrangian KK CFX 5 February 14, 2008 17:48
IcoFoam parallel woes msrinath80 OpenFOAM Running, Solving & CFD 9 July 22, 2007 03:58
Terrible Mistake In Fluid Dynamics History Abhi Main CFD Forum 12 July 8, 2002 10:11


All times are GMT -4. The time now is 04:10.