Many material science, coating, and manufacturing problems involve liquid films where defects that span the film thickness must be removed. Here, we study the surface-tension-driven leveling dynamics of a thin viscous film following closure of an initial hole. The dynamics of the film shape is described by a nonlinear evolution equation, for which we obtain a self-similar solution. The analytical results are verified using time-dependent numerical and experimental results for the profile shapes and the minimum film thickness at the center. The universal behavior we identify can be useful for characterizing the time evolution of the leveling process and estimating material properties from experiments.

In a recent study [Khabaz et al., Phys. Rev. Fluids 2, 093301 (2017)], we showed that jammed soft particle glasses (SPGs) crystallize and order in steady shear flow. Here we investigate the rheology and microstructures of these suspensions in oscillatory shear flow using particle-dynamics simulations. The microstructures in both types of flows are similar, but their evolutions are very different. In both cases the monodisperse and polydisperse suspensions form crystalline and layered structures, respectively, at high shear rates. The crystals obtained in the oscillatory shear flow show fewer defects compared to those in the steady shear. SPGs remain glassy for maximum oscillatory strains less than about the yield strain of the material. For maximum strains greater than the yield strain, microstructural and rheological transitions occur for SPGs. Polydisperse SPGs rearrange into a layered structure parallel to the flow-vorticity plane for sufficiently high maximum shear rates and maximum strains about 10 times greater than the yield strain. Monodisperse suspensions form a face-centered cubic (FCC) structure when the maximum shear rate is low and hexagonal close-packed (HCP) structure when the maximum shear rate is high. In steady shear, the transition from a glassy state to a layered one for polydisperse suspensions included a significant induction strain before the transformation. In oscillatory shear, the transformation begins to occur immediately and with different microstructural changes. A state diagram for suspensions in large amplitude oscillatory shear flow is found to be in close but not exact agreement with the state diagram for steady shear flow. For more modest amplitudes of around one to five times the yield strain, there is a transition from a glassy structure to FCC and HCP crystals, at low and high frequencies, respectively, for monodisperse suspensions. At moderate frequencies, the transition is from glassy to HCP via an intermediate FCC phase.

The so-called wake-moment coefficient ${\tilde{C}}_h$ and lateral wake deflection of three-dimensional windbreaks are explored in the near and far wake. Wind-tunnel experiments were performed to study the functional dependence of ${\tilde{C}}_h$ with windbreak aspect ratio, incidence angle, and the ratio of the windbreak height and surface roughness $\left(h/z_0\right)$. Supported with the data, we also propose basic models for the wake deflection of the windbreak in the near and far fields. The near-wake model is based on momentum conservation considering the drag on the windbreak, whereas the far-wake counterpart is based on existing models for wakes behind surface-mounted obstacles. Results show that ${\tilde{C}}_h$ does not change with windbreak aspect ratios of 10 or greater; however, it may be lower for an aspect ratio of 5. ${\tilde{C}}_h$ is found to change roughly with the cosine of the incidence angle, and to depend strongly on $h/z_0$. The data broadly support the proposed wake-deflection models, though better predictions could be made with improved knowledge of the windbreak drag coefficient.

In this paper we investigate the transition from the vortex to the engulfment regime in a T-mixer when the two entering flows have different viscosity. In particular we consider as working fluid water entering the two inlet channels of the mixer at two different temperatures. Contrary to the isothermal case, at low Reynolds numbers the vortex regime shows only a single reflectional symmetry, due to the nonhomogeneous distribution of the viscosity. Increasing the Reynolds number, a symmetry-breaking bifurcation drives the system to a new steady flow configuration, usually called the engulfment regime, similar to what it is possible to observe in an isothermal case. This flow regime is associated with an increase of the mixing between the two inlet streams. It is shown by direct numerical simulation (DNS) and by stability analysis that the engulfment regime is promoted by the temperature difference. Starting from the DNSs, the resulting flow fields are analyzed in detail considering different temperature jumps between the two inlet boundaries. Furthermore, dedicated linear stability analyses are carried out to investigate the instability mechanism associated with the occurrence of the engulfment regime. In particular, similarly to the case without temperature differences, the onset of engulfment is driven by the momentum equation, and the temperature field does not lead to any additional instability mechanism. However, the existence of a temperature field leads to quantitative changes of the stability characteristics and of the resulting flow fields via a variation of the viscosity coefficient.

