A theoretical model predicting forward blood spatter patterns resulting from a round nose bullet gunshot wound is proposed. The chaotic disintegration of a blood layer located ahead and aside of the bullet is considered in the framework of percolation theory. The size distribution of blood drops is determined, which allows for the prediction of a blood spatter cloud being ejected from the rear side of the target where the bullet exits. Then, droplet trajectories are numerically predicted accounting for gravity and air drag, which is affected by the collective aerodynamic interaction of drops through air. The model predicts the number and area of individual stains, as well as the stain distribution as a function of distance from the region of origin. The theoretical predictions are compared with experimental data acquired in this work from 9 mm Luger copper full metal jacket bullets fired from a handgun. The agreement between the predicted and experimentally measured parameters is found to be good.

Turbulent jets are known to support large-scale vortical wave packets traveling downstream. We show that a propagating helical wave represents a common form of the “optimal” eigenfunction tracking these structures from the near to the far field of a round jet issuing from a pipe. Two first mirror-symmetric modes containing around $5\%$ of the total turbulent kinetic energy capture all significant large-scale events and accurately replicate the full shear-layer dynamics of the azimuthal wave number $m=1$. A family of the most energy-containing traveling waves represents low wave numbers and is described in terms of “empirical” dispersion laws.

Laboratory experiments of the granular collapse of an initially dense-packing column in a fluid are reported. By extracting the temporal evolution of the granular material height profile, both the dynamics of the granular flow and the final deposit are characterized depending on the Stokes number $\text{St}$, based on a dissipative process at the grain scale, the grain-fluid density ratio $r$, and the aspect ratio $a$ of the initial column. A full description of the granular collapse including the transient dynamics of the flow and the characterization of the final shape deposit is proposed. The main contribution of the present study is to provide (i) the $\text{St}$ dependence on the granular collapse beyond the effect of the aspect ratio $a$ already reported in previous studies, (ii) the characterization of granular flow regimes in the $({(d/H_i)}^{1/2}\text{St},{(d/H_i)}^{1/2}r)$ plane, where $d/H_i$ is the particle diameter to initial column height ratio, and (iii) simple correlations to describe the granular collapse which would be of interest for geophysical purposes.

The aim of this paper is to present a fully compressible three-phase (liquid, vapor, and air) model and its application to the simulation of in-nozzle cavitation effects on liquid atomization. The model employs a combination of the homogeneous equilibrium barotropic cavitation model with an implicit sharp interface capturing volume of fluid (VOF) approximation. The numerical predictions are validated against the experimental results obtained for injection of water into the air from a step nozzle, which is designed to produce asymmetric cavitation along its two sides. Simulations are performed for three injection pressures, corresponding to three different cavitation regimes, referred to as cavitation inception, developing cavitation, and hydraulic flip. Model validation is achieved by qualitative comparison of the cavitation, spray pattern, and spray cone angles. The flow turbulence in this study is resolved using the large-eddy simulation approach. The simulation results indicate that the major parameters that influence the primary atomization are cavitation, liquid turbulence, and, to a smaller extent, the Rayleigh-Taylor and Kelvin-Helmholtz aerodynamic instabilities developing on the liquid-air interface. Moreover, the simulations performed indicate that periodic entrainment of air into the nozzle occurs at intermediate cavitation numbers, corresponding to developing cavitation (as opposed to incipient and fully developed cavitation regimes); this transient effect causes a periodic shedding of the cavitation and air clouds and contributes to improved primary atomization. Finally, the cone angle of the spray is found to increase with increased injection pressure but drops drastically when hydraulic flip occurs, in agreement with the relevant experiments.

