On Arnol'd's second nonlinear stability theorem for two-dimensional quasi-geostrophic flow

Arnol'd's second hydrodynamical stability theorem, proven originally for the two-dimensional Euler equations, can establish nonlinear stability of steady flows that are maxima of a suitably chosen energy-Casimir invariant. The usual derivations of this theorem require an assumption of zero disturbance circulation. In the present work an analogue of Arnol'd's second theorem is developed in the more general case of two-dimensional quasi-geostrophic flow, with the important feature that the disturbances are allowed to have non-zero circulation. New nonlinear stability criteria are derived, and explicit bounds are obtained on both the disturbance energy and potential enstrophy which are expressed in terms of the initial disturbance fields. While Arnol'd's stability method relies on the second variation of the energy-Casimir invariant being sign-definite, the new criteria can be applied to cases where the second variation is sign-indefinite because of the disturbance circulations. A version of Andrews' theorem is also established for this problem.

Small time Chung-type LIL for Lévy processes

We prove Chung-type laws of the iterated logarithm for general Lévy processes at zero. In particular, we provide tools to translate small deviation estimates directly into laws of the iterated logarithm. This reveals laws of the iterated logarithm for Lévy processes at small times in many concrete examples. In some cases, exotic norming functions are derived.

A general method for finding extremal states of Hamiltonian dynamical systems, with applications to perfect fluids

In addition to the Hamiltonian functional itself, non-canonical Hamiltonian dynamical systems generally possess integral invariants known as ‘Casimir functionals’. In the case of the Euler equations for a perfect fluid, the Casimir functionals correspond to the vortex topology, whose invariance derives from the particle-relabelling symmetry of the underlying Lagrangian equations of motion. In a recent paper, Vallis, Carnevale & Young (1989) have presented algorithms for finding steady states of the Euler equations that represent extrema of energy subject to given vortex topology, and are therefore stable. The purpose of this note is to point out a very general method for modifying any Hamiltonian dynamical system into an algorithm that is analogous to those of Vallis etal. in that it will systematically increase or decrease the energy of the system while preserving all of the Casimir invariants. By incorporating momentum into the extremization procedure, the algorithm is able to find steadily-translating as well as steady stable states. The method is applied to a variety of perfect-fluid systems, including Euler flow as well as compressible and incompressible stratified flow.

Flow-induced morphological instability of a mushy layer

A morphological instability of a mushy layer due to a forced flow in the melt is analysed. The instability is caused by flow induced in the mushy layer by Bernoulli suction at the crests of a sinusoidally perturbed mush–melt interface. The flow in the mushy layer advects heat away from crests which promotes solidification. Two linear stability analyses are presented: the fundamental mechanism for instability is elucidated by considering the case of uniform flow of an inviscid melt; a more complete analysis is then presented for the case of a parallel shear flow of a viscous melt. The novel instability mechanism we analyse here is contrasted with that investigated by Gilpin et al. (1980) and is found to be more potent for the case of newly forming sea ice.

Frazil dynamics and precipitation in a water column with depth-dependent supercooling

A morphological instability of a mushy layer due to a forced flow in the melt is analysed. The instability is caused by flow induced in the mushy layer by Bernoulli suction at the crests of a sinusoidally perturbed mush–melt interface. The flow in the mushy layer advects heat away from crests which promotes solidification. Two linear stability analyses are presented: the fundamental mechanism for instability is elucidated by considering the case of uniform flow of an inviscid melt; a more complete analysis is then presented for the case of a parallel shear flow of a viscous melt. The novel instability mechanism we analyse here is contrasted with that investigated by Gilpin et al. (1980) and is found to be more potent for the case of newly forming sea ice.

The influence of ocean flow on newly forming sea ice

The heat and mass balance of the Arctic Ocean is very sensitive to the growth and decay of sea ice and the interaction between the heat and salt fields in the oceanic boundary layer. The hydraulic roughness of sea ice controls the detailed nature of turbulent fluxes in the boundary layer and hence is an important ingredient in model parameterizations. We describe a novel mechanism for the generation of corrugations of the sea ice–ocean interface, present a mathematical analysis elucidating the mechanism, and present numerical calculations for geophysically relevant conditions. The mechanism relies on brine flows developing in the sea ice due to Bernoulli suction by flow of ocean past the interface. For oceanic shears at the ice interface of 0.2 s−1, we expect the corrugations to form with a wavelength dependent upon the permeability structure of the sea ice which is described herein. The mechanism should be particularly important during sea ice formation in wind-maintained coastal polynyas and in leads. This paper applies our earlier analyses of the fundamental instability to field conditions and extends it to take account of the anisotropic and heterogeneous permeability of sea ice.

