Pore-scale modeling of two-phase flow instabilities in porous media

Les problèmes d'écoulements multiphasiques en média poreux sont d'un grand intérêt pour de nombreuses applications scientifiques et techniques ; comme la séquestration de C02, l'extraction de pétrole et la dépollution des aquifères. La complexité intrinsèque des systèmes multiphasiques et l'hétérogénéité des formations géologiques sur des échelles multiples représentent un challenge majeur pour comprendre et modéliser les déplacements immiscibles dans les milieux poreux. Les descriptions à l'échelle supérieure basées sur la généralisation de l'équation de Darcy sont largement utilisées, mais ces méthodes sont sujettes à limitations pour les écoulements présentant de l'hystérèse. Les avancées récentes en terme de performances computationnelles et le développement de méthodes précises pour caractériser l'espace interstitiel ainsi que la distribution des phases ont favorisé l'utilisation de modèles qui permettent une résolution fine à l'échelle du pore. Ces modèles offrent un aperçu des caractéristiques de l'écoulement qui ne peuvent pas être facilement observées en laboratoire et peuvent être utilisé pour expliquer la différence entre les processus physiques et les modèles à l'échelle macroscopique existants. L'objet premier de la thèse se porte sur la simulation numérique directe : les équations de Navier-Stokes sont résolues dans l'espace interstitiel et la méthode du volume de fluide (VOF) est employée pour suivre l'évolution de l'interface. Dans VOF, la distribution des phases est décrite par une fonction fluide pour l'ensemble du domaine et des conditions aux bords particulières permettent la prise en compte des propriétés de mouillage du milieu poreux. Dans la première partie de la thèse, nous simulons le drainage dans une cellule Hele-Shaw 2D avec des obstacles cylindriques. Nous montrons que l'approche proposée est applicable même pour des ratios de densité et de viscosité très importants et permet de modéliser la transition entre déplacement stable et digitation visqueuse. Nous intéressons ensuite à l'interprétation de la pression capillaire à l'échelle macroscopique. Nous montrons que les techniques basées sur la moyenne spatiale de la pression présentent plusieurs limitations et sont imprécises en présence d'effets visqueux et de piégeage. Au contraire, une définition basée sur l'énergie permet de séparer les contributions capillaires des effets visqueux. La seconde partie de la thèse est consacrée à l'investigation des effets d'inertie associés aux reconfigurations irréversibles du ménisque causé par l'interface des instabilités. Comme prototype pour ces phénomènes, nous étudions d'abord la dynamique d'un ménisque dans un pore angulaire. Nous montrons que, dans un réseau de pores cubiques, les sauts et reconfigurations sont si fréquents que les effets d'inertie mènent à différentes configurations des fluides. A cause de la non-linéarité du problème, la distribution des fluides influence le travail des forces de pression, qui, à son tour, provoque une chute de pression dans la loi de Darcy. Cela suggère que ces phénomènes devraient être pris en compte lorsque que l'on décrit l'écoulement multiphasique en média poreux à l'échelle macroscopique. La dernière partie de la thèse s'attache à démontrer la validité de notre approche par une comparaison avec des expériences en laboratoire : un drainage instable dans un milieu poreux quasi 2D (une cellule Hele-Shaw avec des obstacles cylindriques). Plusieurs simulations sont tournées sous différentes conditions aux bords et en utilisant différents modèles (modèle intégré 2D et modèle 3D) afin de comparer certaines quantités macroscopiques avec les observations au laboratoire correspondantes. Malgré le challenge de modéliser des déplacements instables, où, par définition, de petites perturbations peuvent grandir sans fin, notre approche numérique apporte de résultats satisfaisants pour tous les cas étudiés. - Problems involving multiphase flow in porous media are of great interest in many scientific and engineering applications including Carbon Capture and Storage, oil recovery and groundwater remediation. The intrinsic complexity of multiphase systems and the multi scale heterogeneity of geological formations represent the major challenges to understand and model immiscible displacement in porous media. Upscaled descriptions based on generalization of Darcy's law are widely used, but they are subject to several limitations for flow that exhibit hysteric and history- dependent behaviors. Recent advances in high performance computing and the development of accurate methods to characterize pore space and phase distribution have fostered the use of models that allow sub-pore resolution. These models provide an insight on flow characteristics that cannot be easily achieved by laboratory experiments and can be used to explain the gap between physical processes and existing macro-scale models. We focus on direct numerical simulations: we solve the Navier-Stokes equations for mass and momentum conservation in the pore space and employ the Volume Of Fluid (VOF) method to track the evolution of the interface. In the VOF the distribution of the phases is described by a fluid function (whole-domain formulation) and special boundary conditions account for the wetting properties of the porous medium. In the first part of this thesis we simulate drainage in a 2-D Hele-Shaw cell filled with cylindrical obstacles. We show that the proposed approach can handle very large density and viscosity ratios and it is able to model the transition from stable displacement to viscous fingering. We then focus on the interpretation of the macroscopic capillary pressure showing that pressure average techniques are subject to several limitations and they are not accurate in presence of viscous effects and trapping. On the contrary an energy-based definition allows separating viscous and capillary contributions. In the second part of the thesis we investigate inertia effects associated with abrupt and irreversible reconfigurations of the menisci caused by interface instabilities. As a prototype of these phenomena we first consider the dynamics of a meniscus in an angular pore. We show that in a network of cubic pores, jumps and reconfigurations are so frequent that inertia effects lead to different fluid configurations. Due to the non-linearity of the problem, the distribution of the fluids influences the work done by pressure forces, which is in turn related to the pressure drop in Darcy's law. This suggests that these phenomena should be taken into account when upscaling multiphase flow in porous media. The last part of the thesis is devoted to proving the accuracy of the numerical approach by validation with experiments of unstable primary drainage in a quasi-2D porous medium (i.e., Hele-Shaw cell filled with cylindrical obstacles). We perform simulations under different boundary conditions and using different models (2-D integrated and full 3-D) and we compare several macroscopic quantities with the corresponding experiment. Despite the intrinsic challenges of modeling unstable displacement, where by definition small perturbations can grow without bounds, the numerical method gives satisfactory results for all the cases studied.

