Source author record

William R. Young

William R. Young appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

8works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

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

8 published item(s)

preprint2022arXiv

Stokes drift and its discontents

The Stokes velocity $\mathbf{u}^\mathrm{S}$, defined approximately by Stokes (1847, Trans. Camb. Philos. Soc., 8, 441-455), and exactly via the Generalized Lagrangian Mean, is divergent even in an incompressible fluid. We show that the Stokes velocity can be naturally decomposed into a solenoidal component, $\mathbf{u}^\mathrm{S}_\mathrm{sol}$, and a remainder that is small for waves with slowly varying amplitudes. We further show that $\mathbf{u}^\mathrm{S}_\mathrm{sol}$ arises as the sole Stokes velocity when the Lagrangian mean flow is suitably redefined to ensure its exact incompressibility. The construction is an application of Soward & Roberts's glm theory (2010, J. Fluid Mech., 661, 45-72) which we specialise to surface gravity waves and implement effectively using a Lie series expansion. We further show that the corresponding Lagrangian-mean momentum equation is formally identical to the Craik-Leibovich equation with $\mathbf{u}^\mathrm{S}_\mathrm{sol}$ replacing $\mathbf{u}^\mathrm{S}$, and we discuss the form of the Stokes pumping associated with both $\mathbf{u}^\mathrm{S}$ and $\mathbf{u}^\mathrm{S}_\mathrm{sol}$.

preprint2020arXiv

The Nusselt numbers of horizontal convection

We consider the problem of horizontal convection in which non-uniform buoyancy, $b_{\rm s}(x,y)$, is imposed on the top surface of a container and all other surfaces are insulating. Horizontal convection produces a net horizontal flux of buoyancy, $\mathbf{J}$, defined by vertically and temporally averaging the interior horizontal flux of buoyancy. We show that $\overline{\mathbf{J}\cdot\mathbf{\nabla}b_{\rm s}}=-κ\langle|\boldsymbol{\nabla}b|^2\rangle$; overbar denotes a space-time average over the top surface, angle brackets denote a volume-time average and $κ$ is the molecular diffusivity of buoyancy $b$. This connection between $\mathbf{J}$ and $κ\langle|\boldsymbol{\nabla}b|^2\rangle$ justifies the definition of the horizontal-convective Nusselt number, $Nu$, as the ratio of $κ\langle|\boldsymbol{\nabla}b|^2\rangle$ to the corresponding quantity produced by molecular diffusion alone. We discuss the advantages of this definition of $Nu$ over other definitions of horizontal-convective Nusselt number currently in use. We investigate transient effects and show that $κ\langle|\boldsymbol{\nabla}b|^2\rangle$ equilibrates more rapidly than other global averages, such as the domain averaged kinetic energy and bottom buoyancy. We show that $κ\langle|\boldsymbol{\nabla} b|^2\rangle$ is essentially the volume-averaged rate of Boussinesq entropy production within the enclosure. In statistical steady state, the interior entropy production is balanced by a flux of entropy through the top surface. This leads to an equivalent "surface Nusselt number", defined as the surface average of vertical buoyancy flux through the top surface times the imposed surface buoyancy $b_{\rm s}(x,y)$. In experiments it is likely easier to evaluate the surface entropy flux, rather than the volume integral of $|\mathbf{\nabla}b|^2$ demanded by $κ\langle|\mathbf{\nabla}b|^2\rangle$.

preprint2020arXiv

Wave-averaged balance: a simple example

In the presence of inertia-gravity waves, the geostrophic and hydrostatic balance that characterises the slow dynamics of rapidly rotating, strongly stratified flows holds in a time-averaged sense and applies to the Lagrangian-mean velocity and buoyancy. We give an elementary derivation of this wave-averaged balance and illustrate its accuracy in numerical solutions of the three-dimensional Boussinesq equations, using a simple configuration in which vertically planar near-inertial waves interact with a barotropic anticylonic vortex. We further use the conservation of the wave-averaged potential vorticity to predict the change in the barotropic vortex induced by the waves.

preprint2019arXiv

Penetration of wind-generated near-inertial waves into a turbulent ocean

An idealized storm scenario is examined in which a wind-generated inertial wave interacts with a turbulent baroclinic quasi-geostrophic flow. The flow is initialized by spinning up a Eady model with a realistic stratification profile. The storm is modeled as an initial value problem for a mixed-layer confined, horizontally-uniform inertial oscillation. The primordial inertial oscillation then evolves under the effects of advection, refraction, dispersion and dissipation. Waves feedback onto the flow by modifying its potential vorticity. In the first few days, refraction dominates and wave energy is attracted (repelled) by regions of negative (positive) vorticity. Wave energy is subsequently drained down anticyclonic pipes. This drainage halts as wave energy encounters weakening vorticity. After a week or two, wave energy accumulates at the bottom of negative vorticity features, i.e. along filamentary structures at shallow depths and in larger anticyclones at greater depths. Wave feedback tends to weaken vortices and thus slow down wave penetration. This effect, however, is found to be weak even for vigorous storms.

preprint2015arXiv

On Gakerkin approximations for the surface-active quasigeostrophic equations

