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

twoPhaseEulerFoam unphysical behaviour

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 2, 2012, 07:06
Default twoPhaseEulerFoam unphysical behaviour
  #1
Senior Member
 
Gerhard Holzinger
Join Date: Feb 2012
Location: Austria
Posts: 165
Rep Power: 14
GerhardHolzinger will become famous soon enoughGerhardHolzinger will become famous soon enough
Hello,

I am trying to simulate a bubble column. My domain is a simple box, partly filled with water.
At the center of a small patch on the bottom air is injected (alpha = 1, Ua = 0.12). Water velocity is zero at the walls and at the inlet.

What puzzles me is the fact that the water levels rises continuously. I would expect a small finite increase of the water level.

To run the case, run:

blockMesh
setFields
twoPhaseEulerFoam

The attached images are taken at T = 0 s, T = 6 s and T = 12 s.

As the amount of water should be constant the water level should not be increasing.
Attached Images
File Type: jpg t0.jpg (11.2 KB, 39 views)
File Type: jpg t6.jpg (13.9 KB, 52 views)
File Type: jpg t12.jpg (11.8 KB, 40 views)
Attached Files
File Type: gz bubbleColumnSeperated.tar.gz (4.2 KB, 21 views)

Last edited by GerhardHolzinger; July 2, 2012 at 07:08. Reason: added description of attached images
GerhardHolzinger is offline   Reply With Quote

Old   July 3, 2012, 02:28
Default update
  #2
Senior Member
 
Gerhard Holzinger
Join Date: Feb 2012
Location: Austria
Posts: 165
Rep Power: 14
GerhardHolzinger will become famous soon enoughGerhardHolzinger will become famous soon enough
I did some reading and found this in H. Rusches thesis on page 116. There he discusses ways to solve the phase continuity equation.

On page 116 he discusses an approach by Spalding. It is about adding a new equation to guarantee boundedness. This new equation is

Code:
alpha = alpha / (alpha + beta);
I am not sure if boundedness was the reason for my problems, because alpha did not exceed unity.

Adding this equation solved the problem, but I have to confess, that I do not know why exactly I do not have this problem anymore.

Code:
{
    word scheme("div(phi,alpha)");
    word schemer("div(phir,alpha)");

    surfaceScalarField phic("phic", phi);
    surfaceScalarField phir("phir", phia - phib);

    if (g0.value() > 0.0)
    {
        surfaceScalarField alphaf(fvc::interpolate(alpha));
        surfaceScalarField phipp(ppMagf*fvc::snGrad(alpha)*mesh.magSf());
        phir += phipp;
        phic += fvc::interpolate(alpha)*phipp;
    }

    for (int acorr=0; acorr<nAlphaCorr; acorr++)
    {
        fvScalarMatrix alphaEqn
        (
             fvm::ddt(alpha)
           + fvm::div(phic, alpha, scheme)
           + fvm::div(-fvc::flux(-phir, beta, schemer), alpha, schemer)
        );

        if (g0.value() > 0.0)
        {
            ppMagf = rUaAf*fvc::interpolate
            (
                (1.0/(rhoa*(alpha + scalar(0.0001))))
               *g0*min(exp(preAlphaExp*(alpha - alphaMax)), expMax)
            );

            alphaEqn -= fvm::laplacian
            (
                (fvc::interpolate(alpha) + scalar(0.0001))*ppMagf,
                alpha,
                "laplacian(alphaPpMag,alpha)"
            );
        }

        alphaEqn.relax();
        alphaEqn.solve();

        #include "packingLimiter.H"

        // addition: option 2 from H. Rusches thesis
        alpha = alpha / (alpha + beta);

        beta = scalar(1) - alpha;

        Info<< "Dispersed phase volume fraction = "
            << alpha.weightedAverage(mesh.V()).value()
            << "  Min(alpha) = " << min(alpha).value()
            << "  Max(alpha) = " << max(alpha).value()
            << endl;
    }
}

