Researcher profile

J. S. Wettlaufer

J. S. Wettlaufer contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
0followers
10topics
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

8 published item(s)

preprint2022arXiv

Analytical solution of stochastic resonance in the nonadiabatic regime

We generalize stochastic resonance to the nonadiabatic limit by treating the double-well potential using two quadratic potentials. We use a singular perturbation method to determine an approximate analytical solution for the probability density function that asymptotically connects local solutions in boundary layers near the two minima with those in the region of the maximum that separates them. The validity of the analytical solution is confirmed numerically. Free from the constraints of the adiabatic limit, the approach allows us to predict the escape rate from one stable basin to another for systems experiencing a more complex periodic forcing.

preprint2022arXiv

Asymptotic interpretation of the Miles mechanism of wind-wave instability

When wind blows over water, ripples are generated on the water surface. These ripples can be regarded as perturbations of the wind field, which is modelled as a parallel inviscid flow. For a given wavenumber $k$, the perturbed streamfunction of the wind field and the complex phase speed are the eigenfunction and the eigenvalue of the so-called Rayleigh equation in a semi-infinite domain. Because of the small air-water density ratio, $ρ_{\rm{a}}/ρ_{\rm{w}}\equivε<<1$, the wind and the ripples are weakly coupled, and the eigenvalue problem can be solved perturbatively. At the leading order, the eigenvalue is equal to the phase speed $c_0$ of surface waves. At order $ε$, the eigenvalue has a finite imaginary part, which implies growth. Miles (1957) showed that the growth rate is proportional to the square modulus of the leading-order eigenfunction evaluated at the so-called critical level $z=z_c$, where the wind speed is equal to $c_0$ and the waves extract energy from the wind. Here, we construct uniform asymptotic approximations of the leading-order eigenfunction for long waves, which we use to calculate the growth rate as a function of $k$. In the strong wind limit, we find that the fastest growing wave is such that the aerodynamic pressure is in phase with the wave slope. The results are confirmed numerically.

preprint2022arXiv

Biolocomotion and premelting in ice

Biota are found in glaciers, ice sheets and permafrost. Ice bound micro-organisms evolve in a complex mobile environment facilitated or hindered by a range of bulk and surface interactions. When a particle is embedded in a host solid near its bulk melting temperature, a melted film forms at the surface of the particle in a process known as interfacial premelting. Under a temperature gradient, the particle is driven by a thermomolecular pressure gradient toward regions of higher temperatures in a process called thermal regelation. When the host solid is ice and the particles are biota, thriving in their environment requires the development of strategies, such as producing exopolymeric substances (EPS) and antifreeze glycoproteins (AFP) that enhance the interfacial water. Therefore, thermal regelation is enhanced and modified by a process we term {\em bio-enhanced premelting}. Additionally, the motion of bioparticles is influenced by chemical gradients influenced by nutrients within the icy host body. We show how the overall trajectory of bioparticles is controlled by a competition between thermal regelation and directed biolocomotion. By re-casting this class of regelation phenomena in the stochastic framework of active Ornstein-Uhlenbeck dynamics, and using multiple scales analysis, we find that for an attractive (repulsive) nutrient source, that thermal regelation is enhanced (suppressed) by biolocomotion. This phenomena is important in astrobiology, the biosignatures of extremophiles and in terrestrial paleoclimatology.

preprint2022arXiv

The combined effects of buoyancy, rotation, and shear on phase boundary evolution

We use well-resolved numerical simulations to study the combined effects of buoyancy, pressure-driven shear and rotation on the melt rate and morphology of a layer of pure solid overlying its liquid phase in three dimensions at a Rayleigh number $Ra = 1.25 \times 10^5$. During thermal convection we find that the rate of melting of the solid phase varies non-monotonically with the strength of the imposed shear flow. In the absence of rotation, depending on whether buoyancy or shear dominates the flow, we observe either domes or ridges aligned in the direction of the shear flow respectively. Furthermore, we show that the geometry of the phase boundary has important effects on the magnitude and evolution of the heat flux in the liquid layer. In the presence of rotation, the strength of which is characterized by the Rossby number, Ro, we observe that for Ro = O(1) the mean flow in the interior is perpendicular to the direction of the constant horizontal applied pressure gradient. As the magnitude of this pressure gradient increases, the geometry of solid-liquid interface evolves from the voids characteristic of melting by rotating convection, to grooves oriented perpendicular or obliquely to the direction of the pressure gradient.

