BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Miguel D. Bustamante

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

ResearcherUnclaimed source record
physics.flu-dynnlin.CDnlin.SIphysics.ao-phmath.APastro-ph.SRCryptography and Securityhep-thmath-phmath.AGmath.MPmath.NTmath.OCmath.RAnlin.PSphysics.data-an

Catalog footprint

What is connected

28works
16topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

A minimal phase-coupling model for intermittency in turbulent systems

Turbulent systems exhibit a remarkable multi-scale complexity, in which spatial structures induce scale-dependent statistics with strong departures from Gaussianity. In Fourier space, this is reflected by pronounced phase synchronization. A quantitative relation between real-space structure, statistics, and phase synchronization is currently missing. Here, we address this problem in the framework of a minimal phase-coupling model, which enables a detailed investigation by means of dynamical systems theory and multi-scale high-resolution simulations. We identify the spectral power-law steepness, which controls the phase coupling, as the control parameter for tuning the non-Gaussian properties of the system. Whereas both very steep and very shallow spectra exhibit close-to-Gaussian statistics, the strongest departures are observed for intermediate slopes comparable to the ones in hydrodynamic and Burgers turbulence. We show that the non-Gaussian regime of the model coincides with a collapse of the dynamical system to a lower-dimensional attractor and the emergence of phase synchronization, thereby establishing a dynamical-systems perspective on turbulent intermittency.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

On the role of continuous symmetries in the solution of the 3D Euler fluid equations and related models

We review and apply the continuous symmetry approach to find the solution of the 3D Euler fluid equations in several instances of interest, via the construction of constants of motion and infinitesimal symmetries, without recourse to Noether's theorem. We show that the vorticity field is a symmetry of the flow, so if the flow admits another symmetry then a Lie algebra of new symmetries can be constructed. For steady Euler flows this leads directly to the distinction of (non-)Beltrami flows: an example is given where the topology of the spatial manifold determines whether extra symmetries can be constructed. Next, we study the stagnation-point-type exact solution of the 3D Euler fluid equations introduced by Gibbon et al. (Physica D, vol.132, 1999, pp.497-510) along with a one-parameter generalisation of it introduced by Mulungye et al. (J. Fluid Mech., vol.771, 2015, pp.468-502). Applying the symmetry approach to these models allows for the explicit integration of the fields along pathlines, revealing a fine structure of blowup for the vorticity, its stretching rate, and the back-to-labels map, depending on the value of the free parameter and on the initial conditions. Finally, we produce explicit blowup exponents and prefactors for a generic type of initial conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Determining when an algebra is an evolution algebra

Evolution algebras are non-associative algebras that describe non-Mendelian hereditary processes and have connections with many other areas. In this paper we obtain necessary and sufficient conditions for a given algebra $A$ to be an evolution algebra. We prove that the problem is equivalent to the so-called $SDC$ $problem$, that is, the $simultaneous$ $diagonalisation$ $via$ $congruence$ of a given set of matrices. More precisely we show that an $n$-dimensional algebra $A$ is an evolution algebra if, and only if, a certain set of $n$ symmetric $n\times n$ matrices $\{M_{1}, \ldots, M_{n}\}$ describing the product of $A$ are $SDC$. We apply this characterisation to show that while certain classical genetic algebras (representing Mendelian and auto-tetraploid inheritance) are not themselves evolution algebras, arbitrarily small perturbations of these are evolution algebras. This is intriguing as evolution algebras model asexual reproduction unlike the classical ones.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Solving the problem of simultaneous diagonalization of complex symmetric matrices via congruence

We provide a solution to the problem of simultaneous $diagonalization$ $via$ $congruence$ of a given set of $m$ complex symmetric $n\times n$ matrices $\{A_{1},\ldots,A_{m}\}$, by showing that it can be reduced to a possibly lower-dimensional problem where the question is rephrased in terms of the classical problem of simultaneous $diagonalization$ $via$ $similarity$ of a new related set of matrices. We provide a procedure to determine in a finite number of steps whether or not a set of matrices is simultaneously diagonalizable by congruence. This solves a long standing problem in the complex case.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Image Encryption Using Elliptic Curves and Rossby/Drift Wave Triads

