zhyang-dev/OpenFOAM-6
1
0
Fork 0
mirror of https://github.com/OpenFOAM/OpenFOAM-6.git synced 2025-12-08 06:57:46 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

constTransport: Handle the case of nMoles = 0 in += and -= operators

Browse Source 
Resolves bug-report http://openfoam.org/mantisbt/view.php?id=1348
This commit is contained in:
Henry
2015-02-25 10:54:50 +00:00
parent d5970175a5
commit 99965b01c4
1 changed files with 14 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

22
src/thermophysicalModels/specie/transport/const/constTransportI.H
View File

@ -153,11 +153,14 @@ inline void Foam::constTransport<Thermo>::operator+=
Thermo::operator+=(st);
molr1 /= this->nMoles();
scalar molr2 = st.nMoles()/this->nMoles();
if (mag(molr1) + mag(st.nMoles()) > SMALL)
{
molr1 /= this->nMoles();
scalar molr2 = st.nMoles()/this->nMoles();
mu_ = molr1*mu_ + molr2*st.mu_;
rPr_ = 1.0/(molr1/rPr_ + molr2/st.rPr_);
mu_ = molr1*mu_ + molr2*st.mu_;
rPr_ = 1.0/(molr1/rPr_ + molr2/st.rPr_);
}
}
@ -171,11 +174,14 @@ inline void Foam::constTransport<Thermo>::operator-=
Thermo::operator-=(st);
molr1 /= this->nMoles();
scalar molr2 = st.nMoles()/this->nMoles();
if (mag(molr1) + mag(st.nMoles()) > SMALL)
{
molr1 /= this->nMoles();
scalar molr2 = st.nMoles()/this->nMoles();
mu_ = molr1*mu_ - molr2*st.mu_;
rPr_ = 1.0/(molr1/rPr_ - molr2/st.rPr_);
mu_ = molr1*mu_ - molr2*st.mu_;
rPr_ = 1.0/(molr1/rPr_ - molr2/st.rPr_);
}
}