Das Ziel der vorliegenden Arbeit ist die Entwicklung von partikelbasierenden Strömungs- und Transportmodellen zur Charakterisierung von kleinskaligen Strömungsprozessen in gesättigten und ungesättigten Poren- und Kluftsystemen. Aufgrund der unzureichenden Prozessbeschreibung von ungesättigter Strömung in Doppelkontinuummodellen mittels der Richardsgleichung und van Genuchten Parametern werden innovative Methoden präsentiert um die zugrunde liegenden hochdynamischen Strömungs- und Transportprozesse zu erfassen. Die Simulation von Strömung und Transport in ungesättigten geklüfteten Aquiferen bildet immer noch ein höchst anspruchsvolles Aufgabenfeld aufgrund von skalenübergreifenden Diskontinuitäten, welche oftmals die Definition eines globalen repräsentativen Einheitsvolumens nicht zulassen. Des Weiteren können die hydraulischen Eigenschaften und potentiellen Parameterräume von geklüfteten Aquiferen oftmals nur durch integrale Ansätze, wie z.B. Pump- und Slugtests, Zeitreihenanalysen von Quellschüttungen und Tracertests ermittelt werden. Doppelkontinuummodelle bieten hierfür einen ausgewogenen Ansatz hinsichtlich der erforderlichen Felddaten und der resultierenden prädiktiven Modellqualität. Der erste Teil dieser Arbeit evaluiert den Doppelkontinuumansatz, welcher die Simulation von Strömung mittels der Richardsgleichung und van Genuchten Parametern in zwei, durch einen linearen Austauschterm gekoppelten, Kontinua ermöglicht. Ganglinien von Karstquellen weisen eine charakteristischen steilen Abfall nach Niederschlagsereignissen auf, der durch das Modell erfolgreich reproduziert werden kann. Das Röhrensystem bildet die hydraulische Brücke zur Karstquelle und nimmt potentialabhängige Wassermengen des geklüfteten Matrixsystems auf. Um die Simulation von schneller Grundwasserinfiltration durch das Röhrenkontinuum innerhalb der ungesättigten Zone zu vermeiden wurde die entsprechende Randbedingung an die untere Grenze des Kontinuums gesetzt. Ein genereller Nachteil des Doppelkontinuumsansatz ist die potentielle Mehrdeutigkeit von Modellergebnissen. Der duale Parameterraum in Kombination mit schwierig zu ermittelnden Parametern, führt zur Existenz von mehr als einem kalibrierten Modell, wie durch mehrdimensionale Sensitivitätsanalysen aufgezeigt wird.  Insbesondere in Karstaquiferen bilden Diskontinuitäten, wie z.B. Lösungsdolinen, Klüfte und Störungssysteme, bevorzugte hydraulische Elemente für schnelle vertikale Grundwasserneubildungsprozesse, die oftmals nicht durch volumeneffektive Modellansätze erfasst werden können. Der Hauptteil dieser Arbeit befasst sich daher mit der Entwicklung von zwei Smoothed Particle Hydrodynamics (SPH) Modellen um ein adäquates numerisches Werkzeug zur partikelbasierenden Simulation von kleinskaligen Strömungen mit freien Oberflächen und Transportprozessen bereitzustellen. SPH Modelle ermöglichen eine Eulersche Beschreibung eines Strömungsfelds auf Basis der Navier-Stokes Gleichung und Partikelbewegung mittels klassischer Newtonscher Mechanik. Der gitterlose Modellansatz ermöglicht flexible Simulationen von hochdynamischen Phasengrenzen in ungesättigten Klüften und Porenräumen. Das erste SPH Modell wird eingesetzt um durch Oberflächenspannung dominierte Tropfen- und Filmströmungen auf glatten und rauhen Kluftoberflächen zu simulieren. Charakteristische dimensionslose Kennzahlen werden über einen weiten Bereich von Benetzungswinkeln und Reynoldszahlen bestimmt. Modellergebnisse weisen einen hervorragende Übereinstimmung mit dimensionslosen Skalierungsfunktionen auf und kritische Kontaktwinkel folgen der zu erwartenden Entnetzungsdynamik. Die Entstehung von adsorbierten Filmen auf trockenen Oberflächen wird für einen breiten Parameterraum bestimmt. Des Weiteren wird der Einfluss von befeuchteten Oberflächen auf die Geschwindigkeitszunahme von Tropfenströmung aufgezeigt und so die Bedeutung der Koexistenz verschiedener Strömungsmodi gezeigt. Der Effekt von Oberflächenrauhigkeit auf Tropfenströmung wird für verschiedene Rauhigkeiten ermittelt und eine deutliche Geschwindigkeitsabnahme demonstriert. Um die makroskopische Kontinuumsbeschreibung der Navier-Stokes Gleichung und atomistische Effekte eines klassischen Partikelsystems der statistischen Mechanik zu kombinieren wurde ein zweites mesoskopisches SPH Modell entwickelt. Diese neue Diskretisation der vollständig gekoppelten Landau-Lifshitz-Navier-Stokes und Advektions- Diffusionsgleichung ermöglicht die Simulation von Strömung und Transport bei gleichzeitiger Berücksichtigung von Fluktuationsdynamiken, welche sich korrekt der Systemskala anpassen. Die Verbindung von klassischer Fickscher Diffusion und thermodynamischen Fluktuationen wird hierbei durch einen effektiven Diffusionskoeffizienten beschrieben. Numerische Experimente zeigen die Präzision des Modells. Grenzflächen zwischen zwei Fluiden unterschiedlicher Konzentration weisen eine korrekte Wellenzahldivergenz entsprechend aktuellen Laborergebnissen auf.

Ce travail a pour objectif de développer et d’appliquer une méthode originale permettant de simuler l’écoulement dans un milieu poreux fracturé. Cette méthode repose sur une approche multicouches double continuum permettant de séparer le comportement des différents aquifères présents sur un site. La résolution des écoulements, basée sur la méthode des Eléments Finis de Crouzeix-Raviart, est associée à une méthode inverse (minimisation de type Quasi-Newton combinée à la méthode de l’état adjoint) et à une paramétrisation multi-échelle.La méthode est appliquée dans un premier temps sur l’aquifère fracturé du site expérimental de Poitiers. Les résultats montrent une bonne restitution du comportement de l’aquifère et aboutissent à des champs de transmissivité plus réguliers par rapport à ceux de l’approche simple continuum. L’application finale est réalisée sur le site de Cadarache (taille plus importante et données d’entrée moins denses). Le calage des deux aquifères présents sur le site est satisfaisant et montre que ceux-ci se comportent globalement de façon indépendante. Ce calage pourra être amélioré localement grâce à données de recharge plus fines. / The aim of this study is to develop and validate a new method for the simulation of flow in fractured porous media. This method is based on a multi-layered and dual continuum approach allowing to discriminate the behavior of different aquifers present on a site. The flow equations are solved using a Crouzeix-Raviart Finite Element method, in association with an inverse method (Quasi-Newton minimization combined with the adjoint state method) and a multi-scale parameterization.The method is first applied and validated on the fractured aquifer of the Hydrogeological Experimental Site of Poitiers. The results closely reproduce the flow behavior of the aquifer and lead to a transmissivity field much more homogeneous than the one obtained with a simple continuum approach. The final application is performed on the site of Cadarache (large scale problem with heterogeneously distributed input data). The model calibration of both aquifers is rather satisfactory and shows that their behavior is globally independent. It could locally be improved if more accurate groundwater recharge data is made available.