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

Nonphysical volume fractions at multiphase inlet

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 17, 2008, 05:28
Default I have a problem with a multip
  #1
Member
 
Roland
Join Date: Mar 2009
Location: Netherlands
Posts: 62
Rep Power: 8
sylvester is on a distinguished road
I have a problem with a multiphase (oil, water and air) inlet of a model of an oil/water separator.



At the inlet (on the left) enters a mix of, for example, 80% oil, 19% water and 1% air. After a couple of timesteps of running multiphaseInterFoam the volume fractions get non-physical (i.e. smaller than 0 and larger than 1), see figures, and after a while the solution explodes.





CheckMesh reports no errors or warnings.

I have perceived little difference between bounded or unbounded schemes, or first or second order schemes that I have tried. For example, one of the fvSchemes settings I have tried is:

ddtSchemes
{
default Euler;
}

gradSchemes
{
default cellLimited Gauss linear 1.0;
}

divSchemes
{
div(rho*phi,U) Gauss linearUpwind cellLimited Gauss linear 1.0;
div(phi,alpha) Gauss limitedLinear01 1.0;
div(phic,alpha) Gauss limitedLinear01 1.0;
}

laplacianSchemes
{
default Gauss linear limited 0.7;
}

interpolationSchemes
{
default linear;
}

snGradSchemes
{
default corrected;
}

fluxRequired
{
default yes;
}

The inlet is specified as fixedValue for U, alphaair, alphaoil and alphawater and zeroGradient for p.

Does anyone know what is going wrong and how to fix it?

Any help on this is highly appreciated!!

Sylvester
sylvester is offline   Reply With Quote

Old   April 17, 2008, 05:35
Default Heya, There have been bugs
  #2
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,758
Rep Power: 21
hjasak will become famous soon enough
Heya,

There have been bugs in the multi-phase interFoam solver, so please make sure you've picked up all bug fixes. Secondly, the error you describe is related to consistency between multiple alpha equations. I would recommend using upwind on div(phi,alpha) and div(phic,alpha) until you sort out the problem. Once everything is running, you can go back to better schemes.

I would guess your problem is with the inlet boundary condition on alphas. Make sure they are all consistent at the inlet and they sum up to 1.

Please keep me posted,

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

Old   April 17, 2008, 05:36
Default Woops, looking at your void fr
  #3
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,758
Rep Power: 21
hjasak will become famous soon enough
Woops, looking at your void fraction in the first picture, it says max(alpha) is 1.5!!! In fact, the SUM of alphas as any point must be 1 - you've got a problem with boundary conditions.

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

Old   April 18, 2008, 08:57
Default Hi Hrvoje, Thanks for your
  #4
Member
 
Roland
Join Date: Mar 2009
Location: Netherlands
Posts: 62
Rep Power: 8
sylvester is on a distinguished road
Hi Hrvoje,

Thanks for your help.

At the start of the run the sum of alphas is 1 everywhere. Actually, the sum always remains 1. But close to the inlet one alpha becomes larger than 1, while the other two become less than 0.

Changing the two div schemes to upwind did help, but not sufficiently. It now just takes a bit longer before alphaoil > 1.

The patches described in this thread (http://www.cfd-online.com/OpenFOAM_D.../126/5412.html) were already applied.

I'm clueless at this point.

regards,
Sylvester
sylvester is offline   Reply With Quote

Old   July 1, 2009, 10:52
Default
  #5
Senior Member
 
Kent Wardle
Join Date: Mar 2009
Location: Illinois, USA
Posts: 195
Rep Power: 10
kwardle is on a distinguished road
I realize this is a stale thread, but was wondering if there is any update on this? Were you able to resolve this issue?
kwardle 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
Volume fractions in muliphase model mange FLUENT 2 May 11, 2015 07:20
volume fractions in homogeneous flow model? Matthias Rohner CFX 1 August 26, 2007 18:42
Mole fractions from mass fractions? Keith CD-adapco 5 September 27, 2002 05:39
Changing inlet mass fractions of CH4 and O2 Prateep Chatterjee FLUENT 4 June 13, 2002 12:06
Higher order for phase volume fractions? George Bergantz CFX 2 December 11, 2000 12:55


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