Observation of the Inverse Energy Cascade in the modified Korteweg de Vries Equation D. Dutykh and E. Tobisch LAMA, UMR 5127 CNRS, Université de Savoie, Campus Scientifique, 73376 Le Bourget-du-Lac Cedex, France Institute for Analysis, J. Kepler UNiversity, Linz; Austria arxiv:1406.3784v1 [physics.flu-dyn] 14 Jun 2014 In this Letter we demonstrate for the first time the formation of the inverse energy cascade in the focusing modified Korteweg de Vries (mkdv) equation. We study numerically the properties of this cascade such as the dependence of the spectrum shape on the initial excitation parameter (amplitude), perturbation magnitude and the size of the spectral domain. Most importantly we found that the inverse cascade is always accompanied by the direct one and they both form a very stable quasi-stationary structure in the Fourier space in the spirit of the FPU-like reoccurrence phenomenon. The formation of this structure is intrinsically related to the development of the nonlinear stage of the Modulational Instability (MI). These results can be used in several fields such as the internal gravity water waves, ion-acoustic waves in plasmas and others. I. INTRODUCTION Currently there exist three main models which describe the energy cascade in nonlinear wave systems. Each model covers its own range of the nonlinearity parameter. Namely, for strong nonlinearities the Kolmogorov cascade with the widely known exponent k 5/3 (k being the vortex scale) introduced for the hydrodynamic turbulence theory [1]. In the opposite limit of weak nonlinearities, the kinetic weak turbulence spectra k α (now k is the wavenumber) were introduced, [2]; in this case the spectrum exponent α is not universal anymore and the nonlinearity parameter ak 10 2, where a is the typical wave amplitude. For the intermediate values of the nonlinearityak = Ø ( 10 1) the Dynamicenergycascade (D-cascade) model was proposed for the equations which feature the Modulational Instability (MI) property [3]. This construction was originally described and studied in details for the NLS-type equations [4, 5]. In particular it was shown [3] that in the focusing (i.e. possessing the MI) NLS equation only one cascading mode is forming due to the intrinsic narrow band spectrum approximation. This prediction was checked numerically and the development of the MI in the NLS is shown on Figure 1. On this bottom panel one can see only one spike in the Fourier spectrum along with the broadening wings. On the other hand, the modified Korteweg de Vries (mkdv) equation does not have any restrictions on the wave spectrum. In this study we consider the following focusing mkdv equation: mkdv(u). = u t +u xxx +6u 2 u x = 0. (1) The formation of the direct D-cascade in this equation has been recently observed and investigated numerically in [6]. The direct cascade forms starting from any base FIG. 1: Modulational instability in the NLS equation. wavenumber k 0 > 0. The main properties of the direct cascade will be briefly described in the next Section. Concerning the inverse cascade, it was not observed yet to our knowledge. In this Letter we aim at presenting solid numerical evidences towards the existence of the inverse D-cascade and studying some of its properties comparing to the direct one. II. DIRECT CASCADE IN THE MKDV EQUATION For our numerical simulations we adopt the set-up used also earlier in [7, 8] for the study of the modulation instability in the physical space, with the following. initial [ condition ] posed on a periodic domain [ l, l] = π/k0,π/k 0 : u(x,0) u 0 (x) = a ( 1+δsin(K 0 x) ) sin(k 0 x), (2) Electronic address: Elena.Tobisch@jku.at where a is the base wave amplitude, δ is the perturbation magnitude and the wavenumbers k 0, K 0 are chosen such
2 Parameter Range Amplitude a = 10 2 2.4 10 1 Perturbation magnitude δ = 5 10 2 5 10 1 Base wavenumber k 0 = 1.8 60.0 Simulation time horizon T = 40 2000 Fourier harmonics N = 1024 131072 TABLE I: Parameters range used in numerical simulations. that theirratiok 0 /K 0 Z. It isimportantto realizethat the amplitude a is the measure of the nonlinearity in the mkdv equation, in the contrast to the NLS equation where the nonlinearity is related to the wave steepness ak 0 defined in terms of the modulational wavenumber K 0. The main characteristics of the direct D-cascade observed in [6] are as follows: (a)a = 0.08, t = 604 (a) the main structure of the cascade is already observable in Fourier spectra from the first instances of the dynamical evolution; (b) the development of the Modulation Instability (MI) in the physical space occurs at the time scale t ε 2 and corresponds to the spectral broadening of the cascading modes; the frequencies and energies of the cascading modes are quasi-stationary; (c) the energy spectrum of the cascade has exponential form, E k exp( αk) which does not depend on the base amplitude a for fixed values of other parameters; (d) the increase of the spectral domain yields the increase of the cascade length but does not changes the exponent α; (e) the increase of the perturbation magnitude δ from 5% to 50% effects the time scale on which the MI occurs in the physical space: higher perturbation accelerates the process; (f) the increase of the base wavenumber k 0 to 2k 0, 3k 0,..., 30k 0 ; changesthe exponent α and yields at some point (beginning with 35k 0 ) the occurrence of the inverse cascade; The appearance of the direct D-cascade has been clearly observed (characteristic cascades are shown in Figure 2), for the wide range of the initial parameters described in Table I. For the sake of illustration, on Figure 2 we show the formation of the direct D-cascade for the initial condition (2) with the following parameters: k 0 = 1.884, m = 0.05, K 0 = 0.00785 and a = 0.08 on the top panel (a) and a = 0.16 on the panel (b). These values of parameters were also used in [8]. However, the numerical resolution was significantly lower in their study (b)a = 0.16, t = 281 FIG. 2: Direct cascade in the mkdv equation. The energy spectrum exponential fit shown on bottom panels is the same for both amplitudes and equal to e 1.45k. and the authors did not report the evolution of the system in the Fourier space. From our simulations it can be easily seen that an increase in the base wave amplitude does not change the shape of the energy spectrum along with the spacing between cascading modes. Moreover, the time instances t = 50 and 90, where these snapshots were taken, indicate that the MI development goes faster with higher amplitudes. The dashed black lines shown on bottom panels show the same exponential fit for both panels. In the next Section we present the results of our studies of the inverse cascade, using the same parameters as abovebutwithbasewavenumbersk 0 shiftedtowardshort waves.
3 III. DOUBLE CASCADE IN THE MKDV EQUATION In order to solve numerically the mkdv equation on a periodic domain we used the classical Fourier-type pseudo-spectral method along with the 2/3-rule for dealising of nonlinear products. The discretization in time was done with the embedded adaptive 5 th order Cash Karp Runge Kutta scheme with the adaptive PI step size control. In most numerical simulations presented below we will use N = 32768 Fourier harmonics, unless a special remark is made. On Figure 3 we show the simultaneous formation of the direct and inverse cascades (which will be referred below as the double cascade) at different times. In order to achievethis result the basewavenumberk 0 was shifted approximatively into the middle of the spectral domain (k 0 = 35 1.884). One can see that starting from t 20 (see Figure 3(d)) the MI is already fully developed in the physical space, which corresponds to the appearance of a quasi-stationary structure double cascade in the Fourier space. For times t 20 the evolution of the system follows the FPU-like pattern, without any significant change in the Fourier domain (see Figure 3(e, f)). As in the case ofthe isolated direct cascade, the double cascade equally has the exponential decay in the Fourier space. However, the two exponents are slightly different in a way that the direct cascade decays generally faster than its inverse counterpart. A. Effect of the amplitude We performed the same simulation as shown on Figure 3, but with twice the initial wave amplitude a = 2 0.08. The results are demonstrated on Figure 4. One can see that the main characteristics of the double cascade remain unaffected, i.e. the number and position of cascading modes along with the exponential shape of the Fourier spectrum. This observation can be made even more precise as we plot the same fitted exponents on bottom panels. One can see on Figure 4(c) a very good agreement between the observed spikes with the observed above exponential decay rate. In contrast to the previous case with a = 0.08, all the processes are accelerated. Most importantly, the development of the MI in the physical space happens earlier with respect to the former computations shown on Figure 3. Consequently, the process enters in the FPU-like recurrence regime earlier showing that the cascade formation is intrinsically linked to the MI development. B. Effect of the perturbation We studied also the effect of the perturbation magnitude δ. We do not report the numerical results here due to the limitations in the Letter length. However, we report that for m = 0.05 0.5 (i.e. 5% to 50% of the base wave amplitude) the same processes described above still take place. Basically, the higher values of the parameter m yield the larger initial amount of the energy in the system which lead to the faster development of the MI. C. Effect of the spectral domain The last series of numerical experiments reported in this Letter aim at the studying the effect of the spectral domain on the energy cascade shape. With this purpose in mind we keep the same initial conditions as above, however we double the size of the spectral domain (consequently, the number of Fourier modes becomes N = 65 536 correspondingly). These results are shown on Figure 5. Our initial intuition, based on the experience with an isolated direct cascade, indicated that theincreaseofthespectraldomainwillleadtothechange of the spectrum exponents by preserving the global structure of the double cascade. Surprisingly, our simulations show that the cascade observed on Figure 3 preservesits shape and location in the Fourierspace(0 k < 120),whileintheadditionalspace (120 k 240) another (apparently slightly smaller) replica of the initial double cascade appears, shifted farther in the Fourier domain. Each replica is an equally stable structure in the Fourier space as the double cascade described above (i.e. with evolution the number and location of cascading modes remain the same along with the spectrum exponents). However, they are different from the isolated double cascade observed for precisely the same initial conditions. The main difference consists, for example, in the number of cascading modes which obviously affects the exponent of the energy spectrum shape. By comparing the development of the MI in the physical space on Figures 3(d,e) and 5(b,c) we can see that in the case of the smaller spectral domain the MI is more developed comparing to the latter. This phenomenon can be explained by the presence of two recurrent coherent structures (instead of one) each of them having a stabilizing effect on the MI. IV. SUMMARY Integrability of a specific evolutionary nonlinear PDE is a mathematical issue related to the spectrum of the corresponding Lax operator [9]. It gives a priori no information about the possible scenarios of the energy transport over the scales (in direct, inverse or both directions) as well as about the shape of the energy spectra. A good illustration is the nonlinear Schrödinger (NLS) equation which is integrable but this fact does not prevent numerous researchers from studying the corresponding energy spectra numerically, experimentally and analytically even
4 (a)t = 5 (b)t = 10 (c)t = 15 (d)t = 20 (e)t = 25 (f)t = 75 FIG. 3: Formation of the direct and inverse cascades in the mkdv equation. nowadays [10]. Let us recall that the complete integrability of the NLS equation was proved by Zakharov & Shabat in 1972 [11], while the observation of the inverse cascade in the NLS was done only 20 years later [12]. The energy transport in the mkdv, also known to be integrable by the inverse scattering transform, is even less studied than the NLS, and the existence of the inverseenergycascadein themkdvisanovelphenomenon which we believe should be studied further on, both numerically and experimentally. In this Letter we presented
5 some preliminary but very promising results on the formation of the inverse energy cascade in the mkdv equation. Briefly our findings can be summarized as: Contrary to the direct D-cascade which can be observed for any k 0 > 0, the inverse cascade starts to appear only for large values of k 0. It is not possible to observe an isolated inverse D-cascade, since the direct cascade forms for any choice of k 0. As a result, our observations show the so-called double cascade which is composed of at least one instance of the direct and inverse cascades. The exponents of the inverse and direct cascades are slightly different. In general, the direct cascade tends to decay faster. Moreover, the values of the exponents were shown not to depend on the initial wave amplitude, but only on the initial wavenumber k 0. The double cascade appears to be a very robust quasi-stationary structure in the Fourier space. In other words, the system at this stage enters into the FPU-like recurrent regime. The formation of the D-cascade is accompanied by the nonlinear development of the MI. The increase of the spectral domain yields the formation of two (and possibly even more) replica of the double cascade, each containing roughly one half of the initial energy. Our conjecture is that the further increase in the spectral domain will lead to the appearance of more and more replica of this double cascade structure. Or in other words, we will see the appearance of more and more recurrent structures sharing the total system energy. This process may be regarded as a mechanism of the MI stabilization. The problem of the MI stabilization was considered by several research groups [13] and the main stabilization mechanism which was put forward has been related to some sort of dissipation. In this study we propose another mechanism which can work in fully conservative systems. The present investigation is restricted to the case of moderate nonlinearities (a Ø ( 10 1) ), since in this case it is straightforward to find an unstable mode k 0. When we go to the field of the finite amplitude waves (a Ø(1)), the stability study becomes a complicated mathematical problem [14]. The investigation of the D- cascade formation in strongly nonlinear regimes, occuring for example in plasma physics [15], is among our next priorities. Acknowledgments This research has been supported by the Austrian Science Foundation (FWF) under projects P22943-N18 and P24671 as well as by the National Science Foundation under Grant No. NSF PHY11-25915. [1] A. Kolmogorov. Proc. R. Soc. Lond. A. 1991, 434, 9-13. [2] Zakharov V.E., Lvov V.S. and Falkovich G. Kolmogorov Spectra of Turbulence I (Springer, 1992) [3] Kartashova E. EPL97(2012): 30004. [4] Kartashova E. Phys. Rev. E86(2012): 041129. [5] Kartashova E. EPL102(2013): 44005. [6] Dutykh D. and Tobisch E. Submitted(2014), E-print: http://hal.archives-ouvertes.fr/hal-00990724/. [7] Grimshaw R., Pelinovsky D., Pelinovsky E. and Talipova T. Phys. D159(1-2)(2001): 35 57. [8] Grimshaw R., Pelinovsky E., Talipova T., Ruderman M. and Erdelyi R. Stud. Appl. Math.114(2)(2005): 189 210. [9] Ablowitz M.J. and Segur H. Solitons and the Inverse Scattering Transform (Society for Industrial & Applied Mathematics, 1981) [10] Kuksin S. and Maiocchi A. E-print: arxiv:1311.6794 [11] Zakharov V.E. and Shabat A.B. Sov. Phys. JETP34(1)(1972): 62 69 [12] Dyachenko S., Newell A., Pushkarev A. and Zakharov V.E. Phys. D57 (1992): 96-160 [13] Segur H., Henderson D., Carter J., Hammack J., Li C.- M., Pheiff D. and Socha K. J. Fluid Mech.539(2005): 229 771. [14] Johnson M., Zumbrun K. Bronski J. Phys D.239 (2010): 2057 2065. [15] Ruderman M., Talipova T. Pelinovsky E. J. Plasma Phys.74(5) (2008): 639 656.
6 (a)t = 5 (b)t = 10 (c)t = 30 FIG. 4: Development of the double cascade for the initial wave amplitude a = 2 0.08.
7 (a)t = 15 (b)t = 20 (c)t = 25 FIG. 5: Formation of two copies of the double inverse and direct cascades in a two times larger spectral domain. The exponents shown with dashed lines were measured in the simulations from Figure 3. The observed discrepancy shows the effect of the spectral domain on the energy spectrum shape.