ENH: Changing Yt to "calculated" to avoid division by zero

This commit is contained in:
Sergio Ferraris
2011-10-05 13:37:07 +01:00
parent 82f569972f
commit 8ede7fdda4
5 changed files with 18 additions and 37 deletions

View File

@ -49,7 +49,8 @@ const ThermoType& Foam::multiComponentMixture<ThermoType>::constructSpeciesData
template<class ThermoType>
void Foam::multiComponentMixture<ThermoType>::correctMassFractions()
{
volScalarField Yt("Yt", Y_[0]);
// It changes Yt patches to "calculated"
volScalarField Yt("Yt", 1.0*Y_[0]);
for (label n=1; n<Y_.size(); n++)
{