Researcher profile

Peter Benner

Peter Benner contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
35works
0followers
14topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

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

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

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

Building this graph slice

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

Published work

35 published item(s)

preprint2025arXiv

Reduced-Order Inference with Structure-Preserving Parametrization for Bending and Rotating Systems

Mechanical systems are often characterized only by their response to certain loads known from experiments or simulations. The obtained data can be used for various purposes: system analysis, design of mathematical models, or construction of reduced-order models for further simulations under different loading conditions. The use of data for reduced-order modeling is an important developing research direction, especially when the high-dimensional system operators are unknown and their low-dimensional approximation is required for accurate but fast simulations. Our goal is to obtain the low-dimensional surrogate model from the available input signal and deformation trajectory data, capturing the correct system behavior for the basic deformation cases, namely bending and rotation, which are present in almost every complex mechanical system. In this work, we propose a methodology to infer the system operators by solving a nonlinear unconstrained optimization problem. The methodology is based on the operator inference approach for second-order systems and includes a parametrization of the unknown operators that preserves their symmetric positive definite or skew-symmetric structure. We demonstrate the performance of the novel approach for three numerical examples that are used to simulate basic bending and rotating.

preprint2024arXiv

Generalized Quadratic Embeddings for Nonlinear Dynamics using Deep Learning

The engineering design process often relies on mathematical modeling that can describe the underlying dynamic behavior. In this work, we present a data-driven methodology for modeling the dynamics of nonlinear systems. To simplify this task, we aim to identify a coordinate transformation that allows us to represent the dynamics of nonlinear systems using a common, simple model structure. The advantage of a common simple model is that customized design tools developed for it can be applied to study a large variety of nonlinear systems. The simplest common model -- one can think of -- is linear, but linear systems often fall short in accurately capturing the complex dynamics of nonlinear systems. In this work, we propose using quadratic systems as the common structure, inspired by the lifting principle. According to this principle, smooth nonlinear systems can be expressed as quadratic systems in suitable coordinates without approximation errors. However, finding these coordinates solely from data is challenging. Here, we leverage deep learning to identify such lifted coordinates using only data, enabling a quadratic dynamical system to describe the system's dynamics. Additionally, we discuss the asymptotic stability of these quadratic dynamical systems. We illustrate the approach using data collected from various numerical examples, demonstrating its superior performance with the existing well-known techniques.

preprint2024arXiv

Guaranteed Stable Quadratic Models and their applications in SINDy and Operator Inference

Scientific machine learning for inferring dynamical systems combines data-driven modeling, physics-based modeling, and empirical knowledge. It plays an essential role in engineering design and digital twinning. In this work, we primarily focus on an operator inference methodology that builds dynamical models, preferably in low-dimension, with a prior hypothesis on the model structure, often determined by known physics or given by experts. Then, for inference, we aim to learn the operators of a model by setting up an appropriate optimization problem. One of the critical properties of dynamical systems is stability. However, this property is not guaranteed by the inferred models. In this work, we propose inference formulations to learn quadratic models, which are stable by design. Precisely, we discuss the parameterization of quadratic systems that are locally and globally stable. Moreover, for quadratic systems with no stable point yet bounded (e.g., chaotic Lorenz model), we discuss how to parameterize such bounded behaviors in the learning process. Using those parameterizations, we set up inference problems, which are then solved using a gradient-based optimization method. Furthermore, to avoid numerical derivatives and still learn continuous systems, we make use of an integral form of differential equations. We present several numerical examples, illustrating the preservation of stability and discussing its comparison with the existing state-of-the-art approach to infer operators. By means of numerical examples, we also demonstrate how the proposed methods are employed to discover governing equations and energy-preserving models.

preprint2023arXiv

Multi-fidelity error estimation accelerates greedy model reduction of complex dynamical systems