Development of a Subgrade Drainage Model for Unpaved Roads, TR-654, 2015

With over 68 thousand miles of gravel roads in Iowa and the importance of these roads within the farm-to-market transportation system, proper water management becomes critical for maintaining the integrity of the roadway materials. However, the build-up of water within the aggregate subbase can lead to frost boils and ultimately potholes forming at the road surface. The aggregate subbase and subgrade soils under these gravel roads are produced with material opportunistically chosen from local sources near the site and, many times, the compositions of these sublayers are far from ideal in terms of proper water drainage with the full effects of this shortcut not being fully understood. The primary objective of this project was to provide a physically-based model for evaluating the drainability of potential subbase and subgrade materials for gravel roads in Iowa. The Richards equation provided the appropriate framework to study the transient unsaturated flow that usually occurs through the subbase and subgrade of a gravel road. From which, we identified that the saturated hydraulic conductivity, Ks, was a key parameter driving the time to drain of subgrade soils found in Iowa, thus being a good proxy variable for accessing roadway drainability. Using Ks, derived from soil texture, we were able to identify potential problem areas in terms of roadway drainage . It was found that there is a threshold for Ks of 15 cm/day that determines if the roadway will drain efficiently, based on the requirement that the time to drain, Td, the surface roadway layer does not exceed a 2-hr limit. Two of the three highest abundant textures (loam and silty clay loam), which cover nearly 60% of the state of Iowa, were found to have average Td values greater than the 2-hr limit. With such a large percentage of the state at risk for the formation of boils due to the soil with relatively low saturated hydraulic conductivity values, it seems pertinent that we propose alternative design and/or maintenance practices to limit the expensive repair work in Iowa. The addition of drain tiles or French mattresses my help address drainage problems. However, before pursuing this recommendation, a comprehensive cost-benefit analysis is needed.

Re-Engineering CTE in Iowa: Aligning Career and Technical Education Reform, Workforce Development and Economic Development, August 2007

Iowa is a relatively small state and is on the rebound economically. It has an overall population that is stable, but which is shifting within the state from more rural areas to suburban and urban centers. There is a very tight labor market with high levels of employment. Iowa now has a time-sensitive opportunity to exert global leadership in renewable energy, while maintaining its leadership in other key industries like finance and agriculture.

Use of LIDAR-based Elevation Data for Highway Drainage Analysis: A Qualitative Assessment, 2003

US Geological Survey (USGS) based elevation data are the most commonly used data source for highway hydraulic analysis; however, due to the vertical accuracy of USGS-based elevation data, USGS data may be too “coarse” to adequately describe surface profiles of watershed areas or drainage patterns. Additionally hydraulic design requires delineation of much smaller drainage areas (watersheds) than other hydrologic applications, such as environmental, ecological, and water resource management. This research study investigated whether higher resolution LIDAR based surface models would provide better delineation of watersheds and drainage patterns as compared to surface models created from standard USGS-based elevation data. Differences in runoff values were the metric used to compare the data sets. The two data sets were compared for a pilot study area along the Iowa 1 corridor between Iowa City and Mount Vernon. Given the limited breadth of the analysis corridor, areas of particular emphasis were the location of drainage area boundaries and flow patterns parallel to and intersecting the road cross section. Traditional highway hydrology does not appear to be significantly impacted, or benefited, by the increased terrain detail that LIDAR provided for the study area. In fact, hydrologic outputs, such as streams and watersheds, may be too sensitive to the increased horizontal resolution and/or errors in the data set. However, a true comparison of LIDAR and USGS-based data sets of equal size and encompassing entire drainage areas could not be performed in this study. Differences may also result in areas with much steeper slopes or significant changes in terrain. LIDAR may provide possibly valuable detail in areas of modified terrain, such as roads. Better representations of channel and terrain detail in the vicinity of the roadway may be useful in modeling problem drainage areas and evaluating structural surety during and after significant storm events. Furthermore, LIDAR may be used to verify the intended/expected drainage patterns at newly constructed highways. LIDAR will likely provide the greatest benefit for highway projects in flood plains and areas with relatively flat terrain where slight changes in terrain may have a significant impact on drainage patterns.

07003-001 Clear Lake Watershed Final Report

Clear Lake, Iowa's third largest natural lake, is a premier natural resource and popular recreational destination in north central Iowa. Despite the lake's already strong recreational use, water quality concerns have not allowed the lake to reach its full potential. Clear Lake is listed on Iowa's 2004 303(d) Impaired Waters List due to excessive levels of phosphorus, bacteria, and turbidity. Urban storm water runoff from the 8,600 acre watershed is a significant contributor to Clear Lake's impairment. Local communities have been working towards the goal of making improvements at all 30 storm water outlets that have a drainage area of five acres or more and have a cost effective solution. Many improvements have already been made, and now there are only seven storm water outlet sites remaining that still need protection in order to meet the goal. The storm water improvements have been very effective in reducing contaminants in urban runoff, achieving reduction levels in the 50-80% range. The proposed Clear Lake Storm Water Improvement Project will address the remaining seven outlet sites and take place over three years. The first year will consist of performing engineering and design of storm water best management practices (BMPs) at the seven outlet sites to determine if a cost effective solution exists for each. Years two and three will consist of installing two storm water improvements each year to implement the most cost effective BMPs at a minimum of four of the seven sites. The grant request addresses one of the main priorities of the Iowa Watershed Improvement Grant: storm water runoff.

05036-013 Joint Drainage District 1-10 Fnial Report

Install new tile system to eliminate the use of 20 agricultural drainage wells (ADW) that have been ruled by EPA to be "vunerable to spills from manure lagoons, direct discharges from septic tanks and accidental releases of materails used in farming". The project includes a tile system that will provide an outlet for all the tile connected to the ADWs and close all 20 ADWs in the drainage district. The pipe ranges in size from 42" to 8" in diameter and approximately 48,000 lineal feet to be installed.

1010-005 Bloody Run Creek Watershed Final Report

The City of Marquette lies in the 65,000 acre Mississippi River watershed, and is surrounded by steep bluffs. Though scenic, controlling water runoff during storm events presents significant challenges. Flash-flooding from the local watershed has plagued the city for decades. The people of Marquette have committed to preserve the water quality of key natural resources in the area including the Bloody Run Creek and associated wetlands by undertaking projects to control the spread of debris and sediment caused by excess runoff during area storm events. Following a July 2007 storm (over 8” of rain in 24 hours) which caused unprecedented flood damage, the City retained an engineering firm to study the area and provide recommendations to eliminate or greatly reduce uncontrolled runoff into the Bloody Run Creek wetland, infrastructure damage and personal property loss. Marquette has received Iowa Great Places designation, and has demonstrated its commitment to wetland preservation with the construction of Phase I of this water quality project. The Bench Area Storm Water Management Plan prepared by the City in 2008 made a number of recommendations to mitigate flash flooding by improving storm water conveyance paths, detention, and infrastructure within the Bench area. Due to steep slopes and rocky geography, infiltration based systems, though desirable, would not be an option over surface based systems. Runoff from the 240 acre watershed comes primarily from large, steep drainage areas to the south and west, flowing to the Bench area down three hillside routes; designated as South East, South Central and South West. Completion of Phase I, which included an increased storage capacity of the upper pond, addressed the South East and South Central areas. The increased upper pond capacity will now allow Phase II to proceed. Phase II will address runoff from the South West drainage area; which engineers have estimated to produce as much water volume as the South Central and South East areas combined. Total costs for Phase I are $1.45 million, of which Marquette has invested $775,000, and IJOBS funding contributed $677,000. Phase II costs are estimated at $617,000. WIRB funding support of $200,000 would expedite project completion, lessen the long term debt impact to the community and aid in the preservation of the Bloody Run Creek and adjoining wetlands more quickly than Marquette could accomplish on its own.

