### Page Content

### to Navigation

## There is no English translation for this web page.

## Absolventen-Seminar • Numerische Mathematik

Verantwortliche Dozenten: | Prof. Dr. Christian Mehl, Prof. Dr. Volker Mehrmann |
---|---|

Koordination: | Benjamin Unger, Dr. Matthias Voigt |

Termine: | Do 10:00-12:00 in MA 376 |

Inhalt: | Vorträge von Bachelor- und Masterstudenten, Doktoranden, Postdocs und manchmal auch Gästen zu aktuellen Forschungsthemen |

Datum | Zeit | Raum | Vortragende(r) | Titel |
---|---|---|---|---|

Do 18.10. | 10:15 Uhr | MA 376 | Vorbesprechung | |

Do 25.10. | 10:15 Uhr | MA 376 | Alexander Grimm | Empirical least-squares fitting of parametrized dynamical systems [abstract] |

Do 01.11. | 10:15 Uhr | MA 376 | no seminar | |

Do 08.11. | 10:15 Uhr | MA 376 | Andrii Dmytryshyn | Geometry of matrix polynomial spaces [abstract] |

Christian Mehl | Generic low-rank perturbations of matrix pencils with symmetry structures [abstract] | |||

Do 15.11. | 10:15 Uhr | MA 376 | no seminar | |

Do 22.11. | 10:15 Uhr | MA 376 | Riccardo Morandin | A new formulation for port-Hamiltonian differential-algebraic equations [abstract] |

Olga Markova | Length realizability question for pairs of quasi-commuting matrices [abstract] | |||

Do 29.11. | 10:15 Uhr | MA 376 | no seminar | |

Do 06.12. | 10:15 Uhr | MA 376 | Arbi Moses Badlyan | - canceled - |

Volker Mehrmann | Computation of Stability Radii for Large-Scale Dissipative Hamiltonian Systems [abstract] | |||

Do 13.12. | 10:15 Uhr | MA 376 | Daniel Bankmann | Adjoint sensitivity equations in optimal control of differential-algebraic equations [abstract] |

Paul Schwerdtner | Robust Control of Delay Systems [abstract] | |||

Do 20.12. | 10:15 Uhr | MA 376 | Arbi Moses Badlyan | Generalized Port-Hamiltonian Systems - Turbulent Pipe Flow [abstract] |

Do 10.01. | 10:15 Uhr | MA 376 | Pia Lutum | Numerical Computation of the Real Structured Stability Radius [abstract] |

Philipp Schulze | - canceled - | |||

Do 17.01. | 10:15 Uhr | MA 376 | Manuel Radons | Semi-automatically optimized calibration of internal combustion engines [abstract] |

Dorothea Hinsen | - canceled - | |||

Do 24.01. | 10:15 Uhr | MA 376 | Rico Berner | Multi-cluster structures in networks of adaptively coupled oscillators [abstract] |

Vander Freitas (National Institute for Space Research INPE, Brazil) | Symmetric circular formations with unitary speed particles [abstract] | |||

Do 31.01. | 10:15 Uhr | MA 376 | Marine Froidevaux | |

Tim Mitchell | ||||

Do 07.02. | 10:15 Uhr | MA 376 | Matthias Voigt | |

Christoph Zimmer | ||||

Do 14.02. | 10:15 Uhr | MA 376 | Benjamin Unger | |

Ines Ahrens |

# Rückblick

- Absolventen Seminar SS 18
- Absolventen Seminar WS 17/18
- Absolventen Seminar SS 17
- Absolventen Seminar WS 16/17
- Absolventen Seminar SS 16
- Absolventen Seminar WS 15/16
- Absolventen Seminar SS 15
- Absolventen Seminar WS 14/15
- Absolventen Seminar SS 14
- Absolventen Seminar WS 13/14
- Absolventen Seminar SS 13
- Absolventen Seminar WS 12/13
- Absolventen Seminar SS 12
- Absolventen Seminar WS 11/12

### Alexander Grimm (Universität Suttgart)

Donnerstag, 25. Oktober 2018

**Empirical least-squares fitting of parametrized dynamical systems**

Given a set of response observations for a parametrized dynamical system, we seek a parametrized dynamical model that will yield uniformly small response error over a range of parameter values yet has low order. Frequently, access to internal system dynamics or equivalently, to realizations of the original system is either not possible or not practical; only response observations over a range of parameter settings might be known. Respecting these typical operational constraints, we propose a two phase approach that first encodes the response data into a high fidelity intermediate model of modest order, followed then by a compression stage that serves to eliminate redundancy in the intermediate model. For the first phase, we extend non-parametric least-squares fitting approaches so as to accommodate parameterized systems. This results in a (discrete) least-squares problem formulated with respect to both frequency and parameter that identifies "local" system response features. The second phase uses an H2-optimal model reduction strategy accommodating the specialized parametric structure of the intermediate model obtained in the first phase. The final compressed model inherits the parametric dependence of the intermediate model and maintains the high fidelity of the intermediate model, while generally having dramatically smaller system order. We provide a variety of numerical examples demonstrating our approach.

### Andrii Dmytryshyn (Umeå University, Sweden)

Donnerstag, 08. November 2018

**Geometry of matrix polynomial spaces**

We study how small perturbations of matrix polynomials may change their elementary divisors and minimal indices by constructing the closure hierarchy graphs (stratifications) of orbits and bundles of matrix polynomial Fiedler linearizations. We show that the stratification graphs do not depend on the choice of Fiedler linearization which means that all the spaces of the matrix polynomial Fiedler linearizations have the same geometry (topology). We also develop the theory for structure preserving stratification of skew-symmetric matrix polynomials. The results are illustrated by examples using the software tool Stratigraphy.

This is a joint work with Stefan Johansson (Umeå University), Bo Kågström (Umeå University), and Paul Van Dooren (Université catholique de Louvain).

### Christian Mehl (TU Berlin)

Donnerstag, 08. November 2018

**Generic low-rank perturbations of matrix pencils with symmetry structures**

In this talk, we discuss the effect that generic structure-preserving low-rank perturbations have on the Jordan structure of a given matrix pencil with symmetry structures. The main result relies on a specific parameterizations of the sets of structured matrix pencils with a given rank.

### Riccardo Morandin (TU Berlin)

Donnerstag, 22. November 2018

**A new formulation for port-Hamiltonian differential-algebraic equations**

Coordinate-based port-Hamiltonian system formulations in the literature are usually either limited to represent ODEs, or to be associated with square systems and quadratic Hamiltonian functions. In this talk, I will present a new formulation for port-Hamiltonian differential-algebraic equations (pHDAEs), that does not have the previous restrictions.

Two equivalent formulations are also provided, that are useful for simplifying the notation, proving some additional results and reaching a better understanding of the new definition. It is proven that this class of pHDAEs is closed under a large class of variable transformations and under linear system interconnection.

A geometric interpretation (Dirac structure) is also associated to the new formulation, and will be useful to extend it to a more abstract context. Collocation-based time integration is analyzed, and with additional assumptions a discrete port-Hamiltonian structure is achieved.

Furthermore, the definition is extended to be applied to infinite-dimensional systems (e.g. PDEs), while preserving most of the previous results.

Finally, a few examples of systems that can be represented with this new formulation are shown.

### Olga Markova (Lomonosov State University, Moscow, Russia)

Donnerstag, 22. November 2018

**Length realizability question for pairs of quasi-commuting matrices**

By the length of a finite system of generators for a finite-dimensional algebra over an arbitrary field we mean the least positive integer k such that the products of length not exceeding k span this algebra (as a vector space). The maximum length for the systems of generators of an algebra is referred to as the length of the algebra. The length evaluation can be a difficult problem, for example, the length of the full matrix algebra is unknown (Paz’s Problem, 1984).

