No Arabic abstract
We report results on the explicit parameterisation of discrete Rossby-wave resonant triads of the Charney-Hasegawa-Mima equation in the small-scale limit (i.e. large Rossby deformation radius), following up from our previous solution in terms of elliptic curves (Bustamante and Hayat, 2013). We find an explicit parameterisation of the discrete resonant wavevectors in terms of two rational variables. We show that these new variables are restricted to a bounded region and find this region explicitly. We argue that this can be used to reduce the complexity of a direct numerical search for discrete triad resonances. Also, we introduce a new direct numerical method to search for discrete resonances. This numerical method has complexity ${mathcal{O}}(N^3)$, where $N$ is the largest wavenumber in the search. We apply this new method to find all discrete irreducible resonant triads in the wavevector box of size $5000$, in a calculation that took about $10.5$ days on a $16$-core machine. Finally, based on our method of mapping to elliptic curves, we discuss some dynamical implications regarding the spread of quadratic invariants across scales via resonant triad interactions, in the form of sharp bounds on the size of the interacting wavevectors.
Linear wave solutions to the Charney-Hasegawa-Mima partial differential equation with periodic boundary conditions have two physical interpretations: Rossby (atmospheric) waves, and drift (plasma) waves in a tokamak. These waves display resonance in triads. In the case of infinite Rossby deformation radius, the set of resonant triads may be described as the set of integer solutions to a particular homogeneous Diophantine equation, or as the set of rational points on a projective surface. We give a rational parametrization of the smooth points on this surface, answering the question: What are all resonant triads? We also give a fiberwise description, yielding a procedure to answer the question: For fixed $r in mathbb{Q}$, what are all wavevectors $(x,y)$ that resonate with a wavevector $(a,b)$ with $a/b = r$?
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.
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.
Kraichnan seminal ideas on inverse cascades yielded new tools to study common phenomena in geophysical turbulent flows. In the atmosphere and the oceans, rotation and stratification result in a flow that can be approximated as two-dimensional at very large scales, but which requires considering three-dimensional effects to fully describe turbulent transport processes and non-linear phenomena. Motions can thus be classified into two classes: fast modes consisting of inertia-gravity waves, and slow quasi-geostrophic modes for which the Coriolis force and horizontal pressure gradients are close to balance. In this paper we review previous results on the strength of the inverse cascade in rotating and stratified flows, and then present new results on the effect of varying the strength of rotation and stratification (measured by the ratio $N/f$ of the Brunt-Vaisala frequency to the Coriolis frequency) on the amplitude of the waves and on the flow quasi-geostrophic behavior. We show that the inverse cascade is more efficient in the range of $N/f$ for which resonant triads do not exist, $1/2 le N/f le 2$. We then use the spatio-temporal spectrum, and characterization of the flow temporal and spatial scales, to show that in this range slow modes dominate the dynamics, while the strength of the waves (and their relevance in the flow dynamics) is weaker.
Implementing multicomponent diffusion models in reacting-flow simulations is computationally expensive due to the challenges involved in calculating diffusion coefficients. Instead, mixture-averaged diffusion treatments are typically used to avoid these costs. However, to our knowledge, the accuracy and appropriateness of the mixture-averaged diffusion models has not been verified for three-dimensional turbulent premixed flames. In this study we propose a fast,efficient, low-memory algorithm and use that to evaluate the role of multicomponent mass diffusion in reacting-flow simulations. Direct numerical simulation of these flames is performed by implementing the Stefan-Maxwell equations in NGA. A semi-implicit algorithm decreases the computational expense of inverting the full multicomponent ordinary diffusion array while maintaining accuracy and fidelity. We first verify the method by performing one-dimensional simulations of premixed hydrogen flames and compare with matching cases in Cantera. We demonstrate the algorithm to be stable, and its performance scales approximately with the number of species squared. Then, as an initial study of multicomponent diffusion, we simulate premixed, three-dimensional turbulent hydrogen flames, neglecting secondary Soret and Dufour effects. Simulation conditions are carefully selected to match previously published results and ensure valid comparison. Our results show that using the mixture-averaged diffusion assumption leads to a 15% under-prediction of the normalized turbulent flame speed for a premixed hydrogen-air flame. This difference in the turbulent flame speed motivates further study into using the mixture-averaged diffusion assumption for DNS of moderate-to-high Karlovitz number flames.