Search

We consider the problem of droplet impact and droplet spreading on a smooth surface in the case of a viscous Newtonian fluid. We revisit the concept of the rim-lamella model, in which the droplet spreading is described by a system of ordinary differential equations (ODEs). We show that these models contain a singularity which needs to be regularized to produce smooth solutions, and we explore the different regularization techniques in detail. We adopt one particular technique for further investigation, and we use differential inequalities to derive upper and lower bounds for the maximum spreading radius $\mathcal{R}_{max}$. Without having to resort to dimensional analysis or scaling arguments, our bounds confirm the leading-order behaviour $\mathcal{R}_{max} \sim \mathrm{Re}^{1/5}$, as well as a correction to the leading-order behaviour involving the combination $\mathrm{We}^{-1/2}\mathrm{Re}^{2/5}$, in agreement with the experimental literature on droplet spread. Our bounds are also consistent with full numerical solutions of the ODEs, as well as with energy-budget calculations, once head loss is accounted for in the latter., Comment: 25 pages, 9 figures

We consider the classical problem of droplet impact and droplet spread on a smooth surface in the case of an ideal inviscid fluid. We revisit the rim-lamella model of Roisman et al. [\textit{Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences}, 458(2022), pp.1411-1430.]. This model comprises a system of ordinary differential equations (ODEs); we present a rigorous theoretical analysis of these ODEs, and derive upper and lower bounds for the maximum spreading radius. Both bounds possess a $\mathrm{We}^{1/2}$ scaling behaviour, and by a sandwich result, the spreading radius itself also possesses this scaling. We demonstrate rigorously that the rim-lamella model is self-consistent: once a rim forms, its height will invariably exceed that of the lamella. We introduce a rational procedure to obtain initial conditions for the rim-lamella model. Our approach to solving the rim-lamella model gives predictions for the maximum droplet spread that are in close agreement with existing experimental studies and direct numerical simulations., Comment: 24 pages, 7 figures

In this paper we investigate analytically and numerically the nonlinear Kelvin lattice, namely a chain of masses and nonlinear springs, as in the alpha-Fermi-Pasta-Ulam-Tsingou (FPUT) chain, where, in addition, each mass is connected to a nonlinear resonator, i.e., a second mass free to oscillate. Both nonlinearities are quadratic in the equations of motion. This setup represents the simplest prototype of nonlinear wave propagation on a nonlinear metamaterial. In the linear case, we diagonalize the system, and the two branches of the dispersion relation can be found. Using this result, we derive in the nonlinear case the equations of motion for the normal variables in Fourier space, obtaining a system governed by triad interactions among the two branches of the dispersion relation. We find that the transfer of energy between these two branches is ruled by three- and four-wave resonant interactions. We perform numerical simulations of the primitive equations of motion and highlight the role of resonances as an efficient mechanism for transferring energy. Moreover, as predicted by the theory, we provide direct evidence that four-wave resonances appear on a time scale that is longer than the time scale for three-wave resonances. We also assess the recurrence behaviour (usual in the FPUT system) for the nonlinear Kelvin lattice, and we show that, while recurrence is observed if all the energy is placed, at time t=0, in the lowest mode of the acoustical branch, a non-recurrent behaviour is observed if the initial energy is located in the optical branch., Comment: Author Accepted Manuscript version (CNSNS)

Laparoscopy is an electrosurgical medical operation often involving an application of high-frequency alternating current to remove undesired biological tissue from the insufflated abdomen accessible through inlet and outlets trocars. One of the main byproducts in this process are the gaseous particles, called surgical smoke, which is found hazardous for both the patient and the operating room staff. The elimination of this hazardous material is an area of active research in the medical community. Thus, understanding dynamics influenced by the underlying flow inside the abdomen is crucial. In this article, we propose a computational fluid dynamics model and analyse the velocity field in an insufflated abdomen shaped domain by identifying the Lagrangian Coherent Structures (LCS) that are responsible for the transportation, mixing and accumulation of the material particles in the flow. By calculating the mixing strength we show that the regions revealed by these material curves are dependent on the angle, positions and number of the outlets and inlets. Hence, a novel utility of LCS in medical surgery is presented that can detail the dynamics of surgical smoke informing better design of effective smoke removal technologies.

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., Comment: Accepted for publication in Philosophical Transactions A

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., Comment: 6 pages, 3 figures. Main changes: Improved introduction and model discussion to emphasize the deterministic system that supports a turbulent attractor. Improved figures and text to introduce a time-dependent order parameter