We propose an image encryption scheme based on quasi-resonant Rossby/drift wave triads (related to elliptic surfaces) and Mordell elliptic curves (MECs). By defining a total order on quasi-resonant triads, at a first stage we construct quasi-resonant triads using auxiliary parameters of elliptic surfaces in order to generate pseudo-random numbers. At a second stage, we employ an MEC to construct a dynamic substitution box (S-box) for the plain image. The generated pseudo-random numbers and S-box are used to provide diffusion and confusion, respectively, in the tested image. We test the proposed scheme against well-known attacks by encrypting all gray images taken from the USC-SIPI image database. Our experimental results indicate the high security of the newly developed scheme. Finally, via extensive comparisons we show that the new scheme outperforms other popular schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Nonholonomic Noetherian symmetries and integrals of the Routh sphere and Chaplygin ball

The dynamics of a spherical body with a non-uniform mass distribution rolling on a plane were discussed by Sergey Chaplygin, whose 150th anniversary we celebrate this year. The Chaplygin top is a non-integrable system, with a colourful range of interesting motions. A special case of this system was studied by Edward Routh, who showed that it is integrable. The Routh sphere has centre of mass offset from the geometric centre, but it has an axis of symmetry through both these points, and equal moments of inertia about all axes orthogonal to the symmetry axis. There are three constants of motion: the total energy and two quantities involving the angular momenta. It is straightforward to demonstrate that these quantities, known as the Jellett and Routh constants, are integrals of the motion. However, their physical significance has not been fully understood. In this paper, we show how the integrals of the Routh sphere arise from Emmy Noether's invariance identity. We derive expressions for the infinitesimal symmetry transformations associated with these constants. We find the finite version of these symmetries and provide their geometrical interpretation. As a further demonstration of the power and utility of this method, we find the Noether symmetries and corresponding Noether integrals for a system introduced recently: the Chaplygin ball on a rotating turntable, confirming that the known integrals are directly obtained from Noether's theorem.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Nonlinear Rossby wave-wave and wave-mean flow theory for long term Solar cycle modulations

The Schwabe cycle of solar activity exhibits modulations and frequency fluctuations on slow time scales of centuries and millennia. Plausible physical explanations for the cause of these long-term variations of the solar cycle are still elusive, with possible theories including stochasticity of alpha effect and fluctuations of the differential rotation. It has been suggested recently in the literature that there exists a possible relation between the spatio-temporal structure of Solar cycle and the nonlinear dynamics of magnetohydrodynamic Rossby waves at the solar tachocline, including both wave-wave and wave-mean flow interactions. Here we extend the nonlinear theory of MHD Rossby waves presented in a previous article to take into account long term modulation effects due to a recently discovered mechanism that allows significant energy transfers throughout different wave triads: the precession resonance mechanism. We have found a large number of Rossby-Haurwitz wave triads whose frequency mismatches are compatible with the solar cycle frequency. Consequently, by analyzing the reduced dynamics of two triads coupled by a single mode (five-wave system), we have demonstrated that in the amplitude regime in which precession resonance occurs, the energy transfer throughout the system yields significant long-term modulations on the main $\sim 11$yr period associated with intra-triad energy exchanges. We further show that such modulations display an inverse relationship between the characteristic wave amplitude and the period of intra-triad energy exchanges, which is consistent with the Waldmeier's law for the solar cycle. In the presence of a constant forcing and dissipation, the five-wave system in the precession resonance regime exhibits irregular amplitude fluctuations with some periods resembling the Grand Minimum states.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

On the convergence of the normal form transformation in discrete Rossby and drift wave turbulence

