Source author record

Colm Connaughton

Colm Connaughton appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

28works
15topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

28 published item(s)

preprint2022arXiv

Forecasting new diseases in low-data settings using transfer learning

Recent infectious disease outbreaks, such as the COVID-19 pandemic and the Zika epidemic in Brazil, have demonstrated both the importance and difficulty of accurately forecasting novel infectious diseases. When new diseases first emerge, we have little knowledge of the transmission process, the level and duration of immunity to reinfection, or other parameters required to build realistic epidemiological models. Time series forecasts and machine learning, while less reliant on assumptions about the disease, require large amounts of data that are also not available in early stages of an outbreak. In this study, we examine how knowledge of related diseases can help make predictions of new diseases in data-scarce environments using transfer learning. We implement both an empirical and a theoretical approach. Using empirical data from Brazil, we compare how well different machine learning models transfer knowledge between two different disease pairs: (i) dengue and Zika, and (ii) influenza and COVID-19. In the theoretical analysis, we generate data using different transmission and recovery rates with an SIR compartmental model, and then compare the effectiveness of different transfer learning methods. We find that transfer learning offers the potential to improve predictions, even beyond a model based on data from the target disease, though the appropriate source disease must be chosen carefully. While imperfect, these models offer an additional input for decision makers during pandemic response.

preprint2021arXiv

Dynamic and interpretable hazard-based models of traffic incident durations

Understanding and predicting the duration or "return-to-normal" time of traffic incidents is important for system-level management and optimisation of road transportation networks. Increasing real-time availability of multiple data sources characterising the state of urban traffic networks, together with advances in machine learning offer the opportunity for new and improved approaches to this problem that go beyond static statistical analyses of incident duration. In this paper we consider two such improvements: dynamic update of incident duration predictions as new information about incidents becomes available and automated interpretation of the factors responsible for these predictions. For our use case, we take one year of incident data and traffic state time-series from the M25 motorway in London. We use it to train models that predict the probability distribution of incident durations, utilising both time-invariant and time-varying features of the data. The latter allow predictions to be updated as an incident progresses, and more information becomes available. For dynamic predictions, time-series features are fed into the Match-Net algorithm, a temporal convolutional hitting-time network, recently developed for dynamical survival analysis in clinical applications. The predictions are benchmarked against static regression models for survival analysis and against an established dynamic technique known as landmarking and found to perform favourably by several standard comparison measures. To provide interpretability, we utilise the concept of Shapley values from the domain of interpretable artificial intelligence to rank the features most relevant to the model predictions at different time horizons. For example, the time of day is always a significantly influential time-invariant feature, whereas the time-series features strongly influence predictions at 5 and 60-minute horizons.

preprint2020arXiv

A non-parametric Hawkes process model of primary and secondary accidents on a UK smart motorway

A self-exciting spatio-temporal point process is fitted to incident data from the UK National Traffic Information Service to model the rates of primary and secondary accidents on the M25 motorway in a 12-month period during 2017-18. This process uses a background component to represent primary accidents, and a self-exciting component to represent secondary accidents. The background consists of periodic daily and weekly components, a spatial component and a long-term trend. The self-exciting components are decaying, unidirectional functions of space and time. These components are determined via kernel smoothing and likelihood estimation. Temporally, the background is stable across seasons with a daily double peak structure reflecting commuting patterns. Spatially, there are two peaks in intensity, one of which becomes more pronounced during the study period. Self-excitation accounts for 6-7% of the data with associated time and length scales around 100 minutes and 1 kilometre respectively. In-sample and out-of-sample validation are performed to assess the model fit. When we restrict the data to incidents that resulted in large speed drops on the network, the results remain coherent.

preprint2020arXiv

Discovering causal factors of drought in Ethiopia

