Source author record

Tatjana Stykel

Tatjana Stykel appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

6works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

6 published item(s)

preprint2022arXiv

Analysis of a quasilinear coupled magneto-quasistatic model: solvability and regularity of solutions

We consider a~quasilinear model arising from dynamical magnetization. This model is described by a~magneto-quasistatic (MQS) approximation of Maxwell's equations. Assuming that the medium consists of a~conducting and a~non-conducting part, the derivative with respect to time is not fully entering, whence the system can be described by an abstract differential-algebraic equation. Furthermore, via magnetic induction, the system is coupled with an equation which contains the induced electrical currents along the associated voltages, which form the input of the system. The aim of this paper is to study well-posedness of the coupled MQS system and regularity of its solutions. Thereby, we rely on the classical theory of gradient systems on Hilbert spaces combined with the concept of $\mathcal{E}$-subgradients using in particular the magnetic energy. The coupled MQS system precisely fits into this general framework.

preprint2022arXiv

Energy-adaptive Riemannian optimization on the Stiefel manifold

This paper addresses the numerical solution of nonlinear eigenvector problems such as the Gross-Pitaevskii and Kohn-Sham equation arising in computational physics and chemistry. These problems characterize critical points of energy minimization problems on the infinite-dimensional Stiefel manifold. To efficiently compute minimizers, we propose a novel Riemannian gradient descent method induced by an energy-adaptive metric. Quantified convergence of the methods is established under suitable assumptions on the underlying problem. A non-monotone line search and the inexact evaluation of Riemannian gradients substantially improve the overall efficiency of the method. Numerical experiments illustrate the performance of the method and demonstrates its competitiveness with well-established schemes.

preprint2022arXiv

Passivity, Port-Hamiltonian Formulation and Solution Estimates for a Coupled Magneto-Quasistatic System

We study a~quasilinear coupled magneto-quasistatic model from a~systems theoretic perspective.} First, by taking the injected voltages as input and the associated currents as output, we prove that the magneto-quasistatic system is passive. Moreover, by defining suitable Dirac and resistive structures, we show that it admits a~representation as a~port-Hamiltonian system. Thereafter, we consider dependence on initial and input data. We show that the current and the magnetic vector potential can be estimated by means of the initial magnetic vector potential and the voltage. We also analyse the free dynamics of the system and study the asymptotic behavior of the solutions for $t\to\infty$.

preprint2021arXiv

Symplectic eigenvalue problem via trace minimization and Riemannian optimization

We address the problem of computing the smallest symplectic eigenvalues and the corresponding eigenvectors of symmetric positive-definite matrices in the sense of Williamson's theorem. It is formulated as minimizing a trace cost function over the symplectic Stiefel manifold. We first investigate various theoretical aspects of this optimization problem such as characterizing the sets of critical points, saddle points, and global minimizers as well as proving that non-global local minimizers do not exist. Based on our recent results on constructing Riemannian structures on the symplectic Stiefel manifold and the associated optimization algorithms, we then propose solving the symplectic eigenvalue problem in the framework of Riemannian optimization. Moreover, a connection of the sought solution with the eigenvalues of a special class of Hamiltonian matrices is discussed. Numerical examples are presented.

preprint2020arXiv

Balanced truncation for parametric linear systems using interpolation of Gramians: a comparison of algebraic and geometric approaches

When balanced truncation is used for model order reduction, one has to solve a pair of Lyapunov equations for two Gramians and uses them to construct a reduced-order model. Although advances in solving such equations have been made, it is still the most expensive step of this reduction method. Parametric model order reduction aims to determine reduced-order models for parameter-dependent systems. Popular techniques for parametric model order reduction rely on interpolation. Nevertheless, the interpolation of Gramians is rarely mentioned, most probably due to the fact that Gramians are symmetric positive semidefinite matrices, a property that should be preserved by the interpolation method. In this contribution, we propose and compare two approaches for Gramian interpolation. In the first approach, the interpolated Gramian is computed as a linear combination of the data Gramians with positive coefficients. Even though positive semidefiniteness is guaranteed in this method, the rank of the interpolated Gramian can be significantly larger than that of the data Gramians. The second approach aims to tackle this issue by performing the interpolation on the manifold of fixed-rank positive semidefinite matrices. The results of the interpolation step are then used to construct parametric reduced-order models, which are compared numerically on two benchmark problems.

preprint2013arXiv

Model Reduction of Descriptor Systems by Interpolatory Projection Methods

In this paper, we investigate interpolatory projection framework for model reduction of descriptor systems. With a simple numerical example, we first illustrate that employing subspace conditions from the standard state space settings to descriptor systems generically leads to unbounded H2 or H-infinity errors due to the mismatch of the polynomial parts of the full and reduced-order transfer functions. We then develop modified interpolatory subspace conditions based on the deflating subspaces that guarantee a bounded error. For the special cases of index-1 and index-2 descriptor systems, we also show how to avoid computing these deflating subspaces explicitly while still enforcing interpolation. The question of how to choose interpolation points optimally naturally arises as in the standard state space setting. We answer this question in the framework of the H2-norm by extending the Iterative Rational Krylov Algorithm (IRKA) to descriptor systems. Several numerical examples are used to illustrate the theoretical discussion.