Multiphase flow and transport in fractured rock is of importance to many practical and engineering applications. In the field of groundwater hydrology an issue of significant environmental concern is the release of dense non-aqueous phase liquids (DNAPLs) which can cause long-term groundwater contamination in fractured aquifers. This study deals with two fundamental processes – fluid displacement and interphase mass transfer – concerning the behavior of the multiphase contaminants in fractured media. The focus of this work has been placed on improving the current understanding of small-scale (single fracture) physics by a combined effort of numerical modeling analysis, laboratory experiments and model development. This thesis contributes to the improved understanding through several aspects. Firstly, the effect of aperture variability, as characterized by geostatistical parameters such as standard deviation and correlation length, on the DNAPL entrapment, dissolution and source-depletion behaviors in single fractures was revealed. Secondly, a novel, generalized approach (adaptive circle fitting approach) to account for the effect of in-plane curvature of fluid-fluid interfaces on immiscible fluid displacement was developed; the new approach has demonstrated good performance when applied to simulate previously published experimental data. Thirdly, the performance of a continuum-based two-phase flow model and an invasion percolation model was compared for modeling fluid displacement in a variable-aperture fracture and the dependence of fracture-scale capillary pressure – saturation relationships on aperture variability was studied. Lastly, through experimental studies and mechanistic numerical modeling of DNAPL dissolution, kinetic mass transfer characteristics of two different entrapment configurations (residual blobs and dead-end pools) were investigated. The obtained understanding from this thesis will be useful for predictive modeling of multiphase contaminant behavior at a larger (fracture network) scale. / Flerfasflöde och ämnestransport i sprickigt berg är av betydelse för många praktiska och tekniska problem. Tunga, svårlösliga organiska vätskor (engelska: dense non-aqueous phase liquids: DNAPLs; t.ex. klorerade lösningsmedel) kan orsaka långvarig förorening av vattenresurser, inklusive akviferer i sprickigt berg, och utgör ett viktigt miljöproblem inom grundvattenhydrologin. Denna studie behandlar två fundamentala processer för spridning av flerfasföroreningar i sprickiga medier – utbredning av den organiska vätskan och massöverföring mellan organisk vätska och vatten. Arbetet har fokuserat på att förbättra nuvarande kunskap om de fysikaliska processerna på liten skala (enskilda sprickor) genom en kombination av numerisk modellering, laboratorieexperiment och modellutveckling. Avhandlingen har bidragit till utökad processförståelse i flera avseenden. För det första har arbetet belyst effekterna av sprickaperturens variabilitet, uttryckt med geostatistiska parametrar som standardavvikelse och rumslig korrelationslängd, på fastläggning och lösning av organiska vätskor i enskilda sprickor, samt utmattningsbeteendet hos dessa källor till grundvattenförorening. För det andra har en ny, generell metod (adaptiva cirkelpassningsmetoden) för att ta hänsyn till effekten av krökningen av gränsytan mellan organisk vätska och vatten i sprickplanet utvecklats; denna metod har visats fungera väl i simuleringar av tidigare publicerade experimentella data. För det tredje, har en jämförelse gjorts mellan en kontinuumbaserad tvåfasflödesmodell och en invasions-perkolationsmodell med avseende på hur väl de kan simulera tvåfasflöde i en spricka med varierande apertur. Här studerades även hur relationen mellan kapillärtryck och mättnadsgrad på sprickplansskala beror av variabiliteten i sprickapertur. Till sist undersöktes lösning av den organiska vätskan i grundvatten för två fastläggningsscenarier (fastläggning i immobila droppar och ansamling i fällor – ”återvändssprickor”) både genom experiment och mekanistisk numerisk modellering. Kunskapen som tagits fram i denna avhandling bedöms vara användbar även för att modellera spridningen av flerfasföroreningar på större (spricknätverks-) skalor.

No description available.

