We study infinite soliton trains solutions of nonlinear Schrödinger equations, i.e. solutions behaving as the sum of infinitely many solitary waves at large time. Assuming the composing solitons have sufficiently large relative speeds, we prove the existence and uniqueness of such a soliton train. We also give a new construction of multi-solitons (i.e. finite trains) and prove uniqueness in an exponentially small neighbourhood, and we consider the case of solutions composed of several solitons and kinks (i.e. solutions with a non-zero background at infinity).

Electromagnetic wave propagation in complex dispersive media is governed by the time dependent Maxwell's equations coupled to equations that describe the evolution of the induced macroscopic polarization. We consider “polydispersive” materials represented by distributions of dielectric parameters in a polarization model. The work focuses on a novel computational framework for such problems involving Polynomial Chaos Expansions as a method to improve the modeling accuracy of the Debye model and allow for easy simulation using the Finite Difference Time Domain (FDTD) method. Stability and dispersion analyzes are performed for the approach utilizing the second order Yee scheme in two spatial dimensions.

In this paper we present a fully discrete A-ø finite element method to solve Maxwell’sequations with a nonlinear degenerate boundary condition, which represents ageneralization of the classical Silver-Müller condition for anon-perfect conductor. The relationship between the normal components of theelectric field E and the magnetic field H obeys a power-law nonlinearity of the type H x n = n x (|E x n|α-1E x n) with α ∈ (0,1]. We prove the existence anduniqueness of the solutions of the proposed A-ø scheme and derive the error estimates. Finally, wepresent some numerical experiments to verify the theoretical result.

This study aims to develop a numerical scheme in collocated Cartesian grids to solve the level set equation together with the incompressible two-phase flow equations. A seventh-order accurate upwinding combined compact difference (UCCD7) scheme has been developed for the approximation of the first-order spatial derivative terms shown in the level set equation. Developed scheme has a higher accuracy with a three-point grid stencil to minimize phase error. To preserve the mass of each phase all the time, the temporal derivative term in the level set equation is approximated by the sixth-order accurate symplectic Runge-Kutta (SRK6) scheme. All the simulated results for the dam-break, Rayleigh-Taylor instability, bubble rising, two-bubble merging, and milkcrown problems in two and three dimensions agree well with the available numerical or experimental results.

High-order and conservative phase space direct solvers that preserve the Euler asymptotic limit of the Boltzmann-BGK equation for modelling rarefied gas flows are explored and studied. The approach is based on the conservative discrete ordinate method for velocity space by using Gauss Hermite or Simpsons quadrature rule and conservation of macroscopic properties are enforced on the BGK collision operator. High-order asymptotic-preserving time integration is adopted and the spatial evolution is performed by high-order schemes including a finite difference weighted essentially non-oscillatory method and correction procedure via reconstruction schemes. An artificial viscosity dissipative model is introduced into the Boltzmann-BGK equation when the correction procedure via reconstruction scheme is used. The effects of the discrete velocity conservative property and accuracy of high-order formulations of kinetic schemes based on BGK model methods are provided. Extensive comparative tests with one-dimensional and two-dimensional problems in rarefied gas flows have been carried out to validate and illustrate the schemes presented. Potentially advantageous schemes in terms of stable large time step allowed and higher-order of accuracy are suggested.

In this work, we present a model for an aerosol (air/particle mixture) in the respiratory system. It consists of the incompressible Navier-Stokes equations for the air and the Vlasov equation for the particles in a fixed or moving domain, coupled through a drag force. We propose a discretization of the model, investigate stability properties of the numerical code and sensitivity to parameter perturbation. We also focus on the influence of the aerosol on the airflow.

We consider the semiclassical Schrödinger equation on a compact negatively curved surface. For any sequence of initial data microlocalized on the unit cotangent bundle, we look at the quantum evolution (below the Ehrenfest time) under small perturbations of the Schrödinger equation, and we prove that, in the semiclassical limit, and for typical perturbations, the solutions become equidistributed on the unit cotangent bundle.