Drought is a costly natural hazard, many aspects of which remain poorly understood. It has many contributory factors, driving its outset, duration, and severity, including land surface, anthropogenic activities, and, most importantly, meteorological anomalies. Prediction plays a crucial role in drought preparedness and risk mitigation. However, this is a challenging task at socio-economically critical lead times (1-2 years), because meteorological anomalies operate at a wide range of temporal and spatial scales. Among them, past studies have shown a correlation between the Sea Surface Temperature (SST) anomaly and the amount of precipitation in various locations in Africa. In its Eastern part, the cooling phase of El Nino-Southern Oscillation (ENSO) and SST anomaly in the Indian ocean are correlated with the lack of rainfall. Given the intrinsic shortcomings of correlation coefficients, we investigate the association among SST modes of variability and the monthly fraction of grid points in Ethiopia, which are in drought conditions in terms of causality. Using the empirical extreme quantiles of precipitation distribution as a proxy for drought, We show that the level of SST second mode of variability in the prior year influences the occurrence of drought in Ethiopia. The causal link between these two variables has a negative coefficient that verifies the conclusion of past studies that rainfall deficiency in the Horn of Africa is associated with ENSO's cooling phase.

preprint2020arXiv

Disease and information spreading at different speeds in multiplex networks

Nowadays, one of the challenges we face when carrying out modeling of epidemic spreading is to develop methods to control disease transmission. In this article we study how the spreading of knowledge of a disease affects the propagation of that disease in a population of interacting individuals. For that, we analyze the interaction between two different processes on multiplex networks: the propagation of an epidemic using the susceptible-infected-susceptible dynamics and the dissemination of information about the disease --and its prevention methods-- using the unaware-aware-unaware dynamics, so that informed individuals are less likely to be infected. Unlike previous related models where disease and information spread at the same time scale, we introduce here a parameter that controls the relative speed between the propagation of the two processes. We study the behavior of this model using a mean-field approach that gives results in good agreement with Monte Carlo simulations on homogeneous complex networks. We find that increasing the rate of information dissemination reduces the disease prevalence, as one may expect. However, increasing the speed of the information process as compared to that of the epidemic process has the counter intuitive effect of increasing the disease prevalence. This result opens an interesting discussion about the effects of information spreading on disease propagation.

preprint2019arXiv

The role of zero-clusters in exchange-driven growth with and without input

The exchange-driven growth model describes the mean field kinetics of a population of composite particles (clusters) subject to pairwise exchange interactions. Exchange in this context means that upon interaction of two clusters, one loses a constituent unit (monomer) and the other gains this unit. Two variants of the exchange-driven growth model appear in applications. They differ in whether clusters of zero size are considered active or passive. In the active case, clusters of size zero can acquire a monomer from clusters of positive size. In the passive case they cannot, meaning that clusters reaching size zero are effectively removed from the system. The large time behaviour is very different for the two variants of the model. We first consider an isolated system. In the passive case, the cluster size distribution tends towards a self-similar evolution and the typical cluster size grows as a power of time. In the active case, we identify a broad class of kernels for which the the cluster size distribution tends to a non-trivial time-independent equilibrium in which the typical cluster size is finite. We next consider a non-isolated system in which monomers are input at a constant rate. In the passive case, the cluster size distribution again attains a self-similar profile in which the typical cluster size grows as a power of time. In the active case, a surprising new behavior is found: the cluster size distribution asymptotes to the same equilibrium profile found in the isolated case but with an amplitude that grows linearly in time.

preprint2016arXiv

Non-equilibrium phase diagram for a model with coalescence, evaporation and deposition

We study a d-dimensional lattice model of diffusing coalescing massive particles, with two parameters controlling deposition and evaporation of monomers. The unique stationary distribution for the system exhibits a phase transition in all dimensions d greater or equal to one between a growing phase, in which the expected mass is infinite at each site, and an exponential phase in which the expected mass is finite. We establish rigorous upper and lower bounds on the critical curve describing the phase transition for this system, and some asymptotics for large or small deposition rates.

preprint2015arXiv

Importance Sampling Variance Reduction for the Fokker-Planck Rarefied Gas Particle Method

Models and methods that are able to accurately and efficiently predict the flows of low-speed rarefied gases are in high demand, due to the increasing ability to manufacture devices at micro and nano scales. One such model and method is a Fokker-Planck approximation to the Boltzmann equation, which can be solved numerically by a stochastic particle method. The stochastic nature of this method leads to noisy estimates of the thermodynamic quantities one wishes to sample when the signal is small in comparison to the thermal velocity of the gas. Recently, Gorji et al have proposed a method which is able to greatly reduce the variance of the estimators, by creating a correlated stochastic process which acts as a control variate for the noisy estimates. However, there are potential difficulties involved when the geometry of the problem is complex, as the method requires the density to be solved for independently. Importance sampling is a variance reduction technique that has already been shown to successfully reduce the noise in direct simulation Monte Carlo calculations. In this paper we propose an importance sampling method for the Fokker-Planck stochastic particle scheme. The method requires minimal change to the original algorithm, and dramatically reduces the variance of the estimates. We test the importance sampling scheme on a homogeneous relaxation, planar Couette flow and a lid-driven-cavity flow, and find that our method is able to greatly reduce the noise of estimated quantities. Significantly, we find that as the characteristic speed of the flow decreases, the variance of the noisy estimators becomes independent of the characteristic speed.

