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.