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

Foam flow through a channel

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 17, 2013, 05:55
Default Foam flow through a channel
  #1
New Member
 
Rogelio Chovet
Join Date: Apr 2013
Posts: 9
Rep Power: 4
rchovet is on a distinguished road
Good day everyone,

I've been working on a project for the last month and it's been really hard to obtain anything from the CFX software. I've read all previous threads with similar problems and tried their solutions, but no one seems to work.

Problem:

Multiphase problem:
Air (Disperse fluid) + Water (Continuous fluid) (Foam Flow)
Square channel (21x1 mm) (2D problem)
Channel length: 315 mm
Mesh element size: 0,5 mm square
Walls at the top and bottom
Symmetry at both side faces
Inlet: 2 cm/s (Air Volume Fraction: 0.7 and Water Volume Fraction: 0.3)
Outlet: Average Static pressure 1 atm

The main objective is to determinate the pressure losses along the channel. I've tried it all: steady state, transient, turbulent, laminar, changing the timestep, prolonging the outlet, openings, specified blend factor, course mesh. The problem doesn't even start it makes 2 iterations and crashes.

If anyone can give me a clue of what to do, you'll relief me from a lot of suffering.

Thanks in advance

RChovet
rchovet is offline   Reply With Quote

Old   April 17, 2013, 06:58
Default
  #2
Senior Member
 
RicochetJ's Avatar
 
Mr CFD
Join Date: Jun 2012
Location: Britain
Posts: 289
Rep Power: 6
RicochetJ is on a distinguished road
When it crashes what does the .out file say?

Have you tried double precision?
RicochetJ is offline   Reply With Quote

Old   April 17, 2013, 07:23
Default
  #3
New Member
 
Rogelio Chovet
Join Date: Apr 2013
Posts: 9
Rep Power: 4
rchovet is on a distinguished road
Thanks for the answer RichochetJ,

I've always used double precision option checked.
In the first iteration it starts to put walls at the outlet (I understand why it does that), then crashes and says:
ERROR #004100018 has occurred in subroutines FINMES
Message:
Fatal overflow in linear solver

Greetings
RCH
rchovet is offline   Reply With Quote

Old   April 17, 2013, 07:54
Default
  #4
Senior Member
 
OJ
Join Date: Apr 2012
Location: United Kindom
Posts: 475
Rep Power: 11
oj.bulmer will become famous soon enough
This is a classic case of divergence. Moreover, the approach you took seems to be like taking all sorts of medicines to see which one gives you relief. Instead, I would first diagnose the disease and then decide the course of medication!

1) It's a multiphase problem, try using a nice hex mesh for faster progress.
2) Try local timescale in steady state and try increasing its value. This should facilitate different timescales through out the domain, depending on the local courant numbers, instead of auto timescale that uses universal timescale. Try value of 5, or increase if necessary.
3) If the solution still diverges, try transient state with adaptive timestep, minimum timestep being 1e-8 s and max being say 1000 s, min max coeff being 5 and 10. Let CFX decide the proper timescale, since you never know whether the timesteps you tried were sufficient, before ruling it out. Make sure you give reasonable convergence target.

OJ
oj.bulmer is offline   Reply With Quote

Old   April 17, 2013, 07:57
Default
  #5
New Member
 
Rogelio Chovet
Join Date: Apr 2013
Posts: 9
Rep Power: 4
rchovet is on a distinguished road
Thanks a lot OJ,

I'll be back as soon as i finish the diagnose.

Greetings
RCH
rchovet is offline   Reply With Quote

Old   April 17, 2013, 08:00
Default
  #6
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,662
Rep Power: 84
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
The overflow is a FAQ: http://www.cfd-online.com/Wiki/Ansys...do_about_it.3F

This is a very high air volume fraction. Are you sure your multiphase model is valid?
ghorrocks is offline   Reply With Quote

Old   April 17, 2013, 08:03
Default
  #7
New Member
 
Rogelio Chovet
Join Date: Apr 2013
Posts: 9
Rep Power: 4
rchovet is on a distinguished road
Glenn,

Thanks for the link, i've already read it. Trust me i've tried it all.
However, how do i know if the multiphase model is valid?

Regards
RCH
rchovet is offline   Reply With Quote

Old   April 17, 2013, 08:16
Default
  #8
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,662
Rep Power: 84
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
Read the CFX documentation on the multiphase models you are using and follow up the references listed - especially the momentum transfer models. These models are only valid for specific ranges of volume fraction.
ghorrocks is offline   Reply With Quote

Old   April 17, 2013, 08:18
Default
  #9
New Member
 
Rogelio Chovet
Join Date: Apr 2013
Posts: 9
Rep Power: 4
rchovet is on a distinguished road
Got it.
I'll be back with the answer
Greetings to all
rchovet is offline   Reply With Quote

Old   May 14, 2013, 06:36
Default Still Stuck!
  #10
New Member
 
Rogelio Chovet
Join Date: Apr 2013
Posts: 9
Rep Power: 4
rchovet is on a distinguished road
Good day everyone,

So I'm back with the results from your previous recommandations.

I change the mesh... Still nothing
I tried a physical timestep and increasing its value... Nothing
Transient state with adaptative timestep... Nope

And finally re-read the whole multiphase model theory and find some interesting things:
There is no volume fraction limit... Only for some drag models (Schiller Nauman and Wen Yu) but this is not the case.
Some time convergence can be obtained by defining a minimum volume fraction bigger than the default (10e-18).... Still it did not work