preprint2014arXiv

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.

preprint2013arXiv

A nonlinear least squares method for the inverse droplet coagulation problem

If the rates, $K(x,y)$, at which particles of size $x$ coalesce with particles of size $y$ is known, then the mean-field evolution of the particle-size distribution of an ensemble of irreversibly coalescing particles is described by the Smoluchowski equation. We study the corresponding inverse problem which aims to determine the coalescence rates, $K(x,y)$ from measurements of the particle size distribution. We assume that $K(x,y)$ is a homogeneous function of its arguments, a case which occurs commonly in practice. The problem of determining, $K(x,y)$, a function to two variables, then reduces to a simpler problem of determining a function of a single variable plus two exponents, $μ$ and $ν$, which characterise the scaling properties of $K(x,y)$. The price of this simplification is that the resulting least squares problem is nonlinear in the exponents $μ$ and $ν$. We demonstrate the effectiveness of the method on a selection of coalescence problems arising in polymer physics, cloud science and astrophysics. The applications include examples in which the particle size distribution is stationary owing to the presence of sources and sinks of particles and examples in which the particle size distribution is undergoing self-similar relaxation in time.

preprint2013arXiv

Modulational Instability in Basic Plasma and Geophysical Models

This is a review of the theory of the modulational instability in idealised fluid models of strongly magnetised plasmas and reduced models of geophysical fluid dynamics, particularly the role it plays in the formation of zonal flows. The discussion focusses on the Charney-Hasegawa-Mima and Hasegawa-Wakatani models. Particular attention is paid to the wave turbulence - zonal flow feedback loop whereby large scale zonal flows which are initially generated by modulational instability of small-scale drift/Rossby waves tend to subsequently suppress these small scale waves by their shearing action. This negative feedback can result in a dynamic equilibrium in which large scale zonal flows grow by drawing energy from small scale turbulence but suppress the small scale turbulence in the process until a balance is reached. In this regime, the level of small scale turbulence is greatly reduced compared to the level one would observe in the absence of the zonal flows.

preprint2013arXiv

Wave turbulence in the two-layer ocean model

This paper looks at the two-layer ocean model from a wave turbulence perspective. A symmetric form of the two-layer kinetic equation for Rossby waves is derived using canonical variables, allowing the turbulent cascade of energy between the barotropic and baroclinic modes to be studied. It turns out that energy is transferred via local triad interactions from the large-scale baroclinic modes to the baroclinic and barotropic modes at the Rossby deformation scale. From there it is then transferred to the large-scale barotropic modes via a nonlocal inverse transfer. Using scale separation a sys- tem of coupled equations were obtained for the small-scale baroclinic component and the large-scale barotropic component. Since the total energy of the small-scale component is not conserved, but the total barotropic plus baroclinic energy is conserved, the baroclinic energy loss at small scales will be compensated by the growth of the barotropic energy at large scales. It is found that this transfer is mostly anisotropic and mostly to the zonal component.

preprint2012arXiv

Collective oscillations in driven coagulation

We present a novel form of collective oscillatory behavior in the kinetics of irreversible coagulation with a constant input of monomers and removal of large clusters. For a broad class of collision rates, this system reaches a non-equilibrium stationary state at large times and the cluster size distribution tends to a universal form characterised by a constant flux of mass through the space of cluster sizes. Universality, in this context, means that the stationary state becomes independent of the cut-off as the cut-off grows. This universality is lost, however, if the aggregation rate between large and small clusters increases sufficiently steeply as a function of cluster sizes. We identify a transition to a regime in which the stationary state vanishes as the cut-off grows. This non-universal stationary state becomes unstable, however, as the cut-off is increased and undergoes a Hopf bifurcation. After this bifurcation, the stationary kinetics are replaced by persistent and periodic collective oscillations. These oscillations carry pulses of mass through the space of cluster sizes. As a result, the average mass flux remains constant. Furthermore, universality is partially restored in the sense that the scaling of the period and amplitude of oscillation is inherited from the dynamical scaling exponents of the universal regime. The implications of this new type of long-time asymptotic behaviour for other driven non-equilibrium systems are discussed.