A new measurement technique called fiber optic reflectometer (FOR) was developed to investigate multiphase flows. The principle and setup of the FOR technique were introduced and applied to various experiments. Based on the coherently mixed signal between the Fresnel reflection off the fiber-liquid interface and the scattered signal off the object, such as a gas bubble, and a solid particle, this single probe technique is capable of simultaneously measuring the velocity of the object with a high accuracy and the phase of the fluid. In addition, bubble diameter, velocity, and void fraction were measured directly. By means of a simple modification of the FOR technique, solute concentration and refractive index change were measured with a greatly improved accuracy. This modified technique was used for measuring of a NaCl concentration in deionized water to validate a new normalization technique. In the second part of this thesis, a plunging breaking wave in deep water has been studied. Using the wave focusing method, a strong plunging breaker was generated with accuracy in the deep water condition in a two-dimensional wave tank. It was possible to describe the breaking process in detail using a high speed camera with a frame rate of 500 or 1000 fps. Four kinds of experimental techniques were employed or developed to investigate the plunging breaker. Bubble image velocimetry (BIV) and particle image velocimetry (PIV) were used to measure the velocity fields. The velocity fields of the highly aerated region were obtained from the BIV measurements. In addition, the modified PIV technique is capable of measuring the velocities in the entire flow field including the aerated region. Mean and turbulent properties were obtained by the ensemble average. The mean velocity, mean vorticity, and mean kinetic energy were examined over the entire flow field. In addition, the Reynolds stresses and turbulent kinetic energy were calculated with high temporal and spatial resolutions. Free surface elevation was obtained from wave gauge measurements. BIV and PIV images were also used to obtain the free surface elevation and the boundary of the aerated region for more accurate results. The FOR technique was used to obtain the void ratio at each splash-up region. Compressibility of the plunging breaker was considered. Mass flux, momentum flux, kinetic energy, and Reynolds stresses at each FOR station were recalculated using the void ratio obtained from the FOR measurements. All terms at the first splash-up region were highly overestimated more than 100 percent unless the void ratio was applied to the calculation of fluxes and energies. Compared with the fully developed first splash-up region, the overestimation at the second and third splash-up was less significant. However, most terms were overestimated by 20~30 percent when the void ratio was not considered.

High explosive charges when detonated ensue in a flow field characterized by several physical phenomena that include blast wave propagation, hydrodynamic instabilities, real gas effects, fluid mixing and afterburn effects. Solid metal particles are often added to explosives to augment the total impulsive loading, either through direct bombardment if inert, or through afterburn energy release if reactive. These multiphase explosive charges, termed as heterogeneous explosives, are of interest from a scientific perspective as they involve the confluence and interplay of various additional physical phenomena such as shock-particle interaction, particle dispersion, ignition, and inter-phase mass, momentum and energy transfer. In the current research effort, chemical explosions in multiphase environments are investigated using a robust, state-of-the-art Eulerian-gas, Lagrangian-solid methodology that can handle both the dense and dilute particle regimes. Explosions into ambient air as well as into aluminum particle clouds are investigated, and hydrodynamic instabilities such as Rayleigh- Taylor and Richtmyer-Meshkov result in a mixing layer where the detonation products mix with the air and afterburn. The particles in the ambient cloud, when present, are observed to pick up significant amounts of momentum and heat from the gas, and thereafter disperse, ignite and burn. The amount of mixing and afterburn are observed to be independent of particle size, but dependent on the particle mass loading and cloud dimensions. Due to fast response times, small particles are observed to cluster as they interact with the vortex rings in the mixing layer, which leads to their preferential ignition/ combustion. The total deliverable impulsive loading from heterogeneous explosive charges containing inert steel particles is estimated for a suite of operating parameters and compared, and it is demonstrated that heterogeneous explosive charges deliver a higher near-field impulse than homogeneous explosive charges containing the same mass of the high explosive. Furthermore, particles are observed to introduce significant amounts of hydrodynamic instabilities in the mixing layer, resulting in augmented fluctuation intensities and fireball size, and different growth rates for heterogeneous explosions compared to homogeneous explosions. For aluminized explosions, the particles are observed to burn in two regimes, and the average particle velocities at late times are observed to be independent of the initial solid volume fraction in the explosive charge. Overall, this thesis provides useful insights on the role played by solid particles in chemical explosions.

An understanding of the motion of soft capsules in microchannels is useful for a number applications. This knowledge can be used to develop devices to sort biological cells based on their size and stiffness. For example, cancer cells have a different stiffness from healthy cells and thus can be readily identified. Additionally, devices can be developed to detect flaws in synthetic particles. Using a 3D hybrid lattice Boltzmann and lattice spring method, the motion of rigid and soft capsules in a pressure-driven microfluidic flow was probed. The effect of inertial drift is evaluated in channels different Reynolds numbers. Other system parameters such as capsule elasticity and channel size are also varied to determine their effect. The equilibrium position of capsules in the channel is also obtained. The equilibrium position of rigid and soft capsules depends on the relative particle size. If the capsule is small, the equilibrium position is found to be closer to the channel wall. Conversely, for larger capsules, the equilibrium position is closer to the channel centerline. The capsule stiffness affects the magnitude of the cross-stream drift velocity. For a given Reynolds number, the equilibrium position of softer capsules is closer to the channel centerline. However, It is found that the equilibrium position of soft capsules is insensitive to the magnitude of the Reynolds number.