We measure and compare the out-of-plane motion in three experimental configurations for approximating two-dimensional flow with electromagnetically driven thin fluid layers. A prior study found that out-of-plane motion grows suddenly when the Reynolds number $\text{Re}$ exceeds a critical value ${\text{Re}}_c$ in a two-layer miscible configuration [Kelley and Ouellette, Phys. Fluids 23, 045103 (2011)]. Here, we confirm that observation; however, we find that a similar onset does not occur in either a single-layer or two-layer immiscible configuration for the ranges $\text{Re}<520$ and $\text{Re}<740$, respectively. Below the critical Reynolds number, the three configurations have out-of-plane motion with similar magnitude. We confirm that two different normalized measures of out-of-plane motion show similar trends among the three configurations as $\text{Re}$ varies. Finally, we compute the vertical velocity profiles using an analytical model of each of the three configurations and provide further evidence suggesting the transition observed in the miscible configuration is due to a shear instability. Our results lead to suggestions for future experimentalists: The single-layer and immiscible configurations most closely approximate two-dimensional flow over a wide range of $\text{Re}$, though the miscible configuration minimizes clumping of tracer particles.

It is well known that when a flexible or flexibly mounted structure is placed perpendicular to the flow of a Newtonian fluid, it can oscillate due to the shedding of separated vortices at high Reynolds numbers. Unlike Newtonian fluids, the flow of viscoelastic fluids can become unstable even at infinitesimal Reynolds numbers due to a purely elastic flow instability that can occur at large Weissenberg numbers. Recent work has shown that these elastic flow instabilities can drive the motion of flexible sheets. The fluctuating fluid forces exerted on the structure from the elastic flow instabilities can lead to a coupling between an oscillatory structural motion and the state of stress in the fluid flow. In this paper, we present the results of an investigation into the flow of a viscoelastic wormlike micelle solution past a flexible circular cylinder. The time variation of the flow field and the state of stress in the fluid are shown using a combination of particle image tracking and flow-induced birefringence images. The static and dynamic responses of the flexible cylinder are presented for a range of flow velocities. The nonlinear dynamics of the structural motion is studied to better understand an observed transition from a symmetric to an asymmetric structural deformation and oscillation behavior.

We consider the rheology of steady two-dimensional granular flows, in different geometries, using discrete element method–based simulations of soft spheres. The flow classification parameter $\left(\psi \right)$, which defines the local flow type (ranging from pure rotation to simple shear to pure extension), varies spatially, to a significant extent, in the flows. We find that the material behaves as a generalized Newtonian fluid. The $\mu \text{−}I$ scaling proposed by Jop et al. [Nature (London) 441, 727 (2006)] is found to be valid in both two-dimensional and unidirectional flows, as observed in previous studies; however, the data for each flow geometry fall on a different curve. The results for the two-dimensional silo flow indicate that the viscosity does not depend directly on the flow type parameter, $\psi$. We find that the scaling based on “granular fluidity” [Zhang and Kamrin, Phys. Rev. Lett. 118, 058001 (2017)] gives good collapse of the data to a single curve for all the geometries. The data for the variation of the solid faction with inertial number show a reasonable collapse for the different geometries.

In this paper, we report an inertial instability that occurs in electro-osmotically driven channel flows. We assume that the charge motion under the influence of an externally applied electric field is confined to a small vicinity of the channel walls that, effectively, drives a bulk flow through a prescribed slip velocity at the boundaries. Here, we study spatially periodic wall velocity modulations in a two-dimensional straight channel numerically. At low slip velocities, the bulk flow consists of a set of vortices along each wall that are left-right symmetric, while at sufficiently high slip velocities, this flow loses its stability through a supercritical bifurcation. Surprisingly, the flow state that bifurcates from a left-right symmetric base flow has a rather strong mean component along the channel, which is similar to pressure-driven velocity profiles. The instability sets in at rather small Reynolds numbers of about 20–30, and we discuss its potential applications in microfluidic devices.

