External loop airlift reactors are widely used for biochemical applications such as syngas fermentation and wastewater treatment. To further understand the inherent gas-liquid flow physics within the reactors, computational modeling and simulations of hydrodynamics for air-water external loop airlift reactors were investigated. The gas-liquid flow dynamics in a bubble column were simulated using a FORTRAN code developed by Los Alamos National Laboratory, CFDLib, which employs an Eulerian-Eulerian ensemble averaged method. A two-dimensional Cartesian coordinate system was used to conduct an extensive grid resolution study; it was found that grid cells smaller than the bubble diameter produced unstable solutions. Next, closure models for drag force and turbulent viscosity were investigated for a simple bubble column geometry. The effects of using a bubble pressure model and two drag coefficient models, the White model and the Schiller-Naumann model, were investigated. The bubble pressure model performed best for homogeneous (low velocity) flows and the Schiller-Naumann model was best for all flow regimes. Based on the studies for bubble column flows, an external loop airlift reactor was simulated using both two- and three-dimensional coordinates and results for gas holdup and riser velocity agreed better with experimental data for the 3D simulations. It was concluded that when performing 2D and 3D simulations, care must be taken when specifying the effective bubble diameter size, especially at high flow rates. Population balance models (PBM) for bubble break-up and coalescence were implemented into CFDLib, validated with experiments, and simulated for the external loop airlift reactor at high inlet superficial gas velocities. The PBM predictions for multiple bubble sizes were comparable with the single bubble size simulations; however, the PBM simulations better predicted the formation of the gas bubble in the downcomer. The 3D PBM simulations also gave better predictions for the average bubble diameter size in the riser. It was concluded that a two-dimensional domain is adequate for gas-liquid flow simulations of a simple bubble column geometry, whereas three-dimensional simulations are required for the complex airlift reactor geometry. To conclude, a two-fluid Eulerian-Eulerian model coupled with a PBM is needed for quantitative as well as physical predictions of gas-liquid external loop airlift reactor flows at high inlet superficial gas velocities. / Ph. D.

