6,458 research outputs found
Information Transmission using the Nonlinear Fourier Transform, Part II: Numerical Methods
In this paper, numerical methods are suggested to compute the discrete and
the continuous spectrum of a signal with respect to the Zakharov-Shabat system,
a Lax operator underlying numerous integrable communication channels including
the nonlinear Schr\"odinger channel, modeling pulse propagation in optical
fibers. These methods are subsequently tested and their ability to estimate the
spectrum are compared against each other. These methods are used to compute the
spectrum of various signals commonly used in the optical fiber communications.
It is found that the layer-peeling and the spectral methods are suitable
schemes to estimate the nonlinear spectra with good accuracy. To illustrate the
structure of the spectrum, the locus of the eigenvalues is determined under
amplitude and phase modulation in a number of examples. It is observed that in
some cases, as signal parameters vary, eigenvalues collide and change their
course of motion. The real axis is typically the place from which new
eigenvalues originate or are absorbed into after traveling a trajectory in the
complex plane.Comment: Minor updates to IEEE Transactions on Information Theory, vol. 60,
no. 7, pp. 4329--4345, July 201
Multi-patch discontinuous Galerkin isogeometric analysis for wave propagation: explicit time-stepping and efficient mass matrix inversion
We present a class of spline finite element methods for time-domain wave
propagation which are particularly amenable to explicit time-stepping. The
proposed methods utilize a discontinuous Galerkin discretization to enforce
continuity of the solution field across geometric patches in a multi-patch
setting, which yields a mass matrix with convenient block diagonal structure.
Over each patch, we show how to accurately and efficiently invert mass matrices
in the presence of curved geometries by using a weight-adjusted approximation
of the mass matrix inverse. This approximation restores a tensor product
structure while retaining provable high order accuracy and semi-discrete energy
stability. We also estimate the maximum stable timestep for spline-based finite
elements and show that the use of spline spaces result in less stringent CFL
restrictions than equivalent piecewise continuous or discontinuous finite
element spaces. Finally, we explore the use of optimal knot vectors based on L2
n-widths. We show how the use of optimal knot vectors can improve both
approximation properties and the maximum stable timestep, and present a simple
heuristic method for approximating optimal knot positions. Numerical
experiments confirm the accuracy and stability of the proposed methods
Error analysis of trigonometric integrators for semilinear wave equations
An error analysis of trigonometric integrators (or exponential integrators)
applied to spatial semi-discretizations of semilinear wave equations with
periodic boundary conditions in one space dimension is given. In particular,
optimal second-order convergence is shown requiring only that the exact
solution is of finite energy. The analysis is uniform in the spatial
discretization parameter. It covers the impulse method which coincides with the
method of Deuflhard and the mollified impulse method of Garc\'ia-Archilla,
Sanz-Serna & Skeel as well as the trigonometric methods proposed by Hairer &
Lubich and by Grimm & Hochbruck. The analysis can also be used to explain the
convergence behaviour of the St\"ormer-Verlet/leapfrog discretization in time.Comment: 25 page
- …