Search

The differential Riccati equation appears in different fields of applied mathematics like control and system theory. Recently Galerkin methods based on Krylov subspaces were developed for the autonomous differential Riccati equation. These methods overcome the prohibitively large storage requirements and computational costs of the numerical solution. In view of memory efficient approximation, we review and extend known solution formulas and identify invariant subspaces for a possibly low-dimensional solution representation. Based on these theoretical findings, we propose a Galerkin projection onto a space related to a low-rank approximation of the algebraic Riccati equation. For the numerical implementation, we provide an alternative interpretation of the modified \emph{Davison-Maki method} via the transformed flow of the differential Riccati equation, which enables us to rule out known stability issues of the method in combination with the proposed projection scheme. We present numerical experiments for large-scale autonomous differential Riccati equations and compare our approach with high-order splitting schemes., Comment: 33 pages, 38 figures

The differential Sylvester equation and its symmetric version, the differential Lyapunov equation, appear in different fields of applied mathematics like control theory, system theory, and model order reduction. The few available straight-forward numerical approaches if applied to large-scale systems come with prohibitively large storage requirements. This shortage motivates us to summarize and explore existing solution formulas for these equations. We develop a unifying approach based on the spectral theorem for normal operators like the Sylvester operator $\mathcal S(X)=AX+XB$ and derive a formula for its norm using an induced operator norm based on the spectrum of $A$ and $B$. In view of numerical approximations, we propose an algorithm that identifies a suitable Krylov subspace using Taylor series and use a projection to approximate the solution. Numerical results for large-scale differential Lyapunov equations are presented in the last sections., Comment: 30 pages, 51 figures

We provide spatial discretizations of nonlinear incompressible Navier-Stokes equations with inputs and outputs in the form of matrices ready to use in any numerical linear algebra package. We discuss the assembling of the system operators and the realization of boundary conditions and inputs and outputs. We describe the two benchmark problems - the driven cavity and the cylinder wake - and provide the corresponding data. The use of the data is illustrated by numerous example setups. The test cases are provided as plain PYTHON or OCTAVE/MATLAB script files for immediate replication.