Geological carbon dioxide (CO2) sequestration (GCS) in deep saline aquifers is viewed as a viable solution to dealing with the impact of anthropogenic CO2 emissions on global warming. The trapping mechanisms that control GCS include capillary trapping, structural trapping, dissolution trapping, and mineral trapping. Wettability and density-driven convection play an important role in GCS, because wettability significantly affects the efficiency of capillary trapping, and density-driven convection greatly decreases the time scale of dissolution trapping. This work focuses on the role of wettability on multiphase flow in porous media, density-driven convection in porous media, and their implications for GCS in deep saline aquifers. Wettability is a critical control over multiphase fluid flow in porous media. However, our understanding on the wettability heterogeneity of a natural rock and its effect on multiphase fluid flow in a natural rock is limited. This work innovatively models the heterogeneous wettability of a rock as a correlated random field. The realistic wetting condition of a natural rock can be reconstructed with in-situ measurements of wettability on the internal surfaces of the rock. A Bentheimer sandstone was used to demonstrate the workflow to model and reconstruct a wettability field. Relative permeability, capillary pressure-water saturation relation are important continuum-scale properties controlling multiphase flow in porous media. This work employed lattice Boltzmann method to simulate the displacement process. We found that pore-scale surface wettability heterogeneity caused noticeable local scCO2 and water redistributions under less water-wet conditions at the pore scale. At the continuum scale, the capillary pressure-water saturation curve under the heterogeneous wetting condition was overall similar to that under the homogeneous wetting condition. This suggested that the impact of local wettability heterogeneity on the capillary pressure-water saturation curve was averaged out at the entire-sample scale. The only difference was that heterogeneous wettability led to a negative entry pressure at the primary drainage stage under the intermediate-wet condition. The impact of pore-scale wettability heterogeneity was more noticeable on the relative permeability curves. Particularly, the variation of the scCO2 relative permeability curve in the heterogeneous wettability scenario was more significant than that in the homogenous wettability scenario. Results showed that higher wettability heterogeneity (i.e., higher standard deviation and higher correlation length) increased the variations in the CO2/brine relative permeability curves. Dissolution of CO2 into brine is a primary mechanism to ensure the long-term security of GCS. CO2 dissolved in brine increases the CO2-brine solution density and thus can cause downward convection. Onset of density-driven instability and onset of convective dissolution are two critical events in the transition process from a diffusion-dominated regime to a convection-dominated regime. In the laboratory, we developed an empirical correlation between light intensity and in-situ solute concentration. Based on the novel and well-controlled experimental methods, we measured the critical Rayleigh-Darcy number and critical times for the onset of density-driven instability and convective dissolution. To further investigate the impact of permeability heterogeneity on density-driven convection, a three-dimensional (3D) fluidics method was proposed to advance the investigation on density-driven convection in porous media. Heterogeneous porous media with desired spatial correlations were efficiently built with 3D-printed elementary porous blocks. In the experiments, methanol-ethylene-glycol (MEG), was used as surrogate fluid to CO2. The heterogeneous porous media were placed in a transparent tank allowing visual observations. Results showed that permeability structure controlled the migration of MEG-rich water. Permeability heterogeneity caused noticeable uncertainty in dissolution rates and uncertainty in dissolution rates increases with correlation length. To sum up, this work comprehensively employed novel experimental methods and large-scale direct simulations to investigate the sequestration of CO2 in saline aquifers at a pore scale and a continuum scale. The findings advanced our understanding on the role of wettability heterogeneity and permeability heterogeneity on GCS in deep saline aquifers. / Doctor of Philosophy / Global warming caused by anthropogenic CO2 emissions is a pressing issue to address of our time. The storage of CO2 in deep saline aquifers is a promising solution because of saline aquifers' vast storage capacity. Property heterogeneity exists extensively in saline aquifers from a continuum scale to a pore scale. The implications of pore-scale wettability heterogeneity and continuum-scale permeability heterogeneity for the storage of CO2 in saline aquifers are not clear. This work is to employ novel experimental methods and powerful simulation tools to investigate the role of wettability heterogeneity and permeability heterogeneity on the storage of CO2 in saline aquifers. This work measured contact angles on the scanned micro-CT images of a Bentheimer sandstone after a CO2 flooding. A correlated lognormal wettability model was put forward with the statistical information of the contact angle measurements. Simulations on the CO2/brine flow in the Bentheimer sandstone were performed. Results showed that the wettability heterogeneity caused noticeable redistributions of CO2/brine compared to scenarios under homogeneous wettability. Impact of wettability on capillary pressure-water saturation curve was not noticeable because the effects were averaged out through the entire rock sample. The standard deviation and correlation length caused variations on the relative permeabilities. This means that we need to take them into consideration in simulating the migration of CO2 in saline aquifers at a reservoir scale. After CO2 pools beneath the impermeable cap rock, dissolution of CO2 into brine dominates the trapping process. Convection caused by CO2 dissolution can greatly accelerate the dissolution rate. The onset of convection is a critical issue and lack of experimental evidence. This work firstly determined the onset time of instability. To further investigate the heterogeneity on the convection, this work proposed a 3D-print-based method to efficiently build heterogeneous porous media with a designed permeability distribution. The experiments were conducted, and results showed that heterogeneity structure of porous media can cause great variations on the dissolution rate of CO2. The findings of this work advanced our understanding on the migration of CO2 in saline aquifers, provided solid basis for assessment and decision on the storage of CO2 into saline aquifers.

