Numerical solution of scattering problems using a Riemann--Hilbert formulation

A fast and accurate numerical method for the solution of scalar and matrix Wiener--Hopf problems is presented. The Wiener--Hopf problems are formulated as Riemann--Hilbert problems on the real line, and a numerical approach developed for these problems is used. It is shown that the known far-field behaviour of the solutions can be exploited to construct numerical schemes providing spectrally accurate results. A number of scalar and matrix Wiener--Hopf problems that generalize the classical Sommerfeld problem of diffraction of plane waves by a semi-infinite plane are solved using the approach

Short- and Long- Time Transport Structures in a Three Dimensional Time Dependent Flow

Lagrangian transport structures for three-dimensional and time-dependent fluid flows are of great interest in numerous applications, particularly for geophysical or oceanic flows. In such flows, chaotic transport and mixing can play important environmental and ecological roles, for examples in pollution spills or plankton migration. In such flows, where simulations or observations are typically available only over a short time, understanding the difference between short-time and long-time transport structures is critical. In this paper, we use a set of classical (i.e. Poincar\'e section, Lyapunov exponent) and alternative (i.e. finite time Lyapunov exponent, Lagrangian coherent structures) tools from dynamical systems theory that analyze chaotic transport both qualitatively and quantitatively. With this set of tools we are able to reveal, identify and highlight differences between short- and long-time transport structures inside a flow composed of a primary horizontal contra-rotating vortex chain, small lateral oscillations and a weak Ekman pumping. The difference is mainly the existence of regular or extremely slowly developing chaotic regions that are only present at short time.Comment: 9 pages, 9 figure

The motion of a buoyant vortex filament

We investigate the motion of a thin vortex filament in the presence of buoyancy. The asymptotic model of Moore & Saffman (Phil. Trans. R. Soc. Lond.A, vol. 272, 1972, pp. 403–429) is extended to take account of buoyancy forces in the force balance on a vortex element. The motion of a buoyant vortex is given by the transverse component of force balance, while the tangential component governs the dynamics of the structure in the core. We show that the local acceleration of axial flow is generated by the external pressure gradient due to gravity. The equations are then solved for vortex rings. An analytic solution for a buoyant vortex ring at a small initial inclination is obtained and asymptotically agrees with the literature

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}}=-\kappa\langle|\boldsymbol{\nabla}b|^2\rangle$; overbar denotes a space-time average over the top surface, angle brackets denote a volume-time average and $\kappa$ is the molecular diffusivity of buoyancy $b$. This connection between $\mathbf{J}$ and $\kappa\langle|\boldsymbol{\nabla}b|^2\rangle$ justifies the definition of the horizontal-convective Nusselt number, $Nu$, as the ratio of $\kappa \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 $\kappa \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 $\kappa\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 $\kappa\langle|\mathbf{\nabla}b|^2\rangle$.Comment: 16 pages, 7 figure

An Arabidopsis rhomboid protease has roles in the chloroplast and in flower development

Increasing numbers of cellular pathways are now recognized to be regulated via proteolytic processing events. The rhomboid family of serine proteases plays a pivotal role in a diverse range of pathways, activating and releasing proteins via regulated intramembrane proteolysis. The prototype rhomboid protease, rhomboid-1 in Drosophila, is the key activator of epidermal growth factor (EGF) receptor pathway signalling in the fly and thus affects multiple aspects of development. The role of the rhomboid family in plants is explored and another developmental phenotype, this time in a mutant of an Arabidopsis chloroplast-localized rhomboid, is reported here. It is confirmed by GFP-protein fusion that this protease is located in the envelope of chloroplasts and of chlorophyll-free plastids elsewhere in the plant. Mutant plants lacking this organellar rhomboid demonstrate reduced fertility, as documented previously with KOM—the one other Arabidopsis rhomboid mutant that has been reported in the literature—along with aberrant floral morphology

Enhanced dispersion of near-inertial waves in an idealized geostrophic flow

This paper presents a simplified model of the process through which a geostrophic flow enhances the vertical propagation of near-inertial activity from the mixed layer into the deeper ocean. The geostrophic flow is idealized as steady and barotropic with a sinusoidal dependence on the north-south coordinate; the corresponding streamfunction takes the form ψ = - Ψ cos (2αy). Near-inertial oscillations are considered in linear theory and disturbances are decomposed into horizontal and vertical normal modes. For this particular flow, the horizontal modes are given in terms of Mathieu functions. The initial-value problem can then be solved by projecting onto this set of normal modes. A detailed solution is presented for the case in which the mixed layer is set into motion as a slab. There is no initial horizontal structure in the model mixed layer; rather, horizontal structure, such as enhanced near-inertial energy in regions of negative vorticity, is impressed on the near-inertial fields by the pre-existing geostrophic flow. Many details of the solution, such as the rate at which near-inertial activity in the mixed layer decays, are controlled by the nondimensional number, Y = 4 Ψf0/H2mixN2mix, where f0 is the inertial frequency, Hmix is the mixed-layer depth, and Nmix is the buoyancy frequency immediately below the base of the mixed layer. When Y is large, near-inertial activity in the mixed layer decays on a time-scale HmixNmix/α2Ψ3/2f01/2. When Y is small, near-inertial activity in the mixed layer decays on a time-scale proportional to N2mixH2mix/α2Ψ2f0