Dispersive wave runup on non-uniform shores

Historically the finite volume methods have been developed for the numerical integration of conservation laws. In this study we present some recent results on the application of such schemes to dispersive PDEs. Namely, we solve numerically a representative of Boussinesq type equations in view of important applications to the coastal hydrodynamics. Numerical results of the runup of a moderate wave onto a non-uniform beach are presented along with great lines of the employed numerical method (see D. Dutykh et al. (2011) for more details).Comment: 8 pages, 6 figures, 18 references. This preprint is submitted to FVCA6 conference proceedings. Other author papers can be downloaded at http://www.lama.univ-savoie.fr/~dutykh

The unified transform: A spectral collocation method for acoustic scattering

This paper employs the unified transform to present a boundary-based spectral collocation method suitable for solving acoustic scattering problems. The method is suitable for both interior and exterior scattering problems, and may be extended to three dimensions. A number of simple two-dimensional examples are presented to illustrate the versatility of this method, and, upon comparison with other spectral methods, or boundary-based methods the approach presented in this paper is shown to be very competitive

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

Evolution of magnetic fields through cosmological perturbation theory

The origin of galactic and extra-galactic magnetic fields is an unsolved problem in modern cosmology. A possible scenario comes from the idea of these fields emerged from a small field, a seed, which was produced in the early universe (phase transitions, inflation, ...) and it evolves in time. Cosmological perturbation theory offers a natural way to study the evolution of primordial magnetic fields. The dynamics for this field in the cosmological context is described by a cosmic dynamo like equation, through the dynamo term. In this paper we get the perturbed Maxwell's equations and compute the energy momentum tensor to second order in perturbation theory in terms of gauge invariant quantities. Two possible scenarios are discussed, first we consider a FLRW background without magnetic field and we study the perturbation theory introducing the magnetic field as a perturbation. The second scenario, we consider a magnetized FLRW and build up the perturbation theory from this background. We compare the cosmological dynamo like equation in both scenarios

Reaction-diffusion systems with constant diffusivities: conditional symmetries and form-preserving transformations

Q-conditional symmetries (nonclassical symmetries) for a general class of two-component reaction-diffusion systems with constant diffusivities are studied. Using the recently introduced notion of Q-conditional symmetries of the first type (R. Cherniha J. Phys. A: Math. Theor., 2010. vol. 43., 405207), an exhaustive list of reaction-diffusion systems admitting such symmetry is derived. The form-preserving transformations for this class of systems are constructed and it is shown that this list contains only non-equivalent systems. The obtained symmetries permit to reduce the reaction-diffusion systems under study to two-dimensional systems of ordinary differential equations and to find exact solutions. As a non-trivial example, multiparameter families of exact solutions are explicitly constructed for two nonlinear reaction-diffusion systems. A possible interpretation to a biologically motivated model is presented

A hybrid analytical-numerical method for solving advection-dispersion problems on a half-line

This paper employs the unified transform, also known as the Fokas method, to solve the advection-dispersion equation on the half-line. This method combines complex analysis with numerics. Compared to classical approaches used to solve linear partial differential equations (PDEs), the unified transform avoids the solution of ordinary differential equations and, more importantly, constructs an integral representation of the solution in the complex plane which is uniformly convergent at the boundaries. As a consequence, such solutions are well suited for numerical computations. Indeed, the numerical evaluation of the solution requires only the computation of a single contour integral involving an integrand which decays exponentially fast for large values of the integration variable. A novel contribution of this paper, with respect to the solution of linear evolution PDEs in general, and the implementation of the unified transform in particular, is the following: using the advection-dispersion equation as a generic example, it is shown that if the transforms of the given data can be computed analytically, then the unified transform yields a fast and accurate method that converges exponentially with the number of evaluations N yet only has complexity O(N). Furthermore, if the transforms are computed numerically using M evaluations, the unified transform gives rise to a method with complexity O(NM). Results are successfully compared to other existing solutions.F.P.J. de B. acknowledges the support from the National Science Foundation Grant No. EAR-1654009. M.J.C. is supported by EPSRC Grant No. EP/L016516/1. A.S.F. is supported by EPSRC, UK, via the senior fellowship