preprint2020arXiv

Precursors to Rare Events in Stochastic Resonance

In stochastic resonance, a periodically forced Brownian particle in a double-well potential jumps between minima at rare increments, the prediction of which poses a major theoretical challenge. Here, we use a path-integral method to find a precursor to these transitions by determining the most probable (or &#34;{optimal}&#34;) space-time path of a particle. We characterize the optimal path using a direct comparison principle between the Langevin and Hamiltonian dynamical descriptions, allowing us to express the jump condition in terms of the accumulation of noise around the stable periodic path. In consequence, as a system approaches a rare event these fluctuations approach one of the deterministic minimizers, thereby providing a precursor for predicting a stochastic transition. We demonstrate the method numerically, which allows us to determine whether a state is following a stable periodic path or will experience an incipient jump with a high probability. The vast range of systems that exhibit stochastic resonance behavior insures broad relevance of our framework, which allows one to extract precursor fluctuations from data.

preprint2020arXiv

Transient Convective Spin-up Dynamics

We study the formation, longevity and breakdown of convective rings during impulsive spin-up in square and cylindrical containers using direct numerical simulations. The rings, which are axisymmetric alternating regions of up- and down-welling flow that can last for O (100) rotation times, were first demonstrated experimentally and arise due to a balance between Coriolis and viscous effects. We study the formation of these rings in the context of the Greenspan-Howard spin-up process, the disruption of which modifies ring formation and evolution. We show that, unless imprinted by boundary geometry, convective rings can only form when the surface providing buoyancy forcing is a free-slip surface, thereby explaining an apparent disagreement between experimental results in the literature. For Prandtl numbers from 1--5 we find that the longest-lived rings occur for intermediate Prandtl numbers, with a Rossby number dependence. Finally, we find that the constant evaporative heat-flux conditions imposed in the experiments are essential in sustaining the rings and in maintaining the vortices that form in consequence of the ring breakdown.

preprint2019arXiv

Nonadiabatic escape and stochastic resonance

We analyze the fluctuation-driven escape of particles from a metastable state under the influence of a weak periodic force. We develop an asymptotic method to solve the appropriate Fokker-Planck equation with mixed natural and absorbing boundary conditions. The approach uses two boundary layers flanking an interior region; most of the probability is concentrated within the boundary layer near the metastable point of the potential and particles transit the interior region before exiting the domain through the other boundary layer, which is near the unstable maximal point of the potential. The dominant processes in each region are given by approximate time-dependent solutions matched to construct the approximate composite solution, which gives the rate of escape with weak periodic forcing. Using reflection we extend the method to a double well potential influenced by white noise and weak periodic forcing, and thereby derive a two-state stochastic model--the simplest treatment of stochastic resonance theory--in the nonadiabatic limit

preprint2015arXiv

Maximal Stochastic Transport in the Lorenz Equations

We calculate the stochastic upper bounds for the Lorenz equations using an extension of the background method. In analogy with Rayleigh-Bénard convection the upper bounds are for heat transport versus Rayleigh number. As might be expected, the stochastic upper bounds are larger than the deterministic counterpart of \citet{Doering15}, but their variation with noise amplitude exhibits interesting behavior. Below the transition to chaotic dynamics the upper bounds increase monotonically with noise amplitude. However, in the chaotic regime this monotonicity depends on the number of realizations in the ensemble; at a particular Rayleigh number the bound may increase or decrease with noise amplitude. The origin of this behavior is the coupling between the noise and unstable periodic orbits, the degree of which depends on the degree to which the ensemble represents the ergodic set. This is confirmed by examining the close returns plots of the full solutions to the stochastic equations and the numerical convergence of the noise correlations. The numerical convergence of both the ensemble and time averages of the noise correlations is sufficiently slow that it is the limiting aspect of the realization of these bounds. Finally, we note that the full solutions of the stochastic equations demonstrate that the effect of noise is equivalent to the effect of chaos.