Kromer, Johannes ; Potyka, Johanna ; Schulte, Kathrin ; Bothe, Dieter (2023)
Efficient sequential PLIC interface positioning for enhanced performance of the three-phase VoF method.
In: Computers & Fluids, 266
doi: 10.1016/j.compfluid.2023.106051
Artikel, Bibliographie
Kurzbeschreibung (Abstract)
This paper presents an efficient algorithm for the sequential positioning, also called nested dissection, of two planes in an arbitrary polyhedron. Two planar interfaces are positioned such that the first plane truncates a given volume from this arbitrary polyhedron and the next plane truncates a second given volume from the residual polyhedron. This is a relevant task in the numerical simulation of three-phase flows when resorting to the geometric Volume-of-Fluid (VoF) method (Hirt and Nichols, 1981) with a Piecewise Linear Interface Calculation (PLIC) (Youngs, 1982). An efficient algorithm for this task significantly speeds up the three-phase PLIC algorithm. The present study describes a method based on a recursive application of the Gaussian divergence theorem, where the fact that the truncated polyhedron shares multiple faces with the original polyhedron can be exploited to reduce the computational effort. A careful choice of the coordinate system origin for the volume computation allows for successive positioning of two planes without reestablishing polyhedron connectivity. Combined with a highly efficient root finding, this results in a significant performance gain in the reconstruction of the three-phase interface configurations.
The performance of the new method is assessed in a series of carefully designed numerical experiments. Compared to a conventional decomposition-based approach, the number of iterations and, thus, of the required truncations was reduced by up to an order of magnitude. The PLIC positioning run-time was reduced by about 90% in our reference implementation. Integrated into the multi-phase flow solver Free Surface 3D (FS3D), an overall performance gain of about 20% was achieved. The reference implementation of the efficient sequential PLIC positioning is available as a Fortran module at https://doi.org/10.18419/darus-2488. Allowing for simple integration into existing numerical schemes, the proposed algorithm is self-contained, requiring no external decomposition libraries.
| Typ des Eintrags: | Artikel |
|---|---|
| Erschienen: | 2023 |
| Autor(en): | Kromer, Johannes ; Potyka, Johanna ; Schulte, Kathrin ; Bothe, Dieter |
| Art des Eintrags: | Bibliographie |
| Titel: | Efficient sequential PLIC interface positioning for enhanced performance of the three-phase VoF method |
| Sprache: | Englisch |
| Publikationsjahr: | 2023 |
| Ort: | Amsterdam |
| Verlag: | Elsevier |
| Titel der Zeitschrift, Zeitung oder Schriftenreihe: | Computers & Fluids |
| Jahrgang/Volume einer Zeitschrift: | 266 |
| DOI: | 10.1016/j.compfluid.2023.106051 |
| URL / URN: | https://www.sciencedirect.com/science/article/abs/pii/S00457... |
| Kurzbeschreibung (Abstract): | This paper presents an efficient algorithm for the sequential positioning, also called nested dissection, of two planes in an arbitrary polyhedron. Two planar interfaces are positioned such that the first plane truncates a given volume from this arbitrary polyhedron and the next plane truncates a second given volume from the residual polyhedron. This is a relevant task in the numerical simulation of three-phase flows when resorting to the geometric Volume-of-Fluid (VoF) method (Hirt and Nichols, 1981) with a Piecewise Linear Interface Calculation (PLIC) (Youngs, 1982). An efficient algorithm for this task significantly speeds up the three-phase PLIC algorithm. The present study describes a method based on a recursive application of the Gaussian divergence theorem, where the fact that the truncated polyhedron shares multiple faces with the original polyhedron can be exploited to reduce the computational effort. A careful choice of the coordinate system origin for the volume computation allows for successive positioning of two planes without reestablishing polyhedron connectivity. Combined with a highly efficient root finding, this results in a significant performance gain in the reconstruction of the three-phase interface configurations. The performance of the new method is assessed in a series of carefully designed numerical experiments. Compared to a conventional decomposition-based approach, the number of iterations and, thus, of the required truncations was reduced by up to an order of magnitude. The PLIC positioning run-time was reduced by about 90% in our reference implementation. Integrated into the multi-phase flow solver Free Surface 3D (FS3D), an overall performance gain of about 20% was achieved. The reference implementation of the efficient sequential PLIC positioning is available as a Fortran module at https://doi.org/10.18419/darus-2488. Allowing for simple integration into existing numerical schemes, the proposed algorithm is self-contained, requiring no external decomposition libraries. |
| Zusätzliche Informationen: | Artikel-ID: 106051 |
| Fachbereich(e)/-gebiet(e): | DFG-Sonderforschungsbereiche (inkl. Transregio) DFG-Sonderforschungsbereiche (inkl. Transregio) > Sonderforschungsbereiche DFG-Sonderforschungsbereiche (inkl. Transregio) > Sonderforschungsbereiche > SFB 1194: Wechselseitige Beeinflussung von Transport- und Benetzungsvorgängen DFG-Sonderforschungsbereiche (inkl. Transregio) > Sonderforschungsbereiche > SFB 1194: Wechselseitige Beeinflussung von Transport- und Benetzungsvorgängen > Projektbereich B: Modellierung und Simulation DFG-Sonderforschungsbereiche (inkl. Transregio) > Sonderforschungsbereiche > SFB 1194: Wechselseitige Beeinflussung von Transport- und Benetzungsvorgängen > Projektbereich B: Modellierung und Simulation > B01: Modellierung und VOF-basierte Simulation der Multiphysik irreversibler thermodynamischer Transferprozesse an dynamischen Kontaktlinien |
| Hinterlegungsdatum: | 07 Dez 2023 13:53 |
| Letzte Änderung: | 07 Dez 2023 13:53 |
| PPN: | |
| Export: | |
| Suche nach Titel in: | TUfind oder in Google |
 |
Frage zum Eintrag |
Optionen (nur für Redakteure)
 |
Redaktionelle Details anzeigen |
Drucken
Impressum