Search

Ocean turbulence at meso- and submesocales affects the propagation of surface waves through refraction and scattering, inducing spatial modulations in significant wave height (SWH). We develop a theoretical framework that relates these modulations to the current that induces them. We exploit the asymptotic smallness of the ratio of typical current speed to wave group speed to derive a linear map -- the U2H map -- between surface current velocity and SWH anomaly. The U2H map is a convolution, non-local in space, expressible as a product in Fourier space by a factor independent of the magnitude of the wavenumber vector. Analytic expressions of the U2H map show how the SWH responds differently to the vortical and divergent parts of the current, and how the anisotropy of the wave spectrum is key to large current-induced SWH anomalies. We implement the U2H map numerically and test its predictions against WAVEWATCH III numerical simulations for both idealised and realistic current configurations.

The refraction of surface gravity waves by currents leads to spatial modulations in the wave field and, in particular, in the significant wave height. We examine this phenomenon in the case of waves scattered by a localised current feature, assuming (i) the smallness of the ratio between current velocity and wave group speed, and (ii) a swell-like, highly directional wave spectrum. We apply matched asymptotics to the equation governing the conservation of wave action in the four-dimensional position--wavenumber space. The resulting explicit formulas show that the modulations in wave action and significant wave height past the localised current are controlled by the vorticity of the current integrated along the primary direction of the swell. We assess the asymptotic predictions against numerical simulations using WAVEWATCH III for a Gaussian vortex. We also consider vortex dipoles to demonstrate the possibility of `vortex cloaking' whereby certain currents have (asymptotically) no impact on the significant wave height. We discuss the role of the ratio of the two small parameters characterising assumptions (i) and (ii) above and show that caustics are only significant for unrealistically large values of this ratio, corresponding to unrealistically narrow directional spectra.

We report a numerical study of horizontal convection (HC) at Prandtl number $Pr = 1$, with both no-slip and free-slip boundary conditions. We obtain 2D and 3D solutions and determine the relation between the Rayleigh number $Ra$ and the Nusselt number $Nu$. In 2D we vary $Ra$ between $0$ and $10^{14}$. In the range $10^6 \le Ra \le 10^{10}$ the $Nu$-$Ra$ relation is $Nu \sim Ra^{1/5}$. With $Ra$ greater than about $10^{11}$ we find a 2D regime with $Nu \sim Ra^{1/4}$ over three decades, up to the highest 2D $Ra$. In 3D, with maximum $Ra = 10^{11.5}$, we find only $Nu \sim Ra^{1/5}$. These results apply to both free slip and no slip boundary conditions. The $Nu \sim Ra^{1/4}$ regime has a double boundary layer (BL): there is a thin BL with thickness $\sim Ra^{-1/4}$ nested inside a thicker BL with thickness $\sim Ra^{-1/5}$. The $Ra^{-1/4}$ BL thickness, which determines $Nu$, coincides with the Kolmogorov and Batchelor scales of HC. Numerical and theoretical results indicate that 3D HC is qualitatively and quantitatively similar to 2D HC. At the same $Ra$, the 3D $Nu$ exceeds the 2D $Nu$ by less than $20$%, i.e., there is very little 3D enhancement of heat transport. Boundary conditions are more important than dimensionality: the 2D free-slip solutions have larger $Nu$ than 3D no-slip solutions. Using the mechanical energy power integral of HC we show that the mean square vorticity of 3D HC is nearly equal to that of 2D HC at the same $Ra$. Thus vorticity amplification by strain-mediated vortex stretching does not operate in 3D HC., Comment: 32 pages, 12 figures, submitted to the Journal of Fluid Mechanics

The coexistence of multiple phytoplankton species despite their reliance on similar resources is often explained with mean-field models assuming mixed populations. In reality, observations of phytoplankton indicate spatial aggregation at all scales, including at the scale of a few individuals. Local spatial aggregation can hinder competitive exclusion since individuals then interact mostly with other individuals of their own species, rather than competitors from different species. To evaluate how microscale spatial aggregation might explain phytoplankton diversity maintenance, an individual-based, multispecies representation of cells in a hydrodynamic environment is required. We formulate a three-dimensional and multispecies individual-based model of phytoplankton population dynamics at the Kolmogorov scale. The model is studied through both simulations and the derivation of spatial moment equations, in connection with point process theory. The spatial moment equations show a good match between theory and simulations. We parameterized the model based on phytoplankters' ecological and physical characteristics, for both large and small phytoplankton. Defining a zone of potential interactions as the overlap between nutrient depletion volumes, we show that local species composition -- within the range of possible interactions -- depends on the size class of phytoplankton. In small phytoplankton, individuals remain in mostly monospecific clusters. Spatial structure therefore favours intra- over inter-specific interactions for small phytoplankton, contributing to coexistence. Large phytoplankton cell neighbourhoods appear more mixed. Although some small-scale self-organizing spatial structure remains and could influence coexistence mechanisms, other factors may need to be explored to explain diversity maintenance in large phytoplankton., Comment: 51 pages; 16 figures

The theory of asymptotic complexity provides an approach to characterizing the behavior of programs in terms of bounds on the number of computational steps executed or use of computational resources. We describe work using ACL2 to prove complexity properties of programs implemented in a simple imperative programming language embedding via an operational semantics in ACL2. We simultaneously prove functional properties of a program and its complexity. We illustrate our approach by describing proofs about a binary search algorithm, proving both that it implements binary search on a sorted list and that it is O(log(n)), where n is the length of the list., Comment: In Proceedings ACL2 2022, arXiv:2205.11103