Abstract
We study the generalization of the dispersionless Kadomtsev  Petviashvili (dKP) equation in dimensions and with nonlinearity of degree , a model equation describing the propagation of weakly nonlinear, quasi one dimensional waves in the absence of dispersion and dissipation, and arising in several physical contexts, like acoustics, plasma physics, hydrodynamics and nonlinear optics. In dimensions and with quadratic nonlinearity, this equation is integrable through a novel IST, and it has been recently shown to be a prototype model equation in the description of the two dimensional wave breaking of localized initial data. In higher dimensions and with higher nonlinearity, the generalized dKP equations are not integrable, but their invariance under motions on the paraboloid allows one to construct in this paper a family of exact solutions describing waves constant on their paraboloidal wave front and breaking simultaneously in all points of it, developing after breaking either multivaluedness or single valued discontinuous shocks. Then such exact solutions are used to build the longtime behavior of the solutions of the Cauchy problem, for small and localized initial data, showing that wave breaking of small initial data takes place in the longtime regime if and only if . At last, the analytic aspects of such a wave breaking are investigated in detail in terms of the small initial data, in both cases in which the solution becomes multivalued after breaking or it develops a discontinuous shock. These results, contained in the 2012 master thesis of one of the authors (FS) [1], generalize those obtained in [2] for the dKP equation in dimensions with quadratic nonlinearity, and are obtained following the same strategy.
On the dispersionless KadomtsevPetviashvili equation with arbitrary nonlinearity and dimensionality:
exact solutions, longtime asymptotics of the Cauchy problem, wave breaking and discontinuous shocks
F. Santucci and P. M. Santini
Istituto Superiore IIS G.Marconi, via Reno snc, 04100 Latina, Italy
Dipartimento di Fisica, Università di Roma ”La Sapienza”, and
Istituto Nazionale di Fisica Nucleare, Sezione di Roma 1
Piazz.le Aldo Moro 2, I00185 Roma, Italy
email: nadir17@libero.it
email: paolo.santini@roma1.infn.it
July 12, 2019
1 Introduction
In this paper we investigate the following generalized dispersionless Kadomtsev  Petviashvili equation
(1.1) 
in dimensions, where is the degree of the nonlinearity, and . Hereafter, we refer to (1.1) as the equation. Equations (1.1) contain, as particular cases, the integrable (through the method of characteristics) Riemann equation [3]
(1.2) 
for , the integrable dispersionless Kadomtsev  Petviashvili (dKP) equation [4, 5, 6, 7]
(1.3) 
for , the nonintegrable Khokhlov  Zobolotskaya (KZ) equation [6]
(1.4) 
for , and the nonintegrable modified dKP (mdKP) equation [8]
(1.5) 
for .
Equation (1.1) describes weakly nonlinear and quasi onedimensional waves, with negligeable dissipation and dispersion, if the linear approximation of the original theory is given by the  dimensional wave equation, at least in some limiting cases like the long wave approximation [2]. In equation (1.1) there are two competing terms: the nonlinear term , responsible for the steepening of the profile, and the term, describing diffraction in the transversal dimensional hyperplane (reminiscence, through the multiscale expansion leading to (1.1), of the diffraction described by the wave operator); therefore diffraction increases as increases [2]. The degree of the nonlinearity originates from expanding the nonlinear terms of the original PDEs in power series, when the first powers of the expansion are absent due, usually, to some symmetry of the problem. In this sense, the degree of universality of (1.1) decreases as the degree of nonlinearity grows (see §2). are relevant in acoustics [4, 5, 6, 2], and also in plasma physics [7, 2, 1] in the long wave approximation; in hydrodynamics, again in the long wave approximation [9, 10]; and in nonlinear optics [1] and in the study of sound waves in antiferromagnets [11].
We remark that equation (1.1) arises as the dispersionless limit of the natural generalization in dimensions and with nonlinearity of degree :
(1.6) 
of the dimensional KadomtsevPetviashvili (KP) equation [7, 12, 9, 10].
Apart from the Riemann equation (1.2), integrable by the method of characteristics, in the family of equations (1.1) only the dKP equation (1.3) is integrable through a novel Inverse Scattering Transform (IST) for integrable dispersionless PDEs [13, 14, 15], recently made rigorous in [16] on the example of the Pavlov equation [17]. This IST allows one to show, in particular, that solutions of dKP depend on through the combination [18]; i.e., these solutions can be written in the characteristic form [18, 19]
(1.7) 
in analogy with the case of the Riemann equation (1.2), for which the dependence of the solution on is through the combination . For this reason, the IST for dKP can be viewed as a significant generalization of the method of characteristics. The formulation (1.7) has allowed one to study in an analytically explicit way the interesting features of the gradient catastrophe of two dimensional waves at finite time [18, 19] and in the longtime regime [18] in terms of the initial data.
The other examples of equations are not integrable; therefore the possibility to investigate a generic wave breaking through equations like (1.7) and the precise form that these equations should take are, in our opinion, challenging open problems; in addition, blow up of the solutions is expected for sufficiently large [8] to complicate the picture. In the recent paper [20], f.i., the formal dependence , motivated by the independent limit, was used to study the generic breaking features of and its dispersive shock formation.
In our paper we give some light on the problem of finding a convenient characteristic form of the type (1.7) for the study of wave breaking of solutions, i) from the construction of a family of exact solutions of exhibiting wave breaking, and ii) from the construction of the longtime behavior of solutions of the Cauchy problem of , for small initial data.
Indeed, after showing in §2 the universality of (1.5) starting from a family of nonlinear wave equations, in §3 we use the invariance of the equations under motions on the paraboloid, to construct a family of exact solutions involving an arbitrary function of one variable, and describing waves constant on their paraboloidal wave front, breaking simultaneously in all points of it, and developing, after breaking, either multivalued overturning profiles or single valued discontinuous shocks. Then we use in §4 such solutions to build a uniform approximation of the solution of the Cauchy problem, in the longtime regime and for small and localized initial data, showing that such a small and localized data evolving according to the equation break, in the longtime regime, iff . In §5 we study the analytic aspects of such a wave breaking, given explicitely in terms of the initial data, providing, in particular, a description of the overturning profile and of the development of a discontinuous shock immediately after breaking; we concentrate, in particular, on the mdKP (1.5) case and on its comparison with the already known dKP (1.3) case. The results of this paper, contained in the 2012 master thesis of one of the authors (FS) [1], generalize analogous ones for the equation in [2], and are obtained following the same strategy.
We end this introduction considering the following consequence of these results: for the family of exact solutions of and for the longtime asymptotics of the solutions of the small data Cauchy problem, the dependence of the solutions on is, f.i., through the combinations , if , and , if , where
(1.8) 
Therefore this dependence involves in a simple way the degree of the nonlinearity and the dimensionality of the problem through the coefficient , and it is consistent with the case of the Riemann equation (1.2), for which , and with the dKP equation (1.3), for which (confirming the apparently mysterious factor in (1.7)); it is not consistent, instead, with the formal dependence used in [20].
In analogy with the dKP case, the results of our paper seem to suggest the following conjecture (admittedly a weak conjecture, due to the absence, in the generic case, of an IST generalizing the method of characteristics).
Conjecture. The solutions of the equation around breaking, for , are described by the following characteristic formulae
(1.9) 
This paper is dedicated to the memory of S. V. Manakov.
2 Universality and applicability of
The universality (and therefore the applicability) of (1.1), through a multiscale expansion, is well described by the following simple model: the family of nonlinear wave equations
(2.1) 
For small amplitudes: , , we substitute in (2.1) the nonlinear term by its Taylor expansion
(2.2) 
obtaining, at , the wave equation
(2.3) 
where is assumed to be positive.
If the waves are quasi onedimensional and we choose as the direction of propagation, the wave lengths in the trasversal directions are small: , where is the transversal wave vector and has to be fixed. Then the dispersion relation becomes [2]
(2.4) 
and the phase of a monochromatic wave reads
(2.5) 
motivating the introduction of the new variables
(2.6) 
Rewriting (2.1) in the new variables and imposing to get the maximal balance, one obtains, at , the equation
(2.7) 
where , and is the tranversal Laplacian given in (1.1b).
If the term vanishes, the maximal balance must involve the cubic term; if also , the maximal balance must involve the quartic term, and so on. In the very special case in which , , and , then the maximal balance must involve the term of order , and is achieved for at the first nontrivial order , obtaining the equation (1.1), with and
(2.8) 
if is odd, and
(2.9) 
with , if is even.
The above considerations explain well why equations (1.1) are less and less universal and applicable, as increases.
3 Exact solutions of
It was observed in [2] that is invariant under the following Lie symmetry group of transformations
(3.1) 
where are the arbitrary parameters of the group, leaving invariant the paraboloid
(3.2) 
Such a symmetry was used to construct a family of exact and implicit solutions of exhibiting wave breaking and playing a relevant role in the longtime regime of the small data Cauchy problem for [2].
Since is also invariant under the transformation (3.1), following the same strategy as in [2], such a symmetry will be used in this section to construct a family of exact and implicit solutions of . In this paper we add, to the exact solutions exhibiting a gradient catastrophe and multivaluedness after breaking, also the exact weak solutions of developing, after breaking, single valued discontinuous shocks of dissipative nature.
Therefore we look for solutions of (1.1) in the form
(3.3) 
reducing (1.1) to the dimensional PDE
(3.4) 
The following change of variables
(3.5) 
where
(3.6) 
are real constants and the coefficient is defined in (1.8), transforms (3.4) into the Riemann equation (1.2) in the variables
(3.7) 
We recall [3] that (3.7) has the general implicit solution
(3.8) 
where is an arbitrary differentiable function of one argument. If, in particular, describes a localized positive hump, the solution breaks first at , on the characteristic , where
(3.9) 
For , the solution becomes multivalued and not acceptable in many physical contexts. Alternatively, for , the regular multivalued solution can be replaced by a weak solution, a single valued discontinuous shock of dissipative nature, whose wave front discontinuity is described by equations [3]
(3.10) 
with the initial conditions
(3.11) 
Behind and ahead the shock, the solution is given by
(3.12) 
where and are respectively the maximum and the minimum among the (three, in the case of a single hump) branches of the implicit equations (3.8).
We remark that in (3.6) is a monotonically increasing function of ; to have it positive, we choose the constants and the intervals as follows:
(3.13) 
Recalling (3.3), (3.5), and (3.13), the exact solutions of in the original variables corresponding to (3.8) read
(3.14) 
where is defined in (1.8).
If, again, is a localized hump, (3.14) describes a wave constant on the paraboloidal wave front (3.2) and simultaneously breaking on it. More precisely, the solution breaks at time defined by
(3.15) 
on the paraboloid
(3.16) 
where are defined in (3.9) (see Figures 1).
Figures 1. Plotting of the exact solution of the mdKP equation (1.5) (), before, at and after breaking. Here we have chosen ; consequently: . If , as , and the wave freezes asymptotically. If , then , and the wave breaks before freezing at
(3.17) 
on the paraboloid (3.16). If, instead, , then and no breaking takes place before the wave freezes.
If , (3.14) reduces to the general solution of the Riemann equation (1.2); if , (3.14) reduces to the class of particular solutions of the dKP equation with quadratic nonlinearity and arbitrary dimensions () constructed in [2]. It is interesting to remark that this type of exact solutions was first derived for the integrable dKP equation (1.3) using its nonlinear Riemnann  Hilbert inverse problem [21].
After breaking: , the regular multivalued solutions (3.14) can be replaced by the discontinuous single valued solutions of obtained using (3.3), (3.5), and (3.13) in (3.10), (3.11), (3.12). More precisely, for , where is defined in (3.15) and (3.17), the solution described by (3.14) for , develops a discontinuous shock on the parabola
(3.18) 
where is characterized by equations (3.10), (3.11), and is defined in (3.13). Behind and ahead the shock, the solution is given by
(3.19) 
where and are respectively the maximum and the minimum among the (three, in the case of a localized bump) branches of the implicit equations (3.14).
4 The Cauchy problem for small and localized initial data and longtime wave breaking
In analogy with [2], now we use the exact solutions of the previous section to construct the longtime behaviour of solutions of the Cauchy problem of (1.1) for small and localized initial data, showing in particular that small and localized initial data evolving according to (1.1) break, in the longtime regime, only if ().
For small and localized initial data of the type
(4.1) 
evolving according to equation (1.1), the behavior of the solution may be approximated by the solution of the linearized equation until nonlinearity becomes relevant.
For initial data, equations (3.7),(3.9) imply that the breaking time is
(4.2) 
so, the nonlinear regime for is characterized by the condition
(4.3) 
where is the inverse of (3.6), while, for , large values of are not compatible with (4.2), and a nonlinear regime never occurs.
So, for , the solution of the Cauchy problem is approximated by the solution of the linearized equation, while for , a matching will be made between the linearized solution and the exact solution (3.14).
In the linear regime, when , the solution is well described by that of the linearized , i.e.:
(4.4) 
where
In the longtime regime
(4.5) 
the multiple integral can be evaluated by the stationary phase method, which gives [2]
(4.6) 
where
(4.7) 
in the spacetime region
(4.8) 
on the paraboloid (3.2), and decays faster outside it. Therefore, small and localized initial data set, in the longtime regime (4.5) governed by the linear theory, on the paraboloid (3.2).
When nonlinearity becomes relevant, i.e., when , the asymptotic solution of is built matching the longtime solution (4.6),(4.7) of the linearized Cauchy problem with the exact solution (3.14). If, f.i., , the first argument of function defined in (4.7) is replaced by , or, equivalently, the arbitrary function in (3.8) acquires the dependence on the second argument , and is identified with . As a result of this matching, in the nonlinear regime (4.2), (4.3) the solution reads
(4.9) 
with for , and for , becoming infinitesimal in the linear regime, as it has to be.
For , the solution is described by its linear form (4.6); therefore wave breaking takes place for the particular values of described by the condition
(4.10) 
i.e., for

