The styles of deformation of the fore-arc wedges along the Chilean convergent margin are observed to vary significantly, despite similar plate kinematic conditions. Here, I focus on the analysis of fore-arc deformation on two regions along the Chilean convergent margin at 20°-24°S and 37°-42°S. Although both regions are subjected to the oblique subduction of the oceanic Nazca plate and backstopped by the Andes mountain chain; they display different patterns of deformation. The northern Chilean study area (20° - 24°S) is characterized by an exceptionally thick crust of about 60 km beneath the Altiplano – Puna plateau, lack of an accretionary wedge in the fore-arc due to hyperarid climate, and consequently a sediment starved trench. Two major margin parallel strike slip faults are observed in this area, the Atacama Fault Zone (AFZ) and the Precordilleran Fault System (PFS). Both strike-slip faults do not exhibit significant recent displacement. The southern study area (37° - 42°S), compared to the northern study area, is characterized by lower topography, high precipitation rates (~2000 mm/yr), and a younger subducted oceanic plate. An active strike-slip fault, the Liquiñe-Ofqui-Fault-Zone (LOFZ), shows ~1 cm/yr recent dextral movement and shapes the surface of this area. Thus, the southern Chilean study area exhibits localized strike-slip motion. Within this area the largest earthquake ever recorded, the 1960 Valdivia earthquake, occurred with a moment magnitude of MW=9.5. I have constructed 2D thermal models and 3D mechanical models for both Chilean study areas to study processes related to active subduction. The applied numerical method is the finite element technique by means of the commercial software package ABAQUS. The thermal models are focused on the thermal conditions along the plate interface. The thermal structure along the plate interface reveals the limits of coupling but also the type of transition from coupled to uncoupled and vice versa. The model results show that shear heating at the plate interface is an important mechanism that should be taken into account. The models also show that the thermal condition at the downdip limit of the coupling zone leads to a sharp decrease of friction along the interface. Due to the different geometries of the two Chilean study areas, such as the slab dip and the thickness of the continental crust, the downdip limit of the southern study area is slightly shallower than that of the northern study area. The results of the 2D thermal models are used to constrain the spatial extent of the coupling zone in the 3D mechanical models. 3D numerical simulations are used to investigate how geometry, rheology and mechanical parameters influence strain partitioning and styles of deformation in the Chilean fore-arc. The general outline of the models is based on the fore-arc geometry and boundary conditions as derived from geophysical and geological field data. I examined the influence of different rheological approaches and varying physical properties of the fore-arc to identify and constrain the parameters controlling the difference in surface deformation between the northern and southern study area. The results of numerical studies demonstrate that a small slab dip, a high coefficient of basal friction, a high obliquity of convergence, and a high Young’s modulus favour localisation of deformation in the fore-arc wedge. This parameter study helped me to constrain preferred models for the two Chilean study areas that fit to first order observations. These preferred models explain the difference in styles of deformation as controlled by the angle of obliquity, the dip of subducting slab, and the strength of wedge material. The difference in styles can be even larger if I apply stronger coupling between plates within the southern area; however, several independent observations indicate opposite tendency showing southward decrease of intensity of coupling. The weaker wedge material of the preferred model for the northern study area is associated with advanced development of the adjacent orogen, the Central Andes. Analysis of world-wide examples of oblique subduction zones supports the conclusion that more mature subduction zones demonstrate less pronounced localization of strike-slip motion. / Die Deformationsmuster der Fore-Arc Keile entlang des chilenischen konvergenten Plattenrandes variieren beachtlich, trotz ähnlicher plattenkinematischer Randbedingungen. In dieser Arbeit konzentriere ich mich auf die Analyse der Deformation des Fore-Arcs in zwei Gebieten entlang des chilenischen konvergenten Plattenrandes zwischen 20°-24°S und 37°-42°S. Obwohl beide Gebiete durch schiefe Subduktion der ozeanischen Nazca Platte und der östlichen Begrenzung durch die Andine Gebirgskette gekennzeichnet sind, zeigen sie unterschiedliche Deformationsmuster an der Oberfläche. Das nördliche chilenische Gebiet (20° - 24°S) ist gekennzeichnet durch eine außergewöhnliche Krustendicke von ~ 60 km unterhalb des Altiplano - Puna Plateaus, dem Fehlen eines akkretionären Prismas im Fore-Arc aufgrund des trockenen Klimas und somit einer nahezu sedimentfreien Tiefseerinne. Zwei große Plattenrand-parallele Strike-Slip Störungen werden in diesem Gebiet beobachtet, die Atacama Fault Zone (AFZ) und das Precordilleran Fault System (PFS). Beide Strike-Slip Störungen zeigen keine signifikanten aktuellen Bewegungsraten. Das südliche Gebiet (37° - 42°S) ist im Vergleich zum nördlichen Gebiet durch eine niedrigere Topographie, hohe Niederschlagsraten (~2000 mm/a) und eine jüngere abtauchende ozeanische Platte gekennzeichnet. Die aktive Strike-Slip Störung, Liquiñe-Ofqui-Fault-Zone (LOFZ), ist gekennzeichnet durch aktuelle dextrale Bewegungsraten von 1 cm/a und prägt die Oberflächenstruktur in dieser Region. Folglich ist der südliche Arbeitsbereich durch lokalisierte Strike-Slip Bewegung charakterisiert. Innerhalb dieses Gebietes ereignete sich das größte instrumentell aufgezeichnete Erdbeben, das 1960 Valdivia Erdbeben, mit einer Stärke von MW=9.5. 2D thermische Modelle und 3D mechanische Modelle wurden für die beiden chilenischen Gebiete konstruiert, um Prozesse im Zusammenhang mit aktiver Subduktion zu untersuchen. Als numerisches Verfahren wurde die Finite Elemente Methode mit Hilfe des kommerziellen Softwarepakets ABAQUS angewandt. Die thermischen Modelle sind auf die thermischen Konditionen entlang der Plattengrenzfläche fokussiert. Die thermische Struktur entlang der Plattengrenzfläche zeigt die Grenzen der Kopplung an aber auch die Art des Überganges von gekoppelt zu nicht gekoppelt und umgekehrt. Die Modellergebnisse zeigen, dass Heizen infolge der Scherung an der Plattengrenzfläche ein wichtiger Faktor ist, der in Betracht gezogen werden sollte. Die Modelle zeigen auch, dass die thermische Struktur an der unteren Begrenzung der Koppelzone zu einer deutlichen Abnahme der Reibung entlang der Grenzfläche führt. Aufgrund der unterschiedlichen Geometrien der zwei chilenischen Untersuchungsgebiete, z.B. Abtauchwinkel der ozeanischen Platte und Krustendicke, ist die untere Begrenzung der Koppelzone des südlichen Untersuchungsgebietes in geringerer Tiefe als die des nördlichen Gebietes. Die Ergebnisse der thermischen 2D Modelle werden genutzt, um die räumliche Ausdehnung der Koppelzone in den mechanischen 3D Modellen festzulegen. Numerische 3D Simulationen werden genutzt, um zu verstehen, wie Geometrien, Rheologien und mechanische Parameter die Verformungspartitionierung und das Deformationsmuster im chilenischen Fore-Arc beeinflussen. Ich habe den Einfluss unterschiedlicher rheologischer Ansätze und unterschiedlicher physikalischer Eigenschaften auf den Fore-Arc untersucht, um Parameter zu identifizieren und zu bestimmen, die den Unterschied des Deformationsmusters zwischen dem nördlichen und südlichen Gebiet steuern. Die Ergebnisse der numerischen Studien stellen heraus, dass ein kleinerer Abtauchwinkel der ozeanischen Platte, ein hoher basaler Reibungskoeffizient, eine hohe Konvergenzschiefe und ein großer Elastizitätsmodul die Lokalisierung der Deformation im Fore-Arc Keil begünstigen. Basierend auf dieser Parameterstudie habe ich Modelle für die beiden chilenischen Gebiete ausgewählt, die in Beobachtungen erster Ordnung übereinstimmen. Diese ausgewählten Modelle erklären die unterschiedlichen Deformationsmuster durch eine größere Konvergenzschiefe, einen kleineren Abtauchwinkel der ozeanischen Platte und ein härteres Keilmaterial für das südliche Untersuchungsgebiet. Der Unterschied bezüglich der Deformationsmuster kann sogar größer sein, wenn ich eine größere Reibung zwischen den Platten im südlichen Gebiet anwende; jedoch zeigen einige unabhängige Beobachtungen eine umgekehrte Tendenz: eine Abnahme der Intensität der Koppelung von Nord nach Süd. Das schwächere Keilmaterial des ausgewählten Modells für das nördliche Untersuchungsgebiet steht im Zusammenhang mit der fortgeschrittenen Entwicklung des angrenzenden Orogens, der zentralen Anden. Die Analyse weltweiter Beispiele von schiefen Subduktionzonen unterstützt die Schlussfolgerung, dass ältere Subduktionzonen weniger ausgeprägte Lokalisierung von Strike-Slip Bewegung aufzeigen.