Kraichnan–Leith–Batchelor similarity theory and two-dimensional inverse cascades

We study the scaling properties and Kraichnan–Leith–Batchelor (KLB) theory of forced inverse cascades in generalized two-dimensional (2D) fluids (α-turbulence models) simulated at resolution 8192x8192. We consider α=1 (surface quasigeostrophic flow), α=2 (2D Euler flow) and α=3. The forcing scale is well resolved, a direct cascade is present and there is no large-scale dissipation. Coherent vortices spanning a range of sizes, most larger than the forcing scale, are present for both α=1 and α=2. The active scalar field for α=3 contains comparatively few and small vortices. The energy spectral slopes in the inverse cascade are steeper than the KLB prediction −(7−α)/3 in all three systems. Since we stop the simulations well before the cascades have reached the domain scale, vortex formation and spectral steepening are not due to condensation effects; nor are they caused by large-scale dissipation, which is absent. One- and two-point p.d.f.s, hyperflatness factors and structure functions indicate that the inverse cascades are intermittent and non-Gaussian over much of the inertial range for α=1 and α=2, while the α=3 inverse cascade is much closer to Gaussian and non-intermittent. For α=3 the steep spectrum is close to that associated with enstrophy equipartition. Continuous wavelet analysis shows approximate KLB scaling ℰ(k)∝k−2 (α=1) and ℰ(k)∝k−5/3 (α=2) in the interstitial regions between the coherent vortices. Our results demonstrate that coherent vortex formation (α=1 and α=2) and non-realizability (α=3) cause 2D inverse cascades to deviate from the KLB predictions, but that the flow between the vortices exhibits KLB scaling and non-intermittent statistics for α=1 and α=2.

Energy consistent discontinuous Galerkin methods for the Navier–Stokes–Korteweg system

We design consistent discontinuous Galerkin finite element schemes for the approximation of the Euler-Korteweg and the Navier-Stokes-Korteweg systems. We show that the scheme for the Euler-Korteweg system is energy and mass conservative and that the scheme for the Navier-Stokes-Korteweg system is mass conservative and monotonically energy dissipative. In this case the dissipation is isolated to viscous effects, that is, there is no numerical dissipation. In this sense the methods are consistent with the energy dissipation of the continuous PDE systems. - See more at: http://www.ams.org/journals/mcom/2014-83-289/S0025-5718-2014-02792-0/home.html#sthash.rwTIhNWi.dpuf

A comparison of duality and energy a posteriori estimates for $\mathrm {L}_{\infty }(0,T;\mathrm {L}_2(\varOmega ))$ in parabolic problems

We use the elliptic reconstruction technique in combination with a duality approach to prove a posteriori error estimates for fully discrete backward Euler scheme for linear parabolic equations. As an application, we combine our result with the residual based estimators from the a posteriori estimation for elliptic problems to derive space-error indicators and thus a fully practical version of the estimators bounding the error in the $ \mathrm {L}_{\infty }(0,T;\mathrm {L}_2(\varOmega ))$ norm. These estimators, which are of optimal order, extend those introduced by Eriksson and Johnson in 1991 by taking into account the error induced by the mesh changes and allowing for a more flexible use of the elliptic estimators. For comparison with previous results we derive also an energy-based a posteriori estimate for the $ \mathrm {L}_{\infty }(0,T;\mathrm {L}_2(\varOmega ))$-error which simplifies a previous one given by Lakkis and Makridakis in 2006. We then compare both estimators (duality vs. energy) in practical situations and draw conclusions.

Memory and burstiness in dynamic networks

A discrete-time random process is described, which can generate bursty sequences of events. A Bernoulli process, where the probability of an event occurring at time t is given by a fixed probability x, is modified to include a memory effect where the event probability is increased proportionally to the number of events that occurred within a given amount of time preceding t. For small values of x the interevent time distribution follows a power law with exponent −2−x. We consider a dynamic network where each node forms, and breaks connections according to this process. The value of x for each node depends on the fitness distribution, \rho(x), from which it is drawn; we find exact solutions for the expectation of the degree distribution for a variety of possible fitness distributions, and for both cases where the memory effect either is, or is not present. This work can potentially lead to methods to uncover hidden fitness distributions from fast changing, temporal network data, such as online social communications and fMRI scans.

