Researcher profile

Howard A. Stone

Howard A. Stone contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
23works
0followers
11topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

23 published item(s)

preprint2022arXiv

Microswimmers near corrugated, periodic surfaces

We explore hydrodynamic interactions between microswimmers and corrugated, or rough, surfaces, as found often in biological systems and microfluidic devices. Using the Lorentz reciprocal theorem for viscous flows we derive exact expressions for the roughness-induced velocities up to first order in the surface-height fluctuations and provide solutions for the translational and angular velocities valid for arbitrary surface shapes. We apply our theoretical predictions to elucidate the impact of a periodic, wavy surface on the velocities of microswimmers modeled in terms of a superposition of Stokes singularities. Our findings, valid in the framework of a far-field analysis, show that the roughness-induced velocities vary non-monotonically with the wavelength of the surface. For wavelengths comparable to the swimmer-surface distance a pusher can experience a repulsive contribution due to the reflection of flow fields at the edge of a surface cavity, which decreases the overall attraction to the wall.

preprint2022arXiv

Possible, Impossible, and Expected Diameters and Production Rates of Droplets in Aerosols and Sprays

Liquid atomization processes generating sprays and aerosols of droplets are used in many delivery and coating systems involving pure solvents, solutions, and suspensions. In our recent experimental work, we introduced a novel liquid atomization process generating micro-sprays and aerosols of submicron-diameter droplets for pure solvents, solutions, and suspensions: gas jets disintegrate thin liquid films that are formed as bubbles approach a liquid surface. Here we develop a theoretical description of droplet sizes and flow rates, using the first principles of conservation of mass and energy, and accounting for the ratios of specific energies and the ratios of specific energy rates provided by the atomizing gas and dissipated by the atomized liquid. We introduce atomization diagrams as a graphical tool to determine possible, impossible, and expected droplet diameters and specific flow rates in aerosols and sprays produced under various conditions. We find a reasonable agreement between the theory and experiments for five different liquids converted into aerosols of submicron-diameter droplets by an atomization process where gas jets disintegrate thin liquid films that are formed as bubbles approach a liquid surface, and also for five traditional pressure nozzles that produce sprays of droplets of hundreds of microns in diameter. Our study explored the overall range of Ohnesorge number between 0.01-100, and the analysis and atomization diagrams contribute to understanding of liquid atomization and can serve as a theoretical framework for comparing different liquid atomization techniques.

preprint2021arXiv

A Geometric Criterion for the Optimal Spreading of Active Polymers in Porous Media

We perform Brownian dynamics simulations of active stiff polymers undergoing run-reverse dynamics, and so mimic bacterial swimming, in porous media. In accord with recent experiments of \emph{Escherichia coli}, the polymer dynamics are characterized by trapping phases interrupted by directed hopping motion through the pores. We find that the effective translational diffusivities of run-reverse agents can be enhanced up to two orders in magnitude, compared to their non-reversing counterparts, and exhibit a non-monotonic behavior as a function of the reversal rate, which we rationalize using a coarse-grained model. Furthermore, we discover a geometric criterion for the optimal spreading, which emerges when their run lengths are comparable to the longest straight path available in the porous medium. More significantly, our criterion unifies results for porous media with disparate pore sizes and shapes and thus provides a fundamental principle for optimal transport of microorganisms and cargo-carriers in densely-packed biological and environmental settings.

preprint2021arXiv

Perturbing ordinary differential equations to generate resonant and repeated root solutions

In the study of ordinary differential equations (ODEs) of the form $\hat{L}[y(x)]=f(x)$, where $\hat{L}$ is a linear differential operator, two related phenomena can arise: resonance, where $f(x)\propto u(x)$ and $\hat{L}[u(x)]=0$, and repeated roots, where $f(x)=0$ and $\hat{L}=\hat{D}^n$ for $n\geq 2$. We illustrate a method to generate exact solutions to these problems by taking a known homogeneous solution $u(x)$, introducing a parameter $ε$ such that $u(x)\rightarrow u(x;ε)$, and Taylor expanding $u(x;ε)$ about $ε= 0$. The coefficients of this expansion $\frac{\partial^k u}{\partialε^k}\big{|}_{ε=0}$ yield the desired resonant or repeated root solutions to the ODE. This approach, whenever it can be applied, is more insightful and less tedious than standard methods such as reduction of order or variation of parameters. While the ideas can be introduced at the undergraduate level, we could not find any elementary or advanced text that illustrates these ideas with appropriate generality.

preprint2021arXiv