Two cases of interferon-alpha-induced sarcoidosis koebnerized along venous drainage lines : new pathogenic insights and review of the literature of interferon-induced sarcoidosis

La sarcoïdose est une affection inflammatoire granuiomateuse systémique d'origine inconnue touchant le plus fréquemment les poumons, le système lymphoïde, le foie, les yeux et la peau. Dans cet article, nous rapportons deux cas de sarcoïdose cutanée touchant les avant-bras de deux patients anciens toxicomanes traités par interféron-a et ribavirine pour une hépatite C chronique. Nous procédons à une revue de la littérature de la sarcoïdose induite par l'interféron et élaborons une nouvelle hypothèse pathogénique de l'effet Koebner dans la sarcoïdose cutanée. Dans le cas des deux patients que nous décrivons, la distribution des lésions cutanées coïncide avec les anciens sites d'injection d'héroïne le long des trajets veineux des deux avant- bras. Cette distribution unique de l'atteinte cutanée suggère que les dommages tissulaires induits par la répétition d'injections percutanées puissent représenter un terrain favorisant au développement local d'une sarcoïdose cutanée. Fait intéressant, il a été récemment démontré que les cellules dendritiques plasmacytoïdes - sous-type de cellules dendritiques généralement absent de la peau - infiltrent rapidement les sites de peau lésée. Ces cellules sont la source d'une production endogène d'interféron-a, cytokine connue pour promouvoir le processus de cicatrisation, mais également pour favoriser le développement de sarcoïdose chez des individus prédisposés. Ainsi, nous postulons que les lésions de sarcoïdose cutanée limitées le long des trajets veineux - sites préalables d'injection percutanée de drogues - peuvent résulter d'une expression locale supplémentaire d'interféron-a. Celle-ci serait en outre favorisée par le traitement de ribavirine dans le cadre de l'hépatite C, connu pour renforcer la production endogène d'interféron-a. L'identification de nombreuses cellules dendritiques plasmacytoïdes circonscrivant l'inflammation granuiomateuse sur la biopsie cutanée de l'un de nos patients semble être un argument dans ce sens, conforté par l'absence de corps étranger détecté en microscopie par lumière polarisée. Cette observation semble pouvoir représenter un point crucial dans la compréhension des mécanismes physiopathologiques à la base de l'infiltration des cicatrices cutanées par la sarcoïdose. Des investigations supplémentaires doivent encore être effectuées afin de confirmer cette hypothèse.

Geologic and Engineering Properties of Pleistocene Materials in Iowa, Iowa Highway Research Board Bulletin No. 20, 1960

This bulletin is a compilation of the reports on completed research done for the Iowa State Highway Research Board Project HR-1, "The Loess and Glacial Till Materials of Iowa; an Investigation of Their Physical and Chemical Properties and Techniques for Processing Them to Increase Their All-Weather Stability for Road Construction.” The research, started in 1950, was done by the Iowa Engineering Experiment Station under its project 283-S. The project was supported by funds from the Iowa State Highway Commission.

Benchmarking project of the Department of Electrical Engineering of Lappeenranta University of Technology

Lappeenrannan teknillisen korkeakoulun sähkötekniikan osasto muutti 1.8.2005 sähkötekniikan tutkinnon kaksiportaiseksi ja vastaamaan näin Bologna-prosessia ja Suomen yliopistolainsäädäntöä. Tutkinnonuudistuksen myötä osasto haluaa varmistaa sähkötekniikan tutkintojen laadun ja vertailtavuuden sekä parantaa opiskelijoidensekä henkilökunnan liikkuvuutta. Tutkintojen laatu ja vertailtavuus osoitetaan sähkötekniikan osaston benchmark-projektilla, jossa kerätään tietoja maisteri- ja tohtorintutkintoa tarjoavista eurooppalaisista yliopistoista. Diplomityö käsittää BM-projektin kolmannen vaiheen suunnittelun ja toteutuksen sisältäen teoriaa benchmark-projekteille tyypillisistä toimintatavoista. Hyväksi havaitut menetelmiä, kuten kyselyt ja matriisit, on tässä työssä otettu soveltuvin osin sähkötekniikan osaston BM-projektin työkaluiksi. Diplomityössä analysoidaan työkalujen avulla BM-kumppaneilta kerättyjä tietoja sekä esitetään ratkaisuja, miten sähkötekniikan osastolla voidaan jatkaa parhaiden toimintatapojen löytämistä.