Void fraction is an essential parameter for understanding the interfacial structure, and heat and mass transfer mechanisms in various gas-liquid flow systems. It becomes critically important to accurately measure void fraction as advanced high fidelity two-phase flow models require high-quality validation data. However, void fraction measurement remains a challenging task to date due to the complexity and rapid-changing characteristic of the gas-liquid boundary flow structure. This study aims to develop an advanced void fraction measurement system based on x-ray and fast line detector technologies. The dissertation has covered the major components necessary to develop a complete measurement system. Spectral analysis of x-ray attenuation in two-phase flow has been performed, and a new void fraction model is developed based on the analysis. The newly developed pixel-to-radial conversion algorithm is capable of converting measured void fraction along with the detector array to the radial distribution in a circular pipe for a wide range of void fraction conditions. The x-ray system attains the radial distributions of key measurable factors such as void fraction and gas velocity. The data are compared with the double-sensor conductivity probe and gas flowmeter for various flow conditions. The results show reasonable agreements between the x-ray and the other measurement techniques. Finally, various 2-D tomography algorithms are implemented for the non-axisymmetric two-phase flow reconstruction. A comprehensive summary of classical absorption tomography for the two-phase flow study is provided. An in-depth sensitivity study is carried out using synthetic bubbles, aiming to investigate the effect of various uncertainty factors such as background noise, off-center shift, void profile effect, etc. The sensitivity study provides a general guideline for the performance of existing 2-D reconstruction algorithms. / Doctor of Philosophy / Gas-liquid flow phenomenon exists in an extensive range of natural and engineering systems, for example, hydraulic pipelines in a nuclear reactor, heat exchanger, pump cavitation, and boilers in the gas-fired power stations. Accurate measurement of the void fraction is essential to understand the behaviors of the two-phase flow phenomenon. However, measuring void fraction distribution in two-phase flow is a difficult task due to its complex and fast-changing interfacial structure. This study developed a comprehensive suite of the non-intrusive x-ray measurement techniques, and a pixel-to-radial conversion algorithm to process the line- and time-averaged void fraction information. The newly developed algorithm, called the Area-based Onion-Peeling (ABOP) method, can convert the pixel measurement to the radial void fraction distribution, which is more useful for studying and modeling axisymmetric flows. Various flow conditions are measured and evaluated for the benchmarking of the algorithm. Finally, classical 2-D reconstruction algorithms are investigated for the void fraction measurement in non-axisymmetric flows. A comprehensive summary of the performance of these algorithms for a two-phase flow study is provided. An in-depth sensitivity study using synthetic bubbles has been performed to examine the effect of uncertainty factors and to benchmark the algorithms for the non-axisymmetric flows.

This thesis presents an investigation into the motions of the horseshoe vortices and the passage vortex, within a plane turbine blade cascade. Fluid motion was determined using a tracer gas technique. Ethylene was injected into the pressure-side and suction-side legs of the horseshoe vortex, near the leading edge of the cascade. Ethylene concentrations were determined at two downstream locations using a flame ionization detector. It was found that the pressure-side leg of the horseshoe vortex moved toward the suction side of the passage, starting the formation of the passage vortex, and was distributed throughout the passage vortex. The suction-side leg of the horseshoe vortex convected once around the periphery of the passage vortex before passing the cascade trailing edge. Downstream of the trailing edge, most of the fluid from the suction-side leg diffused into the passage vortex. However, twice as much fluid from the suction-side leg, as opposed to the pressure-side leg, mixed within the blade wake. At a location 40% of the axial chord downstream of the trailing edge, the passage vortex (shown previously to account for 60% of the overall total pressure losses) contained over 65% of the fluids from both legs. / M.S.