Pressure-driven flow of the viscoelastic Oldroyd-B fluid in narrow non-uniform geometries: analytical results and comparison with simulations

We analyze the pressure-driven flow of a viscoelastic fluid in arbitrarily shaped, narrow channels and present a theoretical framework for calculating the relationship between the flow rate $q$ and pressure drop $Δp$. We utilize the Oldroyd-B model and first identify the characteristic scales and dimensionless parameters governing the flow in the lubrication limit. Employing a perturbation expansion in powers of the Deborah number ($De$), we provide analytical expressions for the velocity, stress, and the $q-Δp$ relation in the weakly viscoelastic limit up to $O(De^2)$. Furthermore, we exploit the reciprocal theorem derived by Boyko $\&$ Stone (Phys. Rev. Fluids, vol. 6, 2021, pp. L081301) to obtain the $q-Δp$ relation at the next order, $O(De^3)$, using only the velocity and stress fields at the previous orders. We validate our analytical results with two-dimensional numerical simulations of the Oldroyd-B fluid in a hyperbolic, symmetric contracting channel and find excellent agreement. For the flow-rate-controlled situation, both our theory and simulations reveal weak dependence of the velocity field on the Deborah number, so that the velocity can be approximated as Newtonian. In contrast to the velocity, the pressure drop strongly depends on the viscoelastic effects and decreases with $De$. Elucidating the relative importance of different terms in the momentum equation contributing to the pressure drop, we identify that a pressure drop reduction for narrow contracting geometries is primarily due to gradients in the viscoelastic shear stresses, while viscoelastic axial stresses have a minor effect on the pressure drop along the symmetry line.

preprint2020arXiv

CO$_2$-driven diffusiophoresis for removal of bacteria

We investigate CO$_2$-driven diffusiophoresis of colloidal particles and bacterial cells in a Hele-Shaw geometry. Combining experiments and a model, we understand the characteristic length and time scales of CO$_2$-driven diffusiophoresis in relation to system dimensions and CO$_2$ diffusivity. Directional migration of wild-type V. cholerae and a mutant lacking flagella, as well as S. aureus and P. aeruginosa, near a dissolving CO$_2$ source shows that diffusiophoresis of bacteria is achieved independent of cell shape and Gram stain. Long-time experiments suggest possible applications for bacterial diffusiophoresis to cleaning systems or anti-biofouling surfaces.

preprint2020arXiv

Harnessing elasticity to generate self-oscillation via an electrohydrodynamic instability

Under a steady DC electric field of sufficient strength, a weakly conducting dielectric sphere in a dielectric solvent with higher conductivity can undergo spontaneous spinning (Quincke rotation) through a pitchfork bifurcation. We design an object composed of a dielectric sphere and an elastic filament. By solving an elasto-electro-hydrodynamic (EEH) problem numerically, we uncover an EEH instability exhibiting diverse dynamic responses. Varying the bending stiffness of the filament, the composite object displays three behaviours: a stationary state, undulatory swimming and steady spinning, where the swimming results from a self-oscillatory instability through a Hopf bifurcation. By conducting a linear stability analysis incorporating an elastohydrodynamic model, we theoretically predict the growth rates and critical conditions, which agree well with the numerical counterparts. We also propose a reduced model system consisting of a minimal elastic structure which reproduces the EEH instability. The elasto-viscous response of the composite structure is able to transform the pitchfork bifurcation into a Hopf bifurcation, leading to self-oscillation. Our results imply a new way of harnessing elastic media to engineer self-oscillations, and more generally, to manipulate and diversify the bifurcations and the corresponding instabilities. These ideas will be useful in designing soft, environmentally adaptive machines.

preprint2020arXiv

Particle motion nearby rough surfaces

We study the hydrodynamic coupling between particles and solid, rough boundaries characterized by random surface textures. Using the Lorentz reciprocal theorem, we derive analytical expressions for the grand mobility tensor of a spherical particle and find that roughness-induced velocities vary nonmonotonically with the characteristic wavelength of the surface. In contrast to sedimentation near a planar wall, our theory predicts continuous particle translation transverse and perpendicular to the applied force. Most prominently, this motion manifests itself in a variance of particle displacements that grows quadratically in time along the direction of the force. This increase is rationalized by surface roughness generating particle sedimentation closer to or farther from the surface, which entails a significant variability of settling velocities.

preprint2020arXiv

Speech can produce jet-like transport relevant to asymptomatic spreading of virus

