From e89c5736db177d9032459de92407f56799d7936d Mon Sep 17 00:00:00 2001 From: Henry Weller Date: Sat, 26 Feb 2022 09:49:46 +0000 Subject: [PATCH] surfaceFilmModels::waxSolventEvaporation: Improved stability where film thickness -> 0 --- .../waxSolventEvaporation.C | 18 +++++++++++------- 1 file changed, 11 insertions(+), 7 deletions(-) diff --git a/src/regionModels/surfaceFilmModels/submodels/thermo/phaseChangeModel/waxSolventEvaporation/waxSolventEvaporation.C b/src/regionModels/surfaceFilmModels/submodels/thermo/phaseChangeModel/waxSolventEvaporation/waxSolventEvaporation.C index c1d8ccc3cf..b6bb1a2958 100644 --- a/src/regionModels/surfaceFilmModels/submodels/thermo/phaseChangeModel/waxSolventEvaporation/waxSolventEvaporation.C +++ b/src/regionModels/surfaceFilmModels/submodels/thermo/phaseChangeModel/waxSolventEvaporation/waxSolventEvaporation.C @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | Website: https://openfoam.org - \\ / A nd | Copyright (C) 2017-2021 OpenFOAM Foundation + \\ / A nd | Copyright (C) 2017-2022 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -222,7 +222,7 @@ void waxSolventEvaporation::correctModel ); // Local surface temperature at which evaporation takes place - scalarField Tloc(dMass.size()); + scalarField Tloc(T); bool filmPresent = false; @@ -320,14 +320,18 @@ void waxSolventEvaporation::correctModel reduce(filmPresent, orOp()); - const dimensionedScalar rho0Bydt + const volScalarField::Internal rho0Bydt ( - "rho0Bydt", - dimDensity/dimTime, - rootVSmall/dt + rho + *max + ( + dimensionedScalar(dimLength, deltaMin_) - delta, + dimensionedScalar(dimLength, 0) + ) + /(dimensionedScalar(dimTime, dt)*film.VbyA()) ); - volScalarField::Internal impingementRate + const volScalarField::Internal impingementRate ( max (