zhyang-dev/OpenFOAM-12
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

filmViscosityModel::thixotropicViscosity: Improved the handling of the implicit shear term

Browse Source
This commit is contained in:
Henry Weller
2019-12-12 16:38:36 +00:00
parent 75f11686ba
commit de9d59ecf1
1 changed files with 1 additions and 10 deletions
Download Patch File Download Diff File Expand all files Collapse all files

11
src/regionModels/surfaceFilmModels/submodels/thermo/filmViscosityModel/thixotropicViscosity/thixotropicViscosity.C
View File

@ -114,7 +114,6 @@ void thixotropicViscosity::correct
const volScalarField alphaRho = film.alpha()*film.rho();
const surfaceScalarField& phiU = film.phiU();
const volScalarField& coverage = film.coverage();
const Time& runTime = this->film().regionMesh().time();
// Shear rate
const volScalarField gDot
@ -123,11 +122,6 @@ void thixotropicViscosity::correct
coverage*mag(U - Uw)/(delta + film.deltaSmall())
);
if (debug && runTime.writeTime())
{
gDot.write();
}
const dimensionedScalar alphaRho0
(
"alphaRho0",
@ -135,9 +129,6 @@ void thixotropicViscosity::correct
rootVSmall
);
const dimensionedScalar c0("c0", dimless/dimTime, rootVSmall);
const volScalarField coeff("coeff", -c_*pow(gDot, d_) + c0);
fvScalarMatrix lambdaEqn
(
fvm::ddt(lambda_)
@ -145,7 +136,7 @@ void thixotropicViscosity::correct
- fvm::Sp(fvc::div(phiU), lambda_)
==
a_*pow((1 - lambda_), b_)
+ fvm::SuSp(coeff, lambda_)
- fvm::Sp(c_*pow(gDot, d_), lambda_)
// Include the effect of the impinging droplets added with lambda = 0
- fvm::Sp