Model order reduction usually consists of two stages: the offline stage and the online stage. The offline stage is the expensive part that sometimes takes hours till the final reduced-order model is derived, especially when the original model is very large or complex. Once the reduced-order model is obtained, the online stage of querying the reduced-order model for simulation is very fast and often real-time capable. This work concerns a strategy to significantly speed up the offline stage of model order reduction for large and complex systems. In particular, it is successful in accelerating the greedy algorithm that is often used in the offline stage for reduced-order model construction. We propose multi-fidelity error estimators and replace the high-fidelity error estimator in the greedy algorithm. Consequently, the computational complexity at each iteration of the greedy algorithm is reduced and the algorithm converges more than 3 times faster without incurring noticeable accuracy loss.

preprint2022arXiv

A Structure-Preserving Divide-and-Conquer Method for Pseudosymmetric Matrices

We devise a spectral divide-and-conquer scheme for matrices that are self-adjoint with respect to a given indefinite scalar product (i.e. pseudosymmetic matrices). The pseudosymmetric structure of the matrix is preserved in the spectral division, such that the method can be applied recursively to achieve full diagonalization. The method is well-suited for structured matrices that come up in computational quantum physics and chemistry. In this application context, additional definiteness properties guarantee a convergence of the matrix sign function iteration within two steps when Zolotarev functions are used. The steps are easily parallelizable. Furthermore, it is shown that the matrix decouples into symmetric definite eigenvalue problems after just one step of spectral division.

preprint2022arXiv

An artificial neural network for surrogate modeling of stress fields in viscoplastic polycrystalline materials

The purpose of this work is the development of an artificial neural network (ANN) for surrogate modeling of the mechanical response of viscoplastic grain microstructures. To this end, a U-Net-based convolutional neural network (CNN) is trained to account for the history dependence of the material behavior. The training data take the form of numerical simulation results for the von Mises stress field under quasi-static tensile loading. The trained CNN (tCNN) can accurately reproduce both the average response as well as the local von Mises stress field. The tCNN calculates the von Mises stress field of grain microstructures not included in the training dataset about 500 times faster than its calculation based on the numerical solution with a spectral solver of the corresponding initial-boundary-value problem. The tCNN is also successfully applied to other types of microstructure morphologies (e.g., matrix-inclusion type topologies) and loading levels not contained in the training dataset.

preprint2022arXiv

Model reduction for second-order systems with inhomogeneous initial conditions

In this paper, we consider the problem of finding surrogate models for large-scale second-order linear time-invariant systems with inhomogeneous initial conditions. For this class of systems, the superposition principle allows us to decompose the system behavior into three independent components. The first behavior corresponds to the transfer between the input and output having zero initial conditions. In contrast, the other two correspond to the transfer between the initial position and the initial velocity conditions having zero input, respectively. Based on this superposition of systems, our goal is to propose model reduction schemes allowing to preserve the second-order structure in the surrogate models. To this aim, we introduce tailored second-order Gramians for each system component and compute them numerically, solving Lyapunov equations. As a consequence, two methodologies are proposed. The first one consists in reducing each of the components independently using a suitable balanced truncation procedure. The sum of these reduced systems provides an approximation of the original system. This methodology allows flexibility on the order of the reduced-order model. The second proposed methodology consists in extracting the dominant subspaces from the sum of Gramians to build the projection matrices leading to a surrogate model. Additionally, we discuss error bounds for the overall output approximation. Finally, the proposed methods are illustrated by means of benchmark problems.

preprint2022arXiv

Neural ODEs with Irregular and Noisy Data

Measurement noise is an integral part while collecting data of a physical process. Thus, noise removal is necessary to draw conclusions from these data, and it often becomes essential to construct dynamical models using these data. We discuss a methodology to learn differential equation(s) using noisy and irregular sampled measurements. In our methodology, the main innovation can be seen in the integration of deep neural networks with the neural ordinary differential equations (ODEs) approach. Precisely, we aim at learning a neural network that provides (approximately) an implicit representation of the data and an additional neural network that models the vector fields of the dependent variables. We combine these two networks by constraining using neural ODEs. The proposed framework to learn a model describing the vector field is highly effective under noisy measurements. The approach can handle scenarios where dependent variables are not available at the same temporal grid. Moreover, a particular structure, e.g., second-order with respect to time, can easily be incorporated. We demonstrate the effectiveness of the proposed method for learning models using data obtained from various differential equations and present a comparison with the neural ODE method that does not make any special treatment to noise.

