Source author record

A. S. Fokas

A. S. Fokas 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

25works
9topics
4close 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

25 published item(s)

preprint2026arXiv

The Focusing NLS Equation on the Half-Line with Periodic Boundary Conditions

We consider the Dirichlet problem for the focusing NLS equation on the half-line, with given Schwartz initial data and boundary data $q(0,t)$ equal to an exponentially decaying perturbation $u(t)$ of the periodic boundary data $ a e^{2iωt + i ε}$ at $x=0.$ It is known from PDE theory that this problem admits a unique solution (for fixed initial data and fixed $u$). On the other hand, the associated inverse scattering transform formalism involves the Neumann boundary value for $x=0$. Thus the implementation of this formalism requires the understanding of the "Dirichlet-to-Neumann" map which characterises the associated Neumann boundary value. We consider this map in an indirect way: we postulate a certain Riemann-Hilbert problem, on a specified contour but with partially unspecified jump data of some generality, and then prove that the solution of the initial-boundary value problem for the focusing NLS constructed through this Riemann-Hilbert problem satisfies all the required properties: the data $q(x,0)$ are Schwartz and $q(0,t)-a e^{2iωt + i ε}$ is exponentially decaying. More specifically, we focus on the case $-3a^2 < ω< a^2.$ By considering a large class of appropriate scattering data for the t-problem, we provide solutions of the above Dirichlet problem such that the data $q_x(0,t)$ is given by an exponentially decaying perturbation of the function $2iab e^{2iωt + i ε},$ where $ω= a^2-2b^2,~~b>0$.

preprint2022arXiv

Time-periodic linear boundary value problems on a finite interval

We study the large time behaviour of the solution of a linear dispersive PDEs posed on a finite interval, when the prescribed boundary conditions are time periodic. We use the approach pioneered in Fokas & Lenells 2012 for nonlinear integrable PDEs. and then applied to linear problems on the half-line in Fokas & van der Weele 2021, to characterise necessary conditions for the solution of such a problem to be periodic, at least in an asymptotic sense. We then fully describe the periodicity properties of the solution in three important illustrative examples, recovering known results for the second-order cases and establishing new ones for the third order one.

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

A quantitative framework for exploring exit strategies from the COVID-19 lockdown

Following the highly restrictive measures adopted by many countries for combating the current pandemic, the number of individuals infected by SARS-CoV-2 and the associated number of deaths is steadily decreasing. This fact, together with the impossibility of maintaining the lockdown indefinitely, raises the crucial question of whether it is possible to design an exit strategy based on quantitative analysis. Guided by rigorous mathematical results, we show that this is indeed possible: we present a robust numerical algorithm which can compute the cumulative number of deaths that will occur as a result of increasing the number of contacts by a given multiple, using as input only the most reliable of all data available during the lockdown, namely the cumulative number of deaths.

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.

preprint2013arXiv

A generating mechanism for higher order rogue waves

We introduce a mechanism for generating higher order rogue waves (HRWs) of the nonlinear Schrödinger(NLS) equation: the progressive fusion and fission of $n$ degenerate breathers associated with a critical eigenvalue $λ_0$, creates an order $n$ HRW. By adjusting the relative phase of the breathers at the interacting area, it is possible to obtain different types of HRWs. The value $λ_0$ is a zero point of the eigenfunction of the Lax pair of the NLS equation and it corresponds to the limit of the period of the breather tending to infinity. By employing this mechanism we prove two conjectures regarding the total number of peaks, as well as a decomposition rule in the circular pattern of an order $n$ HRW.

preprint2013arXiv

Elliptic boundary value problems in convex polygons with low regularity boundary data via the unified method

We use novel integral representations developed by the second author to prove certain rigorous results concerning elliptic boundary value problems in convex polygons. Central to this approach is the so-called global relation, which is a non-local equation in the Fourier space that relates the known boundary data to the unknown boundary values. Assuming that the global relation is satisfied in the weakest possible sense, i.e. in a distributional sense, we prove there exist solutions to Dirichlet, Neumann and Robin boundary value problems with distributional boundary data. We also show that the analysis of the global relation characterises in a straightforward manner the possible existence of both integrable and non-integrable corner-singularities.

preprint2012arXiv

The Laplace equation for the exterior of the Hankel contour and novel identities for hypergeometric functions

By employing conformal mappings, it is possible to express the solution of certain boundary value problems for the Laplace equation in terms of a single integral involving the given boundary data. We show that such explicit formulae can be used to obtain novel identities for special functions. A convenient tool for deriving this type of identities is the so-called \emph{global relation}, which has appeared recently in a wide range of boundary value problems. As a concrete application, we analyze the Neumann boundary value problem for the Laplace equation in the exterior of the so-called Hankel contour, which is the contour that appears in the definition of both the gamma and the Riemann zeta functions. By utilizing the explicit solution of this problem, we derive a plethora of novel identities involving the hypergeometric function.

preprint2011arXiv

A Nonlocal Formulation of Rotational Water Waves

The classical equations of irrotational water waves have recently been reformulated as a system of two equations, one of which is an explicit non-local equation for the wave height and for the velocity potential evaluated on the free surface. Here, in the two dimensional case, (a) we generalise the relevant formulation to the case of constant vorticity, as well as to the case where the free surface is described by a multi-valued function; (b) in the case of travelling waves we derive an upper bound for the free surface; (c) in the case of constant vorticity we construct a sequence of nearly Hamiltonian systems which provide an approximation in the asymptotic limit of certain physical small parameters. In particular, the explicit dependence of the vorticity on the coefficients of the KdV equation is clarified. Also, in the irrotational case we extend the formalism to n>2 dimensions and analyse rigorously the linear limit of these equations.

