Source author record

J. Lenells

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
3topics
1close collaborators

Actions

Connect this record

Log in to claim

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

11 published item(s)

preprint2021arXiv

A new approach to integrable evolution equations on the circle

We propose a new approach for the solution of initial value problems for integrable evolution equations in the periodic setting based on the unified transform. Using the nonlinear Schrödinger equation as a model example, we show that the solution of the initial value problem on the circle can be expressed in terms of the solution of a Riemann-Hilbert problem whose formulation involves quantities which are defined in terms of the initial data alone. Our approach provides an effective solution of the problem on the circle which is conceptually analogous to the solution of the problem on the line via the inverse scattering transform.

preprint2020arXiv

Linearizable boundary value problems for the elliptic sine-Gordon and the elliptic Ernst equations

By employing a novel generalization of the inverse scattering transform method known as the unified transform or Fokas method, it can be shown that the solution of certain physically significant boundary value problems for the elliptic sine-Gordon equation, as well as for the elliptic version of the Ernst equation, can be expressed in terms of the solution of appropriate $2 \times 2$-matrix Riemann--Hilbert (RH) problems. These RH problems are defined in terms of certain functions, called spectral functions, which involve the given boundary conditions, but also unknown boundary values. For arbitrary boundary conditions, the determination of these unknown boundary values requires the analysis of a nonlinear Fredholm integral equation. However, there exist particular boundary conditions, called linearizable, for which it is possible to bypass this nonlinear step and to characterize the spectral functions directly in terms of the given boundary conditions. Here, we review the implementation of this effective procedure for the following linearizable boundary value problems: (a) the elliptic sine-Gordon equation in a semi-strip with zero Dirichlet boundary values on the unbounded sides and with constant Dirichlet boundary value on the bounded side; (b) the elliptic Ernst equation with boundary conditions corresponding to a uniformly rotating disk of dust; (c) the elliptic Ernst equation with boundary conditions corresponding to a disk rotating uniformly around a central black hole; (d) the elliptic Ernst equation with vanishing Neumann boundary values on a rotating disk.

preprint2015arXiv

The nonlinear Schrödinger equation with $t$-periodic data: I. Exact results

We consider the nonlinear Schrödinger equation on the half-line with a given Dirichlet (Neumann) boundary datum which for large $t$ tends to the periodic function $g_0^b(t)$ ($g_1^b(t)$). Assuming that the unknown Neumann (Dirichlet) boundary value tends for large $t$ to a periodic function $g_1^b(t)$ ($g_0^b(t)$), we derive an easily verifiable condition that the functions $g_0^b(t)$ and $g_1^b(t)$ must satisfy. Furthermore, we introduce two different methods, one based on the formulation of a Riemann-Hilbert problem, and one based on a perturbative approach, for constructing $g_1^b(t)$ ($g_0^b(t)$) in terms of $g_0^b(t)$ ($g_1^b(t)$).

preprint2015arXiv

The nonlinear Schrödinger equation with $t$-periodic data: II. Perturbative results

We consider the nonlinear Schrödinger equation on the half-line with a given Dirichlet boundary datum which for large $t$ tends to a periodic function. We assume that this function is sufficiently small, namely that it can be expressed in the form $αg_0^b(t)$, where $α$ is a small constant. Assuming that the Neumann boundary value tends for large $t$ to the periodic function $g_1^b(t)$, we show that $g_1^b(t)$ can be expressed in terms of a perturbation series in $α$ which can be constructed explicitly to any desired order. As an illustration, we compute $g_1^b(t)$ to order $α^8$ for the particular case that $g_0^b(t)$ is the sum of two exponentials. We also show that there exist particular functions $g_0^b(t)$ for which the above series can be summed up, and therefore for these functions $g_1^b(t)$ can be obtained in closed form. The simplest such function is $\exp(iωt)$, where $ω$ is a real constant.

preprint2014arXiv

Perturbative and exact results on the Neumann value for the nonlinear Schrödinger equation on the half-line

