Source author record

Stefan Sauter

Stefan Sauter appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA Numerical Analysis

Catalog footprint

What is connected

6works

2topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

Solvability of Discrete Helmholtz Equations

We study the unique solvability of the discretized Helmholtz problem with Robin boundary conditions using a conforming Galerkin $hp$-finite element method. Well-posedness of the discrete equations is typically investigated by applying a compact perturbation to the continuous Helmholtz problem so that a "sufficiently rich" discretization results in a "sufficiently small" perturbation of the continuous problem and well-posedness is inherited via Fredholm's alternative. The qualitative notion "sufficiently rich", however, involves unknown constants and is only of asymptotic nature. Our paper is focussed on a fully discrete approach by mimicking the tools for proving well-posedness of the continuous problem directly on the discrete level. In this way, a computable criterion is derived which certifies discrete well-posedness without relying on an asymptotic perturbation argument. By using this novel approach we obtain a) new stability results for the $hp$-FEM for the Helmholtz problem b) examples for meshes such that the discretization becomes unstable (stiffness matrix is singular), and c) a simple checking Algorithm MOTZ "marching-of-the-zeros" which guarantees in an a posteriori way that a given mesh is certified for a stable Helmholtz discretization.

preprint2022arXiv

Stabilized leapfrog based local time-stepping method for the wave equation

Local time-stepping methods permit to overcome the severe stability constraint on explicit methods caused by local mesh refinement without sacrificing explicitness. In \cite{DiazGrote09}, a leapfrog based explicit local time-stepping (LF-LTS) method was proposed for the time integration of second-order wave equations. Recently, optimal convergence rates were proved for a conforming FEM discretization, albeit under a CFL stability condition where the global time-step, $Δt$, depends on the smallest elements in the mesh \cite{grote_sauter_1}. In general one cannot improve upon that stability constraint, as the LF-LTS method may become unstable at certain discrete values of $Δt$. To remove those critical values of $Δt$, we apply a slight modification (as in recent work on LF-Chebyshev methods \cite{CarHocStu19}) to the original LF-LTS method which nonetheless preserves its desirable properties: it is fully explicit, second-order accurate, satisfies a three-term (leapfrog like) recurrence relation, and conserves the energy. The new stabilized LF-LTS method also yields optimal convergence rates for a standard conforming FE discretization, yet under a CFL condition where $Δt$ no longer depends on the mesh size inside the locally refined region.

preprint2020arXiv

The Heterogeneous Helmholtz Problem with Spherical Symmetry: Green's Operator and Stability Estimates

We study wave propagation phenomena modelled in the frequency domain by the Helmholtz equation in heterogeneous media with focus on media with discontinuous, highly oscillating wave speed. We restrict to problems with spherical symmetry and will derive explicit representations of the Green's operator and stability estimates which are explicit in the frequency and the wave speed.

preprint2020arXiv

Wavenumber-explicit $hp$-FEM analysis for Maxwell's equations with transparent boundary conditions

The time-harmonic Maxwell equations at high wavenumber $k$ are discretized by edge elements of degree $p$ on a mesh of width $h$. For the case of a ball and exact, transparent boundary conditions, we show quasi-optimality of the Galerkin method under the $k$-explicit scale resolution condition that a) $kh/p$ is sufficient small and b) $p/\log k$ is bounded from below.

preprint2015arXiv

A Posteriori Error Estimation of hp-dG Finite Element Methods for Highly Indefinite Helmholtz Problems (extended version)

In this paper, we will consider an $hp$-finite elements discretization of a highly indefinite Helmholtz problem by some dG formulation which is based on the ultra-weak variational formulation by Cessenat and Deprés. We will introduce an a posteriori error estimator and derive reliability and efficiency estimates which are explicit with respect to the wavenumber and the discretization parameters $h$ and $p$. In contrast to the conventional conforming finite element method for indefinite problems, the dG formulation is unconditionally stable and the adaptive discretization process may start from a very coarse initial mesh. Numerical experiments will illustrate the efficiency and robustness of the method.

preprint2014arXiv

Adaptive Time Discretization for Retarded Potentials

In this paper, we will present advanced discretization methods for solving retarded potential integral equations. We employ a $C^{\infty}$-partition of unity method in time and a conventional boundary element method for the spatial discretization. One essential point for the algorithmic realization is the development of an efficient method for approximation the elements of the arising system matrix. We present here an approach which is based on quadrature for (non-analytic) $C^{\infty}$ functions in combination with certain Chebyshev expansions. Furthermore we introduce an a posteriori error estimator for the time discretization which is employed also as an error indicator for adaptive refinement. Numerical experiments show the fast convergence of the proposed quadrature method and the efficiency of the adaptive solution process.