Da geologische Störungen können als Grundwasserleiter, -Barrieren oder als gemischte leitende /stauende Fluidsysteme wirken. Aufgrund dessen können Störungen maßgeblich den Grundwasserfluss im Untergrund beeinflussen, welcher deutliche Veränderungen des tiefen thermischen Feldes bewirken kann. Grundwasserdynamik und Temperaturveränderungen sind wiederum entscheidende Faktoren für die Exploration geothermischer Energie. Diese Studie untersuchte den Einfluss von Störungen auf das Fluidsystem und das thermische Feld im Untergrund. Sie erforschte die physikalischen Prozesse, welche das Fluidverhalten und die Temperaturverteilung in Störungen und in den umgebenden Gesteinen. Dazu wurden 3D Finite Elemente Simulationen des gekoppelten Fluid und Wärmetransports für synthetische sowie reale Modelszenarien auf unterschiedlichen Skalen durchgeführt. Um den Einfluss einer schräg einfallenden Störung systematisch durch die schrittweise Veränderung der hydraulischen Öffnungsweite und der Permeabilität, zu untersuchen, wurde ein klein-skaliges synthetisches Modell entwickelt. Ein inverser linearer Zusammenhang wurde festgestellt, welcher zeigt, dass sich die Fluidgeschwindigkeit in der Störung jeweils um ~1e-01 m/s verringert, wenn die Öffnungsweite der Störung um jeweils eine Magnitude vergrößert wird. Ein hoher Permeabilitätskontrast zwischen Störung und umgebender Matrix begünstigt die Fluidadvektion hin zur Störung und führt zu ausgeprägten Druck- und Temperaturveränderungen innerhalb und um die Störung herum. Bei geringem Permeabilitätskontrast zwischen Störung und umgebendem Gestein findet hingegen kein Fluidfluss in der Störung statt, wobei das hydrostatische Druck- sowie das Temperaturfeld unverändert bleiben. Auf Grundlage der synthetischen Modellierungsergebnisse wurde der Einfluss von Störungen auf einer größeren Skala anhand eines komplexeren (realen) geologischen Systems analysiert. Dabei handelt es sich um ein 3D Modell des Geothermiestandortes Groß Schönebeck, der ca. 40 km nördlich von Berlin liegt. Die Integration von einer permeablen und drei impermeablen Hauptstörungen, zeigte unterschiedlich starke Einflüsse auf Fluidzirkulation, Temperatur – und Druckfeld. Die modellierte konvektive Zirkulation in der permeablen Störung verändert das thermische Feld stark (bis zu 15 K). In den gering durchlässigen Störungen wird die Wärme ausschließlich durch Diffusion geleitet. Der konduktive Wärmetransport beeinflusst das thermische Feld nicht, bewirkt jedoch lokale Veränderungen des hydrostatischen Druckfeldes. Um den Einfluss großer Störungszonen mit kilometerweitem vertikalen Versatz auf das geothermische Feld der Beckenskala zu untersuchen, wurden gekoppelte Fluid- und Wärmetransportsimulationen für ein 3D Strukturmodell des Gebietes Brandenburg durchgeführt (Noack et al. 2010; 2013). Bezüglich der Störungspermeabilität wurden verschiedene geologische Szenarien modelliert, von denen zwei Endgliedermodelle ausgewertet wurden. Die Ergebnisse zeigten, dass die undurchlässigen Störungen den Fluidfluss nur lokal beeinflussen. Da sie als hydraulische Barrieren wirken, wird der Fluidfluss mir sehr geringen Geschwindigkeiten entlang der Störungen innerhalb eines Bereichs von ~ 1 km auf jeder Seite umgelenkt. Die modellierten lokalen Veränderungen des Grundwasserzirkulationssystems haben keinen beobachtbaren Effekt auf das Temperaturfeld. Hingegen erzeugen permeable Störungszonen eine ausgeprägte thermische Signatur innerhalb eines Einflussbereichs von ~ 2.4-8.8 km in -1000 m Tiefe und ~6-12 km in -3000 m Tiefe. Diese thermische Signatur, in der sich kältere und wärmere Temperaturbereiche abwechseln, wird durch auf- und abwärts gerichteten Fluidfluss innerhalb der Störung verursacht, der grundsätzlich durch existierende Gradienten in der hydraulischen Druckhöhe angetrieben wird. Alle Studien haben gezeigt, dass Störungen einen beachtlichen Einfluss auf den Fluid-, und Wärmefluss haben. Es stellte sich heraus, dass die Permeabilität in der Störung und in den umgebenden geologischen Schichten so wie der spezifische geologische Rahmen entscheidende Faktoren in der Ausbildung verschiedener Wärmetransportmechanismen sind, die sich in Störungen entwickeln können. Die von permeablen Störungen verursachten Temperaturveränderungen können lokal, jedoch groß sein, genauso wie die durch hydraulisch leitende und nichtleitende Störungen hervorgerufenen Veränderungen des Fluidystems. Letztlich haben die Simulationen für die unterschiedlich skalierten Modelle gezeigt, dass die Ergebnisse sich nicht aufeinander übertragen lassen und dass es notwendig ist, jeden geologischen Rahmen hinsichtlich Konfiguration und Größenskala gesondert zu betrachten. Abschließend hat diese Studie demonstriert, dass die Betrachtung von Störungen in 3D Finiten Elementen Modellen für die Simulation von gekoppeltem Fluid- und Wärmetransport auf unterschiedlichen Skalen möglich ist. Da diese Art von numerischen Simulationen sowohl die geologische Struktur des Untergrunds sowie die im Erdinnern ablaufenden physikalischen Prozesse integriert, können sie einen wertvollen Beitrag leisten, indem sie Feld- und Laborgestützte Untersuchungen vervollständigen. / Faults can act as conduits, barriers or mixed conduit/barrier systems to fluid flow. Therefore, faults may significantly influence fluid flow regimes operating in the subsurface, possibly resulting in distinct variations of the deep thermal field. Both, flow dynamics and temperature changes are in turn crucial factors that need to be taken into account for geothermal energy exploration. This study investigated the influence of faults on the subsurface fluid system and thermal field and explored the processes controlling fluid behavior and thermal distribution both within host rocks and faults. For this purpose, 3D finite element simulations of coupled fluid and heat transport have been carried out, both for synthetic and real-case model scenarios on different scales. A small-scale synthetic model was developed to systematically assess the impact of an inclined fault by changing gradually its hydraulic width and its permeability within the simulations. An observed linear inverse relationship revealed that changing the fault width by one order of magnitude results in a fluid velocity decrease (~1e-01 m/s) within the fault. A high permeability contrast between fault and matrix favors fluid advection into the fault and leads to pronounced pressure and temperature changes in and around the same domain. When the permeability contrast between fault domain and host rock is low, however, no fluid flow is observed in the fault, thus resulting in undisturbed hydrostatic pressure and temperature fields. On the basis of these synthetic fault modelling results, the influence of faults on a larger scale have been analyzed within a more complex (real-case) geological setting,- a 3D model of the geothermal site Groß Schönebeck / located ~40 km north of Berlin. The integration of one permeable and three impermeable major faults, resulted in distinct changes observed in the local fluid circulation, thermal and pressure field. Modelled convective circulation within the permeable fault decisively modifies the thermal field (up to 15 K). Within the low permeable faults, heat is transferred only by conduction, inducing no thermal imprint but local deviations of the hydrostatic pressure field. To investigate the impact of major fault zones on the basin-scale geothermal field, coupled fluid and heat transport simulations have been conducted for a 3D structural model for Brandenburg region (Noack et al. 2010; 2013). Different geological scenarios in terms of modelled fault permeability have been carried out of which two end member models are analyzed. The results showed that tight fault zones affect the flow field locally. Acting as hydraulic barriers, fluid flow is deviated with very low velocities along them within a range of ~ 1 km on either sides. The modelled local changes in the groundwater circulation system have no considerable effect on the temperature field. By contrast, permeable fault zones induce a pronounced signature on the thermal field extending over a distance of ~ 2.4-8.8 km at -1000 m depth and ~6-12 km at -3000 m depth. This thermal signature, characterized by alternating cooler and hotter temperature domains, is controlled by up- and downward directed flow within the fault domain, principally driven by existing hydraulic head gradients. All studies demonstrated that faults have a considerable impact on the fluid and heat flow. The permeability in faults and surrounding geological layers as well as the specific geological setting turned out to be crucial factors in controlling the different kinds of heat transfer mechanisms that may evolve in faults. Temperature variations caused by permeable faults may be local but significant as well as changes in fluid dynamics by both conduits and barriers. Thus, the results demonstrated the importance to consider faults in geothermal energy exploration. In the final analysis, the simulations for the small-, regional- and basin-scale models showed that the outcomes cannot be transferred by upscaling and that it is necessary to consider each geological setting separately with respect to its configuration and scale dimension. In summary, this study demonstrated that the consideration of faults in 3D finite element models for coupled fluid and heat transport simulations on different scales is feasible. As these type of numerical simulations integrate both, the structural setting of the subsurface and the physical processes controlling subsurface transport, the outcomes of this thesis may provide positive contributions in that they valuably complement field- and laboratory-based investigations.

Karst aquifers can be conceptualized as dual flow systems comprised of a low-conductive matrix with embedded high-conductive conduits / preferential flow zones. Discharge in conduits ranges from low-velocity laminar flow to high-velocity transitional and turbulent flow. Commonly employed continuum models do not account for the specific behavior of transitional and turbulent flow. In response to this limitation, enhancements have been made to MODFLOW, a commonly used groundwater flow model, by adding a discrete conduit network to the matrix continuum (hybrid model). The Conduit Flow Process (CFP) package is the latest realization of this model approach. CFP Mode 1 (CFPM1) computes laminar and turbulent flow in discrete conduits that are coupled to the laminar continuum model. CFP Mode 2 (CFPM2) accounts for turbulent flow in preferential flow layers by adapting the continuum model. Therefore, laminar hydraulic conduc-tivities are converted into turbulent hydraulic conductivities. CFPM2 was further modified to consider steady turbulent pipe flow. Karst models based on CFPM2 require potentially less input data and computational efforts than karst models based on CFPM1. Furthermore, CFPM2 integrates more easily into MODFLOW versions including e.g. transport models. Parameter studies for a synthetic catchment demonstrates that continuum models with turbulent flow representation and an additional flow barrier between conduits and matrix can represent karst systems similar to hybrid models. For simulation of highly transient flow processes in karst conduit systems, i.e. during flood events, it is crucial to consider dynamics such as free-surface flow, wave propagation, and changes between pressurized and non-pressurized conduit flow. The coupled overland- and groundwater flow model MODBRANCH was therefore enhanced to consider unsteady and non-uniform flow processes in karst conduits. Flow in discrete conduits is simulated using the Saint-Venant-equations for free-surface flow. Contrary to overland flow, the cross sectional area of karst conduits is finite. Accordingly, both pressurized and non-pressurized flow may occur within conduits. To simulate pressurized flow, a hypothetical, narrow, open-top slot (Preissmann slot) is added to the conduit crown, which allows the use of the free-surface flow equations for fully filled conduits. Beyond this, the model features a variable time step to consider wave speed variations, for example due to the transition from free-surface to pressurized flow. Parameter studies for a synthetic catchment demonstrate the significance of free-surface flow representation for variably filled conduits. / Karstgrundwasserleiter können als duale Fließsysteme konzeptionalisiert werden, bestehend aus einer geringdurchlässigen Matrix mit eingebundenen hochdurchlässigen Bereichen, z. B. Karströhren. Der Abfluss in den hochdurchlässigen Bereichen reicht von langsamer laminarer Strömung bis zu schneller turbulenter Strömung. Herkömmliche numerische Grundwasser-strömungsmodelle berücksichtigen nicht die spezifischen Eigenschaften von nicht-laminarer Strömung (Übergangsbereich laminar-turbulent bzw. turbulente Verhältnisse). Ein Ansatz um diese Einschränkung zu umgehen, ist die Erweiterung des laminaren Kontinuums um ein dis-kretes Röhrenmodell, das zustandsabhängig laminare und turbulente Strömung berücksichtigt (Hybridmodell). Eine aktuelle Umsetzung dieses Ansatzes ist Conduit Flow Process (CFP), ein Modul für das weitverbreitete Grundwasserströmungsmodell MODFLOW. CFP Mode 1 (CFPM1) berechnet laminare und turbulente Strömung in diskreten, mit dem Kontinuummodell gekoppelten Röhren. CFP Mode 2 (CFPM2) berücksichtigt nicht-laminare Strömung in hochdurchlässigen Schichten mit einer angepassten hydraulischen Leitfähigkeit des Kontinuummodells. CFPM2 wurde weiter modifiziert, so dass auch turbulente Strömung in Karströhren berechnet werden kann. Dadurch kann möglicherweise der Parameterbedarf sowie der Rechenaufwand gegenüber Hybrid¬modellen reduziert werden. CFPM2 lässt sich einfach in vorhandene MODFLOW Modelle einbinden, z. B. zur Berechnung von Transportprozessen. Parameterstudien für ein idealisiertes Karsteinzugsgebiet zeigen, dass Kontinuummodelle bei Berücksichtigung der turbulenten Strömung sowie des zusätzlichen hydraulischen Widerstand zwischen Röhren und Matrix, Karstsysteme ähnlich wie Hybridmodelle darstellen. Zur Simulation von instationären Prozessen in Karströhren, z. B. ausgeprägte Abflusssignale infolge pulsförmiger Grundwasserneubildung, ist es notwendig, dynamische Prozesse infolge Freispiegelabfluss, Wellenausbreitung sowie Wechsel zwischen Abfluss in teil- und vollgefüllten Röhren zu berücksichtigen. Aus diesem Grund wurde das numerische Modell MODBRANCH, welches ein diskretes Oberflächenwassermodell mit einem Kontinuummodell koppelt, so angepasst, dass instationäre und nichtgleichförmige Abflussprozesse in Karströhren berücksichtigt werden können. Der Abfluss in diskreten Röhren wird dabei mit den Saint-Venant-Gleichungen für Freispiegelabfluss berechnet. Im Gegensatz zu Oberflächengewässern ist der für den Abfluss zur Verfügung stehende Querschnitt in Karströhren limitiert, so dass sowohl Freispiegel- als auch Druckabfluss innerhalb der Röhren auftreten kann. Druckabfluss wird mit Hilfe eines schmalen virtuellen Schlitzes an der Röhrenoberkante simuliert (Preissmann Schlitz), der auch im Fall vollgefüllter Röhren die Anwendung der Gleichungen für Freispiegelabfluss erlaubt. Durch die Verwendung eines variablen Zeitschrittes kann die geänderte Dynamik beim Übergang von Freispiegel- zu Druckabfluss berücksichtigt werden. Parameterstudien für idealisierte, synthetische Karsteinzugsgebiete demonstrieren die Bedeutung der Berücksichtigung von Freispiegelabfluss in teilgefüllter Röhren.