preprint2022arXiv

On an integrated Krylov-ADI solver for large-scale Lyapunov equations

One of the most computationally expensive steps of the low-rank ADI method for large-scale Lyapunov equations is the solution of a shifted linear system at each iteration. We propose the use of the extended Krylov subspace method for this task. In particular, we illustrate how a single approximation space can be constructed to solve all the shifted linear systems needed to achieve a prescribed accuracy in terms of Lyapunov residual norm. Moreover, we show how to fully merge the two iterative procedures in order to obtain a novel, efficient implementation of the low-rank ADI method, for an important class of equations. Many state-of-the-art algorithms for the shift computation can be easily incorporated into our new scheme, as well. Several numerical results illustrate the potential of our novel procedure when compared to an implementation of the low-rank ADI method based on sparse direct solvers for the shifted linear systems.

preprint2022arXiv

The Hamiltonian Extended Krylov Subspace Method

An algorithm for constructing a $J$-orthogonal basis of the extended Krylov subspace $\mathcal{K}_{r,s}=\operatorname{range}\{u,Hu, H^2u,$ $ \ldots, $ $H^{2r-1}u, H^{-1}u, H^{-2}u, \ldots, H^{-2s}u\},$ where $H \in \mathbb{R}^{2n \times 2n}$ is a large (and sparse) Hamiltonian matrix is derived (for $r = s+1$ or $r=s$). Surprisingly, this allows for short recurrences involving at most five previously generated basis vectors. Projecting $H$ onto the subspace $\mathcal{K}_{r,s}$ yields a small Hamiltonian matrix. The resulting HEKS algorithm may be used in order to approximate $f(H)u$ where $f$ is a function which maps the Hamiltonian matrix $H$ to, e.g., a (skew-)Hamiltonian or symplectic matrix. Numerical experiments illustrate that approximating $f(H)u$ with the HEKS algorithm is competitive for some functions compared to the use of other (structure-preserving) Krylov subspace methods.

preprint2021arXiv

A low-rank solution method for Riccati equations with indefinite quadratic terms

Algebraic Riccati equations with indefinite quadratic terms play an important role in applications related to robust controller design. While there are many established approaches to solve these in case of small-scale dense coefficients, there is no approach available to compute solutions in the large-scale sparse setting. In this paper, we develop an iterative method to compute low-rank approximations of stabilizing solutions of large-scale sparse continuous-time algebraic Riccati equations with indefinite quadratic terms. We test the developed approach for dense examples in comparison to other established matrix equation solvers, and investigate the applicability and performance in large-scale sparse examples.

preprint2021arXiv

Reduced basis method for the nonlinear Poisson-Boltzmann equation regularized by the range-separated canonical tensor format

The Poisson-Boltzmann equation (PBE) is a fundamental implicit solvent continuum model for calculating the electrostatic potential of large ionic solvated biomolecules. However, its numerical solution encounters severe challenges arising from its strong singularity and nonlinearity. In [1,2], the effect of strong singularities was eliminated by applying the range-separated (RS) canonical tensor format [3,4] to construct a solution decomposition scheme for the PBE. The RS tensor format allows to derive a smooth approximation to the Dirac delta distribution in order to obtain a regularized PBE (RPBE) model. However, solving the RPBE is still computationally demanding due to its high dimension $\mathcal{N}$, where $\mathcal{N}$ is always in the millions. In this study, we propose to apply the reduced basis method (RBM) and the (discrete) empirical interpolation method ((D)EIM) to the RPBE in order to construct a reduced order model (ROM) of low dimension $N \ll \mathcal{N}$, whose solution accurately approximates the nonlinear RPBE. The long-range potential can be obtained by lifting the ROM solution back to the $\mathcal{N}$-space while the short-range potential is directly precomputed analytically, thanks to the RS tensor format. The sum of both provides the total electrostatic potential. The main computational benefit is the avoidance of computing the numerical approximation of the singular electrostatic potential. We demonstrate in the numerical experiments, the accuracy and efficacy of the reduced basis (RB) approximation to the nonlinear RPBE (NRPBE) solution and the corresponding computational savings over the classical nonlinear PBE (NPBE) as well as over the RBM being applied to the classical NPBE.

