ENH: Updated film netMass() calculation

This commit is contained in:
andy
2012-11-12 10:47:49 +00:00
parent fdcb0f5d0a
commit d1a0319972

View File

@ -179,7 +179,8 @@ inline tmp<volScalarField> kinematicSingleLayer::netMass() const
{
dimensionedScalar d0("SMALL", dimLength, ROOTVSMALL);
return
fvc::surfaceSum(phi_/(fvc::interpolate(delta_) + d0))*time().deltaT()
fvc::surfaceSum(pos(phi_)*phi_/(fvc::interpolate(delta_) + d0))
*time().deltaT()
+ rho_*delta_*magSf();
}