DSpace Collection:

On global regularity of 2D generalized magnetohydrodynamic equations

Tran, Chuong Van — Wed, 15 May 2013 00:00:00 GMT

Abstract: In this article we study the global regularity of 2D generalized magnetohydrodynamic equations (2D GMHD), in which the dissipation terms are $- \nu \left( - \triangle \right)^{\alpha} u$ and $- \kappa \left( - \triangle \right)^{\beta} b$. We show that smooth solutions are global in the following three cases: $\alpha \geqslant 1 / 2, \beta \geqslant 1$; $0 \leqslant \alpha < 1 / 2, 2 \alpha + \beta > 2$; $\alpha \geqslant 2, \beta = 0$. We also show that in the inviscid case $\nu = 0$, if $\beta > 1$, then smooth solutions are global as long as the direction of the magnetic field remains smooth enough.

Sharp global nonlinear stability for a fluid overlying a highly porous material

Hill, Antony A. — Fri, 08 Jan 2010 00:00:00 GMT

Abstract: The stability of convection in a two-layer system in which a layer of fluid with a temperature-dependent viscosity overlies and saturates a highly porous material is studied. Owing to the difficulties associated with incorporating the nonlinear advection term in the Navier-Stokes equations into a stability analysis, previous literature on fluid/porous thermal convection has modelled the fluid using the linear Stokes equations. This paper derives global stability for the full nonlinear system, by utilizing a model proposed by Ladyzhenskaya. The nonlinear stability boundaries are shown to be sharp when compared with the linear instability thresholds.

Nonlinear stability of the one-domain approach to modelling convection in superposed fluid and porous layers

Hill, A A — Wed, 01 Sep 2010 00:00:00 GMT

Abstract: Studies of the nonlinear stability of fluid/porous systems have been developed very recently. A two-domain modelling approach has been adopted in previous works, but was restricted to specific configurations. The extension to the more general case of a Navier–Stokes modelled fluid over a porous material was not achieved for the two-domain approach owing to the difficulties associated with handling the interfacial boundary conditions. This paper addresses this issue by adopting a one-domain approach, where the governing equations for both regions are combined into a unique set of equations that are valid for the entire domain. It is shown that the nonlinear stability bound, in the one-domain approach, is very sharp and hence excludes the possibility of subcritical instabilities. Moreover, the one-domain approach is compared with an equivalent two-domain approach, and excellent agreement is found between the two.

Instability in internal solitary waves with trapped cores

Carr, Magda — Sun, 01 Jan 2012 00:00:00 GMT

Abstract: A numerical method that employs a combination of contour advection and pseudo-spectral techniques is used to investigate instability in internal solitary waves with trapped cores. A three-layer configuration for the background stratification in which the top two layers are linearly stratified and the lower layer is homogeneous is considered throughout. The strength of the stratification in the very top layer is chosen to be sufficient so that waves of depression with trapped cores can be generated. The flow is assumed to satisfy the Dubriel-Jacotin-Long equation both inside and outside of the core region. The Brunt-Vaisala frequency is modelled such that it varies from a constant value outside of the core to zero inside the core over a sharp but continuous transition length. This results in a stagnant core in which the vorticity is zero and the density is homogeneous and approximately equal to that at the core boundary. The time dependent simulations show that instability occurs on the boundary of the core. The instability takes the form of Kelvin-Helmholtz billows. If the instability in the vorticity field is energetic enough, disturbance in the buoyancy field is also seen and fluid exchange takes place across the core boundary. Occurrence of the Kelvin-Helmholtz billows is attributed to the sharp change in the vorticity field at the boundary between the core and the pycnocline. The numerical scheme is not limited by small Richardson number unlike the other alternatives currently available in the literature which appear to be.

Shear induced breaking of large amplitude internal solitary waves

Fructus, D — Sun, 01 Feb 2009 00:00:00 GMT

Convectively induced shear instability in large amplitude internal solitary waves

Carr, Magda — Mon, 01 Dec 2008 00:00:00 GMT

Abstract: Laboratory study has been carried out to investigate the instability of an internal solitary wave of depression in a shallow stratified fluid system. The experimental campaign has been supported by theoretical computations and has focused on a two layered stratification consisting of a homogeneous dense layer below a linearly stratified top layer. The initial background stratification has been varied and it is found that the onset and intensity of breaking are affected dramatically by changes in the background stratification. Manifestations of a combination of shear and convective instability are seen on the leading face of the wave. It is shown that there is an interplay between the two instability types and convective instability induces shear by enhancing isopycnal compression. Variation in the upper boundary condition is also found to have an effect on stability. In particular, the implications for convective instability are shown to be profound and a dramatic increase in wave amplitude is seen for a fixed (as opposed to free) upper boundary condition.

