# Getting a concentration field around a bubble in InterFoam

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

 October 22, 2008, 14:51 Hi all, I have been playin #1 New Member   Azman Join Date: Mar 2009 Location: Aachen, Germany Posts: 3 Rep Power: 8 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.

 March 16, 2012, 08:26 Mass transfer + Interfoam #2 Member   Hanniel Freitas Join Date: Jan 2012 Location: Natal, Brazil Posts: 33 Rep Power: 5 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 .

 January 7, 2013, 14:47 #3 Senior Member     Anton Kidess Join Date: May 2009 Location: Delft, Netherlands Posts: 919 Rep Power: 17 For future reference: InterFoam: Add an equation to only one phase __________________ *On twitter @akidTwit *Spend as much time formulating your questions as you expect people to spend on their answer. *Help define the OpenFOAM stackexchange Q&A site: http://area51.stackexchange.com/prop...oam-technology

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post thibault_pringuey OpenFOAM 0 January 9, 2009 06:59 zoune OpenFOAM Running, Solving & CFD 20 February 4, 2008 19:42 cricke OpenFOAM Running, Solving & CFD 0 December 10, 2007 08:17 helmut OpenFOAM Running, Solving & CFD 2 June 20, 2006 07:31 zhanglei Main CFD Forum 0 January 5, 2001 10:06

All times are GMT -4. The time now is 23:32.