Biomass is considered a biorenewable alternative energy resource that can potentially reduce the use of natural gas and provide low cost power production or process heating needs. Biomass hydrodynamics in a fluidized bed are extremely important to industries that are using biomass material in gasfication processes to yield high quality producer gas. However, biomass particles are typically difficult to fluidize due to their peculiar shape and a second inert material, such as sand, is typically added to the bed. The large differences in size and density between the biomass and inert particles lead to nonuniform distribution of the biomass within the fluidized bed, and particle interactions and mixing become major issues. The main goal of this research was to use CFD as a tool for modeling and analyzing the hydrodynamic behavior of biomassas a single material or as part of a mixture in a fluidized bed. The first part of this research focused on the characterization of biomass particles in a fluidized bed and validation of a numerical model with experimental results obtained from pressure measurements and CT and X-ray radiograph images. For a 2D fluidized bed of glass beads, the pressure drop, void fraction and mean bed height expansion were in quantitative agreement between the experiments and simulations using Syamlal-O'Brien and Gidaspow drag models. It was encouraging that the Gidaspow model predictions were in close agreement because the model does not require knowing the minimum fluidization as an input. Ground walnut shells were used to represent biomass because the material fluidizes uniformly and is classified as a Geldart type B particle. Two-dimensional simulations of ground walnut shells were analyzed to determine parameters that cannot easily be measured experimentally. The parametric study for ground walnut shell indicated that the material can be characterized with a medium sphericity (~0.6) and a relatively large coefficient of restitution (~0.85). In the second part of this work numerical simulations of a ground walnut shell fluidizing bed with side air injection were compared to CT data for the gas-solid distribution to demonstrate the quantitative agreement for bed fluidization. The findings showed that 2D simulations overpredicted the fluidized bed expansion and the results did not demonstrate a uniformly fluidizing bed. The 3D simulations compared well for all cases. This study demonstrates the importance of using a 3D model for a truly 3D flow in order to capture the hydrodynamics of the fluidized bed for a complicated flow and geometry. Finally, CFD modeling of pressure fluctuations was performed on sand and cotton-sand fluidized beds operating at inlet velocities ranging from 1.0-9.0Umf with the objective of predicting characteristic features of bubbling, slugging, and turbulent fluidization regimes. It was determined that the fluidized bed can be modeled using MUSCL discretization and the Ahmadi turbulence model. Three-dimensional sand fluidized beds were simulated for different fluidization regimes. Fluidized beds for all the regimes behaved as second-order dynamic systems. Bubbling fluidized beds showed one broad peak with a maximum at 2.6 Hz while slugging and turbulent showed two distinct peaks. It was observed that the peak at low frequency increased in magnitude as the flow transitioned from a slugging to a turbulent fluidization regime. CFD simulations of fluidized beds with the purpose of studying pressure fluctuations have demonstrated to be a useful tool to obtain hydrodynamic information that will help determine the fluidization regime. Prediction of slugging and turbulent fluidization regimes using CFD have not been reported to date. The work presented here is the first of its kind and can be an important advantage when designing a reactor and evaluating different operation conditions without the need to test them in a pilot plant or a prototype. / Ph. D.

Measuring void fraction distribution in two-phase flow has been a challenging task for many decades because of its complex and fast-changing interfacial structure. In this study, a non-intrusive X-ray measurement system is developed and calibrated to mitigate this challenge. This approach has several advantages over the conventional methods such as the multi-sensor conductivity probe, wire-mesh sensor, impedance void meter, or direct optical imaging. The X-ray densitometry technique is non-intrusive, insensitive to flow regime changes, capable of measuring high temperature or high-pressure flows, and has reasonable penetration depth. With the advancement of detector technology, the system developed in this work can further achieve high spatial resolution (100 micron per pixel) and high temporal resolution (1000 frames per second). This work mainly focuses on the following aspects of the system development: establishing a geometrical model for the line detector system, conducting spectral analysis for X-ray attenuation in two-phase flow, and performing calibration tests. The geometrical model has considered the measurement plane, geometry of the test-section wall and flow channel, relative position of the X-ray source and detector pixels. By assuming axisymmetry, an algorithm has been developed to convert void fraction distribution along the detector pixels to the radial void profile in a circular pipe. The X-ray spectral analysis yielded a novel prediction model for non-chromatic X-rays and non-uniform structure materials such as the internal two-phase flow which contains gas, liquid and solid wall materials. A calibration experiment has been carried out to optimize the detector conversion factor for each detector pixels. Finally, the data measured by the developed X-ray system are compared with the double-sensor conductivity probe and gas flow meter for sample bubbly flow and slug flow conditions. The results show reasonable agreement between these different measuring techniques. / Master of Science

We present 2-D numerical simulations of two phase flow in seafloor hydrothermal systems using the finite control volume numerical scheme FISHES. The FISHES code solves the coupled non-linear equations for mass, momentum, energy, and salt conservation in a NaCl-H2O fluid to model the seafloor hydrothermal processes. These simulations use homogeneous box geometries at a fixed seafloor pressure of 25 MPa with constant bottom temperature boundary conditions that represent a sub-axial magma chamber to explore the effects of permeability, maximum bottom temperature and system depth on the evolution of vent fluid temperature and salinity, and heat output. We also study the temporal and spatial variability in hydrothermal circulation. The two-phase simulation results show that permeability plays an important role in plume structure and heat output of hydrothermal systems, but it has little effect on vent fluid temperature and salinity, given the same bottom temperature. For some permeability values, multiple plumes can vent at the seafloor above the simulated magma chamber. Temporal variability of vent fluid temperature and salinity and the complexity of phase separation suggest that pressure and temperature conditions at the top of the axial magma chamber cannot be easily inferred from vent fluid temperature and salinity alone. Vapor and brine derived fluids can vent at the seafloor simultaneously, even from neighboring locations that are fed by the same plume. / Master of Science

Fluidized bed reactor technology has been widely adopted within the industry as vital component for numerous manufacturing, power generation and gasification processes due to its enhanced mixing characteristics. Computational modeling of fluidized bed hydrodynamics is a significant challenge that has to be tackled for increasing predictive accuracy. The distributor plate of a fluidized bed is typically modeled using a uniform inlet condition, when in reality the inlet is non-uniform inlet. The regions of bed mass that do not fluidize because of the non-uniform inlet conditions form deadzones and remain static between the jets. A new model based on the mass that contributes to the pressure drop is proposed to model a fluidized bed, and has been investigated for a cylindrical reactor for glass beads, ceramic single solids particles, and glass-ceramic, and ceramic-ceramic binary mixtures. The adjusted mass model was shown to accurately predict fluidization characteristics such as pressure drop and minimum fluidization velocity. The effectiveness of the adjusted mass model was further illustrated by applying it to fluidized beds containing coal, switchgrass, poplar wood, and cornstover biomass particles and coal-biomass binary mixtures. The adjusted mass model was further analyzed for bed expansion heights of different mixtures, and for solids distribution by analyzing the solids volume fraction. The effect of increasing the percent biomass in the mixture was also investigated. To further model the non-uniform inlet condition, two different distributor configurations with 5 and 9 jets was considered for a quasi-2D bed, and simulations were performed in both 2D and 3D. Fluidization characteristics and mixing of the bed were analyzed for the simulation. Furthermore, the deadzones formed due to multiple jet configurations of the distributor are quantified and their distributions over the plate were analyzed. / Ph. D.

Escoamentos bifásicos do tipo gás-líquido são encontrados em diversas aplicações. Na indústria petrolífera esses escoamentos são comumente encontrados em tubulações, como por exemplo, na extração e processamento de óleo e gás. No passado, diversos estudos foram realizados para investigar escoamentos bifásicos sobre condições controladas em pequenas e médias plantas de teste, onde os fluidos utilizados são frequentemente água e gás. Esses estudos dão suporte para o desenvolvimento de modelos de escoamento e correlações de engenharia. Nesse trabalho é descrito um sistema de aquisição e gerenciamento de dados experimentais, que possui 4 sensores distribuídos (expansível para até 32) conectados através de um barramento CAN (Controller Area Network). O sistema de gerenciamento foi desenvolvido utilizando tecnologias web, permitindo conexões simultâneas e acesso remoto através do navegador, provendo uma plataforma completa para aquisição, armazenamento, e compartilhamento de dados experimentais. O princípio de funcionamento dos sensores é baseado na medição da variação da condutividade elétrica do meio, permitindo diferenciar as fases líquida e gasosa num escoamento bifásico água-ar. Os sensores resistivos possuem baixo custo de fabricação, são minimamente invasivos, e possuem alta resolução temporal (2 kHz). O protocolo CAN foi escolhido por cobrir grandes distâncias, ter alta velocidade de transferência e baixo custo de implementação. Todos os dados experimentais coletados são armazenados em um banco de dados centralizado, permitindo visualização e acesso de qualquer computador que tenha acesso ao sistema. Resultados iniciais constataram o potencial da plataforma como um sistema flexível e confiável para investigação de escoamentos e gerenciamento de dados experimentais. / Two-phase gas-liquid flows are found in many industrial applications. In the petroleum industry, two phase flow are usually encountered confined to pipes, for instance, in oil and gas extraction and processing. In the past, several studies have been performed to investigate two-phase flows under controlled conditions at small and medium scale test facilities, whereas working fluids are often water and air. Such studies support the development of flow models and engineering correlations. In this work, a data acquisition and management system is introduced, in which 4 distributed sensors (expansible to up 32 sensors) are connected through a CAN bus. The data management system was developed using web technologies, allowing multiple simultaneous connections, providing a complete platform to acquire, store, visualize and share experimental data. The sensors’ operating principle is based on measurement of electrical conductivity of flowing media in order to distinguish liquid and gaseous phase in a two-phase air-water flow. The resistive sensors are low cost, minimally invasive, and have high temporal resolution (2 kHz). CAN bus protocol was chosen due to its robustness and possibility to cover long distances at high data transfer. All collected data are stored in a centralized database, allowing online visualization, and access to experimental data from external network. Initial results have shown the potential of the proposed platform as a flexible and reliable system for two-phase flow investigation and the associated data management.

