59 coeffsDict_(dict.
subDict(typeName +
"Coeffs")),
60 elasticity_(
readScalar(coeffsDict_.lookup(
"elasticity")))
76 const label globalFacei
80 label facei = p.
patchFace(patchi, globalFacei);
84 if (isA<wallPolyPatch>(mesh_.boundaryMesh()[
patchi]))
93 scalar Un = p.
U() & Sf;
97 p.
U() -= (1.0 + elasticity_)*Un*Sf;
105 vector Ub0 = U_.oldTime().boundaryField()[
patchi][facei];
107 scalar dt = spray_.runTime().deltaT().value();
112 vector Cf0 = mesh.
faces()[globalFacei].centre(oldPoints);
114 vector Sf0 = mesh.
faces()[globalFacei].normal(oldPoints);
117 if (
mag(Sf0) > SMALL)
126 scalar magSfDiff =
mag(Sf - Sf0);
130 if (magSfDiff > SMALL)
136 scalar magOmega =
mag(omega);
137 omega /= magOmega+SMALL;
139 scalar phiVel =
::asin(magOmega)/dt;
143 vector vrot = phiVel*(omega ^ (pos - Cf));
147 scalar Un = ((p.
U() - v) & Sfp);
151 p.
U() -= (1.0 + elasticity_)*Un*Sfp;
158 scalar Urn = Ur & Sf;
173 p.
U() -= (1.0 + elasticity_)*Urn*Sf;
182 <<
"bool reflectParcel::wallTreatment(parcel& parcel) const "
183 <<
" parcel has hit a boundary "
184 << mesh_.boundary()[
patchi].type()
185 <<
" which not yet has been implemented."