Number of degrees of freedom and energy spectrum of surface quasi-geostrophic turbulence

Tran, Chuong V. — Sat, 01 Oct 2011 00:00:00 GMT

Abstract: We study both theoretically and numerically surface quasi-geostrophic turbulence regularized by the usual molecular viscosity, with an emphasis on a number of classical predictions. It is found that the system's number of degrees of freedom N, which is defined in terms of local Lyapunov exponents, scales as Re-3/2, where R e is the Reynolds number expressible in terms of the viscosity, energy dissipation rate and system's integral scale. For general power-law energy spectra k(-alpha), a comparison of N with the number of dynamically active Fourier modes, i.e. the modes within the energy inertial range, yields alpha = 5/3. This comparison further renders the scaling Re-1/2 for the exponential dissipation rate at the dissipation wavenumber. These results have been predicted on the basis of Kolmogorov's theory. Our approach thus recovers these classical predictions and is an analytic alternative to the traditional phenomenological method. The implications of the present findings are discussed in conjunction with related results in the literature. Support for the analytic results is provided through a series of direct numerical simulations. Description: L.A.K.B. was supported by an EPSRC post-graduate studentship.

Coronal heating by the partial relaxation of twisted loops

Bareford, Michael — Fri, 01 Feb 2013 00:00:00 GMT

Abstract: Context: Relaxation theory offers a straightforward method for estimating the energy that is released when a magnetic field becomes unstable, as a result of continual convective driving. Aims: We present new results obtained from nonlinear magnetohydrodynamic (MHD) simulations of idealised coronal loops. The purpose of this work is to determine whether or not the simulation results agree with Taylor relaxation, which will require a modified version of relaxation theory applicable to unbounded field configurations. Methods: A three-dimensional (3D) MHD Lagrangian-remap code is used to simulate the evolution of a line-tied cylindrical coronal loop model. This model comprises three concentric layers surrounded by a potential envelope; hence, being twisted locally, each loop configuration is distinguished by a piecewise-constant current profile. Initially, all configurations carry zero-net-current fields and are in ideally unstable equilibrium. The simulation results are compared with the predictions of helicity conserving relaxation theory. Results: For all simulations, the change in helicity is no more than 2% of the initial value; also, the numerical helicities match the analytically-determined values. Magnetic energy dissipation predominantly occurs via shock heating associated with magnetic reconnection in distributed current sheets. The energy release and final field profiles produced by the numerical simulations are in agreement with the predictions given by a new model of partial relaxation theory: the relaxed field is close to a linear force free state; however, the extent of the relaxation region is limited, while the loop undergoes some radial expansion. Conclusions: The results presented here support the use of partial relaxation theory, specifically, when calculating the heating-event distributions produced by ensembles of kink-unstable loops.

Damping of kink waves by mode coupling : I. Analytical treatment

Hood, Alan William — Tue, 01 Jan 2013 00:00:00 GMT

Abstract: Aims. To investigate the spatial damping of propagating kink waves in an inhomogeneous plasma. In the limit of a thin tube surrounded by a thin transition layer, an analytical formulation for kink waves driven in from the bottom boundary of the corona is presented. Methods. The spatial form for the damping of the kink mode was investigated using various analytical approximations. When the density ratio between the internal density and the external density is not too large, a simple di.erential-integral equation was used. Approximate analytical solutions to this equation are presented. Results. For the first time, the form of the spatial damping of the kink mode is shown analytically to be Gaussian in nature near the driven boundary. For several wavelengths, the amplitude of the kink mode is proportional to (1 + exp(-z2 /L2 g))/2, where L2g = 16/ǫκ2 k2 . Although the actual value of 16 in Lg depends on the particular form of the driver, this form is very general and its dependence on the other parameters does not change. For large distances, the damping profile appears to be roughly linear exponential decay. This is shown analytically by a series expansion when the inhomogeneous layer width is small enough.

The influence of a fluid-porous interface on solar pond stability

Hill, A. A — Fri, 01 Feb 2013 00:00:00 GMT

Abstract: The linear instability of the gradient zone of a solar pond containing a fluidporous interface is investigated. It is found that the gradient zone can retain the same stability for lower values of the solute Rayleigh number with the introduction of a porous material compared with a purely fluid layer, whilst maintaining the same lower convective zone temperature. Interestingly, it is also shown that for certain parameter values the penetration of a porous medium into the gradient zone can cause the temperature of the lower convective zone to rise. However, for certain parameter ranges, when the fluid-porous interface is towards the top of the gradient zone, the solar pond can become highly unstable.

A Bayesian approach to fitting Gibbs processes with temporal random effects