We study numerically the region of convergence of the normal form transformation for the case of the Charney-Hasagawa-Mima (CHM) equation to investigate whether certain finite amplitude effects can be described in normal coordinates. We do this by taking a Galerkin truncation of four Fourier modes making part of two triads: one resonant and one non-resonant, joined together by two common modes. We calculate the normal form transformation directly from the equations of motion of our reduced model, successively applying the algorithm to calculate the transformation up to $7^\textrm{th}$ order to eliminate all non-resonant terms, and keeping up to $8$-wave resonances. We find that the amplitudes at which the normal form transformation diverge very closely match with the amplitudes at which a finite-amplitude phenomenon called $precession$ $resonance$ (Bustamante $et$ $al.$ 2014) occurs, characterised by strong energy transfers. This implies that the precession resonance mechanism cannot be explained using the usual methods of normal forms in wave turbulence theory, so a more general theory for intermediate nonlinearity is required.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Exact discrete resonances in the Fermi-Pasta-Ulam-Tsingou system

In systems of N coupled anharmonic oscillators, exact resonant interactions play an important role in the energy exchange between normal modes. In the weakly nonlinear regime, those interactions may facilitate energy equipartition in Fourier space. We consider analytically resonant wave-wave interactions for the celebrated Fermi-Pasta-Ulam-Tsingou (FPUT) system. Using a number-theoretical approach based on cyclotomic polynomials, we show that the problem of finding exact resonances for a system of N particles is equivalent to a Diophantine equation whose solutions depend sensitively on the set of divisors of N. We provide an algorithm to construct all possible resonances, based on two methods: pairing-off and cyclotomic, which we introduce to build up explicit solutions to the 4-, 5- and 6-wave resonant conditions. Our results shed some light in the understanding of the long-standing FPUT paradox, regarding the sensitivity of the resonant manifolds with respect to the number of particles N and the corresponding time scale of the interactions leading to thermalisation. In this light we demonstrate that 6-wave resonances always exist for any N, while 5-wave resonances exist if N is divisible by 3 and N > 6. It is known (for finite N) that 4-wave resonances do not mix energy across the spectrum, so we investigate whether 5-wave resonances can produce energy mixing across a significant region of the Fourier spectrum by analysing the interconnected network of Fourier modes that can interact nonlinearly via resonances. The answer depends on the set of odd divisors of N that are not divisible by 3: the size of this set determines the number of dynamically independent components, corresponding to independent constants of motion (energies). We show that 6-wave resonances connect all these independent components, providing in principle a restoring mechanism for full-scale thermalisation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Derivation of the Biot-Savart equation from the Nonlinear Schrödinger equation

We present a systematic derivation of the Biot-Savart equation from the Nonlinear Schrödinger equation, in the limit when the curvature radius of vortex lines and the inter-vortex distance are much greater than the vortex healing length, or core radius. We derive the Biot-Savart equations in Hamiltonian form with Hamiltonian expressed in terms of vortex lines, $$ H= \frac{κ^2}{8 π}\int_{|{\bf s}-{\bf s}'|>ξ_*} \frac{d{\bf s} \cdot d {\bf s}'}{|{\bf s}-{\bf s}'|} \,,$$ with cut-off length $ξ_* \approx 0.3416293 /\sqrt{ρ_0},$ where $ρ_0$ is the background condensate density far from the vortex lines and $κ$ is the quantum of circulation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Phase and precession evolution in the Burgers equation

We present a phenomenological study of the phase dynamics of the one-dimensional stochastically forced Burgers equation, and of the same equation under a Fourier mode reduction on a fractal set. We study the connection between coherent structures in real space and the evolution of triads in Fourier space. Concerning the one-dimensional case, we find that triad phases show alignments and synchronisations that favour energy fluxes towards small scales --a direct cascade. In addition, strongly dissipative real-space structures are associated with entangled correlations amongst the phase precession frequencies and the amplitude evolution of Fourier triads. As a result, triad precession frequencies show a non-Gaussian distribution with multiple peaks and fat tails, and there is a significant correlation between triad precession frequencies and amplitude growth. Links with dynamical systems approach are briefly discussed, such as the role of unstable critical points in state space. On the other hand, by reducing the fractal dimension $D$ of the underlying Fourier set, we observe: i) a tendency toward a more Gaussian statistics, ii) a loss of alignment of triad phases leading to a depletion of the energy flux, and iii) the simultaneous reduction of the correlation between the growth of Fourier mode amplitudes and the precession frequencies of triad phases.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Robust Quantitative Comparison Criterion of Two Signals based on the Sobolev Norm of Their Difference

