BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

John S. Wettlaufer

John S. Wettlaufer appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.softphysics.flu-dynBiological Physicsphysics.ao-phastro-ph.EPastro-ph.IMastro-ph.SRcond-mat.mtrl-scicond-mat.stat-mechphysics.class-phphysics.geo-ph

Catalog footprint

What is connected

20works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

20 published item(s)

preprint2022arXiv

Premelting controlled active matter in ice

Self-propelled particles can undergo complex dynamics due to a range of bulk and surface interactions. When a particle is embedded in a host solid near its bulk melting temperature, the latter may melt at the surface of the former in a process known as interfacial premelting. The thickness of the melt film depends on the temperature, impurities, material properties and geometry. A temperature gradient is accompanied by a thermomolecular pressure gradient that drives the interfacial liquid from high to low temperatures and hence the particle from low to high temperatures, in a process called thermal regelation. When the host material is ice and the embedded particle is a biological entity, one has a particularly novel form of active matter, which addresses interplay between a wide range of problems, from extremophiles of both terrestrial and exobiological relevance to ecological dynamics in Earth's cryosphere. Of basic importance in all such settings is the combined influence of biological activity and thermal regelation in controlling the redistribution of bioparticles. Therefore, we re-cast this class of regelation phenomena in the stochastic framework of active Ornstein-Uhlenbeck dynamics and make predictions relevant to this and related problems of interest in biological and geophysical problems. We examine how thermal regelation compromises paleoclimate studies in the context of ice core dating and we find that the activity influences particle dynamics during thermal regelation by enhancing the effective diffusion coefficient. Therefore, accurate dating relies on a quantitative treatment of both effects.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Minimal Data Fidelity for Detection of Stellar Features or Companions

Technological advances in instrumentation have led to an exponential increase in exoplanet detection and scrutiny of stellar features such as spots and faculae. While the spots and faculae enable us to understand the stellar dynamics, exoplanets provide us with a glimpse into stellar evolution. While the ubiquity of noise (e.g., telluric, instrumental, or photonic) is unavoidable, combining this with increased spectrographic resolution compounds technological challenges. To account for these noise sources and resolution issues, we use a temporal multifractal framework to study data from the SOAP 2.0 tool, which simulates a stellar spectrum in the presence of a spot, a facula or a planet. Given these controlled simulations, we vary the resolution as well as the signal-to-noise (S/N) ratio to obtain a lower limit on the resolution and S/N required to robustly detect features. We show that a spot and a facula with a 1% coverage of the stellar disk can be robustly detected for a S/N (per pixel) of 35 and 60 respectively, for any spectral resolution above 20,000, while a planet with a radial velocity (RV) of 10 m/s can be detected for a S/N (per pixel) of 600. Rather than viewing noise as an impediment, our approach uses noise as a source of information.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Impurity effects in thermal regelation

When a particle is placed in a material with a lower bulk melting temperature, intermolecular forces can lead to the existence of a premelted liquid film of the lower melting temperature material. Despite the system being below the melting temperatures of both solids, the liquid film is a consequence of thermodynamic equilibrium, controlled by intermolecular, ionic and other interactions. An imposed temperature gradient drives the translation of the particle by a process of melting and refreezing known as thermal regelation. We calculate the rate of regelation of spherical particles surrounded by premelted films that contain ionic impurities. The impurities enhance the rate of motion thereby influencing the dynamics of single particles and distributions of particles, which we describe in addition to the consequences in natural and technological settings.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

The combined effects of shear and buoyancy on phase boundary stability

We study the effects of externally imposed shear and buoyancy driven flows on the stability of a solid-liquid interface. By reanalyzing the data of Gilpin \emph{et al.} [\emph{J. Fluid Mech.}, {\bf 99}(3), 619 (1980)] we show that the instability of the ice-water interface observed in their experiments was affected by buoyancy effects, and that their velocity measurements are more accurately described by Monin-Obukhov theory. A linear stability analysis of shear and buoyancy driven flow of melt over its solid phase shows that buoyancy is the only destabilizing factor and that the regime of shear flow here, by inhibiting vertical motions and hence the upward heat flux, stabilizes the system. It is also shown that all perturbations to the solid-liquid interface decay at a very modest strength of the shear flow. However, at much larger shear, where flow instabilities coupled with buoyancy might enhance vertical motions, a re-entrant instability may arise.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Fluctuation Spectra and Force Generation in Non-equilibrium Systems

