A detailed numerical study of the long time behaviour of dispersive shock waves in solutions to the Kadomtsev–Petviashvili (KP) I equation is presented. It is shown that modulated lump solutions emerge from the dispersive shock waves. For the description of dispersive shock waves, Whitham modulation equations for KP are obtained. It is shown that the modulation equations near the soliton line are hyperbolic for the KPII equation while they are elliptic for the KPI equation leading to a focusing effect and the formation of lumps. Such a behaviour is similar to the appearance of breathers for the focusing nonlinear Schrödinger equation in the semiclassical limit.

VL - 474 UR - https://royalsocietypublishing.org/doi/abs/10.1098/rspa.2017.0458 ER - TY - JOUR T1 - Painlevé IV Critical Asymptotics for Orthogonal Polynomials in the Complex Plane JF - Symmetry, Integrability and Geometry. Methods and Applications Y1 - 2018 A1 - Marco Bertola A1 - José Gustavo Elias Rebelo A1 - Tamara Grava AB -We study the asymptotic behaviour of orthogonal polynomials in the complex plane that are associated to a certain normal matrix model. The model depends on a parameter and the asymptotic distribution of the eigenvalues undergoes a transition for a special value of the parameter, where it develops a corner-type singularity. In the double scaling limit near the transition we determine the asymptotic behaviour of the orthogonal polynomials in terms of a solution of the Painlev´e IV equation. We determine the Fredholm determinant associated to such solution and we compute it numerically on the real line, showing also that the corresponding Painlev´e transcendent is pole-free on a semiaxis.

PB - National Academy of Sciences of Ukraine VL - 14 ER - TY - RPRT T1 - On critical behaviour in systems of Hamiltonian partial differential equations Y1 - 2013 A1 - Boris Dubrovin A1 - Tamara Grava A1 - Christian Klein A1 - Antonio Moro AB -We study the critical behaviour of solutions to weakly dispersive Hamiltonian systems considered as perturbations of elliptic and hyperbolic systems of hydrodynamic type with two components. We argue that near the critical point of gradient catastrophe of the dispersionless system, the solutions to a suitable initial value problem for the perturbed equations are approximately described by particular solutions to the Painlev\'e-I (P$_I$) equation or its fourth order analogue P$_I^2$. As concrete examples we discuss nonlinear Schr\"odinger equations in the semiclassical limit. A numerical study of these cases provides strong evidence in support of the conjecture.