rho = alpha*rhoa + beta*rhob;
The attached images show alpha of the midplane. One can clearly see, that the water level remains where it is.
Attached Images
File Type: jpg t0.jpg (11.2 KB, 16 views)
File Type: jpg t10.jpg (12.8 KB, 19 views)
File Type: jpg t20.jpg (14.0 KB, 26 views)
GerhardHolzinger is offline   Reply With Quote

Old   October 30, 2012, 03:30
Question
  #3
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 664
Rep Power: 8
sharonyue is on a distinguished road
Quote:
Originally Posted by GerhardHolzinger View Post
I did some reading and found this in H. Rusches thesis on page 116. There he discusses ways to solve the phase continuity equation.

On page 116 he discusses an approach by Spalding. It is about adding a new equation to guarantee boundedness. This new equation is



The attached images show alpha of the midplane. One can clearly see, that the water level remains where it is.
Hi Foamers

I am doing the similar work with Gerhard.I wanna simulate a 2D tank partly filled with water. and the air was injected at the bottom.
Actually the twoPhaseEulerFoam‘s parameters are confusing me. so I download the case and change it into a 2D case. except that I didnot do anything.
Just as the picture showes. there is no bubble, but a bulk air like a squid......

how to handle that problem?

run:
fluentMeshToFoam fluent.msh
setFields
twoPhaseEulerFoam
Attached Images
File Type: png QQ??20121030151604.png (4.7 KB, 5 views)
File Type: png QQ??20121030151627.png (12.1 KB, 5 views)
File Type: png QQ??20121030151643.png (29.4 KB, 7 views)
Attached Files
File Type: zip bubbleColumnSeperated2d.zip (72.7 KB, 8 views)
sharonyue is offline   Reply With Quote

Old   November 2, 2012, 22:28
Default
  #4
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 664
Rep Power: 8
sharonyue is on a distinguished road
Quote:
Originally Posted by GerhardHolzinger View Post
I did some reading and found this in H. Rusches thesis on page 116. There he discusses ways to solve the phase continuity equation.

On page 116 he discusses an approach by Spalding. It is about adding a new equation to guarantee boundedness. This new equation is
hi GerhardHolzinger

I simulated lots of different situation and found that the rising water level is relevant to the size of inlet.

But I dont konw why. and Im not sure if this is relevant to the alpha equation.

you can see that the water level remains where it is in this group images.
but the water level rises in the last group image.

does anyone know why that happens?



btw , does the red part represent the big random air bubble?
Attached Images
File Type: jpg 1.jpg (6.9 KB, 3 views)
File Type: jpg 2.jpg (11.1 KB, 7 views)
File Type: png 3.png (36.3 KB, 6 views)

Last edited by sharonyue; November 7, 2012 at 04:13.
sharonyue is offline   Reply With Quote

Old   November 2, 2012, 22:29
Default
  #5
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 664
Rep Power: 8
sharonyue is on a distinguished road
the middle size of inlet

the water level is rising slowly
Attached Images
File Type: png QQ??20121103100907.png (2.9 KB, 5 views)
File Type: png QQ??20121103100943.png (17.7 KB, 5 views)
File Type: png QQ??20121103101005.png (21.5 KB, 3 views)
File Type: png QQ??20121103101841.png (25.4 KB, 3 views)
sharonyue is offline   Reply With Quote

Old   November 2, 2012, 22:30
Default
  #6
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 664
Rep Power: 8
sharonyue is on a distinguished road
the small inlet
level is rising a little bit fast.
Attached Images
File Type: png QQ??20121103101053.png (2.8 KB, 3 views)
File Type: png QQ??20121103101111.png (12.8 KB, 6 views)
File Type: png QQ??20121103101123.png (10.3 KB, 8 views)
File Type: png QQ??20121103101140.png (13.8 KB, 6 views)
sharonyue is offline   Reply With Quote

Old   November 5, 2012, 04:11
Default
  #7
Senior Member
 
