Subsurface Flow Module
Subsurface Flow Module
Zur Analyse von Untergrundströmung
Derzeit werden Endlagerstätten für atomaren Abfall gebaut, in denen verbrauchte Brennstäbe für die nächsten Jahrtausende gelagert werden sollen. Dieses Modell zeigt einen hypothetischen Fall, in dem durch einen Riss, in einem mit Brennstäben gefüllten Kanister, radioaktive Substanzen durch einen Spalt in das umgebende Gestein austreten. Dies führt wiederum zu einer Verfüllung des darüber liegenden Tunnels.
Untergrundströmung in gesättigten und teilgesättigten Medien
Das Subsurface Flow Module ist für Ingenieure und Wissenschaftler geeignet, welche Untergrundströmungen oder Strömungen in anderen porösen Medien simulieren und diese außerdem mit anderen Phänomenen (wie z. B. Poroelastizität, Wärmetransport, chemischen Prozessen und elektromagnetischen Feldern) koppeln möchten. Mithilfe dieses Moduls können Grundwasserströmungen, Ausbreitung von Abfallstoffen und dadurch verursachte Erdreichverschmutzung, Gas- und Erdölströmungen in Richtung der Gas- bzw. Erdölquellen sowie Bodensenkungen durch Grundwasserentnahme modelliert werden. Das Subsurface Flow Module ermöglicht es Modelle von unterirdischen Strömungen durch Kanäle sowie durch gesättigte und teilgesättigte poröse Medien bzw. Felsspalten zu erstellen. Außerdem können diese einfach mit Simulationen von Wärmetransport und Stofftransport sowie von geochemischen Reaktionen und Poroelastizität gekoppelt werden. Viele unterschiedliche Branchen und Industriezweige sehen sich durch komplexe Aufgaben im Bereich der Geophysik und Hydrologie konfrontiert. Ingenieure diverser Sparten (z. B. Hoch- und Tiefbau, Bergbau, Ölindustrie, Agrartechnik, chemischer Industrie, Nukleartechnik und Umwelttechnik) müssen häufig derartige Phänomene in ihre Überlegungen mit einbeziehen. Denn die Industriezweige, in denen sie tätig sind, nehmen unter Umweltgesichtspunkten direkten oder indirekten Einfluss auf unsere Erde.
Untergrundströmungen beeinflussen viele verschiedene geophysikalische Größen
Das Subsurface Flow Module enthält eine Reihe von maßgeschneiderten Interfaces, die sogenannten Physikinterfaces, die bestimmte physikalische Phänomene beschreiben. Mithilfe dieser Interfaces können Sie Strömungen aber auch andere Phänomene im Untergrund modellieren. Diese Interfaces können mit anderen Physikinterfaces des Subsurface Flow Module oder mit einem der anderen Module der COMSOL-Produktpalette kombiniert und direkt gekoppelt werden. Auf diese Weise können Sie zum Beispiel das im Subsurface Flow Module beschriebene poroelastische Verhalten mit nichtlinearen Anwendungen der Festkörpermechanik für Böden und Gestein mit dem Geomechanics Module koppeln.
Berücksichtigung geochemischer und kinetischer Reaktionsraten
COMSOL gibt Ihnen absolute Flexibilität, denn Sie können in die Bearbeitungsfelder der Physikinterfaces des Subsurface Flow Module jede beliebige Gleichung eingeben. Das kann z. B. bei der Definition geochemischer und kinetischer Reaktionsraten in den Interfaces für den Materialtransport sehr praktisch sein. Wenn Sie diese Physikinterfaces zusätzlich noch mit dem Chemical Reaction Engineering Module verbinden, können Sie auch Multi-Spezies-Reaktionen mit Hilfe der benutzerfreundlichen Physikinterfaces zur Definition von chemischen Reaktionen modellieren. Eine Integration und Kombination dieser beiden Modulprodukte kann sehr nützlich sein, z. B. bei Modellierung der vielen Reaktionen, die sich bei der Ausbreitung von radioaktivem Abfall über Tausende von Jahren hinweg ergeben.
Weitere Bilder:
-
Typische geoelektrische Flächenelektrodenanordnung auf der Bodenoberfläche. Die Abbildung zeigt neben dem elektrischen Potenzial auch die Sensitivitätsdichte für die Wenner-Alpha-Konfiguration. -
3D-Darstellung der Verdichtung einer Erdöllagerstätte durch Pumpvorgänge. Hier wurde die Darcy-Gleichung mit den durch Poroelastizität verursachten Verschiebungen gekoppelt.
-
Strömung in einem porösen Block mit einer Spalte, in der die Strömung sehr viel stärker ist als im umgebenden porösen Medium. Der Modellierungsansatz verwendet innere Ränder für die Modellierung der freien Strömung in 2D, während die dreidimensionale Modellierung der Strömung im porösen Medium über Darcy-Gleichungen erfolgt.
-
Geometrien mikroskopischer Poren, die durch das Scannen von Bildern aus dem Elektronenmikroskop erzeugt werden, können in COMSOL importiert werden. Anschließend können dann die Geschwindigkeits- und Druckverteilungen berechnet werden. Modell mit freundlicher Genehmigung von Arturo Keller, University of California, Santa Barbara, USA.
-
Grundwasser fließt in eine relativ trockene Erdbodensäule und transportiert dabei eine Chemikalie. Auf dem Weg durch das unterschiedlich gesättigte Medium lagert sich die Chemikalie an Feststoffpartikeln ab, wodurch sich der Transport der gelösten Stoffe gegenüber dem Wasser verlangsamt. Zudem verringert sich die chemische Konzentration durch biologischen Abbau, und zwar sowohl in der flüssigen als auch in der festen Phase. Die Abbildung zeigt die Konturen von Sättigung und Druckhöhe.
-
Abbau eines Pestizids im ungesättigten Boden. Die Abbildung zeigt den biologischen Abbau des Pestizids sowie die Zwischenprodukte im Verlauf der Zeit.
Physikinterfaces für die Simulation von Untergrundströmung
Das Subsurface Flow Module enthält diverse Physikinterfaces zur Simulation von Untergrundströmung und den damit verbundenen Prozessen:
Strömungen durch poröse Medien
Die Hauptanwendung des Subsurface Flow Module besteht in der Modellierung von Strömungen durch teilgesättigte oder gesättigte poröse Medien. In den Physikinterfaces sind die anwendungsspezifischen Ausdrücke (z. B. für Druckhöhe und hydraulische Höhe) bereits vordefiniert. In teilgesättigten Strömungen verändern sich die hydraulischen Eigenschaften während der Durchströmung des Mediums, da einige Poren gefüllt und andere entleert werden. Um diese Effekte miteinzubeziehen, wird die Richard-Gleichung verwendet. Mit der van Genuchten- und der Brooks-Corey-Formulierung kann außerdem die Retention in Poren berücksichtigt werden. Dichte, dynamische Viskosität, Anteil der gesättigten und restlichen Flüssigkeit, hydraulische Leitfähigkeit sowie Speichermodelle können einfach in die entsprechenden Eingabefelder eingegeben werden.
Strömungen durch gesättigte poröse Medien werden in Abhängigkeit von der Porengröße entweder mittels der Darcy-Gleichung oder anhand der Brinkman-Gleichung modelliert. Sofern viskose Effekte auf die Fluidströmung ignoriert werden können, kommt das Darcy-Gesetz zur Anwendung. Die Strömung wird in diesem Fall ausschließlich über die Druckvariable beschrieben. Ist die Porengröße aber so groß, dass das Fluid Impulsänderungen durch Schereffekte verursachen kann, werden Brinkman-Gleichungen verwendet. Berechnet werden dabei zwar dieselben Variablen wie in den Navier-Stokes-Gleichungen, jedoch sind hier zusätzliche Ausdrücke enthalten, die die Porosität des Mediums berücksichtigen, durch welches das Fluid fließt.
Das Spaltenströmung-Interface löst auf zweidimensionalen Rändern innerhalb einer 3D-Matrix nach dem Druck und ist automatisch mit den Gleichungen zur Simulation der Strömung durch die poröse Umgebungsmatrix gekoppelt. Durch diesen Näherungsansatz müssen die eigentlichen Felsspalten nicht vernetzt werden und der Rechenaufwand minimiert sich erheblich. Strömungen eines Fluids von einem Medium in ein anderes (und wieder zurück) können in selben Modell berücksichtigt werden, da im Subsurface Flow Module alle physikalischen Phänomene zur Beschreibung von porösen Medien automatisch mit den Beschreibungen der freien Strömung gekoppelt sind.
Freie Strömungen
Im Untergrund gibt es Strömungskanäle oder große, miteinander verbundene, porenartige Öffnungen, die sich besser durch die Gleichungen für Fluidströmungen modellieren lassen. Auch Wasser- oder andere Quellen, wie beispielsweise bei der Erdölförderung, lassen sich mithilfe von Strömungsgleichungen beschreiben. Das Subsurface Flow Module unterstützt zwei Arten von freien Strömungen: laminare Strömung und Kriechströmung. Das Laminare Strömung-Interface berechnet die Navier-Stokes-Gleichungen, während das Kriechströmung-Interface eine modifizierte Version dieser Gleichungen berechnet und dabei den Trägheitsterm vernachlässigt. Die Kriechströmung, die auch als Stokes-Strömung bezeichnet wird, eignet sich für Strömungen mit sehr niedrigen Reynoldszahlen.
Materialtransport
Materialtransport, welcher in der Regel durch Konvektion oder Diffusion erfolgt, und Untergrundströmung können einfach miteinander gekoppelt werden. Eigenschaften (wie z. B. Diffusion) werden anhand von Gleichungen beschrieben, die von Variablen abhängig sind (wie z. B. der Konzentration) oder anisotrop sind.
Das Transport, verdünnte Spezies-Interface erweitert den Transportmechanismus um Dispersion und Retardation, und zwar aufgrund von Sorption. Der Transport gelöster Stoffe erfolgt meist in Fließrichtung. Es steht ein Dispersionstensor zur Verfügung, um die Anisotropie miteinbeziehen zu können. Die Sorption beschreibt den Prozess von chemischen Substanzen, die inerhalb des porösen Mediums unterschiedlich schnell absorbiert und dann wieder desorbiert werden. Die Auswirkungen auf den Materialtransport werden entweder durch die Langmuir- oder Freundlich-Isotherme beschrieben, die beide in diesem Interface verfügbar sind, oder auch durch einen benutzerdefinierten Ausdruck. Sorption führt außerdem durch den Retardationsfaktor zur Verlangsamung der Strömung. Für nicht gesättigte Strömungen steht zudem ein Verdampfungs- oder Dispersionsfaktor für die, aus dem gelösten Zustand zu Gas werdenden, chemischen Stoffen zur Verfügung. Darüber hinaus können zur Beschreibung beliebiger Reaktionen, die während des Materialtransports stattfinden, eigenen Gleichungen eingegeben werden. Das Transport, verdünnte Spezies-Interface kann mit jedem anderen geeigneten Physikinterface im CFD Module gekoppelt werden, z. B. wenn Sie Zweiphasenströmungen und die dadurch transportierten gelösten Stoffen modellieren möchten.
Wärmetransport
Der Wärmetransport erfolgt durch Konduktion, Konvektion und Dispersion. Dabei ist die unterschiedliche Leitfähigkeit von flüssiger und fester Phase zu berücksichtigen. In vielen Fällen besteht die feste Phase aus verschiedenen Materialien mit unterschiedlicher Leitfähigkeit, aber auch Fluide können unterschiedlicher Natur sein. Das Wärmetransport in porösen Medien-Interface bietet die Möglichkeit, z. B. durch vordefinierte Mischregeln, effektive Wärmetransporteigenschaften berechnen zu können. Außerdem sind viele weitere Ausdrücke hinterlegt, so z. B. für Wärmedispersion, aufgrund der natürlichen Beschaffenheit des porösen Mediums, oder für geothermische Hintergrunderwärmung. Dispersionseffekte, welche insbesondere aufgrund der gewundenen Transportwege in porösen Medien auftreten, werden durch die alleinige Annahme eines gemittelten konvektiven Terms nicht berücksichtigt.
Poroelastizität
Das leistungsfähige Poroelastizität-Interface ermöglicht die Modellierung von Verdichtung und Absenkung. Dabei wird eine transiente Formulierung des Darcy-Gesetzes mit der linear elastischen Strukturmechanik für das poröse Medium kombiniert. Durch diese Poroelastizitätskopplung kann der Einfluss der Strömung auf die Kompressibilität des porösen Mediums betrachtet werden. Änderungen in den volumetrischen Dehnungen wiederum beeinflussen Impuls- sowie auch Material- und Wärmetransport. Dieses Interface enthält einen Ausdruck für den Spannungstensor, und zwar als Funktion von Dehnungstensor und Biot-Willis-Koeffizient.
Pumping and Injecting from a Single Borehole
Prof. Martin Sauter, Ph.D Student Yulan Jin, and Assoc. Prof. Ekkehard Holzbecher Georg-August-University, Göttingen, Germany
Professor Dr. Martin Sauter leads a group at the University of Göttingen in research on hydrogeology. Cooperating with Hölscher Wasserbau, a leading German dewatering company, the group is studying düsensauginfiltration (DSI), which is a technique for lowering groundwater levels at mining or construction sites. Traditionally, dewatering involves ...
Multiphysics Simulations Help Track Underground Fluid Movements
A. Haas, H. Mahardika, M. Karaoulis, and A. Revil Colorado School of Mines, USA
Combining acoustic and electromagnetic analyses helps scientists to more efficiently measure and survey subterranean regions. Acoustic waves can travel long distances, but have limitations when it comes to providing details about formation properties, and cannot be used to identify the liquids flowing through them. Electromagnetic waves do not ...
Phase Change
This example demonstrates how to model a phase change and predict its impact on a heat transfer analysis. When a material changes phase, for instance from solid to liquid, energy is added to the solid. Instead of creating a temperature rise, the energy alters the material’s molecular structure. Equations for the latent heat of phase changes ...
Pore-Scale Flow
This non-conventional model of porous media flow utilizes creeping (Stokes) flow in the interstices of a porous media. The model comes from the pore-scale flow experiments conducted by Arturo Keller, Maria Auset, and Sanya Sirivithayapakorn of the University of California, Santa Barbara. The geometry used in the model was produced by scanning ...
Biot Poroelasticity
The Poroelasticity interface couples Darcy's law and solid mechanics to assess deformation of porous media that results from fluid withdrawals. The model builds on top of the Terzaghi Compaction example. Results from Terzaghi compaction and Biot poroelasticity analyses are compared to each other and are a good match to published analyses.
Failure of a Multilateral Well
Multilateral wells—those with multiple legs that branch off from a single well—can produce oil efficiently because the legs can tap multiple productive zones and navigate around impermeable ones. Unfortunately, drilling engineers must often mechanically stabilize multilateral wells with a liner or casing, which can cost millions of dollars. ...
Forchheimer Flow
This is a tutorial model of the coupling between flow of a fluid in an open channel and a porous block attached to one of the channel walls. The flow is described by the Navier-Stokes equation in the free region and a Forchheimer-corrected version of the Brinkman equations in the porous region.
Free Convection in a Porous Medium
This example treats the modeling of sub-surface flow where free convection in porous media is analyzed. The results are compared with published literature in the field. The model couples the momentum balance to an energy balance through an equation, dependent on temperature, being directly typed into the source term for the momentum balance. ...
Pesticide Transport and Reaction in Soil
Aldicarb is a commercial pesticide, used on a variety of crops, including cotton, fruits, potatoes, and beans. This arises in the possibility that the general population may be exposed to aldicarb through the ingestion of contaminated water and foods. This example looks at the degradation kinetics of aldicarb and its toxic by-products, ...
Variably Saturated Flow and Transport—Sorbing Solute
In this example water ponded in a ring on the ground moves into a relatively dry soil column and carries a chemical with it. As it moves through the variably saturated soil column, the chemical attaches to solid particles, slowing the solute transport relative to the water. Additionally the chemical concentrations decay from biodegradation in ...
Perforated Well
The ability to describe how fluids funnel into tiny perforations oriented about a well bore is the subject of a rapidly growing number of analyses. In that the perforations are isolated piercings as opposed to rings, the flow field is not suited for axisymmetric analyses - fully 3-D simulations are required. This application allows to simulate ...
Buoyancy Flow with Darcy's Law—The Elder Problem
Density variations can initiate flow even in a still fluid. In earth systems, density variations can arise from naturally occurring salts, subsurface temperature changes, or migrating pollution. This buoyant or density-driven flow factors into fluid movement in salt-lake systems, saline-disposal basins, dense contaminant and leachate plumes, and ...