King, Ruth — Sat, 01 Dec 2012 00:00:00 GMT

Abstract: We consider spatial point pattern data that have been observed repeatedly over a period of time in an inhomogeneous environment. Each spatial point pattern can be regarded as a “snapshot” of the underlying point process at a series of times. Thus, the number of points and corresponding locations of points differ for each snapshot. Each snapshot can be analyzed independently, but in many cases there may be little information in the data relating to model parameters, particularly parameters relating to the interaction between points. Thus, we develop an integrated approach, simultaneously analyzing all snapshots within a single robust and consistent analysis. We assume that sufficient time has passed between observation dates so that the spatial point patterns can be regarded as independent replicates, given spatial covariates. We develop a joint mixed effects Gibbs point process model for the replicates of spatial point patterns by considering environmental covariates in the analysis as fixed effects, to model the heterogeneous environment, with a random effects (or hierarchical) component to account for the different observation days for the intensity function. We demonstrate how the model can be fitted within a Bayesian framework using an auxiliary variable approach to deal with the issue of the random effects component. We apply the methods to a data set of musk oxen herds and demonstrate the increased precision of the parameter estimates when considering all available data within a single integrated analysis.

Collisionless distribution function for the relativistic force-free Harris sheet

Stark, C. R. — Sun, 01 Jan 2012 00:00:00 GMT

Abstract: A self-consistent collisionless distribution function for the relativistic analogue of the force-free Harris sheet is presented. This distribution function is the relativistic generalization of the distribution function for the non-relativistic collisionless force-free Harris sheet recently found by Harrison and Neukirch [Phys. Rev. Lett. 102, 135003 (2009)], as it has the same dependence on the particle energy and canonical momenta. We present a detailed calculation which shows that the proposed distribution function generates the required current density profile (and thus magnetic field profile) in a frame of reference in which the electric potential vanishes identically. The connection between the parameters of the distribution function and the macroscopic parameters such as the current sheet thickness is discussed. (C) 2012 American Institute of Physics. [doi: 10.1063/1.3677268]

Numerical simulation of shear-induced instabilities in internal solitary waves

Carr, Magda — Sun, 25 Sep 2011 00:00:00 GMT

Abstract: A numerical method that employs a combination of contour advection and pseudo-spectral techniques is used to simulate shear-induced instabilities in an internal solitary wave (ISW). A three-layer configuration for the background stratification, in which a linearly stratified intermediate layer is sandwiched between two homogeneous ones, is considered throughout. The flow is assumed to satisfy the inviscid, incompressible, Oberbeck–Boussinesq equations in two dimensions. Simulations are initialized by fully nonlinear, steady-state, ISWs. The results of the simulations show that the instability takes place in the pycnocline and manifests itself as Kelvin–Helmholtz billows. The billows form near the trough of the wave, subsequently grow and disturb the tail. Both the critical Richardson number (Ric) and the critical amplitude required for instability are found to be functions of the ratio of the undisturbed layer thicknesses. It is shown, therefore, that the constant, critical bound for instability in ISWs given in Barad & Fringer (J. Fluid Mech., vol. 644, 2010, pp. 61–95), namely Ric = 0.1 ± 0.01 , is not a sufficient condition for instability. It is also shown that the critical value of Lx/λ required for instability, where Lx is the length of the region in a wave in which Ri < 1/4 and λ is the half-width of the wave, is sensitive to the ratio of the layer thicknesses. Similarly, a linear stability analysis reveals that δiTw (where δi is the growth rate of the instability averaged over Tw, the period in which parcels of fluid are subjected to Ri < 1/4) is very sensitive to the transition between the undisturbed pycnocline and the homogeneous layers, and the amplitude of the wave. Therefore, the alternative tests for instability presented in Fructus et al. (J. Fluid Mech., vol. 620, 2009, pp. 1–29) and Barad & Fringer (J. Fluid Mech., vol. 644, 2010, pp. 61–95), respectively, namely Lx/λ ≥ 0.86 and δiTw > 5 , are shown to be valid only for a limited parameter range. Description: This work was supported by the UK Engineering and Physical Sciences Research Council [grant number EP/F030622/1]

Behind and beyond a theorem on groups related to trivalent graphs

Havas, George — Mon, 01 Dec 2008 00:00:00 GMT

Abstract: In 2006 we completed the proof of a five-part conjecture that was made in 1977 about a family of groups related to trivalent graphs. This family covers all 2-generator, 2-relator groups where one relator specifies that a generator is an involution and the other relator has three syllables. Our proof relies upon detailed but general computations in the groups under question. The proof is theoretical, but based upon explicit proofs produced by machine for individual cases. Here we explain how we derived the general proofs from specific cases. The conjecture essentially addressed only the finite groups in the family. Here we extend the results to infinite groups, effectively determining when members of this family of finitely presented groups are simply isomorphic to a specific quotient.

