Commit 584fed88 authored by Henry Weller's avatar Henry Weller
Browse files

twoPhaseEulerFoam/interfacialModels/dragModels/WenYu: Corrected Re -> Re*alpha_g

Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1794
parent 4fad84b6
......@@ -68,16 +68,16 @@ Foam::tmp<Foam::volScalarField> Foam::dragModels::WenYu::CdRe() const
max(scalar(1) - pair_.dispersed(), pair_.continuous().residualAlpha())
);
volScalarField Re(pair_.Re());
volScalarField CdsRe
volScalarField Res(alpha2*pair_.Re());
volScalarField CdsRes
(
neg(Re - 1000)*24.0*(1.0 + 0.15*pow(Re, 0.687))
+ pos(Re - 1000)*0.44*max(Re, residualRe_)
neg(Res - 1000)*24.0*(1.0 + 0.15*pow(Res, 0.687))
+ pos(Res - 1000)*0.44*max(Res, residualRe_)
);
return
CdsRe
*pow(alpha2, -2.65)
CdsRes
*pow(alpha2, -3.65)
*max(pair_.continuous(), pair_.continuous().residualAlpha());
}
......
......@@ -68,16 +68,16 @@ Foam::tmp<Foam::volScalarField> Foam::dragModels::WenYu::CdRe() const
max(scalar(1) - pair_.dispersed(), pair_.continuous().residualAlpha())
);
volScalarField Re(pair_.Re());
volScalarField CdsRe
volScalarField Res(alpha2*pair_.Re());
volScalarField CdsRes
(
neg(Re - 1000)*24.0*(1.0 + 0.15*pow(Re, 0.687))
+ pos(Re - 1000)*0.44*max(Re, residualRe_)
neg(Res - 1000)*24.0*(1.0 + 0.15*pow(Res, 0.687))
+ pos(Res - 1000)*0.44*max(Res, residualRe_)
);
return
CdsRe
*pow(alpha2, -2.65)
CdsRes
*pow(alpha2, -3.65)
*max(pair_.continuous(), pair_.continuous().residualAlpha());
}
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment