Spectral approach to the scattering map for the semi-classical defocusing Davey–Stewartson II equation

International audience; The inverse scattering approach for the defocusing Davey–Stewartson II equation is given by a system of D-bar equations. We present a numerical approach to semi-classical D-bar problems for real analytic rapidly decreasing potentials. We treat the D-bar problem as a complex linear second order integral equation which is solved with discrete Fourier transforms complemented by a regularization of the singular parts by explicit analytic computation. The resulting algebraic equation is solved either by fixed point iterations or GMRES. Several examples for small values of the semi-classical parameter in the system are discussed.

Shock formation in the dispersionless Kadomtsev-Petviashvili equation

The dispersionless Kadomtsev-Petviashvili (dKP) equation $(u_t+uu_x)_x=u_{yy}$ is one of the simplest nonlinear wave equations describing two-dimensional shocks. To solve the dKP equation we use a coordinate transformation inspired by the method of characteristics for the one-dimensional Hopf equation $u_t+uu_x=0$. We show numerically that the solutions to the transformed equation do not develop shocks. This permits us to extend the dKP solution as the graph of a multivalued function beyond the critical time when the gradients blow up. This overturned solution is multivalued in a lip shape region in the $(x,y)$ plane, where the solution of the dKP equation exists in a weak sense only, and a…

Computation of the topological type of a real Riemann surface

We present an algorithm for the computation of the topological type of a real compact Riemann surface associated to an algebraic curve, i.e., its genus and the properties of the set of fixed points of the anti-holomorphic involution τ \tau , namely, the number of its connected components, and whether this set divides the surface into one or two connected components. This is achieved by transforming an arbitrary canonical homology basis to a homology basis where the A \mathcal {A} -cycles are invariant under the anti-holomorphic involution  τ \tau .

Numerical study of the Kerr solution in rotating coordinates

International audience; The Kerr solution in coordinates corotating with the horizon is studied as a testbed for a spacetime with a helical Killing vector in the Ernst picture. The solution is numerically constructed by solving the Ernst equation with a spectral method and a Newton iteration. We discuss convergence of the iteration for several initial iterates and different values of the Kerr parameters.

Numerical study of blow-up and stability of line solitons for the Novikov-Veselov equation

International audience; We study numerically the evolution of perturbed Korteweg-de Vries solitons and of well localized initial data by the Novikov-Veselov (NV) equation at different levels of the 'energy' parameter E. We show that as |E| -> infinity, NV behaves, as expected, similarly to its formal limit, the Kadomtsev-Petviashvili equation. However at intermediate regimes, i.e. when |E| is not very large, more varied scenarios are possible, in particular, blow-ups are observed. The mechanism of the blow-up is studied.

Erratum to: IST Versus PDE: A Comparative Study

Computational approach to compact Riemann surfaces

International audience; A purely numerical approach to compact Riemann surfaces starting from plane algebraic curves is presented. The critical points of the algebraic curve are computed via a two-dimensional Newton iteration. The starting values for this iteration are obtained from the resultants with respect to both coordinates of the algebraic curve and a suitable pairing of their zeros. A set of generators of the fundamental group for the complement of these critical points in the complex plane is constructed from circles around these points and connecting lines obtained from a minimal spanning tree. The monodromies are computed by solving the defining equation of the algebraic curve on…

Numerical Study of the semiclassical limit of the Davey-Stewartson II equations

We present the first detailed numerical study of the semiclassical limit of the Davey–Stewartson II equations both for the focusing and the defocusing variant. We concentrate on rapidly decreasing initial data with a single hump. The formal limit of these equations for vanishing semiclassical parameter , the semiclassical equations, is numerically integrated up to the formation of a shock. The use of parallelized algorithms allows one to determine the critical time tc and the critical solution for these 2 + 1-dimensional shocks. It is shown that the solutions generically break in isolated points similarly to the case of the 1 + 1-dimensional cubic nonlinear Schrodinger equation, i.e., cubic…

A Study of the Direct Spectral Transform for the Defocusing Davey‐Stewartson II Equation the Semiclassical Limit

International audience; The defocusing Davey-Stewartson II equation has been shown in numerical experiments to exhibit behavior in the semiclassical limit that qualitatively resembles that of its one-dimensional reduction, the defocusing nonlinear Schrodinger equation, namely the generation from smooth initial data of regular rapid oscillations occupying domains of space-time that become well-defined in the limit. As a first step to studying this problem analytically using the inverse scattering transform, we consider the direct spectral transform for the defocusing Davey-Stewartson II equation for smooth initial data in the semiclassical limit. The direct spectral transform involves a sing…

Numerical study of a multiscale expansion of the Korteweg de Vries equation and Painlev\'e-II equation

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…

Nonlinear Dispersive Equations

Nonlinear Dispersive Equations are partial differential equations that naturally arise in physical settings where dispersion dominates dissipation, notably hydrodynamics, nonlinear optics, plasma physics and Bose–Einstein condensates. The topic has traditionally been approached in different ways, from the perspective of modeling of physical phenomena, to that of the theory of partial differential equations, or as part of the theory of integrable systems.This monograph offers a thorough introduction to the topic, uniting the modeling, PDE and integrable systems approaches for the first time in book form. The presentation focuses on three "universal" families of physically relevant equations …

Efficient computation of the branching structure of an algebraic curve

An efficient algorithm for computing the branching structure of a compact Riemann surface defined via an algebraic curve is presented. Generators of the fundamental group of the base of the ramified covering punctured at the discriminant points of the curve are constructed via a minimal spanning tree of the discriminant points. This leads to paths of minimal length between the points, which is important for a later stage where these paths are used as integration contours to compute periods of the surface. The branching structure of the surface is obtained by analytically continuing the roots of the equation defining the algebraic curve along the constructed generators of the fundamental gro…

Numerical study of shock formation in the dispersionless Kadomtsev-Petviashvili equation and dispersive regularizations

The formation of singularities in solutions to the dispersionless Kadomtsev-Petviashvili (dKP) equation is studied numerically for different classes of initial data. The asymptotic behavior of the Fourier coefficients is used to quantitatively identify the critical time and location and the type of the singularity. The approach is first tested in detail in 1+1 dimensions for the known case of the Hopf equation, where it is shown that the break-up of the solution can be identified with prescribed accuracy. For dissipative regularizations of this shock formation as the Burgers' equation and for dispersive regularizations as the Korteweg-de Vries equation, the Fourier coefficients indicate as …

On critical behaviour in generalized Kadomtsev-Petviashvili equations

International audience; An asymptotic description of the formation of dispersive shock waves in solutions to the generalized Kadomtsev–Petviashvili (KP) equation is conjectured. The asymptotic description based on a multiscales expansion is given in terms of a special solution to an ordinary differential equation of the Painlevé I hierarchy. Several examples are discussed numerically to provide strong evidence for the validity of the conjecture. The numerical study of the long time behaviour of these examples indicates persistence of dispersive shock waves in solutions to the (subcritical) KP equations, while in the supercritical KP equations a blow-up occurs after the formation of the disp…

Numerical Study of Blow-Up Mechanisms for Davey-Stewartson II Systems

We present a detailed numerical study of various blow-up issues in the context of the focusing Davey-Stewartson II equation. To this end we study Gaussian initial data and perturbations of the lump and the explicit blow-up solution due to Ozawa. Based on the numerical results it is conjectured that the blow-up in all cases is self similar, and that the time dependent scaling is as in the Ozawa solution and not as in the stable blow-up of standard $L^{2}$ critical nonlinear Schr\"odinger equations. The blow-up profile is given by a dynamically rescaled lump.

IST Versus PDE: A Comparative Study

We survey and compare, mainly in the two-dimensional case, various results obtained by IST and PDE techniques for integrable equations. We also comment on what can be predicted from integrable equations on non integrable ones.

Extremal properties of the determinant of the Laplacian in the Bergman metric on the moduli space of genus two Riemann surfaces

We study extremal properties of the determinant of the Laplacian in the Bergman metric on the moduli space of compact genus two Riemann surfaces. By a combination of analytical and numerical methods we identify four non-degenerate critical points of this function and compute the signature of the Hessian at these points. The curve with the maximal number of automorphisms (the Burnside curve) turns out to be the point of the absolute maximum. Our results agree with the mass formula for orbifold Euler characteristics of the moduli space. A similar analysis is performed for the Bolza's strata of symmetric Riemann surfaces of genus two.

Numerical study of blow-up in solutions to generalized Kadomtsev-Petviashvili equations

We present a numerical study of solutions to the generalized Kadomtsev-Petviashvili equations with critical and supercritical nonlinearity for localized initial data with a single minimum and single maximum. In the cases with blow-up, we use a dynamic rescaling to identify the type of the singularity. We present a discussion of the observed blow-up scenarios.

Numerical study of the stability of the Peregrine solution

International audience; The Peregrine solution to the nonlinear Schrödinger equations is widely discussed as a model for rogue waves in deep water. We present here a detailed fully nonlinear numerical study of high accuracy of perturbations of the Peregrine solution as a solution to the nonlinear Schrödinger (NLS) equations.We study localized and nonlocalized perturbations of the Peregrine solution in the linear and fully nonlinear setting. It is shown that the solution is unstable against all considered perturbations.

Numerical study of the Kadomtsev–Petviashvili equation and dispersive shock waves

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 Schrodinger equation in the semiclassical limit.

Numerical study of soliton stability, resolution and interactions in the 3D Zakharov–Kuznetsov equation

International audience; We present a detailed numerical study of solutions to the Zakharov-Kuznetsov equation in three spatial dimensions. The equation is a three-dimensional generalization of the Korteweg-de Vries equation, though, not completely integrable. This equation is L-2-subcritical, and thus, solutions exist globally, for example, in the H-1 energy space.We first study stability of solitons with various perturbations in sizes and symmetry, and show asymptotic stability and formation of radiation, confirming the asymptotic stability result in Farah et al. (0000) for a larger class of initial data. We then investigate the solution behavior for different localizations and rates of de…

New construction of algebro-geometric solutions to the Camassa-Holm equation and their numerical evaluation

An independent derivation of solutions to the Camassa-Holm equation in terms of multi-dimensional theta functions is presented using an approach based on Fay's identities. Reality and smoothness conditions are studied for these solutions from the point of view of the topology of the underlying real hyperelliptic surface. The solutions are studied numerically for concrete examples, also in the limit where the surface degenerates to the Riemann sphere, and where solitons and cuspons appear.

A numerical approach to Blow-up issues for dispersive perturbations of Burgers' equation

We provide a detailed numerical study of various issues pertaining to the dynamics of the Burgers equation perturbed by a weak dispersive term: blow-up in finite time versus global existence, nature of the blow-up, existence for "long" times, and the decomposition of the initial data into solitary waves plus radiation. We numerically construct solitons for fractionary Korteweg-de Vries equations.

High precision numerical approach for Davey–Stewartson II type equations for Schwartz class initial data

We present an efficient high-precision numerical approach for Davey–Stewartson (DS) II type equa- tions, treating initial data from the Schwartz class of smooth, rapidly decreasing functions. As with previous approaches, the presented code uses discrete Fourier transforms for the spatial dependence and Driscoll’s composite Runge–Kutta method for the time dependence. Since DS equations are non-local, nonlinear Schrödinger equations with a singular symbol for the non-locality, standard Fourier methods in practice only reach accuracy of the order of 10−6or less for typical examples. This was previously demonstrated for the defocusing integrable case by comparison with a numerical approach for …

Numerical study of the transverse stability of the Peregrine solution

We generalise a previously published approach based on a multi-domain spectral method on the whole real line in two ways: firstly, a fully explicit 4th order method for the time integration, based on a splitting scheme and an implicit Runge--Kutta method for the linear part, is presented. Secondly, the 1D code is combined with a Fourier spectral method in the transverse variable both for elliptic and hyperbolic NLS equations. As an example we study the transverse stability of the Peregrine solution, an exact solution to the one dimensional nonlinear Schr\"odinger (NLS) equation and thus a $y$-independent solution to the 2D NLS. It is shown that the Peregrine solution is unstable against all…

Sirolimus Use in Liver Transplant Recipients With Hepatocellular Carcinoma: A Randomized, Multicenter, Open-Label Phase 3 Trial

International audience; BACKGROUND:We investigated whether sirolimus-based immunosuppression improves outcomes in liver transplantation (LTx) candidates with hepatocellular carcinoma (HCC).METHODS:In a prospective-randomized open-label international trial, 525 LTx recipients with HCC initially receiving mammalian target of rapamycin inhibitor-free immunosuppression were randomized 4 to 6 weeks after transplantation into a group on mammalian target of rapamycin inhibitor-free immunosuppression (group A: 264 patients) or a group incorporating sirolimus (group B: 261). The primary endpoint was recurrence-free survival (RFS); intention-to-treat (ITT) analysis was conducted after 8 years. Overal…

Numerical study of blow-up and dispersive shocks in solutions to generalized Korteweg–de Vries equations

Abstract We present a detailed numerical study of solutions to general Korteweg–de Vries equations with critical and supercritical nonlinearity, both in the context of dispersive shocks and blow-up. We study the stability of solitons and show that they are unstable against being radiated away and blow-up. In the L 2 critical case, the blow-up mechanism by Martel, Merle and Raphael can be numerically identified. In the limit of small dispersion, it is shown that a dispersive shock always appears before an eventual blow-up. In the latter case, always the first soliton to appear will blow up. It is shown that the same type of blow-up as for the perturbations of the soliton can be observed whic…

Spectral approach to D-bar problems

We present the first numerical approach to D-bar problems having spectral convergence for real analytic, rapidly decreasing potentials. The proposed method starts from a formulation of the problem in terms of an integral equation that is numerically solved with Fourier techniques. The singular integrand is regularized analytically. The resulting integral equation is approximated via a discrete system that is solved with Krylov methods. As an example, the D-bar problem for the Davey-Stewartson II equations is considered. The result is used to test direct numerical solutions of the PDE.© 2017 Wiley Periodicals, Inc.

Multidomain spectral method for the Gauss hypergeometric function

International audience; We present a multidomain spectral approach for Fuchsian ordinary differential equations in the particular case of the hypergeometric equation. Our hybrid approach uses Frobenius’ method and Moebius transformations in the vicinity of each of the singular points of the hypergeometric equation, which leads to a natural decomposition of the real axis into domains. In each domain, solutions to the hypergeometric equation are constructed via the well-conditioned ultraspherical spectral method. The solutions are matched at the domain boundaries to lead to a solution which is analytic on the whole compactified real line R∪∞, except for the singular points and cuts of the Rie…

A numerical study of the small dispersion limit of the Korteweg-de Vries equation and asymptotic solutions

Abstract We study numerically the small dispersion limit for the Korteweg–de Vries (KdV) equation u t + 6 u u x + ϵ 2 u x x x = 0 for ϵ ≪ 1 and give a quantitative comparison of the numerical solution with various asymptotic formulae for small ϵ in the whole ( x , t ) -plane. The matching of the asymptotic solutions is studied numerically.

On Whitham and Related Equations

The aim of this paper is to study, via theoretical analysis and numerical simulations, the dynamics of Whitham and related equations. In particular, we establish rigorous bounds between solutions of the Whitham and Korteweg–de Vries equations and provide some insights into the dynamics of the Whitham equation in different regimes, some of them being outside the range of validity of the Whitham equation as a water waves model.

Numerical study of the long wavelength limit of the Toda lattice

We present the first detailed numerical study of the Toda equations in $2+1$ dimensions in the limit of long wavelengths, both for the hyperbolic and elliptic case. We first study the formal dispersionless limit of the Toda equations and solve initial value problems for the resulting system up to the point of gradient catastrophe. It is shown that the break-up of the solution in the hyperbolic case is similar to the shock formation in the Hopf equation, a $1+1$ dimensional singularity. In the elliptic case, it is found that the break-up is given by a cusp as for the semiclassical system of the focusing nonlinear Schr\"odinger equation in $1+1$ dimensions. The full Toda system is then studie…

On a nonlinear Schrödinger equation for nucleons in one space dimension

We study a 1D nonlinear Schrödinger equation appearing in the description of a particle inside an atomic nucleus. For various nonlinearities, the ground states are discussed and given in explicit form. Their stability is studied numerically via the time evolution of perturbed ground states. In the time evolution of general localized initial data, they are shown to appear in the long time behaviour of certain cases.

On critical behaviour in systems of Hamiltonian partial differential equations

Abstract 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é-I (P $$_I$$ I ) equation or its fourth-order analogue P $$_I^2$$ I 2 . As concrete examples, we discuss nonlinear Schrödinger equations in the semiclassical limit. A numerical study of these cases provides strong evidence in support of the conjecture.

On the numerical evaluation of algebro-geometric solutions to integrable equations

Physically meaningful periodic solutions to certain integrable partial differential equations are given in terms of multi-dimensional theta functions associated to real Riemann surfaces. Typical analytical problems in the numerical evaluation of these solutions are studied. In the case of hyperelliptic surfaces efficient algorithms exist even for almost degenerate surfaces. This allows the numerical study of solitonic limits. For general real Riemann surfaces, the choice of a homology basis adapted to the anti-holomorphic involution is important for a convenient formulation of the solutions and smoothness conditions. Since existing algorithms for algebraic curves produce a homology basis no…

Algebraic Curves and Riemann Surfaces in Matlab

In the previous chapter, a detailed description of the algorithms for the ‘algcurves’ package in Maple was presented. As discussed there, the package is able to handle general algebraic curves with coefficients given as exact arithmetic expressions, a restriction due to the use of exact integer arithmetic. Coefficients in terms of floating point numbers, i.e., the representation of decimal numbers of finite length on a computer, can in principle be handled, but the floating point numbers have to be converted to rational numbers. This can lead to technical difficulties in practice. One also faces limitations if one wants to study families of Riemann surfaces, where the coefficients in the al…

On universality of critical behavior in the focusing nonlinear Schrödinger equation, elliptic umbilic catastrophe and the Tritronquée solution to the Painlevé-I equation

We argue that the critical behavior near the point of “gradient catastrophe” of the solution to the Cauchy problem for the focusing nonlinear Schrodinger equation \(i\epsilon \varPsi _{t}+\frac{\epsilon^{2}}{2}\varPsi _{xx}+|\varPsi |^{2}\varPsi =0\) , e ≪1, with analytic initial data of the form \(\varPsi (x,0;\epsilon)=A(x)e^{\frac{i}{\epsilon}S(x)}\) is approximately described by a particular solution to the Painleve-I equation.

Multi-domain spectral approach with Sommerfeld condition for the Maxwell equations

We present a multidomain spectral approach with an exterior compactified domain for the Maxwell equations for monochromatic fields. The Sommerfeld radiation condition is imposed exactly at infinity being a finite point on the numerical grid. As an example, axisymmetric situations in spherical and prolate spheroidal coordinates are discussed.