Model reduction for kinetic equations: moment approximations and hierarchical approximate proper orthogonal decomposition
Momentenmodelle sind eine Klasse von spezialisierten approximativen Modellen für kinetische Transportgleichungen. Diese Modelle wandeln die kinetische Gleichung in ein geschwindigkeitsunabhängiges Gleichungssystem für die gewichteten Mittelwerte (Momente) der Lösung um. Je nachdem, welche Gewichtsfu...
| Verfasser: | |
|---|---|
| Weitere Beteiligte: | |
| FB/Einrichtung: | FB 10: Mathematik und Informatik |
| Dokumenttypen: | Dissertation/Habilitation |
| Medientypen: | Text |
| Erscheinungsdatum: | 2021 |
| Publikation in MIAMI: | 18.01.2022 |
| Datum der letzten Änderung: | 18.01.2022 |
| Angaben zur Ausgabe: | [Electronic ed.] |
| Schlagwörter: | Momentenmodelle; Modellreduktion; HAPOD; lineare kinetische Gleichungen; Realisierbarkeit; Entropieminimierung; Finite-Volumen-Verfahren moment models; model reduction; HAPOD; linear kinetic equations; realizability; minimum entropy closure; finite volume methods |
| Fachgebiet (DDC): | 510: Mathematik |
| Lizenz: | InC 1.0 |
| Sprache: | English |
| Hochschulschriftenvermerk: | Zugl.: Münster (Westfalen), Univ., Diss., 2021 |
| Förderung: | Förderer: Deutsche Forschungsgemeinschaft / Projektnummer: 390685587 Förderer: Deutsche Forschungsgemeinschaft / Projektnummer: 194347757 |
| Format: | PDF-Dokument |
| URN: | urn:nbn:de:hbz:6-85019477402 |
| Permalink: | https://nbn-resolving.de/urn:nbn:de:hbz:6-85019477402 |
| Verwandte Dokumente: | |
| Onlinezugriff: | diss_leibner_2021.pdf |
| LEADER | 12189cam a2200421uu 4500 | ||
|---|---|---|---|
| 001 | 9952d183-16d5-4d0a-a60e-188c6d8cc844 | ||
| 003 | miami | ||
| 005 | 20220118 | ||
| 007 | c||||||||||||a| | ||
| 008 | 220118e20220118||||||||||#s||||||||eng|||||| | ||
| 024 | 7 | |a urn:nbn:de:hbz:6-85019477402 |2 urn | |
| 041 | |a eng | ||
| 082 | 0 | |a 510 Mathematik |2 23 | |
| 100 | 1 | |a Leibner, Tobias |0 http://d-nb.info/gnd/1130666093 |0 http://viaf.org/viaf/2061149368877885980006 |4 aut | |
| 110 | 2 | |a Universitäts- und Landesbibliothek Münster |0 http://d-nb.info/gnd/5091030-9 |4 own | |
| 245 | 1 | 0 | |a Model reduction for kinetic equations: moment approximations and hierarchical approximate proper orthogonal decomposition |
| 250 | |a [Electronic ed.] | ||
| 264 | 1 | |c 2021 | |
| 264 | 2 | |b Universitäts- und Landesbibliothek Münster |c 2022-01-18 | |
| 300 | |a xi, 199 | ||
| 502 | |a Zugl.: Münster (Westfalen), Univ., Diss., 2021 | ||
| 505 | 0 | |a Abstract i -- List of Figures ..... ix -- 1 Introduction ..... 1 -- 1.1 Collaborations and previous publications ..... 3 -- 1.2 Notation ..... 4 -- 1.3 Kinetic equations ..... 6 -- 1.3.1 General kinetic transport equations ..... 6 -- 1.3.2 Radiative transfer ..... 8 -- 1.4 Moment models ..... 9 -- 1.4.1 Linear closure ..... 11 -- 1.4.2 Minimum-entropy closures ..... 12 -- 2 First-order finite element moment models ..... 15 -- 2.1 Angular bases in slab geometry ..... 16 -- 2.1.1 Full moments ..... 16 -- 2.1.2 Continuous piecewise linear angular basis (hat functions) ..... 17 -- 2.1.3 Partial moments (discontinuous-Galerkin ansatz) ..... 18 -- 2.2 Angular bases in three dimensions ..... 19 -- 2.2.1 Barycentric-coordinate basis functions (hat functions) ..... 20 -- 2.2.2 Partial moments on the unit sphere ..... 21 -- 2.3 Realizability ..... 22 -- 2.3.1 Structure of the realizable sets ..... 23 -- 2.3.2 Slab geometry (one dimension) ..... 26 -- 2.3.3 Three dimensions ..... 30 -- 2.3.4 Numerically realizable set ..... 31 -- 2.4 Numerical results ..... 34 -- 2.4.1 Slab geometry (one dimension) ..... 35 -- 2.4.2 Three dimensions ..... 39 -- 2.4.3 Computational complexity ..... 41 -- 3 Realizability-preserving finite volume schemes ..... 47 -- 3.1 First-order realizability-preserving finite volume scheme ..... 48 -- 3.1.1 Finite volume grid ..... 48 -- 3.1.2 Finite volume scheme ..... 49 -- 3.2 Second-order realizability-preserving splitting scheme ..... 53 -- 3.2.1 Source system ..... 53 -- 3.2.2 Flux system ..... 54 -- 3.2.3 Time step restriction ..... 56 -- 3.3 Implementation details ..... 59 -- 3.3.1 Solving the optimization problem ..... 59 -- 3.3.2 Solving the eigenvalue problems ..... 65 -- 3.3.3 Realizability limiting ..... 66 -- 3.3.4 Implementation of quadrature rules ..... 70 -- 3.4 Numerical results ..... 71 -- 3.4.1 Slab geometry (one dimension) ..... 72 -- 3.4.2 Three dimensions ..... 76 -- 3.4.3 Quadrature sensitivity ..... 89 -- 3.4.4 Timings ..... 94 -- 4 New numerical scheme based on variable transformation ..... 97 -- 4.1 New scheme in transformed variables ..... 97 -- 4.2 Numerical results ..... 100 -- 4.2.1 Convergence ..... 101 -- 4.2.2 Performance ..... 104 -- 4.2.3 Parallel scaling ..... 111 -- 4.2.4 Masslumping ..... 113 -- 5 Hierarchical approximate proper orthogonal decomposition ..... 117 -- 5.1 Related work ..... 118 -- 5.2 Reduced basis methods ..... 119 -- 5.2.1 General setting ..... 120 -- 5.2.2 Reduced basis method for moment models ..... 121 -- 5.3 Proper orthogonal decomposition ..... 123 -- 5.4 Definition of the HAPOD ..... 126 -- 5.5 Special cases: distributed and incremental HAPOD ..... 127 -- 5.6 Main theorems ..... 128 -- 5.7 Algorithmic benefits ..... 130 -- 5.8 Proofs of main theorems ..... 131 -- 5.9 Reduction of a large linear moment model ..... 136 -- 6 Conclusion and Outlook ..... 143 -- 6.1 First-order finite element-based moment models ..... 143 -- 6.2 Transformed numerical scheme ..... 145 -- 6.3 Reduced-basis methods for moment models ..... 145 -- A Appendix ..... 151 -- A.1 Explicit formulas for integrals ..... 151 -- A.1.1 Mass and flux matrices ..... 151 -- A.1.2 Source-beam boundary values ..... 158 -- A.1.3 Minimum-entropy models in slab geometry ..... 160 -- A.2 Additional figures and data for the variable-transformed scheme ..... 171 -- Bibliography ..... 179 -- Glossaries ..... 191 -- List of Symbols ..... 191 -- List of Acronyms ..... 197. | |
| 506 | 0 | |a free access | |
| 520 | 3 | |a Momentenmodelle sind eine Klasse von spezialisierten approximativen Modellen für kinetische Transportgleichungen. Diese Modelle wandeln die kinetische Gleichung in ein geschwindigkeitsunabhängiges Gleichungssystem für die gewichteten Mittelwerte (Momente) der Lösung um. Je nachdem, welche Gewichtsfunktionen für die Mittelwerte verwendet werden und auf welche Art das Gleichungssystem geschlossen wird, unterscheiden sich die entstehenden Modelle in ihren Eigenschaften. Wird ein linearer Ansatz gewählt, so sind auch die zugehörigen Momentenmodelle linear und mit vergleichsweise geringem Aufwand zu lösen, dafür ist die Lösung aber nicht unbedingt physikalisch sinnvoll und kann zum Beispiel negative Dichten enthalten. Modelle, die auf Entropie-Minimierung beruhen, erhalten hingegen viele der grundlegenden physikalischen Eigenschaften der Lösungen. Dafür liegt der für die Lösung benötigte Rechenaufwand deutlich höher, da ein nichtlineares Optimierungsproblem in jeder Zelle des Zeit-Raum-Gitters gelöst werden muss. Darüber hinaus sind die Optimierungsprobleme nur wohlgestellt, solange die Momente innerhalb eines bestimmten Teilgebiets des Koordinatenraums (der sogenannten Realisierbarkeitsmenge) liegen, was die Entwicklung von numerischen Methoden, insbesondere solcher höherer Ordnung, deutlich erschwert. In dieser Arbeit untersuchen wir deshalb mehrere Ansätze, um die Implementierung und Lösung von (insbesondere Entropie-basierten) Momentenmodellen für lineare kinetische Gleichungen zu vereinfachen. Als erstes konzentrieren wir uns auf die Basisfunktionen (die Gewichtsfunktionen für die Momente). Oft werden diese als Polynome auf dem kompletten Geschwindigkeitsraum gewählt, wohingegen wir uns hier mit stückweise linearen Basisfunktionen beschäftigen. Wir zeigen, dass die entstehenden Modelle für nicht-glatte Lösungen ähnlich gute Approximationseigenschaften haben wie die Standardmodelle, während die Realisierbarkeitsbedingungen sehr viel einfacher und die Rechnungen deutlich schneller sind. Zudem stellen wir ein realisierbarkeitserhaltendes Finite-Volumen-Verfahren zweiter Ordnung vor und zeigen, dass die Implementierung durch die einfachere Struktur der neuen Modelle deutlich weniger komplex ist. Als zweiten Ansatz präsentieren wir ein neues numerisches Verfahren für die Entropie-basierten Momentenmodelle. Durch eine Variablentransformation müssen in dem neuen Verfahren anstelle der nichtlinearen Optimierungsprobleme nur noch positiv definite Matrizen invertiert werden, wodurch die Realisierbarkeitsbedingungen wegfallen und die Lösungen oft deutlich schneller ermittelt werden können als mit Standardverfahren. Schließlich beschäftigen wir uns mit der Möglichkeit einer zusätzlichen Modellreduktion für parameterabhängige Momentenmodelle. Hierzu nutzen wir die Reduzierte-Basis-Methode und berechnen mit der Hauptkomponentenanalyse (Proper Orthogonal Decomposition, POD) ein reduziertes Modell. Wir stellen die hierarchische approximative POD (HAPOD) vor, ein universelles und einfach zu implementierendes Verfahren um eine approximative POD zu berechnen. Die HAPOD nutzt eine fast beliebig an die verfügbaren Ressourcen anpassbare Baumstruktur, so dass jeder Rechenknoten nur die POD einer relativ kleinen Teilmenge der Eingangsvektoren berechnen muss. Wir führen eine ausführliche theoretische Analyse der HAPOD durch und zeigen die Anwendbarkeit auf lineare Momentenmodelle sowie eine gegenüber der POD deutlich verbesserte Effizienz. | |
| 520 | 3 | |a Moment models are a class of specialized approximate models for kinetic transport equations. These models transform the kinetic equation to a system of equations for weighted velocity averages of the solution, called moments, thereby removing the velocity dependency. The properties of the resulting models depend on the chosen weight functions for the moments and on the approach used to close the equations. Closing the system by specifying a linear ansatz function results in linear models that are comparatively easy to solve but may show non-physical behaviour. Minimum-entropy moment closures, on the other hand, conserve many of the fundamental physical properties but require the solution of a non-linear optimization problem for every cell in the space-time grid. In addition, these models are only well-defined if the moments can be kept within a particular subset of the real coordinate space, the so-called realizable set, which is particularly challenging for higher-order numerical solvers. In this thesis, we investigate several approaches to increase usability of (minimum-entropy) moment approximations to linear kinetic equations. First, we focus on the basis functions, i.e. the weight functions for the velocity averages. Typically, these are chosen as polynomials on the whole velocity space. Here, as an alternative, we analyse models based on continuous and discontinuous piecewise linear basis functions. We show that, for non-smooth solutions, approximation quality is comparable while realizability conditions are much simpler and the computations are considerably faster. In addition, we provide second-order finite volume schemes and show that realizability preservation and thus overall implementation is significantly less complex for the new models. Second, we suggest a new numerical method for the minimum-entropy models that is based on a transformation of the moment vectors. This new method replaces the non-linear optimization problems by inversions of positive definite matrices. As a result, it avoids the realizability-related problems and is often several times faster than the standard optimization-based scheme. Finally, we investigate the possibility to further accelerate computations for parameter-dependent moment models by building a reduced model via the reduced basis method utilizing proper orthogonal decomposition (POD). We present the hierarchical approximate POD (HAPOD), a general, easy-to-implement approach to compute an approximate POD based on arbitrary tree hierarchies of worker nodes, where each worker computes a POD of only a small number of input vectors. The tree hierarchy can be freely adapted to optimally suit the available computational resources. We provide rigorous theoretical results and show applicability and performance by applying the HAPOD to a large linear moment model. | |
| 521 | |a specialized | ||
| 536 | |a Förderer: Deutsche Forschungsgemeinschaft / Projektnummer: 390685587 |c 390685587 | ||
| 536 | |a Förderer: Deutsche Forschungsgemeinschaft / Projektnummer: 194347757 |c 194347757 | ||
| 540 | |a InC 1.0 |u https://rightsstatements.org/vocab/InC/1.0/ | ||
| 653 | 0 | |a Momentenmodelle |a Modellreduktion |a HAPOD |a lineare kinetische Gleichungen |a Realisierbarkeit |a Entropieminimierung |a Finite-Volumen-Verfahren | |
| 653 | 0 | |a moment models |a model reduction |a HAPOD |a linear kinetic equations |a realizability |a minimum entropy closure |a finite volume methods | |
| 655 | 7 | |2 DRIVER Types |a Dissertation/Habilitation | |
| 655 | 7 | |2 DCMI Types |a Text | |
| 700 | 1 | |a Ohlberger, Mario |u FB 10: Mathematik und Informatik |0 http://d-nb.info/gnd/1115010972 |0 https://orcid.org/0000-0002-6260-3574 |0 https://www.scopus.com/authid/detail.uri?authorId=6602807638 |0 http://viaf.org/viaf/8208696 |0 http://isni.org/isni/0000000037296360 |0 http://www.wikidata.org/entity/Q59925815 |0 http://id.loc.gov/rwo/agents/n98097714 |4 ths | |
| 787 | 1 | 8 | |i HasPart |o 10.1051/m2an/2021065 |
| 787 | 1 | 8 | |i HasPart |o 10.1016/j.jcp.2020.109547 |
| 856 | 4 | 0 | |3 Zum Volltext |q text/html |u https://nbn-resolving.de/urn:nbn:de:hbz:6-85019477402 |u urn:nbn:de:hbz:6-85019477402 |
| 856 | 4 | 0 | |3 Zum Volltext |q application/pdf |u https://repositorium.uni-muenster.de/document/miami/9952d183-16d5-4d0a-a60e-188c6d8cc844/diss_leibner_2021.pdf |