Refine
Year of publication
Document Type
- ZIB-Report (23)
- Article (1)
- Habilitation (1)
Is part of the Bibliography
- no (25)
Keywords
Institute
- ZIB Allgemein (25) (remove)
In this paper we propose a new finite element realization of the Perfectly Matched Layer method (PML-method). Our approach allows to deal with arbitrary shaped polygonal domains and with certain types of inhomogeneous exterior domains. Among the covered inhomogeneities are open waveguide structures playing an essential role in integrated optics. We give a detailed insight to implementation aspects. Numerical examples show exponential convergence behavior to the exact solution with the thickness of the PML sponge layer.
We establish the relationship between the transparent boundary condition (BPP) of Baskakov and Popov [Wave Motion 14 (1991) 121-128] and Pakpadakis et. al. [J. Acoust. Soc. Am. 92 (1992) 2030-2038] and a second boundary condition (SDY) introduced by Schmidt and Deuflhard [Comp. Math. Appl. 29 (1995) 53-76] and Schmidt and Yevick [J. Compu. Phys. 134 (1997) 96-107], that is explicitly tailored to the form of the underlying numerical propagation scheme. Our analysis demonstrates that if the domain is first discretized in the propagation direction, the SDY expression can be obtained by applying the exact sequence of steps used to derive the BPP procedure. The BPP method is thus an approximate realization of the computationally far simpler and unconditionally stable SDY boundary condition.
We present a general technique for constructing nonlocal transparent boundary conditions for one-dimensional Schrödinger-type equations. Our method supplies boundary conditions for the $\theta$-family of implicit one-step discretizations of Schrödinger's equation in time. The use of Mikusi\'nski's operator approach in time avoids direct and inverse transforms between time and frequency domains and thus implements the boundary conditions in a direct manner.
Adaptive Multigrid Methods for the Vectorial Maxwell Eigenvalue Problem for Optical Waveguide Design
(2000)
This paper has been motivated by the need for a fast robust adaptive multigrid method to solve the vectorial Maxwell eigenvalue problem arising from the design of optical chips. Our nonlinear multigrid methods are based on a previous method for the scalar Helmholtz equation, which must be modified to cope with the null space of the Maxwell operator due to the divergence condition. We present two different approaches. First, we present a multigrid algorithm based on an edge element discretization of time-harmonic Maxwell's equations, including the divergence condition. Second, an explicit elimination of longitudinal magnetic components leads to a nodal discretization known to avoid discrete \emph{spurious modes} also and a vectorial eigenvalue problem, for which we present a multigrid solver. Numerical examples show that the edge element discretization clearly outperforms the nodal element approach.
Transparent Boundary Conditions for Split-Step Pade Approximations of the One-Way Helmholtz Equation
(1999)
In this paper, we generalize the nonlocal discrete transparent boundary condition introduced by Schmidt and Deuflhard {[}Comp. Math. Appl. 29 (1995) 53-76{]} and Schmidt and Yevick {[}J. Comput. Phys. 134 (1997) 96-107{]} to propagation methods based on arbitrary Pad\'e approximations to the two-dimensional one-way Helmholtz equation. Our approach leads to a recursive formula for the coefficients appearing in the nonlocal condition which then yields an unconditionally stable propagation method.
We derive exact discrete nonreflecting boundary conditions for time-harmonic scattering problems modeled by the Helmholtz equation. The main idea is to consider the exterior problem as an initial value problem with initial data given on the boundary of the computational domain. The solution of the exterior problem is obtained via Laplace transformation techniques which supply the boundary conditions in terms of discrete Dirichlet-to-Neumann operators.
The work presents a new approach to the numerical solution of time-harmonic and time-dependent scattering problems. We replace Sommerfeld's radiation condition valid for the Helmholtz equation by a more general concept called pole condition. The pole condition is based on the Laplace transform of the exterior solution and allows a characterization of outgoing waves. Both new insight into the analysis of scattering problems as well as new numerical algorithms are obtained.
This report collects a number of proposals to determine the lowest eigensolutions of the scalar Helmholtz equation. The basic routine of all discussed algorithms is the standard Rayleigh quotient minimization process. The minimization is performed in a direct multilevel manner, and a subspace iteration is used to determine simultaneously a couple of eigensolutions. As smoother the nonlinear Gauß-Seidel, the nonlinear conjugate gradient method and a preconditioned version of this method are compared with respect to their efficiency. The numerical examples are based on realistic 1D and 2D models of integrated optics components.