Search Results

Now showing 1 - 5 of 5
  • Item
    Implementing exact absorbing boundary condition for the linear one-dimensional Schrödinger problem with variable potential by Titchmarsh-Weyl theory
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2009) Ehrhardt, Matthias; Zheng, Chunxiong
    A new approach for simulating the solution of the time-dependent Schrödinger equation with a general variable potential will be proposed. The key idea is to approximate the Titchmarsh-Weyl m-function (exact Dirichlet-to-Neumann operator) by a rational function with respect to a suitable spectral parameter. With the proposed method we can overcome the usual high-frequency restriction for absorbing boundary conditions of general variable potential problems. We end up with a fast computational algorithm for absorbing boundary conditions that are accurate for the full frequency band
  • Item
    Fast numerical methods for waves in periodic media
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2009) Ehrhardt, Matthias; Zheng, Chunxiong
    Periodic media problems widely exist in many modern application areas like semiconductor nanostructures (e.g. quantum dots and nanocrystals), semi-conductor superlattices, photonic crystals (PC) structures, meta materials or Bragg gratings of surface plasmon polariton (SPP) waveguides, etc. Often these application problems are modeled by partial differential equations with periodic coefficients and/or periodic geometries. In order to numerically solve these periodic structure problems efficiently one usually confines the spatial domain to a bounded computational domain (i.e. in a neighborhood of the region of physical interest). Hereby, the usual strategy is to introduce so-called artificial boundaries and impose suitable boundary conditions. For wave-like equations, the ideal boundary conditions should not only lead to w ell-posed problems, but also mimic the perfect absorption of waves traveling out of the computational domain through the artificial boundaries ...
  • Item
    Numerical simulation of waves in periodic structures
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2008) Ehrhardt, Matthias; Han, Houde; Zheng, Chunxiong
    In this work we present a new numerical technique for solving periodic structure problems. This new approach possesses several advantages. First, it allows for a fast evaluation of the Robin-to-Robin operator for periodic array problems. Secondly, this computational method can also be used for bi-periodic structure problems with local defects. Our strategy is an improvement of the recently developed recursive doubling process by Yuan and Lu. In this paper we consider several problems, such as the exterior elliptic problems with strong coercivity, the time-dependent Schrödinger equation and finally the Helmholtz equation with damping.
  • Item
    Evaluation of exact boundary mappings for one-dimensional semiinfinite periodic arrays
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2008) Ehrhardt, Matthias; Sun, Jiguang; Zheng, Chunxiong
    Periodic arrays are structures consisting of geometrically identical subdomains, usually called periodic cells. In this paper, by taking the Helmholtz equation as a model, we consider the definition and evaluation of the exact boundary mappings for general one-dimensional semi-infinite periodic arrays for any real wavenumber. The well-posedness of the Helmholtz equation is established via the limiting absorption principle. An algorithm based on the doubling procedure and extrapolation technique is proposed to derive the exact Sommerfeld-to-Sommerfeld boundary mapping. The advantages of this algorithm are the robustness and simplicity of implementation. But it also suffers from the high computational cost and the resonance wave numbers. To overcome these shortcomings, we propose another algorithm based on a conjecture about the asymptotic behaviour of limiting absorption principle solutions. The price we have to pay is the resolution of two generalized eigenvalue problems, but still the overall computational cost is significantly reduced. Numerical evidences show that this algorithm presents theoretically the same results as the first algorithm. Moreover, some quantitative comparisons between these two algorithms are given.
  • Item
    Exact artificial boundary conditions for problems with period structure
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2008) Ehrhardt, Matthias; Zheng, Chunxiong
    Based on the work of Zheng on the artificial boundary condition for the Schrödinger equation with sinusoidal potentials at infinity, an analytical impedance expression is presented for general second order ODE problems with periodic coefficients and its validity is shown to be strongly supported by numerical evidences. This new expression for the kernel of the Dirichlet-to-Neumann mapping of the artificial boundary conditions is then used for computing the bound states of the Schrödinger operator with periodic potentials at infinity. Other potential applications are associated with the exact artificial boundary conditions for some time-dependent problems with periodic structures. As an example, a two-dimensional hyperbolic equation modeling the TM polarization of the electromagnetic field with a periodic dielectric permittivity is considered.