Many scientific reports document that asymptomatic and presymptomatic individuals contribute to the spread of COVID-19, probably during conversations in social interactions. Droplet emission occurs during speech, yet few studies document the flow to provide the transport mechanism. This lack of understanding prevents informed public health guidance for risk reduction and mitigation strategies, e.g. the &#34;six-foot rule&#34;. Here we analyze flows during breathing and speaking, including phonetic features, using order-of-magnitudes estimates, numerical simulations, and laboratory experiments. We document the spatio-temporal structure of the expelled air flow. Phonetic characteristics of plosive sounds like &#39;P&#39; lead to enhanced directed transport, including jet-like flows that entrain the surrounding air. We highlight three distinct temporal scaling laws for the transport distance of exhaled material including (i) transport over a short distance ($<$ 0.5 m) in a fraction of a second, with large angular variations due to the complexity of speech, (ii) a longer distance, approximately 1 m, where directed transport is driven by individual vortical puffs corresponding to plosive sounds, and (iii) a distance out to about 2 m, or even further, where sequential plosives in a sentence, corresponding effectively to a train of puffs, create conical, jet-like flows. The latter dictates the long-time transport in a conversation. We believe that this work will inform thinking about the role of ventilation, aerosol transport in disease transmission for humans and other animals, and yield a better understanding of linguistic aerodynamics, i.e., aerophonetics.

preprint2020arXiv

Stability of force-driven shear flows in nonequilibrium molecular simulations with periodic boundaries

We analyze the hydrodynamic stability of force-driven parallel shear flows in nonequilibrium molecular simulations with three-dimensional periodic boundary conditions. We show that flows simulated in this way can be linearly unstable, and we derive an expression for the critical Reynolds number as a function of the geometric aspect ratio of the simulation domain. Approximate periodic extensions of Couette and Poiseuille flows are unstable at Reynolds numbers two orders of magnitude smaller than their aperiodic equivalents because the periodic boundaries impose fundamentally different constraints on the flow. This instability has important implications for simulating shear rheology and for designing nonequilibrium simulation methods that are compatible with periodic boundary conditions.

preprint2019arXiv

Non-uniform growth and surface friction determine bacterial biofilm morphology on soft substrates

During development, organisms acquire three-dimensional shapes with important physiological consequences. While the basic mechanisms underlying morphogenesis are known in eukaryotes, it is often difficult to manipulate them in vivo. To circumvent this issue, here we present a study of developing Vibrio cholerae biofilms grown on agar substrates in which the spatiotemporal morphological patterns were altered by varying the agar concentration. Expanding biofilms are initially flat, but later experience a mechanical instability and become wrinkled. Whereas the peripheral region develops ordered radial stripes, the central region acquires a zigzag herringbone-like wrinkle pattern. Depending on the agar concentration, the wrinkles initially appear either in the peripheral region and propagate inward (low agar concentration) or in the central region and propagate outward (high agar concentration). To understand these experimental observations, we developed a model that considers diffusion of nutrients and their uptake by bacteria, bacterial growth/biofilm matrix production, mechanical deformation of both the biofilm and the agar, and the friction between them. Our model demonstrates that depletion of nutrients beneath the central region of the biofilm results in radially-dependent growth profiles, which in turn, produce anisotropic stresses that dictate the morphology of wrinkles. Furthermore, we predict that increasing surface friction (agar concentration) reduces stress anisotropy and shifts the location of the maximum compressive stress, where the wrinkling instability first occurs, toward the center of the biofilm, in agreement with our experimental observations. Our results are broadly applicable to bacterial biofilms with similar morphologies and also provide insight into how other bacterial biofilms form distinct wrinkle patterns.

preprint2019arXiv

Symmetrization of Thin Free-Standing Liquid Films via Capillary-Driven Flow

We present experiments to study the relaxation of a nano-scale cylindrical perturbation at one of the two interfaces of a thin viscous free-standing polymeric film. Driven by capillarity, the film flows and evolves towards equilibrium by first symmetrizing the perturbation between the two interfaces, and eventually broadening the perturbation. A full-Stokes hydrodynamic model is presented which accounts for both the vertical and lateral flows, and which highlights the symmetry in the system. The symmetrization time is found to depend on the membrane thickness, surface tension, and viscosity.

preprint2017arXiv

Collective force generation by groups of migrating bacteria

