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

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 3, 2011, 04:16
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.

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.1 KB, 5 views)
kel85uk 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 23:04.