[kel85uk]

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.

    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