Next: Intermittency and the multi-mode Up: Wave Turbulence "Recent developments Previous: Justification of RPA

# One-mode statistics

We have established above that the one-point statistics is at the heart of the WT theory. All one-point statistical objects can be derived from the one-point amplitude generating function,

which can be obtained from the -point by taking all 's and all 's, except for , equal to zero. Substituting such values to (64) (for the three-wave case) or to (72) (four-wave) we get the following equation for ,

 (74)

where for the three-wave case we have:
 (75) (76)

and in the four-wave case:
 (77) (78)

Here we introduced the wave-action spectrum,

 (79)

Differentiating (80) with respect to we get an equation for the moments :

 (80)

which, for gives the standard kinetic equation,

 (81)

First-order PDE (80) can be easily solved by the method of characteristics. Its steady state solution is

 (82)

which corresponds to the Gaussian values of momenta

 (83)

However, these solutions are invalid at small and high 's because large amplitudes , for which nonlinearity is not weak, strongly contribute in these cases. Because of the integral nature of definitions of and with respect to the , the ranges of amplitudes where WT is applicable are mixed with, and contaminated by, the regions where WT fails. Thus, to clearly separate these regions it is better to work with quantities which are local in , in particular the one-mode probability distribution . Equation for the one-mode PDF can be obtained by applying the inverse Laplace transform to (80). This gives:

 (84)

with is a probability flux in the s-space,

 (85)

Let us consider the steady state solutions, so that

 (86)

Note that in the steady state which follows from kinetic equation (87). The general solution to (92) is

 (87)

where

 (88)

is the general solution to the homogeneous equation (corresponding to ) and is a particular solution,

 (89)

where is the integral exponential function.

At the tail of the PDF, , the solution can be represented as series in ,

 (90)

Thus, the leading order asymptotics of the finite-flux solution is which decays much slower than the exponential (Rayleigh) part and, therefore, describes strong intermittency.

Note that if the weakly nonlinearity assumption was valid uniformly to then we had to put to ensure positivity of and the convergence of its normalisation, . In this case which is a pure Rayleigh distribution corresponding to the Gaussian wave field. However, WT approach fails for for the amplitudes for which the nonlinear time is of the same order or less than the linear wave period and, therefore, we can expect a cut-off of at . Estimate for the value of can be obtained from the dynamical equation (14) by balancing the linear and nonlinear terms and assuming that if the wave amplitude at some happened to be of the critical value then it will also be of similar value for a range of 's of width (i.e. the -modes are strongly correlated when the amplitude is close to critical). This gives:

 (91)

This cutoff can be viewed as a wavebreaking process which does not allow wave amplitudes to exceed their critical value, for . Now the normalisation condition can be satisfied for the finite-flux solutions. Note that in order for our analysis to give the correct description of the PDF tail, the nonlinearity must remain weak for the tail, which means that the breakdown happens far from the PDF core . When one has a strong breakdown predicted in [32] which is hard to describe rigorously due to strong nonlinearity.

Depending on the position in the wavenumber space, the flux can be either positive or negative. As we discussed above, results in an enhanced probability of large wave amplitudes with respect to the Gaussian fields. Positive mean depleted probability and correspond to the wavebreaking value which is closer to the PDF core. When gets into the core, one reaches the wavenumbers at which the breakdown is strong, i.e. of the kind considered in [32]. Consider for example the water surface gravity waves. Analysis of [32] predicts strong breakdown in the high- part of the energy cascade range. According to our picture, these high wavenumbers correspond to the highest positive values of and, therefore, the most depleted PDF tails with respect to the Gaussian distribution. When one moves away from this region toward lower 's, the value of gets smaller and, eventually, changes the sign leading to enhanced PDF tails at low 's. This picture is confirmed by the direct numerical simulations of the water surface equations reported in [23] the results of which are shown in figure 1.

 [width=.5]1.eps

This figure shows that at a high the PDF tail is depleted with respect to the Rayleigh distribution, whereas at a lower it is enhanced which corresponds to intermittency at this scale. Similar conclusion that the gravity wave turbulence is intermittent at low rather than high wavenumbers was reached on the basis of numerical simulations in [33]. To understand the flux reversal leading to intermittency appears in the one-mode statistics, one has to consider fluxes in the multi-mode phase space which will be done in the next section.

Next: Intermittency and the multi-mode Up: Wave Turbulence "Recent developments Previous: Justification of RPA
Dr Yuri V Lvov 2007-01-23