Multiphase flow denotes the simultaneous flow of two or more physically distinct and immiscible substances and it can be widely found in several engineering applications, for instance, power generation, chemical engineering and crude oil extraction and processing. In many of those applications, multiphase flows determine safety and efficiency aspects of processes and plants where they occur. Therefore, the measurement and imaging of multiphase flows has received much attention in recent years, largely driven by a need of many industry branches to accurately quantify, predict and control the flow of multiphase mixtures. Moreover, multiphase flow measurements also form the basis in which models and simulations can be developed and validated. In this work, the use of electrical impedance techniques for multiphase flow measurement has been investigated. Three different impedance sensor systems to quantify and monitor multiphase flows have been developed, implemented and metrologically evaluated. The first one is a complex permittivity needle probe which can detect the phases of a multiphase flow at its probe tip by simultaneous measurement of the electrical conductivity and permittivity at up to 20 kHz repetition rate. Two-dimensional images of the phase distribution in pipe cross section can be obtained by the newly developed capacitance wire-mesh sensor. The sensor is able to discriminate fluids with different relative permittivity (dielectric constant) values in a multiphase flow and achieves frame frequencies of up to 10 000 frames per second. The third sensor introduced in this thesis is a planar array sensor which can be employed to visualize fluid distributions along the surface of objects and near-wall flows. The planar sensor can be mounted onto the wall of pipes or vessels and thus has a minimal influence on the flow. It can be operated by a conductivity-based as well as permittivity-based electronics at imaging speeds of up to 10 000 frames/s. All three sensor modalities have been employed in different flow applications which are discussed in this thesis. The main contribution of this research work to the field of multiphase flow measurement technology is therefore the development, characterization and application of new sensors based on electrical impedance measurement. All sensors present high-speed capability and two of them allow for imaging phase fraction distributions. The sensors are furthermore very robust and can thus easily be employed in a number of multiphase flow applications in research and industry.

Thesis (Ph. D.)--Mississippi State University. Department of Engineering. / Title from title screen. Includes bibliographical references.

During the production and transportation of oil and gas, flow assurance issues may occur due to the solid deposits that are formed and carried by the flowing fluid. Solid deposition may cause serious damage and possible failure to production equipment in the flow lines. The major flow assurance problems that are faced in the fields are concerned with asphaltene, wax and scale deposition, as well as hydrate formations. Hydrates, wax and asphaltene deposition are mostly addressed in deep-water environments, where fluid flows through a long path with a wide range of pressure and temperature variations (Hydrates are generated at high pressure and low temperature conditions). In fact, a large change in the thermodynamic condition of the fluid yields phase instability and triggers solid deposit formations. In contrast, scales are formed in aqueous phase when some incompatible ions are mixed. Among the different flow assurance issues in hydrocarbon reservoirs, asphaltenes are the most complicated one. In fact, the difference in the nature of these molecules with respect to other hydrocarbon components makes this distinction. Asphaltene molecules are the heaviest and the most polar compounds in the crude oils, being insoluble in light n-alkenes and readily soluble in aromatic solvents. Asphaltene is attached to similarly structured molecules, resins, to become stable in the crude oils. Changing the crude oil composition and increasing the light component fractions destabilize asphaltene molecules. For instance, in some field situations, CO₂ flooding for the purpose of enhanced oil recovery destabilizes asphaltene. Other potential parameters that promote asphaltene precipitation in the crude oil streams are significant pressure and temperature variation. In fact, in such situations the entrainment of solid particulates in the flowing fluid and deposition on different zones of the flow line yields serious operational challenges and an overall decrease in production efficiency. The loss of productivity leads to a large number of costly remediation work during a well life cycle. In some cases up to $5 Million per year is the estimated cost of removing the blockage plus the production losses during downtimes. Furthermore, some of the oil and gas fields may be left abandoned prematurely, because of the significance of the damage which may cause loss about $100 Million. In this dissertation, we developed a robust wellbore model which is coupled to our in-house developed compositional reservoir model (UTCOMP). The coupled wellbore/reservoir simulator can address flow restrictions in the wellbore as well as the near-wellbore area. This simulator can be a tool not only to diagnose the potential flow assurance problems in the developments of new fields, but also as a tool to study and design an optimum solution for the reservoir development with different types of flow assurance problems. In addition, the predictive capability of this simulator can prescribe a production schedule for the wells that can never survive from flow assurance problems. In our wellbore simulator, different numerical methods such as, semi-implicit, nearly implicit, and fully implicit schemes along with blackoil and Equation-of-State compositional models are considered. The Equation-of-State is used as state relations for updating the properties and the equilibrium calculation among all the phases (oil, gas, wax, asphaltene). To handle the aqueous phase reaction for possible scales formation in the wellbore a geochemical software package (PHREEQC) is coupled to our simulator as well. The governing equations for the wellbore/reservoir model comprise mass conservation of each phase and each component, momentum conservation of liquid, and gas phase, energy conservation of mixture of fluids and fugacity equations between three phases and wax or asphaltene. The governing equations are solved using finite difference discretization methods. Our simulation results show that scale deposition is mostly initiated from the bottom of the wellbore and near-wellbore where it can extend to the upper part of the well, asphaltene deposition can start in the middle of the well and the wax deposition begins in the colder part of the well near the wellhead. In addition, our simulation studies show that asphaltene deposition is significantly affected by CO₂ and the location of deposition is changed to the lower part of the well in the presence of CO₂. Finally, we applied the developed model for the mechanical remediation and prevention procedures and our simulation results reveal that there is a possibility to reduce the asphaltene deposition in the wellbore by adjusting the well operation condition. / text