Here we consider the problem of shock propagation through a layer of spherical particles. A point particle force model is used to capture the shock-induced aerodynamic force acting upon the particles. The discrete element method (DEM) code liggghts is used to implement the shock-induced force as well as to capture the collisional forces within the system. A volume-fraction-dependent drag correction is applied using Voronoi tessellation to calculate the volume of fluid around each individual particle. A statistically stationary frame is chosen so that spatial and temporal averaging can be performed to calculate ensemble-averaged macroscopic quantities, such as the granular temperature. A parametric study is carried out by varying the coefficient of restitution for three sets of multiphase shock conditions. A self-similar profile is obtained for the granular temperature that is dependent on the coefficient of restitution. A traveling wave structure is observed in the particle concentration downstream of the shock and this instability arises from the volume-fraction-dependent drag force. The intensity of the traveling wave increases significantly as inelastic collisions are introduced. Downstream of the shock, the variance in Voronoi volume fraction is shown to have a strong dependence upon the coefficient of restitution, indicating clustering of particles induced by collisional dissipation. Statistics of the Voronoi volume are computed upstream and downstream of the shock and compared to theoretical results for randomly distributed hard spheres.

Direct numerical simulation (DNS) combined with the Lagrangian point particle model is used to study Rayleigh-Bénard convection in order to understand modifications due to the interaction of inertial, nonisothermal particles with buoyancy-driven turbulence. In this system, turbulence can be altered through direct momentum coupling, as well as through buoyancy modification via thermal coupling between phases. We quantify the effect of the dispersed phase by changes to the total integrated turbulent kinetic energy (TKE) and Nusselt number (Nu). The dispersed particles experience gravitational settling and are introduced at the lower wall so that turbulence must overcome the settling velocity for the particles to vertically distribute throughout the domain. We focus primarily on particle inertia, settling velocity, mass fraction, and the ratio of the particle to fluid specific heat. Furthermore, individual contributions by the momentum coupling and thermal coupling are studied to see which most significantly changes Nu and TKE. Our results show that particles with Stokes number of order unity maximize Nu, corresponding to a peak of clustering and attenuation of TKE. Increased mass fractions lead to a linear increase of Nu and decrease of TKE. With varying specific heat ratio, Nu and TKE exhibit monotonic behaviors, where in the high limit particles become isothermal and depend upon the initialized particle temperature. It is also shown that particles two-way coupled only through momentum attenuate Nu and weaken TKE, while thermal-only coupling also weakens TKE but enhances Nu. When both couplings are present, however, thermal coupling overwhelms the momentum coupling attenuation, and the net result is an enhancement of Nu.

An analytical method is presented for determining the Reynolds shear stress profile in steady, two-dimensional wall-bounded flows using the mean streamwise velocity. The method is then utilized with experimental data to determine the local wall shear stress. The procedure is applicable to flows on smooth and rough surfaces with arbitrary pressure gradients. It is based on the streamwise component of the boundary layer momentum equation, which is transformed into inner coordinates. The method requires velocity profiles from at least two streamwise locations, but the formulation of the momentum equation reduces the dependence on streamwise gradients. The method is verified through application to laminar flow solutions and turbulent DNS results from both zero and nonzero pressure gradient boundary layers. With strong favorable pressure gradients, the method is shown to be accurate for finding the wall shear stress in cases where the Clauser fit technique loses accuracy. The method is then applied to experimental data from the literature from zero pressure gradient studies on smooth and rough walls, and favorable and adverse pressure gradient cases on smooth walls. Data from very near the wall are not required for determination of the wall shear stress. Wall friction velocities obtained using the present method agree with those determined in the original studies, typically to within 2%.

