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

interFoam and Continuous-Species-Transfer (CST) method

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 4, 2011, 05:22
Default interFoam and Continuous-Species-Transfer (CST) method
  #1
New Member
 
Kelvin Loh
Join Date: Mar 2009
Posts: 25
Rep Power: 8
kel85uk is on a distinguished road
Hi all,

I need some help in determining if I am on the right track to implement the jump conditions for species concentration across the air-water interface for a bubble. The base solver I am using is interFoam without AMR/MRF functions. I thought of keeping it simple first before adding more functionality.

Code:
  forAll(Diffusion_cell,cID)
    {
        Diffusion_cell[cID]=D_spec_liq.value()*D_spec_gas.value()/(D_spec_gas.value()*alpha1[cID] + D_spec_liq.value()*(1.0 - alpha1[cID]));
    }

        turbulence->correct();

    forAll(gamma_diff,cIID)
    {
        gamma_diff[cIID] = -Diffusion_cell[cIID]*(1-He)*C[cIID]/(alpha1[cIID] + He*(1.0 - alpha1[cIID]));
    }
        solve(fvm::ddt(C) + fvc::div(phi,C) - fvc::laplacian(Diffusion_cell,C) == fvc::laplacian(gamma_diff,alpha1));
The Diffusion_cell, and gamma_diff variables are of the volScalarField type. I tried a simple cube test case with no flow velocity imposed, but the run was unsuccessful. Also, I noticed that the value of grad(alpha1) is not 0 inside the phases, which cannot be correct.

The equations are based on one of the references (Haroun et al, 2008) in Holger Marschall's 5th OpenFOAM workshop abstract of the same case.

http://www.tfd.chalmers.se/~hani/OFW5/timeSchedule_workshop_online-filer/Abstract_final/HolgerMarschallAbstractOFW5.pdf

Any hints/help whatsoever is much appreciated.

Thank you.

Regards,
kel85uk
Attached Files
File Type: c interFoam_SDWA.C (6.2 KB, 5 views)
kel85uk is offline   Reply With Quote

Old   May 10, 2011, 03:59
Default
  #2
New Member
 
tiantian
Join Date: Sep 2010
Posts: 1
Rep Power: 0
cqworm is on a distinguished road
I am doing the same work.
If possible you could email me for more convenient exchange ideas.

Looking for more professional ones posting here...
cqworm is offline   Reply With Quote

Old   December 22, 2012, 20:05
Default Jump conditions in interFoam
  #3
Member
 
Hanniel Freitas
Join Date: Jan 2012
Location: Natal, Brazil
Posts: 33
Rep Power: 5
hfsf is on a distinguished road
Guys, i've strumbled upon the same question. Recently I've got some help in a post that could be very useful: Equilibrium thru Interface

In that post the user nimasan described and alternative approach to use the jump conditions within the VOF approach, with some interesting reference background.

Best,

hfsf
hfsf 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



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