Research WG Applications of PDEs - Prof. Dr. Christian Engwer

  • Overview and focus

    Many of our applications origin from porous media or biological systems, which exhibit very different kinds of complexity. Complexity can origin from a complicated geometric shapes, which poses a challenge for the numerical solution of PDEs in the complex shaped domain. The other kind of complexity is complexity of the system itself, due to complex couplings between different physical, biological & chemical processes.

    Complex Geometries

    © TODO

     When solving of PDEs on time dependent domains or domains with a complex shape, classic Finite Element Methods pose many problems regarding the construction of the mesh. During the last decade a range of different methods have been developed to decouple the construction of a finite element mesh, i.e. the finite element discretization, from the geometrical details of the domain.

    One approach, our group is working on, is the Unfitted Discontinuous Galerkin method. It offers the possibility to compute simulations with a fine structures on a relatively coarse mesh and was used for the solution of elliptic, parabolic and hyperbolic problems. Using the UDG approach it is easily possible to run simulations directly on image data, e.g. micro-CT images, or to combine it with level-set or phase-field methods to handle moving interfaces.

    Multi-Physics Problems

    © TODO

    In the course of multi-physics applications the efficient coupling of different PDEs on different sub-domains is getting more important. We are working on different aspects of domain decomposition methods and their implementation, either for parallelization and preconditioning, or for the coupling in a multi-physics setting. The latter also includes heterogenous coupling of sub-domains of different dimension.

    Efficient PDE Software

    Dune
    Dune
    © TODO

    We set high value on the development of efficient FEM software. Reusability and interoperability of and with existing software are very important. In this course we are actively participating in the development of the C++ FEM framework DUNE.

    Programming with C++ and using generic programming techniques, allows us to use fine grained interfaces and still employ optimization techniques, like inlining and loop-unrolling. This is the basis for sustainable and efficient software development.

    High Performance Computing

    © TODO

    The speed of a single processor stopped growing in the last years, instead modern chips include many cores to increase the performance. At the same time the architecture of high performance computers like the BlueGene is changing, they include acceleration processors which leaves us with a heterogeneous hardware system. Modern scientific software must cope with these changing requirements. As it is too much a burden to expect scientist to rewrite their code for each new hardware, the software design and the numerical algorithms must be adopted in a way that allows us to port our software with as small work as possible, while still retaining a reasonable performance boost.

  •  

    Latest Publications

    • Kleikamp, Hendrik. . Application of an adaptive model hierarchy to parametrized optimal control problems arXiv. 1st Ed. . doi: 10.48550/arXiv.2402.10708. [submitted / under review]
    • Schembera, Björn; Wübbeling, Frank; Kleikamp, Hendrik; Biedinger, Christine; Fiedler, Jochen; Reidelbach, Marco; Shehu, Aurela; Schmidt, Burkhard; Koprucki, Thomas; Iglezakis, Dorothea; Göddeke, Dominik. . ‘Ontologies for Models and Algorithms in Applied Mathematics and Related Disciplines.’ In Communications in Computer and Information Science.: Springer. [accepted / in Press (not yet published)]

    • Engwer, Christian; Ohlberger, Mario; Renelt, Lukas. . Model order reduction of an ultraweak and optimally stable variational formulation for parametrized reactive transport problems arXiv. doi: 10.48550/arXiv.2310.19674. [submitted / under review]
    • Renelt, Lukas; Ohlberger, Mario; Engwer, Christian. . ‘An optimally stable approximation of reactive transport using discrete test and infinite trial spaces.’ In Finite Volumes for Complex Applications X—Volume 2, Hyperbolic and Related Problems, edited by Franck, Emmanuel; Fuhrmann, Jürgen;, Michel-Dansac, Victor; Navoret, Laurent, 289–298. Cham: Springer. doi: 10.1007/978-3-031-40860-1_30.
    • Kleikamp, Hendrik; Lazar, Martin; Molinari, Cesare. . Be greedy and learn: efficient and certified algorithms for parametrized optimal control problems arXiv. doi: 10.48550/arXiv.2307.15590. [submitted / under review]
    • Haasdonk B, Kleikamp H, Ohlberger M, Schindler F, Wenzel T. . ‘A new certified hierarchical and adaptive RB-ML-ROM surrogate model for parametrized PDEs.’ SIAM Journal on Scientific Computing 45, No. 3: A1039. doi: 10.1137/22M1493318.
    • Medani, Takfarinas; Garcia-Prieto, Juan; Tadel, Francois; Antonakakis, Marios; Erdbrügger, Tim; Höltershinken, Malte; Mead, Wayne; Schrader, Sophie; Joshi, Anand; Engwer, Christian; Wolters, Carsten H.; Mosher, John C.; Leahy, Richard M. . ‘Brainstorm-DUNEuro: An integrated and user-friendly Finite Element Method for modeling electromagnetic brain activity.’ NeuroImage 267: 119851. doi: 10.1016/j.neuroimage.2022.119851.
    • Erdbrügger, T.; Westhoff, A.; Höltershinken, M.; Radecke, J.-O.; Buschermöhle, Y.; Buyx, A.; Wallois, F.; Pursiainen, S.; Gross, J.; Lencer, R.; Engwer, C.; Wolters, C.H. . ‘CutFEM forward modeling for EEG source analysis.’ Frontiers in Human Neuroscience 17: 1216758. doi: 10.3389/fnhum.2023.1216758.
    • Wenzel, Tizian; Haasdonk, Bernard; Kleikamp, Hendrik; Ohlberger, Mario; Schindler, Felix. . Application of Deep Kernel Models for Certified and Adaptive RB-ML-ROM Surrogate Modeling arXiv. doi: 10.48550/ARXIV.2302.14526. [accepted / in Press (not yet published)]

    • Keil T, Kleikamp H, Lorentzen R, Oguntola M, Ohlberger M. . ‘Adaptive machine learning based surrogate modeling to accelerate PDE-constrained optimization in enhanced oil recovery.’ Advances in Computational Mathematics 2022, No. 48: 73. doi: 10.1007/s10444-022-09981-z.

  •  

    Supervised Doctoral Studies

    Wenske, MichaelData-aware methods for the simulation of glioblastoma multiforme
    An unfitted discontinuous Galerkin scheme for a phase-field approximation of pressurized fractures
    Towards Automatic and Reliable Localized Model Order Reduction. Local Training, a Posteriori Error Estimation and Online Enrichment.
    Piastra, Maria CarlaNew Finite Element Methods for MEG and combined EEG/MEG Forward Problem
    Fitted and unfitted finite element methods for solving the EEG forward problem
    Emken, NatalieA coupled bulk-surface reaction-diffusion-advection model for cell polarization
    Vorwerk, JohannesNew Finite Element Methods to Solve the EEG/MEG Forward Problem
    Vorwerk, JohannesNew Finite Element Methods to solve the EEG/MEG Forward Problem
    DTI data based multiscale modelling and simulation of glioma growth