Massiv parallele Simulation von Mehrphasen- und Mehrkomponentenströmungen unter Anwendung des Lattice Boltzmann Verfahrens

This thesis reflects the work mainly performed within the research project FIMOTUM focusing on the determination of transport properties and mechanisms in unsaturated media. The efficient simulation of single- and multiphase flows at the pore scale in highly resolved natural porous media is one of the major topics in this work. For this purpose a simulation kernel which is based on the lattice Boltzmann method (LBM) has been developed and extensively validated. The LBM presented utilizes the Multiple Relaxation Time (MRT) model and fluid/wall boundary conditions of second order accuracy. The model has also been extended to solve multiphase, advection/diffusion and thermal flow problems. Due to the application of an optimized collision model and corresponding boundary conditions, the covered parameter space and the stability of the method could be greatly enhanced. Hence, it was possible to perform simulations in complex geometries at a large scale (2E11+ DoF) which have been obtained with an unprecedented accuracy. A second target of this thesis was the design and implementation of a simulation kernel to perform massively parallel computations with high efficiency. In order to obtain accurate simulation results at reasonable computational effort, a novel grid generation procedure has been developed. The robust and flexible method is based on the decoupling of input geometry and the actual computational grid. It is therefore excellently suited for the grid generation based on natural porous media data sets obtained by CT- or X-ray methods. Aspects concerning the increasing difficulties in pre- and post-processing of large data sets are discussed. Furthermore, special issues in high performance computing environments are highlighted and a tool chain to visualize scientific data in photo-realistic representation is described.Die vorliegende Dissertation gibt im Wesentlichen die Arbeiten wieder, die im Rahmen des FIMOTUM Projektes durchgeführt worden sind, welches sich vornehmlich auf die Untersuchung von Transporteigenschaften in ungesättigten porösen Medien fokussierte. Hierfür wurde ein Software-Prototyp auf Basis der Gitter Boltzmann Methode (LBM) entwickelt und ausführlich validiert. Die vorgestellte LB-Methode basiert auf dem Multiple-Relaxation-Time (MRT) Modell und verwendet Fluid/Wand Randbedingungen mit einer Genauigkeit 2. Ordnung. Das beschriebene Modell wurde zudem für die Simulation von Mehrphasen-, Advektion/Diffusions- und Thermalen Problemen erweitert. Durch die Optimierung des Kollisionsmodells und der entsprechenden Randbedingungen konnte der nutzbare Parameterraum deutlich vergrößert werden, so dass Simulationen in komplexen Geometrien mit mehr als 2.0E+11 Freiheitsgraden möglich wurden. Ein zweites Ziel dieser Arbeit war die Implementierung eines effizienten und hochparallelen Software-Prototypen zur Simulation von fluiddynamischen Problemen. Um möglichst genaue Ergebnisse bei mäßigem Ressourceneinsatz zu erzielen, wurde ein neuartiger Gittergenerierungsprozess entwickelt. Dieses robuste und flexible Verfahren basiert auf der Entkopplung von Eingangsgeometrie und dem eigentlichen Rechengitter. Daher eignet sich dieser Gittergenerator hervorragend für die Erzeugung eines numerischen Gitters aus digitalen Datensätzen natürlicher poröser Medien, wie bspw. Tomographie-Scans. Desweiteren werden, neben allgemeinen Problemen des Hochleistungsrechnens, die zunehmenden Schwierigkeiten bei der Verarbeitung der ständig steigenden Datenmengen im Pre- und Postprocessing diskutiert. Weiterhin wird, unterstützend zur Ergebnisanalyse, eine Prozesskette für die Erzeugung von fotorealistischen Visualisierungen aus Simulationsdaten beschrieben