preprint2021arXiv

Robust output-feedback stabilization for incompressible flows using low-dimensional $\mathcal{H}_{\infty}$-controllers

Output-based controllers are known to be fragile with respect to model uncertainties. The standard $\mathcal{H}_{\infty}$-control theory provides a general approach to robust controller design based on the solution of the $\mathcal{H}_{\infty}$-Riccati equations. In view of stabilizing incompressible flows in simulations, two major challenges have to be addressed: the high-dimensional nature of the spatially discretized model and the differential-algebraic structure that comes with the incompressibility constraint. This work demonstrates the synthesis of low-dimensional robust controllers with guaranteed robustness margins for the stabilization of incompressible flow problems. The performance and the robustness of the reduced-order controller with respect to linearization and model reduction errors are investigated and illustrated in numerical examples.

preprint2020arXiv

A Low-rank Approach for Nonlinear Parameter-dependent Fluid-structure Interaction Problems

Parameter-dependent discretizations of linear fluid-structure interaction problems can be approached with low-rank methods. When discretizing with respect to a set of parameters, the resulting equations can be translated to a matrix equation since all operators involved are linear. If nonlinear FSI problems are considered, a direct translation to a matrix equation is not possible. We present a method that splits the parameter set into disjoint subsets and, on each subset, computes an approximation of the problem related to the upper median parameter by means of the Newton iteration. This approximation is then used as initial guess for one Newton step on a subset of problems.

preprint2020arXiv

A numerical comparison of solvers for large-scale, continuous-time algebraic Riccati equations and LQR problems

In this paper, we discuss numerical methods for solving large-scale continuous-time algebraic Riccati equations. These methods have been the focus of intensive research in recent years, and significant progress has been made in both the theoretical understanding and efficient implementation of various competing algorithms. There are several goals of this manuscript: first, to gather in one place an overview of different approaches for solving large-scale Riccati equations, and to point to the recent advances in each of them. Second, to analyze and compare the main computational ingredients of these algorithms, to detect their strong points and their potential bottlenecks. And finally, to compare the effective implementations of all methods on a set of relevant benchmark examples, giving an indication of their relative performance.

preprint2020arXiv

A rational Even-IRA algorithm for the solution of T-even polynomial eigenvalue problems

In this work we present a rational Krylov subspace method for solving real large-scale polynomial eigenvalue problems with T-even (that is, symmetric/skew-symmetric) structure. Our method is based on the Even-IRA algorithm. To preserve the structure, a sparse T-even linearization from the class of block minimal bases pencils is applied. Due to this linearization, the Krylov basis vectors can be computed in a cheap way. A rational decomposition is derived so that our method explicitly allows for changes of the shift during the iteration. This leads to a method that is able to compute parts of the spectrum of a T-even matrix polynomial in a fast and reliable way.

preprint2020arXiv

Automatic Decoupling and Index-aware Model-Order Reduction for Nonlinear Differential-Algebraic Equations

