37 const scalar molWeight
49 const scalar molWeight
63 molWeight_(st.molWeight_)
99 molWeight_ = st.molWeight_;
105 const scalar sumY = Y_ + st.Y_;
106 if (
mag(sumY) > SMALL)
108 molWeight_ = sumY/(Y_/molWeight_ + st.Y_/st.molWeight_);
125 const scalar sumY =
max(st1.Y_ + st2.Y_, SMALL);
127 if (
mag(sumY) > SMALL)
132 sumY/(st1.Y_/st1.molWeight_ + st2.Y_/st2.molWeight_)
154 scalar diffY = st2.Y_ - st1.Y_;
155 if (
mag(diffY) < SMALL)
160 const scalar diffRW = st2.Y_/st2.molWeight_ - st1.Y_/st1.molWeight_;
165 volatile const bool valid = (
mag(diffRW) > SMALL);
166 const scalar molWeight = valid ? diffY/diffRW : GREAT;
168 scalar molWeight = GREAT;
169 if (
mag(diffRW) > SMALL)
171 molWeight = diffY/diffRW;
175 return specie(diffY, molWeight);
void operator=(const specie &)
Copy assignment, preserve original name.
void operator*=(const scalar)
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
const word & name() const
Name.
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
scalar Y() const
No of moles of this species in mixture.
scalar W() const
Molecular weight [kg/kmol].
Base class of the thermophysical property types.
tmp< faMatrix< Type > > operator+(const faMatrix< Type > &, const faMatrix< Type > &)
word name(const expressions::valueTypeCode typeCode)
A word representation of a valueTypeCode. Empty for INVALID.
tmp< faMatrix< Type > > operator*(const areaScalarField::Internal &, const faMatrix< Type > &)
A class for handling words, derived from Foam::string.
scalar R() const
Gas constant [J/(kg K)].
specie(const specie &)=default
Copy construct.
PtrList< volScalarField > & Y
const scalar RR
Universal gas constant: default in [J/(kmol K)].
tmp< faMatrix< Type > > operator==(const faMatrix< Type > &, const faMatrix< Type > &)
gmvFile<< "tracers "<< particles.size()<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().x()<< " ";}gmvFile<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().y()<< " ";}gmvFile<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().z()<< " ";}gmvFile<< nl;forAll(lagrangianScalarNames, i){ word name=lagrangianScalarNames[i];IOField< scalar > s(IOobject(name, runTime.timeName(), cloud::prefix, mesh, IOobject::MUST_READ, IOobject::NO_WRITE))
void operator+=(const specie &)