Search

Solving optimal control problems for transport-dominated partial differential equations (PDEs) can become computationally expensive, especially when dealing with high-dimensional systems. To overcome this challenge, we focus on developing and deriving reduced-order models that can replace the full PDE system in solving the optimal control problem. Specifically, we explore the use of the shifted proper orthogonal decomposition (POD) as a reduced-order model, which is particularly effective for capturing high-fidelity, low-dimensional representations of transport-dominated phenomena. Furthermore, we propose two distinct frameworks for addressing these problems: one where the reduced-order model is constructed first, followed by optimization of the reduced system, and another where the original PDE system is optimized first, with the reduced-order model subsequently applied to the optimality system. We consider a 1D linear advection equation problem and compare the computational performance of the shifted POD method against the conventional methods like the standard POD when the reduced-order models are used as surrogates within a backtracking line search.

Port-Hamiltonian (pH) systems provide a powerful tool for modeling physical systems. Their energy-based perspective allows for the coupling of various subsystems through energy exchange. Another important class of systems, passive systems, are characterized by their inability to generate energy internally. In this paper, we explore first steps towards understanding the equivalence between passivity and the feasibility of port-Hamiltonian realizations in nonlinear systems. Based on our findings, we present a method to construct port-Hamiltonian representations of a passive system if the dynamics and the Hamiltonian are known.

We discuss structure-preserving time discretization for nonlinear port-Hamiltonian systems with state-dependent mass matrix. Such systems occur, for instance, in the context of structure-preserving nonlinear model order reduction for port-Hamiltonian systems and, in this context, structure-preserving time discretization is crucial for preserving some of the properties of the time-continuous reduced-order model. For this purpose, we introduce a new class of time discretization schemes which is based on so-called discrete gradient pairs and leads to an exact power balance on the time-discrete level. Moreover, for the special case of a pointwise symmetric and positive definite mass matrix, we present an explicit construction of a discrete gradient pair. Finally, we illustrate the theoretical findings by means of a numerical example, where the time-continuous system is a nonlinear reduced-order model for an advection-diffusion problem.

This work focuses on the conservation of quantities such as Hamiltonians, mass, and momentum when solution fields of partial differential equations are approximated with nonlinear parametrizations such as deep networks. The proposed approach builds on Neural Galerkin schemes that are based on the Dirac--Frenkel variational principle to train nonlinear parametrizations sequentially in time. We first show that only adding constraints that aim to conserve quantities in continuous time can be insufficient because the nonlinear dependence on the parameters implies that even quantities that are linear in the solution fields become nonlinear in the parameters and thus are challenging to discretize in time. Instead, we propose Neural Galerkin schemes that compute at each time step an explicit embedding onto the manifold of nonlinearly parametrized solution fields to guarantee conservation of quantities. The embeddings can be combined with standard explicit and implicit time integration schemes. Numerical experiments demonstrate that the proposed approach conserves quantities up to machine precision., Comment: 29 pages, 8 figures

We discuss structure-preserving model order reduction for port-Hamiltonian systems based on an approximation of the full-order state by a linear combination of ansatz functions which depend themselves on the state of the reduced-order model. In recent years, such nonlinear approximation ansatzes have gained more and more attention especially due to their effectiveness in the context of model reduction for transport-dominated systems which are challenging for classical linear model reduction techniques. We demonstrate that port-Hamiltonian reduced-order models can often be obtained by a residual minimization approach where a special weighted norm is used for the residual. Moreover, we discuss sufficient conditions for the resulting reduced-order models to be stable. Finally, the methodology is illustrated by means of two transport-dominated numerical test cases, where the ansatz functions are determined based on snapshot data of the full-order state.

Digital twins have become popular for their ability to monitor and optimize a process or a machine, ideally through its complete life cycle using simulations and sensor data. In this paper, we focus on the challenge of accurate and real-time simulations for digital twins in the context of electrical machines. To build such a digital twin involves not only computational models for the electromagnetic aspects, but also mechanical and thermal effects need to be taken into account. We address mathematical tools that can be employed to carry out the required simulations based on physical laws as well as surrogate or data-driven models. One of those tools is a model hierarchy of very fine to very coarse models as well as reduced order models for obtaining real-time simulations. The required software tools to carry out simulations in the digital twin are also discussed. The simulation models are implemented in a pipeline that allows for the automatic modeling of new machines and the automatic configuration of new digital twins. Finally, the overall implemented digital twin is tested and implemented in a physical demonstrator.

We present a new balancing-based structure-preserving model reduction technique for linear port-Hamiltonian descriptor systems. The proposed method relies on a modification of a set of two dual generalized algebraic Riccati equations that arise in the context of linear quadratic Gaussian balanced truncation for differential algebraic systems. We derive an a priori error bound with respect to a right coprime factorization of the underlying transfer function thereby allowing for an estimate with respect to the gap metric. We further theoretically and numerically analyze the influence of the Hamiltonian and a change thereof, respectively. With regard to this change of the Hamiltonian, we provide a novel procedure that is based on a recently introduced Kalman-Yakubovich-Popov inequality for descriptor systems. Numerical examples demonstrate how the quality of reduced-order models can significantly be improved by first computing an extremal solution to this inequality.

We study an optimization problem related to the approximation of given data by a linear combination of transformed modes. In the simplest case, the optimization problem reduces to a minimization problem well-studied in the context of proper orthogonal decomposition. Allowing transformed modes in the approximation renders this approach particularly useful to compress data with transported quantities, which are prevalent in many flow applications. We prove the existence of a solution to the infinite-dimensional optimization problem. Towards a numerical implementation, we compute the gradient of the cost functional and derive a suitable discretization in time and space. We demonstrate the theoretical findings with three challenging numerical examples., Comment: The only difference to the original version (v1) is a correction of the metadata