Lower-hybrid waves generated by anomalous Doppler resonance in auroral plasmas

Bingham, R. — Sun, 01 Aug 2010 00:00:00 GMT

Abstract: This paper describes sonic aspects of lower-hybrid wave activity in space plasmas. Lower-hybrid waves are particularly important since they can transfer energy efficiently between electrons and ions in a collisionless magnetized plasma. We consider the 'fan' or anomalous Doppler resonance instability driven by energetic electron tails and show that it is responsible for the generation of lower-hybrid waves. We also demonstrate that observations of their intensity are sufficient to drive the modulational instability.

Falling towards forgetfulness : synaptic decay prevents spontaneous recovery of memory

Stone, James V. — Fri, 22 Aug 2008 00:00:00 GMT

Abstract: Long after a new language has been learned and forgotten, relearning a few words seems to trigger the recall of other words. This "free-lunch learning'' (FLL) effect has been demonstrated both in humans and in neural network models. Specifically, previous work proved that linear networks that learn a set of associations, then partially forget them all, and finally relearn some of the associations, show improved performance on the remaining (i.e., nonrelearned) associations. Here, we prove that relearning forgotten associations decreases performance on nonrelearned associations; an effect we call negative free-lunch learning. The difference between free-lunch learning and the negative free-lunch learning presented here is due to the particular method used to induce forgetting. Specifically, if forgetting is induced by isotropic drifting of weight vectors (i.e., by adding isotropic noise), then free-lunch learning is observed. However, as proved here, if forgetting is induced by weight values that simply decay or fall towards zero, then negative free-lunch learning is observed. From a biological perspective, and assuming that nervous systems are analogous to the networks used here, this suggests that evolution may have selected physiological mechanisms that involve forgetting using a form of synaptic drift rather than synaptic decay, because synaptic drift, but not synaptic decay, yields free-lunch learning.

On the relationship between equilibrium bifurcations and ideal MHD instabilities for line-tied coronal loops

Neukirch, T. — Fri, 01 Jan 2010 00:00:00 GMT

Abstract: For axisymmetric models for coronal loops the relationship between the bifurcation points of magnetohydrodynamic (MHD) equilibrium sequences and the points of linear ideal MHD instability is investigated, imposing line-tied boundary conditions. Using a well-studied example based on the Gold -aEuro parts per thousand Hoyle equilibrium, it is demonstrated that if the equilibrium sequence is calculated using the Grad -aEuro parts per thousand Shafranov equation, the instability corresponds to the second bifurcation point and not the first bifurcation point, because the equilibrium boundary conditions allow for modes which are excluded from the linear ideal stability analysis. This is shown by calculating the bifurcating equilibrium branches and comparing the spatial structure of the solutions close to the bifurcation point with the spatial structure of the unstable mode. If the equilibrium sequence is calculated using Euler potentials, the first bifurcation point of the Grad -aEuro parts per thousand Shafranov case is not found, and the first bifurcation point of the Euler potential description coincides with the ideal instability threshold. An explanation of this results in terms of linear bifurcation theory is given and the implications for the use of MHD equilibrium bifurcations to explain eruptive phenomena is briefly discussed.

Automatic presentations and semigroup constructions

Cain, Alan J. — Sun, 01 Aug 2010 00:00:00 GMT

Abstract: An automatic presentation for a relational structure is, informally, an abstract representation of the elements of that structure by means of a regular language such that the relations can all be recognized by finite automata. A structure admitting an automatic presentation is said to be FA-presentable. This paper studies the interaction of automatic presentations and certain semigroup constructions, namely: direct products, free products, finite Rees index extensions and subsemigroups, strong semilattices of semigroups, Rees matrix semigroups, Bruck-Reilly extensions, zero-direct unions, semidirect products, wreath products, ideals, and quotient semigroups. For each case, the closure of the class of FA-presentable semigroups under that construction is considered, as is the question of whether the FA-presentability of the semigroup obtained from such a construction implies the FA-presentability of the original semigroup[s]. Classifications are also given of the FA-presentable finitely generated Clifford semigroups, completely simple semigroups, and completely 0-simple semigroups.

Automatic presentations for semigroups

Cain, Alan J. — Sun, 01 Nov 2009 00:00:00 GMT

Abstract: This paper applies the concept of FA-presentable structures to semigroups. We give a complete classification of the finitely generated FA-presentable cancellative semigroups: namely, a finitely generated cancellative semigroup is FA-presentable if and only if it is a subsemigroup of a virtually abelian group. We prove that all finitely generated commutative semigroups are FA-presentable. We give a complete list of FA-presentable one-relation semigroups and compare the classes of FA-presentable semigroups and automatic semigroups. (C) 2009 Elsevier Inc. All rights reserved. Description: Special Issue: 2nd International Conference on Language and Automata Theory and Applications (LATA 2008)