Trocadores de Calor Solo-Ar (TCSA) são dispositivos usados para melhorar o condicionamento térmico de ambientes construídos. Eles funcionam através da ventilação do ar por um ou mais dutos enterrados, utilizando o solo como fonte ou sumidouro de calor. Em virtude da defasagem entre as temperaturas do ar e das camadas superficiais da Terra, é possível resfriar o ar no verão e aquecê-lo no inverno. Seus princípios de operação baseiam-se na mecânica dos fluidos e transferência de calor, áreas onde a teoria Construtal tem sido usada para melhorar o desempenho, ou reduzir imperfeições, de vários sistemas térmicos, após uma estratégica análise de suas geometrias (ou de seu design), também chamada de método Design Construtal que se fundamenta na lei construtal. Dessa lei, projetos em engenharia devem começar, por exemplo, descobrindo as arquiteturas que facilitam o escoamento entre uma fonte pontual e um volume, ou vice-versa. Explorando esses conceitos para TCSA, onde a corrente é o calor que flui entre os dutos e o solo, esta tese objetiva centralmente avaliar o desempenho térmico desses dispositivos em função de possíveis desenhos de suas estruturas. Particularmente, isso foi feito: (1) usando um duto e variando seu diâmetro e vazão de ar; (2) inserindo novos dutos, até um total de cinco, mantendo a vazão de ar constante, e estudando diversas disposições geométricas. Numericamente, foram empregados dois modelos testados e validados a partir de dados experimentais, que foram simulados no código comercial de dinâmica dos fluidos computacional, FLUENT, o qual baseia-se em métodos de volumes finitos. As geometrias e malhas foram geradas no software GAMBIT. Dos resultados numéricos, esta tese também avançou desenvolvendo novos modelos para o conceito de potencial térmico instantâneo, que é uma medida das diferenças de temperatura entre a saída e a entrada dos dutos de TCSA. Descobriu-se que estes modelos podiam ser simplificados em termos de funções senoidais, facilitando análises e definições analíticas subsequentes, criando um quadro resumo para grandes volumes de dados simulados. Das avaliações dos resultados nos estudos com um duto, foram encontradas faixas para a obtenção de rendimentos térmicos anuais estimados em pelo menos 80%, com mudanças no diâmetro e/ou vazão do ar. Analisando múltiplos dutos, foram reveladas melhorias significativas de rendimento com: (1) o aumento no número de dutos; (2) a redução, até certos limites, da razão entre os espaçamentos verticais e horizontais entre eles; (3) a ampliação da razão entre o volume da instalação e o volume do domínio computacional. Destaca-se, ainda, que em todos os casos foram calculados balanços energéticos, com a estimativa das taxas (e quantidade) de calor trocado, bem como da energia elétrica consumida pelos ventiladores, mostrando sua viabilidade no condicionamento de ar com redução do consumo de energia elétrica. / The Earth-Air Heat Exchangers (EAHE) are devices used to improve the thermal conditions of built environments. They work by blowing the air inside buried ducts to use the soil as heat source or sink. Due to a phase difference between the air and ground temperatures, it is possible to cool the air in the summer and heat it in the winter. These operating principles are based on the areas of fluid mechanics and heat transfer, where the Constructal theory has been used to improve the performance, or reduce imperfections, of several thermal systems, after a strategic evaluation of their geometries (or their design), also called Constructal Design method, which is anchored by the constructal law. From this law, engineering projects should begin, for example, finding the architectures which ease the flow between a source point and a volume, or vice versa. Bringing these concepts to EAHE, where the current is the heat which flows between the ducts and the soil, this thesis mainly aims to assess the thermal performance of EAHE due to possible designs of their layouts. Particularly, this was done by: (1) adopting one duct and varying its diameter and/or the air flow; (2) using up to five ducts, keeping the air flow constant, and studying various geometric shapes. Numerically, two models were used, which have been tested and validated from experimental data, and the simulations were done in the CFD software FLUENT, which is based in the finite volume methods. The geometries and meshes were generated by the software GAMBIT. From the numerical results, this thesis also pursued developing new models for the concept of the instantaneous thermal potential, which is a measure of the temperatures differences between the ducts outlet and inlet. It was found that these models can be simplified in terms of sinusoidal functions, helping analysis, definitions, and creating a framework to summarize large volumes of simulated data. From the evaluation of the data in the studies with one duct, ranges for changes in their diameter and/or air flow were found to obtain at least 80% of estimated annual efficiency. Analyzing multiple ducts, significant improvements in efficiency were also obtained by: (1) increasing the number of ducts; (2) reducing, to some extent, the ratio between their vertical and horizontal spacings; (3) increasing the ratio between the installation volume and the computational domain. It should also be highlighted that in all cases the energy balances were computed, to find the rates (and amounts) of exchanged heat, as well as the electric energy consumed by the fans, showing the economic viability of using EAHE as devices for air conditioning.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / A necessidade de se obter soluções para os diversos problemas estruturais na engenharia, associado ao amplo crescimento e difusão da construção metálica fez surgir, por volta dos anos 60, os perfis tubulares. As ligações soldadas entre perfis tubulares são as ligações de maior difusão nas construções com este tipo de perfil, sendo necessário portanto, estudos mais detalhados que tornem a utilização de ligações soldadas mais seguras e otimizadas. O Eurocode 3 possui um item específico para as ligações entre perfis tubulares, em especial as ligações soldadas devido a suas particularidades. Infelizmente, a norma brasileira NBR8800 não aborda as ligações envolvendo perfis tubulares. Para a execução da análise numérica no presente trabalho, faz-se necessário uma revisão bibliográfica. Com estes resultados, foram feitas modelagens de ligações tubulares com o método dos elementos finitos de forma a otimizar os modelos a serem utilizados em uma análise paramétrica futura. Desta forma, foram desenvolvidos dois modelos numéricos, um considerando ligação tipo T entre perfis tubulares quadrados e outro para uma ligação tipo K entre perfis circulares. Estes modelos foram caracterizados através de elementos de casca com seis graus de liberdade por nó considerando-se adicionalmente o efeito de membrana. A análise não-linear realizada considerou a não-linearidade do material através do critério de plastificação de Von Mises através de uma lei constitutiva tensão versus deformação bi-linear e a não-linearidade geométrica foi atribuída através da formulação de Lagrange atualizado. Dentre as principais conclusões obtidas no presente trabalho, pode-se citar que os resultados para as ligações tipo T, o Eurocode 3, fornece resultados que precisam ser observados com cautela. Todavia, para as ligações do tipo K, os resultados numéricos mostraram-se sempre inferiores aos valores através do Eurocode 3, representando um dimensionamento a favor da segurança. / Tubular structures were conceived in the sixties, to solve several structural engineering problems and were boosted by the fast and continuous growth of the use of steel in the construction industry. The most popular tubular joints are welded. These specific connections still require further investigations to enable a complete understanding of their structural performance aiming to develop optimum and reliable solutions. The Eurocode 3 provides specific guidance for tubular connections, specially the welded joints, due to their particular design peculiarities. Unfortunately, the Brazilian Steel Design Standard, NBR8800, does not provide any recommendation for tubular members joint design. The numerical analysis performed in the current investigation required an updated and comprehensive literature review of tubular joint structural design and response. With this review in hand, numerical models were conceived and optimised with the aid of the Finite Element Method and were used in subsequent parametrical analysis made to enhance their structural response knwolledge. Two numerical models were developed with hollow welded joints of: T square section members and K circular section members. All the models adopted six dof. shell finite elements including the membrane action. The material response was non-linear using a bilinear stress vs. strain curve and a Von Mises yield criterion. The geometric non-linearity was considered through the use of an Updated Lagragean formulation. The most significant conclusions of the present investigation were related to the analyzed T joints results that indicated that the Eurocode 3 provisions should be adopted with caution. On the other hand, for the K joints, the numerical results were always lower than the Eurocode 3 results, and could be considered a lower bound for the investigated problem.

