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();