MULTISCALED SOLITARY WAVES

It is analytically shown how competing nonlinearities yield multiscaled structures for internal solitary waves in stratified shallow fluids. These solitary waves only exist for large amplitudes beyond the limit of applicability of the KdV equation or its usual extensions. Multiscaling phenomenon exists or does not exist for almost identical density profiles. Trapped core inside the wave prevents appearance of such multiple scales within the core area. The structural stability of waves of large amplitude is briefly discussed. Waves of large amplitude displaying both quadratic, cubic and higher order nonlinear terms have 5 stable and unstable branches. Multiscaled waves without vortex core are shown to be structurally unstable. It is anticipated that multiscaling phenomena exist for solitary waves in various physical contexts.


Introduction
The typical horizontal scale (or scales) is a major characteristic of a plane disturbance propagating in a nonuniform medium.
Usually, in an ideal density stratified shallow fluid, a wave of small, albeit finite amplitude has one typical scale resulting from the (local) balance between nonlinearity and dispersion like in the realm of Korteveg-de Vries (KdV) equation (Helfrich and Melville , 2006).Solitary waves of permanent form for which capillary dispersion is of the same order as the gravitational one may have oscillatory outskirts as predicted by Benjamin (1992).When viscosity is taken into account, transient effects leading to various length scales are discussed for the KdV type equation with cubic nonlinearity, for example by Grimshaw et al. (2003).In the present note it is shown that for the gravitational dispersion, ignoring all other earlier mentioned effects, solitary waves with multiple scales are possible.These solutions exist only for disturbances of finite amplitude exceeding the range of applicability of the extended KdV model, which incorporates both quadratic and cubic nonlinearities.Higher nonlinearity in the existing small amplitude KdV or mKdV models leads to the correction of the wave length scale without generation of multiscaling.For appearance of multiscaling, the various competitive nonlinearities should be of the same order and that order needs to be higher than the cubic one as analytically discussed below.This effect was initially noticed by Derzho and Borisov (1990) in the Russian journal but the result was not widely disseminated.Recently Dunphy et al. (2011) presented numerical procedure that provides fast calculations for gravitational waves between rigid lids .This model is able to work with fine density stratifications.Dunphy et al. (2011) reported two-humped and usual one-humped solitary internal waves solutions for nearly identical density profiles in a two-pycnocline density stratification.Lamb and Wan (1998) have numerically shown that in some stratifications with two pycnoclines three conjugate flow solutions leading to two-humped solitary waves were present.Makarenko et al. (2009) theoretically considered continuous stratification in order to characterize the role of vertical structure of the fluid density in the context of waves close to the limiting amplitude.To the best of author knowledge, neither specific nonlinearity in terms of power series of wave amplitude necessary to reveal a two-humped structure nor regions of density profiles with single pycnocline at which such structures exist were not examined in the literature.Kurkina et al. (2011) derived KdV like equation with quadratic and quartic nonlinear terms for interfacial transient waves for the specific three layer geometry.Assumption on small albeit finite wave amplitude was essential to balance nonlinearity and dispersion in that study.
Table-top limiting solutions were reported and they were stable within the accuracy of their numerical scheme.In the current paper, an asymptotic model presented in earlier papers by the author addresses multiscaling phenomena for internal solitary waves under free surface in the frame of Dubreil-Jacotin-Long equation (DJL) (Long , 1965).Special attention is given to the case of complicated nonlinearity involving both quadratic, cubic and quartic nonlinear terms for the case of continuous stratification with a single pycnocline .Solitary waves of permanent form, their existence and structural stability are discussed.
It is worth noting that family of solutions is richer than two-humped structures.It is expectable that such multiscaled solitary waves exist in other physical systems where complicated competitive nonlinearities are balanced by dispersion.

