[vbchris]

I just wanted to check some code with the community to see if its doing what I think it should.

I want to Implement the following Drag model (based on Roghair, 2011).
http://fcre.tnw.utwente.nl/people/ni...oghair2011.pdf

Cd/[Cinf*(1-alpha)] = ( 1 + 18/Eo*alpha )

Where

Eo is the Eotvos # (1.4 for the system in question)
Cinf is the drag coefficient for a single bubble in pure liquid
Cd is the swarm corrected coefficient
alpha is the gas phase fraction (gas-liquid system)

This is for twoPhaseEulerFoam (OF version 2.2.0)

    volScalarField Re(max(Ur*phase1_.d()/phase2_.nu(), scalar(1.0e-3)));

    volScalarField Cds
    (
        neg(Re - 1000)*(24.0*(1.0 + 0.15*pow(Re, 0.687))/Re)
      + pos(Re - 1000)*0.44
    );

    volScalarField Cde
    (
      	(Cds + Cds*(18/1.4)*alpha1_)*(scalar(1)-alpha1_)
    );

    return 0.75*Cde*phase2_.rho()*Ur/phase1_.d();

[hester]

Hello Chris,

I know the answer is a bit late.

I think your implementation is not correct. I implemented the model after Roghair myself in OF 2.3.0.
Roghair calculates the Cd for a single bubble after Eq. 12- 14 in his paper. You used the Schiller-Naumann model instead. I don't think this will work. Also, its' a bad idea to put a number for the Eötvös number in the equation, it's much better to look it up. This way you can use the drag model for different systems.

This is how I implemented the model. I use values for residualRe and residualEo, which I define in my constant/ directory, to make sure Re and Eo never become zero.

    residualRe_("residualRe", dimless, dict.lookup("residualRe")),
    residualEo_("residualEo", dimless, dict.lookup("residualEo"))

Then I calculate CdRe() using the following code.

Foam::tmp<Foam::volScalarField> Foam::dragModels::Roghair::CdRe() const
{
    volScalarField Eo(max(pair_.Eo(), residualEo_));
    volScalarField Re(max(pair_.Re(), residualRe_));
    volScalarField alpha(pair_.dispersed());

    volScalarField CdRe = 16./Re * (1. + 2./(1. + 16./Re + 3.315/sqrt(Re)));
    volScalarField CdEo = 4.* Eo/(Eo + 9.5);
    volScalarField Cdsingle = sqrt(pow(CdRe, 2.) + pow(CdEo, 2.));
    volScalarField Cd = Cdsingle * (1. - alpha) * (1. + 18./Eo *alpha);

    return Cd * Re;
}

I'm not sure if this would look the same in OF 2.2.0.

Regards,
hester