Neste trabalho, modela-se a ocorrência de contaminação de um aquífero hipotético por derramamento de substância tóxica e analisa-se a solução de descontaminação baseada na retirada do contaminante através de bombeamento feito por poços de extração. O projeto do sistema de remediação envolve a escolha do número de poços a serem instalados, suas localizações e vazões de modo a maximizar a quantidade de poluente extraída e ao mesmo tempo minimizar o custo total do sistema. A busca da solução ótima é feita de forma automática, através de um sistema paralelo de otimização por simulação numérica, composto por três subsistemas: (a) simulador numérico - resolve numericamente o modelo matemático do aquífero contaminado; (b) otimizador automático - implementa o método dos algoritmos genéticos para busca das localizações e vazões ótimas dos poços de extração; (c)sistema computacional distribuído - gerencia a distribuição e a execução paralela das simulações numéricas. Foram feitos experimentos em vários ambientes computacionais: homogêneo, heterogêneo, em grande escala, usando máquinas não dedicadas, interligadas por rede local e ambiente de grade, e seus resultados demonstram a aplicabilidade da metodologia. / In this dissertation, a hypothetical aquifer that has been contaminated by the dumping of toxic substances is modeled. The remediation strategy considered is based on withdrawal, which requires the removal of contaminated groundwater from the aquifer by pumping. The design of such a system involves the choice of the number of extracting wells to be installed, their locations and pumping rates,with the goal of maximizing the amount of contaminant extracted, while minimizing the cost of the system. To find the optimal solution, a numerical simulation optimization parallel system is used, which is composed by three subsystems: (a) numerical simulator - numerically solves the mathematical model ofthe contaminated aquifer; (b) optimizer - implements the genetic algorithm method to search for optimal locations and pumping rates for the extracting wells; (c)distributed computing system - manages the distribuition and parallel execution of the numerical simulations. Experiments were done in many different computational environments: homogeneous, heterogeneous, in large scale, using non dedicated computers, connected via local network, and computational grids, and their results demonstrate the methodologys applicability.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / A necessidade de se obter soluções para os diversos problemas estruturais na engenharia, associado ao amplo crescimento e difusão da construção metálica fez surgir, por volta dos anos 60, os perfis tubulares. As ligações soldadas entre perfis tubulares são as ligações de maior difusão nas construções com este tipo de perfil, sendo necessário portanto, estudos mais detalhados que tornem a utilização de ligações soldadas mais seguras e otimizadas. O Eurocode 3 possui um item específico para as ligações entre perfis tubulares, em especial as ligações soldadas devido a suas particularidades. Infelizmente, a norma brasileira NBR8800 não aborda as ligações envolvendo perfis tubulares. Para a execução da análise numérica no presente trabalho, faz-se necessário uma revisão bibliográfica. Com estes resultados, foram feitas modelagens de ligações tubulares com o método dos elementos finitos de forma a otimizar os modelos a serem utilizados em uma análise paramétrica futura. Desta forma, foram desenvolvidos dois modelos numéricos, um considerando ligação tipo T entre perfis tubulares quadrados e outro para uma ligação tipo K entre perfis circulares. Estes modelos foram caracterizados através de elementos de casca com seis graus de liberdade por nó considerando-se adicionalmente o efeito de membrana. A análise não-linear realizada considerou a não-linearidade do material através do critério de plastificação de Von Mises através de uma lei constitutiva tensão versus deformação bi-linear e a não-linearidade geométrica foi atribuída através da formulação de Lagrange atualizado. Dentre as principais conclusões obtidas no presente trabalho, pode-se citar que os resultados para as ligações tipo T, o Eurocode 3, fornece resultados que precisam ser observados com cautela. Todavia, para as ligações do tipo K, os resultados numéricos mostraram-se sempre inferiores aos valores através do Eurocode 3, representando um dimensionamento a favor da segurança. / Tubular structures were conceived in the sixties, to solve several structural engineering problems and were boosted by the fast and continuous growth of the use of steel in the construction industry. The most popular tubular joints are welded. These specific connections still require further investigations to enable a complete understanding of their structural performance aiming to develop optimum and reliable solutions. The Eurocode 3 provides specific guidance for tubular connections, specially the welded joints, due to their particular design peculiarities. Unfortunately, the Brazilian Steel Design Standard, NBR8800, does not provide any recommendation for tubular members joint design. The numerical analysis performed in the current investigation required an updated and comprehensive literature review of tubular joint structural design and response. With this review in hand, numerical models were conceived and optimised with the aid of the Finite Element Method and were used in subsequent parametrical analysis made to enhance their structural response knwolledge. Two numerical models were developed with hollow welded joints of: T square section members and K circular section members. All the models adopted six dof. shell finite elements including the membrane action. The material response was non-linear using a bilinear stress vs. strain curve and a Von Mises yield criterion. The geometric non-linearity was considered through the use of an Updated Lagragean formulation. The most significant conclusions of the present investigation were related to the analyzed T joints results that indicated that the Eurocode 3 provisions should be adopted with caution. On the other hand, for the K joints, the numerical results were always lower than the Eurocode 3 results, and could be considered a lower bound for the investigated problem.

Trocadores de Calor Solo-Ar (TCSA) são dispositivos usados para melhorar o condicionamento térmico de ambientes construídos. Eles funcionam através da ventilação do ar por um ou mais dutos enterrados, utilizando o solo como fonte ou sumidouro de calor. Em virtude da defasagem entre as temperaturas do ar e das camadas superficiais da Terra, é possível resfriar o ar no verão e aquecê-lo no inverno. Seus princípios de operação baseiam-se na mecânica dos fluidos e transferência de calor, áreas onde a teoria Construtal tem sido usada para melhorar o desempenho, ou reduzir imperfeições, de vários sistemas térmicos, após uma estratégica análise de suas geometrias (ou de seu design), também chamada de método Design Construtal que se fundamenta na lei construtal. Dessa lei, projetos em engenharia devem começar, por exemplo, descobrindo as arquiteturas que facilitam o escoamento entre uma fonte pontual e um volume, ou vice-versa. Explorando esses conceitos para TCSA, onde a corrente é o calor que flui entre os dutos e o solo, esta tese objetiva centralmente avaliar o desempenho térmico desses dispositivos em função de possíveis desenhos de suas estruturas. Particularmente, isso foi feito: (1) usando um duto e variando seu diâmetro e vazão de ar; (2) inserindo novos dutos, até um total de cinco, mantendo a vazão de ar constante, e estudando diversas disposições geométricas. Numericamente, foram empregados dois modelos testados e validados a partir de dados experimentais, que foram simulados no código comercial de dinâmica dos fluidos computacional, FLUENT, o qual baseia-se em métodos de volumes finitos. As geometrias e malhas foram geradas no software GAMBIT. Dos resultados numéricos, esta tese também avançou desenvolvendo novos modelos para o conceito de potencial térmico instantâneo, que é uma medida das diferenças de temperatura entre a saída e a entrada dos dutos de TCSA. Descobriu-se que estes modelos podiam ser simplificados em termos de funções senoidais, facilitando análises e definições analíticas subsequentes, criando um quadro resumo para grandes volumes de dados simulados. Das avaliações dos resultados nos estudos com um duto, foram encontradas faixas para a obtenção de rendimentos térmicos anuais estimados em pelo menos 80%, com mudanças no diâmetro e/ou vazão do ar. Analisando múltiplos dutos, foram reveladas melhorias significativas de rendimento com: (1) o aumento no número de dutos; (2) a redução, até certos limites, da razão entre os espaçamentos verticais e horizontais entre eles; (3) a ampliação da razão entre o volume da instalação e o volume do domínio computacional. Destaca-se, ainda, que em todos os casos foram calculados balanços energéticos, com a estimativa das taxas (e quantidade) de calor trocado, bem como da energia elétrica consumida pelos ventiladores, mostrando sua viabilidade no condicionamento de ar com redução do consumo de energia elétrica. / The Earth-Air Heat Exchangers (EAHE) are devices used to improve the thermal conditions of built environments. They work by blowing the air inside buried ducts to use the soil as heat source or sink. Due to a phase difference between the air and ground temperatures, it is possible to cool the air in the summer and heat it in the winter. These operating principles are based on the areas of fluid mechanics and heat transfer, where the Constructal theory has been used to improve the performance, or reduce imperfections, of several thermal systems, after a strategic evaluation of their geometries (or their design), also called Constructal Design method, which is anchored by the constructal law. From this law, engineering projects should begin, for example, finding the architectures which ease the flow between a source point and a volume, or vice versa. Bringing these concepts to EAHE, where the current is the heat which flows between the ducts and the soil, this thesis mainly aims to assess the thermal performance of EAHE due to possible designs of their layouts. Particularly, this was done by: (1) adopting one duct and varying its diameter and/or the air flow; (2) using up to five ducts, keeping the air flow constant, and studying various geometric shapes. Numerically, two models were used, which have been tested and validated from experimental data, and the simulations were done in the CFD software FLUENT, which is based in the finite volume methods. The geometries and meshes were generated by the software GAMBIT. From the numerical results, this thesis also pursued developing new models for the concept of the instantaneous thermal potential, which is a measure of the temperatures differences between the ducts outlet and inlet. It was found that these models can be simplified in terms of sinusoidal functions, helping analysis, definitions, and creating a framework to summarize large volumes of simulated data. From the evaluation of the data in the studies with one duct, ranges for changes in their diameter and/or air flow were found to obtain at least 80% of estimated annual efficiency. Analyzing multiple ducts, significant improvements in efficiency were also obtained by: (1) increasing the number of ducts; (2) reducing, to some extent, the ratio between their vertical and horizontal spacings; (3) increasing the ratio between the installation volume and the computational domain. It should also be highlighted that in all cases the energy balances were computed, to find the rates (and amounts) of exchanged heat, as well as the electric energy consumed by the fans, showing the economic viability of using EAHE as devices for air conditioning.

