Refine
Year of publication
Language
- English (182)
Keywords
Project
Application Area
- D (182) (remove)
To approximate convolutions which occur in evolution equations with memory terms, a variable-stepsize algorithm is presented for which advancing $N$ steps requires only $O(N\log N)$ operations and $O(\log N)$ active memory, in place of $O(N^2)$ operations and $O(N)$ memory for a direct implementation. A basic feature of the fast algorithm is the reduction, via contour integral representations, to differential equations which are solved numerically with adaptive step sizes. Rather than the kernel itself, its Laplace transform is used in the algorithm. The algorithm is illustrated on three examples: a blow-up example originating from a Schrödinger equation with concentrated nonlinearity, chemical reactions with inhibited diffusion, and viscoelasticity with a fractional order constitutive law.
A new approach to derive transparent boundary conditions (TBCs) for wave, Schrödinger, heat and drift-diffusion equations is presented. It relies on the pole condition and distinguishes between physical reasonable and unreasonable solutions by the location of the singularities of the spatial Laplace transform of the exterior solution. To obtain a numerical algorithm, a Möbius transform is applied to map the Laplace transform onto the unit disc. In the transformed coordinate the solution is expanded into a power series. Finally, equations for the coefficients of the power series are derived. These are coupled to the equation in the interior, and yield transparent boundary conditions.
Numerical results are presented in the last section, showing that the error introduced by the new approximate TBCs decays exponentially in the number of coefficients.
In this review article we discuss different techniques to solve numerically the
time-dependent Schrödinger equation on unbounded domains.
We present in detail the most recent approaches and describe briefly alternative ideas pointing out the relations between these works.
We conclude with several numerical examples from
different application areas to compare the presented techniques. We mainly focus on the one-dimensional problem but also touch upon the situation in two space dimensions and the cubic nonlinear case.
We consider hybrid systems of differential-algebraic equations and present
a general framework for general nonlinear over- and underdetermined hybrid
systems that allows the
analysis of existence and uniqueness and the application of index reduction
methods for hybrid differential-algebraic systems.
A particular difficulty in the numerical simulation of hybrid systems is
(numerical) chattering, i.e., fast oscillations between modes of operations.
A regularization technique using sliding modes allows to regularize the
system behavior in the case of chattering.
Further, we show how chattering behavior during the numerical solution can
be prevented using sliding mode simulation. The advantage of the sliding mode
simulation is illustrated by numerical examples.
The nonlinear Schrödinger equation based on the Taylor approximation of the material dispersion can become invalid for ultrashort and few-cycle optical pulses. Instead, we use a rational fit to the dispersion function such that the resonances are naturally accounted for. This approach allows us to derive a simple non-envelope model for short pulses propagating in one spatial dimension. This model is further investigated numerically and analytically.
A new concept is introduced for the adaptive finite element discretization of partial differential equations that have a sparsely
representable solution. Motivated by recent work on compressed sensing, a recursive mesh refinement procedure is presented that uses linear programming to find a good approximation to the sparse solution on a given refinement level. Then only those parts of the mesh are refined that belong to nonzero expansion coefficients. Error estimates for this procedure are refined and the behavior of the procedure is demonstrated via some simple elliptic model problems.
We consider linear discrete-time descriptor systems, i.e., systems of linear equations of the form $E_{k+1} x^{k+1} = A_k x^k + f^k$ for $k \in \IZ$, where all $E_k$ and $A_k$ are matrices, $f_k$ are vectors and $x_k$ are the vectors of the solution we are looking for. We study the existence and uniqueness of solutions. A strangeness index is defined for such systems. Compared to the continuous-time case, it turns out, that in the discrete-time case it makes a difference, if one has an initial condition and one wants a solution in the future or if one has an initial condition and one wants a solution into the past and future at the same time.
Studying high-dimensional Hamiltonian systems with microstructure, it is an important and challenging problem to identify reduced macroscopic models that describe some effective dynamics on large spatial and temporal scales. This paper concerns the question how reasonable macroscopic Lagrangian and Hamiltonian structures can by derived from the microscopic system. In the first part we develop a general approach to this problem by considering non-canonical Hamiltonian structures on the tangent bundle. This approach can be applied to all Hamiltonian lattices (or Hamiltonian PDEs) and involves three building blocks: (i) the embedding of the microscopic system, (ii) an invertible two-scale transformation that encodes the underlying scaling of space and time, (iii) an elementary model reduction that is based on a Principle of Consistent Expansions. In the second part we exemplify the reduction approach and derive various reduced PDE models for the atomic chain. The reduced equations are either related to long wave-length motion or describe the macroscopic modulation of an oscillatory microstructure.
We study both theoretically and experimentally typical operation
regimes of 40 GHz monolithic mode-locked lasers. The underlying Traveling Wave Equation model reveals quantitative agreement for characteristics of the fundamental mode-locking as pulse width and repetition frequency tuning, as well as qualitative agreement with the experiments for other dynamic regimes. Especially the appearance of stable harmonic mode-locking at 80 GHz
has been predicted theoretically and confirmed by measurements.
Furthermore, we derive and apply a simplified Delay-Differential-Equation model
which guides us to a qualitative analysis of bifurcations responsible for the appearance
and the breakup of different mode-locking regimes. Higher harmonics of mode-locking are predicted by this model as well.
We explore the concept of passive-feedback lasers for direct signal
modulation at 40 Gbit/s. Based on numerical simulation and bifurcation
analysis, we explain the main mechanisms in these devices which are
crucial for modulation at high speed. The predicted effects are
demonstrated experimentally by means of correspondingly designed
devices. In particular a significant improvement of the modulation
bandwidth at low injection currents can be demonstrated.