-
Hierarchical Super-Localized Orthogonal Decomposition Method
Authors:
Jose C. Garay,
Hannah Mohr,
Daniel Peterseim,
Christoph Zimmer
Abstract:
We present the construction of a sparse-compressed operator that approximates the solution operator of elliptic PDEs with rough coefficients. To derive the compressed operator, we construct a hierarchical basis of an approximate solution space, with superlocalized basis functions that are quasi-orthogonal across hierarchy levels with respect to the inner product induced by the energy norm. The sup…
▽ More
We present the construction of a sparse-compressed operator that approximates the solution operator of elliptic PDEs with rough coefficients. To derive the compressed operator, we construct a hierarchical basis of an approximate solution space, with superlocalized basis functions that are quasi-orthogonal across hierarchy levels with respect to the inner product induced by the energy norm. The superlocalization is achieved through a novel variant of the Super-Localized Orthogonal Decomposition method that is built upon corrections of basis functions arising from the Localized Orthogonal Decomposition method. The hierarchical basis not only induces a sparse compression of the solution space but also enables an orthogonal multiresolution decomposition of the approximate solution operator, decoupling scales and solution contributions of each level of the hierarchy. With this decomposition, the solution of the PDE reduces to the solution of a set of independent linear systems per level with mesh-independent condition numbers that can be computed simultaneously. We present an accuracy study of the compressed solution operator as well as numerical results illustrating our theoretical findings and beyond, revealing that desired optimal error rates with well-behaved superlocalized basis functions can still be attained even in the challenging case of coefficients with high-contrast channels.
△ Less
Submitted 26 July, 2024;
originally announced July 2024.
-
Positivity preserving finite element method for the Gross-Pitaevskii ground state: discrete uniqueness and global convergence
Authors:
Moritz Hauck,
Yizhou Liang,
Daniel Peterseim
Abstract:
We propose a positivity preserving finite element discretization for the nonlinear Gross-Pitaevskii eigenvalue problem. The method employs mass lumping techniques, which allow to transfer the uniqueness up to sign and positivity properties of the continuous ground state to the discrete setting. We further prove that every non-negative discrete excited state up to sign coincides with the discrete g…
▽ More
We propose a positivity preserving finite element discretization for the nonlinear Gross-Pitaevskii eigenvalue problem. The method employs mass lumping techniques, which allow to transfer the uniqueness up to sign and positivity properties of the continuous ground state to the discrete setting. We further prove that every non-negative discrete excited state up to sign coincides with the discrete ground state. This allows one to identify the limit of fully discretized gradient flows, which are typically used to compute the discrete ground state, and thereby establish their global convergence. Furthermore, we perform a rigorous a priori error analysis of the proposed non-standard finite element discretization, showing optimal orders of convergence for all unknowns. Numerical experiments illustrate the theoretical results of this paper.
△ Less
Submitted 27 May, 2024;
originally announced May 2024.
-
Hierarchical Rank-One Sequence Convexification for the Relaxation of Variational Problems with Microstructures
Authors:
Maximilian Köhler,
Timo Neumeier,
Malte. A. Peter,
Daniel Peterseim,
Daniel Balzani
Abstract:
This paper presents an efficient algorithm for the approximation of the rank-one convex hull in the context of nonlinear solid mechanics. It is based on hierarchical rank-one sequences and simultaneously provides first and second derivative information essential for the calculation of mechanical stresses and the computational minimization of discretized energies. For materials, whose microstructur…
▽ More
This paper presents an efficient algorithm for the approximation of the rank-one convex hull in the context of nonlinear solid mechanics. It is based on hierarchical rank-one sequences and simultaneously provides first and second derivative information essential for the calculation of mechanical stresses and the computational minimization of discretized energies. For materials, whose microstructure can be well approximated in terms of laminates and where each laminate stage achieves energetic optimality with respect to the current stage, the approximate envelope coincides with the rank-one convex envelope. Although the proposed method provides only an upper bound for the rank-one convex hull, a careful examination of the resulting constraints shows a decent applicability in mechanical problems. Various aspects of the algorithm are discussed, including the restoration of rotational invariance, microstructure reconstruction, comparisons with other semi-convexification algorithms, and mesh independency. Overall, this paper demonstrates the efficiency of the algorithm for both, well-established mathematical benchmark problems as well as nonconvex isotropic finite-strain continuum damage models in two and three dimensions. Thereby, for the first time, a feasible concurrent numerical relaxation is established for an incremental, dissipative large-strain model with relevant applications in engineering problems.
△ Less
Submitted 27 May, 2024;
originally announced May 2024.
-
A simple collocation-type approach to numerical stochastic homogenization
Authors:
Moritz Hauck,
Hannah Mohr,
Daniel Peterseim
Abstract:
This paper proposes a novel collocation-type numerical stochastic homogenization method for prototypical stochastic homogenization problems with random coefficient fields of small correlation lengths. The presented method is based on a recently introduced localization technique that enforces a super-exponential decay of the basis functions relative to the underlying coarse mesh, resulting in consi…
▽ More
This paper proposes a novel collocation-type numerical stochastic homogenization method for prototypical stochastic homogenization problems with random coefficient fields of small correlation lengths. The presented method is based on a recently introduced localization technique that enforces a super-exponential decay of the basis functions relative to the underlying coarse mesh, resulting in considerable computational savings during the sampling phase. More generally, the collocation-type structure offers a particularly simple and computationally efficient construction in the stochastic setting with minimized communication between the patches where the basis functions of the method are computed. An error analysis that bridges numerical homogenization and the quantitative theory of stochastic homogenization is performed. In a series of numerical experiments, we study the effect of the correlation length and the discretization parameters on the approximation quality of the method.
△ Less
Submitted 2 April, 2024;
originally announced April 2024.
-
Quantum Realization of the Finite Element Method
Authors:
Matthias Deiml,
Daniel Peterseim
Abstract:
This paper presents a quantum algorithm for the solution of prototypical second-order linear elliptic partial differential equations discretized by $d$-linear finite elements on Cartesian grids of a bounded $d$-dimensional domain. An essential step in the construction is a BPX preconditioner, which transforms the linear system into a sufficiently well-conditioned one, making it amenable to quantum…
▽ More
This paper presents a quantum algorithm for the solution of prototypical second-order linear elliptic partial differential equations discretized by $d$-linear finite elements on Cartesian grids of a bounded $d$-dimensional domain. An essential step in the construction is a BPX preconditioner, which transforms the linear system into a sufficiently well-conditioned one, making it amenable to quantum computation. We provide a constructive proof demonstrating that our quantum algorithm can compute suitable functionals of the solution to a given tolerance $\texttt{tol}$ with a complexity linear in $\texttt{tol}^{-1}$ for a fixed dimension $d$, neglecting logarithmic terms. This complexity is proportional to that of its one-dimensional counterpart and improves previous quantum algorithms by a factor of order $\texttt{tol}^{-2}$. We also detail the design and implementation of a quantum circuit capable of executing our algorithm, and present simulator results that support the quantum feasibility of the finite element method in the near future, paving the way for quantum computing approaches to a wide range of PDE-related challenges.
△ Less
Submitted 28 March, 2024;
originally announced March 2024.
-
Mixed finite elements for the Gross-Pitaevskii eigenvalue problem: a priori error analysis and guaranteed lower energy bound
Authors:
Dietmar Gallistl,
Moritz Hauck,
Yizhou Liang,
Daniel Peterseim
Abstract:
We establish an a priori error analysis for the lowest-order Raviart-Thomas finite element discretisation of the nonlinear Gross-Pitaevskii eigenvalue problem. Optimal convergence rates are obtained for the primal and dual variables as well as for the eigenvalue and energy approximations. In contrast to conformal approaches, which naturally imply upper energy bounds, the proposed mixed discretisat…
▽ More
We establish an a priori error analysis for the lowest-order Raviart-Thomas finite element discretisation of the nonlinear Gross-Pitaevskii eigenvalue problem. Optimal convergence rates are obtained for the primal and dual variables as well as for the eigenvalue and energy approximations. In contrast to conformal approaches, which naturally imply upper energy bounds, the proposed mixed discretisation provides a guaranteed and asymptotically exact lower bound for the ground state energy. The theoretical results are illustrated by a series of numerical experiments.
△ Less
Submitted 9 February, 2024;
originally announced February 2024.
-
Reconstructing microstructures from statistical descriptors using neural cellular automata
Authors:
Paul Seibert,
Alexander Raßloff,
Yichi Zhang,
Karl Kalina,
Paul Reck,
Daniel Peterseim,
Markus Kästner
Abstract:
The problem of generating microstructures of complex materials in silico has been approached from various directions including simulation, Markov, deep learning and descriptor-based approaches. This work presents a hybrid method that is inspired by all four categories and has interesting scalability properties. A neural cellular automaton is trained to evolve microstructures based on local informa…
▽ More
The problem of generating microstructures of complex materials in silico has been approached from various directions including simulation, Markov, deep learning and descriptor-based approaches. This work presents a hybrid method that is inspired by all four categories and has interesting scalability properties. A neural cellular automaton is trained to evolve microstructures based on local information. Unlike most machine learning-based approaches, it does not directly require a data set of reference micrographs, but is trained from statistical microstructure descriptors that can stem from a single reference. This means that the training cost scales only with the complexity of the structure and associated descriptors. Since the size of the reconstructed structures can be set during inference, even extremely large structures can be efficiently generated. Similarly, the method is very efficient if many structures are to be reconstructed from the same descriptor for statistical evaluations. The method is formulated and discussed in detail by means of various numerical experiments, demonstrating its utility and scalability.
△ Less
Submitted 28 September, 2023;
originally announced September 2023.
-
Super-localised wave function approximation of Bose-Einstein condensates
Authors:
Daniel Peterseim,
Johan Wärnegård,
Christoph Zimmer
Abstract:
This paper presents a novel spatial discretisation method for the reliable and efficient simulation of Bose-Einstein condensates modelled by the Gross-Pitaevskii equation and the corresponding nonlinear eigenvector problem. The method combines the high-accuracy properties of numerical homogenisation methods with a novel super-localisation approach for the calculation of the basis functions. A rigo…
▽ More
This paper presents a novel spatial discretisation method for the reliable and efficient simulation of Bose-Einstein condensates modelled by the Gross-Pitaevskii equation and the corresponding nonlinear eigenvector problem. The method combines the high-accuracy properties of numerical homogenisation methods with a novel super-localisation approach for the calculation of the basis functions. A rigorous numerical analysis demonstrates superconvergence of the approach compared to classical polynomial and multiscale finite element methods, even in low regularity regimes. Numerical tests reveal the method's competitiveness with spectral methods, particularly in capturing critical physical effects in extreme conditions, such as vortex lattice formation in fast-rotating potential traps. The method's potential is further highlighted through a dynamic simulation of a phase transition from Mott insulator to Bose-Einstein condensate, emphasising its capability for reliable exploration of physical phenomena.
△ Less
Submitted 21 September, 2023;
originally announced September 2023.
-
Computational polyconvexification of isotropic functions
Authors:
Timo Neumeier,
Malte A. Peter,
Daniel Peterseim,
David Wiedemann
Abstract:
Based on the characterization of the polyconvex envelope of isotropic functions by their signed singular value representations, we propose a simple algorithm for the numerical approximation of the polyconvex envelope. Instead of operating on the $d^2$-dimensional space of matrices, the algorithm requires only the computation of the convex envelope of a function on a $d$-dimensional manifold, which…
▽ More
Based on the characterization of the polyconvex envelope of isotropic functions by their signed singular value representations, we propose a simple algorithm for the numerical approximation of the polyconvex envelope. Instead of operating on the $d^2$-dimensional space of matrices, the algorithm requires only the computation of the convex envelope of a function on a $d$-dimensional manifold, which is easily realized by standard algorithms. The significant speedup associated with the dimensional reduction from $d^2$ to $d$ is demonstrated in a series of numerical experiments.
△ Less
Submitted 28 July, 2023;
originally announced July 2023.
-
Riemannian Newton methods for energy minimization problems of Kohn-Sham type
Authors:
R. Altmann,
D. Peterseim,
T. Stykel
Abstract:
This paper is devoted to the numerical solution of constrained energy minimization problems arising in computational physics and chemistry such as the Gross-Pitaevskii and Kohn-Sham models. In particular, we introduce the Riemannian Newton methods on the infinite-dimensional Stiefel and Grassmann manifolds. We study the geometry of these two manifolds, its impact on the Newton algorithms, and pres…
▽ More
This paper is devoted to the numerical solution of constrained energy minimization problems arising in computational physics and chemistry such as the Gross-Pitaevskii and Kohn-Sham models. In particular, we introduce the Riemannian Newton methods on the infinite-dimensional Stiefel and Grassmann manifolds. We study the geometry of these two manifolds, its impact on the Newton algorithms, and present expressions of the Riemannian Hessians in the infinite-dimensional setting, which are suitable for variational spatial discretizations. A series of numerical experiments illustrates the performance of the methods and demonstrates its supremacy compared to other well-established schemes such as the self-consistent field iteration and gradient descent schemes.
△ Less
Submitted 27 June, 2024; v1 submitted 25 July, 2023;
originally announced July 2023.
-
A reduced basis super-localized orthogonal decomposition for reaction-convection-diffusion problems
Authors:
Francesca Bonizzoni,
Moritz Hauck,
Daniel Peterseim
Abstract:
This paper presents a method for the numerical treatment of reaction-convection-diffusion problems with parameter-dependent coefficients that are arbitrary rough and possibly varying at a very fine scale. The presented technique combines the reduced basis (RB) framework with the recently proposed super-localized orthogonal decomposition (SLOD). More specifically, the RB is used for accelerating th…
▽ More
This paper presents a method for the numerical treatment of reaction-convection-diffusion problems with parameter-dependent coefficients that are arbitrary rough and possibly varying at a very fine scale. The presented technique combines the reduced basis (RB) framework with the recently proposed super-localized orthogonal decomposition (SLOD). More specifically, the RB is used for accelerating the typically costly SLOD basis computation, while the SLOD is employed for an efficient compression of the problem's solution operator requiring coarse solves only. The combined advantages of both methods allow one to tackle the challenges arising from parametric heterogeneous coefficients. Given a value of the parameter vector, the method outputs a corresponding compressed solution operator which can be used to efficiently treat multiple, possibly non-affine, right-hand sides at the same time, requiring only one coarse solve per right-hand side.
△ Less
Submitted 28 November, 2022;
originally announced November 2022.
-
Multidimensional rank-one convexification of incremental damage models at finite strains
Authors:
Daniel Balzani,
Maximilian Köhler,
Timo Neumeier,
Malte A. Peter,
Daniel Peterseim
Abstract:
This paper presents computationally feasible rank-one relaxation algorithms for the efficient simulation of a time-incremental damage model with nonconvex incremental stress potentials in multiple spatial dimensions. While the standard model suffers from numerical issues due to the lack of convexity, the relaxation by rank-one convexification prevents non-existence of minimizers and mesh dependenc…
▽ More
This paper presents computationally feasible rank-one relaxation algorithms for the efficient simulation of a time-incremental damage model with nonconvex incremental stress potentials in multiple spatial dimensions. While the standard model suffers from numerical issues due to the lack of convexity, the relaxation by rank-one convexification prevents non-existence of minimizers and mesh dependence of the solutions of finite element discretizations. By the combination, modification and parallelization of the underlying convexification algorithms, the novel approach becomes computationally feasible. A descent method and a Newton scheme enhanced by step-size control prevent stability issues related to local minima in the energy landscape and the computation of derivatives. Numerical techniques for the construction of continuous derivatives of the approximated rank-one convex envelope are discussed. A series of numerical experiments demonstrates the ability of the computationally relaxed model to capture softening effects and the mesh independence of the computed approximations. An interpretation in terms of microstructural damage evolution is given, based on the rank-one lamination process.
△ Less
Submitted 9 February, 2023; v1 submitted 24 November, 2022;
originally announced November 2022.
-
Computational multiscale methods for nondivergence-form elliptic partial differential equations
Authors:
Philip Freese,
Dietmar Gallistl,
Daniel Peterseim,
Timo Sprekeler
Abstract:
This paper proposes novel computational multiscale methods for linear second-order elliptic partial differential equations in nondivergence-form with heterogeneous coefficients satisfying a Cordes condition. The construction follows the methodology of localized orthogonal decomposition (LOD) and provides operator-adapted coarse spaces by solving localized cell problems on a fine scale in the spiri…
▽ More
This paper proposes novel computational multiscale methods for linear second-order elliptic partial differential equations in nondivergence-form with heterogeneous coefficients satisfying a Cordes condition. The construction follows the methodology of localized orthogonal decomposition (LOD) and provides operator-adapted coarse spaces by solving localized cell problems on a fine scale in the spirit of numerical homogenization. The degrees of freedom of the coarse spaces are related to nonconforming and mixed finite element methods for homogeneous problems. The rigorous error analysis of one exemplary approach shows that the favorable properties of the LOD methodology known from divergence-form PDEs, i.e., its applicability and accuracy beyond scale separation and periodicity, remain valid for problems in nondivergence-form.
△ Less
Submitted 24 June, 2023; v1 submitted 24 November, 2022;
originally announced November 2022.
-
A Super-Localized Generalized Finite Element Method
Authors:
Philip Freese,
Moritz Hauck,
Tim Keil,
Daniel Peterseim
Abstract:
This paper presents a novel multi-scale method for elliptic partial differential equations with arbitrarily rough coefficients. In the spirit of numerical homogenization, the method constructs problem-adapted ansatz spaces with uniform algebraic approximation rates. Localized basis functions with the same super-exponential localization properties as the recently proposed Super-Localized Orthogonal…
▽ More
This paper presents a novel multi-scale method for elliptic partial differential equations with arbitrarily rough coefficients. In the spirit of numerical homogenization, the method constructs problem-adapted ansatz spaces with uniform algebraic approximation rates. Localized basis functions with the same super-exponential localization properties as the recently proposed Super-Localized Orthogonal Decomposition enable an efficient implementation. The method's basis stability is enforced using a partition of unity approach. A natural extension to higher order is presented, resulting in higher approximation rates and enhanced localization properties. We perform a rigorous a priori and a posteriori error analysis and confirm our theoretical findings in a series of numerical experiments. In particular, we demonstrate the method's applicability for challenging high-contrast channeled coefficients.
△ Less
Submitted 17 November, 2022;
originally announced November 2022.
-
Neural network approximation of coarse-scale surrogates in numerical homogenization
Authors:
Fabian Kröpfl,
Roland Maier,
Daniel Peterseim
Abstract:
Coarse-scale surrogate models in the context of numerical homogenization of linear elliptic problems with arbitrary rough diffusion coefficients rely on the efficient solution of fine-scale sub-problems on local subdomains whose solutions are then employed to deduce appropriate coarse contributions to the surrogate model. However, in the absence of periodicity and scale separation, the reliability…
▽ More
Coarse-scale surrogate models in the context of numerical homogenization of linear elliptic problems with arbitrary rough diffusion coefficients rely on the efficient solution of fine-scale sub-problems on local subdomains whose solutions are then employed to deduce appropriate coarse contributions to the surrogate model. However, in the absence of periodicity and scale separation, the reliability of such models requires the local subdomains to cover the whole domain which may result in high offline costs, in particular for parameter-dependent and stochastic problems. This paper justifies the use of neural networks for the approximation of coarse-scale surrogate models by analyzing their approximation properties. For a prototypical and representative numerical homogenization technique, the Localized Orthogonal Decomposition method, we show that one single neural network is sufficient to approximate the coarse contributions of all occurring coefficient-dependent local sub-problems for a non-trivial class of diffusion coefficients up to arbitrary accuracy. We present rigorous upper bounds on the depth and number of non-zero parameters for such a network to achieve a given accuracy. Further, we analyze the overall error of the resulting neural network enhanced numerical homogenization surrogate model.
△ Less
Submitted 6 September, 2022;
originally announced September 2022.
-
Super-localized orthogonal decomposition for convection-dominated diffusion problems
Authors:
Francesca Bonizzoni,
Philip Freese,
Daniel Peterseim
Abstract:
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 pro…
▽ More
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.
△ Less
Submitted 4 June, 2022;
originally announced June 2022.
-
Super-localized Orthogonal Decomposition for high-frequency Helmholtz problems
Authors:
Philip Freese,
Moritz Hauck,
Daniel Peterseim
Abstract:
We propose a novel variant of the Localized Orthogonal Decomposition (LOD) method for time-harmonic scattering problems of Helmholtz type with high wavenumber $κ$. On a coarse mesh of width $H$, the proposed method identifies local finite element source terms that yield rapidly decaying responses under the solution operator. They can be constructed to high accuracy from independent local snapshot…
▽ More
We propose a novel variant of the Localized Orthogonal Decomposition (LOD) method for time-harmonic scattering problems of Helmholtz type with high wavenumber $κ$. On a coarse mesh of width $H$, the proposed method identifies local finite element source terms that yield rapidly decaying responses under the solution operator. They can be constructed to high accuracy from independent local snapshot solutions on patches of width $\ell H$ and are used as problem-adapted basis functions in the method. In contrast to the classical LOD and other state-of-the-art multi-scale methods, the localization error decays super-exponentially as the oversampling parameter $\ell$ is increased. This implies that optimal convergence is observed under the substantially relaxed oversampling condition $\ell \gtrsim (\log \tfracκ{H})^{(d-1)/d}$ with $d$ denoting the spatial dimension. Numerical experiments demonstrate the significantly improved offline and online performance of the method also in the case of heterogeneous media and perfectly matched layers.
△ Less
Submitted 21 December, 2021;
originally announced December 2021.
-
Energy-adaptive Riemannian optimization on the Stiefel manifold
Authors:
Robert Altmann,
Daniel Peterseim,
Tatjana Stykel
Abstract:
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…
▽ More
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.
△ Less
Submitted 16 April, 2022; v1 submitted 22 August, 2021;
originally announced August 2021.
-
Super-localization of elliptic multiscale problems
Authors:
Moritz Hauck,
Daniel Peterseim
Abstract:
Numerical homogenization aims to efficiently and accurately approximate the solution space of an elliptic partial differential operator with arbitrarily rough coefficients in a $d$-dimensional domain. The application of the inverse operator to some standard finite element space defines an approximation space with uniform algebraic approximation rates with respect to the mesh size parameter $H$. Th…
▽ More
Numerical homogenization aims to efficiently and accurately approximate the solution space of an elliptic partial differential operator with arbitrarily rough coefficients in a $d$-dimensional domain. The application of the inverse operator to some standard finite element space defines an approximation space with uniform algebraic approximation rates with respect to the mesh size parameter $H$. This holds even for under-resolved rough coefficients. However, the true challenge of numerical homogenization is the localized computation of a localized basis for such an operator-dependent approximate solution space. This paper presents a novel localization technique that enforces the super-exponential decay of the basis relative to $H$. This shows that basis functions with supports of width $\mathcal O(H|\log H|^{(d-1)/d})$ are sufficient to preserve the optimal algebraic rates of convergence in $H$ without pre-asymptotic effects. A sequence of numerical experiments illustrates the significance of the new localization technique when compared to the so far best localization to supports of width $\mathcal O(H|\log H|)$.
△ Less
Submitted 23 November, 2022; v1 submitted 28 July, 2021;
originally announced July 2021.
-
Operator Compression with Deep Neural Networks
Authors:
Fabian Kröpfl,
Roland Maier,
Daniel Peterseim
Abstract:
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…
▽ More
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.
△ Less
Submitted 28 March, 2022; v1 submitted 25 May, 2021;
originally announced May 2021.
-
Multi-resolution Localized Orthogonal Decomposition for Helmholtz problems
Authors:
Moritz Hauck,
Daniel Peterseim
Abstract:
We introduce a novel multi-resolution Localized Orthogonal Decomposition (LOD) for time-harmonic acoustic scattering problems that can be modeled by the Helmholtz equation. The method merges the concepts of LOD and operator-adapted wavelets (gamblets) and proves its applicability for a class of complex-valued, non-hermitian and indefinite problems. It computes hierarchical bases that block-diagona…
▽ More
We introduce a novel multi-resolution Localized Orthogonal Decomposition (LOD) for time-harmonic acoustic scattering problems that can be modeled by the Helmholtz equation. The method merges the concepts of LOD and operator-adapted wavelets (gamblets) and proves its applicability for a class of complex-valued, non-hermitian and indefinite problems. It computes hierarchical bases that block-diagonalize the Helmholtz operator and thereby decouples the discretization scales. Sparsity is preserved by a novel localization strategy that improves stability properties even in the elliptic case. We present a rigorous stability and a-priori error analysis of the proposed method for homogeneous media. In addition, we investigate the fast solvability of the blocks by a standard iterative method. A sequence of numerical experiments illustrates the sharpness of the theoretical findings and demonstrates the applicability to scattering problems in heterogeneous media.
△ Less
Submitted 23 November, 2022; v1 submitted 22 April, 2021;
originally announced April 2021.
-
Localization and delocalization of ground states of Bose-Einstein condensates under disorder
Authors:
Robert Altmann,
Patrick Henning,
Daniel Peterseim
Abstract:
This paper studies the localization behaviour of Bose-Einstein condensates in disorder potentials, modeled by a Gross-Pitaevskii eigenvalue problem on a bounded interval. In the regime of weak particle interaction, we are able to quantify exponential localization of the ground state, depending on statistical parameters and the strength of the potential. Numerical studies further show delocalizatio…
▽ More
This paper studies the localization behaviour of Bose-Einstein condensates in disorder potentials, modeled by a Gross-Pitaevskii eigenvalue problem on a bounded interval. In the regime of weak particle interaction, we are able to quantify exponential localization of the ground state, depending on statistical parameters and the strength of the potential. Numerical studies further show delocalization if we leave the identified parameter range, which is in agreement with experimental data. These mathematical and numerical findings allow the prediction of physically relevant regimes where localization of ground states may be observed experimentally.
△ Less
Submitted 8 October, 2021; v1 submitted 1 June, 2020;
originally announced June 2020.
-
Reconstruction of quasi-local numerical effective models from low-resolution measurements
Authors:
Alfonso Caiazzo,
Roland Maier,
Daniel Peterseim
Abstract:
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 repres…
▽ More
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.
△ Less
Submitted 4 May, 2020;
originally announced May 2020.
-
A priori error analysis of a numerical stochastic homogenization method
Authors:
Julian Fischer,
Dietmar Gallistl,
Daniel Peterseim
Abstract:
This paper provides an a~priori error analysis of a localized orthogonal decomposition method (LOD) for the numerical stochastic homogenization of a model random diffusion problem. If the uniformly elliptic and bounded random coefficient field of the model problem is stationary and satisfies a quantitative decorrelation assumption in form of the spectral gap inequality, then the expected $L^2$ err…
▽ More
This paper provides an a~priori error analysis of a localized orthogonal decomposition method (LOD) for the numerical stochastic homogenization of a model random diffusion problem. If the uniformly elliptic and bounded random coefficient field of the model problem is stationary and satisfies a quantitative decorrelation assumption in form of the spectral gap inequality, then the expected $L^2$ error of the method can be estimated, up to logarithmic factors, by $H+(\varepsilon/H)^{d/2}$; $\varepsilon$ being the small correlation length of the random coefficient and $H$ the width of the coarse finite element mesh that determines the spatial resolution. The proof bridges recent results of numerical homogenization and quantitative stochastic homogenization.
△ Less
Submitted 2 December, 2020; v1 submitted 25 December, 2019;
originally announced December 2019.
-
The $J$-method for the Gross-Pitaevskii eigenvalue problem
Authors:
Robert Altmann,
Patrick Henning,
Daniel Peterseim
Abstract:
This paper studies the $J$-method of [E. Jarlebring, S. Kvaal, W. Michiels. SIAM J. Sci. Comput. 36-4:A1978-A2001, 2014] for nonlinear eigenvector problems in a general Hilbert space framework. This is the basis for variational discretization techniques and a mesh-independent numerical analysis. A simple modification of the method mimics an energy-decreasing discrete gradient flow. In the case of…
▽ More
This paper studies the $J$-method of [E. Jarlebring, S. Kvaal, W. Michiels. SIAM J. Sci. Comput. 36-4:A1978-A2001, 2014] for nonlinear eigenvector problems in a general Hilbert space framework. This is the basis for variational discretization techniques and a mesh-independent numerical analysis. A simple modification of the method mimics an energy-decreasing discrete gradient flow. In the case of the Gross-Pitaevskii eigenvalue problem, we prove global convergence towards an eigenfunction for a damped version of the $J$-method. More importantly, when the iterations are sufficiently close to an eigenfunction, the damping can be switched off and we recover a local linear convergence rate previously known from the discrete setting. This quantitative convergence analysis is closely connected to the~$J$-method's unique feature of sensitivity with respect to spectral shifts. Contrary to classical gradient flows, this allows both the selective approximation of excited states as well as the amplification of convergence beyond linear rates in the spirit of the Rayleigh quotient iteration for linear eigenvalue problems. These advantageous convergence properties are demonstrated in a series of numerical experiments involving exponentially localized states under disorder potentials and vortex lattices in rotating traps.
△ Less
Submitted 9 December, 2020; v1 submitted 1 August, 2019;
originally announced August 2019.
-
Localized computation of eigenstates of random Schrödinger operators
Authors:
Robert Altmann,
Daniel Peterseim
Abstract:
This paper concerns the numerical approximation of low-energy eigenstates of the linear random Schrödinger operator. Under oscillatory high-amplitude potentials with a sufficient degree of disorder it is known that these eigenstates localize in the sense of an exponential decay of their moduli. We propose a reliable numerical scheme which provides localized approximations of such localized states.…
▽ More
This paper concerns the numerical approximation of low-energy eigenstates of the linear random Schrödinger operator. Under oscillatory high-amplitude potentials with a sufficient degree of disorder it is known that these eigenstates localize in the sense of an exponential decay of their moduli. We propose a reliable numerical scheme which provides localized approximations of such localized states. The method is based on a preconditioned inverse iteration including an optimal multigrid solver which spreads information only locally. The practical performance of the approach is illustrated in various numerical experiments in two and three space dimensions and also for a non-linear random Schrödinger operator.
△ Less
Submitted 7 November, 2019; v1 submitted 22 March, 2019;
originally announced March 2019.
-
Computational high frequency scattering from high contrast heterogeneous media
Authors:
Daniel Peterseim,
Barbara Verfürth
Abstract:
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 presen…
▽ More
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.
△ Less
Submitted 28 January, 2020; v1 submitted 26 February, 2019;
originally announced February 2019.
-
Computational Multiscale Methods for Linear Poroelasticity with High Contrast
Authors:
Shubin Fu,
Robert Altmann,
Eric T. Chung,
Roland Maier,
Daniel Peterseim,
Sai-Mang Pun
Abstract:
In this work, we employ the Constraint Energy Minimizing Generalized Multiscale Finite Element Method (CEM-GMsFEM) to solve the problem of linear heterogeneous poroelasticity with coefficients of high contrast. The proposed method makes use of the idea of energy minimization with suitable constraints in order to generate efficient basis functions for the displacement and the pressure. These basis…
▽ More
In this work, we employ the Constraint Energy Minimizing Generalized Multiscale Finite Element Method (CEM-GMsFEM) to solve the problem of linear heterogeneous poroelasticity with coefficients of high contrast. The proposed method makes use of the idea of energy minimization with suitable constraints in order to generate efficient basis functions for the displacement and the pressure. These basis functions are constructed by solving a class of local auxiliary optimization problems based on eigenfunctions containing local information on the heterogeneity. Techniques of oversampling are adapted to enhance the computational performance. Convergence of first order is shown and illustrated by a number of numerical tests.
△ Less
Submitted 10 December, 2018;
originally announced December 2018.
-
Sobolev gradient flow for the Gross-Pitaevskii eigenvalue problem: global convergence and computational efficiency
Authors:
Patrick Henning,
Daniel Peterseim
Abstract:
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.…
▽ More
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.
△ Less
Submitted 2 April, 2020; v1 submitted 3 December, 2018;
originally announced December 2018.
-
From Domain Decomposition to Homogenization Theory
Authors:
Daniel Peterseim,
Dora Varga,
Barbara Verfürth
Abstract:
This paper rediscovers a classical homogenization result for a prototypical linear elliptic boundary value problem with periodically oscillating diffusion coefficient. Unlike classical analytical approaches such as asymptotic analysis, oscillating test functions, or two-scale convergence, the result is purely based on the theory of domain decomposition methods and standard finite elements techniqu…
▽ More
This paper rediscovers a classical homogenization result for a prototypical linear elliptic boundary value problem with periodically oscillating diffusion coefficient. Unlike classical analytical approaches such as asymptotic analysis, oscillating test functions, or two-scale convergence, the result is purely based on the theory of domain decomposition methods and standard finite elements techniques. The arguments naturally generalize to problems far beyond periodicity and scale separation and we provide a brief overview on such applications.
△ Less
Submitted 15 November, 2018;
originally announced November 2018.
-
Sparse Compression of Expected Solution Operators
Authors:
Michael Feischl,
Daniel Peterseim
Abstract:
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…
▽ More
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.
△ Less
Submitted 16 March, 2020; v1 submitted 4 July, 2018;
originally announced July 2018.
-
Quantitative Anderson localization of Schrödinger eigenstates under disorder potentials
Authors:
Robert Altmann,
Patrick Henning,
Daniel Peterseim
Abstract:
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…
▽ More
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.
△ Less
Submitted 9 February, 2020; v1 submitted 27 March, 2018;
originally announced March 2018.
-
Explicit Computational Wave Propagation in Micro-Heterogeneous Media
Authors:
Roland Maier,
Daniel Peterseim
Abstract:
Explicit time stepping schemes are popular for linear acoustic and elastic wave propagation due to their simple nature which does not require sophisticated solvers for the inversion of the stiffness matrices. However, explicit schemes are only stable if the time step size is bounded by the mesh size in space subject to the so-called CFL condition. In micro-heterogeneous media, this condition is ty…
▽ More
Explicit time stepping schemes are popular for linear acoustic and elastic wave propagation due to their simple nature which does not require sophisticated solvers for the inversion of the stiffness matrices. However, explicit schemes are only stable if the time step size is bounded by the mesh size in space subject to the so-called CFL condition. In micro-heterogeneous media, this condition is typically prohibitively restrictive because spatial oscillations of the medium need to be resolved by the discretization in space. This paper presents a way to reduce the spatial complexity in such a setting and, hence, to enable a relaxation of the CFL condition. This is done using the Localized Orthogonal Decomposition method as a tool for numerical homogenization. A complete convergence analysis is presented with appropriate, weak regularity assumptions on the initial data.
△ Less
Submitted 15 October, 2018; v1 submitted 21 March, 2018;
originally announced March 2018.
-
Computational multiscale methods for linear heterogeneous poroelasticity
Authors:
Robert Altmann,
Eric Chung,
Roland Maier,
Daniel Peterseim,
Sai-Mang Pun
Abstract:
We consider a strongly heterogeneous medium saturated by an incompressible viscous fluid as it appears in geomechanical modeling. This poroelasticity problem suffers from rapidly oscillating material parameters, which calls for a thorough numerical treatment. In this paper, we propose a method based on the local orthogonal decomposition technique and motivated by a similar approach used for linear…
▽ More
We consider a strongly heterogeneous medium saturated by an incompressible viscous fluid as it appears in geomechanical modeling. This poroelasticity problem suffers from rapidly oscillating material parameters, which calls for a thorough numerical treatment. In this paper, we propose a method based on the local orthogonal decomposition technique and motivated by a similar approach used for linear thermoelasticity. Therein, local corrector problems are constructed in line with the static equations, whereas we propose to consider the full system. This allows to benefit from the given saddle point structure and results in two decoupled corrector problems for the displacement and the pressure. We prove the optimal first-order convergence of this method and verify the result by numerical experiments.
△ Less
Submitted 22 December, 2018; v1 submitted 2 January, 2018;
originally announced January 2018.
-
Thermo-optical interactions in a dye-microcavity photon Bose-Einstein condensate
Authors:
Hadiseh Alaeian,
Mira Schedensack,
Clara Bartels,
Daniel Peterseim,
Martin Weitz
Abstract:
Superfluidity and Bose-Einstein condensation are usually considered as two closely related phenomena. Indeed, in most macroscopic quantum systems, like liquid helium, ultracold atomic Bose gases, and exciton-polaritons, condensation and superfluidity occur in parallel. In photon Bose-Einstein condensates realized in the dye microcavity system, thermalization does not occur by direct interaction of…
▽ More
Superfluidity and Bose-Einstein condensation are usually considered as two closely related phenomena. Indeed, in most macroscopic quantum systems, like liquid helium, ultracold atomic Bose gases, and exciton-polaritons, condensation and superfluidity occur in parallel. In photon Bose-Einstein condensates realized in the dye microcavity system, thermalization does not occur by direct interaction of the condensate particles as in the above described systems, i.e. photon-photon interactions, but by absorption and re-emission processes on the dye molecules, which act as a heat reservoir. Currently, there is no experimental evidence for superfluidity in the dye microcavity system, though effective photon interactions have been observed from thermo-optic effects in the dye medium. In this work, we theoretically investigate the implications of effective thermo-optic photon interactions, a temporally delayed and spatially non-local effect, on the photon condensate, and derive the resulting Bogoliubov excitation spectrum. The calculations suggest a linear photon dispersion at low momenta, fulfilling the Landau's criterion of superfluidity . We envision that the temporally delayed and long-range nature of the thermo-optic photon interaction offer perspectives for novel quantum fluid phenomena.
△ Less
Submitted 4 October, 2017;
originally announced October 2017.
-
Numerical Homogenization of Heterogeneous Fractional Laplacians
Authors:
Donald L. Brown,
Joscha Gedicke,
Daniel Peterseim
Abstract:
In this paper, we develop a numerical multiscale method to solve the fractional Laplacian with a heterogeneous diffusion coefficient. When the coefficient is heterogeneous, this adds to the computational costs. Moreover, the fractional Laplacian is a nonlocal operator in its standard form, however the Caffarelli-Silvestre extension allows for a localization of the equations. This adds a complexity…
▽ More
In this paper, we develop a numerical multiscale method to solve the fractional Laplacian with a heterogeneous diffusion coefficient. When the coefficient is heterogeneous, this adds to the computational costs. Moreover, the fractional Laplacian is a nonlocal operator in its standard form, however the Caffarelli-Silvestre extension allows for a localization of the equations. This adds a complexity of an extra spacial dimension and a singular/degenerate coefficient depending on the fractional order. Using a sub-grid correction method, we correct the basis functions in a natural weighted Sobolev space and show that these corrections are able to be truncated to design a computationally efficient scheme with optimal convergence rates. A key ingredient of this method is the use of quasi-interpolation operators to construct the fine scale spaces. Since the solution of the extended problem on the critical boundary is of main interest, we construct a projective quasi-interpolation that has both $d$ and $d+1$ dimensional averages over subsets in the spirit of the Scott-Zhang operator. We show that this operator satisfies local stability and local approximation properties in weighted Sobolev spaces. We further show that we can obtain a greater rate of convergence for sufficient smooth forces, and utilizing a global $L^2$ projection on the critical boundary. We present some numerical examples, utilizing our projective quasi-interpolation in dimension $2+1$ for analytic and heterogeneous cases to demonstrate the rates and effectiveness of the method.
△ Less
Submitted 3 September, 2017;
originally announced September 2017.
-
Numerical stochastic homogenization by quasilocal effective diffusion tensors
Authors:
Dietmar Gallistl,
Daniel Peterseim
Abstract:
This paper proposes a numerical upscaling procedure for elliptic boundary value problems with diffusion tensors that vary randomly on small scales. The resulting effective deterministic model is given through a quasilocal discrete integral operator, which can be further compressed to an effective partial differential operator. Error estimates consisting of a priori and a posteriori terms are provi…
▽ More
This paper proposes a numerical upscaling procedure for elliptic boundary value problems with diffusion tensors that vary randomly on small scales. The resulting effective deterministic model is given through a quasilocal discrete integral operator, which can be further compressed to an effective partial differential operator. Error estimates consisting of a priori and a posteriori terms are provided that allow one to quantify the impact of uncertainty in the diffusion coefficient on the expected effective response of the process.
△ Less
Submitted 22 January, 2019; v1 submitted 28 February, 2017;
originally announced February 2017.
-
An analysis of a class of variational multiscale methods based on subspace decomposition
Authors:
Ralf Kornhuber,
Daniel Peterseim,
Harry Yserentant
Abstract:
Numerical homogenization tries to approximate the solutions of elliptic partial differential equations with strongly oscillating coefficients by functions from modified finite element spaces. We present in this paper a class of such methods that are very closely related to the method of Målqvist and Peterseim [Math. Comp. 83, 2014]. Like the method of Målqvist and Peterseim, these methods do not m…
▽ More
Numerical homogenization tries to approximate the solutions of elliptic partial differential equations with strongly oscillating coefficients by functions from modified finite element spaces. We present in this paper a class of such methods that are very closely related to the method of Målqvist and Peterseim [Math. Comp. 83, 2014]. Like the method of Målqvist and Peterseim, these methods do not make explicit or implicit use of a scale separation. Their compared to that in the work of Målqvist and Peterseim strongly simplified analysis is based on a reformulation of their method in terms of variational multiscale methods and on the theory of iterative methods, more precisely, of additive Schwarz or subspace decomposition methods.
△ Less
Submitted 21 January, 2018; v1 submitted 14 August, 2016;
originally announced August 2016.
-
Crank-Nicolson Galerkin approximations to nonlinear Schrödinger equations with rough potentials
Authors:
Patrick Henning,
Daniel Peterseim
Abstract:
This paper analyses the numerical solution of a class of non-linear Schrödinger equations by Galerkin finite elements in space and a mass- and energy conserving variant of the Crank-Nicolson method due to Sanz-Serna in time. The novel aspects of the analysis are the incorporation of rough, discontinuous potentials in the context of weak and strong disorder, the consideration of some general class…
▽ More
This paper analyses the numerical solution of a class of non-linear Schrödinger equations by Galerkin finite elements in space and a mass- and energy conserving variant of the Crank-Nicolson method due to Sanz-Serna in time. The novel aspects of the analysis are the incorporation of rough, discontinuous potentials in the context of weak and strong disorder, the consideration of some general class of non-linearities, and the proof of convergence with rates in $L^{\infty}(L^2)$ under moderate regularity assumptions that are compatible with discontinuous potentials. For sufficiently smooth potentials, the rates are optimal without any coupling condition between the time step size and the spatial mesh width.
△ Less
Submitted 17 June, 2017; v1 submitted 7 August, 2016;
originally announced August 2016.
-
Computation of quasilocal effective diffusion tensors and connections to the mathematical theory of homogenization
Authors:
Dietmar Gallistl,
Daniel Peterseim
Abstract:
This paper aims at bridging existing theories in numerical and analytical homogenization. For this purpose the multiscale method of Målqvist and Peterseim [Math. Comp. 2014], which is based on orthogonal subspace decomposition, is reinterpreted by means of a discrete integral operator acting on standard finite element spaces. The exponential decay of the involved integral kernel motivates the use…
▽ More
This paper aims at bridging existing theories in numerical and analytical homogenization. For this purpose the multiscale method of Målqvist and Peterseim [Math. Comp. 2014], which is based on orthogonal subspace decomposition, is reinterpreted by means of a discrete integral operator acting on standard finite element spaces. The exponential decay of the involved integral kernel motivates the use of a diagonal approximation and, hence, a localized piecewise constant coefficient. In a periodic setting, the computed localized coefficient is proved to coincide with the classical homogenization limit. An a priori error analysis shows that the local numerical model is appropriate beyond the periodic setting when the localized coefficient satisfies a certain homogenization criterion, which can be verified a posteriori. The results are illustrated in numerical experiments.
△ Less
Submitted 23 October, 2017; v1 submitted 6 August, 2016;
originally announced August 2016.
-
Error analysis of a variational multiscale stabilization for convection-dominated diffusion equations in 2d
Authors:
Guanglian Li,
Daniel Peterseim,
Mira Schedensack
Abstract:
We formulate a stabilized quasi-optimal Petrov-Galerkin method for singularly perturbed convection-diffusion problems based on the variational multiscale method. The stabilization is of Petrov-Galerkin type with a standard finite element trial space and a problem-dependent test space based on pre-computed fine-scale correctors. The exponential decay of these correctors and their localisation to lo…
▽ More
We formulate a stabilized quasi-optimal Petrov-Galerkin method for singularly perturbed convection-diffusion problems based on the variational multiscale method. The stabilization is of Petrov-Galerkin type with a standard finite element trial space and a problem-dependent test space based on pre-computed fine-scale correctors. The exponential decay of these correctors and their localisation to local patch problems, which depend on the direction of the velocity field and the singular perturbation parameter, is rigorously justified. Under moderate assumptions, this stabilization guarantees stability and quasi-optimal rate of convergence for arbitrary mesh Péclet numbers on fairly coarse meshes at the cost of additional inter-element communication.
△ Less
Submitted 15 June, 2016;
originally announced June 2016.
-
Adaptive Mesh Refinement Strategies in Isogeometric Analysis - A Computational Comparison
Authors:
Paul Hennig,
Markus Kästner,
Philipp Morgenstern,
Daniel Peterseim
Abstract:
We explain four variants of an adaptive finite element method with cubic splines and compare their performance in simple elliptic model problems. The methods in comparison are Truncated Hierarchical B-splines with two different refinement strategies, T-splines with the refinement strategy introduced by Scott et al. in 2012, and T-splines with an alternative refinement strategy introduced by some o…
▽ More
We explain four variants of an adaptive finite element method with cubic splines and compare their performance in simple elliptic model problems. The methods in comparison are Truncated Hierarchical B-splines with two different refinement strategies, T-splines with the refinement strategy introduced by Scott et al. in 2012, and T-splines with an alternative refinement strategy introduced by some of the authors. In four examples, including singular and non-singular problems of linear elasticity and the Poisson problem, the H1-errors of the discrete solutions, the number of degrees of freedom as well as sparsity patterns and condition numbers of the discretized problem are compared.
△ Less
Submitted 3 May, 2016;
originally announced May 2016.
-
Efficient implementation of the Localized Orthogonal Decomposition method
Authors:
Christian Engwer,
Patrick Henning,
Axel Målqvist,
Daniel Peterseim
Abstract:
In this paper we present algorithms for an efficient implementation of the Localized Orthogonal Decomposition method (LOD). The LOD is a multiscale method for the numerical simulation of partial differential equations with a continuum of inseparable scales. We show how the method can be implemented in a fairly standard Finite Element framework and discuss its realization for different types of pro…
▽ More
In this paper we present algorithms for an efficient implementation of the Localized Orthogonal Decomposition method (LOD). The LOD is a multiscale method for the numerical simulation of partial differential equations with a continuum of inseparable scales. We show how the method can be implemented in a fairly standard Finite Element framework and discuss its realization for different types of problems, such as linear elliptic problems with rough coefficients and linear eigenvalue problems.
△ Less
Submitted 20 February, 2019; v1 submitted 4 February, 2016;
originally announced February 2016.
-
Robust Numerical Upscaling of Elliptic Multiscale Problems at High Contrast
Authors:
Daniel Peterseim,
Robert Scheichl
Abstract:
We present a new approach to the numerical upscaling for elliptic problems with rough diffusion coefficient at high contrast. It is based on the localizable orthogonal decomposition of $H^1$ into the image and the kernel of some novel stable quasi-interpolation operators with local $L^2$-approximation properties, independent of the contrast. We identify a set of sufficient assumptions on these qua…
▽ More
We present a new approach to the numerical upscaling for elliptic problems with rough diffusion coefficient at high contrast. It is based on the localizable orthogonal decomposition of $H^1$ into the image and the kernel of some novel stable quasi-interpolation operators with local $L^2$-approximation properties, independent of the contrast. We identify a set of sufficient assumptions on these quasi-interpolation operators that guarantee in principle optimal convergence without pre-asymptotic effects for high-contrast coefficients. We then give an example of a suitable operator and establish the assumptions for a particular class of high-contrast coefficients. So far this is not possible without any pre-asymptotic effects, but the optimal convergence is independent of the contrast and the asymptotic range is largely improved over other discretisation schemes. The new framework is sufficiently flexible to allow also for other choices of quasi-interpolation operators and the potential for fully robust numerical upscaling at high contrast.
△ Less
Submitted 25 January, 2016;
originally announced January 2016.
-
Relaxing the CFL condition for the wave equation on adaptive meshes
Authors:
Daniel Peterseim,
Mira Schedensack
Abstract:
The Courant-Friedrichs-Lewy (CFL) condition guarantees the stability of the popular explicit leapfrog method for the wave equation. However, it limits the choice of the time step size to be bounded by the minimal mesh size in the spatial finite element mesh. This essentially prohibits any sort of adaptive mesh refinement that would be required to reveal optimal convergence rates on domains with re…
▽ More
The Courant-Friedrichs-Lewy (CFL) condition guarantees the stability of the popular explicit leapfrog method for the wave equation. However, it limits the choice of the time step size to be bounded by the minimal mesh size in the spatial finite element mesh. This essentially prohibits any sort of adaptive mesh refinement that would be required to reveal optimal convergence rates on domains with re-entrant corners. This paper shows how a simple subspace projection step inspired by numerical homogenisation can remove the critical time step restriction so that the CFL condition and approximation properties are balanced in an optimal way, even in the presence of spatial singularities.
△ Less
Submitted 24 February, 2017; v1 submitted 19 January, 2016;
originally announced January 2016.
-
Multiscale Petrov-Galerkin Method for High-Frequency Heterogeneous Helmholtz Equations
Authors:
Donald L. Brown,
Dietmar Gallistl,
Daniel Peterseim
Abstract:
This paper presents a multiscale Petrov-Galerkin finite element method for time-harmonic acoustic scattering problems with heterogeneous coefficients in the high-frequency regime. We show that the method is pollution- free also in the case of heterogeneous media provided that the stability bound of the continuous problem grows at most polynomially with the wave number k. By generalizing classical…
▽ More
This paper presents a multiscale Petrov-Galerkin finite element method for time-harmonic acoustic scattering problems with heterogeneous coefficients in the high-frequency regime. We show that the method is pollution- free also in the case of heterogeneous media provided that the stability bound of the continuous problem grows at most polynomially with the wave number k. By generalizing classical estimates of [Melenk, Ph.D. Thesis 1995] and [Hetmaniuk, Commun. Math. Sci. 5 (2007)] for homogeneous medium, we show that this assumption of polynomially wave number growth holds true for a particular class of smooth heterogeneous material coefficients. Further, we present numerical examples to verify our stability estimates and implement an example in the wider class of discontinuous coefficients to show computational applicability beyond our limited class of coefficients.
△ Less
Submitted 30 November, 2015;
originally announced November 2015.
-
Generalized finite element methods for quadratic eigenvalue problems
Authors:
Axel Målqvist,
Daniel Peterseim
Abstract:
We consider a large-scale quadratic eigenvalue problem (QEP), formulated using P1 finite elements on a fine scale reference mesh. This model describes damped vibrations in a structural mechanical system. In particular we focus on problems with rapid material data variation, e.g., composite materials. We construct a low dimensional generalized finite element (GFE) space based on the localized ortho…
▽ More
We consider a large-scale quadratic eigenvalue problem (QEP), formulated using P1 finite elements on a fine scale reference mesh. This model describes damped vibrations in a structural mechanical system. In particular we focus on problems with rapid material data variation, e.g., composite materials. We construct a low dimensional generalized finite element (GFE) space based on the localized orthogonal decomposition (LOD) technique. The construction involves the (parallel) solution of independent localized linear Poisson-type problems. The GFE space is then used to compress the large-scale algebraic QEP to a much smaller one with a similar modeling accuracy. The small scale QEP can then be solved by standard techniques at a significantly reduced computational cost. We prove convergence with rate for the proposed method and numerical experiments confirm our theoretical findings.
△ Less
Submitted 20 October, 2015;
originally announced October 2015.
-
Complexity of hierarchical refinement for a class of admissible mesh configurations
Authors:
Annalisa Buffa,
Carlotta Giannelli,
Philipp Morgenstern,
Daniel Peterseim
Abstract:
An adaptive isogeometric method based on $d$-variate hierarchical spline constructions can be derived by considering a refine module that preserves a certain class of admissibility between two consecutive steps of the adaptive loop [6]. In this paper we provide a complexity estimate, i.e., an estimate on how the number of mesh elements grows with respect to the number of elements that are marked f…
▽ More
An adaptive isogeometric method based on $d$-variate hierarchical spline constructions can be derived by considering a refine module that preserves a certain class of admissibility between two consecutive steps of the adaptive loop [6]. In this paper we provide a complexity estimate, i.e., an estimate on how the number of mesh elements grows with respect to the number of elements that are marked for refinement by the adaptive strategy. Our estimate is in the line of the similar ones proved in the finite element context, [3,24].
△ Less
Submitted 18 September, 2015;
originally announced September 2015.
-
Variational Multiscale Stabilization and the Exponential Decay of Fine-scale Correctors
Authors:
Daniel Peterseim
Abstract:
This paper addresses the variational multiscale stabilization of standard finite element methods for linear partial differential equations that exhibit multiscale features. The stabilization is of Petrov-Galerkin type with a standard finite element trial space and a problem-dependent test space based on pre-computed fine-scale correctors. The exponential decay of these correctors and their localis…
▽ More
This paper addresses the variational multiscale stabilization of standard finite element methods for linear partial differential equations that exhibit multiscale features. The stabilization is of Petrov-Galerkin type with a standard finite element trial space and a problem-dependent test space based on pre-computed fine-scale correctors. The exponential decay of these correctors and their localisation to local cell problems is rigorously justified. The stabilization eliminates scale-dependent pre-asymptotic effects as they appear for standard finite element discretizations of highly oscillatory problems, e.g., the poor $L^2$ approximation in homogenization problems or the pollution effect in high-frequency acoustic scattering.
△ Less
Submitted 20 October, 2015; v1 submitted 28 May, 2015;
originally announced May 2015.
-
Stable Multiscale Petrov-Galerkin Finite Element Method for High Frequency Acoustic Scattering
Authors:
Dietmar Gallistl,
Daniel Peterseim
Abstract:
We present and analyze a pollution-free Petrov-Galerkin multiscale finite element method for the Helmholtz problem with large wave number $κ$ as a variant of [Peterseim, ArXiv:1411.1944, 2014]. We use standard continuous $Q_1$ finite elements at a coarse discretization scale $H$ as trial functions, whereas the test functions are computed as the solutions of local problems at a finer scale $h$. The…
▽ More
We present and analyze a pollution-free Petrov-Galerkin multiscale finite element method for the Helmholtz problem with large wave number $κ$ as a variant of [Peterseim, ArXiv:1411.1944, 2014]. We use standard continuous $Q_1$ finite elements at a coarse discretization scale $H$ as trial functions, whereas the test functions are computed as the solutions of local problems at a finer scale $h$. The diameter of the support of the test functions behaves like $mH$ for some oversampling parameter $m$. Provided $m$ is of the order of $\log(κ)$ and $h$ is sufficiently small, the resulting method is stable and quasi-optimal in the regime where $H$ is proportional to $κ^{-1}$. In homogeneous (or more general periodic) media, the fine scale test functions depend only on local mesh-configurations. Therefore, the seemingly high cost for the computation of the test functions can be drastically reduced on structured meshes. We present numerical experiments in two and three space dimensions.
△ Less
Submitted 14 July, 2015; v1 submitted 17 March, 2015;
originally announced March 2015.