Experiments were performed in a vertical channel to study the properties of a bubbly flow composed of two distinct bubble size species. Bubbles were produced using a capillary bank with tubes with two distinct inner diameters; the flow through each capillary size was controlled such that the amount of large or small bubbles could be controlled. Using water and water-glycerin mixtures, a wide range of Reynolds and Weber number ranges were investigated. The gas volume fraction ranged between 0.5% and 6%. The measurements of the mean bubble velocity of each species and the liquid velocity variance were obtained and contrasted with the monodisperse flows with equivalent gas volume fractions. We found that the bidispersity can induce a reduction of the mean bubble velocity of the large species; for the small size species, the bubble velocity can be increased, decreased, or remain unaffected depending of the flow conditions. The liquid velocity variance of the bidisperse flows is, in general, bound by the values of the small and large monodisperse values; interestingly, in some cases, the liquid velocity fluctuations can be larger than either monodisperse case. A simple model for the liquid agitation for bidisperse flows is proposed, with good agreement with the experimental measurements.

It is known that rapidly rotating turbulent flows are characterized by the emergence of simultaneous upscale and downscale energy transfer. Indeed, both numerics and experiments show the formation of large-scale anisotropic vortices together with the development of small-scale dissipative structures. However the organization of interactions leading to this complex dynamics remains unclear. Two different mechanisms are known to be able to transfer energy upscale in a turbulent flow. The first is characterized by two-dimensional interactions among triads lying on the two-dimensional, three-component (2D3C)/slow manifold, namely on the Fourier plane perpendicular to the rotation axis. The second mechanism is three-dimensional and consists of interactions between triads with the same sign of helicity (homochiral). Here, we present a detailed numerical study of rotating flows using a suite of high-Reynolds-number direct numerical simulations (DNS) within different parameter regimes to analyze both upscale and downscale cascade ranges. We find that the upscale cascade at wave numbers close to the forcing scale is generated by increasingly dominant homochiral interactions which couple the three-dimensional bulk and the 2D3C plane. This coupling produces an accumulation of energy in the 2D3C plane, which then transfers energy to smaller wave numbers thanks to the two-dimensional mechanism. In the forward cascade range, we find that the energy transfer is dominated by heterochiral triads and is dominated primarily by interaction within the fast manifold where $k_z\ne 0$. We further analyze the energy transfer in different regions in the real-space domain. In particular, we distinguish high-strain from high-vorticity regions and we uncover that while the mean transfer is produced inside regions of strain, the rare but extreme events of energy transfer occur primarily inside the large-scale column vortices.

Contact between particles and motile cells underpins a wide variety of biological processes, from nutrient capture and ligand binding to grazing, viral infection, and cell-cell communication. The window of opportunity for these interactions depends on the basic mechanism determining contact time, which is currently unknown. By combining experiments on three different species—Chlamydomonas reinhardtii, Tetraselmis subcordiforms, and Oxyrrhis marina—with simulations and analytical modeling, we show that the fundamental physical process regulating proximity to a swimming microorganism is hydrodynamic particle entrainment. The resulting distribution of contact times is derived within the framework of Taylor dispersion as a competition between advection by the cell surface and microparticle diffusion, and predicts the existence of an optimal tracer size that is also observed experimentally. Spatial organization of flagella, swimming speed, and swimmer and tracer size influence entrainment features and provide tradeoffs that may be tuned to optimize the estimated probabilities for microbial interactions like predation and infection.

The Aplysia, commonly referred to as the “sea hare,” is a marine mollusc that swims using large-amplitude flapping of its wide, winglike parapodia. In this study, flow simulations with a relatively simple kinematical model are used to gain insights into the vortex dynamics, thrust generation, and energetics of locomotion for this animal. A unique vortex pattern characterized by three distinct trains of vortex ringlike structures is observed in the wake of this animal. These vortex rings are associated with a positive momentum flux in the wake that counteracts the drag generated by the body. Simulations indicate propulsive efficiencies of up to 24% and terminal swimming speeds of about 0.9 body length per cycle. Swimming speeds are found to increase with increasing parapodial flapping amplitude as well as wavelength of undulation.