In this paper we present a methodology for data accesses when solving batches of Tridiagonal and Pentadiagonal matrices that all share the same left-hand-side (LHS) matrix. The intended application is to the numerical solution of Partial Differential Equations via the finite-difference method, although the methodology is applicable more broadly. By only storing one copy of this matrix, a significant reduction in storage overheads is obtained, together with a corresponding decrease in compute time. Taken together, these two performance enhancements lead to an overall more efficient implementation over the current state of the art algorithms cuThomasBatch and cuPentBatch, allowing for a greater number of systems to be solved on a single GPU. We demonstrate the methodology in the case of the Diffusion Equation, Hyperdiffusion Equation, and the Cahn--Hilliard Equation, all in one spatial dimension. In this last example, we demonstrate how the method can be used to perform $2^{20}$ independent simulations of phase separation in one dimension. In this way, we build up a robust statistical description of the coarsening phenomenon which is the defining behavior of phase separation. We anticipate that the method will be of further use in other similar contexts requiring statistical simulation of physical systems., Comment: arXiv admin note: substantial text overlap with arXiv:1909.04539

The goal of this study is to analyze the fine structure of nonlinear modal interactions in different 1D Burgers and 3D Navier-Stokes flows. This analysis is focused on preferential alignments characterizing the phases of Fourier modes participating in triadic interactions, which are key to determining the nature of energy fluxes between different scales. We develop novel diagnostic tools designed to probe the level of coherence among triadic interactions realizing different flow scenarios. We consider extreme 1D viscous Burgers flows and 3D Navier-Stokes flows which are complemented by singularity-forming inviscid Burgers flows as well as viscous Burgers flows and Navier-Stokes flows corresponding to generic turbulent and simple unimodal initial data, such as the Taylor-Green vortex. The main finding is that while the extreme viscous Burgers and Navier-Stokes flows reveal the same relative level of enstrophy amplification by nonlinear effects, this behaviour is realized via modal interactions with vastly different levels of coherence. In the viscous Burgers flows the flux-carrying triads have phase values which saturate the nonlinearity thereby maximizing the energy flux towards small scales. On the other hand, in 3D Navier-Stokes flows with the extreme initial data the energy flux to small scales is realized by a very small subset of helical triads. The second main finding concerns the role of initial coherence. Comparison of the flows resulting from the extreme and generic initial conditions shows striking similarities between these two types of flows, for the 1D viscous Burgers equation as well as the 3D Navier-Stokes equation., Comment: 57 pages, 31 figures

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.

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., Comment: Accepted and published version (Entropy 2020, 22, 454)

In this paper we present a new methodology for data accesses when solving batches of Tridiagonal and Pentadiagonal matrices that all share the same LHS matrix. By only storing one copy of this matrix there is a significant reduction in storage overheads and the authors show that there is also a performance increase in terms of compute time. These two results combined lead to an overall more efficient implementation over the current state of the art algorithms cuThomasBatch and cuPentBatch, allowing for a greater number of systems to be solved on a single GPU.

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.

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.

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.

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.

We apply Wave Turbulence theory to describe the dynamics on nonlinear one-dimensional chains. We consider $\alpha$ and $\beta$ Fermi-Pasta-Ulam-Tsingou (FPUT) systems, and the discrete nonlinear Klein-Gordon chain. We demonstrate that resonances are responsible for the irreversible transfer of energy among the Fourier modes. We predict that all the systems thermalize for large times, and that the equipartition time scales as a power-law of the strength of the nonlinearity. Our methodology is not limited to only these systems and can be applied to the case of a finite number of modes, such as in the original FPUT experiment, or to the thermodynamic limit, i.e. when the number of modes approach infinity. In the latter limit, we perform state of the art numerical simulations and show that the results are consistent with theoretical predictions. We suggest that the route to thermalization, based only on the presence of exact resonance, has universal features. Moreover, a by-product of our analysis is the asymptotic integrability, up to four wave interactions, of the discrete nonlinear Klein-Gordon chain.

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.

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.

