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

Getting a concentration field around a bubble in InterFoam

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

Like Tree1Likes
  • 1 Post By azman

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 22, 2008, 14:51
Default Hi all, I have been playin
  #1
New Member
 
Azman
Join Date: Mar 2009
Location: Aachen, Germany
Posts: 3
Rep Power: 17
azman is on a distinguished road
Hi all,

I have been playing around with the InterFOAM solver for about 2 weeks now to simulate a rising bubble (my bubbles rise! all thanks to Sebastian Gatze's contribution to this Forum), and was toying with the idea of solving the concentration field of oxygen bubble rising in water using InterFOAM. A OpenFoam user has been sharing his work on rising bubbles in another thread using InterFOAM, and I was thinking whether it makes sense to get a concentration field around a bubble.

The first strategy that came to my mind was ..just adding the scalarTransport equation (just like how it was defined in the ScalarTransportFoam) in my InterFoam solver, after the gamma equation and velocity equation is solved.

solve
(
fvm::ddt(C)
+ fvm::div(U, C)
- fvm::laplacian(DT, C)
);
with DT = gamma*D02_liquid + (1-gamma)*DO2_air

I am however not sure how I can include a jump condition at the interface, such that CO2_liquidinterface = HenrysKoeff*C02_gasinterface. I fear that this is not possible with the VOF method since VOF doesnt track the interface per se, but just the vol. fraction of each cell, and magically reconstruct the interface. Am I right?

I recently came across a paper from Bothe et al., about direct numerical simulation of mass transfer between rising bubbles and the surrounding liquid that can be accessed through the link below.

http://chemie.uni-paderborn.de/filea...ecke/Literatur /bubblyflows.pdf

In their work with a self-built fvm code, the VOF method was used too, with the PLIC method to reconstruct the interface. The scalar transport equation was solved as follows:

Inside each phases:
dC'/dt + div(C'*u) = DT.grad(C')

whereby C' = C_liq at the liquid phase
C' = C_gas/H at the gas phase
DT (as above)

at the interphase:
C'_L = C'_G
D02_liquid*grad(C').n = H*D02_Gas*grad(C').n

Does anyone know what to include in the solver, such that the condition at the interphase is fulfilled.
raj kumar saini likes this.
azman is offline   Reply With Quote

Old   March 16, 2012, 07:26
Arrow Mass transfer + Interfoam
  #2
Member
 
Hanniel Freitas
Join Date: Jan 2012
Location: Natal, Brazil
Posts: 34
Rep Power: 14
hfsf is on a distinguished road
Hi, azman.
Could you send me your case of mass transfer coupled with interFoam?
My email is: hanniel88(at)gmail(dot)com
Replace (at) for @
And (dot) for .
hfsf is offline   Reply With Quote

Old   January 7, 2013, 13:47
Default
  #3
Senior Member
 
akidess's Avatar
 
Anton Kidess
Join Date: May 2009
Location: Germany
Posts: 1,377
Rep Power: 30
akidess will become famous soon enough
For future reference: http://www.cfd-online.com/Forums/ope...one-phase.html
__________________
*On twitter @akidTwit
*Spend as much time formulating your questions as you expect people to spend on their answer.
akidess is offline   Reply With Quote

Old   June 7, 2022, 04:21
Smile
  #4
Member
 
sadra mahmoudi
Join Date: Feb 2021
Location: Austria
Posts: 39
Rep Power: 5
sadra2003 is on a distinguished road
Quote:
Originally Posted by azman View Post
Hi all,

I have been playing around with the InterFOAM solver for about 2 weeks now to simulate a rising bubble (my bubbles rise! all thanks to Sebastian Gatze's contribution to this Forum), and was toying with the idea of solving the concentration field of oxygen bubble rising in water using InterFOAM. A OpenFoam user has been sharing his work on rising bubbles in another thread using InterFOAM, and I was thinking whether it makes sense to get a concentration field around a bubble.

The first strategy that came to my mind was ..just adding the scalarTransport equation (just like how it was defined in the ScalarTransportFoam) in my InterFoam solver, after the gamma equation and velocity equation is solved.

solve
(
fvm::ddt(C)
+ fvm::div(U, C)
- fvm::laplacian(DT, C)
);
with DT = gamma*D02_liquid + (1-gamma)*DO2_air

I am however not sure how I can include a jump condition at the interface, such that CO2_liquidinterface = HenrysKoeff*C02_gasinterface. I fear that this is not possible with the VOF method since VOF doesnt track the interface per se, but just the vol. fraction of each cell, and magically reconstruct the interface. Am I right?

I recently came across a paper from Bothe et al., about direct numerical simulation of mass transfer between rising bubbles and the surrounding liquid that can be accessed through the link below.

http://chemie.uni-paderborn.de/filea...ecke/Literatur /bubblyflows.pdf

In their work with a self-built fvm code, the VOF method was used too, with the PLIC method to reconstruct the interface. The scalar transport equation was solved as follows:

Inside each phases:
dC'/dt + div(C'*u) = DT.grad(C')

whereby C' = C_liq at the liquid phase
C' = C_gas/H at the gas phase
DT (as above)

at the interphase:
C'_L = C'_G
D02_liquid*grad(C').n = H*D02_Gas*grad(C').n

Does anyone know what to include in the solver, such that the condition at the interphase is fulfilled.



Hello azman,

I hope you are doing well. I would have a similar question.
I am trying to simulate a single bubble movement in a solution of water and sugar with interIsoFoam solver, OF2112. I modified the solver and coupled the density, surface tension and viscosity of solution to the concentration of sugar which is different in various parts of the domain. In order to solve the distribution of sugar (a passive scalar) in the geometry, I added a new equation to the solver as below:

fvScalarMatrix CEqn
(
fvm::ddt(C)
+ fvm::div(phi, C)
- fvm::laplacian(dc,C)
==
fvOptions(C)
);


CEqn.relax();
fvOptions.constrain(CEqn);
CEqn.solve();
fvOptions.correct(C);

Now, the problem is, sugar concentration penetrates inside the bubble which is not correct. I would like to know how can I prevent sugar entering the bubble? I would be more than happy if you share your opinion with me.
sadra2003 is offline   Reply With Quote

Reply

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
How to modify interFoam to keep a constant velocity field thibault_pringuey OpenFOAM 0 January 9, 2009 05:59
Problem with the pressure field using interFoam zoune OpenFOAM Running, Solving & CFD 20 February 4, 2008 18:42
Problems calculating field gh with interFoam cricke OpenFOAM Running, Solving & CFD 0 December 10, 2007 07:17
Putting submesh field values into field on parent mesh helmut OpenFOAM Running, Solving & CFD 2 June 20, 2006 07:31
field by field decomposition vs. component-wise zhanglei Main CFD Forum 0 January 5, 2001 09:06


All times are GMT -4. The time now is 15:45.