Simulation of two-phase flow processes in heterogeneous porous media with adaptive methods

In recent years the demand for numerical simulations which help to support the work of engineers has constantly gained weight. In the same way the computational power advances, the complexity of the problems and demands posed to the numerical simulation tools increases, too. This thesis expands the simulation tool MUFTE-UG for two-phase flow processes in porous media with adaptive methods The approaches presented here can be divided into two parts: On the one hand, a space adaptive method is introduced. Here, the element-size resolution throughout the discretization mesh is changed in certain areas from time step to time step. A marking of the relevant elements which need to be refined or coarsened is realized by an empirically derived error indicator. While at first three different indicators are compared, the final investigations of the various test cases are performed with an indicator locating a steep gradient of the saturation distribution in the system. As it shows that the applied h-adaptive strategy is not mass conservative due to mesh manipulations, an algorithm is developed which ensures the mass preservation for coarsening as well as for refinement. On the other hand, the discretization method is adaptively adjusted inside the domain. Since advection-dominated processes require an other numerical treatment than diffusiondominated processes (in the here presented case realized by the application of a ’fully upwinding’ or a ’centrally weighted’ scheme), at first a suitable indicator needs to be found which accounts for these processes. For this, the two-phase (element) Pecletnumber is derived which describes the ratio between advection and diffusion. The developed methods are applied to homogeneous and heterogeneous test cases. For the choice of the homogeneous test cases it is considered that the here introduced methods need to be capable of handling purely advection-dominant problems (e.g. Buckley-Leverett problem), purely diffusion-dominant problems (e.g. McWhorter problem) and problems, where both effect appear in the domain at the same time (e.g. Sandbox problem). The heterogeneous test case resembles the Sandbox problem with a lense. Overall it can be said that the deployed space adaptive methods work very well for all the investigated test cases. The results obtained by the adaptive choice of discretization method are only mildly satisfactory. Here it shows, especially for the heterogeneous case, that the switch to the centrally weighted scheme needs a very careful adjusting. In den vergangenen Jahren hat die Nachfrage nach numerischen Simulationsprogrammen, die eine Unterstutzung der taglichen Arbeit von Ingenieuren bieten, standig an Gewicht zugenommen. Im gleichen Mase in dem die Rechnerleistung zunimmt, steigen auch die Komplexitat der Fragestellungen und die Anforderungen an die Simulationsprogramme. Diese Arbeit erweitert das Simulationsprogramm MUFTE-UG fur Zwei-Phasen Stromungsprozesse in porosen Medien um adaptive Methoden. Die hier vorgestellten Vorgehensweisen konnen in zwei Bereiche unterteilt werden: Auf der einen Seite wird eine orts-adaptive Methode vorgestellt. Dabei wird die Elementgrose in bestimmten Gebieten des Diskretisierungsnetzes von Zeitschritt zu Zeitschritt angepasst. Das Markieren der relevanten Elemente, die verfeinert oder vergrobert werden sollen, wird mittels eines empirischen Fehlerindikators realisiert. Zunachst werden drei verschiedene Indikatoren verglichen. Die letztendlichen Untersuchungen werden mit einem Indikator durchgefuhrt, der einen steilen Gradienten der Sattigung im Gebiet lokalisiert. Da sich zeigt, dass das verwendete h-adaptive Verfahren aufgrund von Netz-Manipulationen nicht massenkonservativ ist, wird ein Algorithmus entwickelt der die Erhaltung der Masse fur das Vergrobern und das Verfeinern sicherstellt. Auf der anderen Seite findet ein adaptives Anpassen der Diskretisierungsmethode innerhalb des Gebietes statt. Da advektions-dominante Prozesse eine andere numerische Behandlung als diffusions-dominante Prozesse benotigen (was in dieser Arbeit durch Anwendung der ’fully upwinding’ oder ’zentral gewichteten’ Methode geschieht), muss zunachst ein geeigneter Indikator gefunden werden, der diese Prozesse identifiziert. Hierfur wird die Zwei-Phasen (element-basierte) Peclet-Zahl hergeleitet, die das Verhaltnis zwischen Advektion und Diffusion beschreibt. Die entwickelten Methoden werden fur homogene und heterogene Testfalle eingesetzt. Bei der Wahl der homogenen Falle wurde berucksichtigt, dass die vorgestellten Methoden in der Lage sein sollen sowohl rein advektions-dominante Probleme (z.B. Buckley-Leverett Problem) und rein diffusions-dominante Problem (z.B. McWhorterProblem) als auch Probleme behandeln konnen soll, bei denen beide Effekte gleichzeitig im Gebiet auftreten (z.B. Sandbox Problem). Als heterogenes Test wird das Sandbox Problem mit einer Linse gewahlt. Zusammenfassend kann gesagt werden, dass die Ergebnisse der verwendeten ortsadaptiven Methoden fur alle Testfalle sehr gut sind. Die Ergebnisse, die mit der adaptiven Anpassung der Diskretisierungsmethode erreicht wurden sind nur masig zufriedenstellend. Hier zeigt sich, dass, insbesondere fur den heterogenen Fall, der Ubergang zur zentral gewichteten Methode sehr sorgfaltig gewahlt werden muss.