From de9d59ecf15c8f21000ea49693021e9a53cf714b Mon Sep 17 00:00:00 2001 From: Henry Weller Date: Thu, 12 Dec 2019 16:38:36 +0000 Subject: [PATCH] filmViscosityModel::thixotropicViscosity: Improved the handling of the implicit shear term --- .../thixotropicViscosity/thixotropicViscosity.C | 11 +---------- 1 file changed, 1 insertion(+), 10 deletions(-) diff --git a/src/regionModels/surfaceFilmModels/submodels/thermo/filmViscosityModel/thixotropicViscosity/thixotropicViscosity.C b/src/regionModels/surfaceFilmModels/submodels/thermo/filmViscosityModel/thixotropicViscosity/thixotropicViscosity.C index ff2b5171f3..e4329ddb8e 100644 --- a/src/regionModels/surfaceFilmModels/submodels/thermo/filmViscosityModel/thixotropicViscosity/thixotropicViscosity.C +++ b/src/regionModels/surfaceFilmModels/submodels/thermo/filmViscosityModel/thixotropicViscosity/thixotropicViscosity.C @@ -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