Cancellative and Malcev presentations for finite Rees index subsemigroups and extensions

Cain, Alan James — Fri, 01 Feb 2008 00:00:00 GMT

Abstract: It is known that, for semigroups, the property of admitting a finite presentation is preserved on passing to subsemigroups and extensions of finite Rees index. The present paper shows that the same holds true for Malcev, cancellative, left-cancellative and right-cancellative presentations. (A Malcev (respectively, cancellative, left-cancellative, right-cancellative) presentation is a presentation of a special type that can be used to define any group-embeddable (respectively, cancellative, left-cancellative, right-cancellative) semigroup.).

The steady-state form of large-amplitude internal solitary waves

King, Stuart Edward — Mon, 10 Jan 2011 00:00:00 GMT

Abstract: A new numerical scheme for obtaining the steady-state form of an internal solitary wave of large amplitude is presented. A stratified inviscid two-dimensional fluid under the Boussinesq approximation flowing between horizontal rigid boundaries is considered. The stratification is stable, and buoyancy is continuously differentiable throughout the domain of the flow. Solutions are obtained by tracing the buoyancy frequency along streamlines from the undisturbed far field. From this the vorticity field can be constructed and the streamfunction may then be obtained by inversion of Laplace's operator. The scheme is presented as an iterative solver, where the inversion of Laplace's operator is performed spectrally. The solutions agree well with previous results for stratification in which the buoyancy frequency is a discontinuous function. The new numerical scheme allows significantly larger amplitude waves to be computed than have been presented before and it is shown that waves with Richardson numbers as low as 0.062 can be computed straightforwardly. The method is also extended to deal in a novel way with closed streamlines when they occur in the domain. The new solutions are tested in independent fully nonlinear time-dependent simulations and are verified to be steady. Waves with regions of recirculation are also discussed.

Impeded inverse energy transfer in the Charney--Hasegawa--Mima model of quasi-geostrophic flows

Tran, Chuong Van — Sat, 25 Mar 2006 00:00:00 GMT

Abstract: The behaviour of turbulent flows within the single-layer quasi-geostrophic (Charney-Hasegawa-Mima) model is shown to be strongly dependent on the Rossby deformation wavenumber lambda (or free-surface elasticity). Herein, we derive a bound oil the inverse energy transfer, specifically on the growth rate dl/dt of the characteristic length scale e representing the energy centroid. It is found that dl/dt <= 2 parallel to q parallel to(infinity)/(l(s)lambda(2)), where parallel to q parallel to(infinity) is the supremum of the potential vorticity and l(s) represents the potential enstrophy centroid of the reservoir, both invariant. This result implies that in the potential-energy-dominated regime (l >= l(s) >> lambda(-1)) the inverse energy transfer is strongly impeded, in the sense that under the usual time scale no significant transfer of energy to larger scales occurs. The physical implication is that the elasticity of the free surface impedes turbulent energy transfer in wavenumber space, effectively rendering large-scale vortices long-lived and inactive. Results from numerical simulations of forced-dissipative turbulence confirm this prediction.

Vanishing enstrophy dissipation in two-dimensional Navier--Stokes turbulence in the inviscid limit

Tran, Chuong Van — Tue, 25 Jul 2006 00:00:00 GMT

Abstract: Batchelor (Phys. Fluids, vol. 12, 1969, p. 233) developed a theory of two-dimensional turbulence based on the assumption that the dissipation of enstrophy (mean-square vorticity) tends to a finite non-zero constant in the limit of infinite Reynolds number Re. Here, by assuming power-law spectra, including the one predicted by Batchelor's theory, we prove that the maximum dissipation of enstrophy is in fact zero in this limit. Specifically, as Re -> infinity, the dissipation approaches zero no slower than (ln Re)(-1/2). The physical reason behind this result is that the decrease of viscosity enhances the production of both palinstrophy (mean-square vorticity gradients) and its dissipation - but in such a way that the net growth of palinstrophy is less rapid than the decrease of viscosity, resulting in vanishing enstrophy dissipation. This result generalizes to a rich class of quasi-geostrophic models as well as to the case of a passive tracer in layerwise-two-dimensional turbulent flows having bounded enstrophy.

Quasi-geostrophic vortices in compressible atmospheres

Scott, Richard Kirkness — Tue, 10 May 2005 00:00:00 GMT

