Colm Connaughton

Condensation of Classical Nonlinear Waves

We study the formation of a large-scale coherent structure (a condensate) in classical wave equations by considering the defocusing nonlinear Schrödinger equation as a representative model. We formulate a thermodynamic description of the classical condensation process by using a wave turbulence theory with ultraviolet cutoff. In three dimensions the equilibrium state undergoes a phase transition for sufficiently low energy density, while no transition occurs in two dimensions, in complete analogy with standard Bose-Einstein condensation in quantum systems. On the basis of a modified wave turbulence theory, we show that the nonlinear interaction makes the transition to condensation subcritical. The theory is in quantitative agreement with the numerical integration of the nonlinear Schrödinger equation.

Modulational instability of Rossby and drift waves and generation of zonal jets

We study the modulational instability of geophysical Rossby and plasma drill waves within the Charney-Hasegawa-Mima (CH M) model both theoretically, using truncated (four-mode and three-mode) models, and numerically, using direct simulations of CHM equation in the Fourier space. We review the linear theory of Gill (Geophys. Fluid Dyn., vol. 6, 1974, p. 29) and extend it to show that for strong primary waves the most unstable modes are perpendicular to the primary wave, which correspond to generation of a zonal flow if the primary wave is purely meridional. For weak waves, the maximum growth occurs for off-zonal inclined modulations that are close to being in three-wave resonance with the primary wave. Our numerical simulations confirm the theoretical predictions of the linear theory as well as the nonlinear jet pinching predicted by Manin & Nazarenko (Pit vs. Fluids, vol. 6, 1994, p. 1158). We find that, for strong primary waves, these narrow zonal jets further roll up into Karman-like vortex streets, and at this moment the truncated models fail. For weak primary waves, the growth of the unstable mode reverses and the system oscillates between a dominant jet and a dominate primary wave, so that the truncated description holds for longer. The two-dimensional vortex streets appear to be more stable than purely one-dimensional zonal jets, and their zonal-averaged speed can reach amplitudes much stronger than is allowed by the Rayleigh-Kuo instability criterion for the one-dimensional case. In the long term, the system transitions to turbulence helped by the vortex-pairing instability (for strong waves) and the resonant wave wave interactions (for weak waves).

Dimensional analysis and weak turbulence

Abstract In the study of weakly turbulent wave systems possessing incomplete self-similarity, it is possible to use dimensional arguments to derive the scaling exponents of the Kolmogorov–Zakharov spectra, provided the order of the resonant wave interactions responsible for nonlinear energy transfer is known. Furthermore, one can easily derive conditions for the breakdown of the weak turbulence approximation. It is found that for incompletely self-similar systems dominated by three wave interactions, the weak turbulence approximation usually cannot break down at small scales. It follows that such systems cannot exhibit small scale intermittency. For systems dominated by four wave interactions, the incomplete self-similarity property implies that the scaling of the interaction coefficient depends only on the physical dimension of the system. These results are used to build a complete picture of the scaling properties of the surface wave problem where both gravity and surface tension play a role. We argue that, for large values of the energy flux, there should be two weakly turbulent scaling regions matched together via a region of strongly nonlinear turbulence.

Non-stationary spectra of local wave turbulence

Abstract The evolution of the Kolmogorov–Zakharov (K–Z) spectrum of weak turbulence is studied in the limit of strongly local interactions where the usual kinetic equation, describing the time evolution of the spectral wave-action density, can be approximated by a PDE. If the wave action is initially compactly supported in frequency space, it is then redistributed by resonant interactions producing the usual direct and inverse cascades, leading to the formation of the K–Z spectra. The emphasis here is on the direct cascade. The evolution proceeds by the formation of a self-similar front which propagates to the right leaving a quasi-stationary state in its wake. This front is sharp in the sense that the solution remains compactly supported until it reaches infinity. If the energy spectrum has infinite capacity, the front takes infinite time to reach infinite frequency and leaves the K–Z spectrum in its wake. On the other hand, if the energy spectrum has finite capacity, the front reaches infinity within a finite time, t ∗ , and the wake is steeper than the K–Z spectrum. For this case, the K–Z spectrum is set up from the right after the front reaches infinity. The slope of the solution in the wake can be related to the speed of propagation of the front. It is shown that the anomalous slope in the finite capacity case corresponds to the unique front speed which ensures that the front tip contains a finite amount of energy as the connection to infinity is made. We also introduce, for the first time, the notion of entropy production in wave turbulence and show how it evolves as the system approaches the stationary K–Z spectrum.

Feedback of zonal flows on wave turbulence driven by small-scale instability in the Charney-Hasegawa-Mima model

We demonstrate theoretically and numerically the zonal-flow/drift-wave feedback mechanism in an idealised 2-dimensional model of plasma turbulence driven by a small-scale instability. Zonal flows are generated by a secondary modulational instability of the modes which are directly driven by the primary instability. The zonal flows then suppress the small scales thereby arresting the energy injection into the system, a process which can be described using nonlocal (in scale) wave turbulence theory. Finally, the arrest of the energy input results in saturation of the zonal flows at a level which can be estimated from the theory and the system reaches stationarity without large-scale dissipation.

