[kaifu]

Hi,

I notice there is one line in LISAAtomization.C#L148 in Openfoam 7.0

    vector diff = pos - injectionPos;
    scalar pWalk = mag(diff);
    scalar traveledTime = pWalk/Urel;
    scalar h = diff & injectorDirection_;
    scalar delta = sqrt(sqr(pWalk) - sqr(h));
    scalar hSheet = volFlowRate/(constant::mathematical::pi*delta*Urel);

where,
diff is the displacement vector of current position to the injection position of the parcel,
pWalk is the distance,
h, injectorDirection_ and delta see the following picture.
hSheet is the thickness of liquid sheet. [m]
volFlowRate is volume rate of the injector [m3/s],
Urel is the velocity of parcel [m/s]



I am curious about the hSheet. Should not it be calculated by the following sentence?

hSheet = volFlowRate/(constant::mathematical::twoPi*delta*Urel);