preprint2012arXiv

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\leq 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 \leq 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.

preprint2012arXiv

Externally forced triads of resonantly interacting waves: boundedness and integrability properties

We revisit the problem of a triad of resonantly interacting nonlinear waves driven by an external force applied to the unstable mode of the triad. The equations are Hamiltonian, and can be reduced to a dynamical system for 5 real variables with 2 conservation laws. If the Hamiltonian, $H$, is zero we reduce this dynamical system to the motion of a particle in a one-dimensional time-independent potential and prove that the system is integrable. Explicit solutions are obtained for some particular initial conditions. When explicit solution is not possible we present a novel numerical/analytical method for approximating the dynamics. Furthermore we show analytically that when H=0 the motion is generically bounded. That is to say the waves in the forced triad are bounded in amplitude for all times for any initial condition with the single exception of one special choice of initial condition for which the forcing is in phase with the nonlinear oscillation of the triad. This means that the energy in the forced triad generically remains finite for all time despite the fact that there is no dissipation in the system. We provide a detailed characterisation of the dependence of the period and maximum energy of the system on the conserved quantities and forcing intensity. When $H \neq 0$ we reduce the problem to the motion of a particle in a one-dimensional time-periodic potential. Poincaré sections of this system provide strong evidence that the motion remains bounded when $H \neq 0$ and is typically quasi-periodic although periodic orbits can certainly be found. Throughout our analyses, the phases of the modes in the triad play a crucial role in understanding the dynamics.

preprint2012arXiv

Percolation transition in the kinematics of nonlinear resonance broadening in Charney-Hasegawa-Mima model of Rossby wave turbulence

We study the kinematics of nonlinear resonance broadening of interacting Rossby waves as modelled by the Charney-Hasegawa-Mima equation on a biperiodic domain. We focus on the set of wave modes which can interact quasi-resonantly at a particular level of resonance broadening and aim to characterise how the structure of this set changes as the level of resonance broadening is varied. The commonly held view that resonance broadening can be thought of as a thickening of the resonant manifold is misleading. We show that in fact the set of modes corresponding to a single quasi-resonant triad has a nontrivial structure and that its area in fact diverges for a finite degree of broadening. We also study the connectivity of the network of modes which is generated when quasi-resonant triads share common modes. This network has been argued to form the backbone for energy transfer in Rossby wave turbulence. We show that this network undergoes a percolation transition when the level of resonance broadening exceeds a critical value. Below this critical value, the largest connected component of the quasi-resonant network contains a negligible fraction of the total number of modes in the system whereas above this critical value a finite fraction of the total number of modes in the system are contained in the largest connected component. We argue that this percolation transition should correspond to the transition to turbulence in the system.

preprint2011arXiv

Instantaneous Gelation in Smoluchowski's Coagulation Equation Revisited

We study the solutions of the Smoluchowski coagulation equation with a regularisation term which removes clusters from the system when their mass exceeds a specified cut-off size, M. We focus primarily on collision kernels which would exhibit an instantaneous gelation transition in the absence of any regularisation. Numerical simulations demonstrate that for such kernels with monodisperse initial data, the regularised gelation time decreases as M increases, consistent with the expectation that the gelation time is zero in the unregularised system. This decrease appears to be a logarithmically slow function of M, indicating that instantaneously gelling kernels may still be justifiable as physical models despite the fact that they are highly singular in the absence of a cut-off. We also study the case when a source of monomers is introduced in the regularised system. In this case a stationary state is reached. We present a complete analytic description of this regularised stationary state for the model kernel, K(m_1,m_2)=max{m_1,m_2}^v, which gels instantaneously when M tends to infinity if v>1. The stationary cluster size distribution decays as a stretched exponential for small cluster sizes and crosses over to a power law decay with exponent v for large cluster sizes. The total particle density in the stationary state slowly vanishes as (Log M^(v-1))^-1/2 when M gets large. The approach to the stationary state is non-trivial : oscillations about the stationary state emerge from the interplay between the monomer injection and the cut-off, M, which decay very slowly when M is large. A quantitative analysis of these oscillations is provided for the addition model which describes the situation in which clusters can only grow by absorbing monomers.

