From e9df26cbb8599a74b360a0af310dd5d4dfd1a4a5 Mon Sep 17 00:00:00 2001 From: Henry Weller Date: Sun, 22 Jan 2023 11:47:39 +0000 Subject: [PATCH] surfaceFilmModels::momentumSurfaceFilm: Minor code improvements --- .../momentumSurfaceFilm/momentumSurfaceFilm.C | 12 +++--------- 1 file changed, 3 insertions(+), 9 deletions(-) diff --git a/src/surfaceFilmModels/momentumSurfaceFilm/momentumSurfaceFilm.C b/src/surfaceFilmModels/momentumSurfaceFilm/momentumSurfaceFilm.C index 774daab590..21e1472478 100644 --- a/src/surfaceFilmModels/momentumSurfaceFilm/momentumSurfaceFilm.C +++ b/src/surfaceFilmModels/momentumSurfaceFilm/momentumSurfaceFilm.C @@ -394,15 +394,9 @@ void Foam::momentumSurfaceFilm::solveAlpha ) ); - surfaceScalarField phid + const surfaceScalarField phid ( "phid", - // constrainFilmField - // ( - // rhof - // *constrainPhiHbyA(fvc::flux(HbyA) - alpharAUf*phiu, U_, alpha_), - // 0 - // ) rhof*constrainPhiHbyA(fvc::flux(HbyA) - alpharAUf*phiu, U_, alpha_) ); @@ -412,8 +406,6 @@ void Foam::momentumSurfaceFilm::solveAlpha alphaf*rhof*alpharAUf*rhogf ); - mesh().schemes().setFluxRequired(alpha_.name()); - while (pimple_.correctNonOrthogonal()) { // Film thickness equation @@ -833,6 +825,8 @@ Foam::momentumSurfaceFilm::momentumSurfaceFilm { alpha_ == delta_/VbyA(); + mesh().schemes().setFluxRequired(alpha_.name()); + if (readFields) { transferPrimaryRegionThermoFields();