research/OpenFOAM-2.4.x
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

constTransport: Avoid /0 in + and - operators averaging 1/Pr when the number of moles = 0

Browse source 
Resolves bug-report http://www.openfoam.org/mantisbt/view.php?id=1348
This commit is contained in:
Henry 2015-02-09 21:25:36 +00:00
parent 24f0273b92
commit 3087fcbb1e
1 changed files with 40 additions and 16 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

56
src/thermophysicalModels/specie/transport/const/constTransportI.H
View file

@ -203,15 +203,27 @@ inline Foam::constTransport<Thermo> Foam::operator+
static_cast<const Thermo&>(ct1) + static_cast<const Thermo&>(ct2)
);
scalar molr1 = ct1.nMoles()/t.nMoles();
scalar molr2 = ct2.nMoles()/t.nMoles();
if (mag(ct1.nMoles()) + mag(ct2.nMoles()) < SMALL)
{
return constTransport<Thermo>
(
t,
0,
ct1.rPr_
);
}
else
{
scalar molr1 = ct1.nMoles()/t.nMoles();
scalar molr2 = ct2.nMoles()/t.nMoles();
return constTransport<Thermo>
(
t,
molr1*ct1.mu_ + molr2*ct2.mu_,
1.0/(molr1/ct1.rPr_ + molr2/ct2.rPr_)
);
return constTransport<Thermo>
(
t,
molr1*ct1.mu_ + molr2*ct2.mu_,
1.0/(molr1/ct1.rPr_ + molr2/ct2.rPr_)
);
}
}
@ -227,15 +239,27 @@ inline Foam::constTransport<Thermo> Foam::operator-
static_cast<const Thermo&>(ct1) - static_cast<const Thermo&>(ct2)
);
scalar molr1 = ct1.nMoles()/t.nMoles();
scalar molr2 = ct2.nMoles()/t.nMoles();
if (mag(ct1.nMoles()) + mag(ct2.nMoles()) < SMALL)
{
return constTransport<Thermo>
(
t,
0,
ct1.rPr_
);
}
else
{
scalar molr1 = ct1.nMoles()/t.nMoles();
scalar molr2 = ct2.nMoles()/t.nMoles();
return constTransport<Thermo>
(
t,
molr1*ct1.mu_ - molr2*ct2.mu_,
1.0/(molr1/ct1.rPr_ - molr2/ct2.rPr_)
);
return constTransport<Thermo>
(
t,
molr1*ct1.mu_ - molr2*ct2.mu_,
1.0/(molr1/ct1.rPr_ - molr2/ct2.rPr_)
);
}
}