preprint2011arXiv

Mixed flux-equipartition solutions of a diffusion model of nonlinear cascades

We present a parametric study of a nonlinear diffusion equation which generalises Leith's model of a turbulent cascade to an arbitrary cascade having a single conserved quantity. There are three stationary regimes depending on whether the Kolmogorov exponent is greater than, less than or equal to the equilibrium exponent. In the first regime, the large scale spectrum scales with the Kolmogorov exponent. In the second regime, the large scale spectrum scales with the equilibrium exponent so the system appears to be at equilibrium at large scales. Furthermore, in this equilibrium-like regime, the amplitude of the large-scale spectrum depends on the small scale cut-off. This is interpreted as an analogue of cascade nonlocality. In the third regime, the equilibrium spectrum acquires a logarithmic correction. An exact analysis of the self-similar, non-stationary problem shows that time-evolving cascades have direct analogues of these three regimes.

preprint2011arXiv

Some remarks on the inverse Smoluchowski problem for cluster-cluster aggregation

It is proposed to revisit the inverse problem associated with Smoluchowski's coagulation equation. The objective is to reconstruct the functional form of the collision kernel from observations of the time evolution of the cluster size distribution. A regularised least squares method originally proposed by Wright and Ramkrishna (1992) based on the assumption of self-similarity is implemented and tested on numerical data generated for a range of different collision kernels. This method expands the collision kernel as a sum of orthogonal polynomials and works best when the kernel can be expressed exactly in terms of these polynomials. It is shown that plotting an "L-curve" can provide an a-priori understanding of the optimal value of the regularisation parameter and the reliability of the inversion procedure. For kernels which are not exactly expressible in terms of the orthogonal polynomials it is found empirically that the performance of the method can be enhanced by choosing a more complex regularisation function.

preprint2010arXiv

Developing Homogeneous Isotropic Turbulence

We investigate the self-similar evolution of the transient energy spectrum which precedes the establishment of the Kolmogorov spectrum in homogeneous isotropic turbulence in three dimensions using the EDQNM closure model. The transient evolution exhibits self-similarity of the second kind and has a non-trivial dynamical scaling exponent which results in the transient spectrum having a scaling which is steeper than the Kolmogorov $k^{-5/3}$ spectrum. Attempts to detect a similar phenomenon in DNS data are inconclusive owing to the limited range of scales available.

preprint2010arXiv

Nonlocal wave turbulence in the Charney-Hasegawa-Mima equation: a short review

Rossby wave turbulence, as modelled by the Charney-Hasegawa-Mima (CHM) equation, is nonlocal in scale. As a result, the formal stationary Kolmogorov-Zakharov solutions of the Rossby wave kinetic equation, which describe local cascades, are not valid. Rather the solution of the kinetic equation is dominated by interactions between the large and small scales. This suggests an alternative analytic approach based on an expansion of the collision integral in a small parameter obtained from scale separation. This expansion approximates the integral collision operator in the kinetic equation by anisotropic diffusion operators in wavenumber space as first shown in a series of papers by Balk, Nazarenko and Zakharov in the early 1990's. In this note we summarize the foundations of this theory and provide the technical details which were absent from the original papers.

preprint2010arXiv

On the Non-equilibrium Phase Transition in Evaporation-Deposition Models

We study a system of diffusing-aggregating particles with deposition and evaporation of monomers. By combining theoretical and numerical methods, we establish a clearer understanding of the non-equilibrium phase transition known to occur in such systems. The transition is between a growing phase in which the total mass increases for all time and a non-growing phase in which the total mass is bounded. In addition to deriving rigorous bounds on the position of the transition point, we show that the growing phase is in the same universality class as diffusion-aggregation models with deposition but no evaporation. In this regime, the flux of mass in mass space becomes asymptotically constant (as a function of mass) at large times. The magnitude of this flux depends on the evaporation rate but the fact that it is asymptotically constant does not. The associated constant flux relation exactly determines the scaling of the two-point mass correlation function with mass in all dimensions while higher-order mass correlation functions exhibit nonlinear multi-scaling in dimension less than 2. If the deposition rate is below some critical value, a different stationary state is reached at large times characterised by a global balance between evaporation and deposition with a scale-by-scale balance between the mass fluxes due to aggregation and evaporation. Both the mass distribution and the flux decay exponentially in this regime. Finally, we develop a scaling theory of the model near the critical point, which yields non-trivial scaling laws for the critical two-point mass correlation function with mass. These results are well supported by numerical measurements.

