### Page Content

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

## Absolventen-Seminar • Numerische Mathematik

Verantwortliche Dozenten: |
Prof. Dr. Christian Mehl [1], Prof. Dr. Volker Mehrmann
[2] |
---|---|

Koordination: | Ines
Ahrens [3] |

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 11.04. | 10:15 Uhr | MA
376 | no seminar | |

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

Benjamin Unger | Feedback regularization of DAEs via delays
[abstract] | |||

Do
25.04. | 10:15 Uhr | MA
376 | Volker Mehrmann [4] | Numerical analysis of finite element systems modeling
elastic stents [abstract] |

Do 02.05. | 10:15 Uhr | MA 376 | Christoph Zimmer | On
the solvability of port-Hamiltonian partial differential equation with
linear constraints [abstract] |

Paul Van Dooren | Role modeling using a low rank similarity matrix
[abstract] | |||

Do
09.05. | 10:15 Uhr | MA 376 | Michal
Wojtylak | The gap distance between linear
pencils [abstract] |

Paul Schwerdtner [5] | Structure
Preserving or Realization Independent H-infinity Approximation
[abstract] | |||

Do 16.05. | 10:15 Uhr | MA 376 | Ninoslav Truhar
[6] | On an Eigenvector-Dependent Nonlinear
Eigenvalue Problem from Perspective of Relative Perturbation Theory
[abstract] |

Dorothea
Hinsen | A modeling of a power network with the
telegraph equations [abstract] | |||

Do 23.05. | 10:15 Uhr | MA 376 | Jennifer Przybilla | Model
Reduction of Differential-Algebraic Systems by Parameter-Dependent
Balanced Truncation [abstract] |

Riccardo Morandin [7] | Energy-based hierarchical modeling of power networks
[abstract] | |||

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

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

Michelle Stahl | Analysis of optimization methods for nonlinear functions
to improve the geometry of injectors used for excavator ships
[abstract] | |||

Do
13.06. | 10:15 Uhr | MA 376 | Serhiy Yanchuk
[8] | On absolute stability of delay
differential equations with one discrete delay [abstract] |

Do
20.06. | 10:15 Uhr | MA 376 | Ute Kandler | Inexact methods for the solution of large scale Hermitian
eigenvalue problems [abstract] |

Philipp Krah [9] | Towards
mode-based model order reduction for flows with sharp fronts and
complex topolgies [abstract] | |||

Do 27.06. | 10:15 Uhr | MA 376 | ||

Do 04.07. | 10:15 Uhr | MA
376 | Ines Ahrens [10] | The Pantelides Algorithm for DAEs with delay
[abstract] |

Felix
Black | Computation of reduced order models for
transport phenomena via shifted proper orthogonal decomposition
[abstract] | |||

Do
11.07. | 10:15 Uhr | MA 376 | Daniel
Bankmann | Optimizing condition numbers in
robust stabilization [abstract] |

Marine Froidevaux | Structure-preservation in non-linear eigenvalue problems
arising from simulations of photonic crystals
[abstract] |

# Rückblick

- Absolventen Seminar WS 18/19 [11]
- Absolventen Seminar SS 18 [12]
- Absolventen Seminar WS 17/18 [13]
- Absolventen Seminar SS 17 [14]
- Absolventen Seminar WS 16/17 [15]
- Absolventen Seminar SS 16 [16]
- Absolventen Seminar WS 15/16 [17]
- Absolventen Seminar SS 15 [18]
- Absolventen Seminar WS 14/15 [19]
- Absolventen Seminar SS 14 [20]
- Absolventen Seminar WS 13/14 [21]
- Absolventen Seminar SS 13 [22]
- Absolventen Seminar WS 12/13 [23]
- Absolventen Seminar SS 12 [24]
- Absolventen Seminar WS 11/12 [25]

### Daniel Bankmann (TU Berlin)

Donnerstag, 11. Juli 2019

**Optimizing condition numbers
in robust stabilization**