From biofilm and colony formation in bacteria to wound healing and embryonic development in multicellular organisms, groups of living cells must often move collectively. While considerable study has probed the biophysical mechanisms of how eukaryotic cells generate forces during migration, little such study has been devoted to bacteria, in particular with regard to the question of how bacteria generate and coordinate forces during collective motion. This question is addressed here for the first time using traction force microscopy. We study two distinct motility mechanisms of Myxococcus xanthus, namely twitching and gliding. For twitching, powered by type-IV pilus retraction, we find that individual cells exert local traction in small hotspots with forces on the order of 50 pN. Twitching of bacterial groups also produces traction hotspots, however with amplified forces around 100 pN. Although twitching groups migrate slowly as a whole, traction fluctuates rapidly on timescales <1.5 min. Gliding, the second motility mechanism, is driven by lateral transport of substrate adhesions. When cells are isolated, gliding produces low average traction on the order of 1 Pa. However, traction is amplified in groups by a factor of ~5. Since advancing protrusions of gliding cells push on average in the direction of motion, we infer a long-range compressive load sharing among sub-leading cells. Together, these results show that the forces generated during twitching and gliding have complementary characters and both forces are collectively amplified in groups.

preprint2014arXiv

A pinned or free-floating plate on a thin viscous film. Part 1: A rigid plate

A pinned or free-floating rigid plate lying on the free surface of a thin film of viscous fluid, which itself lies on top of a horizontal substrate that is moving to the right at a constant speed is considered. The focus of the present work is to describe how the competing effects of the speed of the substrate, surface tension, viscosity, and, in the case of a pinned plate, the prescribed pressure in the reservoir of fluid at its upstream end, determine the possible equilibrium positions of the plate, the free surface, and the flow within the film. The present problems are of interest both in their own right as paradigms for a range of fluid-structure interaction problems in which viscosity and surface tension both play an important role, and as a first step towards the study of elastic effects.

preprint2013arXiv

Buckling dynamics of a solvent-stimulated stretched elastomeric sheet

When stretched uniaxially, a thin elastic sheet may exhibit buckling. The occurrence of buckling depends on the geometrical properties of the sheet and the magnitude of the applied strain. Here we show that an elastomeric sheet initially stable under uniaxial stretching can destabilize when exposed to a solvent that swells the elastomer. We demonstrate experimentally and computationally that the features of the buckling pattern depend on the magnitude of stretching, and this observation offers a new way for controlling the shape of a swollen homogeneous thin sheet.

preprint2012arXiv

Source-like solution for radial imbibition into a homogeneous semi-infinite porous medium

We describe the imbibition process from a point source into a homogeneous semi-infinite porous material. When body forces are negligible, the advance of the wetting front is driven by capillary pressure and resisted by viscous forces. With the assumption that the wetting front assumes a hemispherical shape, our analytical results show that the absorbed volume flow rate is approximately constant with respect to time, and that the radius of the wetting evolves in time as r \approx t^(1/3). This cube-root law for the long-time dynamics is confirmed by experiments using a packed cell of glass microspheres with average diameter of 42 μm. This result complements the classical one-dimensional imbibition result where the imbibition length l \approx t^(1/2), and studies in axisymmetric porous cones with small opening angles where l \approx t^(1/4) at long times.

preprint2011arXiv

Dynamical thickening transition in plate coating with concentrated surfactant solutions

We present a large range of experimental data concerning the influence of surfactants on the well-known Landau-Levich-Derjaguin experiment where a liquid film is generated by pulling a solid plate out of a bath. The thickness h of the film was measured as a function of the pulling velocity V for different kind of surfactant and at various concentrations. Measuring the thickening factor $α=h/h_{LLD}$, where hLLD is obtained for a pure liquid, in a wide range of capillary ($Ca=ηV/γ$), two regimes of constant thickening can be identified: at small capillary number, $α$ is large due to a confinement and surface elasticity (or Marangoni) effects and at large Ca, $α$ is slightly higher than unity, due to surface viscous effects. At intermediate Ca, $α$ decreases as Ca increases along a &#34;dynamic transition&#34;. In the case of non-ionic surfactants, the dynamic transition occurs at a fixed Ca, independently of the surfactant concentration, while for ionic surfactants, the dynamic transition depends on the concentration due to the existence of an electrostatic barrier. The control of physico-chemical parameters allowed us to elucidate the nature of the dynamic transition and to relate it to surface rheology.

preprint2010arXiv

Semi-permeable vesicles composed of natural clay