preprint2011arXiv

A Transform Method for Evolution PDEs on a Finite Interval

We study initial boundary value problems for linear scalar partial differential equations with constant coefficients, with spatial derivatives of {\em arbitrary order}, posed on the domain $\{t>0, 0<x<L\}$. We first show that by analysing the so-called {\em global relation}, which is an algebraic relation defined in the complex $k$-plane coupling all boundary values of the solution, it is possible to identify how many boundary conditions must be prescribed at each end of the space interval in order for the problem to be well posed. We then show that the solution can be expressed as an integral in the complex $k$-plane. This integral is defined in terms of an $x$-transform of the initial condition and a $t$-transform of the boundary conditions. For particular cases, such as the case of periodic boundary conditions, or the case of boundary value problems for {\em second} order PDEs, the integral can be rewritten as an infinite series. However, there exist initial boundary value problems for which the only representation is an integral which {\em cannot} be written as an infinite series. An example of such a problem is provided by the linearised version of the KdV equation. Thus, contrary to common belief, the solution of many linear initial boundary value problems on a finite interval {\em cannot} be expressed in terms of an infinite series.

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

A Novel Approach to Elastodynamics: I. The Two-Dimensional Case

We introduce a new approach to constructing analytic solutions of the linear PDEs describing elastodynamics. This approach is illustrated for the case of a homogeneous isotropic half-plane body satisfying arbitrary initial conditions and Lamb's boundary conditions. A particular case of this problem, namely the case of homogeneous initial conditions, was first solved by Lamb using the Fourier-Laplace transform. The solution of the general problem can also be expressed in terms of the Fourier transform, but this representation involves transforms of unknown boundary values. This necessitates the formulation and solution of a cumbersome auxiliary problem, which expresses the unknown boundary values in terms of the Laplace transform of the given boundary data. The new approach, which is applicable to arbitrary initial and boundary conditions, bypasses the above auxiliary problem and expresses the solutions directly in terms of the given initial and boundary data.

preprint2010arXiv

A Novel Approach to Elastodynamics: II. The Three-Dimensional Case

A new approach was recently introduced by the authors for constructing analytic solutions of the linear PDEs describing elastodynamics. Here, this approach is applied to the case of a homogeneous isotropic half-space body satisfying arbitrary initial conditions and Lamb's boundary conditions. A particular case of this problem, namely the case of homogeneous initial conditions and normal point load boundary conditions, was first solved by Lamb using the Fourier-Laplace transform. The general problem solved here can also be analysed via the Fourier transform, but in this case, the solution representation involves transforms of \textit{unknown} boundary values; this necessitates the formulation and solution of a cumbersome auxiliary problem, which expresses the unknown boundary values in terms of the Laplace transform of the given boundary data. The new approach, which is applicable to arbitrary initial and boundary conditions, bypasses the above auxiliary problem and expresses the solutions directly in terms of the given initial and boundary conditions.

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 Class of Integrable ODEs Related to the Painlevé Equations

One of the authors has recently introduced the concept of conjugate Hamiltonian systems: the solution of the equation $h=H(p,q,t),$ where $H$ is a given Hamiltonian containing $t$ explicitly, yields the function $t=T(p,q,h)$, which defines a new Hamiltonian system with Hamiltonian $T$ and independent variable $h.$ By employing this construction and by using the fact that the classical Painlevé equations are Hamiltonian systems, it is straightforward to associate with each Painlevé equation two new integrable ODEs. Here, we investigate the conjugate Painlevé II equations. In particular, for these novel integrable ODEs, we present a Lax pair formulation, as well as a class of implicit solutions. We also construct conjugate equations associated with Painlevé I and Painlevé IV equations.

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.

preprint2010arXiv

The Kadomtsev-Petviashvili II Equation on the Half-Plane

The KPII equation is an integrable nonlinear PDE in 2+1 dimensions (two spatial and one temporal), which arises in several physical circumstances, including fluid mechanics where it describes waves in shallow water. It provides a multidimensional generalisation of the renowned KdV equation. In this work, we employ a novel approach recently introduced by one of the authors in connection with the Davey-Stewartson equation \cite{FDS2009}, in order to analyse the initial-boundary value problem for the KPII equation formulated on the half-plane. The analysis makes crucial use of the so-called d-bar formalism, as well as of the so-called global relation. A novel feature of boundary as opposed to initial-value problems in 2+1 is that the d-bar formalism now involves a function in the complex plane which is discontinuous across the real axis.

preprint2009arXiv

A Novel Method of Solution for the Fluid Loaded Plate

We study the Cauchy problem associated with the equations governing a fluid loaded plate formulated on either the line or the half-line. We show that in both cases the problem can be solved by employing the unified approach to boundary value problems introduced by on of the authors in the late 1990s. The problem on the full line was analysed by Crighton et. al. using a combination of Laplace and Fourier transforms. The new approach avoids the technical difficulty of the a priori assumption that the amplitude of the plate is in $L^1_{dt}(R^+)$ and furthermore yields a simpler solution representation which immediately implies the problem is well-posed. For the problem on the half-line, a similar analysis yields a solution representation, but this formula involves two unknown functions. The main difficulty with the half-line problem is the characterisation of these two functions. By employing the so-called global relation, we show that the two functions can be obtained via the solution of a complex valued integral equation of the convolution type. This equation can be solved in closed form using the Laplace transform. By prescribing the initial data $η_0$ to be in $H^3(R^+)$, we show that the solution depends continuously on the initial data, and hence, the problem is well-posed.