Recently the generation of eddy currents by interacting surface waves was observed experimentally. The phenomenon provides the possibility for manipulation of particles which are immersed in the fluid. The analysis shows that the amplitude of the established eddy currents produced by stationary surface waves does not depend on the fluid viscosity in the free surface case. The currents become parametrically larger, being inversely proportional to the square root of the fluid viscosity in the case when the fluid surface is covered by an almost incompressible thin liquid (i.e., shear elasticity is zero) film formed by an insoluble agent with negligible internal viscous losses as compared to the dissipation in the fluid bulk. Here we extend the theory for a thin insoluble film with zero shear elasticity and small shear and dilational viscosities on the case of an arbitrary elastic compression modulus. We find both contributions into the Lagrangian motion of passive tracers, which are the advection by the Eulerian vertical vorticity and the Stokes drift. Whereas the Stokes drift contribution preserves its value for the free surface case outside a thin viscous sublayer, the Eulerian vertical vorticity strongly depends on the fluid viscosity at high values of the film compression modulus. The Stokes drift acquires a strong dependence on the fluid viscosity inside the viscous sublayer; however, the change is compensated by an opposite change in the Eulerian vertical vorticity. As a result, the vertical dependence of the intensity of eddy currents is given by a sum of two decaying exponents with both decrements being of the order of the wave number. The decrements are numerically different, so the Eulerian contribution becomes dominant at some depth for the surface film with any compression modulus.

At atmospheric pressure, a drop of ethanol impacting on a solid surface produces a splash. Reducing the ambient pressure below its atmospheric value suppresses this splash. The origin of this so-called pressure effect is not well understood, and this study presents an in-depth comparison between various theoretical models that aim to predict splashing and simulations. In this paper, the pressure effect is explored numerically by resolving the Navier-Stokes equations at a 3-nm resolution. In addition to reproducing numerous experimental observations, it is found that different models all provide elements of what is observed in the simulations. The skating droplet model correctly predicts the existence and scaling of a gas film under the droplet, the lamella formation theory is able to correctly predict the scaling of the lamella ejection velocity as a function of the impact velocity for liquids with different viscosity, and lastly, the dewetting theory's hypothesis of a lift force acting on the liquid sheet after ejection is consistent with our results.

Ionic wind refers to the acceleration of partially ionized air between two high-voltage electrodes. We study the momentum transfer from ions to air, resulting from ionic wind created by two asymmetric electrodes and producing a net thrust. This electrohydrodynamic (EHD) thrust, has already been measured in previous studies with digital scales. In this study, we provide more insights into the electrohydrodynamic momentum transfer for a wire-to-cylinder(s) positive dc corona discharge. We provide a simple and general theoretical derivation for EHD thrust, which is proportional to the current/mobility ratio and also to an effective distance integrated on the surface of the electrodes. By considering various electrode configurations, our investigation brings out the physical origin of previously obtained optimal configurations, associated with a better tradeoff between Coulomb forcing, friction occurring at the collector, and wake interactions. By measuring two-dimensional velocity fields using particle image velocimetry (PIV), we are able to evaluate the resulting local net force, including the pressure gradient. It is shown that the contribution of velocity fluctuations in the wake of the collecting electrode(s) must be taken into account to recover the net thrust. We confirm the proportionality between the EHD force and the current/mobility ratio experimentally, and evaluate the ion mobility from PIV measurements. A spectral analysis of the velocity fluctuations indicates a dominant frequency corresponding to a Strouhal number of 0.3 based on the ionic wind velocity and the collector size. Finally, the effective mobility of charge carriers is estimated by a PIV based method inside the drift region.

The physics of nucleation in water is an important issue for many areas, ranging from biomedical to engineering applications. Within the present study, we investigate numerically homogeneous nucleation in a microchannel induced by shock reflection to gain a better understanding of the mechanism of homogeneous nucleation. The liquid expands due to the reflected shock and homogeneous cavitation nuclei are generated. An Eulerian-Lagrangian approach is employed for modeling this process in a microchanel. Two-dimensional axisymmetric Euler equations are solved for obtaining the time evolution of shock, gas bubble, and the ambient fluid. The dynamics of dispersed vapor bubbles is coupled with the surrounding fluid in a Lagrangian framework, describing bubble location and bubble size variation. Our results reproduce nuclei distributions at different stages of homogeneous nucleation and are in good agreement with experimental results. We obtain numerical data for the negative pressure that water can sustain under the process of homogeneous nucleation. An energy transformation description for the homogeneous nucleation inside a microchannel flow is derived and analyzed in detail.