We extend the index-aware model-order reduction method to systems of nonlinear differential-algebraic equations with a special nonlinear term f(Ex), where E is a singular matrix. Such nonlinear differential-algebraic equations arise, for example, in the spatial discretization of the gas flow in pipeline networks. In practice, mathematical models of real-life processes pose challenges when used in numerical simulations, due to complexity and system size. Model-order reduction aims to eliminate this problem by generating reduced-order models that have lower computational cost to simulate, yet accurately represent the original large-scale system behavior. However, direct reduction and simulation of nonlinear differential-algebraic equations is difficult due to hidden constraints which affect the choice of numerical integration methods and model-order reduction techniques. We propose an extension of index-aware model-order reduction methods to nonlinear differential-algebraic equations without any kind of linearization. The proposed model-order reduction approach involves automatic decoupling of nonlinear differential-algebraic equations into nonlinear ordinary differential equations and algebraic equations. This allows applying standard model-order reduction techniques to both parts without worrying about the index. The same procedure can also be used to simulate nonlinear differential-algebraic equations using standard integration schemes. We illustrate the performance of our proposed method for nonlinear differential-algebraic equations arising from gas flow models in pipeline networks.

preprint2020arXiv

Clustering-Based Model Order Reduction for Nonlinear Network Systems

Clustering by projection has been proposed as a way to preserve network structure in linear multi-agent systems. Here, we extend this approach to a class of nonlinear network systems. Additionally, we generalize our clustering method which restores the network structure in an arbitrary reduced-order model obtained by projection. We demonstrate this method on a number of examples.

preprint2020arXiv

Data-Driven Learning of Reduced-order Dynamics for a Parametrized Shallow Water Equation

This paper discusses a non-intrusive data-driven model order reduction method that learns low-dimensional dynamical models for a parametrized shallow water equation. We consider the shallow water equation in non-traditional form (NTSWE). We focus on learning low-dimensional models in a non-intrusive way. That means, we assume not to have access to a discretized form of the NTSWE in any form. Instead, we have snapshots that are obtained using a black-box solver. Consequently, we aim at learning reduced-order models only from the snapshots. Precisely, a reduced-order model is learnt by solving an appropriate least-squares optimization problem in a low-dimensional subspace. Furthermore, we discuss computational challenges that particularly arise from the optimization problem being ill-conditioned. Moreover, we extend the non-intrusive model order reduction framework to a parametric case where we make use of the parameter dependency at the level of the partial differential equation. We illustrate the efficiency of the proposed non-intrusive method to construct reduced-order models for NTSWE and compare it with an intrusive method (proper orthogonal decomposition). We furthermore discuss the predictive capabilities of both models outside the range of the training data.

preprint2020arXiv

Data-Driven Snapshot Calibration via Monotonic Feature Matching

Snapshot matrices of hyperbolic equations have a slow singular value decay, resulting in inefficient reduced-order models. We develop on the idea of inducing a faster singular value decay by computing snapshots on a transformed spatial domain, or the so-called snapshot calibration/transformation. We are particularly interested in problems involving shock collision, shock rarefaction-fan collision, shock formation, etc. For such problems, we propose a realizable algorithm to compute the spatial transform using monotonic feature matching. We consider discontinuities and kinks as features, and by carefully partitioning the parameter domain, we ensure that the spatial transform has properties that are desirable both from a theoretical and an implementation standpoint. We use these properties to prove that our method results in a fast m-width decay of a so-called calibrated manifold. A crucial observation we make is that due to calibration, the m-width does not only depend on m but also on the accuracy of the full order model, which is in contrast to elliptic and parabolic problems that do not need calibration. The method we propose only requires the solution snapshots and not the underlying partial differential equation (PDE) and is therefore, data-driven. We perform several numerical experiments to demonstrate the effectiveness of our method.

preprint2020arXiv

Efficient and Accurate Algorithms for Solving the Bethe-Salpeter Eigenvalue Problem for Crystalline Systems

Optical properties of materials related to light absorption and scattering are explained by the excitation of electrons. The Bethe-Salpeter equation is the state-of-the-art approach to describe these processes from first principles (ab initio), i.e. without the need for empirical data in the model. To harness the predictive power of the equation, it is mapped to an eigenvalue problem via an appropriate discretization scheme. The eigenpairs of the resulting large, dense, structured matrix can be used to compute dielectric properties of the considered crystalline or molecular system. The matrix always shows a $2\times 2$ block structure. Additionally, certain definiteness properties typically hold. One form can be acquired for crystalline systems, another one is more general and can for example be used to study molecules. In this work, we present new theoretical results characterizing the structure of the two forms in the language of non-standard scalar products. These results enable us to develop a new perspective on the state-of-the-art solution approach for crystalline systems. This new viewpoint is used to develop two new methods for solving the eigenvalue problem. One requires less computational effort while providing the same degree of accuracy. The other one improves the expected accuracy, compared to methods currently in use, with a comparable performance. Both methods are well suited for high performance environments and only rely on basic numerical linear algebra building blocks.

preprint2020arXiv

GR decompositions and their relations to Cholesky-like factorizations

For a given matrix, we are interested in computing GR decompositions $A=GR$, where $G$ is an isometry with respect to given scalar products. The orthogonal QR decomposition is the representative for the Euclidian scalar product. For a signature matrix, a respective factorization is given as the hyperbolic QR decomposition. Considering a skew-symmetric matrix leads to the symplectic QR decomposition. The standard approach for computing GR decompositions is based on the successive elimination of subdiagonal matrix entries. For the hyperbolic and symplectic case, this approach does in general not lead to a satisfying numerical accuracy. An alternative approach computes the QR decomposition via a Cholesky factorization, but also has bad stability. It is improved by repeating the procedure a second time. In the same way, the hyperbolic and the symplectic QR decomposition are related to the $LDL^T$ and a skew-symmetric Cholesky-like factorization. We show that methods exploiting this connection can provide better numerical stability than elimination-based approaches.

preprint2020arXiv

High Performance Solution of Skew-symmetric Eigenvalue Problems with Applications in Solving the Bethe-Salpeter Eigenvalue Problem

We present a high-performance solver for dense skew-symmetric matrix eigenvalue problems. Our work is motivated by applications in computational quantum physics, where one solution approach to solve the so-called Bethe-Salpeter equation involves the solution of a large, dense, skew-symmetric eigenvalue problem. The computed eigenpairs can be used to compute the optical absorption spectrum of molecules and crystalline systems. One state-of-the art high-performance solver package for symmetric matrices is the ELPA (Eigenvalue SoLvers for Petascale Applications) library. We extend the methods available in ELPA to skew-symmetric matrices. This way, the presented solution method can benefit from the optimizations available in ELPA that make it a well-established, efficient and scalable library, such as GPU support. We compare performance and scalability of our method to the only available high-performance approach for skew-symmetric matrices, an indirect route involving complex arithmetic. In total, we achieve a performance that is up to 3.67 higher than the reference method using Intel's ScaLAPACK implementation. The runtime to solve the Bethe-Salpeter-Eigenvalue problem can be improved by a factor of 10. Our method is freely available in the current release of the ELPA library.

preprint2020arXiv

Low-Rank and Total Variation Regularization and Its Application to Image Recovery

In this paper, we study the problem of image recovery from given partial (corrupted) observations. Recovering an image using a low-rank model has been an active research area in data analysis and machine learning. But often, images are not only of low-rank but they also exhibit sparsity in a transformed space. In this work, we propose a new problem formulation in such a way that we seek to recover an image that is of low-rank and has sparsity in a transformed domain. We further discuss various non-convex non-smooth surrogates of the rank function, leading to a relaxed problem. Then, we present an efficient iterative scheme to solve the relaxed problem that essentially employs the (weighted) singular value thresholding at each iteration. Furthermore, we discuss the convergence properties of the proposed iterative method. We perform extensive experiments, showing that the proposed algorithm outperforms state-of-the-art methodologies in recovering images.

preprint2020arXiv

Low-rank Linear Fluid-structure Interaction Discretizations

