示例#1
0
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);
}