static double hyij(long i, long j) { return (Delta/12.0) * (12.0*(alphaC(i)*ad(j)) +4.0*( betaC(i)*ad(j) + alphaC(i)*bd(j) + gammaC(i)*ad(j) + alphaC(i)*cd(j)) +1.0*( betaC(i)*cd(j) + gammaC(i)*bd(j)) +2.0*( betaC(i)*bd(j) + gammaC(i)*cd(j)) ); }
void Foam::relativeVelocityModel::update() { tmp<volVectorField> URel(Ur()); tmp<volScalarField> betaC(alphaC_*rhoC_); tmp<volScalarField> betaD(alphaD_*rhoD_); tmp<volScalarField> rhoM(betaC() + betaD()); tmp<volVectorField> Udm = URel()*betaC()/rhoM; tmp<volVectorField> Ucm = Udm() - URel; Udm_ = Udm(); tau_ = betaD*sqr(Udm) + betaC*sqr(Ucm); }