Stabilization is a well-established task for linear control
systems, in which one tries to push all the eigenvalues of the system
matrix via a feedback to a region contained in the complex half plane
while optimizing certain condition numbers as the condition number of
the closed loop matrix.

The feedback gain can either be
computed by pole-placement techniques or by solving an optimal control
problem with the same resulting feedback.

It has been
shown, that for the numerically robust computation of these feedback
matrices, it is necessary to compute the solutions of the equivalent
optimal control problem via structure-exploiting computations of
deflating subspaces of an associated even matrix pencil. However, also
in that case, the computation of the feedback may be ill-conditioned
for a bad choice of weight parameters in the optimal control
problem.

In this talk we present an approach for detecting
when this computation is ill-conditioned. Also, we construct relations
that only keep well-conditioned parameters and decrease the dimension
of the parameter set. Finally, we show with some small study
example, that with this smaller parameter set, the results of
the overall optimization may improve in
accuracy.

### Marine Froidevaux (TU Berlin)

Donnerstag, 11. Juli 2019

**Structure-preservation in
non-linear eigenvalue problems arising from simulations of photonic
crystals**

Band-structure diagrams corresponding to photonic crystals can be numerically approximated by combining and discretizing the famous Maxwell equations of electromagnetism. In this talk we will focus on 2D photonic crystals where the electric permittivity is approximated by the rational Drude-Lorentz model. A finite element discretization of the system leads to parametric non-linear eigenvalue problems (NLEVPs) that are typically large, sparse and possess a symmetry in the spectrum.

Possible approaches to solve the NLEVPs are the use of contour integral methods, Newton-type methods or some linearization “by hand” followed by an iterative Krylov method. We will focus on the third approach and recall how to reduce the large-scale matrix equation with the help of the reduced basis method. We will see how to preserve the structure of the problem through the linearization and model reduction steps, as well as during the solution of the reduced system. We will also see that the extension of an existing structure-preserving Krylov method to large-scale complex even problems contains complications.

### Ines Ahrens (TU Berlin)

Donnerstag, 04. Juli 2019

**The Pantelides Algorithm for DAEs with
delay**

The strangeness index for DAEs is based on the derivative array,
which consists of the system itself plus its time derivatives. Index
reduction is performed by selecting certain important equations from
the derivative array. In a large-scale setting with high index, this
might become computationally infeasible. However, if it is known a
priori which equations of the original system need to be
differentiated, then the computational cost can be reduced. One way to
determine these equations is the Pantelides algorithm.

If
the DAE features in addition a delay term, taking derivatives might
not be sufficient, and instead, the derivative array must be shifted
in time, thus increasing the computational complexity even further. In
this talk I will explain how one can modify the Pantelides algorithm
such that it determines the equations which need to be shifted and/or
differentiated. This is joint work with Benjamin
Unger.

### Felix Black (TU Berlin)

Donnerstag, 04. Juli 2019

**Computation of reduced order
models for transport phenomena via shifted proper orthogonal
decomposition**

Many classical model order reduction methods are formulated in a
projection framework, building the reduced order model (ROM) by
projecting the full order model (FOM) onto a suitable subspace. If the
FOM exhibits advective transport, classical methods often fail to
produce low-dimensional models with a small approximation error. One
strategy to remedy this problem is the shifted proper orthogonal
decomposition (shifted POD), which extends the classical POD by
introducing additional transformation operators associated with the
modes. The transformation operators are parametrized by paths in a
suitable parameter space, thus allowing the transformed modes to cope
with the convection. In contrast to classical methods that project
onto a fixed subspace, the ROM of the shifted POD is thus obtained by
projecting onto a time- and/or state-dependent subspace that adapts
itself to the problem. On the one hand this approach is very flexible,
on the other hand it introduces additional complexity to the online
stage, since in addition to the time-dependent coefficients, also the
paths need to be computed from the ROM. In this talk, we present the
online stage of the shifted POD and discuss how to build a ROM from
which we may compute the coefficients and the paths. This is joint
work with Philipp Schulze and Benjamin Unger.

### Ute Kandler (Fraunhofer Institut)