The closest I've got to any result is using the multiphase homogeneous model... Because my bubbles diameter is 0.5mm so I assume is small enough to apply it. Is this correct? Anyways, it does not converges and results are far away from reality.

So any other suggestions?

Thanks in advance for any kind of idea... I'm getting desperate jejejeje.

Greetings

RCH
rchovet is offline   Reply With Quote

Old   May 14, 2013, 08:26
Default
  #11
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,662
Rep Power: 84
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
Rather than just try multiphase models at random, how you done some basic analysis to determine what is going on in your flow so you can choose an appropriate multiphase model? For instance, how much slip is there between the phases? Does each different phase require a different velocity field or can a single velocity field describe it? Or a velocity field with a defined slip? There are different multiphase models for each of these approaches.

We will worry about convergence after you have selected the correct multiphase model.
ghorrocks is offline   Reply With Quote

Old   May 14, 2013, 08:37
Default
  #12
New Member
 
Rogelio Chovet
Join Date: Apr 2013
Posts: 9
Rep Power: 4
rchovet is on a distinguished road
Thanks for the commitment glenn,
I truly appreciate it...

To answer the questions:

The horizontal foam flow should drain and a 2 mm slip liquid layer should form at the bottom while the top becomes "dryer". So there should be foam all along the conduct except at the bottom.

Because of the high air volume fraction, the air bubbles and the continuous liquid present the same velocity field. However, the velocity field of the liquid film at the bottom can present different values, depending on the total flow.

I tough that the eulerian-eulerian model should work fine.

Let me know what you think about it.

Thanks again

RCH
rchovet is offline   Reply With Quote

Old   May 14, 2013, 20:03
Default
  #13
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,662
Rep Power: 84
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
Yes, you will need to do this with a Eularian-Eularian model. But which one? MUSIG? Algebraic slip? homogenous/inhomogenous?
ghorrocks is offline   Reply With Quote

Old   May 15, 2013, 08:27
Default Foam characteristics
  #14
New Member
 
Rogelio Chovet
Join Date: Apr 2013
Posts: 9
Rep Power: 4
rchovet is on a distinguished road
Thanks again for the answer Glenn,

So I read the theory (again, jejeje) and I'm a little confuse. We already stablished that it was a Eulerian-Eulerian model. Going through the specific model I got this:
MUSIG - depends on the bubbles diameter and their capacity to change its form
ASM - depends on the slip velocity equilibrium
Inhomogeneous and Homogeneous -are the classic ones that I've been trying to use.

From the description, It seems that all of them might work. Can you provide me a little of help here?

My foam properties are:
bubble diameter: 0.5 mm
dispersed phase: Air at 25
continuos phase: Water
Surface tension: 31.5 mN.m

I'll attach a picture of the foam.

Thanks again
Regards

RCH
Attached Images
File Type: jpg Foam Picture.jpg (53.3 KB, 12 views)
rchovet is offline   Reply With Quote

Old   May 15, 2013, 08:37
Default
  #15
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 10,662
Rep Power: 84
ghorrocks has a spectacular aura aboutghorrocks has a spectacular aura aboutghorrocks has a spectacular aura about
Those photographs convince me that it is unlikely you are going to get this to work. The foam is showing some tricky rheology - the right hand image shows peaks, and the left image shows a short "beam" of foam hanging out the end of the duct. This means the foam has elastic/plastic properties for it to hold these shapes against its own weight. CFX does not have any multiphase models (that I am aware of) which can model elastic/plastic deformation. This elastic/plastic behaviour is caused by the bubbles acting like little rubber balls glued together - from the gas in the bubbles being compressed and the surface tension holding it together tightly.

So you have some choices:
* Model it anyway using CFX and you will probably be miles off as there is no appropriate physical model
* Model it using another software which has a good foam model. I have no idea if such a software exists.
* Simplify it down to some form of representative fluid with a representative viscosity, density, and possibly non-Newtonian properties to match the true multiphase flow. This will also probably be difficult to get anywhere near accurate.
* Give up and model something easier. Use experimental results to guide your designs as simulations are too difficult in this type of flow.

I had to do some foam modelling in CFX a few years ago and I chose the last option.
ghorrocks is offline   Reply With Quote

Old   May 15, 2013, 08:41
Default
  #16
New Member
 
Rogelio Chovet
Join Date: Apr 2013
Posts: 9
Rep Power: 4
rchovet is on a distinguished road
Thank you very much Glenn,

I knew it was going to be tricky... Anyways, I'll follow your suggestion and go for something easier. What really matters is the non-Newtonian properties and the pressure losses this type of fluid can creates.

Best regards
RCH
rchovet is offline   Reply With Quote

Reply

Tags
convergence problem, foam flow, multiphase flow, pressure losses

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
open channel flow with open foam windyday OpenFOAM 4 September 26, 2011 17:04
channelFoam for a 3D pipe AlmostSurelyRob OpenFOAM 3 June 24, 2011 13:06
Yacht in Open Channel Flow andreimour FLUENT 1 October 14, 2010 23:54
BlockMesh FOAM warning gaottino OpenFOAM Native Meshers: blockMesh 7 July 19, 2010 14:11
Flow within a thin channel al0Vera FLUENT 0 June 15, 2009 10:20


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