Many biological systems are appropriately viewed as passive inclusions immersed in an active bath: from proteins on active membranes to microscopic swimmers confined by boundaries. The non-equilibrium forces exerted by the active bath on the inclusions or boundaries often regulate function, and such forces may also be exploited in artificial active materials. Nonetheless, the general phenomenology of these active forces remains elusive. We show that the fluctuation spectrum of the active medium, the partitioning of energy as a function of wavenumber, controls the phenomenology of force generation. We find that for a narrow, unimodal spectrum, the force exerted by a non-equilibrium system on two embedded walls depends on the width and the position of the peak in the fluctuation spectrum, and oscillates between repulsion and attraction as a function of wall separation. We examine two apparently disparate examples: the Maritime Casimir effect and recent simulations of active Brownian particles. A key implication of our work is that important non-equilibrium interactions are encoded within the fluctuation spectrum. In this sense the noise becomes the signal.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Large deformations of a soft porous material

Compressing a porous material will decrease the volume of the pore space, driving fluid out. Similarly, injecting fluid into a porous material can expand the pore space, distorting the solid skeleton. This poromechanical coupling has applications ranging from tissue mechanics to hydrogeology. The classical theory of linear poroelasticity captures this coupling by combining Darcy's law with Terzaghi's effective stress and linear elasticity in a linearized kinematic framework. Linear poroelasticity is a good model for very small deformations, but it becomes increasingly inappropriate for moderate to large deformations, which are common in the context of phenomena such as swelling and damage, and for soft materials such as gels and tissues. The well-known theory of large-deformation poroelasticity combines Darcy's law with Terzaghi's effective stress and nonlinear elasticity in a rigorous kinematic framework. This theory has been used extensively in biomechanics to model large elastic deformations in soft tissues, and in geomechanics to model large elastoplastic deformations in soils. Here, we first provide an overview and discussion of this theory with an emphasis on the physics of poromechanical coupling. We present the large-deformation theory in an Eulerian framework to minimize the mathematical complexity, and we show how this nonlinear theory simplifies to linear poroelasticity under the assumption of small strain. We then compare the predictions of linear poroelasticity with those of large-deformation poroelasticity in the context of two uniaxial model problems: Fluid outflow driven by an applied mechanical load (the consolidation problem) and compression driven by a steady fluid throughflow. We use these problems to explore the steady and dynamical errors associated with the linear model in both situations, as well as the impact of introducing a deformation-dependent permeability.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Statistical Mechanics and the Climatology of the Arctic Sea Ice Thickness Distribution

We study the seasonal changes in the thickness distribution of Arctic sea ice, $g(h)$, under climate forcing. Our analytical and numerical approach is based on a Fokker-Planck equation for $g(h)$ (Toppaladoddi \& Wettlaufer \emph{Phys. Rev. Lett.} {\bf 115}, 148501, 2015), in which the thermodynamic growth rates are determined using observed climatology. In particular, the Fokker-Planck equation is coupled to the observationally consistent thermodynamic model of Eisenman \& Wettlaufer (\emph{Proc. Natl. Acad. Sci. USA} {\bf 106}, pp. 28-32, 2009). We find that due to the combined effects of thermodynamics and mechanics, $g(h)$ spreads during winter and contracts during summer. This behavior is in agreement with recent satellite observations from CryoSat-2 (Kwok \& Cunningham, \emph{Phil. Trans. R. Soc. A} {\bf 373}, 20140157, 2015). Because $g(h)$ is a probability density function, we quantify all of the key moments (e.g., mean thickness, fraction of thin/thick ice, mean albedo, relaxation time scales) as greenhouse-gas radiative forcing, $ΔF_0$, increases. The mean ice thickness decays exponentially with $ΔF_0$, but {\em much slower} than do solely thermodynamic models. This exhibits the crucial role that ice mechanics plays in maintaining the ice cover, by redistributing thin ice to thick ice--far more rapidly than can thermal growth alone.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Fluid-driven deformation of a soft granular material