Numerical solutions of the isotropic 3-wave kinetic equation

We show that the isotropic 3-wave kinetic equation is equivalent to the mean-field rate equations for an aggregation-fragmentation problem with an unusual fragmentation mechanism. This analogy is used to write the theory of 3-wave turbulence almost entirely in terms of a single scaling parameter. A new numerical method for solving the kinetic equation over a large range of frequencies is developed by extending Lees method for solving aggregation equations. The new algorithm is validated against some analytical calculations of the Kolmogorov-Zakharov (KZ) constant for some families of model interaction coefficients. The algorithm is then applied to study some wave turbulence problems in which the finiteness of the dissipation scale is an essential feature. Firstly, it is shown that for finite capacity cascades, the dissipation of energy becomes independent of the cut-off frequency as this cut-off is taken to infinity. This is an explicit indication of the presence of a dissipative anomaly. Secondly, a preliminary numerical study is presented of the so-called bottleneck effect in a wave turbulence context. It is found that the structure of the bottleneck, depends non-trivially on the interaction coefficient. Finally, some results are presented on the complementary phenomenon of thermalisation in closed wave systems which demonstrates explicitly for the first time the existence of so-called mixed solutions of the kinetic equation which exhibit aspects of both KZ and equilibrium equipartition spectra.

Breakdown of Kolmogorov scaling in models of cluster aggregation

We describe a model of cluster aggregation with a source which provides a rare example of an analytically tractable turbulent system. The steady state is characterized by a constant mass flux from small masses to large. Thus it can be studied using a phenomenological theory, inspired by Kolmogorovs 1941 theory, which assumes constant flux and self-similarity. We prove that such self-similarity is violated in dimensions less than or equal to two. We then use dynamical renormalization group techniques to show that the scaling of multipoint correlation functions implies nontrivial multifractality. The analytical results are supported by Monte Carlo simulations.

Rossby and drift wave turbulence and zonal flows : The Charney–Hasegawa–Mima model and its extensions

A detailed study of the Charney-Hasegawa-Mima model and its extensions is presented. These simple nonlinear partial differential equations suggested for both Rossby waves in the atmosphere and also drift waves in a magnetically-confined plasma exhibit some remarkable and nontrivial properties, which in their qualitative form survive in more realistic and complicated models, and as such form a conceptual basis for understanding the turbulence and zonal flow dynamics in real plasma and geophysical systems. Two idealised scenarios of generation of zonal flows by small-scale turbulence are explored: a modulational instability and turbulent cascades. A detailed study of the generation of zonal flows by the modulational instability reveals that the dynamics of this zonal flow generation mechanism differ widely depending on the initial degree of nonlinearity. A numerical proof is provided for the extra invariant in Rossby and drift wave turbulence -zonostrophy and the invariant cascades are shown to be characterised by the zonostrophy pushing the energy to the zonal scales. A small scale instability forcing applied to the model demonstrates the well-known drift wave - zonal flow feedback loop in which the turbulence which initially leads to the zonal flow creation, is completely suppressed and the zonal flows saturate. The turbulence spectrum is shown to diffuse in a manner which has been mathematically predicted. The insights gained from this simple model could provide a basis for equivalent studies in more sophisticated plasma and geophysical fluid dynamics models in an effort to fully understand the zonal flow generation, the turbulent transport suppression and the zonal flow saturation processes in both the plasma and geophysical contexts as well as other wave and turbulence systems where order evolves from chaos.

Driven Brownian coagulation of polymers

We present an analysis of the mean-field kinetics of Brownian coagulation of droplets and polymers driven by input of monomers which aims to characterize the long time behavior of the cluster size distribution as a function of the inverse fractal dimension, a, of the aggregates. We find that two types of long time behavior are possible. For 0≤a<1/2 the size distribution reaches a stationary state with a power law distribution of cluster sizes having exponent 3/2. The amplitude of this stationary state is determined exactly as a function of a. For 1/2<a≤1, the cluster size distribution never reaches a stationary state. Instead a bimodal distribution is formed in which a narrow population of small clusters near the monomer scale is separated by a gap (where the cluster size distribution is effectively zero) from a population of large clusters which continue to grow for all time by absorbing small clusters. The marginal case, a=1/2, is difficult to analyze definitively, but we argue that the cluster size distribution becomes stationary and there is a logarithmic correction to the algebraic tail.

Combining Gaussian processes, mutual information and a genetic algorithm for multi-target optimization of expensive-to-evaluate functions

A novel approach to multi-target optimization of expensive-to-evaluate functions is explored that is based on a combined application of Gaussian processes, mutual information and a genetic algorithm. The aim of the approach is to find an approximation to the optimal solution (or the Pareto optimal solutions) within a small budget. The approach is shown to compare favourably with a surrogate based online evolutionary algorithm on two synthetic problems.