The most challenging problem in the implementation of the so-called \textit{unified transform} to the analysis of the nonlinear Schrödinger equation on the half-line is the characterization of the unknown boundary value in terms of the given initial and boundary conditions. For the so-called \textit{linearizable} boundary conditions this problem can be solved explicitly. Furthermore, for non-linearizable boundary conditions which decay for large $t$, this problem can be largely bypassed in the sense that the unified transform yields useful asymptotic information for the large $t$ behavior of the solution. However, for the physically important case of periodic boundary conditions it is necessary to characterize the unknown boundary value. Here, we first present a perturbative scheme which can be used to compute explicitly the asymptotic form of the Neumann boundary value in terms of the given $τ$-periodic Dirichlet datum to any given order in a perturbation expansion. We then discuss briefly an extension of the pioneering results of Boutet de Monvel and co-authors which suggests that if the Dirichlet datum belongs to a large class of particular $τ$-periodic functions, which includes $\{a \exp(i ωt) \, | \, a>0, \, ω\geq a^2\}$, then the large $t$ behavior of the Neumann value is given by a $τ$-periodic function which can be computed explicitly.

preprint2014arXiv

The Unified Transform for the modified Helmholtz equation in the exterior of a square

The Unified Transform provides a novel method for analyzing boundary value problems for linear and for integrable nonlinear PDEs. The numerical implementation of this method to linear elliptic PDEs formulated in the {\it interior} of a polygon has been investigated by several authors (see the article by Iserles, Smitheman, and one of the authors in this book). Here, we show that the Unified Transform also yields a novel numerical technique for computing the solution of linear elliptic PDEs in the {\it exterior} of a polygon. One of the advantages of this new technique is that it actually yields directly the scattering amplitude. Details are presented for the modified Helmholtz equation in the exterior of a square.

preprint2011arXiv

The Unified Method: I Non-Linearizable Problems on the Half-Line

Boundary value problems for integrable nonlinear evolution PDEs formulated on the half-line can be analyzed by the unified method introduced by one of the authors and used extensively in the literature. The implementation of this general method to this particular class of problems yields the solution in terms of the unique solution of a matrix Riemann-Hilbert problem formulated in the complex $k$-plane (the Fourier plane), which has a jump matrix with explicit $(x,t)$-dependence involving four scalar functions of $k$, called spectral functions. Two of these functions depend on the initial data, whereas the other two depend on all boundary values. The most difficult step of the new method is the characterization of the latter two spectral functions in terms of the given initial and boundary data, i.e. the elimination of the unknown boundary values. For certain boundary conditions, called linearizable, this can be achieved simply using algebraic manipulations. Here, we present an effective characterization of the spectral functions in terms of the given initial and boundary data for the general case of non-linearizable boundary conditions. This characterization is based on the analysis of the so-called global relation, on the analysis of the equations obtained from the global relation via certain transformations leaving the dispersion relation of the associated linearized PDE invariant, and on the computation of the large $k$ asymptotics of the eigenfunctions defining the relevant spectral functions.

preprint2011arXiv

The Unified Method: II NLS on the Half-Line with $t$-Periodic Boundary Conditions

Boundary value problems for integrable nonlinear evolution PDEs formulated on the half-line can be analyzed by the unified method introduced by one of the authors and used extensively in the literature. The implementation of this general method to this particular class of problems yields the solution in terms of the unique solution of a matrix Riemann-Hilbert problem formulated in the complex $k$-plane (the Fourier plane), which has a jump matrix with explicit $(x,t)$-dependence involving four scalar functions of $k$, called spectral functions. Two of these functions depend on the initial data, whereas the other two depend on all boundary values. The most difficult step of the new method is the characterization of the latter two spectral functions in terms of the given initial and boundary data, i.e. the elimination of the unknown boundary values. For certain boundary conditions, called linearizable, this can be achieved simply using algebraic manipulations. Here, we first present an effective characterization of the spectral functions in terms of the given initial and boundary data for the general case of non-linearizable boundary conditions. This characterization is based on the analysis of the so-called global relation and on the introduction of the so-called Gelfand-Levitan-Marchenko representations of the eigenfunctions defining the spectral functions. We then concentrate on the physically significant case of $t$-periodic Dirichlet boundary data. After presenting certain heuristic arguments which suggest that the Neumann boundary values become periodic as $t\to\infty$, we show that for the case of the NLS with a sine-wave as Dirichlet data, the asymptotics of the Neumann boundary values can be computed explicitly at least up to third order in a perturbative expansion and indeed at least up to this order are asymptotically periodic.