We outline the derivation of a two-phase continuum theory for grains, jumping above a bed of sand, while accelerated by a turbulent shearing flow, colliding with the bed, rebounding, and, perhaps, generating other grains. Relations between the shear and normal stresses and vertical derivatives of components of the average particle velocity are determined by averaging the dynamical equations for the particle trajectories. This provides the closure for the system of differential equations that govern the behavior of the wind and particles above the bed. Boundary conditions are obtained by averaging the results of experiments on rebound and ejection of particles from a particle bed. We solve the resulting system of equations subject to the derived boundary conditions for steady, uniform flows over both particle and rigid beds, and obtain unsteady, uniform solutions and steady, nonuniform solutions that provide information regarding saturation times and lengths, respectively.

In this paper, we consider how the statistical moments of the separation between two fluid particles grow with time when their separation lies in the dissipation range of turbulence. In this range, the fluid velocity field varies smoothly and the relative velocity of two fluid particles depends linearly upon their separation. While this may suggest that the rate at which fluid particles separate is exponential in time, this is not guaranteed because the strain rate governing their separation is a strongly fluctuating quantity in turbulence. Indeed, Afik and Steinberg [Nat. Commun. 8, 468 (2017)] argue that there is no convincing evidence that the moments of the separation between fluid particles grow exponentially with time in the dissipation range of turbulence. Motivated by this, we use direct numerical simulations (DNS) to compute the moments of particle separation over very long periods of time in a statistically stationary, isotropic turbulent flow to see if we ever observe evidence for exponential separation. Our results show that if the initial separation between the particles is infinitesimal, the moments of the particle separation first grow as power laws in time, but we then observe convincing evidence that at sufficiently long times the moments do grow exponentially. However, this exponential growth is only observed after extremely long times $\gtrsim 200{\tau }_{\eta }$, where ${\tau }_{\eta }$ is the Kolmogorov time scale. This is due to fluctuations in the strain rate about its mean value measured along the particle trajectories, the effect of which on the moments of the particle separation persists for very long times. We also consider the backward-in-time (BIT) moments of the article separation, and observe that they too grow exponentially in the long-time regime. However, a dramatic consequence of the exponential separation is that at long times the difference between the rate of the particle separation forward in time (FIT) and BIT grows exponentially in time, leading to incredibly strong irreversibility in the dispersion. This is in striking contrast to the irreversibility of their relative dispersion in the inertial range, where the difference between FIT and BIT is constant in time according to Richardson's phenomenology.

Many common engineering problems involve the study of turbulence interaction with other physical processes. For many such physical processes, solutions are expressed most naturally in physical space, necessitating the use of physical space solutions. For simulating isotropic turbulence in physical space, linear forcing is a commonly used strategy because it produces realistic turbulence in an easy-to-implement formulation. However, the method resolves a smaller range of scales on the same mesh than spectral forcing. We propose an alternative approach for turbulence forcing in physical space that uses the low-pass filtered velocity field as the basis of the forcing term. This method is shown to double the range of scales captured by linear forcing while maintaining the flexibility and low computational cost of the original method. This translates to a 60% increase of the Taylor microscale Reynolds number on the same mesh. An extension is made to scalar mixing wherein a scalar field is forced to have an arbitrarily chosen, constant variance. Filtered linear forcing of the scalar field allows for control over the length scale of scalar injection, which could be important when simulating scalar mixing.

Time-resolved stereoscopic particle image velocimetry measurements were conducted of a wing-tip vortex decaying in free-stream turbulence. The objective of the research was to experimentally investigate the mechanism causing the increased rate of decay of the vortex in the presence of turbulence. It was observed that the circulation of the vortex core experienced periods of rapid loss and recovery when immersed in free-stream turbulence. These events were not observed when the vortex was in a laminar free stream. A connection was made between these events and distortion of the vortex, coinciding with stripping of core fluid from the vortex core. Specifically, vortex stripping events were connected to asymmetry in the vortex core, and this asymmetry was associated with instances of rapid circulation loss. The increased rate of decay of the vortex in turbulence coincided with the formation of secondary vortical structures which wrapped azimuthally around the primary vortex.

