Minor changes to the momentum source and particle velocity

This commit is contained in:
Achuth1992
2020-03-04 13:06:56 +01:00
parent 74052b2319
commit 2af9d37307
2 changed files with 2 additions and 2 deletions

View File

@ -5,7 +5,7 @@ fvVectorMatrix UEqn
+ turbulence->divDevReff(U)
==
fvOptions(U)
+ (lambda*(1-voidfractionNext)/U.mesh().time().deltaT())*(Us - fvm::Sp(1,U))
+ (lambda*(1-voidfractionNext)/U.mesh().time().deltaT())*(fvc::Sp(1,Us)-fvm::Sp(1,U))
);
UEqn.relax();

View File

@ -177,7 +177,7 @@ void cfdemCloudIBmodified::calcForcingTerm(volVectorField& Us)
// calc particle velocity
//for(int i=0;i<3;i++) rVec[i]=Us.mesh().C()[cell][i]-position(index)[i];
for(int i=0;i<3;i++) uP[i] = velocities()[index][i];
Us[cell] = uP;//(1-voidfractions_[index][subCell])*uP;
Us[cell] = (1-voidfractions_[index][subCell])*uP;
}
}
}