Gerhard Holzinger
Join Date: Feb 2012
Location: Austria
Posts: 165
Rep Power: 14
GerhardHolzinger will become famous soon enoughGerhardHolzinger will become famous soon enough
Did you add the equation I mentioned?

In OpenFOAM 2.1.x (some time after 2.1.1 came out) the two-phase solvers were adapted to use MULES. Now there is no problem with free water surfaces.
GerhardHolzinger is offline   Reply With Quote

Old   November 7, 2012, 04:12
Default
  #8
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 664
Rep Power: 8
sharonyue is on a distinguished road
Quote:
Originally Posted by GerhardHolzinger View Post
Did you add the equation I mentioned?

In OpenFOAM 2.1.x (some time after 2.1.1 came out) the two-phase solvers were adapted to use MULES. Now there is no problem with free water surfaces.
yeah, I upgrade it, its fixed.

I have a little question.

Alberto said in twoPhaseEulerFoam, there would be no bubbles. but why there are bubbles in the tutorials/bed2/ and bed/?
sharonyue is offline   Reply With Quote

Old   November 7, 2012, 04:28
Default
  #9
Senior Member
 
Gerhard Holzinger
Join Date: Feb 2012
Location: Austria
Posts: 165
Rep Power: 14
GerhardHolzinger will become famous soon enoughGerhardHolzinger will become famous soon enough
I think Alberto meant, the difference between an Euler-Euler and a VOF-approach.

In Euler-Euler all phases are distributed. You do not calculate single bubbles but a gaseous phase.
So when simulating a bubble column in an Eulerian-Framework, there are no bubbles. There are only phase 'a' and 'b'. One of them is water and the other one is air.
To consider the bubbly nature of air inside water, you have to provide a bubble diameter in your case setup. This is needed e.g. for the calculation of drag.
You solve the momentum equation for both water and air in all your domain. Additionally you solve a phase continuity equation that tells you, where your water and where your air is located. The field alpha is the volume fraction of the disperse phase.
So, if you are having a bubble column, alpha=0 means there is only water and alpha=1 means there is only air.
Air bubbles rising through water mean a value of alpha between 0 and 1. So, there are no actual bubbles. There is only a volume fraction in an Eulerian simulation.


You can play the same game with solid particles and air, then you can call it fluidized bed.
GerhardHolzinger is offline   Reply With Quote

Old   November 7, 2012, 05:09
Default
  #10
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 664
Rep Power: 8
sharonyue is on a distinguished road
Well thank you for your detailed answer and the information you always provide.Gerhard Holzinger. but now I am still confused.


Quote:
Originally Posted by GerhardHolzinger View Post
Air bubbles rising through water mean a value of alpha between 0 and 1.
I dont know why the alpha is between 0 and 1 in a bubble, should it be 0?

Quote:
Originally Posted by GerhardHolzinger View Post
So when simulating a bubble column in an Eulerian-Framework, there are no bubbles.
its very weird that there are no bubbles in a bubble column, Im really lost......confused....

I want to simulate air injected in a stirred tank filled with high viscosity (about 2000cp)non-Newtonian fluid(xanthan gum), in experiment, I find there are numerous small air bubble (1mm) cause of the high viscosity,and many big bubbles breakup like boiling in the free water level.
see it here Coalescence and breakup in twoPhaseEulerFoam

So in reality , we can see lots of big bubbles and numerous small bubbles.
just as you said, can i perceive that two-fluid model can not simulate big bubble?
sharonyue is offline   Reply With Quote

Old   November 7, 2012, 05:25
Default
  #11
Senior Member
 
Gerhard Holzinger
Join Date: Feb 2012
Location: Austria
Posts: 165
Rep Power: 14
GerhardHolzinger will become famous soon enoughGerhardHolzinger will become famous soon enough
Quote:
Originally Posted by sharonyue View Post

I dont know why the alpha is between 0 and 1 in a bubble, should it be 0?

