42 scalar time0 = time0_;
54 const label nHoles = it->
nHoles();
57 scalar mass = it->
mass(time0, time, twoD_, angleOfWedge_);
64 scalar mp = mass/Np/nHoles;
68 scalar constT =
max(0.0, it->
tsoi() - time0);
81 for(label j=0; j<Np; j++)
84 scalar toi = time0 + constT + deltaT*j/scalar(Np);
86 for(label n=0; n<nHoles; n++)
107 scalar symComponent = direction & axisOfSymmetry_;
108 vector normal = direction - symComponent*axisOfSymmetry_;
109 normal /=
mag(normal);
115 label injectorCell = mesh_.
findCell(injectionPosition);
119 if (injectorCell >= 0)
121 scalar liquidCore = 1.0;
147 injectedLiquidKE_ += 0.5*pPtr->
m()*
magSqr(U);
149 scalar dt = time - toi;
155 bool keepParcel = pPtr->
move(*
this);