Parallel multi-scale simulations with octrees and coupled applications

Klimach, Harald; Behr, Marek; Resch, Michael; Roller, Sabine P.

doi:HT019121639

Parallel multi-scale simulations with octrees and coupled applications = Parallele Multi-Skalen Simulationen mit Octrees und gekoppelte Anwendungen

Klimach, Harald^{GRS Aachen*}

2016

VerantwortlichkeitsangabeHarald Klimach

Ausgabe1st edition

ImpressumSiegen : universi - Universitätsverlag Siegen 2016

Umfang1 Online-Ressource (xv, 185 Seiten) : Illustrationen, Diagramme

ISBN978-3-936533-82-8

ReiheSimulation Techniques in Siegen ; 1

Dissertation, Rheinisch-Westfälische Technische Hochschule Aachen, 2016

Auch veröffentlicht auf dem Publikationsserver der RWTH Aachen University

Genehmigende Fakultät
Fak04

Hauptberichter/Gutachter
Roller, Sabine P. (Thesis advisor) ; Behr, Marek (Thesis advisor) ; Resch, Michael (Thesis advisor)

Tag der mündlichen Prüfung/Habilitation
2016-01-28

Online
URN: urn:nbn:de:hbz:82-rwth-2016-073855
URL: http://publications.rwth-aachen.de/record/670169/files/670169.pdf
URL: http://publications.rwth-aachen.de/record/670169/files/670169.pdf?subformat=pdfa

Einrichtungen

Inhaltliche Beschreibung (Schlagwörter)
Hochleistungsrechnen (frei) ; Oktalbaum (frei) ; numerische Strömungssimulation (frei) ; Gittererzeugung (frei) ; Strömungsakustik (frei)

Thematische Einordnung (Klassifikation)
DDC: 620

Kurzfassung
Simulationen physikalischer Gegebenheiten müssen oft das Zusammenspiel vieler Phänomene und Skalen berücksichtigen. Bei aeroakustischen Problemen zum Beispiel, muss beides berücksichtigt werden, sowohl die Strömung,die den Lärm erzeugt, als auch der Transport der Schallwellen. In dieser Arbeit werden numerische Ansätze für solche Probleme auf großen, verteilt parallelen Rechensystemen untersucht. Das Kopplungsframework KOP wird, soweit wie möglich, parallelisiert und ein neues Framework (APES) wird entwickelt um fundamentale Beschränkungen der Skalierbarkeit zu überwinden. In beiden Implementierungen werden Verfahren hoher Ordnung eingesetzt, da diese eine hochauflösende Simulation mit weniger Freiheitsgraden ermöglichen, als Verfahren niedrigerer Ordnung. Diese Eigenschaft von Verfahren hoher Ordnung ist ein wichtiger Vorteil auf modernen Supercomputersystemen, da der Speicher, der benötigt wird um die Freiheitsgradeabzubilden eine knappe Ressource darstellt. Die vorgestellten Methoden ermöglichen die transiente Simulation von Mehrskalenproblemen, allerdings werden für detaillierte Simulationen noch immer große Mengen an Rechenressourcen benötigt. Im Rahmen dieser Arbeit wird deshalb ein Fokus auf die effiziente Nutzung moderner Rechensysteme gelegt. KOP verwendet diskrete Punkte um die Kopplung zu realisieren. Dies erlaubt die Interaktion zwischen Gebieten mit unterschiedlicher Diskretisierung aber auch verschiedenen Gleichungen. Beide Implementierungen verwenden eine explizite Zeitintegration um die zeitabhängigen Simulationen aufzulösen. Das Kopplungsframework erlaubt von Gebiet zu Gebiet variierende Zeitschritte. Diverse Erhaltungsgleichungen, von linearisierten Euler Gleichungen und den Maxwell Gleichung, bis hin zu den Navier-Stokes Gleichungen, können mit den vorgestellten Verfahren gelöst werden. Ein vollständig verteilter Kopplungsmechanismus wird im Rahmen der Arbeit entwickelt, der es ermöglicht, diese Gleichungen in großen Simulationen gekoppelt zu verwenden. APES erlaubt die Verwendung spektraler Diskretisierungen. Dazu bringt es eine eigene Toolchain mit, die eine skalierbare Ausführung der gesamten Simulation sicherstellt. Insbesondere, umfasst dies auch einen Gittergenerator, der Geometrien mit Polynomen hoher Ordnung darstellen kann. Das robuste Verfahren, das hier zum Einsatz kommt ermöglicht es, ingenieurtechnische Fragestellungen auch mit solchen Verfahren hoher Ordnung in Angriff zu nehmen.

Physical simulations often require the consideration of many phenomena and scales. For example in aeroacoustic problems, both, the flow generatingthe noise and the sound wave propagation needs to be considered. This work investigates numerical approaches to such problems on large distributed and parallel computing systems. The coupling framework KOP is parallelized as far as possible and to overcome fundamental scalability limits a new framework APES is developed. Both implementations utilize high-order discretizations, as these allow for accurate simulations with less degrees of freedoms than lower order methods. This property of high-order methodsis an important feature for modern supercomputing systems, as memory to represent degrees of freedom in a simulation is a scarce resource. The presented methods enable the transient simulation of multi-scale setups but detailed resolutions still require large amounts of computational resources. A focus is put on the efficient utilization of modern computing systems to address this need. Besides the scalability of the implementations, the importance of single core optimization and vectorization is illustrated. KOP uses discrete points to realize the coupling and allows for the interaction between domains with differing discretizations and solved equation systems. Arbitrary mesh configurations are supported and both, structuredand unstructured mesh solvers are available in the framework. In both framworksexplicit time integration methods are deployed to resolve the timedependent simulations. The coupling allows for a varying time step widthover the participating domains by a sub-cycling method. Various conservationlaws can be solved by the presented frameworks ranging from Maxwell’s equations and linearized Euler equations to full compressible Navier-Stokesequations. A fully distributed coupling approach is developed that allows for coupling of those in a large-scale simulation to solve, for example, aeroacousticproblems.APES enables high-order discretizations in the spectral regime. It involvesa fully scalable toolchain for mesh-based simulations featuring a meshgeneration and a post-processing tool to support the solvers. The common foundation of these tools is an Octree representation for the mesh, andthis work specifically covers the generation of high-order geometry approximationsin the developed mesh generator Seeder. This robust mechanismworks for arbitrarily complex surfaces and offers a practical way to tackleengineering tasks with spectral element discretizations.

OpenAccess:
PDF PDF (PDFA)
(additional files)