“…It is well known that steady and unsteady jet/puff systems can be well approximated by analytic models that assume axial symmetry and a self similar profile for the average centerline and radial components of the velocity field in cylindrical coordinates [94, 95, 103, 104] (see figure 4) where f, g are empiric Gaussian or polynomial functions of the self similar variable η = r/z and the centerline velocity is U c = U z for r = 0 along the z axis, hence f ( η ), g ( η ) must satisfy U z = U c and U r = 0 at r = 0 (see examples in [77, 94, 95, 96, 97, 98, 99, 100, 101]). An axially symmetric self similar jet/puff system fulfills the conservation of linear specific momentum Q = V U c (puff) and force where V is the penetration volume [94, 100, 101], hence for an initial time t = t 0 .…”