More and more experimental evidence demonstrates that the slip boundary condition plays an important role in the study of nano- or micro-scale fluid. We propose a homogenization approach to study the effective slippage problem. We show that the effective slip length obtained by homogenization agrees with the results obtained by the traditional method in the literature for the simplest Stokes flow; then we use our approach to deal with two examples which seem quite hard by other analytical methods. We also include some numerical results to validate our analytical results.

We consider the cubic nonlinear Schrödinger equation posed on the spatial domain $\mathbb{R}\times \mathbb{T}^{d}$. We prove modified scattering and construct modified wave operators for small initial and final data respectively ($1\leqslant d\leqslant 4$). The key novelty comes from the fact that the modified asymptotic dynamics are dictated by the resonant system of this equation, which sustains interesting dynamics when $d\geqslant 2$. As a consequence, we obtain global strong solutions (for $d\geqslant 2$) with infinitely growing high Sobolev norms $H^{s}$.

In this paper, we numerically study the ground and first excited states of the fractional Schrödinger equation in an infinite potential well. Due to the nonlocality of the fractional Laplacian, it is challenging to find the eigenvalues and eigenfunctions of the fractional Schrödinger equation analytically. We first introduce a normalized fractional gradient flow and then discretize it by a quadrature rule method in space and the semi-implicit Euler method in time. Our numerical results suggest that the eigenfunctions of the fractional Schrödinger equation in an infinite potential well differ from those of the standard (non-fractional) Schrödinger equation. We find that the strong nonlocal interactions represented by the fractional Laplacian can lead to a large scattering of particles inside of the potential well. Compared to the ground states, the scattering of particles in the first excited states is larger. Furthermore, boundary layers emerge in the ground states and additionally inner layers exist in the first excited states of the fractional nonlinear Schrödinger equation. Our simulated eigenvalues are consistent with the lower and upper bound estimates in the literature.

We discuss the development, verification, and performance of a GPU accelerated discontinuous Galerkin method for the solutions of two dimensional nonlinear shallow water equations. The shallow water equations are hyperbolic partial differential equations and are widely used in the simulation of tsunami wave propagations. Our algorithms are tailored to take advantage of the single instruction multiple data (SIMD) architecture of graphic processing units. The time integration is accelerated by local time stepping based on a multi-rate Adams-Bashforthscheme. A total variational bounded limiter is adopted for nonlinear stability of the numerical scheme. This limiter is coupled with a mass and momentum conserving positivity preserving limiter for the special treatment of a dry or partially wet element in the triangulation. Accuracy, robustness and performance are demonstrated with the aid of test cases. Furthermore, we developed a unified multi-threading model OCCA. The kernels expressed in OCCA model can be cross-compiled with multi-threading models OpenCL, CUDA, and OpenMP. We compare the performance of the OCCA kernels when cross-compiled with these models.

We define and investigate, via numerical analysis, a one dimensional toy-model of a cloud chamber. An energetic quantum particle, whose initial state is a superposition of two identical wave packets with opposite average momentum, interacts during its evolution and exchanges (small amounts of) energy with an array of localized spins. Triggered by the interaction with the environment, the initial superposition state turns into an incoherent sum of two states describing the following situation: or the particle is going to the left and a large number of spins on the left side changed their states, or the same is happening on the right side. This evolution is reminiscent of what happens in a cloud chamber where a quantum particle, emitted as a spherical wave by a radioactive source, marks its passage inside a supersaturated vapour-chamber in the form of a sequence of small liquid bubbles arranging themselves around a possible classical trajectory of the particle.

We compare six fixed-stepsize fourth-order numerical methods for a number of test problems described by a system of coupled Korteweg–de Vries equations. Particular attention is paid to the ability of these methods to preserve fixed points (solitary waves) and the invariants of the system, and establishing to what extent the conservation of integral invariants is indicative of the solution error for these methods.