Donnerstag, 20. Juni 2019

**Inexact methods for the solution of large
scale Hermitian eigenvalue problems**

This talk focuses on the solution of high dimensional Hermitian eigenproblems in situations where vector operations cannot be carried out exactly. To this end an inexact Arnoldi method with the aim to approximate extreme eigenvalues and eigenvectors is introduced. This method is particularly well-suited for large scale problems as it efficiently reduces the storage and computational requirements by constructing an orthonormal basis of the associated Krylov subspace. In particular, we investigate the influence of the inexactness of the vector operations on the convergence behavior of the inexact Arnoldi method and more general of inexact Krylov subspace methods. By using the concept of the angle of inclusion, a bound for the distance of the exact invariant subspace to an inexact Krylov subspace is analysed. To avoid the sensitivity of the bound to the gap between the desired and the remaining eigenvalues the concept of nested subspaces is used.

### Philipp Krah (TU Berlin)

Donnerstag, 20. Juni 2019

**Towards mode-based model order reduction
for flows with sharp fronts and complex topolgies**

I will present a model order reduction approach based on snapshots of 1D/2D flow fields, where the flow features sharp fronts. The approach parametrizes the flow field with the help of a levelset function which is zero at the location of the front and a nonlinear function which approximates the jump. I will explain how the freedom of choice of the levelset function can help to improve mode-based model order reduction. Therefore I will consider a 1D and 2D example of synthetic data and a 2D example of a simulated burning flame.

### Serhiy Yanchuk (TU Berlin)

Donnerstag, 13. Juni 2019

**On absolute stability of
delay differential equations with one discrete delay**

In this talk I discuss the properties of the eigenvalues of linear systems of DDEs of the form x'(t)=Ax(t) + B x(t-τ) with one discrete delay τ and constant matrices A and B. In particular, I show how the asymptotic spectrum for large delay determines the absolute stability, i.e. the stability for all positive delays.

### Pia Lutum (TU Berlin)

Donnerstag, 06. Juni 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 perturbations. In this talk we develop an algorithm to compute the real structured stability radius for large-scale linear time-invariant systems in the continuous-time case. Therefore, we first use model order reduction, especially rational interpolation. In the second step, we use a level-set method to compute the real structured stability radius for reduced systems. We improve this method by combining it with an algorithm which calculates the eigenvalues of an arising Hamiltonian matrix more precisely and faster. We illustrate the efficiency of the developed method using several numerical examples.

### Michelle Stahl (TU Berlin)

Donnerstag, 06. Juni 2019

**Analysis of optimization
methods for nonlinear functions to improve the geometry of injectors
used for excavator ships**

In this talk we will deal with numerical optimization methods for
nonlinear functions to improve the geometry of injectors applied in
dredging.

The reason is a working ship that is used as an
hydraulic excavator. It uses 36 equal nozzles for water injection but
it is known, that using nozzles of different sizes lead to better
results.

For improving the geometries as a multidimensional
problem I will introduce two optimization methods that work without
use of derivatives. The first one is a simplex algorithm, called
“Nelder-Mead method” or “Downhill- Simplex method”. Then I
will show the idea and the application of the evolution strategy by
Ingo Rechenberg, a genetic algorithm that mimics the biological
evolution.

### Jennifer Przybilla (TU Berlin)

Donnerstag, 23. Mai 2019

**Model Reduction of
Differential-Algebraic Systems by Parameter-Dependent Balanced
Truncation**

In this talk we will deal with the application of balanced truncation for parameter-dependent differential-algebraic systems. For this, we have to solve parameter-dependent Lyapunov equations, which we do with the help of the reduced-basis method. In order to use this method, we first have to deal with the algebraic parts of the system and make the system strictly dissipative in order to apply error estimators in the reduced basis method.

### Riccardo Morandin (TU Berlin)

Donnerstag, 23. Mai 2019

**Energy-based hierarchical modeling of
power networks**