Escoamentos bifásicos líquido-gás estão presentes na natureza e em muitas atividades industriais. Neste tipo de escoamento, as fases líquida e gasosa podem assumir diferentes configurações espaciais dentro da tubulação, chamadas padrões de escoamento. O escoamento bifásico líquido-gás em golfadas é o padrão de escoamento mais frequente nas aplicações industriais, ocorrendo em uma ampla faixa de velocidades das fases segundo os estudos de diversos autores. A modelagem matemática para o escoamento em golfadas compreende desde modelos simples em regime estacionário até modelos mais complexos, em regime transiente. E, para solução destes modelos são necessárias correlações empíricas e distribuições estatísticas dos parâmetros característicos do escoamento. Assim, no presente trabalho, vários modelos baseados em redes neurais artificiais são apresentados como suporte à caracterização dos parâmetros do escoamento bifásico em golfadas em função das séries temporais de fração de vazio obtidas experimentalmente. As séries temporais de fração de vazio são medidas com um par de sensores de malha de eletrodos instalado na seção de testes de uma planta experimental do NUEMUTFPR e descritas em Castillo (2015). A partir das séries temporais de fração de vazio medidas são calculados os parâmetros de interesse para o escoamento em estudo: comprimento da bolha alongada de gás, comprimento do pistão de líquido, velocidade de translação da bolha alongada e desvios padrões para essas variáveis. Essas variáveis medidas e calculadas são utilizadas para a obtenção de um conjunto de modelos baseados em rede neural artificial. Após obtenção dos modelos é realizado um estudo de simulação no qual esses modelos são usados para estimar os parâmetros que caracterizam o escoamento bifásico em golfadas. Análises detalhadas dos resultados mostraram que as variáveis relacionadas à fase gasosa são estimadas com maior acurácia que as variáveis relacionadas à fase líquida. Como aplicação imediata do modelo obtido, apresenta-se sua utilização como uma ferramenta de cálculo das condições iniciais para um modelo matemático fenomenológico de escoamento bifásico em golfadas com leve mudança de inclinação baseado no método de seguimento de pistões. O diferencial do presente trabalho está na predição da característica intermitente do escoamento bifásico líquido-gás em golfadas a partir do modelo neural, além da estimação de parâmetros médios para as variáveis de interesse com taxas de incerteza variando entre 10% e 16%. / Gas-liquid two-phase flows are present in nature and in different industrial activities. In this type of flow, the liquid and gas phases assume different spatial configurations inside the pipe, called flow patterns. Slug flow is one of the most frequent flow patterns in industrial applications, occurring over a wide range of phase velocities according to studies presented by several authors. The mathematical modelling of slug flow comprises from simple steady state models to more complex models for transient regimes. Those models require closure relationships, e.g. empirical correlations and statistical distributions of characteristic flow parameters. In this work, several models based on artificial neural networks are presented as a support to the characterization of the two-phase slug flow parameters that depend on experimentally obtained void fraction time series. The void fraction time series are measured with a pair of wiremesh sensors installed in a test section of an experimental rig in the premises of the NUEM/UTFPR labs and described in Castillo (2015). From the time series of void fraction measurements relevant parameters to the flow under consideration are computed: the length of the elongated gas bubble, the liquid slug length, the translational velocity of the elongated bubble and the standard deviations for those variables. Those measured and calculated variables are used to obtain a set of artificial neural network-based models. After obtaining such models, a simulation study in which those models are used to estimate the parameters that characterize the two-phase slug flows is carried out. Detailed analysis of the results showed that the variables related to the gas phase are estimated with greater accuracy than the ones related to the liquid phase. As an immediate application of the obtained model, its use as a tool to calculate the initial conditions for a phenomenological mathematical model of twophase slug flow with a slight change of inclination based on a slug tracking method is presented. The differential of this study is to predict the intermittent features of the twophase slug flow by means of a neural model, as well as the estimation of average parameters for the variables of interest with uncertainly rates ranging between 10% and 16%.