Lock-release experiments are performed focusing upon the evolution of near-pure glycerol flowing into fresh water. If the lock height is sufficiently tall, the current is found to propagate for many lock lengths close to the speed predicted for energy-conserving moderately non-Boussinesq gravity currents. The current then slows to a near stop as the current head ceases to be elevated relative to its tail and the current as a whole forms a wedge shape. By contrast, an experiment of near-pure glycerol advancing under air exhibits the well-known slowing of the current such that the front position increases as a one-fifth power of time. The evolution of a viscous gravity current in water is also qualitatively different from that for a high-Reynolds number gravity current which transitions smoothly from a constant speed to self-similar to viscous regime. The reason a viscous gravity current flowing under water moves initially at near-constant speed is not due to a lubrication layer forming below the current. Rather it is due to the return flow of water into the lock establishing a current with an elevated head that is taller than the viscous boundary layer depth near the current nose. The flow near the top of the head advances to the nose where it comes into contact with the tank bottom. Meanwhile the ambient fluid is pushed up and over the head rather than being drawn underneath it. The front slows rapidly to a near stop as the head height reduces to that comparable to the boundary layer depth underneath the head. The initial speed and entrainment into the current are shown to depend upon the ratio, $R_{\ell }$, of the starting current height to the characteristic boundary layer depth. In particular, entrainment via the turbulent shear flow over the head is found to increase the volume by less than $10\%$ during its evolution if $R_{\ell }\lesssim 10$ but increases by as much as $100\%$ for high-Reynolds number gravity currents. A conceptual model is developed that captures the transition from an inertially driven current to its sudden near stop by viscous forces.

The evaporation of droplets within turbulent sprays involves unsteady, multiscale, and multiphase processes which make its comprehension and modeling capabilities still limited. The present work aims to investigate the dynamics of droplet vaporization within a turbulent spatial developing jet in dilute, nonreacting conditions. We address the problem considering a turbulent jet laden with acetone droplets and using the direct numerical simulation framework based on a hybrid Eulerian-Lagrangian approach and the point droplet approximation. A detailed statistical analysis of both phases is presented. In particular, we show how crucial is the preferential sampling of the vapor phase induced by the inhomogeneous localization of the droplets through the flow. Strong droplet preferential segregation develops suddenly downstream from the inflow section both within the turbulent core and the jet mixing layer. Two distinct mechanisms have been found to drive this phenomenon: the inertial small-scale clustering in the jet core and the intermittent dynamics of droplets across the turbulent-nonturbulent interface in the mixing layer, where dry air entrainment occurs. These phenomenologies strongly affect the overall vaporization process and lead to an impressive widening of the droplet size and vaporization rate distributions in the downstream evolution of the turbulent spray.

We theoretically study the self-propulsion of a thin (slender) colloid driven by asymmetric chemical reactions on its surface at vanishing Reynolds number. Using the method of matched asymptotic expansions, we obtain the colloid self-propulsion velocity as a function of its shape and surface physicochemical properties. The mechanics of self-phoresis for rod-like swimmers has a richer spectrum of behaviors than spherical swimmers due to the presence of two small length scales, the slenderness of the rod and the width of the slip layer. This leads to subtleties in taking the limit of vanishing slenderness. As a result, even for very thin rods, the distribution of curvature along the surface of the swimmer, namely, its shape, plays a surprising role in determining the efficiency of propulsion. We find that thin cylindrical self-phoretic swimmers with blunt ends move faster than thin prolate spheroid shaped swimmers with the same aspect ratio.

The effects of changing Strouhal number on the three-dimensional wake produced by a rigid, bioinspired trapezoidal pitching panel are analyzed through the use of stereoscopic particle image velocimetry over a Strouhal number range of 0.17–0.56. The results show that for all cases, at least some section of the wake comprises an alternating series of interacting vortex rings. The behavior of the flows induced by these vortex rings is consistent with the wake phenomena of spanwise compression and transverse expansion. Increases in Strouhal number correspond to an increased rate of spanwise compression, a greater amount of transverse expansion, and the movement of the location of wake breakdown onset upstream.

We experimentally investigate the bursting dynamics of a confined liquid film suspended in air and find a viscous dynamics distinctly different from the unconfined counterpart due to lack of circular symmetry in the shape of the expanding hole: the confined-viscous bursting proceeds at a constant speed and a rim formed at the bursting tip does not grow. We find a confined-viscous to confined-inertial crossover, as well as an unconfined-inertial to confined-inertial crossover, at which bursting speed does not change although the circular symmetry in the hole shape breaks dynamically.