Researcher profile

Miguel D. Bustamante

Miguel D. Bustamante contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
11topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

11 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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.