During the 125th European Study Group with Industry held in Limassol, Cyprus, 5-9 December 2016, one of the participating companies, Engino.net Ltd, posed a very interesting challenge to the members of the study group. Engino.net Ltd is a Cypriot company, founded in 2004, that produces a series of toy sets -- the Engino$^{\circledR}$ toy sets -- consisting of a number of building blocks which can be assembled by pupils to compose toy models. Depending on the contents of a particular toy set, the company has developed a number of models that can be built utilizing the blocks present in the set, however the production of a step-by-step assembly manual for each model could only be done manually. The goal of the challenge posed by the company was to implement a procedure to automatically generate the assembly instructions for a given toy. In the present paper we propose a graph-theoretic approach to model the problem and provide a series of results to solve it by employing modified versions of well established algorithms in graph theory. An algorithmic procedure to obtain a hierarchical, physically feasible decomposition of a given toy model, from which the assembly instructions can be recovered, is proposed., Comment: Final version, accepted for publication in the European Journal of Applied Mathematics

We present a theoretical and numerical study of Fourier space triad phase dynamics in one-dimensional stochastically forced Burgers equation at Reynolds number $\mathrm{Re} \approx 2.7 \times 10^4$. We demonstrate that Fourier triad phases over the inertial range display a collective behaviour characterised by intermittent periods of synchronisation and alignment, reminiscent of Kuramoto model (Kuramoto 1984) and directly related to collisions of shocks in physical space. These periods of synchronisation favour efficient energy fluxes across the inertial range towards small scales, resulting in strong bursts of dissipation and enhanced coherence of Fourier energy spectrum. The fast time scale of the onset of synchronisation relegates energy dynamics to a passive role: this is further examined using a reduced system with the Fourier amplitudes fixed in time -- a phase-only model. We show that intermittent triad phase dynamics persists without amplitude evolution and we broadly recover many of the characteristics of the full Burgers system. In addition, for both full Burgers and phase-only systems the physical space velocity statistics reveals that triad phase alignment is directly related to the non-Gaussian statistics typically associated with structure-function intermittency in turbulent systems., Comment: Accepted for publication in Journal of Fluid Mechanics

We describe the theory and present numerical evidence for a new type of nonlinear resonant interaction between gravity waves on the surface of deep water. The resonance constitutes a generalisation of the usual 'exact' resonance as we show that exchanges of energy between the waves can be enhanced when the interaction is three-wave rather than four and the linear frequency mismatch, or detuning, is non-zero i.e. $\omega_1\pm\omega_2\pm\omega_3 \neq0.$ This is possible because the resonance condition is now a match between the so-called 'precession frequency' of a given $\textit{triad interaction}$ and an existent nonlinear frequency in the system. In the limit of weak nonlinearity this precession frequency is simply due to the linear 'drift' of the triad phase; therefore, it tends toward the detuning. This means precession resonance of this type can occur at finite amplitudes, with nonlinear corrections contributing to the resonance. We report energy transfer efficiencies of up to 40%, depending on the model options. To the authors' knowledge this represents the first new type of nonlinear resonance in surface gravity waves since the seminal work of Benjamin & Feir (1967).

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.

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}.$

We present a systematic derivation of the Biot-Savart equation from the Nonlinear Schr\"odinger 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{\kappa^2}{8 \pi}\int_{|{\bf s}-{\bf s}'|>\xi_*} \frac{d{\bf s} \cdot d {\bf s}'}{|{\bf s}-{\bf s}'|} \,,$$ with cut-off length $\xi_* \approx 0.3416293 /\sqrt{\rho_0},$ where $\rho_0$ is the background condensate density far from the vortex lines and $\kappa$ is the quantum of circulation.

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.

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.

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., Comment: This paper supersedes arXiv:1305.5517 [physics.flu-dyn]

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.

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

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., Comment: Improved version