Recent developments in the energy market require new mathematical models and suitable algorithms for the efficient usage of the existing networks. Starting from these real-world problems, it is necessary to concentrate on methods for large-scale energy networks and, in particular, address optimization and stability analysis, model predictive control, model-order reduction, uncertainty quantification, and related topics on energy networks. The abstract setting allows for consideration of applications arising from both, gas and power networks.

Our framework of choice is the one of port-Hamitlonian
descriptor systems, or pHDAEs. These are energy-based equations with a
special structure, that guarantees several beneficial properties, e.g.
inherent stability and passivity, physical interpretation of
variables, simple interconnection, structure-preserving model
reduction, robust numerical integration and simplification of feedback
stabilization.

In this talk we present a model hierarchy for some components
of a power grid, specifically synchronous generators and transmission
lines, where all equations presented are pHDAEs. Furthermore, we show
how these components can be easily interconnected to form a larger
network, while preserving the port-Hamiltonian structure. To do so, we
apply Kirchhoff's laws and exploit the common structure of the
electrical components, in a way that can be easily extended to other
devices and more complex models.

### Dorothea Hinsen (TU Berlin)

Donnerstag, 16. Mai 2019

**A modeling of a power network with the
telegraph equations **

In recent years energy transition and the increasing electricity
demand have led to growing interest in modeling power networks, which
withstand unexpected occurrences as voltage or transient
instabilities.

One way to approach modeling power networks
is with port-Hamiltonian systems.

The power networks we are
dealing with consist of generators, loads and transmission lines.

In this talk we discuss an approach to a power network model,
where the generators and the loads are described with port-Hamiltonian
equations.

However, the transmission lines are depicted
with the telegraph equations. This leads us to a PDAE model, which we
will discuss.

### Ninoslav Truhar (Josip Juraj Strossmayer University of Osijek)

Donnerstag, 16. Mai 2019

**On an Eigenvector-Dependent Nonlinear
Eigenvalue Problem from Perspective of Relative Perturbation Theory
**

The talk contains two part. I the first part we consider the
eigenvector-dependent nonlinear eigenvalue problem (NEPv) H(V) V = V
Λ, where H(V) ∈ C^{n,n} is an Hermitian matrix-valued
function of V ∈ C^{n,k} with orthonormal columns, i.e.,
V^{H} V = I_{k}, k < n (usually k ≪ n). We
present the conditions on existence and uniqueness for the solvability
of NEPv using the well known results of the relative perturbation
theory. Our results are motiveted by the results on NEPv presented in
Y. Cai, L.-H. Zhang, Z. Bai, and R.-C. Li, *On an
eigenvector-dependent nonlinear eigenvalue problem*, SIAM J.
Matrix Anal. Appl. 2018, where among the other results one can find
conditions for existence and uniqueness for the solvability of an
NEPv. These results are based on well-known standard perturbation
theory for Hermitian matrices. The differences between so called
standard perturbation theory approach, and our new (relative
perturbation theory) approach have been illustrated in several
numerical examples.

In the second part we present an upper
and a lower bound for the the Frobenius norm of the matrix sin(Θ), of
the sines of canonical angles between unperturbed and perturbed
eigenspaces of a regular generalized Hermitian eigenvalue problem A x
= λ B x where A and B are Hermitian n × n matrices, under a
feasible non Hermitian perturbation. As one application of the
obtained bounds we present the corresponding upper and the lower
bounds for eigenspaces of a matrix pair (A,B) obtained by a
linearization of regular quadratic eigenvalue problem ( λ^{2}
M + λ D + K ) u = 0, where M is positive definite and D and K are
semidefinite. We also apply obtained upper and lower bounds to the
important problem which considers the influence of adding a damping on
mechanical systems. The new results show that for certain additional
damping the upper bound can be too pessimistic, but the lower bound
can reflect a behaviour of considered eigenspaces
properly.

### Michal Wojtylak (Jagiellonian University, Krakow)

Donnerstag, 09. Mai 2019

**The gap distance between linear
pencils**

