From b333a3118b9d7bca49a94f8f7ea379eecdfe0320 Mon Sep 17 00:00:00 2001 From: andy Date: Tue, 3 May 2011 15:54:41 +0100 Subject: [PATCH] ENH: Updated film enthalpy source treatment --- .../SurfaceFilmModel/ThermoSurfaceFilm/ThermoSurfaceFilm.C | 2 +- .../surfaceFilmModels/thermoSingleLayer/thermoSingleLayer.C | 5 +---- 2 files changed, 2 insertions(+), 5 deletions(-) diff --git a/src/lagrangian/intermediate/submodels/Thermodynamic/SurfaceFilmModel/ThermoSurfaceFilm/ThermoSurfaceFilm.C b/src/lagrangian/intermediate/submodels/Thermodynamic/SurfaceFilmModel/ThermoSurfaceFilm/ThermoSurfaceFilm.C index 541e983155..cfa6b785ca 100644 --- a/src/lagrangian/intermediate/submodels/Thermodynamic/SurfaceFilmModel/ThermoSurfaceFilm/ThermoSurfaceFilm.C +++ b/src/lagrangian/intermediate/submodels/Thermodynamic/SurfaceFilmModel/ThermoSurfaceFilm/ThermoSurfaceFilm.C @@ -467,7 +467,7 @@ void Foam::ThermoSurfaceFilm::splashInteraction // Apply correction to velocity for 2-D cases meshTools::constrainDirection(mesh, mesh.solutionD(), pPtr->U()); -Info<< "NEW PARTICLE: " << *pPtr << endl; + // Add the new parcel this->owner().addParticle(pPtr); diff --git a/src/regionModels/surfaceFilmModels/thermoSingleLayer/thermoSingleLayer.C b/src/regionModels/surfaceFilmModels/thermoSingleLayer/thermoSingleLayer.C index e6ce90d3af..48c6762b39 100644 --- a/src/regionModels/surfaceFilmModels/thermoSingleLayer/thermoSingleLayer.C +++ b/src/regionModels/surfaceFilmModels/thermoSingleLayer/thermoSingleLayer.C @@ -273,15 +273,12 @@ void thermoSingleLayer::solveEnergy() updateSurfaceTemperatures(); - dimensionedScalar hs0("SMALL", hs_.dimensions(), SMALL); - solve ( fvm::ddt(deltaRho_, hs_) + fvm::div(phi_, hs_) == -// - hsSp_ - - fvm::Sp(hsSp_/(hs_ + hs0), hs_) + - hsSp_ + q(hs_) + radiation_->Shs() - fvm::SuSp(rhoSp_, hs_)