42 sum += injectors_[i].properties()->injectedMass(t);
55 sum += injectors_[i].properties()->mass();
68 label Nf = fuels_->components().size();
72 scalar
T = injectors_[i].properties()->T(time);
75 scalar hl = fuels_->hl(pi, T, X);
76 scalar Wl = fuels_->W(X);
79 for(label j=0; j<Nf; j++)
81 label
k = liquidToGasIndex_[j];
85 sum += injectors_[i].properties()->injectedMass(time)*(hg-hl);
135 for(label j=0; j<Nf; j++)
137 label
k = liquidToGasIndex_[j];
144 scalar
h = hg - hlat;
171 label celli =
elmnt().cell();
173 scalar
pc =
p()[celli];
179 for(label j=0; j<Nf; j++)
181 label
k = liquidToGasIndex_[j];
189 scalar
h = hg - hlat + (pc - psat)/rho;
215 sum +=
elmnt().m()*ke;
232 return injectedLiquidKE_;
249 label nHoles = injectors_[nozzlei].properties()->nHoles();
253 for(label i=0;i<nHoles;i++)
255 ip += injectors_[nozzlei].properties()->position(i);
261 ip = injectors_[nozzlei].properties()->position(0);
284 dist[n] =
mag(first().position() - ip);
295 scalar de =
mag(
elmnt().position() - ip);
296 scalar me =
elmnt().m();
306 while ( ( i < n-1 ) && ( !found ) )
311 for(label j=n; j>i; j--)
334 scalar mLimit = prc*
mTot;
335 scalar mOff = (1.0 - prc)*mTot;
339 if (mLimit > mTot - m[Np-1])
345 scalar mOffSum = 0.0;
348 while ((mOffSum < mOff) && (i>0))
374 scalar numerator = 0.0, denominator = VSMALL;
383 label celli =
elmnt().cell();
384 scalar Pc =
p()[celli];
389 numerator += tmp*
elmnt().d();
396 return numerator/denominator;