In this manuscript we present a method for the quantitative comparison of two surfaces, applicable to temporal and/or spatial extent in one or two dimensions. Often surface comparisons are simply overlaid graphs of results from different methodologies that are qualitative at best; it is the purpose of this work to facilitate quantitative evaluation. The surfaces can be analytical, numerical, and/or experimental, and the result returned by the method, termed surface similarity parameter or normalized error, has been normalized so that its value lies between zero and one. When the parameter has a value of zero, the surfaces are in perfect agreement, whereas a value of one is indicative of perfect disagreement. To provide insight regarding the magnitude of the parameter, several canonical cases are presented, followed by results from breaking water wave experimental measurements with numerical simulations, and by a comparison of a prescribed, periodic, square-wave surface profile and the subsequent manufactured surface.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Atypical late-time singular regimes accurately diagnosed in stagnation-point-type solutions of 3D Euler flows

We revisit, both numerically and analytically, the finite-time blowup of the infinite-energy solution of 3D Euler equations of stagnation-point-type introduced by Gibbon et al. (1999). By employing the method of mapping to regular systems, presented in Bustamante (2011) and extended to the symmetry-plane case by Mulungye et al. (2015), we establish a curious property of this solution that was not observed in early studies: before but near singularity time, the blowup goes from a fast transient to a slower regime that is well resolved spectrally, even at mid-resolutions of $512^2.$ This late-time regime has an atypical spectrum: it is Gaussian rather than exponential in the wavenumbers. The analyticity-strip width decays to zero in a finite time, albeit so slowly that it remains well above the collocation-point scale for all simulation times $t < T^* - 10^{-9000}$, where $T^*$ is the singularity time. Reaching such a proximity to singularity time is not possible in the original temporal variable, because floating point double precision ($\approx 10^{-16}$) creates a `machine-epsilon' barrier. Due to this limitation on the \emph{original} independent variable, the mapped variables now provide an improved assessment of the relevant blowup quantities, crucially with acceptable accuracy at an unprecedented closeness to the singularity time: $T^*- t \approx 10^{-140}.$

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Symmetry-plane model of 3D Euler flows and mapping to regular systems to improve blowup assessment using numerical and analytical solutions

Motivated by the work on stagnation-point type exact solutions (with infinite energy) of 3D Euler fluid equations by Gibbon et al. (1999) and the subsequent demonstration of finite-time blowup by Constantin (2006) we introduce a one-parameter family of models of the 3D Euler fluid equations on a 2D symmetry plane. Our models are seen as a deformation of the 3D Euler equations which respects the variational structure of the original equations so that explicit solutions can be found for the supremum norms of the basic fields. The value of the model's parameter determines if there is finite-time blowup, and the singularity time can be computed explicitly in terms of the initial conditions and the model's parameter. We use a representative parameter value, for which the solution blows up in finite-time, as a benchmark for the systematic study of errors in numerical solutions. We compare numerical integrations of our "original" model with a "mapped" version of these equations. The mapped version is a globally regular (in time) system of equations, obtained via a bijective nonlinear mapping of time and fields from the original model equations (Bustamante 2011). We show that the mapped system's numerical solution increases accuracy in estimates of supremum norms and singularity time while entailing only a small additional computational cost. We study the Fourier spectrum of the model's numerical solution and find that the analyticity-strip width (a measure of the solution's analyticity) tends to zero as a power law in a finite time. This is in agreement with the finite-time blowup of the supremum norms, in the light of rigorous bounds stemming from the bridge (Bustamante & Brachet 2012) between the analyticity-strip method and the Beale-Kato-Majda type of theorems. We conclude by discussing the implications of this research on the analysis of numerical solutions of the 3D Euler fluid equations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Robust energy transfer mechanism via precession resonance in nonlinear turbulent wave systems

A robust energy transfer mechanism is found in nonlinear wave systems, which favours transfers towards modes interacting via triads with nonzero frequency mismatch, applicable in meteorology, nonlinear optics and plasma wave turbulence. We introduce the concepts of truly dynamical degrees of freedom and triad precession. Transfer efficiency is maximal when the triads' precession frequencies resonate with the system's nonlinear frequencies, leading to a collective state of synchronised triads with strong turbulent cascades at intermediate nonlinearity. Numerical simulations confirm analytical predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Complete classification of discrete resonant Rossby/drift wave triads on periodic domains

We consider the set of Diophantine equations that arise in the context of the barotropic vorticity equation on periodic domains, when nonlinear wave interactions are studied to leading order in the amplitudes. The solutions to this set of Diophantine equations are of interest in atmosphere (Rossby waves) and Tokamak plasmas (drift waves), because they provide the values of the spectral wavevectors that interact resonantly via three-wave interactions. These come in "triads", i.e., groups of three wavevectors. We provide the full solution to the Diophantine equations in the case of infinite Rossby deformation radius. The method is completely new, and relies on mapping the unknown variables to rational points on quadratic forms of "Minkowski" type. Classical methods invented centuries ago by Fermat, Euler, Lagrange and Minkowski, are used to classify all solutions to our original Diophantine equations, thus providing a computational method to generate numerically all the resonant triads in the system. Our method has a clear computational advantage over brute-force numerical search: on a 10000^2 grid, the brute-force search would take 15 years using optimised C++ codes, whereas our method takes about 40 minutes. The method is extended to generate quasi-resonant triads, which are defined by relaxing the resonant condition on the frequencies, allowing for a small mismatch. Quasi-resonances are robust with respect to physical perturbations, unlike exact resonances. Therefore, the new method is really valuable in practical terms. We show that the set of quasi-resonances form an intricate network of clusters of connected triads, whose structure depends on the value of the allowed mismatch. We provide some quantitative comparison between the clusters' structure and the onset of fully nonlinear turbulence in the barotropic vorticity equation, and provide perspectives for new research.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Counting of discrete Rossby/drift wave resonant triads

The purpose of this note is to remove the confusion about counting of resonant wave triads for Rossby and drift waves in the context of the Charney-Hasegawa-Mima equation. In particular, we aim to point out a major error of over-counting of triads in the paper "Discrete exact and quasi-resonances of Rossby/drift waves on beta-plane with periodic boundary conditions", by Kartashov and Kartashova, arXiv:1307.8272v1 [physics.flu-dyn] (2013).

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Counting of discrete Rossby/drift wave resonant triads (again)

The purpose of our earlier note (arXiv:1309.0405 [physics.flu-dyn]) was to remove the confusion over counting of resonant wave triads for Rossby and drift waves in the context of the Charney-Hasegawa-Mima equation. A comment by Kartashov and Kartashova (arXiv:1309.0992v1 [physics.flu-dyn]) on that note has further confused the situation. The present note aims to remove this obfuscation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Quadratic invariants for discrete clusters of weakly interacting waves

We consider discrete clusters of quasi-resonant triads arising from a Hamiltonian three-wave equation. A cluster consists of N modes forming a total of M connected triads. We investigate the problem of constructing a functionally independent set of quadratic constants of motion. We show that this problem is equivalent to an underlying basic linear problem, consisting of finding the null space of a rectangular (M times N) matrix A with entries 1, -1 and 0. In particular, we prove that the number of independent quadratic invariants is equal to J = N - M* (greater than or equal to N - M), where M* is the number of linearly independent rows in A. Thus, the problem of finding all independent quadratic invariants is reduced to a linear algebra problem in the Hamiltonian case. We establish that the properties of the quadratic invariants (e.g., locality) are related to the topological properties of the clusters (e.g., types of linkage). To do so, we formulate an algorithm for decomposing large clusters into smaller ones and show how various invariants are related to certain parts of a cluster, including the basic structures leading to M* < M. We illustrate our findings by presenting examples from the Charney-Hasegawa-Mima wave model, and by showing a classification of small (up to three-triad) clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Robust energy transfer mechanism and critically balanced turbulence via non-resonant triads in nonlinear wave systems

A robust energy transfer mechanism is found in nonlinear wave systems, which favours transfers towards modes interacting via non-resonant triads, applicable in meteorology, nonlinear optics and plasma wave turbulence. Transfer efficiency is maximal when the frequency mismatch of the non-resonant triad balances the system's nonlinear frequency: at intermediate levels of oscillation amplitudes an instability is triggered that explores unstable manifolds of periodic orbits, so turbulent cascades are most efficient at intermediate nonlinearity. Numerical simulations confirm analytical predictions.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Interplay between the Beale-Kato-Majda theorem and the analyticity-strip method to investigate numerically the incompressible Euler singularity problem

Numerical simulations of the incompressible Euler equations are performed using the Taylor-Green vortex initial conditions and resolutions up to $4096^3$. The results are analyzed in terms of the classical analyticity strip method and Beale, Kato and Majda (BKM) theorem. A well-resolved acceleration of the time-decay of the width of the analyticity strip $δ(t)$ is observed at the highest resolution for $3.7<t<3.85$ while preliminary 3D visualizations show the collision of vortex sheets. The BKM criterium on the power-law growth of supremum of the vorticity, applied on the same time-interval, is not inconsistent with the occurrence of a singularity around $t \simeq 4$. These new findings lead us to investigate how fast the analyticity strip width needs to decrease to zero in order to sustain a finite-time singularity consistent with the BKM theorem. A new simple bound of the supremum norm of vorticity in terms of the energy spectrum is introduced and used to combine the BKM theorem with the analyticity-strip method. It is shown that a finite-time blowup can exist only if $δ(t)$ vanishes sufficiently fast at the singularity time. In particular, if a power law is assumed for $δ(t)$ then its exponent must be greater than some critical value, thus providing a new test that is applied to our $4096^3$ Taylor-Green numerical simulation. Our main conclusion is that the numerical results are not inconsistent with a singularity but that higher-resolution studies are needed to extend the time-interval on which a well-resolved power-law behavior of $δ(t)$ takes place, and check whether the new regime is genuine and not simply a crossover to a faster exponential decay.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

3D Euler equations and ideal MHD mapped to regular systems: probing the finite-time blowup hypothesis

We prove by an explicit construction that solutions to incompressible 3D Euler equations defined in the periodic cube can be mapped bijectively to a new system of equations whose solutions are globally regular. We establish that the usual Beale-Kato-Majda criterion for finite-time singularity (or blowup) of a solution to the 3D Euler system is equivalent to a condition on the corresponding \emph{regular} solution of the new system. In the hypothetical case of Euler finite-time singularity, we provide an explicit formula for the blowup time in terms of the regular solution of the new system. The new system is amenable to being integrated numerically using similar methods as in Euler equations. We propose a method to simulate numerically the new regular system and describe how to use this to draw robust and reliable conclusions on the finite-time singularity problem of Euler equations, based on the conservation of quantities directly related to energy and circulation. The method of mapping to a regular system can be extended to any fluid equation that admits a Beale-Kato-Majda type of theorem, e.g. 3D Navier-Stokes, 2D and 3D magnetohydrodynamics, and 1D inviscid Burgers. We discuss briefly the case of 2D ideal magnetohydrodynamics. In order to illustrate the usefulness of the mapping, we provide a thorough comparison of the analytical solution versus the numerical solution in the case of 1D inviscid Burgers equation.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Resonance clustering in wave turbulent regimes: Integrable dynamics

Two fundamental facts of the modern wave turbulence theory are 1) existence of power energy spectra in $k$-space, and 2) existence of "gaps" in this spectra corresponding to the resonance clustering. Accordingly, three wave turbulent regimes are singled out: \emph{kinetic}, described by wave kinetic equations and power energy spectra; \emph{discrete}, characterized by resonance clustering; and \emph{mesoscopic}, where both types of wave field time evolution coexist. In this paper we study integrable dynamics of resonance clusters appearing in discrete and mesoscopic wave turbulent regimes. Using a novel method based on the notion of dynamical invariant we establish that some of the frequently met clusters are integrable in quadratures for arbitrary initial conditions and some others -- only for particular initial conditions. We also identify chaotic behaviour in some cases. Physical implications of the results obtained are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

3D Euler about a 2D Symmetry Plane

Initial results from new calculations of interacting anti-parallel Euler vortices are presented with the objective of understanding the origins of singular scaling presented by Kerr (1993) and the lack thereof by Hou and Li (2006). Core profiles designed to reproduce the two results are presented, new more robust analysis is proposed, and new criteria for when calculations should be terminated are introduced and compared with classical resolution studies and spectral convergence tests. Most of the analysis is on a $512 \times 128 \times 2048$ mesh, with new analysis on a just completed $1024 \times 256 \times 2048$ used to confirm trends. One might hypothesize that there is a finite-time singularity with enstrophy growth like $Ω\sim (T_c-t)^{-γ_Ω}$ and vorticity growth like $||ω||_\infty \sim (T_c-t)^{-γ}$. The new analysis would then support $γ_Ω\approx 1/2$ and $γ> 1$. These represent modifications of the conclusions of Kerr (1993). Issues that might arise at higher resolution are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Generalized Eulerian-Lagrangian description of Navier-Stokes and resistive MHD dynamics

New generalized equations of motion for the Weber-Clebsch potentials that describe both the Navier-Stokes and MHD dynamics are derived. These depend on a new parameter, which has dimensions of time for Navier-Stokes and inverse velocity for MHD. Direct numerical simulations are performed. For Navier-Stokes, the generalized formalism captures the intense reconnection of vortices of the Boratav, Pelz and Zabusky flow, in agreement with the previous study by Ohkitani and Constantin. For MHD, the new formalism is used to detect magnetic reconnection in several flows: the 3D Arnold, Beltrami and Childress (ABC) flow and the (2D and 3D) Orszag-Tang vortex. It is concluded that periods of intense activity in the magnetic enstrophy are correlated with periods of increasingly frequent resettings. Finally, the positive correlation between the sharpness of the increase in resetting frequency and the spatial localization of the reconnection region is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Multi-Lagrangians, Hereditary Operators and Lax Pairs for the Korteweg-de Vries Positive and Negative Hierarchies

We present an approach to the construction of action principles for differential equations, and apply it to field theory in order to construct systematically, for integrable equations which are based on a Nijenhuis (or hereditary) operator, a ladder of action principles which is complementary to the well-known multi-Hamiltonian formulation. We work out results for the Korteweg-de Vries (KdV) equation, which is a member of the positive hierarchy related to a hereditary operator. Three negative hierarchies of (negative) evolution equations are defined naturally from the hereditary operator as well, in the context of field theory. The Euler-Lagrange equations arising from the action principles are equivalent to the original evolution equation + deformations, which are obtained in terms of the positive and negative evolution vectors. We recognize the Liouville, Sinh-Gordon, Hunter-Zheng and Camassa-Holm equations as negative equations. The ladder for KdV is directly mappable to a ladder for any of these negative equations and other positive equations (e.g., the Harry-Dym and a special case of the Krichever-Novikov equations): a new nonlocal action principle for the deformed system Sinh-Gordon + spatial translation vector is presented. Several nonequivalent, nonlocal time-reparametrization invariant action principles for KdV are constructed. Hamiltonian and Symplectic operators are obtained in factorized form. Alternative Lax pairs for all negative flows are constructed, using the flows and the hereditary operator as only input. From this result we prove that all positive and negative equations in the hierarchies share the same sets of local and nonlocal constants of the motion for KdV, which are explicitly obtained using the local and nonlocal action principles for KdV.

0 citations0 reviewsNo code linkedOpen access