(4.11) 
For , reduces to the general solution of the Riemann equation (1.2) [3]. The cases and , the dKP and the KZ equations, have been already investigated in [18] and [2] respectively; therefore, in the following section we mainly focus on the case , corresponding to the nonintegrable modified dKP equation (1.5) and to its asymptotic solution
(4.12) 
comparing the results with the dKP case investigated in [18], with its asymptotic solution
(4.13) 
Remark 1. The estimate of the first correction to the asymptotics (4.9) reads
(4.14) 
Indeed, if , one has
(4.15) 
where and is expressed in terms of ( is defined in the second of equations (5.1)), through the PDE
(4.16) 
Since, from (4.3), , equation (4.15) yields (4.14). Similarly, if ,
(4.17) 
where now solves the PDE
(4.18) 
Then equation (4.17) yields equation (4.14) too. Remark 2. As is was observed in [2], for a gaussian initial condition
(4.19) 
where is constant, the asymptotic solution can be written in terms of special functions, obtaining for :
(4.20) 
where
(4.21) 
is the Euler gamma function and is the Kummer confluent hypergeometric function [22, 23]. For and , it yields, respectively:
(4.22) 
Remark 3. The fact that small data break only when the inequality (4.10) is satisfied was expected from physical considerations. Indeed, as already mentioned in the introduction, in equation (1.1) there are two competing terms: the nonlinear term , responsible for the steepening of the profile, and the term, describing diffraction in the transversal dimensional hyperplane (reminiscence, through the multiscale expansion leading to (1.1), of the wave operator); therefore diffraction increases as increases [2]. For small initial data, the term is initially smaller and smaller increasing , and the solution evolves in a linear way for a longer and longer time, diffracting transversally through the diffraction channels. So, diffraction, increasing with , acts on a very long time, increasing with , before the nonlinear regime could become relevant, and if and/or are sufficiently large, one expects that the solution would be diffracted away almost completely, before reaching the nonlinear regime, and will not break.
We have, in particular, that breaking takes place in the following longtime regimes (see (4.3)):
(the Riemann equation with quadratic nonlinearity): ;
(the dKP equation): ;
(the KZ equation): ;
(the Riemann equation with cubic nonlinearity): ;
(the mdKP equation): .
It is interesting to remark that, as it was observed in [2], if (the physically more relevant case), wave breaking takes place only for ; i.e., in physical space!
5 Analytic aspects of the wave breaking
In this section we use the longtime behavior of solutions of the small data Cauchy problem for to show explicitly the analytic aspects of the solution in the neighbourhood of the breaking time in terms of the initial data, represented by function G defined in (4.7), as it was already done in [18, 2].
One solves the second of equations (5.1) with respect to the parameter , obtaining and replaces it into the first, to obtain the solution . The inversion of the second of equations (5.1) is possible iff its derivative is different from zero. Therefore the singularity manifold (SM) of the solution is the  dimensional manifold characterized by the equation
(5.4) 
Since
(5.5) 
the slope of the localized wave becomes infinity (the socalled gradient catastrophe) on the SM, and the dimensional wave “breaks” except on the direction represented by the vector field [2], [19], on which:
(5.6) 
Since the KZ case has been studied in detail in [2], in the following we mainly concentrate on the mdKP case , and we compare it with the dKP case investigated in [18].
From (5.4), the first breaking time and the corresponding characteristic parameters are defined by the global minimum of a function of the variables :