88 scalar molr1 = this->nMoles();
92 molr1 /= this->nMoles();
93 scalar molr2 = ico.
nMoles()/this->nMoles();
95 rho_ = molr1*rho_ + molr2*ico.rho_;
100 scalar molr1 = this->nMoles();
104 molr1 /= this->nMoles();
105 scalar molr2 = ico.
nMoles()/this->nMoles();
107 rho_ = molr1*rho_ - molr2*ico.rho_;
124 scalar nMoles = ico1.
nMoles() + ico2.nMoles();
125 scalar molr1 = ico1.nMoles()/nMoles;
126 scalar molr2 = ico2.nMoles()/nMoles;
130 static_cast<const specie&>(ico1)
131 + static_cast<const specie&>(ico2),
132 molr1*ico1.rho_ + molr2*ico2.rho_
138 const incompressible& ico1,
139 const incompressible& ico2
142 scalar nMoles = ico1.
nMoles() + ico2.nMoles();
143 scalar molr1 = ico1.nMoles()/nMoles;
144 scalar molr2 = ico2.nMoles()/nMoles;
146 return incompressible
148 static_cast<const specie&>(ico1)
149 - static_cast<const specie&>(ico2),
150 molr1*ico1.rho_ - molr2*ico2.rho_
157 const incompressible& ico
160 return incompressible(s*static_cast<const specie&>(ico), ico.rho_);
165 const incompressible& ico1,
166 const incompressible& ico2