We report a simple route to form robust, inorganic, semi-permeable compartments composed of montmorillonite, a natural plate-like clay mineral that occurs widely in the environment. Mechanical forces due to shear in a narrow gap assemble clay nanoplates from an aqueous suspension onto air bubbles. Translucent vesicles suspended in a single-phase liquid are produced when the clay-covered air bubbles are exposed to a variety of water-miscible organic liquids. These vesicles of clay are mechanically robust and are stable in water and other liquids. The formation of clay vesicles can be described by a physical mechanism that recognizes changes in the wetting characteristics of clay-covered air bubbles in organic liquids. The clay vesicles are covered with small pores and so intrinsically exhibit size-selective permeability, which allows spontaneous compartmentalization of self-assembling molecules in aqueous environments. The results we report here expand our understanding of potential paths to micro-compartmentalization in natural settings and are of relevance to theories of colloidal aggregation, mineral cycles, and the origins of life.

preprint2010arXiv

Surfactant-assisted spreading of an oil-in-water emulsion on the surface of a liquid bath

This fluid dynamics video shows how an oil-in-water emulsion stabilized by an ionic surfactant spreads on the free surface of a layer of pure water. The spreading shows two intriguing features: a transparent area surrounding the source of oil droplets, and a fast retraction of the layer of oil droplets on itself once the source has emptied. We show that the dynamics of spreading are strongly connected to the interfacial/bulk properties of the surfactant.

preprint2007arXiv

Dissolution Arrest and Stability of Armored Bubbles

Dissolving armored bubbles stabilize with nonspherical shapes by jamming the initially Brownian particles adsorbed on their interfaces. In a gas-saturated solution, these shapes are characterized by planar facets or folds for decreasing ratios of the particle to bubble radii. We perform numerical simulations that mimic dissolution, and show that the faceted shape represents a local minimum of energy during volume reduction. This minimum is marked by the vanishing of the Laplace overpressure $ΔP$, which together with the existence of a $V$-interval where $dΔP/dV>0$ guarantees stability against dissolution. The reduction of $ΔP$ is due to the saddle-shape deformation of most of the interface which accompanies the reduction in the mean curvature of the interface.

preprint2006arXiv

Mechanics of Interfacial Composite Materials

Recent experiments and simulations have demonstrated that particle-covered interfaces can exist in stable non-spherical shapes as a result of the steric jamming of the interfacially trapped particles, which confers the interface with solid-like properties. We provide an experimental and theoretical characterization of the mechanical properties of these armored objects, with attention given to the two-dimensional granular state of the interface. Small inhomogeneous stresses produce a plastic response while homogeneous stresses produce a weak elastic response. Shear-driven particle-scale rearrangements explain the basic threshold needed to obtain the near-perfect plastic deformation that is observed. Furthermore, the inhomogeneous stress state of the interface is exhibited experimentally by using surfactants to destabilize the particles on the surface. Since the interfacially trapped particles retain their individual characteristics, armored interfaces can be recognized as a kind of composite material with distinct chemical, structural and mechanical properties.

preprint2006arXiv

Microstructure, morphology and lifetime of armored bubbles exposed to surfactants

We report the behavior of particle-stabilized bubbles (armored bubbles) when exposed to various classes and concentrations of surfactants. The bubbles are non-spherical due to the jamming of the particles on the interface and are stable to dissolution prior to the addition of surfactant. We find that the dissolving bubbles exhibit distinct morphological, microstructural, and lifetime changes, which correlate with the concentration of surfactant employed. For low concentrations of surfactant an armored bubble remains non-spherical while dissolving, while for concentrations close to and above the surfactant CMC a bubble reverts to a spherical shape before dissolving. We propose a microstructural interpretation, supported by our experimental observations of particle dynamics on the bubble interface, that recognizes the role of interfacial jamming and stresses in particle stabilization and surfactant-mediated destabilization of armored bubbles.

preprint2006arXiv

Relaxation time of the topological T1 process in a two-dimensional foam

The elementary topological T1 process in a two-dimensional foam corresponds to the &#34;flip&#34; of one soap film with respect to the geometrical constraints. From a mechanical point of view, this T1 process is an elementary relaxation process through which the entire structure of an out-of-equilibrium foam evolves. The dynamics of this elementary relaxation process has been poorly investigated and is generally neglected during simulations of foams. We study both experimentally and theoretically the T1 dynamics in a dry two-dimensional foam. We show that the dynamics is controlled by the surface viscoelastic properties of the soap films (surface shear plus dilatational viscosity, ms+k, and Gibbs elasticity e), and is independent of the shear viscosity of the bulk liquid. Moreover, our approach illustrates that the dynamics of T1 relaxation process provides a convenient tool for measuring the surface rheological properties: we obtained e = 32+/-8 mN/m and ms+k = 1.3+/-0.7 mPa.m.s for SDS, and e = 65+/-12 mN/m and ms+k = 31+/-12 mPa.m.s for BSA, in good agreement with values reported in the literature.