Abstract: This paper considers the effect of an exponential variation in the background density field (as exists in compressible atmospheres) on the structure and dynamics of the quasi-geostrophic system, and compares the results with the corresponding Boussinesq limit in which background density variations are assumed small. The behaviour of the compressible system is understood via a closed-form analytic expression for the Green's function of the inversion operator relating potential vorticity and streamfunction. This expression makes explicit the anisotropy of the Green's function, inherited from the density profile, which has a slow, algebraic decay directly above the source and an exponential decay in all other directions. An immediate consequence for finite-volume vortices is a differential rotation of upper and lower levels that results in counterintuitive behaviour during the nonlinear evolution of ellipsoidal vortices, in which vortex destruction is confined to the lower vortex and wave activity is seen to propagate downwards. This is in contrast to the Boussinesq limit, which exhibits symmetric destruction of the upper and lower vortex, and in contrast to naive expectations based on a consideration of the mass distribution alone, which would lead to greater destruction of the upper vortex. Finally, the presence of a horizontal lower boundary introduces a strong barotropic component that is absent in the unbounded case (the presence of an upper boundary has almost no effect). The lower boundary also alters the differential rotation in the lower vortex with important consequences for the nonlinear evolution: for very small separation between the lower boundary and the vortex, the differential rotation is reversed leading to strong deformations of the middle vortex; for a critical separation, the vortex is stabilized by the reduction of the differential rotation, and remains coherent over remarkably long times.

Subsemigroups of virtually free groups : finite Malcev presentations and testing for freeness

Cain, AJ — Sat, 01 Jul 2006 00:00:00 GMT

Abstract: This paper shows that, given a finite subset X of a finitely generated virtually free group F, the freeness of the subsemigroup of F generated by X can be tested algorithmically. (A group is virtually free if it contains a free subgroup of finite index.) It is then shown that every finitely generated subsemigroup, of F has a finite Malcev presentation (a type of semigroup presentation which can be used to define any semigroup that embeds in a group), and that such a presentation can be effectively found from any finite generating set.

The critical merger distance between two co-rotating quasi-geostrophic vortices

Reinaud, Jean Noel — Mon, 10 Jan 2005 00:00:00 GMT

Abstract: This paper examines the critical merger or strong interaction distance between two equal-potential-vorticity quasi-geostrophic vortices. The interaction between the two vortices depends on five parameters: their volume ratio, their height-to-width aspect ratios and their vertical and horizontal offsets. Due to the size of the parameter space, a direct investigation solving the full quasi-geostrophic equations is impossible. We instead determine the critical merger distance approximately using an asymptotic approach. We associate the merger distance with the margin of stability for a family of equilibrium states having prescribed aspect and volume ratios, and vertical offset. The equilibrium states are obtained using an asymptotic solution method which models vortices by ellipsoids. The margin itself is determined by a linear stability analysis. We focus on the interaction between oblate to moderately prolate vortices, the shapes most commonly found in turbulence. Here, a new unexpected instability is found and discussed for prolate vortices which is manifested by the tilting of vortices toward each other. It implies than tall vortices may merge starting from greater separation distances than previously thought.

The shape of vortices in quasi-geostrophic turbulence

Reinaud, Jean Noel — Fri, 10 Jan 2003 00:00:00 GMT

Abstract: The present work discusses the most commonly occurring shape of the coherent vortical structures in rapidly rotating stably stratified turbulence, under the quasi-geostrophic approximation. In decaying turbulence, these vortices-coherent regions of the materially-invariant potential vorticity-dominate the flow evolution, and indeed the flow evolution is governed by their interactions. An analysis of several exceptionally high-resolution simulations of quasi-geostrophic turbulence is performed. The results indicate that the population of vortices exhibits a mean height-to-width aspect ratio less than unity, in fact close to 0.8. This finding is justified here by a simple model, in which vortices are taken to be ellipsoids of uniform potential vorticity. The model focuses on steady ellipsoids within a uniform background strain flow. This background flow approximates the effects of surrounding vortices in a turbulent flow on a given vortex. It is argued that the vortices which are able to withstand the highest levels of strain are those most likely to be found in the actual turbulent flow. Our calculations confirm that the optimal height-to-width aspect ratio is close to 0.8 for a wide range of background straining flows.

The quasi-geostrophic ellipsoidal vortex model

Dritschel, David Gerard — Sun, 25 Apr 2004 00:00:00 GMT