Orientador: Maria Cecília Amorim Teixeira da Silva / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Civil, Arquitetura e Urbanismo / Made available in DSpace on 2018-08-17T07:55:36Z (GMT). No. of bitstreams: 1 Pires_SusanadeLima_M.pdf: 919425 bytes, checksum: d13d24300a7bb026a2c9352fd83b1c35 (MD5) Previous issue date: 2006 / Resumo: Este trabalho apresenta um modelo numérico para o estudo de pilares isolados de concreto armado submetidos à flexão normal composta onde as não linearidades física e geométrica são levadas em conta de maneira rigorosa. O procedimento numérico desenvolvido para o cálculo dos deslocamentos é baseado no Método dos Elementos Finitos. Para a inclusão das não linearidades física e geométrica é utilizado o Método de Newton-Raphson Modificado que é um procedimento iterativo para a determinação de raízes de funções não lineares. A relação tensão-deformação do concreto é representada pelo diagrama parábola-retângulo. O aço é considerado um material elastoplástico perfeito. É admitida perfeita aderência entre o concreto e o aço. São desenvolvidos algoritmos para o dimensionamento de pilares e para o estudo da estabilidade do pilar. É desenvolvido um programa de computador utilizando os algoritmos mencionados acima. A eficiência do programa é verificada por meio da comparação com trabalhos experimentais e com o programa comercial CAD/TQS. Os resultados do dimensionamento feito com o programa proposto apresentaram boa concordância com os dados experimentais selecionados. Finalmente é feito um estudo comparativo entre os resultados do dimensionamento de pilares feito pelo programa elaborado, com os resultados obtidos pelas simplificações admitidas pela Norma Brasileira NBR 6118/2003. Esta comparação mostra que, para os exemplos apresentados, a NBR 6118/2003 superestima o dimensionamento de pilares esbeltos realizados através dos métodos simplificados. Para pilares com índice de esbeltez superiores a 90, a norma está certa em impedir o dimensionamento com os métodos simplificados pois podem apresentar resultados contra a segurança. Os efeitos causados pela deformação lenta não foram incluídos neste trabalho / Abstract: This work presents a numerical model for the study of reinforced concrete columns subjected to combined axial load and bending where the geometric and material nonlinearities are rigorously considered. The numerical procedure developed to calculate the displacements is based on the Finite Element Method. An iterative procedure for the determination of roots of nonlinear functions, the Modified Newton-Raphson Method, is used for including nonlinearities. The parabolic-rectangular diagram is adopted as the concrete stress-strain relationship. The steel is considered a perfect elastoplastic material, and a perfect concrete-to-steel bond is admitted. Algorithms are developed for dimensioning and studying the stability of columns. A computer program is developed using the algorithms mentioned above. The efficiency of the program is verified against available experimental data and commercial programs such as CAD/TQS. The values obtained from the program have presented good agreement with experimental data. Finally a comparative study is made among the results of the designed columns with values obtained from Brazilian Code NBR 6118/2003 approaches. This comparison shows that, for the presented examples, NBR 6118/2003 overestimates the slender columns design when using the approaches, and for columns with slenderness index over 90 the code is correct in preventing the design through the simplified methods, for it can lead to unsafe results. The effects caused by the creep are not enclosed in this work / Mestrado / Estruturas / Mestre em Engenharia Civil

