Source author record

Achim Schädle

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

math.NA math-ph math.MP

Catalog footprint

Research overview

8works

3topics

4close collaborators

Actions

Is this your work?

Open graph Browse papers

preprint2016arXiv

Fast and oblivious algorithms for dissipative and 2D wave equations

The use of time-domain boundary integral equations has proved very effective and efficient for three dimensional acoustic and electromagnetic wave equations. In even dimensions and when some dissipation is present, time-domain boundary equations contain an infinite memory tail. Due to this, computation for longer times becomes exceedingly expensive. In this paper we show how oblivious quadrature, initially designed for parabolic problems, can be used to significantly reduce both the cost and the memory requirements of computing this tail. We analyse Runge-Kutta based quadrature and conclude the paper with numerical experiments.

preprint2012arXiv

Transparent boundary conditions based on the Pole Condition for time-dependent, two-dimensional problems

The pole condition approach for deriving transparent boundary conditions is extended to the time-dependent, two-dimensional case. Non-physical modes of the solution are identified by the position of poles of the solution's spatial Laplace transform in the complex plane. By requiring the Laplace transform to be analytic on some problem dependent complex half-plane, these modes can be suppressed. The resulting algorithm computes a finite number of coefficients of a series expansion of the Laplace transform, thereby providing an approximation to the exact boundary condition. The resulting error decays super-algebraically with the number of coefficients, so relatively few additional degrees of freedom are sufficient to reduce the error to the level of the discretization error in the interior of the computational domain. The approach shows good results for the Schrödinger and the drift-diffusion equation but, in contrast to the one-dimensional case, exhibits instabilities for the wave and Klein-Gordon equation. Numerical examples are shown that demonstrate the good performance in the former and the instabilities in the latter case.

preprint2011arXiv

High order Curl-conforming Hardy space infinite elements for exterior Maxwell problems

A construction of prismatic Hardy space infinite elements to discretize wave equations on unbounded domains $Ω$ in $H^1_{loc}(Ω)$ , $H_{loc}(curl;Ω)$ and $H_{loc}(div;Ω)$ is presented. As our motivation is to solve Maxwell's equations we take care that these infinite elements fit into the discrete de Rham diagram, i.e. they span discrete spaces, which together with the exterior derivative form an exact sequence. Resonance as well as scattering problems are considered in the examples. Numerical tests indicate super-algebraic convergence in the number of additional unknowns per degree of freedom on the coupling boundary that are required to realize the Dirichlet to Neumann map.

preprint2010arXiv

Hardy space infinite elements for Helmholtz-type problems with unbounded inhomogeneities

This paper introduces a class of approximate transparent boundary conditions for the solution of Helmholtz-type resonance and scattering problems on unbounded domains. The computational domain is assumed to be a polygon. A detailed description of two variants of the Hardy space infinite element method which relays on the pole condition is given. The method can treat waveguide-type inhomogeneities in the domain with non-compact support. The results of the Hardy space infinite element method are compared to a perfectly matched layer method. Numerical experiments indicate that the approximation error of the Hardy space decays exponentially in the number of Hardy space modes.

preprint2006arXiv

Domain Decomposition Method for Maxwell's Equations: Scattering off Periodic Structures

We present a domain decomposition approach for the computation of the electromagnetic field within periodic structures. We use a Schwarz method with transparent boundary conditions at the interfaces of the domains. Transparent boundary conditions are approximated by the perfectly matched layer method (PML). To cope with Wood anomalies appearing in periodic structures an adaptive strategy to determine optimal PML parameters is developed. We focus on the application to typical EUV lithography line masks. Light propagation within the multi-layer stack of the EUV mask is treated analytically. This results in a drastic reduction of the computational costs and allows for the simulation of next generation lithography masks on a standard personal computer.

preprint2005arXiv

A spectral order method for inverting sectorial Laplace transforms

Laplace transforms which admit a holomorphic extension to some sector strictly containing the right half plane and exhibiting a potential behavior are considered. A spectral order, parallelizable method for their numerical inversion is proposed. The method takes into account the available information about the errors arising in the evaluations. Several numerical illustrations are provided.

preprint2005arXiv

Fast and oblivious convolution quadrature

We give an algorithm to compute $N$ steps of a convolution quadrature approximation to a continuous temporal convolution using only $O(N \log N)$ multiplications and $O(\log N)$ active memory. The method does not require evaluations of the convolution kernel, but instead $O(\log N)$ evaluations of its Laplace transform, which is assumed sectorial. The algorithm can be used for the stable numerical solution with quasi-optimal complexity of linear and nonlinear integral and integro-differential equations of convolution type. In a numerical example we apply it to solve a subdiffusion equation with transparent boundary conditions.

preprint2005arXiv

Fast Runge-Kutta approximation of inhomogeneous parabolic equations

The result after $N$ steps of an implicit Runge-Kutta time discretization of an inhomogeneous linear parabolic differential equation is computed, up to accuracy $ε$ , by solving only $O\Big(\log N \log \frac1ε\Big)$ linear systems of equations. We derive, analyse, and numerically illustrate this fast algorithm.

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Achim Schädle

Research overview

Is this your work?

8 publications

Fast and oblivious algorithms for dissipative and 2D wave equations

Transparent boundary conditions based on the Pole Condition for time-dependent, two-dimensional problems

High order Curl-conforming Hardy space infinite elements for exterior Maxwell problems

Hardy space infinite elements for Helmholtz-type problems with unbounded inhomogeneities

Domain Decomposition Method for Maxwell's Equations: Scattering off Periodic Structures

A spectral order method for inverting sectorial Laplace transforms

Fast and oblivious convolution quadrature

Fast Runge-Kutta approximation of inhomogeneous parabolic equations

See the researcher in context

Building this map preview