Foam::tmp<Foam::volScalarField> Foam::dragModels::SyamlalOBrien::CdRe() const { volScalarField alpha2 ( max(scalar(1) - pair_.dispersed(), pair_.continuous().residualAlpha()) ); volScalarField A(pow(alpha2, 4.14)); volScalarField B ( neg(alpha2 - 0.85)*(0.8*pow(alpha2, 1.28)) + pos0(alpha2 - 0.85)*(pow(alpha2, 2.65)) ); volScalarField Re(pair_.Re()); volScalarField Vr ( 0.5 *( A - 0.06*Re + sqrt(sqr(0.06*Re) + 0.12*Re*(2.0*B - A) + sqr(A)) ) ); volScalarField CdsRe(sqr(0.63*sqrt(Re) + 4.8*sqrt(Vr))); return CdsRe *max(pair_.continuous(), pair_.continuous().residualAlpha()) /sqr(Vr); }
Foam::tmp<Foam::volScalarField> Foam::dragModels::WenYu::CdRe() const { volScalarField alpha2(max(scalar(1) - pair_.dispersed(), residualAlpha_)); volScalarField Re(pair_.Re()); volScalarField CdsRe ( neg(Re - 1000)*24.0*(1.0 + 0.15*pow(Re, 0.687)) + pos(Re - 1000)*0.44*max(Re, residualRe_) ); return CdsRe *pow(alpha2, -2.65) *max(pair_.continuous(), residualAlpha_); }
Foam::tmp<Foam::volScalarField> Foam::dragModels::GidaspowSchillerNaumann::CdRe() const { volScalarField alpha2 ( max(scalar(1) - pair_.dispersed(), pair_.continuous().residualAlpha()) ); volScalarField Re(alpha2*pair_.Re()); volScalarField CdsRe ( neg(Re - 1000)*24.0*(1.0 + 0.15*pow(Re, 0.687))/alpha2 + pos(Re - 1000)*0.44*max(Re, residualRe_) ); return CdsRe *pow(alpha2, -2.65) *max(pair_.continuous(), pair_.continuous().residualAlpha()); }