Fluid-structure interaction models involve parameters that describe the solid and the fluid behavior. In simulations, there often is a need to vary these parameters to examine the behavior of a fluid-structure interaction model for different solids and different fluids. For instance, a shipping company wants to know how the material, a ship's hull is made of, interacts with fluids at different Reynolds and Strouhal numbers before the building process takes place. Also, the behavior of such models for solids with different properties is considered before the prototype phase. A parameter-dependent linear fluid-structure interaction discretization provides approximations for a bundle of different parameters at one step. Such a discretization with respect to different material parameters leads to a big block-diagonal system matrix that is equivalent to a matrix equation as discussed in [KressnerTobler 2011]. The unknown is then a matrix which can be approximated using a low-rank approach that represents the iterate by a tensor. This paper discusses a low-rank GMRES variant and a truncated variant of the Chebyshev iteration. Bounds for the error resulting from the truncation operations are derived. Numerical experiments show that such truncated methods applied to parameter-dependent discretizations provide approximations with relative residual norms smaller than $10^{-8}$ within a twentieth of the time used by individual standard approaches.

preprint2020arXiv

Matrix Equations, Sparse Solvers: M-M.E.S.S.-2.0.1 -- Philosophy, Features and Application for (Parametric) Model

Matrix equations are omnipresent in (numerical) linear algebra and systems theory. Especially in model order reduction (MOR) they play a key role in many balancing based reduction methods for linear dynamical systems. When these systems arise from spatial discretizations of evolutionary partial differential equations, their coefficient matrices are typically large and sparse. Moreover, the numbers of inputs and outputs of these systems are typically far smaller than the number of spatial degrees of freedom. Then, in many situations the solutions of the corresponding large-scale matrix equations are observed to have low (numerical) rank. This feature is exploited by M-M.E.S.S. to find successively larger low-rank factorizations approximating the solutions. This contribution describes the basic philosophy behind the implementation and the features of the package, as well as its application in the model order reduction of large-scale linear time-invariant (LTI) systems and parametric LTI systems.

preprint2020arXiv

On Error Estimation for Reduced-order Modeling of Linear Non-parametric and Parametric Systems

Motivated by a recently proposed error estimator for the transfer function of the reduced-order model of a given linear dynamical system, we further develop more theoretical results in this work. Furthermore, we propose several variants of the error estimator, and compare those variants with the existing ones both theoretically and numerically. It has been shown that some of the proposed error estimators perform better than or equally well as the existing ones. All the error estimators considered can be easily extended to estimate output error of reduced-order modeling for steady linear parametric systems.

preprint2020arXiv

On the Eigenvalue Distribution for a Beam with Attached Masses

We study a mathematical model of a hinged flexible beam with piezoelectric actuators and electromagnetic shaker in this paper. The shaker is modelled as a mass and spring system attached to the beam. To analyze free vibrations of this mechanical system, we consider the corresponding spectral problem for a fourth-order differential operator with interface conditions that characterize the shaker dynamics. The characteristic equation is studied analytically, and asymptotic estimates of eigenvalues are obtained. The eigenvalue distribution is also illustrated by numerical simulations under a realistic choice of mechanical parameters.

preprint2020arXiv

On the Solution of the Nonsymmetric T-Riccati Equation

The nonsymmetric T-Riccati equation is a quadratic matrix equation where the linear part corresponds to the so-called T-Sylvester or T-Lyapunov operator that has previously been studied in the literature. It has applications in macroeconomics and policy dynamics. So far, it presents an unexplored problem in numerical analysis, and both, theoretical results and computational methods, are lacking in the literature. In this paper we provide some sufficient conditions for the existence and uniqueness of a nonnegative minimal solution and its efficient computation is deeply analyzed. Both the small-scale and the large-scale setting are addressed and Newton-Kleinman-like methods are derived. The convergence of these procedures to the minimal solution is proved and several numerical results illustrate the computational efficiency of the proposed methods.

preprint2020arXiv

Space and Chaos-Expansion Galerkin POD Low-order Discretization of PDEs for Uncertainty Quantification

The quantification of multivariate uncertainties in partial differential equations can easily exceed any computing capacity unless proper measures are taken to reduce the complexity of the model. In this work, we propose a multidimensional Galerkin Proper Orthogonal Decomposition that optimally reduces each dimension of a tensorized product space. We provide the analytical framework and results that define and quantify the low-dimensional approximation. We illustrate its application for uncertainty modeling with Polynomial Chaos Expansions and show its efficiency in a numerical example.