Spectral theory of some non-selfadjoint linear differential operators

We give a characterisation of the spectral properties of linear differential operators with constant coefficients, acting on functions defined on a bounded interval, and determined by general linear boundary conditions. The boundary conditions may be such that the resulting operator is not selfadjoint. We associate the spectral properties of such an operator $S$ with the properties of the solution of a corresponding boundary value problem for the partial differential equation $\partial_t q \pm iSq=0$. Namely, we are able to establish an explicit correspondence between the properties of the family of eigenfunctions of the operator, and in particular whether this family is a basis, and the existence and properties of the unique solution of the associated boundary value problem. When such a unique solution exists, we consider its representation as a complex contour integral that is obtained using a transform method recently proposed by Fokas and one of the authors. The analyticity properties of the integrand in this representation are crucial for studying the spectral theory of the associated operator.Comment: 1 figur

Zero-dispersion limit for integrable equations on the half-line with linearisable data

We study the zero-dispersion limit for certain initial boundary value problems for the defocusing nonlinear Schrödinger (NLS) equationand for the Korteweg-de Vries (KdV)equation with dominant surface tension. These problems are formulated on the half-line and they involve linearisable boundaryconditions.Peer Reviewe

Gross-Neveu Models, Nonlinear Dirac Equations, Surfaces and Strings

Recent studies of the thermodynamic phase diagrams of the Gross-Neveu model (GN2), and its chiral cousin, the NJL2 model, have shown that there are phases with inhomogeneous crystalline condensates. These (static) condensates can be found analytically because the relevant Hartree-Fock and gap equations can be reduced to the nonlinear Schr\"odinger equation, whose deformations are governed by the mKdV and AKNS integrable hierarchies, respectively. Recently, Thies et al have shown that time-dependent Hartree-Fock solutions describing baryon scattering in the massless GN2 model satisfy the Sinh-Gordon equation, and can be mapped directly to classical string solutions in AdS3. Here we propose a geometric perspective for this result, based on the generalized Weierstrass spinor representation for the embedding of 2d surfaces into 3d spaces, which explains why these well-known integrable systems underlie these various Gross-Neveu gap equations, and why there should be a connection to classical string theory solutions. This geometric viewpoint may be useful for higher dimensional models, where the relevant integrable hierarchies include the Davey-Stewartson and Novikov-Veselov systems.Comment: 27 pages, 1 figur

The implementation of the unified transform to the nonlinear Schrödinger equation with periodic initial conditions

Funder: Göran Gustafssons Stiftelse för Naturvetenskaplig och Medicinsk Forskning; doi: http://dx.doi.org/10.13039/501100003426Funder: Ruth and Nils-Erik Stenbäck FoundationFunder: Engineering and Physical Sciences Research Council; doi: http://dx.doi.org/10.13039/501100000266AbstractThe unified transform method (UTM) provides a novel approach to the analysis of initial boundary value problems for linear as well as for a particular class of nonlinear partial differential equations called integrable. If the latter equations are formulated in two dimensions (either one space and one time, or two space dimensions), the UTM expresses the solution in terms of a matrix Riemann–Hilbert (RH) problem with explicit dependence on the independent variables. For nonlinear integrable evolution equations, such as the celebrated nonlinear Schrödinger (NLS) equation, the associated jump matrices are computed in terms of the initial conditions and all boundary values. The unknown boundary values are characterized in terms of the initial datum and the given boundary conditions via the analysis of the so-called global relation. In general, this analysis involves the solution of certain nonlinear equations. In certain cases, called linearizable, it is possible to bypass this nonlinear step. In these cases, the UTM solves the given initial boundary value problem with the same level of efficiency as the well-known inverse scattering transform solves the initial value problem on the infinite line. We show here that the initial boundary value problem on a finite interval with x-periodic boundary conditions (which can alternatively be viewed as the initial value problem on a circle) belongs to the linearizable class. Indeed, by employing certain transformations of the associated RH problem and by using the global relation, the relevant jump matrices can be expressed explicitly in terms of the so-called scattering data, which are computed in terms of the initial datum. Details are given for NLS, but similar considerations are valid for other well-known integrable evolution equations, including the Korteweg–de Vries (KdV) and modified KdV equations.</jats:p