Uniform factorial decay estimate for the remainder of rough Taylor expansion

We establish an uniform factorial decay estimate for the Taylor approximation of solutions to controlled differential equations. Its proof requires a factorial decay estimate for controlled paths which is interesting in its own right.

The problem of two fixed centers: bifurcation diagram for positive energies

We give a comprehensive analysis of the Euler-Jacobi problem of motion in the field of two fixed centers with arbitrary relative strength and for positive values of the energy. These systems represent nontrivial examples of integrable dynamics and are analysed from the point of view of the energy-momentum mapping from the phase space to the space of the integration constants. In this setting, we describe the structure of the scattering trajectories in phase space and derive an explicit description of the bifurcation diagram, i.e., the set of critical value of the energy-momentum map.

On the structure of infinity-harmonic maps

Let H ∈ C 2(ℝ N×n ), H ≥ 0. The PDE system arises as the Euler-Lagrange PDE of vectorial variational problems for the functional E ∞(u, Ω) = ‖H(Du)‖ L ∞(Ω) defined on maps u: Ω ⊆ ℝ n → ℝ N . (1) first appeared in the author's recent work. The scalar case though has a long history initiated by Aronsson. Herein we study the solutions of (1) with emphasis on the case of n = 2 ≤ N with H the Euclidean norm on ℝ N×n , which we call the “∞-Laplacian”. By establishing a rigidity theorem for rank-one maps of independent interest, we analyse a phenomenon of separation of the solutions to phases with qualitatively different behaviour. As a corollary, we extend to N ≥ 2 the Aronsson-Evans-Yu theorem regarding non existence of zeros of |Du| and prove a maximum principle. We further characterise all H for which (1) is elliptic and also study the initial value problem for the ODE system arising for n = 1 but with H(·, u, u′) depending on all the arguments.

Optimal infinity-quasiconformal immersions

For a Hamiltonian K ∈ C2(RN × n) and a map u:Ω ⊆ Rn − → RN, we consider the supremal functional (1) The “Euler−Lagrange” PDE associated to (1)is the quasilinear system (2) Here KP is the derivative and [ KP ] ⊥ is the projection on its nullspace. (1)and (2)are the fundamental objects of vector-valued Calculus of Variations in L∞ and first arose in recent work of the author [N. Katzourakis, J. Differ. Eqs. 253 (2012) 2123–2139; Commun. Partial Differ. Eqs. 39 (2014) 2091–2124]. Herein we apply our results to Geometric Analysis by choosing as K the dilation function which measures the deviation of u from being conformal. Our main result is that appropriately defined minimisers of (1)solve (2). Hence, PDE methods can be used to study optimised quasiconformal maps. Nonconvexity of K and appearance of interfaces where [ KP ] ⊥ is discontinuous cause extra difficulties. When n = N, this approach has previously been followed by Capogna−Raich ? and relates to Teichmüller’s theory. In particular, we disprove a conjecture appearing therein.

A METHOD FOR THE STUDY OF ACCRETION DISK EMISSION IN CATACLYSMIC VARIABLES. I. THE MODEL

We have developed a spectrum synthesis method for modeling the ultraviolet (UV) emission from the accretion disk from cataclysmic variables (CVs). The disk is separated into concentric rings, with an internal structure from the Wade & Hubeny disk-atmosphere models. For each ring, a wind atmosphere is calculated in the comoving frame with a vertical velocity structure obtained from a solution of the Euler equation. Using simple assumptions, regarding rotation and the wind streamlines, these one-dimensional models are combined into a single 2.5-dimensional model for which we compute synthetic spectra. We find that the resulting line and continuum behavior as a function of the orbital inclination is consistent with the observations, and verify that the accretion rate affects the wind temperature, leading to corresponding trends in the intensity of UV lines. In general, we also find that the primary mass has a strong effect on the P Cygni absorption profiles, the synthetic emission line profiles are strongly sensitive to the wind temperature structure, and an increase in the mass-loss rate enhances the resonance line intensities. Synthetic spectra were compared with UV data for two high orbital inclination nova-like CVs-RW Tri and V347 Pup. We needed to include disk regions with arbitrary enhanced mass loss to reproduce reasonably well widths and line profiles. This fact and a lack of flux in some high ionization lines may be the signature of the presence of density-enhanced regions in the wind, or alternatively, may result from inadequacies in some of our simplifying assumptions.