BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Mario Ohlberger appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
We present an integrated approach for the use of simulated data from full order discretization as well as projection-based Reduced Basis reduced order models for the training of machine learning approaches, in particular Kernel Methods, in order to achieve fast, reliable predictive models for the chemical conversion rate in reactive flows with varying transport regimes.
In this contribution, we develop an efficient surrogate modeling framework for simulation-based optimization of enhanced oil recovery, where we particularly focus on polymer flooding. The computational approach is based on an adaptive training procedure of a neural network that directly approximates an input-output map of the underlying PDE-constrained optimization problem. The training process thereby focuses on the construction of an accurate surrogate model solely related to the optimization path of an outer iterative optimization loop. True evaluations of the objective function are used to finally obtain certified results. Numerical experiments are given to evaluate the accuracy and efficiency of the approach for a heterogeneous five-spot benchmark problem.
We consider machine-learning of time-dependent quantities of interest derived from solution trajectories of parabolic partial differential equations. For large-scale or long-time integration scenarios, where using a full order model (FOM) to generate sufficient training data is computationally prohibitive, we propose an adaptive hierarchy of intermediate Reduced Basis reduced order models (ROM) to augment the FOM training data by certified ROM training data required to fit a kernel model.
In many applications, for instance when describing dynamics of fluids or gases, hyperbolic conservation laws arise naturally in the modeling of conserved quantities of a system, like mass or energy. These types of equations exhibit highly nonlinear behaviors like shock formation or shock interaction. In the case of parametrized hyperbolic equations, where, for instance, varying transport velocities are considered, these nonlinearities and strong transport effects result in a highly nonlinear solution manifold. This solution manifold cannot be approximated properly by linear subspaces. To this end, nonlinear approaches for model order reduction of hyperbolic conservation laws are required. We propose a new method for nonlinear model order reduction that is especially well-suited for hyperbolic equations with discontinuous solutions. The approach is based on a space-time discretization and employs diffeomorphic transformations of the underlying space-time domain to align the discontinuities. To derive a reduced model for the diffeomorphisms, the Lie group structure of the diffeomorphism group is used to associate diffeomorphisms with corresponding velocity fields via the exponential map. In the linear space of velocity fields, standard model order reduction techniques, such as proper orthogonal decomposition, can be applied to extract a reduced subspace. For a parametrized Burgers' equation with two merging shocks, numerical experiments show the potential of the approach.
Projection based model order reduction has become a mature technique for simulation of large classes of parameterized systems. However, several challenges remain for problems where the solution manifold of the parameterized system cannot be well approximated by linear subspaces. While the online efficiency of these model reduction methods is very convincing for problems with a rapid decay of the Kolmogorov n-width, there are still major drawbacks and limitations. Most importantly, the construction of the reduced system in the offline phase is extremely CPU-time and memory consuming for large scale and multi scale systems. For practical applications, it is thus necessary to derive model reduction techniques that do not rely on a classical offline/online splitting but allow for more flexibility in the usage of computational resources. A promising approach with this respect is model reduction with adaptive enrichment. In this contribution we investigate Petrov-Galerkin based model reduction with adaptive basis enrichment within a Trust Region approach for the solution of multi scale and large scale PDE constrained parameter optimization.
In this contribution we propose and rigorously analyze new variants of adaptive Trust-Region methods for parameter optimization with PDE constraints and bilateral parameter constraints. The approach employs successively enriched Reduced Basis surrogate models that are constructed during the outer optimization loop and used as model function for the Trust-Region method. Each Trust-Region sub-problem is solved with the projected BFGS method. Moreover, we propose a non-conforming dual (NCD) approach to improve the standard RB approximation of the optimality system. Rigorous improved a posteriori error bounds are derived and used to prove convergence of the resulting NCD-corrected adaptive Trust-Region Reduced Basis algorithm. Numerical experiments demonstrate that this approach enables to reduce the computational demand for large scale or multi-scale PDE constrained optimization problems significantly.
In this contribution we device and analyze improved variants of the non-conforming dual approach for trust-region reduced basis (TR-RB) approximation of PDE-constrained parameter optimization that has recently been introduced in [Keil et al.. A non-conforming dual approach for adaptive Trust-Region Reduced Basis approximation of PDE-constrained optimization. arXiv:2006.09297, 2020]. The proposed methods use model order reduction techniques for parametrized PDEs to significantly reduce the computational demand of parameter optimization with PDE constraints in the context of large-scale or multi-scale applications. The adaptive TR approach allows to localize the reduction with respect to the parameter space along the path of optimization without wasting unnecessary resources in an offline phase. The improved variants employ projected Newton methods to solve the local optimization problems within each TR step to benefit from high convergence rates. This implies new strategies in constructing the RB spaces, together with an estimate for the approximation of the hessian. Moreover, we present a new proof of convergence of the TR-RB method based on infinite-dimensional arguments, not restricted to the particular case of an RB approximation and provide an a posteriori error estimate for the approximation of the optimal parameter. Numerical experiments demonstrate the efficiency of the proposed methods.
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.
In this paper, we suggest a new Heterogeneous Multiscale Method (HMM) for the Helmholtz equation with high contrast. The method is constructed for a setting as in Bouchitté and Felbacq (C.R. Math. Acad. Sci. Paris 339(5):377--382, 2004), where the high contrast in the parameter leads to unusual effective parameters in the homogenized equation. We revisit existing homogenization approaches for this special setting and analyze the stability of the two-scale solution with respect to the wavenumber and the data. This includes a new stability result for solutions to the Helmholtz equation with discontinuous diffusion matrix. The HMM is defined as direct discretization of the two-scale limit equation. With this approach we are able to show quasi-optimality and an a priori error estimate under a resolution condition that inherits its dependence on the wavenumber from the stability constant for the analytical problem. Numerical experiments confirm our theoretical convergence results and examine the resolution condition. Moreover, the numerical simulation gives a good insight and explanation of the physical phenomenon of frequency band gaps.
In this article we introduce a procedure, which allows to recover the potentially very good approximation properties of tensor-based model reduction procedures for the solution of partial differential equations in the presence of interfaces or strong gradients in the solution which are skewed with respect to the coordinate axes. The two key ideas are the location of the interface either by solving a lower-dimensional partial differential equation or by using data functions and the subsequent removal of the interface of the solution by choosing the determined interface as the lifting function of the Dirichlet boundary conditions. We demonstrate in numerical experiments for linear elliptic equations and the reduced basis-hierarchical model reduction approach that the proposed procedure locates the interface well and yields a significantly improved convergence behavior even in the case when we only consider an approximation of the interface.
In this paper we introduce a new hierarchical model reduction framework for the Fokker-Planck equation. We reduce the dimension of the equation by a truncated basis expansion in the velocity variable, obtaining a hyperbolic system of equations in space and time. Unlike former methods like the Legendre moment models, the new framework generates a suitable problem-dependent basis of the reduced velocity space that mimics the shape of the solution in the velocity variable. To that end, we adapt the framework of [M. Ohlberger and K. Smetana. A dimensional reduction approach based on the application of reduced basis methods in the framework of hierarchical model reduction. SIAM J. Sci. Comput., 36(2):A714-A736, 2014] and derive initially a parametrized elliptic partial differential equation (PDE) in the velocity variable. Then, we apply ideas of the Reduced Basis method to develop a greedy-algorithm that selects the basis from solutions of the parametrized PDE. Numerical experiments demonstrate the potential of this new method.
Parametric model order reduction using reduced basis methods can be an effective tool for obtaining quickly solvable reduced order models of parametrized partial differential equation problems. With speedups that can reach several orders of magnitude, reduced basis methods enable high fidelity real-time simulations of complex systems and dramatically reduce the computational costs in many-query applications. In this contribution we analyze the methodology, mainly focussing on the theoretical aspects of the approach. In particular we discuss what is known about the convergence properties of these methods: when they succeed and when they are bound to fail. Moreover, we highlight some recent approaches employing nonlinear approximation techniques which aim to overcome the current limitations of reduced basis methods.
In this paper, we suggest a new heterogeneous multiscale method (HMM) for the time-harmonic Maxwell equations in locally periodic media. The method is constructed by using a divergence-regularization in one of the cell problems. This allows us to introduce fine-scale correctors that are not subject to a cumbersome divergence-free constraint and which can hence easily be implemented. To analyze the method, we first revisit classical homogenization theory for time-harmonic Maxwell equations and derive a new homogenization result that makes use of the divergence-regularization in the two-scale homogenized equation. We then show that the HMM is equivalent to a discretization of this equation. In particular, writing both problems in a fully coupled two-scale formulation is the crucial starting point for a corresponding numerical analysis of the method. With this approach we are able to prove rigorous a priori error estimates in the $\mathbf{H}(\mbox{curl})$- and the $H^{-1}$-norm and we derive reliable and efficient localized residual-based a posteriori error estimates.
The cross gramian matrix is a tool for model reduction and system identification, but it is only computable for square control systems. For symmetric systems the cross gramian possesses a useful relation to the system's associated Hankel singular values. Yet, many real-life models are neither square nor symmetric. In this work, concepts from decentralized control are used to approximate a cross gramian for non-symmetric and non-square systems. To illustrate this new non-symmetric cross gramian, it is applied in the context of model order reduction.
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.
In this contribution we present the first formulation of a heterogeneous multiscale method for an incompressible immiscible two-phase flow system with degenerate permeabilities. The method is in a general formulation which includes oversampling. We do not specify the discretization of the derived macroscopic equation, but we give two examples of possible realizations, suggesting a finite element solver for the fine scale and a vertex centered finite volume method for the effective coarse scale equations. Assuming periodicity, we show that the method is equivalent to a discretization of the homogenized equation. We provide an a-posteriori estimate for the error between the homogenized solutions of the pressure and saturation equations and the corresponding HMM approximations. The error estimate is based on the results recently achieved in [C. Canc{è}s, I. S. Pop, and M. Vohral\'ık. An a posteriori error estimate for vertex-centered finite volume discretizations of immiscible incompressible two-phase flow. Math. Comp., 2014].
This work introduces the empirical cross gramian for multiple-input-multiple-output systems. The cross gramian is a tool for reducing the state space of control systems, which conjoins controllability and observability information into a single matrix and does not require balancing. Its empirical gramian variant extends the application of the cross gramian to nonlinear systems. Furthermore, for parametrized systems, the empirical gramians can also be utilized for sensitivity analysis or parameter identification and thus for parameter reduction. This work also introduces the empirical joint gramian, which is derived from the empirical cross gramian. The joint gramian not only allows a reduction of the parameter space, but also the combined state and parameter space reduction, which is tested on a linear and a nonlinear control system. Controllability- and observability-based combined reduction methods are also presented, which are benchmarked against the joint gramian.
In this contribution we present an accelerated optimization-based approach for combined state and parameter reduction of a parametrized linear control system which is then used as a surrogate model in a Bayesian inverse setting. Following the basic ideas presented in [Lieberman, Willcox, Ghattas. Parameter and state model reduction for large-scale statistical inverse settings, SIAM J. Sci. Comput., 32(5):2523-2542, 2010], our approach is based on a generalized data-driven optimization functional in the construction process of the surrogate model and the usage of a trust-region-type solution strategy that results in an additional speed-up of the overall method. In principal, the model reduction procedure is based on the offline construction of appropriate low-dimensional state and parameter spaces and an online inversion step based on the resulting surrogate model that is obtained through projection of the underlying control system onto the reduced spaces. The generalization and enhancements presented in this work are shown to decrease overall computational time and increase accuracy of the reduced order model and thus allow an application to extreme-scale problems. Numerical experiments for a generic model and a fMRI connectivity model are presented in order to compare the computational efficiency of our improved method with the original approach.
The compact Discontinuous Galerkin 2 (CDG2) method was successfully tested for elliptic problems, scalar convection-diffusion equations and compressible Navier-Stokes equations. In this paper we use the newly developed DG method to solve a mathematical model for early stages of atherosclerotic plaque formation. Atherosclerotic plaque is mainly formed by accumulation of lipid-laden cells in the arterial walls which leads to a heart attack in case the artery is occluded or a thrombus is built through a rupture of the plaque. After describing a mathematical model and the discretization scheme, we present some benchmark tests comparing the CDG2 method to other commonly used DG methods. Furthermore, we take parallelization and higher order discretization schemes into account.
In this work, we propose a novel model order reduction approach for two-phase flow in porous media by introducing a formulation in which the mobility, which realizes the coupling between phase saturations and phase pressures, is regarded as a parameter to the pressure equation. Using this formulation, we introduce the Localized Reduced Basis Multiscale method to obtain a low-dimensional surrogate of the high-dimensional pressure equation. By applying ideas from model order reduction for parametrized partial differential equations, we are able to split the computational effort for solving the pressure equation into a costly offline step that is performed only once and an inexpensive online step that is carried out in every time step of the two-phase flow simulation, which is thereby largely accelerated. Usage of elements from numerical multiscale methods allows us to displace the computational intensity between the offline and online step to reach an ideal runtime at acceptable error increase for the two-phase flow simulation.
A common approach in model reduction is balanced truncation, which is based on gramian matrices classifiying certain attributes of states or parameters of a given dynamic system. Initially restricted to linear systems, the empirical gramians not only extended this concept to nonlinear systems, but also provide a uniform computational method. This work introduces a unified software framework supplying routines for six types of empirical gramians. The gramian types will be discussed and applied in a model reduction framework for multiple-input-multiple-output (MIMO) systems.
We recently introduced the joint gramian for combined state and parameter reduction [C. Himpe and M. Ohlberger. Cross-Gramian Based Combined State and Parameter Reduction for Large-Scale Control Systems. arXiv:1302.0634, 2013], which is applied in this work to reduce a parametrized linear time-varying control system modeling a hyperbolic network. The reduction encompasses the dimension of nodes and parameters of the underlying control system. Networks with a hyperbolic structure have many applications as models for large-scale systems. A prominent example is the brain, for which a network structure of the various regions is often assumed to model propagation of information. Networks with many nodes, and parametrized, uncertain or even unknown connectivity require many and individually computationally costly simulations. The presented model order reduction enables vast simulations of surrogate networks exhibiting almost the same dynamics with a small error compared to full order model.
We analytically and numerically analyze groundwater flow in a homogeneous soil described by the Richards equation, coupled to surface water represented by a set of ordinary differential equations (ODE's) on parts of the domain boundary, and with nonlinear outflow conditions of Signorini's type. The coupling of the partial differential equation (PDE) and the ODE's is given by nonlinear Robin boundary conditions. This article provides two major new contributions regarding these infiltration conditions. First, an existence result for the continuous coupled problem is established with the help of a regularization technique. Second, we analyze and validate a solver-friendly discretization of the coupled problem based on an implicit-explicit time discretization and on finite elements in space. The discretized PDE leads to convex spatial minimization problems which can be solved efficiently by monotone multigrid. Numerical experiments are provided using the DUNE numerics framework.