CFD Online Discussion Forums - interFoam and Continuous-Species-Transfer (CST) method

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