Initial boundary value problems for the nonlinear Schrödinger equation

A new spectral method for solving initial boundary value problems for linear and integrable nonlinear partial differential equations in two independent variables is applied to the nonlinear Schrödinger equation and to its linearized version in the domain {x≥l(t), t≥0}. We show that there exist two cases: (a) if l″(t)<0, then the solution of the linear or nonlinear equations can be obtained by solving the respective scalar or matrix Riemann-Hilbert problem, which is defined on a time-dependent contour; (b) if l″(t)>0, then the Riemann-Hilbert problem is replaced by a respective scalar or matrix problem on a time-independent domain. In both cases, the solution is expressed in a spectrally decomposed form.

Lie symmetry analysis and reductions of a two-dimensional integrable generalization of the Camassa-Holm equation

In this Letter we investigate Lie symmetries of a (2 + 1)-dimensional integrable generalization of the Camassa-Holm (CH) equation. Through the similarity reductions we obtain four different (1 + 1)-dimensional systems of partial differential equations in which one of them turns out to be a (1 + 1)-dimensional CH equation. We establish their integrability by providing the Lax pair for all of them. Further, we present a brief analysis for some types of particular solutions which include the cuspon, peakon and soliton solutions for the two-dimensional generalization of the CH equation. (C) 2000 Published by Elsevier B.V. B.V.

Symmetry analysis of an integrable reaction-diffusion equation

In this paper, we investigate the invariance and integrability properties of an integrable two-component reaction-diffusion equation. We perform Painleve analysis for both the reaction-diffusion equation modelled by a coupled nonlinear partial differential equations and its general similarity reduced ordinary differential equation and confirm its integrability. Further, we perform Lie symmetry analysis for this model. Interestingly our investigations reveals a rich variety of particular solutions, which have not been reported in the literature, for this model. (C) 2000 Elsevier B.V. Ltd. All rights reserved.

On the integrable perturbations of the Camassa-Holm equation

We present an investigation of the nonlinear partial differential equations (PDE) which are asymptotically representable as a linear combination of the equations from the Camassa-Holm hierarchy. For this purpose we use the infinitesimal transformations of dependent and independent variables of the original PDE. This approach is helpful for the analysis of the systems of the PDE which can be asymptotically represented as the evolution equations of polynomial structure. © 2000 American Institute of Physics.

Lamé differential equations and electrostatics

The problem of existence and uniqueness of polynomial solutions of the Lamé differential equation A(x)y″ + 2B(x)y′ + C(x)y = 0, where A(x),B(x) and C(x) are polynomials of degree p + 1,p and p - 1, is under discussion. We concentrate on the case when A(x) has only real zeros aj and, in contrast to a classical result of Heine and Stieltjes which concerns the case of positive coefficients rj in the partial fraction decomposition B(x)/A(x) = ∑j p=0 rj/(x - aj), we allow the presence of both positive and negative coefficients rj. The corresponding electrostatic interpretation of the zeros of the solution y(x) as points of equilibrium in an electrostatic field generated by charges rj at aj is given. As an application we prove that the zeros of the Gegenbauer-Laurent polynomials are the points of unique equilibrium in a field generated by two positive and two negative charges. © 2000 American Mathematical Society.

On partial derivative-Neumann's problem in generalized functions framework

In this paper, we show that the equation delta u/delta (z) over bar + Gu = f, where the elements involved are in generalized functions context, has a local solution in the generalized functions context.

Virasoro Action on Imaginary Verma Modules and the Operator Form of the KZ-Equation

We define the Virasoro algebra action on imaginary Verma modules for affine and construct an analogue of the Knizhnik-Zamolodchikov equation in the operator form. Both these results are based on a realization of imaginary Verma modules in terms of sums of partial differential operators.

A uniqueness theorem for the determination of sources in the Germain–Lagrange plate equation

We prove a uniqueness result related to the Germain–Lagrange dynamic plate differential equation. We consider the equation {∂2u∂t2+△2u=g⊗f,in ]0,+∞)×R2,u(0)=0,∂u∂t(0)=0, where uu stands for the transverse displacement, ff is a distribution compactly supported in space, and g∈Lloc1([0,+∞)) is a function of time such that g(0)≠0g(0)≠0 and there is a T0>0T0>0 such that g∈C1[0,T0[g∈C1[0,T0[. We prove that the knowledge of uu over an arbitrary open set of the plate for any interval of time ]0,T[]0,T[, 0

Fractional Fokker-Planck-Kolmogorov type Equations and their Associated Stochastic Differential Equations

MSC 2010: 26A33, 35R11, 35R60, 35Q84, 60H10 Dedicated to 80-th anniversary of Professor Rudolf Gorenflo

On Modified Method of Simplest Equation for Obtaining Exact Solutions of Nonlinear PDEs: Case of Elliptic Simplest Equation

2010 Mathematics Subject Classification: 74J30, 34L30.

A Simple Example of Localized Parametric Resonance for the Wave Equation

2000 Mathematics Subject Classification: 35L05, 35P25, 47A40.

Codomain Rigidity of the Dirichlet to Neumann Operator for the Riemannian Wave Equation

We study the Dirichlet to Neumann operator for the Riemannian wave equation on a compact Riemannian manifold. If the Riemannian manifold is modelled as an elastic medium, this operator represents the data available to an observer on the boundary of the manifold when the manifold is set into motion through boundary vibrations. We study the Dirichlet to Neumann operator when vibrations are imposed and data recorded on disjoint sets, a useful setting for applications. We prove that this operator determines the Dirichlet to Neumann operator where sources and observations are on the same set, provided a spectral condition on the Laplace-Beltrami operator for the manifold is satisfied. We prove this by providing an implementable procedure for determining a portion of the Riemannian manifold near the area where sources are applied. Drawing on established results, an immediate corollary is that a compact Riemannian manifold can be reconstructed from the Dirichlet to Neumann operator where sources and observations are on disjoint sets.

Codomain Rigidity of the Dirichlet to Neumann Operator for the Riemannian Wave Equation

We study the Dirichlet to Neumann operator for the Riemannian wave equation on a compact Riemannian manifold. If the Riemannian manifold is modelled as an elastic medium, this operator represents the data available to an observer on the boundary of the manifold when the manifold is set into motion through boundary vibrations. We study the Dirichlet to Neumann operator when vibrations are imposed and data recorded on disjoint sets, a useful setting for applications. We prove that this operator determines the Dirichlet to Neumann operator where sources and observations are on the same set, provided a spectral condition on the Laplace-Beltrami operator for the manifold is satisfied. We prove this by providing an implementable procedure for determining a portion of the Riemannian manifold near the area where sources are applied. Drawing on established results, an immediate corollary is that a compact Riemannian manifold can be reconstructed from the Dirichlet to Neumann operator where sources and observations are on disjoint sets.

Implicit difference approximation for the time fractional diffusion equation

In this paper, we consider a time fractional diffusion equation on a finite domain. The equation is obtained from the standard diffusion equation by replacing the first-order time derivative by a fractional derivative (of order $0<\alpha<1$ ). We propose a computationally effective implicit difference approximation to solve the time fractional diffusion equation. Stability and convergence of the method are discussed. We prove that the implicit difference approximation (IDA) is unconditionally stable, and the IDA is convergent with $O(\tau+h^2)$, where $\tau$ and $h$ are time and space steps, respectively. Some numerical examples are presented to show the application of the present technique.

Krylov subspace methods for approximating functions of symmetric positive definite matrices with applications to applied statistics and anomalous diffusion

Matrix function approximation is a current focus of worldwide interest and finds application in a variety of areas of applied mathematics and statistics. In this thesis we focus on the approximation of A^(-α/2)b, where A ∈ ℝ^(n×n) is a large, sparse symmetric positive definite matrix and b ∈ ℝ^n is a vector. In particular, we will focus on matrix function techniques for sampling from Gaussian Markov random fields in applied statistics and the solution of fractional-in-space partial differential equations. Gaussian Markov random fields (GMRFs) are multivariate normal random variables characterised by a sparse precision (inverse covariance) matrix. GMRFs are popular models in computational spatial statistics as the sparse structure can be exploited, typically through the use of the sparse Cholesky decomposition, to construct fast sampling methods. It is well known, however, that for sufficiently large problems, iterative methods for solving linear systems outperform direct methods. Fractional-in-space partial differential equations arise in models of processes undergoing anomalous diffusion. Unfortunately, as the fractional Laplacian is a non-local operator, numerical methods based on the direct discretisation of these equations typically requires the solution of dense linear systems, which is impractical for fine discretisations. In this thesis, novel applications of Krylov subspace approximations to matrix functions for both of these problems are investigated. Matrix functions arise when sampling from a GMRF by noting that the Cholesky decomposition A = LL^T is, essentially, a `square root' of the precision matrix A. Therefore, we can replace the usual sampling method, which forms x = L^(-T)z, with x = A^(-1/2)z, where z is a vector of independent and identically distributed standard normal random variables. Similarly, the matrix transfer technique can be used to build solutions to the fractional Poisson equation of the form ϕn = A^(-α/2)b, where A is the finite difference approximation to the Laplacian. Hence both applications require the approximation of f(A)b, where f(t) = t^(-α/2) and A is sparse. In this thesis we will compare the Lanczos approximation, the shift-and-invert Lanczos approximation, the extended Krylov subspace method, rational approximations and the restarted Lanczos approximation for approximating matrix functions of this form. A number of new and novel results are presented in this thesis. Firstly, we prove the convergence of the matrix transfer technique for the solution of the fractional Poisson equation and we give conditions by which the finite difference discretisation can be replaced by other methods for discretising the Laplacian. We then investigate a number of methods for approximating matrix functions of the form A^(-α/2)b and investigate stopping criteria for these methods. In particular, we derive a new method for restarting the Lanczos approximation to f(A)b. We then apply these techniques to the problem of sampling from a GMRF and construct a full suite of methods for sampling conditioned on linear constraints and approximating the likelihood. Finally, we consider the problem of sampling from a generalised Matern random field, which combines our techniques for solving fractional-in-space partial differential equations with our method for sampling from GMRFs.