Abstract: We present a simple approximate model for studying general aspects of vortex interactions in a rotating stably-stratified fluid. The model idealizes vortices by ellipsoidal volumes of uniform potential vorticity, a materially conserved quantity in an inviscid, adiabatic fluid. Each vortex thus possesses 9 degrees of freedom, 3 for the centroid and 6 for the shape and orientation. Here, we develop equations for the time evolution of these quantities for a general system of interacting vortices. An isolated ellipsoidal vortex is well known to remain ellipsoidal in a fluid with constant background rotation and uniform stratification, as considered here. However, the interaction between any two ellipsoids in general induces weak non-ellipsoidal perturbations. We develop a unique projection method, which follows directly from the Hamiltonian structure of the system, that effectively retains just the part of the interaction which preserves ellipsoidal shapes. This method does not use a moment expansion, e.g. local expansions of the flow in a Taylor series. It is in fact more general, and consequently more accurate. Comparisons of the new model with the full equations of motion prove remarkably close.

The merger of vertically offset quasi-geostrophic vortices

Reinaud, Jean Noel — Fri, 25 Oct 2002 00:00:00 GMT

Abstract: We examine the critical merging distance between two equal-volume, equal-potential-vorticity quasi-geostrophic vortices. We focus on how this distance depends on the vertical offset between the two vortices, each having a unit mean height-to-width aspect ratio. The vertical direction is special in the quasi-geostrophic model (used to capture the leading-order dynamical features of stably stratified and rapidly rotating geophysical flows) since vertical advection is absent. Nevertheless vortex merger may still occur by horizontal advection. In this paper, we first investigate the equilibrium states for the two vortices as a function of their vertical and horizontal separation. We examine their basic properties together with their linear stability. These findings are next compared to numerical simulations of the nonlinear evolution of two spheres of potential vorticity. Three different regimes of interaction are identified, depending on the vertical offset. For a small offset, the interaction differs little from the case when the two vortices are horizontally aligned. On the other hand, when the vertical offset is comparable to the mean vortex radius, strong interaction occurs for greater horizontal gaps than in the horizontally aligned case, and therefore at significantly greater full separation distances. This perhaps surprising result is consistent with the linear stability analysis and appears to be a consequence of the anisotropy of the quasi-geostrophic equations. Finally, for large vertical offsets, vortex merger results in the formation of a metastable tilted dumbbell vortex.

The persistence of balance in geophysical flows

Dritschel, David Gerard — Wed, 10 Jan 2007 00:00:00 GMT

Abstract: Rotating stably stratified geophysical flows can exhibit a near 'balanced' evolution controlled by the conservative advection of a single scalar quantity, the potential vorticity (PV). This occurs frequently in the Earth's atmosphere and oceans where motions tend to be weak compared with the background planetary rotation and where stratification greatly inhibits vertical motion. Under these circumstances, both high-frequency acoustic waves and lower-frequency inertia-gravity waves (IGWs) contribute little to the flow evolution compared with the even-lower-frequency advection of PV. Moreover, this 'slow' PV-controlled balanced evolution appears unable to excite these higher-frequency waves in any significant way-i.e. balance persists. The present work pushes the limits of balance by systematically exploring the evolution of a range of highly nonlinear flows in which motions are comparable with the background rotation. These flows do not possess a frequency separation between PV advection and IGWs. Nonetheless, the flows exhibit a remarkable persistence of balance. Even when flows are not initialized to minimize the amount of IGWs initially present, and indeed even when flows are deliberately seeded with significant IGW amplitudes, the flow evolution-over many inertial periods (days)-remains strongly controlled by PV advection. Description: This paper introduces a novel, powerful way to understand the why geophysical flows are largely under the control of a single scalar field, the potential vorticity, a materially conserved dynamical tracer in the absence of viscous and diabatic effects.

Bending and twisting instabilities of columnar elliptical vortices in a rotating strongly stratified fluid

Billant, Paul — Fri, 25 Aug 2006 00:00:00 GMT