We study the representation of solutions of the three-dimensional quasigeostrophic (QG) equations using Galerkin series with standard vertical modes, with particular attention to the incorporation of active surface buoyancy dynamics. We extend two existing Galerkin approaches (A and B) and develop a new Galerkin approximation (C). Approximation A, due to \cite{flierl1978}, represents the streamfunction as a truncated Galerkin series and defines the potential vorticity (PV) that satisfies the inversion problem exactly. Approximation B, due to \cite{tulloch_smith2009b}, represents the PV as a truncated Galerkin series and calculates the streamfunction that satisfies the inversion problem exactly. Approximation C, the true Galerkin approximation for the QG equations, represents both streamfunction and PV as truncated Galerkin series, but does not satisfy the inversion equation exactly. The three approximations are fundamentally different unless the boundaries are isopycnal surfaces. We discuss the advantages and limitations of approximations A, B, and C in terms of mathematical rigor and conservation laws, and illustrate their relative efficiency by solving linear stability problems with nonzero surface buoyancy. With moderate number of modes, B and C have have superior accuracy than A at high wavenumbers. Because B lacks conservation of energy, we recommend approximation C for constructing solutions to the surface-active QG equations using Galerkin series with standard vertical modes.

preprint2014arXiv

Mixing by microorganisms in stratified fluids

We examine the vertical mixing induced by the swimming of microorganisms at low Reynolds and Péclet numbers in a stably stratified ocean, and show that the global contribution of oceanic microswimmers to vertical mixing is negligible. We propose two approaches to estimating the mixing efficiency, $η$, or the ratio of the rate of potential energy creation to the total rate-of-working on the ocean by microswimmers. The first is based on scaling arguments and estimates $η$ in terms of the ratio between the typical organism size, $a$, and an intrinsic length scale for the stratified flow, $\ell = \left ( νκ/ N^2 \right )^{1/4}$, where $ν$ is the kinematic viscosity, $κ$ the diffusivity, and $N$ the buoyancy frequency. In particular, for small organisms in the relevant oceanic limit, $a / \ell \ll 1$, we predict the scaling $η\sim (a / \ell)^3$. The second estimate of $η$ is formed by solving the full coupled flow-stratification problem by modeling the swimmer as a regularized force dipole, and computing the efficiency numerically. Our computational results, which are examined for all ratios $a/\ell$, validate the scaling arguments in the limit $a / \ell \ll 1$ and further predict $η\approx 1.2 \left ( a / \ell \right )^3$ for vertical swimming and $η\approx 0.15 \left ( a / \ell \right )^3$ for horizontal swimming. These results, relevant for any stratified fluid rich in biological activity, imply that the mixing efficiency of swimming microorganisms in the ocean is at very most 8\% and is likely smaller by at least two orders of magnitude.

preprint2014arXiv

Refraction of swell by surface currents

Using recordings of swell from pitch-and-roll buoys, we have reproduced the classic observations of long-range surface wave propagation originally made by Munk et al. (1963) using a triangular array of bottom pressure measurements. In the modern data, the direction of the incoming swell fluctuates by about $\pm 10^\circ$ on a time scale of one hour. But if the incoming direction is averaged over the duration of an event then, in contrast with the observations by Munk et al. (1963), the sources inferred by great-circle backtracking are most often in good agreement with the location of large storms on weather maps of the Southern Ocean. However there are a few puzzling failures of great-circle backtracking e.g., in one case, the direct great-circle route is blocked by the Tuamoto Islands and the inferred source falls on New Zealand. Mirages like this occur more frequently in the bottom-pressure observations of Munk et al. (1963), where several inferred sources fell on the Antarctic continent. Using spherical ray tracing we investigate the hypothesis that the refraction of waves by surface currents produces the mirages. With reconstructions of surface currents inferred from satellite altimetry, we show that mesoscale vorticity significantly deflects swell away from great-circle propagation so that the source and receiver are connected by a bundle of many rays, none of which precisely follow a great circle. The $\pm 10^\circ$ directional fluctuations at the receiver result from the arrival of wave packets that have travelled along the different rays within this multipath. The occasional failure of great-circle backtracking, and the associated mirages, probably results from partial topographic obstruction of the multipath, which biases the directional average at the receiver.

preprint2011arXiv

The advection-condensation model and water vapour PDFs

The statistically steady humidity distribution resulting from an interaction of advection, modeled as an uncorrelated random walk of moist parcels on an isentropic surface, and a vapour sink, modeled as immediate condensation whenever the specific humidity exceeds a specified saturation humidity, is explored with theory and simulation. A source supplies moisture at the deep-tropical southern boundary of the domain, and the saturation humidity is specified as a monotonically decreasing function of distance from the boundary. The boundary source balances the interior condensation sink, so that a stationary spatially inhomogeneous humidity distribution emerges. An exact solution of the Fokker-Planck equation delivers a simple expression for the resulting probability density function (PDF) of the water vapour field and also of the relative humidity. This solution agrees completely with a numerical simulation of the process, and the humidity PDF exhibits several features of interest, such as bimodality close to the source and unimodality further from the source. The PDFs of specific and relative humidity are broad and non-Gaussian. The domain averaged relative humidity PDF is bimodal with distinct moist and dry peaks, a feature which we show agrees with middleworld isentropic PDFs derived from the ERA interim dataset.