Model for internal waves
Let us consider the two-dimensional steady motion of an ideal stratified fluid in a frame of reference moving with phase speed of wave c .The approach is asymptotic being based on the DJL equation for waves without a priori limitation on amplitude.This approach has started from the pioneering work by Benney and Ko (1978).Let us consider the stratification in the form where σ denotes Boussinesq parameter.In (Derzho and Velarde, 1995) it was shown that for this case the dimensionless (primed) streamfunction ψ = −ψ/cH of a solitary disturbance obeys the equation where µ is the aspect ratio H/L and λ = σgH c 2 .In Eq.( 2) z denotes the vertical axis, taken positive upwards and x corresponds to the horizontal axis; z and x are scaled with H and L, the given vertical and horizontal scales respectively.Expecting no confusion we have, for simplicity, dropped the primes in Eq.( 2).Let us locate the bottom and the surface at the dimensionless heights z = −0.5 and z = 0.5 + η(x), respectively, where η(x) denotes surface displacement.The boundary conditions at the bottom and surface are The solution of Eqs.(2-5) is sought in the form where zeroth order variables are of order unity.Below we shall provide solution for the first mode, which is most frequently observed in nature.The analysis for the higher modes is similar.In the zeroth order where the amplitude function A(x) is to be determined at a higher order.For the solution to the first order equation to exist the solvability condition (Fredholm alternative) demands In order to (locally) balance nonlinearity and dispersion we have to require max(σ/µ 2 , δ/µ 2 ) ∼ 1 thus determining L. Benney and Ko (1978) suggested to consider the nonlinear terms as power series in the Boussinesq parameter instead of the small amplitude parameter.Derzho and Velarde (1995) somewhat extended this idea to account a more general undisturbed flow state.After straightforward integrations for remaining solitary wave solution, Eqs.(8-9) can be reduced to The Weierstrass approximation theorem states that every continuous function defined on a closed interval can be uniformly approximated as closely as desired by a polynomial function.Recent account on the topic is reviewed in Hazewinkel (2001).
Thus, the integral below can be represented with the help of some N th-order polynomial according to the Weierstrass approximation theorem.In the current study only polynomial formula for stratification is considered, thus it directly leads to nonlinearities in the polynomial form.
For the wave of amplitude A 0 Eq.(10) yields (13) Equations (12-14) determine completely both profile and phase velocity of a solitary wave with amplitude A 0 .

Multiscaling
The function f in the form of a M th-order polynomial generates P N with the index N = M − 1.The power index of Φ is thus max(1, M − 2).The condition for Eq.( 8) to possess a multiscaled solution reduces to the condition that Φ(A, A 0 ) must be sign-defined with several extrema within [0, A 0 ] .Thus it must have more than two imaginary roots on that interval.It determines M ≥ 4 , i.e. for a stratification in the form of cubic polynomial or if wave amplitude is small enough to neglect A 4 and higher order nonlinearities, multiscaled solitary waves do not exist because f has no imaginary roots for this case.This is why classical KdV or mKdV can not provide multiscaled solitary waves over flat bottom.Let us consider wave structures for the density stratification in the form, which produces quadratic, cubic and quartic terms in Eq.( 8).Thus Eq. ( 12) for this case of stratification becomes Two humped solitary waves for the stratification given by Eq.( 17) exist in the domain shown in Fig. 1.
Indeed, the maximum derivative on x in the dimensionless coordinates is of order unity.However, the wave has a pronounced two-scale structure with typical length scales, which are much larger than the length L used to scale the derivative.A solitary wave with three typical length scales ( three-humped one) is shown in Fig. 4.
For this case the stratification profile is which produces in Eq. ( 8) nonlinear terms up to A 6 .Generally, one can expect at most M/2 different scales for a stratification in the form of polynomial with even power index M , and (M − 1)/2 otherwise.Further, we wish to examine the structure of solitary waves of permanent form for the stratification given by Eq.( 16).We only consider the case α = −1.39 and focus on the waves of permanent form under free surface, their domain of existence, 5 Fig. 5 shows that the phase velocity is an increasing function for 0 < A 0 < A 2 and A 0 > A 1 .For A 1 < A 0 < A 2 the phase velocity decreases with amplitude and there are no steady solitary wave solutions.When 0 < A 0 < A 2 solitary waves are  widening as amplitude increases with a table top limiting shape with a local maximum for the wave velocity as shown in Fig. 5 and Fig. 6.
Such waves are structurally stable according to Bona et al. (1987) as both the wave energy E = ∞ −∞ A 2 dx and the wave velocity increase as amplitude increases.

