Order reduction of semilinear differential matrix and tensor equations

Kirsten, Gerhardus Petrus (2021) Order reduction of semilinear differential matrix and tensor equations, [Dissertation thesis], Alma Mater Studiorum Università di Bologna. Dottorato di ricerca in Matematica, 34 Ciclo. DOI 10.48676/unibo/amsdottorato/10010.
Documenti full-text disponibili:
[img] Documento PDF (English) - Richiede un lettore di PDF come Xpdf o Adobe Acrobat Reader
Disponibile con Licenza: Salvo eventuali più ampie autorizzazioni dell'autore, la tesi può essere liberamente consultata e può essere effettuato il salvataggio e la stampa di una copia per fini strettamente personali di studio, di ricerca e di insegnamento, con espresso divieto di qualunque utilizzo direttamente o indirettamente commerciale. Ogni altro diritto sul materiale è riservato.
Download (5MB)


In this thesis, we are interested in approximating, by model order reduction, the solution to large-scale matrix- or tensor-valued semilinear Ordinary Differential Equations (ODEs). Under specific hypotheses on the linear operators and the considered domain, these ODEs often stem from the space discretization on a tensor basis of semilinear Partial Differential Equations (PDEs) with a dimension greater than or equal to two. The bulk of this thesis is devoted to the case where the discrete system is a matrix equation. We consider separately the cases of general Lipschitz continuous nonlinear functions and the Differential Riccati Equation (DRE) with a quadratic nonlinear term. In both settings, we construct a pair of left-right approximation spaces that leads to a reduced semilinear matrix differential equation with the same structure as the original problem, which can be more rapidly integrated with matrix-oriented integrators. For the DRE, under certain assumptions on the data, we show that a reduction process onto rational Krylov subspaces obtains significant computational and memory savings as opposed to current approaches. In the more general setting, a challenging difference lies in selecting and constructing the two approximation bases to handle the nonlinear term effectively. In addition, the nonlinear term also needs to be approximated for efficiency. To this end, in the framework of the Proper Orthogonal Decomposition (POD) methodology and the Discrete Empirical Interpolation Method (DEIM), we derive a novel matrix-oriented reduction process leading to a practical, structure-aware low order approximation of the original problem. In the final part of the thesis, we consider the multidimensional setting. Here we extend the matrix-oriented POD-DEIM algorithm to the tensor setting and illustrate how we can apply it to systems of such equations. Moreover, we discuss how to integrate the reduced-order model and, in particular, how to solve the resulting tensor-valued linear systems.

Tipologia del documento
Tesi di dottorato
Kirsten, Gerhardus Petrus
Dottorato di ricerca
Settore disciplinare
Settore concorsuale
Parole chiave
Rational Krylov subspace method, Differential Riccati equation Proper orthogonal decomposition, Discrete empirical interpolation method, Semilinear matrix differential equations, Semilinear tensor differential equations, Exponential integrators
Data di discussione
2 Dicembre 2021

Altri metadati

Statistica sui download

Gestione del documento: Visualizza la tesi