Compressing a porous, fluid-filled material will drive the interstitial fluid out of the pore space, as when squeezing water out of a kitchen sponge. Inversely, injecting fluid into a porous material can deform the solid structure, as when fracturing a shale for natural gas recovery. These poromechanical interactions play an important role in geological and biological systems across a wide range of scales, from the propagation of magma through the Earth's mantle to the transport of fluid through living cells and tissues. The theory of poroelasticity has been largely successful in modeling poromechanical behavior in relatively simple systems, but this continuum theory is fundamentally limited by our understanding of the pore-scale interactions between the fluid and the solid, and these problems are notoriously difficult to study in a laboratory setting. Here, we present a high-resolution measurement of injection-driven poromechanical deformation in a system with granular microsctructure: We inject fluid into a dense, confined monolayer of soft particles and use particle tracking to reveal the dynamics of the multi-scale deformation field. We find that a continuum model based on poroelasticity theory captures certain macroscopic features of the deformation, but the particle-scale deformation field exhibits dramatic departures from smooth, continuum behavior. We observe particle-scale rearrangement and hysteresis, as well as petal-like mesoscale structures that are connected to material failure through spiral shear banding.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Interfacial tension and a three-phase generalized self-consistent theory of non-dilute soft composite solids

In the dilute limit Eshelby's inclusion theory captures the behavior of a wide range of systems and properties. However, because Eshelby's approach neglects interfacial stress, it breaks down in soft materials as the inclusion size approaches the elastocapillarity length $L$. Here, we use a three-phase generalized self-consistent method to calculate the elastic moduli of composites comprised of an isotropic, linear-elastic compliant solid hosting a spatially random monodisperse distribution of spherical liquid droplets. As opposed to similar approaches, we explicitly capture the liquid-solid interfacial stress when it is treated as an isotropic, strain-independent surface tension. Within this framework, the composite stiffness depends solely on the ratio of the elastocapillarity length $L$ to the inclusion radius $R$. Independent of inclusion volume fraction, we find that the composite is stiffened by the inclusions whenever $R < 3L/2$. Over the same range of parameters, we compare our results with alternative approaches (dilute and Mori-Tanaka theories that include surface tension). Our framework can be easily extended to calculate the composite properties of more general soft materials where surface tension plays a role.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the nature of the sea ice albedo feedback in simple models

We examine the nature of the ice-albedo feedback in a long standing approach used in the dynamic-thermodynamic modeling of sea ice. The central issue examined is how the evolution of the ice area is treated when modeling a partial ice cover using a two-category-thickness scheme; thin sea ice and open water in one category and "thick" sea ice in the second. The problem with the scheme is that the area-evolution is handled in a manner that violates the basic rules of calculus, which leads to a neglected area-evolution term that is equivalent to neglecting a leading-order latent heat flux. We demonstrate the consequences by constructing energy balance models with a fractional ice cover and studying them under the influence of increased radiative forcing. It is shown that the neglected flux is particularly important in a decaying ice cover approaching the transitions to seasonal or ice-free conditions. Clearly, a mishandling of the evolution of the ice area has leading-order effects on the ice-albedo feedback. Accordingly, it may be of considerable importance to re-examine the relevant climate model schemes and to begin the process of converting them to fully resolve the sea ice thickness distribution in a manner such as remapping, which does not in principle suffer from the pathology we describe.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Surface tension and the Mori-Tanaka theory of non-dilute soft composite solids

Eshelby's theory is the foundation of composite mechanics, allowing calculation of the effective elastic moduli of composites from a knowledge of their microstructure. However it ignores interfacial stress and only applies to very dilute composites -- i.e. where any inclusions are widely spaced apart. Here, within the framework of the Mori-Tanaka multiphase approximation scheme, we extend Eshelby's theory to treat a composite with interfacial stress in the non-dilute limit. In particular we calculate the elastic moduli of composites comprised of a compliant, elastic solid hosting a non-dilute distribution of identical liquid droplets. The composite stiffness depends strongly on the ratio of the droplet size, $R$, to an elastocapillary length scale, $L$. Interfacial tension substantially impacts the effective elastic moduli of the composite when $R/L\lesssim 100$. When $R < 3L/2$ ($R=3L/2$) liquid inclusions stiffen (cloak the far-field signature of) the solid.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Tailoring boundary geometry to optimize heat transport in turbulent convection