In this talk we discuss the length evaluation problem for quasi-commuting pairs of matrices (we say that matrices A and B quasi-commute if the products AB and BA are linearly dependent). We provide sharp upper and lower bounds for the length of such pairs and show how the interval between these extremal values is divided into intervals of realizable values for the length and ``gaps'', i.e. non-realizable values.

The talk is based on a joint work with Alexander Guterman and Volker Mehrmann.

### Arbi Moses Badlyan (TU Berlin)

Donnerstag, 06. Dezember 2018

**Generalized Port-Hamiltonian Systems – Turbulent Pipe Flow **

In this talk I present the notion of a generalized port-Hamiltonian system based on modified, explicitly state dependent Stokes-Dirac structures.

I will demonstrate the mathematical framework by taking the example of a turbulent pipe flow. The mathematical model of this pipe flow is given as system of coupled partial-differential equations complemented by closure relations and contains a non-local heat-conduction term. The presented results are part of a joint work with Prof. Volker Mehrmann (TU Berlin).

### Volker Mehrmann (TU Berlin)

Donnerstag, 06. Dezember 2018

**Computation of Stability Radii for Large-Scale Dissipative Hamiltonian Systems**

A linear time-invariant dissipative Hamiltonian (DH) system is always Lyapunov stable and under weak further conditions even asymptotically stable. In various applications there is uncertainty on the system matrices, and it is desirable to know whether the system remains asymptotically stable uniformly against all possible uncertainties within a given perturbation set. Such robust stability considerations motivate the concept of stability radius for DH systems, i.e., what is the maximal perturbation permissible to the coefficients, while preserving the asymptotic stability. We consider two stability radii, the unstructured one where the coefficients are subject to unstructured perturbation, and the structured one where the perturbations preserve the DH structure.

We propose new algorithms to compute these stability radii for large scale problems by tailoring subspace frameworks that are interpolatory and guaranteed to converge at a super-linear rate in theory. At every iteration, they first solve a reduced problem and then expand the subspaces in order to attain certain Hermite interpolation properties between the full and reduced problems. The reduced problems are solved by means of the adaptations of existing level-set algorithms for H-infinity-norm computation in the unstructured case, while, for the structured radii, we benefit from algorithms that approximate the objective eigenvalue function with a piece-wise quadratic global underestimator. The performance of the new approaches is illustrated with several examples including a system that arises from a finite-element modeling of an industrial disk brake.

Joint work with N. Aliyev, V. Mehrmann, and E. Mengi

### Daniel Bankmann (TU Berlin)

Donnerstag, 13. Dezember 2018

**Adjoint sensitivity equations in optimal control of differential-algebraic equations **

Optimal control problems for differential-algebraic equations appear in a variety of applications coming from electrical or mechanical engineering.

Sometimes, the system's descriptions are also depending on parameters and one is not only interested in computing the solution of such a problem. In addition, we are interested in computing how the solution changes for small changes in the parameters, i. e. the sensitivities of the optimal solution with respect to the parameters.

Solutions of the optimal control problem can be characterized by so-called necessary conditions. These constitute a boundary value problem for differential-algebraic equations.

In the current literature, only sensitivity analysis for initial value problems of differential-algebraic equations or boundary value problems for ordinary differential equations has been carried out.

We close this gap, by formulating the correct boundary conditions based on the standard adjoint equations for differential-algebraic equations of index 1. This includes jump conditions on certain variables which pose new challenges on solvability and existence of solutions.

In this talk, we present an adjoint sensitivity boundary value problem for index 1 differential-algebraic equations with boundary values. Based on a flow formulation, we analyze properties like existence and uniqueness and apply these results to the necessary conditions coming from the optimal control problem.

### Paul Schwerdtner (TU Berlin)

Donnerstag, 13. Dezember 2018

**Robust Control of Delay Systems**

We review the shortcomings of an algorithm that computes the H-infinity norm of possibly irrational and large scale transfer functions using a subspace projection based approach. After that, we show how these can be overcome by employing rational interpolation in the Loewner matrix framework and compare our new approach to state-of-the-art methods in various numerical examples.