preprint2020arXiv

Stabilization of Crystallization Models Governed by Hyperbolic Systems

This paper deals with mathematical models of continuous crystallization described by hyperbolic systems of partial differential equations coupled with ordinary and integro-differential equations. The considered systems admit nonzero steady-state solutions with constant inputs. To stabilize these solutions, we present an approach for constructing control Lyapunov functionals based on quadratic forms in weighted L2-spaces. It is shown that the proposed control design scheme guarantees exponential stability of the closed-loop system.

preprint2019arXiv

Analysis of switching strategies for the optimization of periodic chemical reactions with controlled flow-rate

An isoperimetric optimal control problem with non-convex cost is considered for a class of nonlinear control systems with periodic boundary conditions. This problem arises in chemical engineering as the maximization of the product of non-isothermal reactions by consuming a fixed amount of input reactants. It follows from the Pontryagin maximum principle that the optimal controls are piecewise constant in the considered case. We focus on a parametrization of optimal controls in terms of switching times in order to estimate the cost under different switching strategies. We exploit the Chen-Fliess functional expansion of solutions to the considered nonlinear system with bang-bang controls to satisfy the boundary conditions and evaluate the cost analytically for small periods. In contrast to the previous results in this area, the system under consideration is not control-affine, and the integrand of the cost depends on the state. This approach is applied to non-isothermal chemical reactions with simultaneous modulation of the input concentration and the volumetric flow-rate.

preprint2018arXiv

Hankel-Norm Approximation of Large-Scale Descriptor Systems

The Hankel-norm approximation is a model reduction method which provides the best approximation in the Hankel semi-norm. In this paper the computation of the optimal Hankel-norm approximation is generalized to the case of linear time-invariant continuous-time descriptor systems. An efficient algorithm is developed by refining the generalized balanced truncation square root method. For a wide practical usage, adaptations of the introduced algorithm towards stable computations and sparse systems are made as well as an approach for a projection-free algorithm. To show the approximation behavior of the introduced method, numerical examples are presented.

preprint2016arXiv

Exponential Stability and Stabilization of Extended Linearizations via Continuous Updates of Riccati Based Feedback

Many recent works on stabilization of nonlinear systems target the case of locally stabilizing an unstable steady state solutions against small perturbation. In this work we explicitly address the goal of driving a system into a nonattractive steady state starting from a well developed state for which the linearization based local approaches will not work. Considering extended linearizations or state-dependent coefficient representations of nonlinear systems, we develop sufficient conditions for stability of solution trajectories. We find that if the coefficient matrix is uniformly stable in a sufficiently large neighborhood of the current state, then the state will eventually decay. Based on these analytical results we propose an update scheme that is designed to maintain the stabilization property of Riccati based feedback constant during a certain period of the state evolution. We illustrate the general applicability of the resulting algorithm for setpoint stabilization of nonlinear autonomous systems and its numerical efficiency in two examples.

preprint2016arXiv

Space-time Galerkin POD with application in optimal control of semi-linear parabolic partial differential equations

In the context of Galerkin discretizations of a partial differential equation (PDE), the modes of the classical method of Proper Orthogonal Decomposition (POD) can be interpreted as the ansatz and trial functions of a low-dimensional Galerkin scheme. If one also considers a Galerkin method for the time integration, one can similarly define a POD reduction of the temporal component. This has been described earlier but not expanded upon -- probably because the reduced time discretization globalizes time which is computationally inefficient. However, in finite-time optimal control systems, time \textit{is} a global variable and there is no disadvantage from using a POD reduced Galerkin scheme in time. In this paper, we provide a newly developed generalized theory for space-time Galerkin POD, prove its optimality in the relevant function spaces, show its application for the optimal control of nonlinear PDEs, and, by means of a numerical example with Burgers' equation, discuss the competitiveness by comparing to standard approaches.