Abstract: In this paper, we investigate the three-dimensional stability of the Moore-Saffman elliptical vortex in a rotating stratified fluid. By means of an asymptotic analysis for long vertical wavelength perturbations and small Froude number, we study the effects of Rossby number, external strain, and ellipticity of the vortex on the stability of azimuthal modes m = 1 (corresponding to a bending instability) and m = 2 (corresponding to a twisting instability). In the case of a quasi-geostrophic fluid (small Rossby number), the asymptotic results are in striking agreement with previous numerical stability analyses even for vertical wavelengths of order one. For arbitrary Rossby number, the key finding is that the Rossby number has no effect on the domains of long-wavelength instability of these two modes: the two-dimensional or three-dimensional nature of the instabilities is controlled only by the background strain rate gamma and by the rotation rate Omega of the principal axes of the elliptical vortex relative to the rotating frame of reference. For the m = 1 mode, it is shown that when Omega < -gamma, the vortex is stable to any long-wavelength disturbances, when -gamma < Omega less than or similar to 0, two-dimensional perturbations are most unstable, when 0 less than or similar to Omega < gamma, long-wavelength three-dimensional disturbances are the most unstable, and finally when gamma < Omega, short-wavelength three-dimensional perturbations are the most unstable. Similarly, the m = 2 instability is two-dimensional or three-dimensional depending only on gamma and Omega, independent of the Rossby number. This means that if a long-wavelength three-dimensional instability exists for a given elliptical vortex in a quasi-geostrophic fluid, a similar instability should be observed for any other Rossby number, in particular for infinite Rossby number (strongly stratified fluids). This implies that the planetary rotation plays a minor role in the nature of the instabilities observed in rotating strongly stratified fluids. The present results for the azimuthal mode m = 1 suggest that the vortex-bending instabilities observed previously in quasi-geostrophic fluids (tall-column instability) and in strongly stratified fluids (zigzag instability) are fundamentally related. Description: This is a comprehensive analysis of the linear stability of columnar elliptical vortices subject to two-dimensional strain in a rotating, stratified fluid. It is the culmination of two lines of research, one started by Dritschel involving the tall-column instability, and another started by Billant and Chomaz involving the zigzag instability. Our joint work unifies these instabilities, and shows that they exist over a vast parameter space. This work represents over 7 years of collaborative effort.

Revisiting Batchelor's theory of two-dimensional turbulence

Dritschel, David Gerard — Sun, 25 Nov 2007 00:00:00 GMT

Abstract: Recent mathematical results have shown that a central assumption in the theory of two-dimensional turbulence proposed by Batchelor (Phys. Fluids, vol. 12, 1969, p. 233) is false. That theory, which predicts a X-2/3 k(-1) enstrophy spectrum in the inertial range of freely-decaying turbulence, and which has evidently been successful in describing certain aspects of numerical simulations at high Reynolds numbers Re, assumes that there is a finite, non-zero enstrophy dissipation X in the limit of infinite Re. This, however, is not true for flows having finite vorticity. The enstrophy dissipation in fact vanishes. We revisit Batchelor's theory and propose a simple modification of it to ensure vanishing X in the limit Re -> infinity. Our proposal is supported by high Reynolds number simulations which confirm that X decays like 1/ln Re, and which, following the time of peak enstrophy dissipation, exhibit enstrophy spectra containing an increasing proportion of the total enstrophy (omega(2))/2 in the inertial range as Re increases. Together with the mathematical analysis of vanishing X, these observations motivate a straightforward and, indeed, alarmingly simple modification of Batchelor's theory: just replace Batchelor's enstrophy spectrum X(2/3)k(-1) with (omega(2))k(-1)(In Re)(-1).

A balanced approach to modelling rotating stably stratified geophysical flows

Dritschel, David Gerard — Sun, 10 Aug 2003 00:00:00 GMT

Abstract: We describe a new approach to modelling three-dimensional rotating stratified flows under the Boussinesq approximation. This approach is based on the explicit conservation of potential vorticity, and exploits the underlying leading-order geostrophic and hydrostratic balances inherent in these equations in the limit of small Froude and Rossby numbers. These balances are not imposed, but instead are used to motivate the use of a pair of new variables expressing the departure from geostrophic and hydrostratic balance. These new variables are the ageostrophic horizontal vorticity components, i.e. the vorticity not directly associated with the displacement of isopycnal surfaces. The use of potential vorticity and ageostrophic horizontal vorticity, rather than the usual primitive variables of velocity and density, reveals a deep mathematical structure and appears to have advantages numerically. This change of variables results in a diagnostic equation, of Monge-Amp re type, for one component of a vector potential phi, and two Poisson equations for the other two components. The curl of phi gives the velocity field while the divergence of phi is proportional to the displacement of isopycnal surfaces. This diagnostic equation makes transparent the conditions for both static and inertial stability, and may change form from (spatially) elliptic to (spatially) hyperbolic even when the flow is statically and inertially stable. A numerical method based on these new variables is developed and used to examine the instability of a horizontal elliptical shear zone (modelling a jet streak). The basic-state flow is in exact geostrophic and hydrostratic balance. Given a small perturbation however, the shear zone destabilizes by rolling up into a street of vortices and radiating inertia-gravity waves. Description: This work was the first to show how one can rewrite the equations for a rotating stratified fluid in a way which makes potential vorticity conservation explicit. Potential vorticity is linked closely to balance, a state void of high-frequency gravity waves. The mathematical transformation reveals a deep underlying mathematical structure, including explicit conditions for inertial and static stability as well as a new double Monge-Ampere equation. This work forms the cornerstone of much subsequent research into the fundamental nature of rotating stratified fluids.

Penetrative convection in a superposed porous-medium–fluid layer via internal heating

Carr, Magda — Tue, 01 Jun 2004 00:00:00 GMT