Commit b333a311 authored by andy's avatar andy
Browse files

ENH: Updated film enthalpy source treatment

parent 6fea2cd9
...@@ -467,7 +467,7 @@ void Foam::ThermoSurfaceFilm<CloudType>::splashInteraction ...@@ -467,7 +467,7 @@ void Foam::ThermoSurfaceFilm<CloudType>::splashInteraction
// Apply correction to velocity for 2-D cases // Apply correction to velocity for 2-D cases
meshTools::constrainDirection(mesh, mesh.solutionD(), pPtr->U()); meshTools::constrainDirection(mesh, mesh.solutionD(), pPtr->U());
Info<< "NEW PARTICLE: " << *pPtr << endl;
// Add the new parcel // Add the new parcel
this->owner().addParticle(pPtr); this->owner().addParticle(pPtr);
......
...@@ -273,15 +273,12 @@ void thermoSingleLayer::solveEnergy() ...@@ -273,15 +273,12 @@ void thermoSingleLayer::solveEnergy()
updateSurfaceTemperatures(); updateSurfaceTemperatures();
dimensionedScalar hs0("SMALL", hs_.dimensions(), SMALL);
solve solve
( (
fvm::ddt(deltaRho_, hs_) fvm::ddt(deltaRho_, hs_)
+ fvm::div(phi_, hs_) + fvm::div(phi_, hs_)
== ==
// - hsSp_ - hsSp_
- fvm::Sp(hsSp_/(hs_ + hs0), hs_)
+ q(hs_) + q(hs_)
+ radiation_->Shs() + radiation_->Shs()
- fvm::SuSp(rhoSp_, hs_) - fvm::SuSp(rhoSp_, hs_)
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment