BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

C. Klein

C. Klein appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
nlin.SImath.APmath-phmath.MPmath.NAmath.AGNumerical Analysisastro-ph.HEComputational Geometrynlin.PSastro-ph.COastro-ph.EPgr-qcmath.CAphysics.flu-dynquant-ph

Catalog footprint

What is connected

35works
16topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

35 published item(s)

preprint2025arXiv

All-sky neutrino point-source search with IceCube combined track and cascade data

Despite extensive efforts, discovery of high-energy astrophysical neutrino sources remains elusive. We present an event-level simultaneous maximum likelihood analysis of tracks and cascades using IceCube data collected from 04/06/2008 to 05/23/2022 to search the whole sky for neutrino sources and, using a source catalog, for coincidence of neutrino emission with gamma-ray emission. This is the first time a simultaneous fit of different detection channels is used to conduct a time-integrated all-sky scan with IceCube. Combining all-sky tracks, with superior pointing-power and sensitivity in the northern sky, with all-sky cascades, with good energy-resolution and sensitivity in the southern sky, we have developed the most sensitive point-source search to date by IceCube which targets the entire sky. The most significant point in the northern sky aligns with NGC 1068, a Seyfert II galaxy, which, from the catalog search, shows a 3.5$σ$ excess over background after accounting for trials. The most significant point in the southern sky does not align with any source in the catalog and is not significant after accounting for trials. A search for the single most significant Gaussian flare at the locations of NGC 1068, PKS 1424+240, and the southern highest significance point shows results consistent with expectations for steady emission. Notably, this is the first time that a flare shorter than four years has been excluded as being responsible for NGC 1068's emergence as a neutrino source. Our results show that combining tracks and cascades when conducting neutrino source searches improves sensitivity and can lead to new discoveries.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Multi-domain spectral approach for the Hilbert transform on the real line

A multi-domain spectral method is presented to compute the Hilbert transform on the whole compactified real line, with a special focus on piece-wise analytic functions and functions with algebraic decay towards infinity. Several examples of these and other types of functions are discussed. As an application solitons to generalized Benjamin-Ono equations are constructed.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Numerical scattering for the defocusing Davey-Stewartson II equation for initial data with compact support

In this work we present spectral algorithms for the numerical scattering for the defocusing Davey-Stewartson (DS) II equation with initial data having compact support on a disk, i.e., for the solution of d-bar problems. Our algorithms use polar coordinates and implement a Chebychev spectral scheme for the radial dependence and a Fourier spectral method for the azimuthal dependence. The focus is placed on the construction of complex geometric optics (CGO) solutions which are needed in the scattering approach for DS. We discuss two different approaches: The first constructs a fundamental solution to the d-bar system and applies the CGO conditions on the latter. This is especially efficient for small values of the modulus of the spectral parameter $k$. The second approach uses a fixed point iteration on a reformulated d-bar system containing the spectral parameter explicitly, a price paid to have simpler asymptotics. The approaches are illustrated for the example of the characteristic function of the disk and are shown to exhibit spectral convergence, i.e., an exponential decay of the numerical error with the number of collocation points. An asymptotic formula for large $|k|$ is given for the reflection coefficient.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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ödinger (NLS) equation and thus a $y$-independent solution to the 2D NLS. It is shown that the Peregrine solution is unstable against all standard perturbations, and that some perturbations can even lead to a blow-up for the elliptic NLS equation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Numerical study of the Kerr solution in rotating coordinates

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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 shock front develops. A local expansion reveals the universal scaling structure of the shock, which after a suitable change of coordinates corresponds to a generic cusp catastrophe. We provide a heuristic derivation of the shock front position near the critical point for the solution of the dKP equation, and study the solution of the dKP equation when a small amount of dissipation is added. Using multiple-scale analysis, we show that in the limit of small dissipation and near the critical point of the dKP solution, the solution of the dissipative dKP equation converges to a Pearcey integral. We test and illustrate our results by detailed comparisons with numerical simulations of both the regularized equation, the dKP equation, and the asymptotic description given in terms of the Pearcey integral.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Computational approach to compact Riemann surfaces

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 collocation points along these contours and by analytically continuing the roots. The collocation points are chosen to correspond to Chebychev collocation points for an ensuing Clenshaw-Curtis integration of the holomorphic differentials which gives the periods of the Riemann surface with spectral accuracy. At the singularities of the algebraic curve, Puiseux expansions computed by contour integration on the circles around the singularities are used to identify the holomorphic differentials. The Abel map is also computed with the Clenshaw-Curtis algorithm and contour integrals. As an application of the code, solutions to the Kadomtsev-Petviashvili equation are computed on non-hyperelliptic Riemann surfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Numerical study of the stability of the Peregrine breather

The Peregrine breather is widely discussed as a model for rogue waves in deep water. We present here a detailed numerical study of perturbations of the Peregrine breather as a solution to the nonlinear Schrödinger (NLS) equations. We first address the modulational instability of the constant modulus solution to NLS. Then we study numerically localized and nonlocalized perturbations of the Peregrine breather in the linear and fully nonlinear setting. It is shown that the solution is unstable against all considered perturbations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the tritronquée solutions of P$_I^2$

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,...,6$, called {\em tritronquée}, 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ée solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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 which is solved with Fourier techniques. The singular integrand is regularized analytically. The resulting integral equation is approximated via a discrete system which is solved with Krylov methods. As an example, the D-bar problem for the Davey-Stewartson II equations is solved. The result is used to test direct numerical solutions of the PDE.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Multidomain spectral method for Schrödinger equations

A multidomain spectral method with compactified exterior domains combined with stable second and fourth order time integrators is presented for Schrödinger equations. The numerical approach allows high precision numerical studies of solutions on the whole real line. At examples for the linear and cubic nonlinear Schrödinger equation, this code is compared to transparent boundary conditions and perfectly matched layers approaches. The code can deal with asymptotically non vanishing solutions as the Peregrine breather being discussed as a model for rogue waves. It is shown that the Peregrine breather can be numerically propagated with essentially machine precision, and that localized perturbations of this solution can be studied.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Numerical study of blow-up in solutions to generalized Korteweg-de Vries equations

We present a detailed numerical study of solutions to general Korteweg-de Vries equations with critical and supercritical nonlinearity. 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 Raphaël 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 which indicates that the theory by Martel, Merle and Raphaël is also applicable to initial data with a mass much larger than the soliton mass. We study the scaling of the blow-up time $t^{*}$ in dependence of the small dispersion parameter $ε$ and find an exponential dependence $t^{*}(ε)$ and that there is a minimal blow-up time $t^{*}_{0}$ greater than the critical time of the corresponding Hopf solution for $ε\to0$. To study the cases with blow-up in detail, we apply the first dynamic rescaling for generalized Korteweg-de Vries equations. This allows to identify the type of the singularity.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Numerical study of fractional Nonlinear Schrödinger equations

Using a Fourier spectral method, we provide a detailed numerically investigation of dispersive Schrödinger type equations involving a fractional Laplacian. By an appropriate choice of the dispersive exponent, both mass and energy sub- and supercritical regimes can be computed in one spatial dimension, only. This allows us to study the possibility of finite time blow-up versus global existence, the nature of the blow-up, the stability and instability of nonlinear ground states, and the long time dynamics of solutions. The latter is also studied in a semiclassical setting. Moreover, we numerically construct ground state solutions to the fractional nonlinear Schrödinger equation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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ödinger equation in $1+1$ dimensions. The full Toda system is then studied for finite small values of the dispersion parameter $ε$ in the vicinity of the shocks of the dispersionless Toda equations. We determine the scaling in $ε$ of the difference between the Toda solution for small $ε$ and the singular solution of the dispersionless Toda system. In the hyperbolic case, the same scaling proportional to $ε^{2/7}$ is found as in the small dispersion limit of the Korteweg-de Vries and the defocusing nonlinear Schrödinger equations. In the elliptic case, we obtain the same scaling proportional to $ε^{2/5}$ as in the semiclassical limit for the focusing nonlinear Schrödinger equation. We also study the formation of dispersive shocks for times much larger than the break-up time in the hyperbolic case. In the elliptic case, an $L_{\infty}$ blow-up is observed instead of a dispersive shock for finite times greater than the break-up time. The $ε$-dependence of the blow-up time is determined.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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, are numerically integrated up to the formation of a shock. The use of parallelized algorithms allows to determine the critical time $t_{c}$ 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 Schrödinger equation, i.e., cubic singularities in the defocusing case and square root singularities in the focusing case. For small values of $ε$, the full Davey-Stewartson II equations are integrated for the same initial data up to the critical time $t_{c}$. The scaling in $ε$ of the difference between these solutions is found to be the same as in the $1+1$ dimensional case, proportional to $ε^{2/7}$ for the defocusing case and proportional to $ε^{2/5}$ in the focusing case. We document the Davey-Stewartson II solutions for small $ε$ for times much larger than the critical time $t_{c}$. It is shown that zones of rapid modulated oscillations are formed near the shocks of the solutions to the semiclassical equations. For smaller $ε$, the oscillatory zones become smaller and more sharply delimited to lens shaped regions. Rapid oscillations are also found in the focusing case for initial data where the singularities of the solution to the semiclassical equations do not coincide.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On critical behaviour in systems of Hamiltonian partial differential equations

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$) equation or its fourth order analogue P$_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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On the integrable elliptic cylindrical Kadomtsev-Petviashvili equation

There exist two versions of the Kadomtsev-Petviashvili equation, related to the Cartesian and cylindrical geometries of the waves. In this paper we derive and study a new version, related to the elliptic cylindrical geometry. The derivation is given in the context of surface waves, but the derived equation is a universal integrable model applicable to generic weakly-nonlinear weakly-dispersive waves. We also show that there exist nontrivial transformations between all three versions of the KP equation associated with the physical problem formulation, and use them to obtain new classes of approximate solutions for water waves.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Asteroid rotation periods from the Palomar Transient Factory survey

The Palomar Transient Factory (PTF) is a synoptic survey designed to explore the transient and variable sky in a wide variety of cadences. We use PTF observations of fields that were observed multiple times (>=10) per night, for several nights, to find asteroids, construct their lightcurves and measure their rotation periods. Here we describe the pipeline we use to achieve these goals and present the results from the first four (overlapping) PTF fields analyzed as part of this program. These fields, which cover an area of 21 deg^2, were observed on four nights with a cadence of ~20 min. Our pipeline was able to detect 624 asteroids, of which 145 (~20%) were previously unknown. We present high quality rotation periods for 88 main-belt asteroids and possible period or lower limit on the period for an additional 85 asteroids. For the remaining 451 asteroids, we present lower limits on their photometric amplitudes. Three of the asteroids have lightcurves that are characteristic of binary asteroids. We estimate that implementing our search for all existing high-cadence PTF data will provide rotation periods for about 10,000 asteroids mainly in the magnitude range ~14 to ~20.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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 $τ$, 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 $\mathcal{A}$-cycles are invariant under the anti-holomorphic involution $τ$.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

WKB analysis of Bohmian dynamics

We consider a semi-classically scaled Schrödinger equation with WKB initial data. We prove that in the classical limit the corresponding Bohmian trajectories converge (locally in measure) to the classical trajectories before the appearance of the first caustic. In a second step we show that after caustic onset this convergence in general no longer holds. In addition, we provide numerical simulations of the Bohmian trajectories in the semiclassical regime which illustrate the above results.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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 group.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Fourth order time-stepping for Kadomtsev-Petviashvili and Davey-Stewartson equations

Purely dispersive partial differential equations as the Korteweg-de Vries equation, the nonlinear Schrödinger equation and higher dimensional generalizations thereof can have solutions which develop a zone of rapid modulated oscillations in the region where the corresponding dispersionless equations have shocks or blow-up. To numerically study such phenomena, fourth order time-stepping in combination with spectral methods is beneficial to resolve the steep gradients in the oscillatory region. We compare the performance of several fourth order methods for the Kadomtsev-Petviashvili and the Davey-Stewartson equations, two integrable equations in 2+1 dimensions: exponential time-differencing, integrating factors, time-splitting, implicit Runge-Kutta and Driscoll's composite Runge-Kutta method. The accuracy in the numerical conservation of integrals of motion is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Multiwavelength campaign on Mrk 509. I. Variability and spectral energy distribution

(Abridged) Active galactic nuclei show a wealth of interesting physical processes, some of which are poorly understood. We want to address a number of open questions, including the location and physics of the outflow from AGN, the nature of the continuum emission, the geometry and physical state of the X-ray broad emission line region, the Fe-K line complex, the metal abundances of the nucleus and finally the interstellar medium of our own Galaxy. We study one of the best targets for these aims, the Seyfert 1 galaxy Mrk 509 with a multiwavelength campaign using five satellites (XMM-Newton, INTEGRAL, Chandra, HST and Swift) and two ground-based facilities (WHT and PAIRITEL). Our observations cover more than five decades in frequency, from 2 um to 200 keV. The combination of high-resolution spectroscopy and time variability allows us to disentangle and study the different components. Our campaign covers 100 days from September to December 2009, and is centred on a simultaneous set of deep XMM-Newton and INTEGRAL observations with regular time intervals, spanning seven weeks. We obtain a continuous light curve in the X-ray and UV band, showing a strong, up to 60% flux increase in the soft X-ray band during the three weeks in the middle of our deepest monitoring campaign, and which is correlated with an enhancement of the UV flux. This allows us to study the time evolution of the continuum and the outflow. By stacking the observations, we have also obtained one of the best X-ray and UV spectra of a Seyfert galaxy ever obtained. In this paper we also study the effects of the spectral energy distribution (SED) that we obtained on the photo-ionisation equilibrium. Thanks to our broad-band coverage, uncertainties on the SED do not strongly affect the determination of this equilibrium.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Numerical Study of Blowup in the Davey-Stewartson System

Nonlinear dispersive partial differential equations such as the nonlinear Schrödinger equations can have solutions that blow-up. We numerically study the long time behavior and potential blowup of solutions to the focusing Davey-Stewartson II equation by analyzing perturbations of the lump and the Ozawa exact solutions. It is shown in this way that the lump is unstable to both blowup and dispersion, and that blowup in the Ozawa solution is generic.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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 not related to automorphisms of the curve, we study symplectic transformations to an adapted basis and give explicit formulae for M-curves. As examples we discuss solutions of the Davey-Stewartson and the multi-component nonlinear Schrödinger equations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Transverse stability of periodic traveling waves in Kadomtsev-Petviashvili equations: A numerical study

We numerically investigate transverse stability and instability of so-called cnoidal waves, i.e., periodic traveling wave solutions of the Korteweg-de Vries equation, under the time-evolution of the Kadomtsev-Petviashvili equation. In particular, we find that in KP-I small amplitude cnoidal waves are stable (at least for spatially localized perturbations) and only become unstable above a certain threshold. In contrast to that, KP-II is found to be stable for all amplitudes, or, equivalently, wave speeds. This is in accordance with recent analytical results for solitary waves given in \cite{RT1, RT2}.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Numerical Study of breakup in generalized Korteweg-de Vries and Kawahara equations

This article is concerned with a conjecture by one of the authors on the formation of dispersive shocks in a class of Hamiltonian dispersive regularizations of the quasilinear 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 behaviour of their solution near the point of gradient catastrophe for the transport equation is described locally by a special solution of a Painlevé-type equation. This local description holds also for solutions to equations where blow up 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 $ε^2$ where $ε^2$ is the small dispersion parameter. Corrections up to order $ε^4 $ are obtained and tested numerically.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On linear degeneracy of integrable quasilinear systems in higher dimensions

We investigate $(d+1)$-dimensional quasilinear systems which are integrable by the method of hydrodynamic reductions. In the case $d\geq 3$ we formulate a conjecture that any such system with an irreducible dispersion relation must be linearly degenerate. We prove this conjecture in the 2-component case, providing a complete classification of multi-dimensional integrable systems in question. In particular, our results imply the non-existence of 2-component integrable systems of hydrodynamic type for $d\geq 6$. In the second half of the paper we discuss a numerical and analytical evidence for the impossibility of the breakdown of smooth initial data for linearly degenerate systems in 2+1 dimensions.

0 citations0 reviewsNo code linkedOpen access