5
For A 0 > A 1 wave profiles are shown in Fig. 7. Waves change from the table top solution to solitary waves with a single scale via multiscaled structures.For the particular stratification considered here waves are structurally unstable (Bona et al., 1987) since wave energy decreases as shown in Fig. 8 but the wave velocity increases with the increasing of wave amplitude.Interesting observation is these waves of sufficiently large amplitude could be stable as the energy is eventually increased as shown in Fig. 8.For the stratification considered here, it does not matter because solution with vortex core appears at lower amplitude when energy is still the decreasing function of amplitude.However it leads to the interesting phenomenon -waves with vortex core could stabilise the wave.The idea is that the vortex core leads to widening of wave (Derzho and Grimshaw , 1997) and consequently to the increase of its energy, thus the structural stability criterion will be satisfied.For the considered particular stratification waves with vortex core are initially unstable as increase of energy due to the vortex core and associated widening does not compensate the decrease of energy in the wave outside the vortex core.Nonetheless, above some amplitude waves become structurally stable.When wave amplitude further increases the permanent wave of limiting amplitude becomes infinitely wide as shown by Derzho and Grimshaw (1997).
The theory described above is valid for wave amplitudes below A * , a certain amplitude at which a vortex core started to appear inside the wave.For nearly linear density profile A * = 1/π.Derzho and Grimshaw (1997) have shown that where ν is the supercriticality parameter defined such that B varies from zero to one as wave amplitude does from A * to the maximum value allowed predicted there.R(A * ) depends on stratification profile and is fixed.It is straightforward to notice that B(x) is monotonic and therefore multiscaling in the vortex core area does not exist when A > A * .
Multiscaling effects similar to the discussed above, could be observed in various physical media.Derzho and Grimshaw (2005) reported that solitary Rossby waves in channels obey the same KdV type equation with complicated nonlinearity due to the mean shear variations.Coriolis force for Rossby waves plays the same role as gravitational force for the internal gravity waves.The results on multiscaling for Rossby waves with and without trapped core will be reported elsewhere.

Conclusion
For the particular case of a nonlinear dispersive medium like a shallow stratified fluid embedded in the gravity field, we have addressed multiscaled solitary waves which are predicted when there exists competition of several different types of nonlinearity.The mechanism leading to these solutions differs from the mechanism of multiscaling due to the competition of different types of dispersion or effects due to the dissipation.We have shown that the length used to scale the x-derivative does not simply coincide with the typical length scale of the wave, as for KdV.Moreover, multiscaled (multi humped) disturbances exist for sufficiently large amplitudes, at least terms in fourth order of waves amplitude should be accounted.Multiscaling (multi humped) phenomenon exists or does not exist for almost identical density profiles, two pycnoclines case studied earlier is not necessary for the existence of multiscaling.Continuous stratification given by Eq.( 16) was studied in more detail.The structure of permanent solitary waves and how multiscaling appeared were presented.Structural stability was examined using the criterion proposed by Bona et al. (1987).It was shown that both stable and unstable solutions of the KdV type equation with quadratic, cubic and quartic nonlinearities are available.Multiscaled waves without trapped core belong to the unstable solutions.Trapped core inside the wave prevents appearance of such multiple scales within the core area.However, trapped core could stabilise the multiscaled solution in the sense of structural stability.The case when trapped core and multiscaling are combined together is beyond the scope of the present study and will be presented elsewhere.It is noted that multiscaling phenomena exist for solitary waves in various physical origins, for example, for Rossby waves on a shear flow (Derzho and Grimshaw, 2005) or inertial waves in swirling flows (Derzho and Grimshaw, 2002).