An electrochemical high-speed wall shear raster microscope is presented. It involves chronoamperometric measurements on a microelectrode that is flush-mounted in a submerged test specimen. Wall shear rates are derived from the measured microelectrode signal by numerically solving a convection-diffusion equation with an optimization approach. This way, the unsteady wall shear rates from the collapse of a laser pulse seeded cavitation bubble close to a substrate are measured. By planar scanning, they are resolved in high spatial resolution. The wall shear rates are related to the bubble dynamics via synchronized high-speed imaging of the bubble shape.

The flow field of a five-bladed marine propeller operating at design condition, obtained using large eddy simulation, is used to calculate the resulting far-field sound. The results of three acoustic formulations are compared, and the effects of the underlying assumptions are quantified. The integral form of the Ffowcs-Williams and Hawkings (FW-H) equation is solved on the propeller surface, which is discretized into a collection of $N$ radial strips. Further assumptions are made to reduce FW-H to a Curle acoustic analogy and a point-force dipole model. Results show that although the individual blades are strongly tonal in the rotor plane, the propeller is acoustically compact at low frequency and the tonal sound interferes destructively in the far field. The propeller is found to be acoustically compact for frequencies up to 100 times the rotation rate. The overall far-field acoustic signature is broadband. Locations of maximum sound of the propeller occur along the axis of rotation both up and downstream. The propeller hub is found to be a source of significant sound to observers in the rotor plane, due to flow separation and interaction with the blade-root wakes. The majority of the propeller sound is generated by localized unsteadiness at the blade tip, which is caused by shedding of the tip vortex. Tonal blade sound is found to be caused by the periodic motion of the loaded blades. Turbulence created in the blade boundary layer is convected past the blade trailing edge leading to generation of broadband noise along the blade. Acoustic energy is distributed among higher frequencies as local Reynolds number increases radially along the blades. Sound source correlation and spectra are examined in the context of noise modeling.

The nonlinear dynamics of a swirling wake flow stemming from a Graboswksi-Berger vortex [Grabowski and Berger, J. Fluid Mech. 75, 525 (1976)] in a semi-infinite domain is addressed at low Reynolds numbers for a fixed swirl number $S=1.095$, defined as the ratio between the characteristic tangential velocity and the centerline axial velocity. In this system, only pure hydrodynamic instabilities develop and interact through the quadratic nonlinearities of the Navier-Stokes equations. Such interactions lead to the onset of chaos at a Reynolds value of $\text{Re}=220$. This chaotic state is reached by following a Ruelle-Takens-Newhouse scenario, which is initiated by a Hopf bifurcation (the spiral vortex breakdown) as the Reynolds number increases. At larger Reynolds value, a frequency synchronization regime appears followed by a chaotic state again. This scenario is corroborated by nonlinear time series analyses. Stability analysis around the time-average flow and temporal-azimuthal Fourier decomposition of the nonlinear flow distributions both identify successfully the developing vortices and provide deeper insight into the development of the flow patterns leading to this route to chaos. Three single-helical vortices are involved: the primary spiral associated with the spiral vortex breakdown, a downstream spiral, and a near-wake spiral. As the Reynolds number increases, the frequencies of these vortices become closer, increasing their interactions by nonlinearity to eventually generate a strong chaotic axisymmetric oscillation.

