Das Projekt "ROCKFLOW - Numerische Simulation von Stroemung, Stoff- und Waermetransport im Festgestein" wird vom Umweltbundesamt gefördert und von Universität Hannover, Institut für Strömungsmechanik und Elektronisches Rechnen im Bauwesen durchgeführt. Fuer zahlreiche Fragestellungen im Zusammenhang mit der Deponierung von Abfaellen, der Altlastensanierung, der Grundwassergewinnung und der Nutzung geothermischer Energie sind Modelluntersuchungen im Kluftgestein durchzufuehren. Dabei sind einerseits geringleitende Formationen fuer die Deponierung von Interesse, die eine wirksame geologische Barriere zur Isolierung von Schadstoffen bilden, andererseits geht es um die Bewirtschaftung von Aquiferen zur Grund- und Thermalwassergewinnung. Die Simulation von Stroemungs- und Transportprozessen in klueftig-poroesen Grundwasserleitern und Grundwassergeringleitern stellt spezifische Anforderungen an die modelltechnische Umsetzung, die aus der signifikanten Inhomogenitaet des klueftigen Untergrunds erwachsen. Zur numerischen Simulation solcher Vorgaenge ist das Finite-Elemente-Programmsystem ROCKFLOW entwickelt worden. Das Programm besteht aus einer Reihe von FE-Rechnenkernen (Kernels), welche die prozessspezifischen Differentialgleichungen mittels Galerkin-FEM approximieren. Diese Rechenkerne sind miteinander verknuepfbar (Models), so dass gekoppelte Prozesse (z.B. Tracertransport durch eine Gasstroemung) simuliert werden koennen. Physikalische Prozesse: Folgende physikalische Prozesse sind modellierbar: - Grundwasserstroemung (Sicker- und Kluftstroemung) - Gasstroemung (kompressible Fluide) - Mehrphasenstroemungen (Systeme aus in- und kompressiblen Fluiden) - nicht- (Forchheimer) und liniare Fliessgesetze (Darcy) - hydrodynamische Dispersion (Scheidegger-Ansatz) - Zerfallreaktionen - nicht- (Freundlich, Langmuir) und lineare Gleichgewichtssorption (Henry) - Dichtestroemungen. Numerik: ROCKFLOW ist ein Finite-Elemente-Simulator, wobei verschieden-dimensionale isoparametrische Elemente beliebig im Raum koppelbar sind. Auf der Basis der Methode der gewichteten Residuen wird eine zur prozessbeschreibenden Differentialgleichung aequivalente sog. 'schwache' Integralformulierung abgeleitet. Es stehen verschiedene Loeser zur Verfuegung (Gauss, BiCGSTAB, QMRCGSTAB), um die resultierenden albebraischen Gleichungssysteme zu loesen. Nichtlineare Probleme werden mit Picard- oder Newton-Verfahren behandelt. Gitteradaption: Ab der dritten Version stehen Methoden fuer eine problemangepasste Gitteradaption zur Verguegung. Der Algorithmus zur Gitteradaption basiert auf einem hierarischen Konzept zur Verfeinerung und Vergroeberung gekoppelter verschieden-dimensionaler Elemente. Diskretisierungsfehler koennen entweder mit heuristischen Indikatoren oder einem analytischem Estimator lokalisiert und quantifiziert werden.