Many boundary value problems occur in a natural way while studying fluid flow problems in a channel. The solutions of two such boundary value problems are obtained and analysed in the context of flow problems involving three layers of fluids of different constant densities in a channel, associated with an impermeable bottom that has a small undulation. The top surface of the channel is either bounded by a rigid lid or free to the atmosphere. The fluid in each layer is assumed to be inviscid and incompressible, and the flow is irrotational and two-dimensional. Only waves that are stationary with respect to the bottom profile are considered in this paper. The effect of surface tension is neglected. In the process of obtaining solutions for both the problems, regular perturbation analysis along with a Fourier transform technique is employed to derive the first-order corrections of some important physical quantities. Two types of bottom topography, such as concave and convex, are considered to derive the profiles of the interfaces. We observe that the profiles are oscillatory in nature, representing waves of variable amplitude with distinct wave numbers propagating downstream and with no wave upstream. The observations are presented in tabular and graphical forms.

We introduce a novel solution concept, denoted ${\it\alpha}$-dissipative solutions, that provides a continuous interpolation between conservative and dissipative solutions of the Cauchy problem for the two-component Camassa–Holm system on the line with vanishing asymptotics. All the ${\it\alpha}$-dissipative solutions are global weak solutions of the same equation in Eulerian coordinates, yet they exhibit rather distinct behavior at wave breaking. The solutions are constructed after a transformation into Lagrangian variables, where the solution is carefully modified at wave breaking.

This paper presents a model of a 1D–1D dynamic multi-structure, supporting propagation of a transition wave. It is used to explain the recent phenomenon of the collapse of the San Saba bridge. An analytical model is supplied with illustrative numerical simulations.

We study the propagation of wave packets for a one-dimensional system of two coupled Schrödinger equations with a cubic nonlinearity, in the semiclassical limit. Couplings are induced by the nonlinearity and by the potential, whose eigenvalues present an avoided crossing: at one given point, the gap between them reduces as the semiclassical parameter becomes smaller. For data which are coherent states polarized along an eigenvector of the potential, we prove that when the wave function propagates through the avoided crossing point there are transitions between the eigenspaces at leading order. We analyze the nonlinear effects, which are noticeable away from the crossing point, but see that in a small time interval around this point the nonlinearity’s role is negligible at leading order, and the transition probabilities can be computed with the linear Landau–Zener formula.

This paper is concerned with the Cauchy problem for a nonlinear Schrödinger equation with a harmonic potential and exponential growth nonlinearity in two space dimensions. In the defocusing case, global well-posedness is obtained. In the focusing case, existence of nonglobal solutions is discussed via potential-well arguments.

In this paper, we consider the existence of multi-soliton structures for the nonlinear Klein–Gordon (NLKG) equation in $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}\mathbb{R}^{1+d}$ . We prove that, independently of the unstable character of NLKG solitons, it is possible to construct a $N$ -soliton family of solutions to the NLKG equation, of dimension $2N$ , globally well defined in the energy space $H^1\times L^2$ for all large positive times. The method of proof involves the generalization of previous works on supercritical Nonlinear Schrödinger (NLS) and generalized Korteweg–de Vries (gKdV) equations by Martel, Merle, and the first author [R. Côte, Y. Martel and F. Merle, Rev. Mat. Iberoam. 27 (1) (2011), 273–302] to the wave case, where we replace the unstable mode associated to the linear NLKG operator by two generalized directions that are controlled without appealing to modulation theory. As a byproduct, we generalize the linear theory described in Grillakis, Shatah, and Strauss [J. Funct. Anal. 74 (1) (1987), 160–197] and Duyckaerts and Merle [Int. Math. Res. Pap. IMRP (2008), Art ID rpn002] to the case of boosted solitons, and provide new solutions to be studied using the recent work of Nakanishi and Schlag [Zurich Lectures in Advanced Mathematics, vi+253 pp (European Mathematical Society (EMS), Zürich, 2011)] theory.

In this paper, we obtain the well posedness of the linear stochastic Korteweg–de Vries equation by the Galerkin method, and then establish the Carleman estimate, leading to the unique continuation property (UCP) for the linear stochastic Korteweg–de Vries equation. This UCP cannot be obtained from the classical Holmgren uniqueness theorem.