Karst aquifers can be conceptualized as dual flow systems comprised of a low-conductive matrix with embedded high-conductive conduits / preferential flow zones. Discharge in conduits ranges from low-velocity laminar flow to high-velocity transitional and turbulent flow. Commonly employed continuum models do not account for the specific behavior of transitional and turbulent flow. In response to this limitation, enhancements have been made to MODFLOW, a commonly used groundwater flow model, by adding a discrete conduit network to the matrix continuum (hybrid model). The Conduit Flow Process (CFP) package is the latest realization of this model approach. CFP Mode 1 (CFPM1) computes laminar and turbulent flow in discrete conduits that are coupled to the laminar continuum model. CFP Mode 2 (CFPM2) accounts for turbulent flow in preferential flow layers by adapting the continuum model. Therefore, laminar hydraulic conduc-tivities are converted into turbulent hydraulic conductivities. CFPM2 was further modified to consider steady turbulent pipe flow. Karst models based on CFPM2 require potentially less input data and computational efforts than karst models based on CFPM1. Furthermore, CFPM2 integrates more easily into MODFLOW versions including e.g. transport models. Parameter studies for a synthetic catchment demonstrates that continuum models with turbulent flow representation and an additional flow barrier between conduits and matrix can represent karst systems similar to hybrid models. For simulation of highly transient flow processes in karst conduit systems, i.e. during flood events, it is crucial to consider dynamics such as free-surface flow, wave propagation, and changes between pressurized and non-pressurized conduit flow. The coupled overland- and groundwater flow model MODBRANCH was therefore enhanced to consider unsteady and non-uniform flow processes in karst conduits. Flow in discrete conduits is simulated using the Saint-Venant-equations for free-surface flow. Contrary to overland flow, the cross sectional area of karst conduits is finite. Accordingly, both pressurized and non-pressurized flow may occur within conduits. To simulate pressurized flow, a hypothetical, narrow, open-top slot (Preissmann slot) is added to the conduit crown, which allows the use of the free-surface flow equations for fully filled conduits. Beyond this, the model features a variable time step to consider wave speed variations, for example due to the transition from free-surface to pressurized flow. Parameter studies for a synthetic catchment demonstrate the significance of free-surface flow representation for variably filled conduits. / Karstgrundwasserleiter können als duale Fließsysteme konzeptionalisiert werden, bestehend aus einer geringdurchlässigen Matrix mit eingebundenen hochdurchlässigen Bereichen, z. B. Karströhren. Der Abfluss in den hochdurchlässigen Bereichen reicht von langsamer laminarer Strömung bis zu schneller turbulenter Strömung. Herkömmliche numerische Grundwasser-strömungsmodelle berücksichtigen nicht die spezifischen Eigenschaften von nicht-laminarer Strömung (Übergangsbereich laminar-turbulent bzw. turbulente Verhältnisse). Ein Ansatz um diese Einschränkung zu umgehen, ist die Erweiterung des laminaren Kontinuums um ein dis-kretes Röhrenmodell, das zustandsabhängig laminare und turbulente Strömung berücksichtigt (Hybridmodell). Eine aktuelle Umsetzung dieses Ansatzes ist Conduit Flow Process (CFP), ein Modul für das weitverbreitete Grundwasserströmungsmodell MODFLOW. CFP Mode 1 (CFPM1) berechnet laminare und turbulente Strömung in diskreten, mit dem Kontinuummodell gekoppelten Röhren. CFP Mode 2 (CFPM2) berücksichtigt nicht-laminare Strömung in hochdurchlässigen Schichten mit einer angepassten hydraulischen Leitfähigkeit des Kontinuummodells. CFPM2 wurde weiter modifiziert, so dass auch turbulente Strömung in Karströhren berechnet werden kann. Dadurch kann möglicherweise der Parameterbedarf sowie der Rechenaufwand gegenüber Hybrid¬modellen reduziert werden. CFPM2 lässt sich einfach in vorhandene MODFLOW Modelle einbinden, z. B. zur Berechnung von Transportprozessen. Parameterstudien für ein idealisiertes Karsteinzugsgebiet zeigen, dass Kontinuummodelle bei Berücksichtigung der turbulenten Strömung sowie des zusätzlichen hydraulischen Widerstand zwischen Röhren und Matrix, Karstsysteme ähnlich wie Hybridmodelle darstellen. Zur Simulation von instationären Prozessen in Karströhren, z. B. ausgeprägte Abflusssignale infolge pulsförmiger Grundwasserneubildung, ist es notwendig, dynamische Prozesse infolge Freispiegelabfluss, Wellenausbreitung sowie Wechsel zwischen Abfluss in teil- und vollgefüllten Röhren zu berücksichtigen. Aus diesem Grund wurde das numerische Modell MODBRANCH, welches ein diskretes Oberflächenwassermodell mit einem Kontinuummodell koppelt, so angepasst, dass instationäre und nichtgleichförmige Abflussprozesse in Karströhren berücksichtigt werden können. Der Abfluss in diskreten Röhren wird dabei mit den Saint-Venant-Gleichungen für Freispiegelabfluss berechnet. Im Gegensatz zu Oberflächengewässern ist der für den Abfluss zur Verfügung stehende Querschnitt in Karströhren limitiert, so dass sowohl Freispiegel- als auch Druckabfluss innerhalb der Röhren auftreten kann. Druckabfluss wird mit Hilfe eines schmalen virtuellen Schlitzes an der Röhrenoberkante simuliert (Preissmann Schlitz), der auch im Fall vollgefüllter Röhren die Anwendung der Gleichungen für Freispiegelabfluss erlaubt. Durch die Verwendung eines variablen Zeitschrittes kann die geänderte Dynamik beim Übergang von Freispiegel- zu Druckabfluss berücksichtigt werden. Parameterstudien für idealisierte, synthetische Karsteinzugsgebiete demonstrieren die Bedeutung der Berücksichtigung von Freispiegelabfluss in teilgefüllter Röhren.