By tailoring the geometry of the upper boundary in turbulent Rayleigh-Bénard convection we manipulate the boundary layer -- interior flow interaction, and examine the heat transport using the Lattice Boltzmann method. For fixed amplitude and varying boundary wavelength $λ$, we find that the exponent $β$ in the Nusselt-Rayleigh scaling relation, $Nu-1 \propto Ra^β$, is maximized at $λ\equiv λ_{\text{max}} \approx (2 π)^{-1}$, but decays to the planar value in both the large ($λ\gg λ_{\text{max}}$) and small ($λ\ll λ_{\text{max}}$) wavelength limits. The changes in the exponent originate in the nature of the coupling between the boundary layer and the interior flow. We present a simple scaling argument embodying this coupling, which describes the maximal convective heat flux.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Stiffening solids with liquid inclusions

From bone and wood to concrete and carbon fibre, composites are ubiquitous natural and engineering materials. Eshelby's inclusion theory describes how macroscopic stress fields couple to isolated microscopic inclusions, allowing prediction of a composite's bulk mechanical properties from a knowledge of its microstructure. It has been extended to describe a wide variety of phenomena from solid fracture to cell adhesion. Here, we show experimentally and theoretically that Eshelby's theory breaks down for small liquid inclusions in a soft solid. In this limit, an isolated droplet's deformation is strongly size-dependent with the smallest droplets mimicking the behaviour of solid inclusions. Furthermore, in opposition to the predictions of conventional composite theory, we find that finite concentrations of small liquid inclusions enhance the stiffness of soft solids. A straight-forward extension of Eshelby's theory, accounting for the surface tension of the solid-liquid interface, explains our experimental observations. The counterintuitive effect of liquid-stiffening of solids is expected whenever droplet radii are smaller than an elastocapillary length, given by the ratio of the surface tension to Young's modulus of the solid matrix.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Surface tension and the mechanics of liquid inclusions in compliant solids

Eshelby's theory of inclusions has wide-reaching implications across the mechanics of materials and structures including the theories of composites, fracture, and plasticity. However, it does not include the effects of surface stress, which has recently been shown to control many processes in soft materials such as gels, elastomers and biological tissue. To extend Eshelby's theory of inclusions to soft materials, we consider liquid inclusions within an isotropic, compressible, linear-elastic solid. We solve for the displacement and stress fields around individual stretched inclusions, accounting for the bulk elasticity of the solid and the surface tension (\textit{i.e.} isotropic strain-independent surface stress) of the solid-liquid interface. Surface tension significantly alters the inclusion's shape and stiffness as well as its near- and far-field stress fields. These phenomenon depend strongly on the ratio of inclusion radius, $R$, to an elastocapillary length, $L$. Surface tension is significant whenever inclusions are smaller than $100L$. While Eshelby theory predicts that liquid inclusions generically reduce the stiffness of an elastic solid, our results show that liquid inclusions can actually stiffen a solid when $R<3L/2$. Intriguingly, surface tension cloaks the far-field signature of liquid inclusions when $R=3L/2$. These results are have far-reaching applications from measuring local stresses in biological tissue, to determining the failure strength of soft composites.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Hertz Beyond Belief

We examine the validity of Hertz's linear elastic theory for central collisions using a viscoelastic model. This model explains why Hertz's theory is accurate in predicting the collision time and maximum contact area even when $40\%$ of the kinetic energy is lost due to viscous dissipation. The main reason is that both the collision time and maximum contact area have a very weak dependence on the impact velocity. Moreover, we show that colliding objects exhibit an apparent size dependent yield strength, which results from larger objects dissipating less energy at a given impact velocity.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Patterning droplets with durotaxis

Numerous cell-types have shown a remarkable ability to detect and move along gradients in stiffness of an underlying substrate -- a process known as durotaxis. The mechanisms underlying durotaxis are still unresolved, but generally believed to involve active sensing and locomotion. Here, we show that simple liquid droplets also undergo durotaxis. By modulating substrate stiffness, we obtain fine control of droplet position on soft, flat substrates. Unlike other control mechanisms, droplet durotaxis works without imposing chemical, thermal, electrical or topographical gradients. This enables a new approach to large-scale droplet patterning and is potentially useful for many applications, such as microfluidics, thermal control and microfabrication.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Surface tension and contact with soft elastic solids

Johnson-Kendall-Robert (JKR) theory is the basis of modern contact mechanics. It describes how two deformable objects adhere together, driven by adhesion energy and opposed by elasticity. However, it does not include solid surface tension, which also opposes adhesion by acting to flatten the surface of soft solids. We tested JKR theory to see if solid surface tension affects indentation behaviour. Using confocal microscopy, we characterised the indentation of glass particles into soft, silicone substrates. While JKR theory held for particles larger than a critical, elastocapillary lengthscale, it failed for smaller particles. Instead, adhesion of small particles mimicked the adsorption of particles at a fluid interface, with a size-independent contact angle between the undeformed surface and the particle given by a generalised version of Young's law. A simple theory quantitatively captures this behaviour, and explains how solid surface tension dominates elasticity for small-scale indentation of soft materials.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Finite-sample-size effects on convection in mushy layers