its very weird that there are no bubbles in a bubble column, Im really lost......confused....
?
What I meant with "There are no bubbles" is, that you don't see any bubbles when you look on your results in paraView.

I attached some pictures. Two of them are from paraView, showing the same image, just with a different scale on the color-bar. The third picture is an image of nearly the same situation in an experiment.

You see in both cases a bubble plume rising. In the photo you hardly see any actual bubbles because the exposure time is too long. So rising bubbles result in a somehow blurred image of the mean bubble path.

The two images out of paraView show the mean values of alpha - the volume fraction of the air. In this images there are no bubbles, only values of alpha between 0 and 1. If in a cell the value of alpha is 0, then there are no bubbles in this cell.
If in a cell the value of alpha is 0.1, then there are bubbles rising through this cell, and the overall bubble volume is 0.1 times the cell volume.
Attached Images
File Type: jpg alphaMidPlaneT10rescaledRange.jpg (14.2 KB, 36 views)
File Type: jpg alphaMidPlaneT10fullRange.jpg (12.3 KB, 29 views)
File Type: jpg Img0678.jpg (55.9 KB, 46 views)
GerhardHolzinger is offline   Reply With Quote

Old   November 7, 2012, 06:13
Smile
  #12
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 664
Rep Power: 8
sharonyue is on a distinguished road
Quote:
Originally Posted by GerhardHolzinger View Post
What I meant with "There are no bubbles" is, that you don't see any bubbles when you look on your results in paraView.
wow! thank you for that image...now I know why there is no bubble,

one last question.

in your experiment, the bubble is so small. normally if air was injected in water,there would be lots of big bubbles than yours(right?or I just make a common sense mistake?). did u add grille in the inlet?
if u dont add grille in the inlet, or enlagrge the inlet, just as the image with a small scale on the color-bar I attached.

maybe its bubbles?...

anyway,thanks for your extremely patience...I own you a favor~
Attached Images
File Type: jpg ???.jpg (11.2 KB, 12 views)
File Type: jpg ???1.jpg (13.0 KB, 16 views)
sharonyue is offline   Reply With Quote

Old   November 7, 2012, 07:52
Default
  #13
Senior Member
 
Gerhard Holzinger
Join Date: Feb 2012
Location: Austria
Posts: 165
Rep Power: 14
GerhardHolzinger will become famous soon enoughGerhardHolzinger will become famous soon enough
The bubble size when air is injected into water depends on how you inject it. In my experiment air enters through a porous plate. Therefore, the bubbles are relatively fine. If air enters through a perforated plate with 1mm holes, then the bubbles are bigger.
GerhardHolzinger is offline   Reply With Quote

Old   April 14, 2014, 21:02
Default
  #14
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 664
Rep Power: 8
sharonyue is on a distinguished road
It has been a long time,

To those who is confused to this unclear boundary of alpha in twoPhaseEulerFoam with interFoam.

Because in twoPhaseEulerFoam it uses two-fluid model. It means, in every cell, you can have air and water's velocity respectively. So these two phases are permeating with each other. Apart from this, because this two-fluid model is not surface tracking method like VOF. So it can not have a sharp phase boundary.

But why there are some "big bubbles"(in red) in this case? This is just becuase air's volume fraction is big such as 0.6 or 0.7. I mean, This is not big bubble. It can be seen as bubble cluster.

If Im wrong, please correct me.
sharonyue 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
twoPhaseEulerFoam - unphysical behaviour Lydia OpenFOAM Running, Solving & CFD 14 June 14, 2012 03:21
twoPhaseEulerFoam - weird behaviour grjmell OpenFOAM Running, Solving & CFD 5 May 22, 2012 04:07
twoPhaseEulerFoam freemankofi OpenFOAM 0 May 23, 2011 16:24
Unphysical behaviour in dynamic mesh/vof model xaero FLUENT 0 April 24, 2011 07:50
Unstable behaviour after long period of stablility plunge11 FLUENT 1 April 6, 2011 09:15


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