Beyond the Spectral Theorem: Spectrally Decomposing Arbitrary Functions of Nondiagonalizable Operators

Nonlinearities in finite dimensions can be linearized by projecting them into infinite dimensions. Unfortunately, often the linear operator techniques that one would then use simply fail since the operators cannot be diagonalized. This curse is well known. It also occurs for finite-dimensional linear operators. We circumvent it by developing a meromorphic functional calculus that can decompose arbitrary functions of nondiagonalizable linear operators in terms of their eigenvalues and projection operators. It extends the spectral theorem of normal operators to a much wider class, including circumstances in which poles and zeros of the function coincide with the operator spectrum. By allowing the direct manipulation of individual eigenspaces of nonnormal and nondiagonalizable operators, the new theory avoids spurious divergences. As such, it yields novel insights and closed-form expressions across several areas of physics in which nondiagonalizable dynamics are relevant, including memoryful stochastic processes, open non unitary quantum systems, and far-from-equilibrium thermodynamics. The technical contributions include the first full treatment of arbitrary powers of an operator. In particular, we show that the Drazin inverse, previously only defined axiomatically, can be derived as the negative-one power of singular operators within the meromorphic functional calculus and we give a general method to construct it. We provide new formulae for constructing projection operators and delineate the relations between projection operators, eigenvectors, and generalized eigenvectors. By way of illustrating its application, we explore several, rather distinct examples.Comment: 29 pages, 4 figures, expanded historical citations; http://csc.ucdavis.edu/~cmg/compmech/pubs/bst.ht

On Some Classes of mKdV Periodic Solutions

We obtain exact periodic solutions of the positive and negative modified Kortweg-de Vries (mKdV) equations. We examine the dynamical stability of these solitary wave lattices through direct numerical simulations. While the positive mKdV breather lattice solutions are found to be unstable, the two-soliton lattice solution of the same equation is found to be stable. Similarly, a negative mKdV lattice solution is found to be stable. We also touch upon the implications of these results for the KdV equation.Comment: 8 pages, 3 figures, to appear in J. Phys.

Linear and nonlinear dynamics in stratified shear flows

Stably stratified shear flows, in which a less dense layer of fluid lies above and moves counter to a more dense layer below, are ubiquitous in geophysical fluid dynamics. These are often found to be unstable if the non-dimensional Richardson number Ri, quantifying the strength of stratification to shear, is sufficiently low. This is of particular importance in oceanography, where shear instabilities are conjectured to be important in the generation of turbulence in the deep ocean, an area of huge uncertainty in contemporary climate models. The Miles-Howard theorem tells us that for a steady, inviscid, parallel shear flow, if the local Richardson number is everywhere greater than one quarter, the flow is stable to infinitesimal perturbations. Though an important result, the strong restrictions in the applicability of this theorem mean care must be used when applying the criterion of Ri > 1/4 for stability. This thesis explores some of these limitations, beginning with an overview in chapter 1. Chapter 2 explores the infinitesimal restriction of the Miles-Howard theorem, by asking whether finite-amplitude perturbations could lead to significant nonlinear behaviour, in a so-called subcritical instability. It is found that while the classical Kelvin-Helmholtz instability does indeed exhibit subcriticality, nonlinear steady states are found only just above Ri = 1/4. Chapter 3 investigates in detail a hitherto unknown linear instability, which was discovered in chapter 2. Behaving similarly to the classic Holmboe instability, it exists for Ri > 1/4 when viscosity is introduced, and reveals new insights into the possible physical interpretations of stratified shear instability. Chapter 4 revisits the results of chapter 2 but considers two cases of the Prandtl number Pr, the ratio of diffusivity of the momentum to density. When Pr = 0.7, as is approximately the case for air, a simple supercritical instability is found. However, for Pr = 7, corresponding approximately to water, strong subcritical behaviour is observed, and it is demonstrated that finite-amplitude perturbations can trigger Kelvin-Helmholtz-like behaviour well above Ri = 1/4. Chapter 5 considers the time-varying, non-parallel flow of an oblique internal gravity wave incident on a shear layer. Using direct-adjoint looping, it is shown that the disturbances which maximise energy after a certain time, so-called linear optimal perturbations, can be convective-like rolls in the spanwise direction, rather than a shear instability, calling into question the relevance of the classical shear instabilities in oceanography. Chapter 6 concludes the thesis with a discussion of the implications of the results

On the Exponential Stability of Stochastic Perturbed Singular Systems in Mean Square

The approach of Lyapunov functions is one of the most efficient ones for the investigation of the stability of stochastic systems, in particular, of singular stochastic systems. The main objective of the paper is the analysis of the stability of stochastic perturbed singular systems by using Lyapunov techniques under the assumption that the initial conditions are consistent. The uniform exponential stability in mean square and the practical uniform exponential stability in mean square of solutions of stochastic perturbed singular systems based on Lyapunov techniques are investigated. Moreover, we study the problem of stability and stabilization of some classes of stochastic singular systems. Finally, an illustrative example is given to illustrate the effectiveness of the proposed results

Linear dynamics of weakly viscous accretion disks: A disk analog of Tollmien-Schlichting waves

This paper discusses new perspectives and approaches to the problem of disk dynamics where, in this study, we focus on the effects of viscous instabilities influenced by boundary effects. The Boussinesq approximation of the viscous large shearing box equations is analyzed in which the azimuthal length scale of the disturbance is much larger than the radial and vertical scales. We examine the stability of a non-axisymmetric potential vorticity mode, i.e. a PV-anomaly. in a configuration in which buoyant convection and the strato-rotational instability do not to operate. We consider a series of boundary conditions which show the PV-anomaly to be unstable both on a finite and semi-infinite radial domains. We find these conditions leading to an instability which is the disk analog of Tollmien-Schlichting waves. When the viscosity is weak, evidence of the instability is most pronounced by the emergence of a vortex sheet at the critical layer located away from the boundary where the instability is generated. For some boundary conditions a necessary criterion for the onset of instability for vertical wavelengths that are a sizable fraction of the layer's thickness and when the viscosity is small is that the appropriate Froude number of the flow be greater than one. This instability persists if more realistic boundary conditions are applied, although the criterion on the Froude number is more complicated. The unstable waves studied here share qualitative features to the instability seen in rotating Blasius boundary layers. The implications of these results are discussed. An overall new strategy for exploring and interpreting disk instability mechanisms is also suggested.Comment: Accepted for publication in Astronomy and Astrophysics. 18 pages. This version 3 with corrected style fil