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

[swak4Foam] calculate bubble velocity

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

Like Tree3Likes
  • 2 Post By gschaider
  • 1 Post By arsalan.dryi

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 1, 2012, 13:28
Default calculate bubble velocity
  #1
Senior Member
 
Nima Samkhaniani
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,265
Blog Entries: 1
Rep Power: 23
nimasam is on a distinguished road
hi former
i guess, it should be possible to access a bubble velocity with swak4Foam
but i dont know how the procedure can be?
the procedure should be like that:

1) select all cells with alpha < 0.5
2) calculate the gravity center in each time step

any comment or suggestion?
nimasam is offline   Reply With Quote

Old   March 1, 2012, 15:10
Default
  #2
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 4,224
Rep Power: 49
gschaider will become famous soon enoughgschaider will become famous soon enough
Quote:
Originally Posted by nimasam View Post
hi former
i guess, it should be possible to access a bubble velocity with swak4Foam
but i dont know how the procedure can be?
the procedure should be like that:

1) select all cells with alpha < 0.5
2) calculate the gravity center in each time step

any comment or suggestion?
Why use a cutoff? that way you're going to get jumps in your result. Something like

"(pos()*vol()*(1-alpha1))/sum(vol()*(1-alpha1))"

(with an accumulation sum) might give you the center of the "non-fluid".Have a look at my presentation from the last workshop (you'll find it on the swak4Foam-page on the Wiki). Slide 76 has a similar application.

BTW: if you're interested in the velocity of the liquid interface then you might want to have a look at slide 155 where it s demonstrated how to calculate that with sampledSurfaces
nimasam and arsalan.dryi like this.
gschaider is online now   Reply With Quote

Old   May 24, 2016, 15:56
Default
  #3
Member
 
Arsalan
Join Date: Jul 2014
Posts: 74
Rep Power: 10
arsalan.dryi is on a distinguished road
Hi Foamers,

I'm doing a 3D simulation of two and three bubble rising using a modified interFoam solver and I need to bubbles centre position, velocity and surface area.

For a single bubble rising I used swak4Foam expressions for example for bubble centre position in Y as follows :
Code:
    bubbleCentreY
    {
        type swakExpression;
        valueType internalField;
        verbose true;
    variables (
    "Vol= sum (alpha1 < 0.5 ? vol() : 0);"
    "VolY= sum (alpha1 < 0.5 ? pos().y*vol() : 0);"    
    );        
    expression "VolY/Vol";
        accumulations (
        min        
        );

     
    }
Is there a way to compute two or three bubble properties in this manner?

Thanks in advance,
Best Regards,
Arsalan.
zeynab hoseini likes this.
arsalan.dryi is offline   Reply With Quote

Old   March 8, 2022, 04:14
Smile
  #4
Member
 
sadra mahmoudi
Join Date: Feb 2021
Location: Austria
Posts: 34
Rep Power: 3
sadra2003 is on a distinguished road
Hello all,

Thank you for the info you provided here.
I would be appreciative if you let me know your opinion. (I have attached my case)
I am also working on the terminal velocity of bubbles. Using the paraView, I measure the center of a bubble location in two successive time steps, then by deviding the displacement of the center of bubble to the time difference, I want to calculate the velocity. The problem is, the value that I gain is 50% lower than the reported values ​​in the literature. I am using a 2D mesh in openFoam 8 using interFoam.
The contents of my 0 folder are:

U file:

dimensions [0 1 -1 0 0 0 0];

internalField uniform (0 0 0);

boundaryField
{
bottom
{
type noSlip;
}
outlet
{
type noSlip;
}
walls
{
type slip;
}
defaultFaces
{
type empty;
}
}


p_rgh file:
dimensions [1 -1 -2 0 0 0 0];

internalField uniform 0;

boundaryField
{
bottom
{
type zeroGradient;
}

outlet
{
type zeroGradient;
}

walls
{
type zeroGradient;
}

defaultFaces
{
type empty;
}
}

alpha file:


dimensions [0 0 0 0 0 0 0];

internalField uniform 0;

boundaryField
{
bottom
{
type zeroGradient;
}

outlet
{
type zeroGradient;
}

walls
{
type zeroGradient;
}

defaultFaces
{
type empty;
}
}


transportProperties file:

phases (air water);

air
{
transportModel Newtonian;
nu nu [ 0 2 -1 0 0 0 0 ] 1.5E-5;
rho rho [ 1 -3 0 0 0 0 0 ] 1.18;

}

water
{
transportModel Newtonian;
nu nu [ 0 2 -1 0 0 0 0 ] 7.22E-7;
rho rho [ 1 -3 0 0 0 0 0 ] 995.7;

}

Thanks a lot.
Attached Files
File Type: zip case.zip (8.6 KB, 2 views)
sadra2003 is offline   Reply With Quote

Reply

Tags
bubble, swak4foam

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
VOF bubble drift velocity niriosrat STAR-CCM+ 1 October 18, 2016 10:01
How to calculate velocity of rising bubble in Ansys fluent using VOF model BUSHRA KHATOON Fluent Multiphase 0 July 19, 2016 02:18
Calculate particle velocity over space giack OpenFOAM Post-Processing 0 April 20, 2013 11:23
How to calculate swirl velocity drsrinivasan Main CFD Forum 0 November 23, 2012 00:25
velocity of air Robert Main CFD Forum 4 January 22, 2007 19:42


All times are GMT -4. The time now is 06:59.