-
Construction of local reduced spaces for Friedrichs' systems via randomized training
Authors:
Christian Engwer,
Mario Ohlberger,
Lukas Renelt
Abstract:
This contribution extends the localized training approach, traditionally employed for multiscale problems and parameterized partial differential equations (PDEs) featuring locally heterogeneous coefficients, to the class of linear, positive symmetric operators, known as Friedrichs' operators. Considering a local subdomain with corresponding oversampling domain we prove the compactness of the trans…
▽ More
This contribution extends the localized training approach, traditionally employed for multiscale problems and parameterized partial differential equations (PDEs) featuring locally heterogeneous coefficients, to the class of linear, positive symmetric operators, known as Friedrichs' operators. Considering a local subdomain with corresponding oversampling domain we prove the compactness of the transfer operator which maps boundary data to solutions on the interior domain. While a Caccioppoli-inequality quantifying the energy decay to the interior holds true for all Friedrichs' systems, showing a compactness result for the graph-spaces hosting the solution is additionally necessary. We discuss the mixed formulation of a convection-diffusion-reaction problem where the necessary compactness result is obtained by the Picard-Weck-Weber theorem. Our numerical results, focusing on a scenario involving heterogeneous diffusion fields with multiple high-conductivity channels, demonstrate the effectiveness of the proposed method.
△ Less
Submitted 29 April, 2024;
originally announced April 2024.
-
Model order reduction of an ultraweak and optimally stable variational formulation for parametrized reactive transport problems
Authors:
Christian Engwer,
Mario Ohlberger,
Lukas Renelt
Abstract:
This contribution introduces a model order reduction approach for an advection-reaction problem with a parametrized reaction function. The underlying discretization uses an ultraweak formulation with an $L^2$-like trial space and an 'optimal' test space as introduced by Demkowicz et al. This ensures the stability of the discretization and in addition allows for a symmetric reformulation of the pro…
▽ More
This contribution introduces a model order reduction approach for an advection-reaction problem with a parametrized reaction function. The underlying discretization uses an ultraweak formulation with an $L^2$-like trial space and an 'optimal' test space as introduced by Demkowicz et al. This ensures the stability of the discretization and in addition allows for a symmetric reformulation of the problem in terms of a dual solution which can also be interpreted as the normal equations of an adjoint least-squares problem. Classic model order reduction techniques can then be applied to the space of dual solutions which also immediately gives a reduced primal space. We show that the necessary computations do not require the reconstruction of any primal solutions and can instead be performed entirely on the space of dual solutions. We prove exponential convergence of the Kolmogorov $N$-width and show that a greedy algorithm produces quasi-optimal approximation spaces for both the primal and the dual solution space. Numerical experiments based on the benchmark problem of a catalytic filter confirm the applicability of the proposed method.
△ Less
Submitted 28 June, 2024; v1 submitted 30 October, 2023;
originally announced October 2023.
-
Domain of Dependence stabilization for the acoustic wave equation on 2D cut-cell meshes
Authors:
Gunnar Birke,
Christian Engwer,
Sandra May,
Florian Streitbürger
Abstract:
Cut-cell meshes are an attractive alternative to avoid common mesh generation problems. For hyperbolic problems they pose additional challenges, as elements can become arbitrarily small, leading to prohibitive time step restrictions for explicit time stepping methods. To alleviate this small cell problem we consider a particular stabilization method, the Domain of Dependence (DoD) method. So far,…
▽ More
Cut-cell meshes are an attractive alternative to avoid common mesh generation problems. For hyperbolic problems they pose additional challenges, as elements can become arbitrarily small, leading to prohibitive time step restrictions for explicit time stepping methods. To alleviate this small cell problem we consider a particular stabilization method, the Domain of Dependence (DoD) method. So far, while posessing many favorable theoretical properties, in two dimensions the DoD method was essentially restricted to the transport equation. In this work we extend the DoD method to the acoustic wave equation in two dimensions and provide numerical results for validation.
△ Less
Submitted 9 April, 2023;
originally announced April 2023.
-
An optimally stable approximation of reactive transport using discrete test and infinite trial spaces
Authors:
Lukas Renelt,
Christian Engwer,
Mario Ohlberger
Abstract:
In this contribution we propose an optimally stable ultraweak Petrov-Galerkin variational formulation and subsequent discretization for stationary reactive transport problems. The discretization is exclusively based on the choice of discrete approximate test spaces, while the trial space is a priori infinite dimensional. The solution in the trial space or even only functional evaluations of the so…
▽ More
In this contribution we propose an optimally stable ultraweak Petrov-Galerkin variational formulation and subsequent discretization for stationary reactive transport problems. The discretization is exclusively based on the choice of discrete approximate test spaces, while the trial space is a priori infinite dimensional. The solution in the trial space or even only functional evaluations of the solution are obtained in a post-processing step. We detail the theoretical framework and demonstrate its usage in a numerical experiment that is motivated from modeling of catalytic filters.
△ Less
Submitted 29 May, 2023; v1 submitted 28 March, 2023;
originally announced March 2023.
-
The Localized Subtraction Approach For EEG and MEG Forward Modeling
Authors:
Malte B. Höltershinken,
Pia Lange,
Tim Erdbrügger,
Yvonne Buschermöhle,
Fabrice Wallois,
Alena Buyx,
Sampsa Pursiainen,
Johannes Vorwerk,
Christian Engwer,
Carsten H. Wolters
Abstract:
In FEM-based EEG and MEG source analysis, the subtraction approach has been proposed to simulate sensor measurements generated by neural activity. While this approach possesses a rigorous foundation and produces accurate results, its major downside is that it is computationally prohibitively expensive in practical applications. To overcome this, we developed a new approach, called the localized su…
▽ More
In FEM-based EEG and MEG source analysis, the subtraction approach has been proposed to simulate sensor measurements generated by neural activity. While this approach possesses a rigorous foundation and produces accurate results, its major downside is that it is computationally prohibitively expensive in practical applications. To overcome this, we developed a new approach, called the localized subtraction approach. This approach is designed to preserve the mathematical foundation of the subtraction approach, while also leading to sparse right-hand sides in the FEM formulation, making it efficiently computable. We achieve this by introducing a cut-off into the subtraction, restricting its influence to the immediate neighborhood of the source. In this work, this approach will be presented, analyzed, and compared to other state-of-the-art FEM right-hand side approaches. Furthermore, we discuss how to arrive at an efficient and stable implementation. We perform validation in multi-layer sphere models where analytical solutions exist. There, we demonstrate that the localized subtraction approach is vastly more efficient than the subtraction approach. Moreover, we find that for the EEG forward problem, the localized subtraction approach is less dependent on the global structure of the FEM mesh when compared to the subtraction approach. Additionally, we show the localized subtraction approach to rival, and in many cases even surpass, the other investigated approaches in terms of accuracy. For the MEG forward problem, we show the localized subtraction approach and the subtraction approach to produce highly accurate approximations of the volume currents close to the source. The localized subtraction approach thus reduces the computational cost of the subtraction approach to an extent that makes it usable in practical applications without sacrificing rigorousness and accuracy.
△ Less
Submitted 24 February, 2023;
originally announced February 2023.
-
DoD stabilization for higher-order advection in two dimensions
Authors:
Florian Streitbürger,
Gunnar Birke,
Christian Engwer,
Sandra May
Abstract:
When solving time-dependent hyperbolic conservation laws on cut cell meshes one has to overcome the small cell problem: standard explicit time stepping is not stable on small cut cells if the time step is chosen with respect to larger background cells. The domain of dependence (DoD) stabilization is designed to solve this problem in a discontinuous Galerkin framework. It adds a penalty term to the…
▽ More
When solving time-dependent hyperbolic conservation laws on cut cell meshes one has to overcome the small cell problem: standard explicit time stepping is not stable on small cut cells if the time step is chosen with respect to larger background cells. The domain of dependence (DoD) stabilization is designed to solve this problem in a discontinuous Galerkin framework. It adds a penalty term to the space discretization that restores proper domains of dependency. In this contribution we introduce the DoD stabilization for solving the advection equation in 2d with higher order. We show an $L^2$ stability result for the stabilized semi-discrete scheme for arbitrary polynomial degrees $p$ and provide numerical results for convergence tests indicating orders of $p+1$ in the $L^1$ norm and between $p+\frac 1 2$ and $p+1$ in the $L^{\infty}$ norm.
△ Less
Submitted 7 January, 2023;
originally announced January 2023.
-
DoD Stabilization of linear hyperbolic PDEs on general cut-cell meshes
Authors:
Gunnar Birke,
Christian Engwer,
Sandra May,
Florian Streitbürger
Abstract:
Standard numerical methods for hyperbolic PDEs require for stability a CFL-condition which implies that the time step size depends on the size of the elements of the mesh. On cut-cell meshes, elements can become arbitrarily small and thus the time step size cannot take the size of small cut-cells into account but has to be chosen based on the background mesh elements. A remedy for this is the so c…
▽ More
Standard numerical methods for hyperbolic PDEs require for stability a CFL-condition which implies that the time step size depends on the size of the elements of the mesh. On cut-cell meshes, elements can become arbitrarily small and thus the time step size cannot take the size of small cut-cells into account but has to be chosen based on the background mesh elements. A remedy for this is the so called DoD (domain of dependence) stabilization for which several favorable theoretical and numerical properties have been shown in one and two space dimensions. Up to now the method is restricted to stabilization of cut-cells with exactly one inflow and one outflow face, i.e. triangular cut-cells with a no-flow face. We extend the DoD stabilization to cut-cells with multiple in- and out-flow faces by properly considering the flow distribution inside the cut-cell. We further prove L2-stability for the semi-discrete formulation in space and present numerical results to validate the proposed extension.
△ Less
Submitted 6 January, 2023;
originally announced January 2023.
-
CutFEM forward modeling for EEG source analysis
Authors:
Tim Erdbrügger,
Andreas Westhoff,
Malte Hoeltershinken,
Jan-Ole Radecke,
Yvonne Buschermoehle,
Alena Buyx,
Fabrice Wallois,
Sampsa Pursiainen,
Joachim Gross,
Rebekka Lencer,
Christian Engwer,
Carsten Wolters
Abstract:
Source analysis of Electroencephalography (EEG) data requires the computation of the scalp potential induced by current sources in the brain. This so-called EEG forward problem is based on an accurate estimation of the volume conduction effects in the human head, represented by a partial differential equation which can be solved using the finite element method (FEM). FEM offers flexibility when mo…
▽ More
Source analysis of Electroencephalography (EEG) data requires the computation of the scalp potential induced by current sources in the brain. This so-called EEG forward problem is based on an accurate estimation of the volume conduction effects in the human head, represented by a partial differential equation which can be solved using the finite element method (FEM). FEM offers flexibility when modeling anisotropic tissue conductivities but requires a volumetric discretization, a mesh, of the head domain. Structured hexahedral meshes are easy to create in an automatic fashion, while tetrahedral meshes are better suited to model curved geometries. Tetrahedral meshes thus offer better accuracy, but are more difficult to create. Methods: We introduce CutFEM for EEG forward simulations to integrate the strengths of hexahedra and tetrahedra. It belongs to the family of unfitted finite element methods, decoupling mesh and geometry representation. Following a description of the method, we will employ CutFEM in both controlled spherical scenarios and the reconstruction of somatosensory evoked potentials. Results: CutFEM outperforms competing FEM approaches with regard to numerical accuracy, memory consumption and computational speed while being able to mesh arbitrarily touching compartments. Conclusion: CutFEM balances numerical accuracy, computational efficiency and a smooth approximation of complex geometries that has previously not been available in FEM-based EEG forward modeling.
△ Less
Submitted 30 November, 2022;
originally announced November 2022.
-
A Hardware-aware and Stable Orthogonalization Framework
Authors:
Nils-Arne Dreier,
Christian Engwer
Abstract:
The orthogonalization process is an essential building block in Krylov space methods, which takes up a large portion of the computational time. Commonly used methods, like the Gram-Schmidt method, consider the projection and normalization separately and store the orthogonal base explicitly. We consider the problem of orthogonalization and normalization as a QR decomposition problem on which we app…
▽ More
The orthogonalization process is an essential building block in Krylov space methods, which takes up a large portion of the computational time. Commonly used methods, like the Gram-Schmidt method, consider the projection and normalization separately and store the orthogonal base explicitly. We consider the problem of orthogonalization and normalization as a QR decomposition problem on which we apply known algorithms, namely CholeskyQR and TSQR. This leads to methods that solve the orthogonlization problem with reduced communication costs, while maintaining stability and stores the orthogonal base in a locally orthogonal representation. Furthermore, we discuss the novel method as a framework which allows us to combine different orthogonalization algorithms and use the best algorithm for each part of the hardware. After the formulation of the methods, we show their advantageous performance properties based on a performance model that takes data transfers within compute nodes as well as message passing between compute nodes into account. The theoretic results are validated by numerical experiments.
△ Less
Submitted 28 April, 2022;
originally announced April 2022.
-
Realistic head modeling of electromagnetic brain activity: An integrated Brainstorm pipeline from MRI data to the FEM solution
Authors:
Takfarinas Medani,
Juan Garcia-Prieto,
Francois Tadel,
Sophie Schrader,
Anand Joshi,
Christian Engwer,
Carsten H. Wolters,
John C. Mosher,
Richard M. Leahy
Abstract:
Human brain activity generates scalp potentials (electroencephalography EEG), intracranial potentials (iEEG), and external magnetic fields (magnetoencephalography MEG), all capable of being recorded, often simultaneously, for use in research and clinical applications. The so-called forward problem is the modeling of these fields at their sensors for a given putative neural source configuration. Wh…
▽ More
Human brain activity generates scalp potentials (electroencephalography EEG), intracranial potentials (iEEG), and external magnetic fields (magnetoencephalography MEG), all capable of being recorded, often simultaneously, for use in research and clinical applications. The so-called forward problem is the modeling of these fields at their sensors for a given putative neural source configuration. While early generations modeled the head as a simple set of isotropic spheres, today s ubiquitous magnetic resonance imaging (MRI) data allows detailed descriptions of head compartments with assigned isotropic and anisotropic conductivities. In this paper, we present a complete pipeline, integrated into the Brainstorm software, that allows users to generate an individual and accurate head model from the MRI and then calculate the electromagnetic forward solution using the finite element method (FEM). The head model generation is performed by the integration of the latest tools for MRI segmentation and FEM mesh generation. The final head model is divided into five main compartments: white matter, grey matter, CSF, skull, and scalp. For the isotropic compartments, widely-used default conductivity values are assigned. For the brain tissues, we use the process of the effective medium approach (EMA) to estimate anisotropic conductivity tensors from diffusion-weighted imaging (DWI) data. The FEM electromagnetic calculations are performed by the DUNEuro library, integrated into Brainstorm and accessible with a user-friendly graphical interface. This integrated pipeline, with full tutorials and example data sets freely available on the Brainstorm website, gives the neuroscience community easy access to advanced tools for electromagnetic modeling using FEM.
△ Less
Submitted 2 November, 2020;
originally announced November 2020.
-
Modeling glioma invasion with anisotropy- and hypoxia-triggered motility enhancement: from subcellular dynamics to macroscopic PDEs with multiple taxis
Authors:
G. Corbin,
C. Engwer,
A. Klar,
J. Nieto,
J. Soler,
C. Surulescu,
M. Wenske
Abstract:
We deduce a model for glioma invasion making use of DTI data and accounting for the dynamics of brain tissue being actively degraded by tumor cells via excessive acidity production, but also according to the local orientation of tissue fibers. Our approach has a multiscale character: we start with a microscopic description of single cell dynamics including biochemical and/or biophysical effects of…
▽ More
We deduce a model for glioma invasion making use of DTI data and accounting for the dynamics of brain tissue being actively degraded by tumor cells via excessive acidity production, but also according to the local orientation of tissue fibers. Our approach has a multiscale character: we start with a microscopic description of single cell dynamics including biochemical and/or biophysical effects of the tumor microenvironment, translated on the one hand into cell stress and corresponding forces and on the other hand into receptor binding dynamics; these lead on the mesoscopic level to kinetic equations involving transport terms w.r.t. all kinetic variables and eventually, by appropriate upscaling, to a macroscopic reaction-diffusion equation for glioma density with multiple taxis, coupled to (integro-)differential equations characterizing the evolution of acidity and macro- and mesoscopic tissue. Our approach also allows for a switch between fast and slower moving regimes, %diffusion- and drift-dominated regimes, according to the local tissue anisotropy. We perform numerical simulations to investigate the behavior of solutions w.r.t. various scenarios of tissue dynamics and the dominance of each of the tactic terms, also suggesting how the model can be used to perform a numerical necrosis-based tumor grading or support radiotherapy planning by dose painting. We also provide a discussion about alternative ways of including cell level environmental influences in such multiscale modeling approach, ultimately leading in the macroscopic limit to (multiple) taxis.
△ Less
Submitted 22 June, 2020;
originally announced June 2020.
-
Estimating the extent of glioblastoma invasion
Authors:
Christian Engwer,
Michael Wenske
Abstract:
Glioblastoma Multiforme is a malignant brain tumor with poor prognosis. There have been numerous attempts to model the invasion of tumorous glioma cells via partial differential equations in the form of advection-diffusion-reaction equations. The patient-wise parametrisation of these models, and their validation via experimental data has been found to be difficult, as time sequence measurements ar…
▽ More
Glioblastoma Multiforme is a malignant brain tumor with poor prognosis. There have been numerous attempts to model the invasion of tumorous glioma cells via partial differential equations in the form of advection-diffusion-reaction equations. The patient-wise parametrisation of these models, and their validation via experimental data has been found to be difficult, as time sequence measurements are generally missing. Also the clinical interest lies in the actual (invisible) tumor extent for a particular MRI/DTI scan and not in a predictive estimate. Therefore we propose a stationalised approach to estimate the extent of glioblastoma (GBM) invasion at the time of a given MRI/DTI scan. The underlying dynamics can be derived from an instationary GBM model, falling into the wide class of advection-diffusion-reaction equations. The stationalisation is introduced via an analytical solution of the Fisher-KPP equation, the simplest model in the considered model class. We investigate the applicability in 1D and 2D, in the presence of inhomogeneous diffusion coefficients and on a real 3D DTI-dataset.
△ Less
Submitted 15 January, 2020;
originally announced January 2020.
-
Monotonicity considerations for stabilized DG cut cell schemes for the unsteady advection equation
Authors:
Florian Streitbürger,
Christian Engwer,
Sandra May,
Andreas Nüßing
Abstract:
For solving unsteady hyperbolic conservation laws on cut cell meshes, the so called small cell problem is a big issue: one would like to use a time step that is chosen with respect to the background mesh and use the same time step on the potentially arbitrarily small cut cells as well. For explicit time stepping schemes this leads to instabilities. In a recent preprint [arXiv:1906.05642], we propo…
▽ More
For solving unsteady hyperbolic conservation laws on cut cell meshes, the so called small cell problem is a big issue: one would like to use a time step that is chosen with respect to the background mesh and use the same time step on the potentially arbitrarily small cut cells as well. For explicit time stepping schemes this leads to instabilities. In a recent preprint [arXiv:1906.05642], we propose penalty terms for stabilizing a DG space discretization to overcome this issue for the unsteady linear advection equation. The usage of the proposed stabilization terms results in stable schemes of first and second order in one and two space dimensions. In one dimension, for piecewise constant data in space and explicit Euler in time, the stabilized scheme can even be shown to be monotone. In this contribution, we will examine the conditions for monotonicity in more detail.
△ Less
Submitted 26 December, 2019;
originally announced December 2019.
-
Strategies for the vectorized Block Conjugate Gradients method
Authors:
Nils-Arne Dreier,
Christian Engwer
Abstract:
Block Krylov methods have recently gained a lot of attraction. Due to their increased arithmetic intensity they offer a promising way to improve performance on modern hardware. Recently Frommer et al. presented a block Krylov framework that combines the advantages of block Krylov methods and data parallel methods. We review this framework and apply it on the Block Conjugate Gradients method,to sol…
▽ More
Block Krylov methods have recently gained a lot of attraction. Due to their increased arithmetic intensity they offer a promising way to improve performance on modern hardware. Recently Frommer et al. presented a block Krylov framework that combines the advantages of block Krylov methods and data parallel methods. We review this framework and apply it on the Block Conjugate Gradients method,to solve linear systems with multiple right hand sides. In this course we consider challenges that occur on modern hardware, like a limited memory bandwidth, the use of SIMD instructions and the communication overhead. We present a performance model to predict the efficiency of different Block CG variants and compare these with experimental numerical results.
△ Less
Submitted 26 December, 2019;
originally announced December 2019.
-
Dynamic and weighted stabilizations of the $L$-scheme applied to a phase-field model for fracture propagation
Authors:
Christian Engwer,
Iuliu Sorin Pop,
Thomas Wick
Abstract:
We consider a phase-field fracture propagation model, which consists of two (nonlinear) coupled partial differential equations. The first equation describes the displacement evolution, and the second is a smoothed indicator variable, describing the crack position. We propose an iterative scheme, the so-called $L$-scheme, with a dynamic update of the stabilization parameters during the iterations.…
▽ More
We consider a phase-field fracture propagation model, which consists of two (nonlinear) coupled partial differential equations. The first equation describes the displacement evolution, and the second is a smoothed indicator variable, describing the crack position. We propose an iterative scheme, the so-called $L$-scheme, with a dynamic update of the stabilization parameters during the iterations. Our algorithmic improvements are substantiated with two numerical tests. The dynamic adjustments of the stabilization parameters lead to a significant reduction of iteration numbers in comparison to constant stabilization values.
△ Less
Submitted 15 December, 2019;
originally announced December 2019.
-
Exa-Dune -- Flexible PDE Solvers, Numerical Methods and Applications
Authors:
Peter Bastian,
Mirco Altenbernd,
Nils-Arne Dreier,
Christian Engwer,
Jorrit Fahlke,
René Fritze,
Markus Geveler,
Dominik Göddeke,
Oleg Iliev,
Olaf Ippisch,
Jan Mohring,
Steffen Müthing,
Mario Ohlberger,
Dirk Ribbrock,
Nikolay Shegunov,
Stefan Turek
Abstract:
In the Exa-Dune project we have developed, implemented and optimised numerical algorithms and software for the scalable solution of partial differential equations (PDEs) on future exascale systems exhibiting a heterogeneous massively parallel architecture. In order to cope with the increased probability of hardware failures, one aim of the project was to add flexible, application-oriented resilien…
▽ More
In the Exa-Dune project we have developed, implemented and optimised numerical algorithms and software for the scalable solution of partial differential equations (PDEs) on future exascale systems exhibiting a heterogeneous massively parallel architecture. In order to cope with the increased probability of hardware failures, one aim of the project was to add flexible, application-oriented resilience capabilities into the framework. Continuous improvement of the underlying hardware-oriented numerical methods have included GPU-based sparse approximate inverses, matrix-free sum-factorisation for high-order discontinuous Galerkin discretisations as well as partially matrix-free preconditioners. On top of that, additional scalability is facilitated by exploiting massive coarse grained parallelism offered by multiscale and uncertainty quantification methods where we have focused on the adaptive choice of the coarse/fine scale and the overlap region as well as the combination of local reduced basis multiscale methods and the multilevel Monte-Carlo algorithm. Finally, some of the concepts are applied in a land-surface model including subsurface flow and surface runoff.
△ Less
Submitted 6 November, 2019; v1 submitted 4 November, 2019;
originally announced November 2019.
-
The DUNE Framework: Basic Concepts and Recent Developments
Authors:
Peter Bastian,
Markus Blatt,
Andreas Dedner,
Nils-Arne Dreier,
Christian Engwer,
René Fritze,
Carsten Gräser,
Christoph Grüninger,
Dominic Kempf,
Robert Klöfkorn,
Mario Ohlberger,
Oliver Sander
Abstract:
This paper presents the basic concepts and the module structure of the Distributed and Unified Numerics Environment and reflects on recent developments and general changes that happened since the release of the first Dune version in 2007 and the main papers describing that state [1, 2]. This discussion is accompanied with a description of various advanced features, such as coupling of domains and…
▽ More
This paper presents the basic concepts and the module structure of the Distributed and Unified Numerics Environment and reflects on recent developments and general changes that happened since the release of the first Dune version in 2007 and the main papers describing that state [1, 2]. This discussion is accompanied with a description of various advanced features, such as coupling of domains and cut cells, grid modifications such as adaptation and moving domains, high order discretizations and node level performance, non-smooth multigrid methods, and multiscale methods. A brief discussion on current and future development directions of the framework concludes the paper.
△ Less
Submitted 22 June, 2020; v1 submitted 30 September, 2019;
originally announced September 2019.
-
A stabilized DG cut cell method for discretizing the linear transport equation
Authors:
Christian Engwer,
Sandra May,
Andreas Nüßing,
Florian Streitbürger
Abstract:
We present new stabilization terms for solving the linear transport equation on a cut cell mesh using the discontinuous Galerkin (DG) method in two dimensions with piecewise linear polynomials. The goal is to allow for explicit time stepping schemes, despite the presence of cut cells. Using a method of lines approach, we start with a standard upwind DG discretization for the background mesh and ad…
▽ More
We present new stabilization terms for solving the linear transport equation on a cut cell mesh using the discontinuous Galerkin (DG) method in two dimensions with piecewise linear polynomials. The goal is to allow for explicit time stepping schemes, despite the presence of cut cells. Using a method of lines approach, we start with a standard upwind DG discretization for the background mesh and add penalty terms that stabilize the solution on small cut cells in a conservative way. Then, one can use explicit time stepping, even on cut cells, with a time step length that is appropriate for the background mesh. In one dimension, we show monotonicity of the proposed scheme for piecewise constant polynomials and total variation diminishing in the means stability for piecewise linear polynomials. We also present numerical results in one and two dimensions that support our theoretical findings.
△ Less
Submitted 13 June, 2019;
originally announced June 2019.
-
DUNEuro -- A software toolbox for forward modeling in bioelectromagnetism
Authors:
Sophie Schrader,
Andreas Westhoff,
Maria Carla Piastra,
Tuuli Miinalainen,
Sampsa Pursiainen,
Johannes Vorwerk,
Heinrich Brinck,
Carsten H. Wolters,
Christian Engwer
Abstract:
Accurate and efficient source analysis in electro- and magnetoencephalography using sophisticated realistic head geometries requires advanced numerical approaches. This paper presents DUNEuro, a free and open source C++ software toolbox for forward modeling in bioelectromagnetism. Building upon the DUNE framework, it provides implementations of modern fitted and unfitted finite element methods to…
▽ More
Accurate and efficient source analysis in electro- and magnetoencephalography using sophisticated realistic head geometries requires advanced numerical approaches. This paper presents DUNEuro, a free and open source C++ software toolbox for forward modeling in bioelectromagnetism. Building upon the DUNE framework, it provides implementations of modern fitted and unfitted finite element methods to efficiently solve the forward problems in electro- and magnetoencephalography. The user can choose between a variety of different source models that are implemented. The software's aim is to provide interfaces that are extendible and easy-to-use. In order to enable a closer integration into existing analysis pipelines, interfaces to Python and Matlab are provided. The practical use is demonstrated by a source analysis example of somatosensory evoked potentials using a realistic six compartment head model.
△ Less
Submitted 14 January, 2021; v1 submitted 9 January, 2019;
originally announced January 2019.
-
Function space bases in the dune-functions module
Authors:
Christian Engwer,
Carsten Gräser,
Steffen Müthing,
Oliver Sander
Abstract:
The dune-functions Dune module provides interfaces for functions and function space bases. It forms one abstraction level above grids, shape functions, and linear algebra, and provides infrastructure for full discretization frameworks like dune-pdelab and dune-fem. This document describes the function space bases provided by dune-functions. These are based on an abstract description of bases for p…
▽ More
The dune-functions Dune module provides interfaces for functions and function space bases. It forms one abstraction level above grids, shape functions, and linear algebra, and provides infrastructure for full discretization frameworks like dune-pdelab and dune-fem. This document describes the function space bases provided by dune-functions. These are based on an abstract description of bases for product spaces as trees of simpler bases. From this description, many different numberings of degrees of freedom by multi-indices can be derived in a natural way. We describe the abstract concepts, document the programmer interface, and give a complete example program that solves the stationary Stokes equation using Taylor-Hood elements.
△ Less
Submitted 25 June, 2018;
originally announced June 2018.
-
A high-level C++ approach to manage local errors, asynchrony and faults in an MPI application
Authors:
Christian Engwer,
Mirco Altenbernd,
Nils-Arne Dreier,
Dominik Göddeke
Abstract:
C++ advocates exceptions as the preferred way to handle unexpected behaviour of an implementation in the code. This does not integrate well with the error handling of MPI, which more or less always results in program termination in case of MPI failures. In particular, a local C++ exception can currently lead to a deadlock due to unfinished communication requests on remote hosts. At the same time,…
▽ More
C++ advocates exceptions as the preferred way to handle unexpected behaviour of an implementation in the code. This does not integrate well with the error handling of MPI, which more or less always results in program termination in case of MPI failures. In particular, a local C++ exception can currently lead to a deadlock due to unfinished communication requests on remote hosts. At the same time, future MPI implementations are expected to include an API to continue computations even after a hard fault (node loss), i.e. the worst possible unexpected behaviour.
In this paper we present an approach that adds extended exception propagation support to C++ MPI programs. Our technique allows to propagate local exceptions to remote hosts to avoid deadlocks, and to map MPI failures on remote hosts to local exceptions. A use case of particular interest are asynchronous 'local failure local recovery' resilience approaches. Our prototype implementation uses MPI-3.0 features only. In addition we present a dedicated implementation, which integrates seamlessly with MPI-ULFM, i.e. the most prominent proposal for extending MPI towards fault tolerance.
Our implementation is available at https://gitlab.dune-project.org/christi/test-mpi-exceptions .
△ Less
Submitted 13 April, 2018; v1 submitted 12 April, 2018;
originally announced April 2018.
-
ArbiLoMod: Local Solution Spaces by Random Training in Electrodynamics
Authors:
Andreas Buhr,
Christian Engwer,
Mario Ohlberger,
Stephan Rave
Abstract:
The simulation method ArbiLoMod has the goal to provide users of Finite Element based simulation software with quick re-simulation after localized changes to the model under consideration. It generates a Reduced Order Model (ROM) for the full model without ever solving the full model. To this end, a localized variant of the Reduced Basis method is employed, solving only small localized problems in…
▽ More
The simulation method ArbiLoMod has the goal to provide users of Finite Element based simulation software with quick re-simulation after localized changes to the model under consideration. It generates a Reduced Order Model (ROM) for the full model without ever solving the full model. To this end, a localized variant of the Reduced Basis method is employed, solving only small localized problems in the generation of the reduced basis. The key to quick re-simulation lies in recycling most of the localized basis vectors after a localized model change. In this publication, ArbiLoMod's local training algorithm is analyzed numerically for the non-coercive problem of time harmonic Maxwell's equations in 2D, formulated in H(curl).
△ Less
Submitted 5 December, 2016; v1 submitted 20 June, 2016;
originally announced June 2016.
-
A Mixed Finite Element Method to Solve the EEG Forward Problem
Authors:
Johannes Vorwerk,
Christian Engwer,
Sampsa Pursiainen,
Carsten H. Wolters
Abstract:
Finite element methods have been shown to achieve high accuracies in numerically solving the EEG forward problem and they enable the realistic modeling of complex geometries and important conductive features such as anisotropic conductivities. To date, most of the presented approaches rely on the same underlying formulation, the continuous Galerkin (CG)-FEM. In this article, a novel approach to so…
▽ More
Finite element methods have been shown to achieve high accuracies in numerically solving the EEG forward problem and they enable the realistic modeling of complex geometries and important conductive features such as anisotropic conductivities. To date, most of the presented approaches rely on the same underlying formulation, the continuous Galerkin (CG)-FEM. In this article, a novel approach to solve the EEG forward problem based on a mixed finite element method (Mixed-FEM) is introduced. To obtain the Mixed-FEM formulation, the electric current is introduced as an additional unknown besides the electric potential. As a consequence of this derivation, the Mixed-FEM is, by construction, current preserving, in contrast to the CG-FEM. Consequently, a higher simulation accuracy can be achieved in certain scenarios, e.g., when the diameter of thin insulating structures, such as the skull, is in the range of the mesh resolution.
A theoretical derivation of the Mixed-FEM approach for EEG forward simulations is presented, and the algorithms implemented for solving the resulting equation systems are described. Subsequently, first evaluations in both sphere and realistic head models are presented, and the results are compared to previously introduced CG-FEM approaches. Additional visualizations are shown to illustrate the current preserving property of the Mixed-FEM.
Based on these results, it is concluded that the newly presented Mixed-FEM can at least complement and in some scenarios even outperform the established CG-FEM approaches, which motivates a further evaluation of the Mixed-FEM for applications in bioelectromagnetism.
△ Less
Submitted 23 August, 2017; v1 submitted 19 April, 2016;
originally announced April 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.
-
An unfitted discontinuous Galerkin scheme for conservation laws on evolving surfaces
Authors:
Christian Engwer,
Thomas Ranner,
Sebastian Westerheide
Abstract:
Motivated by considering partial differential equations arising from conservation laws posed on evolving surfaces, a new numerical method for an advection problem is developed and simple numerical tests are performed. The method is based on an unfitted discontinuous Galerkin approach where the surface is not explicitly tracked by the mesh which means the method is extremely flexible with respect t…
▽ More
Motivated by considering partial differential equations arising from conservation laws posed on evolving surfaces, a new numerical method for an advection problem is developed and simple numerical tests are performed. The method is based on an unfitted discontinuous Galerkin approach where the surface is not explicitly tracked by the mesh which means the method is extremely flexible with respect to geometry. Furthermore, the discontinuous Galerkin approach is well-suited to capture the advection driven by the evolution of the surface without the need for a space-time formulation, back-tracking trajectories or streamline diffusion. The method is illustrated by a one-dimensional example and numerical results are presented that show good convergence properties for a simple test problem.
△ Less
Submitted 2 February, 2016;
originally announced February 2016.
-
The Unfitted Discontinuous Galerkin Method for Solving the EEG Forward Problem
Authors:
Andreas Nüßing,
Carsten H. Wolters,
Heinrich Brinck,
Christian Engwer
Abstract:
Objective: The purpose of this study is to introduce and evaluate the unfitted discontinuous Galerkin finite element method (UDG-FEM) for solving the electroencephalography (EEG) forward problem. Methods: This new approach for source analysis does not use a geometry conforming volume triangulation, but instead uses a structured mesh that does not resolve the geometry. The geometry is described usi…
▽ More
Objective: The purpose of this study is to introduce and evaluate the unfitted discontinuous Galerkin finite element method (UDG-FEM) for solving the electroencephalography (EEG) forward problem. Methods: This new approach for source analysis does not use a geometry conforming volume triangulation, but instead uses a structured mesh that does not resolve the geometry. The geometry is described using level set functions and is incorporated implicitly in its mathematical formulation. As no triangulation is necessary, the complexity of a simulation pipeline and the need for manual interaction for patient specific simulations can be reduced and is comparable with that of the FEM for hexahedral meshes. In addition, it maintains conservation laws on a discrete level. Here, we present the theory for UDG-FEM forward modeling, its verification using quasi-analytical solutions in multi-layer sphere models and an evaluation in a comparison with a discontinuous Galerkin (DG-FEM) method on hexahedral and on conforming tetrahedral meshes. We furthermore apply the UDG-FEM forward approach in a realistic head model simulation study. Results: The given results show convergence and indicate a good overall accuracy of the UDG-FEM approach. UDG-FEM performs comparable or even better than DG-FEM on a conforming tetrahedral mesh while providing a less complex simulation pipeline. When compared to DG-FEM on hexahedral meshes, an overall better accuracy is achieved. Conclusion: The UDG-FEM approach is an accurate, flexible and promising method to solve the EEG forward problem. Significance: This study shows the first application of the UDG-FEM approach to the EEG forward problem.
△ Less
Submitted 2 May, 2016; v1 submitted 28 January, 2016;
originally announced January 2016.
-
Geometric Integration Over Irregular Domains with topologic Guarantees
Authors:
Christian Engwer,
Andreas Nüßing
Abstract:
Implicitly described domains are a well established tool in the simulation of time dependent problems, e.g. using level-set methods. In order to solve partial differential equations on such domains, a range of numerical methods was developed, e.g. the Immersed Boundary method, Unfitted Finite Element or Unfitted discontinuous Galerkin methods, eXtended or Generalised Finite Element methods, just t…
▽ More
Implicitly described domains are a well established tool in the simulation of time dependent problems, e.g. using level-set methods. In order to solve partial differential equations on such domains, a range of numerical methods was developed, e.g. the Immersed Boundary method, Unfitted Finite Element or Unfitted discontinuous Galerkin methods, eXtended or Generalised Finite Element methods, just to name a few. Many of these methods involve integration over cut-cells or their boundaries, as they are described by sub-domains of the original level-set mesh. We present a new algorithm to geometrically evaluate the integrals over domains described by a first-order, conforming level-set function. The integration is based on a polyhedral reconstruction of the implicit geometry, following the concepts of the Marching Cubes algorithm. The algorithm preserves various topological properties of the implicit geometry in its polyhedral reconstruction, making it suitable for Finite Element computations. Numerical experiments show second order accuracy of the integration. An implementation of the algorithm is available as free software, which allows for an easy incorporation into other projects. The software is in productive use within the DUNE framework.
△ Less
Submitted 14 January, 2016;
originally announced January 2016.
-
ArbiLoMod, a Simulation Technique Designed for Arbitrary Local Modifications
Authors:
Andreas Buhr,
Christian Engwer,
Mario Ohlberger,
Stephan Rave
Abstract:
Engineers manually optimizing a structure using Finite Element based simulation software often employ an iterative approach where in each iteration they change the structure slightly and resimulate. Standard Finite Element based simulation software is usually not well suited for this workflow, as it restarts in each iteration, even for tiny changes. In settings with complex local microstructure, w…
▽ More
Engineers manually optimizing a structure using Finite Element based simulation software often employ an iterative approach where in each iteration they change the structure slightly and resimulate. Standard Finite Element based simulation software is usually not well suited for this workflow, as it restarts in each iteration, even for tiny changes. In settings with complex local microstructure, where a fine mesh is required to capture the geometric detail, localized model reduction can improve this workflow. To this end, we introduce ArbiLoMod, a method which allows fast recomputation after arbitrary local modifications. It employs a domain decomposition and a localized form of the Reduced Basis Method for model order reduction. It assumes that the reduced basis on many of the unchanged domains can be reused after a localized change. The reduced model is adapted when necessary, steered by a localized error indicator. The global error introduced by the model order reduction is controlled by a robust and efficient localized a posteriori error estimator, certifying the quality of the result. We demonstrate ArbiLoMod for a coercive, parameterized example with changing structure.
△ Less
Submitted 29 November, 2016; v1 submitted 24 December, 2015;
originally announced December 2015.
-
The interface for functions in the dune-functions module
Authors:
Christian Engwer,
Carsten Gräser,
Steffen Müthing,
Oliver Sander
Abstract:
The dune-functions dune module introduces a new programmer interface for discrete and non-discrete functions. Unlike the previous interfaces considered in the existing dune modules, it is based on overloading operator(), and returning values by-value. This makes user code much more readable, and allows the incorporation of newer C++ features such as lambda expressions. Run-time polymorphism is imp…
▽ More
The dune-functions dune module introduces a new programmer interface for discrete and non-discrete functions. Unlike the previous interfaces considered in the existing dune modules, it is based on overloading operator(), and returning values by-value. This makes user code much more readable, and allows the incorporation of newer C++ features such as lambda expressions. Run-time polymorphism is implemented not by inheritance, but by type erasure, generalizing the ideas of the std::function class from the C++11 standard library. We describe the new interface, show its possibilities, and measure the performance impact of type erasure and return-by-value.
△ Less
Submitted 18 December, 2015;
originally announced December 2015.
-
A discontinuous Galerkin Method for the EEG Forward Problem using the Subtraction Approach
Authors:
Christian Engwer,
Johannes Vorwerk,
Jakob Ludewig,
Carsten H. Wolters
Abstract:
In order to perform electroencephalography (EEG) source reconstruction, i.e., to localize the sources underlying a measured EEG, the electric potential distribution at the electrodes generated by a dipolar current source in the brain has to be simulated, which is the so-called EEG forward problem. To solve it accurately, it is necessary to apply numerical methods that are able to take the individu…
▽ More
In order to perform electroencephalography (EEG) source reconstruction, i.e., to localize the sources underlying a measured EEG, the electric potential distribution at the electrodes generated by a dipolar current source in the brain has to be simulated, which is the so-called EEG forward problem. To solve it accurately, it is necessary to apply numerical methods that are able to take the individual geometry and conductivity distribution of the subject's head into account. In this context, the finite element method (FEM) has shown high numerical accuracy with the possibility to model complex geometries and conductive features, e.g., white matter conductivity anisotropy. In this article, we introduce and analyze the application of a discontinuous Galerkin (DG) method, a finite element method that includes features of the finite volume framework, to the EEG forward problem. The DG-FEM approach fulfills the conservation property of electric charge also in the discrete case, making it attractive for a variety of applications. Furthermore, as we show, this approach can alleviate modeling inaccuracies that might occur in head geometries when using classical FE methods, e.g., so-called "skull leakage effects", which may occur in areas where the thickness of the skull is in the range of the mesh resolution. Therefore, we derive a DG formulation of the FEM subtraction approach for the EEG forward problem and present numerical results that highlight the advantageous features and the potential benefits of the proposed approach.
△ Less
Submitted 15 November, 2016; v1 submitted 16 November, 2015;
originally announced November 2015.
-
A numerically stable a posteriori error estimator for reduced basis approximations of elliptic equations
Authors:
Andreas Buhr,
Christian Engwer,
Mario Ohlberger,
Stephan Rave
Abstract:
The Reduced Basis (RB) method is a well established method for the model order reduction of problems formulated as parametrized partial differential equations. One crucial requirement for the application of RB schemes is the availability of an a posteriori error estimator to reliably estimate the error introduced by the reduction process. However, straightforward implementations of standard residu…
▽ More
The Reduced Basis (RB) method is a well established method for the model order reduction of problems formulated as parametrized partial differential equations. One crucial requirement for the application of RB schemes is the availability of an a posteriori error estimator to reliably estimate the error introduced by the reduction process. However, straightforward implementations of standard residual based estimators show poor numerical stability, rendering them unusable if high accuracy is required. In this work we propose a new algorithm based on representing the residual with respect to a dedicated orthonormal basis, which is both easy to implement and requires little additional computational overhead. A numerical example is given to demonstrate the performance of the proposed algorithm.
△ Less
Submitted 30 July, 2014;
originally announced July 2014.