preprint2010arXiv

Scaling properties of one-dimensional cluster-cluster aggregation with Levy diffusion

We present a study of the scaling properties of cluster-cluster aggregation with a source of monomers in the stationary state when the spatial transport of particles occurs by Levy flights. We show that the transition from mean-field statistics to fluctuation-dominated statistics which, for the more commonly considered case of diffusive transport, occurs as the spatial dimension of the system is tuned through two from above, can be mimicked even in one dimension by varying the characteristic exponent, beta, of the the Levy jump length distribution. We also show that the two-point mass correlation function, responsible for the flux of mass in the stationary state, is strongly universal: its scaling exponent is given by the mean field value independent of the spatial dimension and independent of the value of beta. Finally we study numerically the two point spatial correlation function which characterises the structure of the depletion zone around heavy particles in the diffusion limited regime. We find that this correlation function vanishes with a non-trivial fractional power of the separation between particles as this separation goes to zero. We provide a scaling argument for the value of this exponent which is in reasonable agreement with the numerical measurements.

preprint2010arXiv

The life-cycle of drift-wave turbulence driven by small scale instability

We demonstrate theoretically and numerically the zonal-flow/drift-wave feedback mechanism for the LH transition in an idealised 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 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 damping of the large scales.

preprint2009arXiv

Dynamical Scaling and the Finite Capacity Anomaly in 3-Wave Turbulence

We present a systematic study of the dynamical scaling process leading to the establishment of the Kolmogorov--Zakharov (KZ) spectrum in weak 3-wave turbulence. In the finite capacity case, in which the transient spectrum reaches infinite frequency in finite time, the dynamical scaling exponent is anomalous in the sense that it cannot be determined from dimensional considerations. As a consequence, the transient spectrum preceding the establishment of the steady state is steeper than the KZ spectrum. Constant energy flux is actually established from right to left in frequency space after the singularity of the transient solution. From arguments based on entropy production, a steeper transient spectrum is heuristically plausible.

preprint2009arXiv

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

We study the modulational instability of geophysical Rossby and plasma drift waves within the Charney-Hasegawa-Mima (CHM) model both theoretically, using truncated (four-mode and three-mode) models, and numerically, using direct simulations of CHM equation in the Fourier space. The linear theory predicts instability for any amplitude of the primary wave. 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 weaker waves, the maximum growth occurs for off-zonal inclined modulations. For very weak primary waves the unstable waves are close to being in three-wave resonance with the primary wave. The nonlinear theory predicts that the zonal flows generated by the linear instability experience pinching into narrow zonal jets. Our numerical simulations confirm the theoretical predictions of the linear theory as well as of the nonlinear pinching. We find that, for strong primary waves, these narrow zonal jets further roll up into Karman-like vortex streets. On the other hand, for weak primary waves, the growth of the unstable mode reverses and the system oscillates between a dominant jet and a dominate primary wave. The 2D vortex streets appear to be more stable than purely 1D zonal jets, and their zonal-averaged speed can reach amplitudes much stronger than is allowed by the Rayleigh-Kuo instability criterion for the 1D case. We find that the truncation models work well for both the linear stage and and often even for the medium-term nonlinear behavior. In the long term, the system transitions to turbulence helped by the vortex-pairing instability (for strong waves) and by the resonant wave-wave interactions (for weak waves).

preprint2007arXiv

Constant flux relation for driven dissipative systems

Conservation laws constrain the stationary state statistics of driven dissipative systems because the average flux of a conserved quantity between driving and dissipation scales should be constant. This requirement leads to a universal scaling law for flux-measuring correlation functions, which generalizes the 4/5-th law of Navier-Stokes turbulence. We demonstrate the utility of this simple idea by deriving new exact scaling relations for models of aggregating particle systems in the fluctuation-dominated regime and for energy and wave action cascades in models of strong wave turbulence.