After that, we utilize the new computation method for H-infinity norms in an optimization loop to compute robust fixed order controllers for time delay systems. Specifically, we alter the controller parameters of a given fixed order controller to optimize the H-infinity norm of the resulting closed loop transfer function.

Furthermore, we show how the realization independence of the H-infinity norm computation method can be exploited to facilitate imposing performance requirements on the closed loop transfer function.

Finally, we describe future work that is concerned with an efficient check for asymptotic stability of delay systems, since that is currently the bottleneck in our optimization based approach.

### Arbi Moses Badlyan (TU Berlin)

Donnerstag, 20. Dezember 2018

**Generalized Port-Hamiltonian Systems - Turbulent Pipe Flow**

In this talk I will present the notion of a generalized port-Hamiltonian system based on explicitly state dependent Stokes-Dirac structures.

I will demonstrate the mathematical framework by taking the example of a turbulent pipe flow. The mathematical model of this pipe flow is given as system of coupled partial-differential equations complemented by closure relations and contains a non-local transport term. The presented results are part of a joint work with Prof. Christopher Beattie (Virginia Tech) and Prof. Volker Mehrmann (TU Berlin).

### Pia Lutum (TU Berlin)

Donnerstag, 10. Januar 2019

**Numerical Computation of the Real Structured Stability Radius**

The real structured stability radius is a measurement of the robustness of a stable matrix under certain real perturbations. In this talk we briefly discuss an iterative algorithm to compute the stability radius. Based on this algorithm, we further discuss an iterative method to compute the real structured stability radius for large-scale matrices by projections.

### Manuel Radons (TU Berlin)

Donnerstag, 17. Januar 2019

**Semi-automatically optimized calibration of internal combustion engines**

Modern combustion engines incorporate a number of actuators and sensors that can be used to control and optimize the performance and emissions. We describe a semi-automatic method to simultaneously measure and calibrate the actuator settings and the resulting behavior of the engine. The method includes an adaptive process for refining the measurements, a data cleaning step, and an optimization procedure. The optimization works in a discretized space and incorporates the conditions to describe the dependence between the actuators and the engine behavior as well as emission bounds. We demonstrate our method on practical examples.

### Rico Berner (TU Berlin)

Donnerstag, 24. Januar 2019

**Multi-cluster structures in networks of adaptively coupled oscillators**

Dynamical systems on networks with adaptive couplings appear naturally in real-world systems such as power grid networks, social networks as well as neuronal networks. We investigate collective behaviour in a paradigmatic network of adaptively coupled phase oscillators. The coupling topology of the network changes slowly depending on the dynamics of the oscillators. We show that such a system gives rise to numerous complex dynamics, including relative equilibria and hierarchical multi-cluster states. An analytic treatment for equilibria and multi-cluster solutions as well as the existence of continuous families of these states is presented and parameter regimes of high multi-stability are found. In addition, we give an interpretation for equilibria as functional units which are building blocks in multi-cluster structures. Our results contribute to the understanding of mechanisms for pattern formation in adaptive networks, such as the emergence of multi-layer structure in neural systems.

### Vander Freitas (National Institute for Space Research INPE, Brazil)

Donnerstag, 24. Januar 2019

**Symmetric circular formations with unitary speed particles**

Several living beings such as birds, fish, flies and bees perform collective motion when foraging, migrating or escaping from predators. They are able to flock to a certain direction, rotate around a common center, or swarm within an area. Biologists try to understand the local rules that lead single individuals to macroscopic group behaviors. Inspired on that, it is possible to develop models for mobile vehicles, whose applications include data collection, surveillance, monitoring, etc. Motivated by such problems, we study a model of particles with phase-coupled oscillators dynamics, focusing on the special case in which they rotate around a common center and are able to group into clusters. Those formations arise from an optimization procedure of certain potentials. We then investigate the parameter space, the impact of adding and removing particles, and how to exchange from one configuration to another.