preprint2011arXiv

The Unified Method: III Non-Linearizable Problems on the Interval

Boundary value problems for integrable nonlinear evolution PDEs formulated on the finite interval can be analyzed by the unified method introduced by one of the authors and used extensively in the literature. The implementation of this general method to this particular class of problems yields the solution in terms of the unique solution of a matrix Riemann-Hilbert problem formulated in the complex $k$-plane (the Fourier plane), which has a jump matrix with explicit $(x,t)$-dependence involving six scalar functions of $k$, called spectral functions. Two of these functions depend on the initial data, whereas the other four depend on all boundary values. The most difficult step of the new method is the characterization of the latter four spectral functions in terms of the given initial and boundary data, i.e. the elimination of the unknown boundary values. Here, we present an effective characterization of the spectral functions in terms of the given initial and boundary data. We present two different characterizations of this problem. One is based on the analysis of the so-called global relation, on the analysis of the equations obtained from the global relation via certain transformations leaving the dispersion relation of the associated linearized PDE invariant, and on the computation of the large $k$ asymptotics of the eigenfunctions defining the relevant spectral functions. The other is based on the analysis of the global relation and on the introduction of the so-called Gelfand-Levitan-Marchenko representations of the eigenfunctions defining the relevant spectral functions. We also show that these two different characterizations are equivalent and that in the limit when the length of the interval tends to infinity, the relevant formulas reduce to the analogous formulas obtained recently for the case of boundary value problems formulated on the half-line.

preprint2010arXiv

Explicit soliton asymptotics for the Korteweg-de Vries equation on the half-line

Integrable PDEs on the line can be analyzed by the so-called Inverse Scattering Transform (IST) method. A particularly powerful aspect of the IST is its ability to predict the large $t$ behavior of the solution. Namely, starting with initial data $u(x,0)$, IST implies that the solution $u(x,t)$ asymptotes to a collection of solitons as $t \to \infty$, $x/t = O(1)$; moreover the shapes and speeds of these solitons can be computed from $u(x,0)$ using only {\it linear} operations. One of the most important developments in this area has been the generalization of the IST from initial to initial-boundary value (IBV) problems formulated on the half-line. It can be shown that again $u(x,t)$ asymptotes into a collection of solitons, where now the shapes and the speeds of these solitons depend both on $u(x,0)$ and on the given boundary conditions at $x = 0$. A major complication of IBV problems is that the computation of the shapes and speeds of the solitons involves the solution of a {\it nonlinear} Volterra integral equation. However, for a certain class of boundary conditions, called linearizable, this complication can be bypassed and the relevant computation is as effective as in the case of the problem on the line. Here, after reviewing the general theory for KdV, we analyze three different types of linearizable boundary conditions. For these cases, the initial conditions are: (a) restrictions of one and two soliton solutions at $t = 0$; (b) profiles of certain exponential type; (c) box-shaped profiles. For each of these cases, by computing explicitly the shapes and the speeds of the asymptotic solitons, we elucidate the influence of the boundary.

preprint2010arXiv

On a novel integrable generalization of the sine-Gordon equation

We consider an integrable generalization of the sine-Gordon (sG) equation that was earlier derived by one of the authors using bi-Hamiltonian methods. This equation is related to the sG equation in the same way that the Camassa-Holm equation is related to the KdV equation. In this paper we: (a) Derive a Lax pair. (b) Use the Lax pair to solve the initial value problem on the line. (c) Analyze solitons. (d) Show that the generalized sG and sG equations are related by a Liouville transformation. (e) Derive conservation laws. (f) Analyze traveling-wave solutions.