ENH: Thermo film - added vapour recoil pressure effect

This commit is contained in:
andy 2012-06-11 14:27:23 +01:00
parent af6c3269b3
commit bc4824dcd9

View File

@ -255,6 +255,9 @@ void thermoSingleLayer::updateSubmodels()
// Update source fields
hsSp_ += primaryEnergyPCTrans_/magSf()/time().deltaT();
rhoSp_ += primaryMassPCTrans_/magSf()/time().deltaT();
// Vapour recoil pressure
pSp_ -= sqr(primaryMassPCTrans_/magSf()/time_.deltaT())/2.0/rhoPrimary_;
}