PB - SISSA U1 - 7280 U2 - Mathematics U4 - 1 U5 - MAT/07 FISICA MATEMATICA ER - TY - RPRT T1 - On the tritronquée solutions of P$_I^2$ Y1 - 2013 A1 - Tamara Grava A1 - Andrey Kapaev A1 - Christian Klein AB -For equation P$_I^2$, the second member in the P$_I$ hierarchy, we prove existence of various degenerate solutions depending on the complex parameter $t$ and evaluate the asymptotics in the complex $x$ plane for $|x|\to\infty$ and $t=o(x^{2/3})$. Using this result, we identify the most degenerate solutions $u^{(m)}(x,t)$, $\hat u^{(m)}(x,t)$, $m=0,\dots,6$, called {\em tritronqu\'ee}, describe the quasi-linear Stokes phenomenon and find the large $n$ asymptotics of the coefficients in a formal expansion of these solutions. We supplement our findings by a numerical study of the tritronqu\'ee solutions.

PB - SISSA U1 - 7282 U2 - Mathematics U4 - 1 U5 - MAT/07 FISICA MATEMATICA ER - TY - JOUR T1 - The KdV hierarchy: universality and a Painleve transcendent JF - International Mathematics Research Notices, vol. 22 (2012) , page 5063-5099 Y1 - 2012 A1 - Tom Claeys A1 - Tamara Grava KW - Small-Dispersion limit AB - We study the Cauchy problem for the Korteweg-de Vries (KdV) hierarchy in the small dispersion limit where $\e\to 0$. For negative analytic initial data with a single negative hump, we prove that for small times, the solution is approximated by the solution to the hyperbolic transport equation which corresponds to $\e=0$. Near the time of gradient catastrophe for the transport equation, we show that the solution to the KdV hierarchy is approximated by a particular Painlev\'e transcendent. This supports Dubrovins universality conjecture concerning the critical behavior of Hamiltonian perturbations of hyperbolic equations. We use the Riemann-Hilbert approach to prove our results. PB - Oxford University Press UR - http://hdl.handle.net/1963/6921 N1 - This article was published in "International Mathematics Research Notices, vol. 22 (2012) , page 5063-5099 U1 - 6902 U2 - Mathematics U4 - 1 ER - TY - JOUR T1 - Numerical study of the small dispersion limit of the Korteweg-de Vries equation and asymptotic solutions JF - Physica D 241, nr. 23-24 (2012): 2246-2264 Y1 - 2012 A1 - Tamara Grava A1 - Christian Klein KW - Korteweg-de Vries equation AB - We study numerically the small dispersion limit for the Korteweg-de Vries (KdV) equation $u_t+6uu_x+\epsilon^{2}u_{xxx}=0$ for $\epsilon\ll1$ and give a quantitative comparison of the numerical solution with various asymptotic formulae for small $\epsilon$ in the whole $(x,t)$-plane. The matching of the asymptotic solutions is studied numerically. PB - Elsevier U1 - 7069 U2 - Physics U4 - -1 ER - TY - JOUR T1 - Numerical Study of breakup in generalized Korteweg-de Vries and Kawahara equations JF - SIAM J. Appl. Math. 71 (2011) 983-1008 Y1 - 2011 A1 - Boris Dubrovin A1 - Tamara Grava A1 - Christian Klein AB - This article is concerned with a conjecture in [B. Dubrovin, Comm. Math. Phys., 267 (2006), pp. 117–139] on the formation of dispersive shocks in a class of Hamiltonian dispersive regularizations of the quasi-linear transport equation. The regularizations are characterized by two arbitrary functions of one variable, where the condition of integrability implies that one of these functions must not vanish. It is shown numerically for a large class of equations that the local behavior of their solution near the point of gradient catastrophe for the transport equation is described by a special solution of a Painlevé-type equation. This local description holds also for solutions to equations where blowup can occur in finite time. Furthermore, it is shown that a solution of the dispersive equations away from the point of gradient catastrophe is approximated by a solution of the transport equation with the same initial data, modulo terms of order $\\\\epsilon^2$, where $\\\\epsilon^2$ is the small dispersion parameter. Corrections up to order $\\\\epsilon^4$ are obtained and tested numerically. PB - SIAM UR - http://hdl.handle.net/1963/4951 U1 - 4732 U2 - Mathematics U3 - Mathematical Physics U4 - -1 ER - TY - RPRT T1 - Numerical Solution of the Small Dispersion Limit of the Camassa-Holm and Whitham Equations and Multiscale Expansions Y1 - 2010 A1 - Simonetta Abenda A1 - Tamara Grava A1 - Christian Klein AB - The small dispersion limit of solutions to the Camassa-Holm (CH) equation is characterized by the appearance of a zone of rapid modulated oscillations. An asymptotic description of these oscillations is given, for short times, by the one-phase solution to the CH equation, where the branch points of the corresponding elliptic curve depend on the physical coordinates via the Whitham equations. We present a conjecture for the phase of the asymptotic solution. A numerical study of this limit for smooth hump-like initial data provides strong evidence for the validity of this conjecture.... UR - http://hdl.handle.net/1963/3840 U1 - 487 U2 - Mathematics U3 - Mathematical Physics ER - TY - JOUR T1 - Painlevé II asymptotics near the leading edge of the oscillatory zone for the Korteweg-de Vries equation in the small-dispersion limit JF - Comm. Pure Appl. Math. 63 (2010) 203-232 Y1 - 2010 A1 - Tom Claeys A1 - Tamara Grava AB - In the small dispersion limit, solutions to the Korteweg-de Vries equation develop an interval of fast oscillations after a certain time. We obtain a universal asymptotic expansion for the Korteweg-de Vries solution near the leading edge of the oscillatory zone up to second order corrections. This expansion involves the Hastings-McLeod solution of the Painlev\\\\\\\'e II equation. We prove our results using the Riemann-Hilbert approach. PB - Wiley UR - http://hdl.handle.net/1963/3799 U1 - 527 U2 - Mathematics U3 - Mathematical Physics ER - TY - JOUR T1 - Solitonic asymptotics for the Korteweg-de Vries equation in the small dispersion limit JF - SIAM J. Math. Anal. 42 (2010) 2132-2154 Y1 - 2010 A1 - Tamara Grava A1 - Tom Claeys AB - We study the small dispersion limit for the Korteweg-de Vries (KdV) equation $u_t+6uu_x+\\\\epsilon^{2}u_{xxx}=0$ in a critical scaling regime where $x$ approaches the trailing edge of the region where the KdV solution shows oscillatory behavior. Using the Riemann-Hilbert approach, we obtain an asymptotic expansion for the KdV solution in a double scaling limit, which shows that the oscillations degenerate to sharp pulses near the trailing edge. Locally those pulses resemble soliton solutions of the KdV equation. UR - http://hdl.handle.net/1963/3839 U1 - 488 U2 - Mathematics U3 - Mathematical Physics ER - TY - JOUR T1 - Initial value problem of the Whitham equations for the Camassa-Holm equation JF - Physica D 238 (2009) 55-66 Y1 - 2009 A1 - Tamara Grava A1 - Virgil U. Pierce A1 - Fei-Ran Tian AB - We study the Whitham equations for the Camassa-Holm equation. The equations are neither strictly hyperbolic nor genuinely nonlinear. We are interested in the initial value problem of the Whitham equations. When the initial values are given by a step function, the Whitham solution is self-similar. When the initial values are given by a smooth function, the Whitham solution exists within a cusp in the x-t plane. On the boundary of the cusp, the Whitham equation matches the Burgers solution, which exists outside the cusp. PB - Elsevier UR - http://hdl.handle.net/1963/3429 U1 - 906 U2 - Mathematics U3 - Mathematical Physics ER - TY - JOUR T1 - On universality of critical behaviour in the focusing nonlinear Schrödinger equation, elliptic umbilic catastrophe and the \\\\it tritronquée solution to the Painlevé-I equation JF - J. Nonlinear Sci. 19 (2009) 57-94 Y1 - 2009 A1 - Boris Dubrovin A1 - Tamara Grava A1 - Christian Klein AB - We argue that the critical behaviour near the point of ``gradient catastrophe\\\" of the solution to the Cauchy problem for the focusing nonlinear Schr\\\\\\\"odinger equation $ i\\\\epsilon \\\\psi_t +\\\\frac{\\\\epsilon^2}2\\\\psi_{xx}+ |\\\\psi|^2 \\\\psi =0$ with analytic initial data of the form $\\\\psi(x,0;\\\\epsilon) =A(x) e^{\\\\frac{i}{\\\\epsilon} S(x)}$ is approximately described by a particular solution to the Painlev\\\\\\\'e-I equation. UR - http://hdl.handle.net/1963/2525 U1 - 1593 U2 - Mathematics U3 - Mathematical Physics ER - TY - JOUR T1 - Universality of the break-up profile for the KdV equation in the small dispersion limit using the Riemann-Hilbert approach JF - Comm. Math. Phys. 286 (2009) 979-1009 Y1 - 2009 A1 - Tamara Grava A1 - Tom Claeys AB - We obtain an asymptotic expansion for the solution of the Cauchy problem for the Korteweg-de Vries (KdV) equation in the small dispersion limit near the point of gradient catastrophe (x_c,t_c) for the solution of the dispersionless equation.\\nThe sub-leading term in this expansion is described by the smooth solution of a fourth order ODE, which is a higher order analogue to the Painleve I equation. This is in accordance with a conjecture of Dubrovin, suggesting that this is a universal phenomenon for any Hamiltonian perturbation of a hyperbolic equation. Using the Deift/Zhou steepest descent method applied on the Riemann-Hilbert problem for the KdV equation, we are able to prove the asymptotic expansion rigorously in a double scaling limit. UR - http://hdl.handle.net/1963/2636 U1 - 1487 U2 - Mathematics U3 - Mathematical Physics ER - TY - JOUR T1 - Numerical study of a multiscale expansion of the Korteweg-de Vries equation and Painlevé-II equation JF - Proc. R. Soc. A 464 (2008) 733-757 Y1 - 2008 A1 - Tamara Grava A1 - Christian Klein AB - The Cauchy problem for the Korteweg de Vries (KdV) equation with small dispersion of order $\\\\e^2$, $\\\\e\\\\ll 1$, is characterized by the appearance of a zone of rapid modulated oscillations. These oscillations are approximately described by the elliptic solution of KdV where the amplitude, wave-number and frequency are not constant but evolve according to the Whitham equations. Whereas the difference between the KdV and the asymptotic solution decreases as $\\\\epsilon$ in the interior of the Whitham oscillatory zone, it is known to be only of order $\\\\epsilon^{1/3}$ near the leading edge of this zone. To obtain a more accurate description near the leading edge of the oscillatory zone we present a multiscale expansion of the solution of KdV in terms of the Hastings-McLeod solution of the Painlev\\\\\\\'e-II equation. We show numerically that the resulting multiscale solution approximates the KdV solution, in the small dispersion limit, to the order $\\\\epsilon^{2/3}$. UR - http://hdl.handle.net/1963/2592 U1 - 1530 U2 - Mathematics U3 - Mathematical Physics ER - TY - RPRT T1 - Numerical solution of the small dispersion limit of Korteweg de Vries and Whitham equations Y1 - 2007 A1 - Tamara Grava A1 - Christian Klein AB - The Cauchy problem for the Korteweg de Vries (KdV) equation with small dispersion of order $\\\\epsilon^2$, is characterized by the appearance of a zone of rapid modulated oscillations of wave-length of order $\\\\epsilon$. These oscillations are approximately described by the elliptic solution of KdV where the amplitude, wave-number and frequency are not constant but evolve according to the Whitham equations. In this manuscript we give a quantitative analysis of the discrepancy between the numerical solution of the KdV equation in the small dispersion limit and the corresponding approximate solution for values of $\\\\epsilon$ between $10^{-1}$ and $10^{-3}$. The numerical results are compatible with a difference of order $\\\\epsilon$ within the `interior\\\' of the Whitham oscillatory zone, of order $\\\\epsilon^{1/3}$ at the left boundary outside the Whitham zone and of order $\\\\epsilon^{1/2}$ at the right boundary outside the Whitham zone. JF - Comm. Pure Appl. Math. 60 (2007) 1623-1664 UR - http://hdl.handle.net/1963/1788 U1 - 2756 U2 - Mathematics U3 - Mathematical Physics ER - TY - RPRT T1 - Numerical study of a multiscale expansion of KdV and Camassa-Holm equation Y1 - 2007 A1 - Tamara Grava A1 - Christian Klein AB - We study numerically solutions to the Korteweg-de Vries and Camassa-Holm equation close to the breakup of the corresponding solution to the dispersionless equation. The solutions are compared with the properly rescaled numerical solution to a fourth order ordinary differential equation, the second member of the Painlev\\\\\\\'e I hierarchy. It is shown that this solution gives a valid asymptotic description of the solutions close to breakup. We present a detailed analysis of the situation and compare the Korteweg-de Vries solution quantitatively with asymptotic solutions obtained via the solution of the Hopf and the Whitham equations. We give a qualitative analysis for the Camassa-Holm equation UR - http://hdl.handle.net/1963/2527 U1 - 1591 U2 - Mathematics U3 - Mathematical Physics ER - TY - JOUR T1 - Reciprocal transformations and flat metrics on Hurwitz spaces JF - J. Phys. A 40 (2007) 10769-10790 Y1 - 2007 A1 - Simonetta Abenda A1 - Tamara Grava AB - We consider hydrodynamic systems which possess a local Hamiltonian structure of Dubrovin-Novikov type. To such a system there are also associated an infinite number of nonlocal Hamiltonian structures. We give necessary and sufficient conditions so that, after a nonlinear transformation of the independent variables, the reciprocal system still possesses a local Hamiltonian structure of Dubrovin-Novikov type. We show that, under our hypotheses, bi-hamiltonicity is preserved by the reciprocal transformation. Finally we apply such results to reciprocal systems of genus g Whitham-KdV modulation equations. UR - http://hdl.handle.net/1963/2210 U1 - 2034 U2 - Mathematics U3 - Mathematical Physics ER - TY - RPRT T1 - Large Parameter Behavior of Equilibrium Measures Y1 - 2006 A1 - Tamara Grava A1 - Fei-Ran Tian AB - We study the equilibrium measure for a logarithmic potential in the presence of an external field V*(x) + tp(x), where t is a parameter, V*(x) is a smooth function and p(x) a monic polynomial. When p(x) is of an odd degree, the equilibrium measure is shown to be supported on a single interval as |t| is sufficiently large. When p(x) is of an even degree, the equilibrium measure is supported on two disjoint intervals as t is negatively large; it is supported on a single interval for convex p(x) as t is positively large and is likely to be supported on multiple disjoint intervals for non-convex p(x). JF - Commun. Math. Sci. 4 (2006) 551-573 UR - http://hdl.handle.net/1963/1789 U1 - 2755 U2 - Mathematics U3 - Mathematical Physics ER - TY - JOUR T1 - Thomae type formulae for singular Z_N curves JF - Lett. Math. Phys. 76 (2006) 187-214 Y1 - 2006 A1 - Victor Z. Enolski A1 - Tamara Grava AB - We give an elementary and rigorous proof of the Thomae type formula for singular $Z_N$ curves. To derive the Thomae formula we use the traditional variational method which goes back to Riemann, Thomae and Fuchs. An important step of the proof is the use of the Szego kernel computed explicitly in algebraic form for non-singular 1/N-periods. The proof inherits principal points of Nakayashiki\\\'s proof [31], obtained for non-singular ZN curves. UR - http://hdl.handle.net/1963/2125 U1 - 2118 U2 - Mathematics U3 - Mathematical Physics ER - TY - JOUR T1 - Modulation of the Camassa-Holm equation and reciprocal transformations JF - Ann. Inst. Fourier (Grenoble) 55 (2005) 1803-1834 Y1 - 2005 A1 - Simonetta Abenda A1 - Tamara Grava AB - We derive the modulation equations or Whitham equations for the Camassa-Holm (CH) equation. We show that the modulation equations are hyperbolic and admit bi-Hamiltonian structure. Furthermore they are connected by a reciprocal transformation to the modulation equations of the first negative flow of the Korteweg de Vries (KdV) equation. The reciprocal transformation is generated by the Casimir of the second Poisson bracket of the KdV averaged flow. We show that the geometry of the bi-Hamiltonian structure of the KdV and CH modulation equations is quite different: indeed the KdV averaged bi-Hamiltonian structure can always be related to a semisimple Frobenius manifold while the CH one cannot. UR - http://hdl.handle.net/1963/2305 U1 - 1711 U2 - Mathematics U3 - Mathematical Physics ER - TY - JOUR T1 - Singular Z_N curves, Riemann-Hilbert problem and modular solutions of the Schlesinger equation JF - Int. Math. Res. Not. 2004, no. 32, 1619-1683 Y1 - 2004 A1 - Victor Z. Enolski A1 - Tamara Grava AB - We are solving the classical Riemann-Hilbert problem of rank N>1 on the extended complex plane punctured in 2m+2 points, for NxN quasi-permutation monodromy matrices. Following Korotkin we solve the Riemann-Hilbert problem in terms of the Szego kernel of certain Riemann surfaces branched over the given 2m+2 points. These Riemann surfaces are constructed from a permutation representation of the symmetric group S_N to which the quasi-permutation monodromy representation has been reduced. The permutation representation of our problem generates the cyclic subgroup Z_N. For this reason the corresponding Riemann surfaces of genus N(m-1) have Z_N symmetry. This fact enables us to write the matrix entries of the solution of the NxN Riemann-Hilbert problem as a product of an algebraic function and theta-function quotients. The algebraic function turns out to be related to the Szego kernel with zero characteristics. From the solution of the Riemann- Hilbert problem we automatically obtain a particular solution of the Schlesinger system. The tau-function of the Schlesinger system is computed explicitly. The rank 3 problem with four singular points (0,t,1,\\\\infty) is studied in detail. The corresponding solution of the Riemann-Hilbert problem and the Schlesinger system is given in terms of Jacobi\\\'s theta-function with modulus T=T(t), Im(T)>0. The function T=T(t) is invertible if it belongs to the Siegel upper half space modulo the subgroup \\\\Gamma_0(3) of the modular group. The inverse function t=t(T) generates a solution of a general Halphen system. UR - http://hdl.handle.net/1963/2540 U1 - 1579 U2 - Mathematics U3 - Mathematical Physics ER - TY - THES T1 - On the Cauchy Problem for the Whitham Equations Y1 - 1998 A1 - Tamara Grava KW - Korteweg de Vries equation PB - SISSA UR - http://hdl.handle.net/1963/5555 U1 - 5382 U2 - Mathematics U3 - Mathematical Physics U4 - -1 ER -