The effect of particle edge geometry on the descent motion of freely falling planar particles is examined through experiments. Various planar particles, such as disk and polygons, with identical frontal areas ($A_p$) and different number of edges (or perimeter) are used. All particles are designed such that their values of Galileo number ($G$) and dimensionless moment of inertia ($I^{*}$) correspond to the previously identified fluttering regime of particle motion. Several modes of secondary motion are observed for the same particle and conditions, and these are not equally probable. This probability depends on the particle shape. Disks and heptagons were found to prefer a “planar zig-zag” behavior. These planar motions are composed of gliding sweeps and turning sections. As the number of sides in the polygon decreases, i.e., for hexagons and pentagons, the trajectory transitions to a more three-dimensional form. These trajectories were found to be restricted to one plane per swing, but the subsequent swings are in other planes. Further decrease in number of sides to a square results in the trajectories having a severe out-of-plane motion. These subregimes of particle motion within the fluttering regime are consistent with those reported for disks in previous studies. Based on this information, a length scale that accounts for the frontal area of the particles and its edge geometry (i.e., perimeter) is proposed. This length scale represents an approach to determine an equivalent disk for planar particles such that the phase diagram in the Reynolds number ($\text{Re}$) and dimensionless moment of inertia ($I^{*}$) domain can be used to characterize the motion of planar particles with different frontal geometries. However, further experiments covering other domains of the regime map are needed to verify its universality.

A series of laboratory experiments were performed to study the Peregrine Breather (PB) evolution in a wave flume of finite depth and deep water. Experimental cases were selected with water depths $k_{0}h$ ($k_0$ is the wave number and $h$ is the water depth) varying from 3.11 to 8.17 and initial steepness $k_0{a}_0$ ($a_0$ is the background wave amplitude) in the range 0.06 to 0.12, and the corresponding initial Ursell number in the range 0.03 to 0.061. Experimental results indicate that the water depth plays an important role in the formation of the extreme waves in finite depth; the maximum wave amplification of the PB packets is also strongly dependent on the initial Ursell number. For experimental cases with the initial Ursell number larger than 0.05, the maximum crest amplification can exceed three. If the initial Ursell number is nearly 0.05, a shorter propagation distance is needed for maximum amplification of the height in deeper water. A time-frequency analysis using the wavelet transform reveals that the energy of the higher harmonics is almost in-phase with the carrier wave. The contribution of the higher harmonics to the extreme wave is significant for the cases with initial Ursell number larger than 0.05 in water depth $k_{0}h<$ 5.0. Additionally, the experimental results are compared with computations based on both the nonlinear Schrödinger (NLS) equation and the Dysthe equation, both with a dissipation term. It is found that both models with a dissipation term can predict the maximum amplitude amplification of the primary waves. However, the Dysthe equation also can predict the group horizontal asymmetry.

We study the growth of concentration fluctuations of weakly inertial particles in the turbulent channel flow starting with a smooth initial distribution. The steady-state concentration is singular and multifractal so the growth describes the increasingly rugged structure of the distribution. We demonstrate that inhomogeneity influences the growth of concentration fluctuations profoundly. For homogeneous turbulence the growth is exponential and is fully determined by Kolmogorov scale eddies.We derive lognormality of the statistics in this case. The growth exponents of the moments are proportional to the sum of Lyapunov exponents, which is quadratic in the small inertia of the particles. In contrast, for inhomogeneous turbulence the growth is linear in inertia. It involves correlations of inertial range and viscous scale eddies that turn the growth into a stretched exponential law with exponent three halves. We demonstrate using direct numerical simulations that the resulting growth rate can differ by orders of magnitude over channel height. This strong variation might have relevance in the planetary boundary layer.

Locomotion in Stokes flow is an intensively studied problem because it describes important biological phenomena such as the motility of many species' sperm, bacteria, algae, and protozoa. Numerical computations can be challenging, particularly in three dimensions, due to the presence of moving boundaries and complex geometries; methods which combine ease of implementation and computational efficiency are therefore needed. A recently proposed method to discretize the regularized Stokeslet boundary integral equation without the need for a connected mesh is applied to the inertialess locomotion problem in Stokes flow. The mathematical formulation and key aspects of the computational implementation in matlab or GNU Octave are described, followed by numerical experiments with biflagellate algae and multiple uniflagellate sperm swimming between no-slip surfaces, for which both swimming trajectories and flow fields are calculated. These computational experiments required minutes of time on modest hardware; an extensible implementation is provided in a GitHub repository. The nearest-neighbor discretization dramatically improves convergence and robustness, a key challenge in extending the regularized Stokeslet method to complicated three-dimensional biological fluid problems.