Researcher profile

Daniel Peterseim

Daniel Peterseim contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

8 published item(s)

preprint2022arXiv

Energy-adaptive Riemannian optimization on the Stiefel manifold

This paper addresses the numerical solution of nonlinear eigenvector problems such as the Gross-Pitaevskii and Kohn-Sham equation arising in computational physics and chemistry. These problems characterize critical points of energy minimization problems on the infinite-dimensional Stiefel manifold. To efficiently compute minimizers, we propose a novel Riemannian gradient descent method induced by an energy-adaptive metric. Quantified convergence of the methods is established under suitable assumptions on the underlying problem. A non-monotone line search and the inexact evaluation of Riemannian gradients substantially improve the overall efficiency of the method. Numerical experiments illustrate the performance of the method and demonstrates its competitiveness with well-established schemes.

preprint2022arXiv

Operator Compression with Deep Neural Networks

This paper studies the compression of partial differential operators using neural networks. We consider a family of operators, parameterized by a potentially high-dimensional space of coefficients that may vary on a large range of scales. Based on existing methods that compress such a multiscale operator to a finite-dimensional sparse surrogate model on a given target scale, we propose to directly approximate the coefficient-to-surrogate map with a neural network. We emulate local assembly structures of the surrogates and thus only require a moderately sized network that can be trained efficiently in an offline phase. This enables large compression ratios and the online computation of a surrogate based on simple forward passes through the network is substantially accelerated compared to classical numerical upscaling approaches. We apply the abstract framework to a family of prototypical second-order elliptic heterogeneous diffusion operators as a demonstrating example.

preprint2022arXiv

Super-localized orthogonal decomposition for convection-dominated diffusion problems

This paper presents a multi-scale method for convection-dominated diffusion problems in the regime of large Péclet numbers. The application of the solution operator to piecewise constant right-hand sides on some arbitrary coarse mesh defines a finite-dimensional coarse ansatz space with favorable approximation properties. For some relevant error measures, including the $L^2$-norm, the Galerkin projection onto this generalized finite element space even yields $\varepsilon$-independent error bounds, $\varepsilon$ being the singular perturbation parameter. By constructing an approximate local basis, the approach becomes a novel multi-scale method in the spirit of the Super-Localized Orthogonal Decomposition (SLOD). The error caused by basis localization can be estimated in an a-posteriori way. In contrast to existing multi-scale methods, numerical experiments indicate $\varepsilon$-independent convergence without preasymptotic effects even in the under-resolved regime of large mesh Péclet numbers.

preprint2020arXiv

Computational high frequency scattering from high contrast heterogeneous media

This article considers the computational (acoustic) wave propagation in strongly heterogeneous structures beyond the assumption of periodicity. A high contrast between the constituents of microstructured multiphase materials can lead to unusual wave scattering and absorption, which are interesting and relevant from a physical viewpoint, for instance, in the case of crystals with defects. We present a computational multiscale method in the spirit of the Localized Orthogonal Decomposition and provide its rigorous a priori error analysis for two-phase diffusion coefficients that vary between $1$ and very small values. Special attention is paid to the extreme regimes of high frequency, high contrast, and their previously unexplored coexistence. A series of numerical experiments confirms the theoretical results and demonstrates the ability of the multiscale approach to efficiently capture relevant physical phenomena.

preprint2020arXiv

Quantitative Anderson localization of Schrödinger eigenstates under disorder potentials

This paper concerns spectral properties of linear Schrödinger operators under oscillatory high-amplitude potentials on bounded domains. Depending on the degree of disorder, we prove the existence of spectral gaps amongst the lowermost eigenvalues and the emergence of exponentially localized states. We quantify the rate of decay in terms of geometric parameters that characterize the potential. The proofs are based on the convergence theory of iterative solvers for eigenvalue problems and their optimal local preconditioning by domain decomposition.

preprint2020arXiv

Reconstruction of quasi-local numerical effective models from low-resolution measurements

We consider the inverse problem of reconstructing an effective model for a prototypical diffusion process in strongly heterogeneous media based on coarse measurements. The approach is motivated by quasi-local numerical effective forward models that are provably reliable beyond periodicity assumptions and scale separation. The goal of this work is to show that an identification of the matrix representation related to these effective models is possible. On the one hand, this provides a reasonable surrogate in cases where a direct reconstruction is unfeasible due to a mismatch between the coarse data scale and the microscopic quantities to be reconstructed. On the other hand, the approach allows us to investigate the requirement for a certain non-locality in the context of numerical homogenization. Algorithmic aspects of the inversion procedure and its performance are illustrated in a series of numerical experiments.

preprint2020arXiv

Sobolev gradient flow for the Gross-Pitaevskii eigenvalue problem: global convergence and computational efficiency

We propose a new normalized Sobolev gradient flow for the Gross-Pitaevskii eigenvalue problem based on an energy inner product that depends on time through the density of the flow itself. The gradient flow is well-defined and converges to an eigenfunction. For ground states we can quantify the convergence speed as exponentially fast where the rate depends on spectral gaps of a linearized operator. The forward Euler time discretization of the flow yields a numerical method which generalizes the inverse iteration for the nonlinear eigenvalue problem. For sufficiently small time steps, the method reduces the energy in every step and converges globally in $H^1$ to an eigenfunction. In particular, for any nonnegative starting value, the ground state is obtained. A series of numerical experiments demonstrates the computational efficiency of the method and its competitiveness with established discretizations arising from other gradient flows for this problem.

preprint2020arXiv

Sparse Compression of Expected Solution Operators

We show that the expected solution operator of prototypical linear elliptic partial differential operators with random coefficients is well approximated by a computable sparse matrix. This result is based on a random localized orthogonal multiresolution decomposition of the solution space that allows both the sparse approximate inversion of the random operator represented in this basis as well as its stochastic averaging. The approximate expected solution operator can be interpreted in terms of classical Haar wavelets. When combined with a suitable sampling approach for the expectation, this construction leads to an efficient method for computing a sparse representation of the expected solution operator.