We propose a new hyper-reduction method for a recently introduced nonlinear model reduction framework based on dynamically transformed basis functions and especially well-suited for advection-dominated systems. Furthermore, we discuss applying this new method to a wildland fire model whose dynamics feature traveling combustion waves and local ignition and is thus challenging for classical model reduction schemes based on linear subspaces. The new hyper-reduction framework allows us to construct parameter-dependent reduced-order models (ROMs) with efficient offline/online decomposition. The numerical experiments demonstrate that the ROMs obtained by the novel method outperform those obtained by a classical approach using the proper orthogonal decomposition and the discrete empirical interpolation method in terms of run time and accuracy.

We study linear quadratic Gaussian (LQG) control design for linear port-Hamiltonian systems. To this end, we exploit the freedom in choosing the weighting matrices and propose a specific choice which leads to an LQG controller which is port-Hamiltonian and, thus, in particular stable and passive. Furthermore, we construct a reduced-order controller via balancing and subsequent truncation. This approach is closely related to classical LQG balanced truncation and shares a similar a priori error bound with respect to the gap metric. By exploiting the non-uniqueness of the Hamiltonian, we are able to determine an optimal pH representation of the full-order system in the sense that the error bound is minimized. In addition, we discuss consequences for pH-preserving balanced truncation model reduction which results in two different classical H-infinity-error bounds. Finally, we illustrate the theoretical findings by means of two numerical examples.

We propose a new model reduction framework for problems that exhibit transport phenomena. As in the moving finite element method (MFEM), our method employs time-dependent transformation operators and, especially, generalizes MFEM to arbitrary basis functions. The new framework is suitable to obtain a low-dimensional approximation with small errors even in situations where classical model order reduction techniques require much higher dimensions for a similar approximation quality. Analogously to the MFEM framework, the reduced model is designed to minimize the residual, which is also the basis for an a-posteriori error bound. Moreover, since the dependence of the transformation operators on the reduced state is nonlinear, the resulting reduced order model is obtained by projecting the original evolution equation onto a nonlinear manifold. Furthermore, for a special case, we show a connection between our approach and the method of freezing, which is also known as symmetry reduction. Besides the construction of the reduced order model, we also analyze the problem of finding optimal basis functions based on given data of the full order solution. Especially, we show that the corresponding minimization problem has a solution and reduces to the proper orthogonal decomposition of transformed data in a special case. Finally, we demonstrate the effectiveness of our method with several analytical and numerical examples.

We present a framework for constructing a structured realization of a linear time-invariant dynamical system solely from a discrete sampling of an input and output trajectory of the system. We estimate the transfer function of the original model at selected frequencies using a modification of the empirical transfer function estimation that was recently presented in [Peherstorfer, Gugercin, Willcox, SIAM J. Sci. Comput., 39(5):2152-2178, 2017]. Our realization interpolates the transfer function estimates and can be seen as a generalization of the Loewner framework to structured systems. We demonstrate the presented framework by means of a delay example.

We propose a new algorithm to compute a shifted proper orthogonal decomposition (sPOD) for systems dominated by multiple transport velocities. The sPOD is a recently proposed mode decomposition technique which overcomes the poor performance of classical methods like the proper orthogonal decomposition (POD) for transport-dominated phenomena. This is achieved by identifying the transport directions and velocities and by shifting the modes in space to track the transports. Our new algorithm carries out a residual minimization in which the main computational cost arises from solving a nonlinear optimization problem scaling with the snapshot dimension. We apply the algorithm to snapshot data from the simulation of a pulsed detonation combuster and observe that very few sPOD modes are sufficient to obtain a good approximation. For the same accuracy, the common POD needs ten times as many modes and, in contrast to the sPOD modes, the POD modes do not reflect the moving front profiles properly.

We introduce a novel model order reduction method for large-scale linear switched systems (LSS) where the coefficient matrices are affected by a low-rank switching. The key idea is to replace the LSS by a non-switched system with extended input and output vectors - called the envelope system - which is able to reproduce the dynamical behavior of the original LSS by applying a certain feedback law. The envelope system can be reduced using standard model order reduction schemes and then transformed back to an LSS. Furthermore, we present an upper bound for the output error of the reduced-order LSS and show how to preserve quadratic Lyapunov stability. The approach is tested by means of various numerical examples demonstrating the efficacy of the presented method.

We present a framework for constructing structured realizations of linear dynamical systems having transfer functions of the form $C(\sum_{k=1}^K h_k(s)A_k)^{-1}B$ where $h_1,h_2,\ldots,h_K$ are prescribed functions that specify the surmised structure of the model. Our construction is data-driven in the sense that an interpolant is derived entirely from measurements of a transfer function. Our approach extends the Loewner realization framework to more general system structure that includes second-order (and higher) systems as well as systems with internal delays. Numerical examples demonstrate the advantages of this approach.

Transport-dominated phenomena provide a challenge for common mode-based model reduction approaches. We present a model reduction method, which is suited for these kind of systems. It extends the proper orthogonal decomposition (POD) by introducing time-dependent shifts of the snapshot matrix. The approach, called shifted proper orthogonal decomposition (sPOD), features a determination of the {\it multiple} transport velocities and a separation of these. One- and two-dimensional test examples reveal the good performance of the sPOD for transport-dominated phenomena and its superiority in comparison to the POD., Comment: extended and clarified version