We report theoretical and experimental investigations of the flow instability responsible for the mushy-layer mode of convection and the formation of chimneys, drainage channels devoid of solid, during steady-state solidification of aqueous ammonium chloride. Under certain growth conditions a state of steady mushy-layer growth with no flow is unstable to the onset of convection, resulting in the formation of chimneys. We present regime diagrams to quantify the state of the flow as a function of the initial liquid concentration, the porous-medium Rayleigh number, and the sample width. For a given liquid concentration, increasing both the porous-medium Rayleigh number and the sample width caused the system to change from a stable state of no flow to a different state with the formation of chimneys. Decreasing the concentration ratio destabilized the system and promoted the formation of chimneys. As the initial liquid concentration increased, onset of convection and formation of chimneys occurred at larger values of the porous-medium Rayleigh number, but the critical cell widths for chimney formation are far less sensitive to the liquid concentration. At the highest liquid concentration, the mushy-layer mode of convection did not occur in the experiment. The formation of multiple chimneys and the morphological transitions between these states are discussed. The experimental results are interpreted in terms of a previous theoretical analysis of finite amplitude convection with chimneys, with a single value of the mushy-layer permeability consistent with the liquid concentrations considered in this study.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Nonlinear mushy-layer convection with chimneys: stability and optimal solute fluxes

We model buoyancy-driven convection with chimneys -- channels of zero solid fraction -- in a mushy layer formed during directional solidification of a binary alloy in two-dimensions. A large suite of numerical simulations is combined with scaling analysis in order to study the parametric dependence of the flow. Stability boundaries are calculated for states of finite-amplitude convection with chimneys, which for a narrow domain can be interpreted in terms of a modified Rayleigh number criterion based on the domain width and mushy-layer permeability. For solidification in a wide domain with multiple chimneys, it has previously been hypothesised that the chimney spacing will adjust to optimise the rate of removal of potential energy from the system. For a wide variety of initial liquid concentration conditions, we consider the detailed flow structure in this optimal state and derive scaling laws for how the flow evolves as the strength of convection increases. For moderate mushy-layer Rayleigh numbers these flow properties support a solute flux that increases linearly with Rayleigh number. This behaviour does not persist indefinitely, however, with porosity-dependent flow saturation resulting in sub-linear growth of the solute flux for sufficiently large Rayleigh numbers. Finally, we consider the influence of the porosity dependence of permeability, with a cubic function and a Carmen-Kozeny permeability yielding qualitatively similar system dynamics and flow profiles for the optimal states.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Ice-lens formation and geometrical supercooling in soils and other colloidal materials

We present a new, physically-intuitive model of ice-lens formation and growth during the freezing of soils and other dense, particulate suspensions. Motivated by experimental evidence, we consider the growth of an ice-filled crack in a freezing soil. At low temperatures, ice in the crack exerts large pressures on the crack walls that will eventually cause the crack to split open. We show that the crack will then propagate across the soil to form a new lens. The process is controlled by two factors: the cohesion of the soil, and the geometrical supercooling of the water in the soil; a new concept introduced to measure the energy available to form a new ice lens. When the supercooling exceeds a critical amount (proportional to the cohesive strength of the soil) a new ice lens forms. This condition for ice-lens formation and growth does not appeal to any ad hoc, empirical assumptions, and explains how periodic ice lenses can form with or without the presence of a frozen fringe. The proposed mechanism is in good agreement with experiments, in particular explaining ice-lens pattern formation, and surges in heave rate associated with the growth of new lenses. Importantly for systems with no frozen fringe, ice-lens formation and frost heave can be predicted given only the unfrozen properties of the soil. We use our theory to estimate ice-lens growth temperatures obtaining quantitative agreement with the limited experimental data that is currently available. Finally we suggest experiments that might be performed in order to verify this theory in more detail. The theory is generalizable to complex natural-soil scenarios, and should therefore be useful in the prediction of macroscopic frost heave rates.

0 citations0 reviewsNo code linkedOpen access