We introduce a new distance between linear pencils, which is based
on their graphs ker[A,-E]. Several basic properties of this distance
will be shown. Next, we will formulate the problem of distance to
singularity and compare it with the original one given by Byers, He
and Mehrmann in 1998. The talk is based on a joint work: Berger,
Gernandt, Trunk, Winkler, MW, LAA 2019.

### Paul Schwerdtner (TU Berlin)

Donnerstag, 09. Mai 2019

**Structure Preserving or Realization
Independent H-infinity Approximation**

In this talk we revisit linorm_subsp, an algorithm that is used to
compute the H-infinity norm of transfer functions of large-scale
dynamical systems and show how it can be made more efficient for
irrational transfer functions using an inner loop outer loop
strategy.

Then, we present a greedy interpolation approach to
address the H-infinity model order reduction problem that is using the
evaluation of the H-infinity norm computation provided by
linorm_subsp.

Starting from an initial reduced order model, we
compute the H-infinity norm of the error transfer function and then
place an new interpolation point where this H-infinity norm is
attained. In this way, we construct a sequence of reduced order models
aiming at minimizing the H-infinity norm of the error transfer
function.

Using interpolation for this allows to either preserve
the given model structure or construct surrogate models with a
structure defined by the user to approximate the given system,
respectively.

### Christoph Zimmer (TU Berlin)

Donnerstag, 02. Mai 2019

**On the solvability of
port-Hamiltonian partial differential equation with linear
constraints **

Through its natural approach of energy as underlying quantity and its many nice properties, port-Hamiltonan (pH) systems rose as a popular modeling tool for dynamical system like electrical, mechanical, and hydrodynamical ones. Therefore, it is not surprising that the concept was extended to non-autonomous systems, to descriptor systems and to infinite dimensional systems. In this talk, we investigate non-autonomous semi-explicit pH descriptor systems. Mainly focusing on the existence of solution, we will start with finite dimensional systems. Afterwards, we will transfer the used tricks to the infinite dimensional case and prove the existence of solutions as well as their smoothness.

### Paul Van Dooren (Catholic University of Louvain)

Donnerstag, 02. Mai 2019

**Role modeling using a low rank similarity
matrix**

Community detection is a popular approach used analyze large networks and obtain relevant statistical properties by finding community structures in networks. However, community detection algorithms cannot be used to find non-community structures in networks, such as cyclic graph structures which appear in food web networks. The role extraction problem represents large networks by a smaller, general graph structure, called role graphs. In this presentation, we analyze the neighborhood pattern similarity measure used to solve the role extraction problem. We show theoretically that under certain assumptions the role graphs can be recovered from a low-rank factorization of the similarity matrix due to the relationship between the dominant eigenvalues of the similarity matrix and the number of roles. We also show how it applies to the detection of overlapping communities and bipartite communities.

### Volker Mehrmann (TU Berlin)

Donnerstag, 25. April 2019

** Numerical analysis of finite element
systems modeling elastic stents**

A new model description for the numerical simulation of elastic stents is proposed. Based on the new formulation an inf-sup inequality for the finite element discretization is proved and the proof of the inf-sup inequality for the continuous problem is simplified. The new formulation also leads to faster simulation times despite an increased number of variables. The techniques also simplify the analysis and numerical solution of the evolution problem describing the movement of the stent under external forces. The results are illustrated via numerical examples, see [1].

[1] L. Grubišić, M. Ljulj, V. Mehrmann, and J. Tambača, Modeling and discretization methods for the numerical simulation of elastic stents, arxiv.org/1812.10096 [26], Preprint 01-2019, Institute of Mathematics, TU Berlin, submitted for publication, 2019.

### Benjamin Unger (TU Berlin)

Donnerstag, 18. April 2019

**Feedback regularization of
DAEs via delays**

We study linear time-invariant delay differential-algebraic equations (DDAEs). Such equations arise naturally, if a feedback controller is applied to a descriptor system, since the controller requires some time to measure the state and to compute the feedback resulting in the time-delay. We present an existence and uniqueness result for DDAEs within the space of piecewise smooth distributions and an algorithm to determine whether a DDAE is delay-regular. As a consequence, we show that a DAE can be regularized by a feedback if and only if it can be regularized by a delayed feedback. This is joint work with Stephan Trenn (University of Groningen).