Hydrocarbon reserves are generally produced through wells drilled into reservoir pay zones. During production, gas liberation from the oil phase occurs due to pressure decline in the wellbore. Thus, we expect multiphase flow in some sections of the wellbore. As a multi-phase/multi-component gas-oil mixture flows from the reservoir to the surface, pressure, temperature, composition, and liquid holdup distributions are interrelated. Modeling these multiphase flow parameters is important to design production strategies such as artificial lift procedures. A wellbore fluid flow model can also be used for pressure transient test analysis and interpretation. Considering heat exchange in the wellbore is important to compute fluid flow parameters accurately. Modeling multiphase fluid flow in the wellbore becomes more complicated due to heat transfer between the wellbore fluids and the surrounding formations. Due to mass, momentum, and energy exchange between the wellbore and the reservoir, the wellbore model should be coupled with a numerical reservoir model to simulate fluid flow accurately. This model should be non-isothermal to consider the effect of temperature. Our research shows that, in some cases, ignoring compositional effects may lead to errors in pressure profile prediction for the wellbore. Nearly all multiphase wellbore simulations are currently performed using the "black oil" approach. The primary objective of this study was to develop a non-isothermal wellbore simulator to model transient fluid flow and temperature and couple the model to a reservoir simulator called General Purpose Adaptive Simulator (GPAS). The coupled wellbore/reservoir simulator can be applied to steady state problems, such as production from, or injection to a reservoir as well as during transient phenomena such as well tests to accurately model wellbore effects. Fluid flow in the wellbore may be modeled either using the blackoil approach or the compositional approach, as required by the complexity of the fluids. The simulation results of the new model were compared with field data for pressure gradients and temperature distribution obtained from wireline conveyed pressure recorder and acoustic fluid level measurements for a gas/oil producer well during a buildup test. The model results are in good agreement with the field data. Our simulator gave us further insights into the wellbore dynamics that occur during transient problems such as phase segregation and counter-current multiphase flow. We show that neglecting these multiphase flow dynamics would lead to unreliable results in well testing analysis.

The phenomenon of particle retention in granular materials has a wide range of implications. For agricultural operations, these particles can be contaminants transported through the ground that can eventually reach to aquifers, consequently contaminating the water. In oil reservoirs, these particles can be clays that get detached from the rock and migrate with the flow after a change of pressure, plugging the reservoir with the consequent reduction in permeability. These particles can also be traceable nanoparticles, introduced in the reservoir with the purpose of identifying bypassed oil. For all these reasons it is important to understand the mechanisms that contribute to the transport and retention of these particles. In this dissertation the retention of micro and nano size particles was investigated. In saturated model sediments (sphere packs), we analyzed the retention of particles by the mechanism of straining (size exclusion). The analysis focused on experiments reported in the literature in which particles smaller than the smallest pore throats were retained in the sediment. The analysis yields a mechanistic explanation of these observations, by indentifying the retention sites as gaps between pairs of sediment grains. A predictive model was developed that yields a relationship between the straining rate constant and particle size in agreement with the experimental observations. In unsaturated granular materials, the relative contributions of grain surfaces, interfacial areas and contact lines between phases to the retention of colloidal size particles were investigated. An important part of this analysis was the identification and calculation of the length of the contact lines between phases. This estimation of contact line lengths in porous media is the first of its kind. The algorithm developed to compute contact line length yielded values consistent with observations from beads pack and real rocks, which were obtained independently from analysis of high resolution images. Additionally, the predictions of interfacial areas in granular materials were consistent with an established thermodynamic theory of multiphase flow in porous media. Since there is a close relationship between interfacial areas and contact lines this supports the accuracy of the contact line length estimations. Predictions of contact line length and interfacial area in model sediments, combined with experimental values of retention of colloidal size particles in columns of glass beads suggested that it is plausible for interfacial area and contact line to contribute in the same proportion to the retention of particles. The mechanism of retention of surface treated nanoparticles in sedimentary rocks was also investigated, where it was found that retention is reversible and dominated by attractive van der Waals forces between the particles and the rock’s grain surfaces. The intricate combination of factors that affect retention makes the clear identification of the mechanism responsible for trapping a complex task. The work presented in this dissertation provides significant insight into the retention mechanisms in relevant scenarios. / text