diff/fg_numerische_mathematik/v_menue/mitarbeiter/chris

tian_mehl/parameter/en/minhilfe/

diff/fg_numerische_mathematik/v_menue/mitarbeiter/mehrm

ann/volker_mehrmann_home/parameter/en/minhilfe/

diff/fg_numerische_mathematik/v_menue/mitarbeiter/ines_

ahrens/home/parameter/en/minhilfe/

diff/fg_numerische_mathematik/v_menue/mitarbeiter/mehrm

ann/volker_mehrmann_home/parameter/en/minhilfe/

diff/fg_numerische_mathematik/v_menue/mitarbeiter/paul_

schwerdtner/home/parameter/en/minhilfe/

vnici-i-suradnici/120

diff/fg_numerische_mathematik/v_menue/mitarbeiter/ricca

rdo_morandin/home/parameter/en/minhilfe/

diff/gp_angewandte_dynamische_systeme/v_menue/yanchuk/v

_menue/homepage_of_serhiy_yanchuk/parameter/en/minhilfe

/

mechanik_und_technische_akustik_ista/numerische_fluiddy

namik/menue/about_us/people/philipp_krah/parameter/en/m

inhilfe/

mdiff/fg_numerische_mathematik/v_menue/mitarbeiter/ines

_ahrens/home/parameter/en/minhilfe/

mdiff/fg_numerische_mathematik/v_menue/vortraege/absolv

enten_seminar/absolventen_seminar_ws_1819/parameter/en/

minhilfe/

mdiff/fg_numerische_mathematik/v_menue/vortraege/absolv

enten_seminar/absolventen_seminar_ss_18/parameter/en/mi

nhilfe/

mdiff/fg_numerische_mathematik/v_menue/vortraege/absolv

enten_seminar/absolventen_seminar_ws_1718/parameter/en/

minhilfe/

mdiff/fg_numerische_mathematik/v_menue/vortraege/absolv

enten_seminar_main/absolventen_seminar_ss_17/parameter/

en/minhilfe/

mdiff/fg_numerische_mathematik/v_menue/vortraege/absolv

enten_seminar_main/absolventen_seminar_ws_1617/paramete

r/en/minhilfe/

mdiff/fg_numerische_mathematik/v_menue/vortraege/absolv

enten_seminar_main/absolventen_seminar_ss_16/parameter/

en/minhilfe/

mdiff/fg_numerische_mathematik/v_menue/vortraege/absolv

enten_seminar/absolventen_seminar_ws_1516/parameter/en/

minhilfe/

mdiff/fg_numerische_mathematik/v_menue/vortraege/absolv

enten_seminar_main/absolventen_seminar_ss_15/parameter/

en/minhilfe/

mdiff/fg_numerische_mathematik/v_menue/vortraege/absolv

enten_seminar_main/absolventen_seminar_ws_1415/paramete

r/en/minhilfe/

mdiff/fg_numerische_mathematik/v_menue/vortraege/absolv

enten_seminar_main/absolventen_seminar_ws_1314/paramete

r/en/minhilfe/

mdiff/fg_numerische_mathematik/v_menue/vortraege/absolv

enten_seminar_main/absolventen_seminar_ws_1314/paramete

r/en/minhilfe/

mdiff/fg_numerische_mathematik/v_menue/vortraege/absolv

enten_seminar_main/absolventen_seminar_ss_13/parameter/

en/minhilfe/

mdiff/fg_numerische_mathematik/v_menue/vortraege/absolv

enten_seminar_main/absolventen_seminar_ws_1213/paramete

r/en/minhilfe/

mdiff/fg_numerische_mathematik/v_menue/vortraege/absolv

enten_seminar_main/absolventen_seminar_ss_12/parameter/

en/minhilfe/

mdiff/fg_numerische_mathematik/v_menue/vortraege/absolv

enten_seminar_main/absolventen_seminar_ws_1112/paramete

r/en/minhilfe/