No description available.

No description available.

Hydraulic fracturing is a process for the extraction of hydrocarbons from underground formations. It involves pumping a specialized fluid into the wellbore under high pressures to form and support fractures in the rock. Fracturing stimulates the well to increase the production of oil and the natural gas which are the pillars of the energy economy. Key to this process is the use of proppants, which are solid materials used to keep the fractures open. Understanding the transport of proppant particles through a fluid is important to improve the efficiency and reduce environmental impact of fracturing. An increase of the settling velocity for instance, will impede the hydraulic fracturing process by reducing well productivity, or necessitate use of chemical additives. This thesis presents a theoretical investigation of the settling velocity of proppant particles. The effect of different parameters on the settling velocity were studied by manipulating the main factors that can influence particle transport. These include size of the particle (300 μm- 2000 μm), sphericity, density (1200 kg/m3-3500 kg/m3) and concentration. These typical values were obtained from commercially available proppants currently used in industry. Various correlations were investigated, assuming the carrier (fracturing) fluid to be an ideal Newtonian and as a power law (non-Newtonian) fluid. This will help predict the settling velocity for proppant particles in order to increase well productivity, and improve hydraulic fracturing efficiency. The models show that changing the carrier fluid viscosity and particle properties such as diameter, density, sphericity, and concentration leads to a significant change in the proppant settling velocity. For instance, reduction in particle size, density, and sphericity tend to reduce the settling velocity, while increasing the concentration of the particles and the fluid viscosity reduce the settling velocity.

With the depletion of oil reserves and increase in oil price, the enhanced oil recovery methods such as polymer flooding to increase oil production from water flooded fields are becoming more attractive. Effective design of these processes is challenging because the polymer chemistry has a strong effect on reaction and fluid rheology, which in turn has a strong effect on fluid transport. We have implemented a well-established polymer model within the Implicit Parallel Accurate Reservoir Simulator (IPARS), which enables parallel simulation of non-Newtonian fluid flow through porous media. The following properties of polymer solution are modeled in this work: 1) polymer adsorption; 2) polymer viscosity as a function of salinity, hardness, polymer concentration, and shear rate; 3) permeability reduction; 4) inaccessible pore volume. IPARS enables field-scale polymer flooding simulation with its parallel computation capability. In this thesis, several numerical examples are presented. The result of polymer module is verified by UTCHEM, a three-dimensional chemical flood simulator developed at the University of Texas at Austin. The parallel capability is also tested. The influence of different shear rate calculations is investigated in homogeneous and heterogeneous reservoirs. We observed that the wellbore velocity calculation instead of Darcy velocity reduces the grid effect for coarse mesh. We noted that the injection bottom hole pressure is very sensitive to the shear rate calculation. However, cumulative oil recovery and overall oil saturation appear to not be sensitive to grid and shear rate calculation for same reservoir. There are two models to model the ground water infiltration in vadose zone. One is Richard’s Equation (RE) model. And the other is two-phase flow model. In this work, we compare the two-phase model with an RE model to ascertain, under common scenarios such as infiltration or injection of water into initially dry soils, the similarities and differences in solutions behaviors, the ability of each model to simulate such infiltration processes under realistic scenarios, and to investigate the numerical efficiencies and difficulties which arise in these models. Six different data sets were assembled as benchmark infiltration problems in the unsaturated zone. The comparison shows that two-phase model holds for general porous media and is not limited by several assumptions that must be made for the RE formulation, while RE is applicable only for shallow regions (vadose) that are only several meters in depth and a fully saturated bottom boundary condition must be assumed. / text

Even though fluid mechanics is well developed as a science, there are many physical phenomena that we do not yet fully understand. One of these is the deformation rates and fluid stresses generated in a boundary layer for a non-Newtonian fluid. One such non-Newtonian fluid would be a waxy crude oil flowing in a centrifugal pump. This type of flow can be numerically modeled by a rotating disk system, in combination with an appropriate constitutive equation, such as the relation for a Bingham fluid. A Bingham fluid does not begin to flow until the stress magnitude exceeds the yield stress. However, experimental measurements are also required to serve as a database against which the results of the numerical simulation can be interpreted and validated. The purpose of the present research is to gain a better understanding of the behavior of a Bingham fluid in the laminar boundary layer on a rotating disk. For this project, two different techniques were employed: numerical simulation, and laboratory investigations using Particle Image Velocimetry (PIV) and flow visualization. Both methods were applied to the flow of a Bingham fluid over a rotating disk. In the numerical investigations, the flow was characterized by the dimensionless yield stress Bingham number, By, which is the ratio of the yield and viscous stresses. Using von Kármáns similarity transformation, and introducing the rheological behavior of the fluid into the conservation equations, the corresponding nonlinear two-point boundary value problem was formulated. A solution to the problem under investigation was obtained by numerical integration of the set of Ordinary Differential Equations (ODEs) using a multiple shooting method. The influence of the Bingham number on the flow behavior was identified. It decreases the magnitude of the radial and axial velocity components, and increases the magnitude of the tangential velocity component, which has a pronounced effect on the moment coefficient, CM, and the volume flow rate, Q. In the laboratory investigations, since the waxy crude oils are naturally opaque, an ambitious experimental plan to create a transparent oil that was rheologically similar to the Amna waxy crude oil from Libya was developed. The simulant was used for flow visualization experiments, where a transparent fluid was required. To fulfill the demand of the PIV system for a higher degree of visibility, a second Bingham fluid was created and rheologically investigated. The PIV measurements were carried out for both filtered tap water and the Bingham fluid in the same rotating disk apparatus that was used for the flow visualization experiments. Both the axial and radial velocity components in the (r-z) plane were measured for various rotational speeds. Comparison between the numerical and experimental results for the axial and radial velocity profiles for water was found to be satisfactory. Significant discrepancies were found between numerical results and measured values for the Bingham fluid, especially at low rotational speeds, mostly relating to the formation of a yield surface within the tank. Even though the flow in a pump is in some ways different from that of a disk rotating in a tank, some insight about the behavior of the pump flow can be drawn. One conclusion is that the key difference between the flow of a Bingham fluid in rotating equipment from that of a Newtonian fluid such as water relates to the yield surface introduced by the yield stress of the material, which causes an adverse effect on the performance and efficiency of such equipment.

Tissue engineering is a multidisciplinary subject aimed at producing the immunologically tolerant artificial tissues/organs to repair or replace damaged ones. In this field, tissue scaffold plays a key role to support cell growth and new tissue regeneration. For fabrication of tissue scaffolds with individual external geometry and predefined inner structure, rapid prototyping (RP) systems based on fluid dispensing techniques have proved to be very promising. The present research conducted a comprehensive study on the dispensing-based fabrication process. <p>First of all, the scaffold materials are characterized in terms of their biocompatibility and flow behaviour. The biocompatibility of biomaterials of PLLA, PCL, collagen, chitosan, and gelatine is evaluated in terms of supporting neuron cells adhesion and outgrowth. Chitosan solution (2% w/v) in acetic acid is shown to be the most promising among the examined biomaterials for the fabrication of nerve tissue scaffolds. Its non-Newtonian flow behaviour is identified by using a commercial rheometer. <p>In the fabrication process, the flow rate of biomaterials dispensed, the profile of strand cross-sections, and the scaffold porosity are very important and must be precisely controlled. A model is developed to represent the flow rate of biomaterials dispensed under the assumptions that the flow is incompressible, steady, laminar, and axisymmetric. Also, the profile and size of line strands at different layers and portions are modeled based on the Young-Laplace equation. Thus the dispensing-based fabrication process can be predicted in terms of the flow rate and the scaffold porosity. <p>The effects of operation conditions on the fabrication result are identified theoretically and experimentally. Simulation result shows that a higher driving pressure, a higher temperature, and a larger needle diameter will result in a larger size of the strand cross-sections and lower scaffold porosity. The change pattern, however, is nonlinear, which is affected by the fluid surface tension and non-Newtonian flow behaviour of scaffold biomaterials. <p>To verify the effectiveness of the developed models, experiments were carried out on a commercial dispensing system (C-720, Asymtek, USA). To avoid the possible error derived from the temperature difference between the dispensing system and the rheometer, a new method is presented to characterize the fluid properties used for model predictions. Experimental results illustrate that the developed models, combined with the new identification method, are very promising to predict the dispensing-based fabrication process.

Even though fluid mechanics is well developed as a science, there are many physical phenomena that we do not yet fully understand. One of these is the deformation rates and fluid stresses generated in a boundary layer for a non-Newtonian fluid. One such non-Newtonian fluid would be a waxy crude oil flowing in a centrifugal pump. This type of flow can be numerically modeled by a rotating disk system, in combination with an appropriate constitutive equation, such as the relation for a Bingham fluid. A Bingham fluid does not begin to flow until the stress magnitude exceeds the yield stress. However, experimental measurements are also required to serve as a database against which the results of the numerical simulation can be interpreted and validated. The purpose of the present research is to gain a better understanding of the behavior of a Bingham fluid in the laminar boundary layer on a rotating disk. For this project, two different techniques were employed: numerical simulation, and laboratory investigations using Particle Image Velocimetry (PIV) and flow visualization. Both methods were applied to the flow of a Bingham fluid over a rotating disk. In the numerical investigations, the flow was characterized by the dimensionless yield stress Bingham number, By, which is the ratio of the yield and viscous stresses. Using von Kármáns similarity transformation, and introducing the rheological behavior of the fluid into the conservation equations, the corresponding nonlinear two-point boundary value problem was formulated. A solution to the problem under investigation was obtained by numerical integration of the set of Ordinary Differential Equations (ODEs) using a multiple shooting method. The influence of the Bingham number on the flow behavior was identified. It decreases the magnitude of the radial and axial velocity components, and increases the magnitude of the tangential velocity component, which has a pronounced effect on the moment coefficient, CM, and the volume flow rate, Q. In the laboratory investigations, since the waxy crude oils are naturally opaque, an ambitious experimental plan to create a transparent oil that was rheologically similar to the Amna waxy crude oil from Libya was developed. The simulant was used for flow visualization experiments, where a transparent fluid was required. To fulfill the demand of the PIV system for a higher degree of visibility, a second Bingham fluid was created and rheologically investigated. The PIV measurements were carried out for both filtered tap water and the Bingham fluid in the same rotating disk apparatus that was used for the flow visualization experiments. Both the axial and radial velocity components in the (r-z) plane were measured for various rotational speeds. Comparison between the numerical and experimental results for the axial and radial velocity profiles for water was found to be satisfactory. Significant discrepancies were found between numerical results and measured values for the Bingham fluid, especially at low rotational speeds, mostly relating to the formation of a yield surface within the tank. Even though the flow in a pump is in some ways different from that of a disk rotating in a tank, some insight about the behavior of the pump flow can be drawn. One conclusion is that the key difference between the flow of a Bingham fluid in rotating equipment from that of a Newtonian fluid such as water relates to the yield surface introduced by the yield stress of the material, which causes an adverse effect on the performance and efficiency of such equipment.

Tissue engineering is a multidisciplinary subject aimed at producing the immunologically tolerant artificial tissues/organs to repair or replace damaged ones. In this field, tissue scaffold plays a key role to support cell growth and new tissue regeneration. For fabrication of tissue scaffolds with individual external geometry and predefined inner structure, rapid prototyping (RP) systems based on fluid dispensing techniques have proved to be very promising. The present research conducted a comprehensive study on the dispensing-based fabrication process. <p>First of all, the scaffold materials are characterized in terms of their biocompatibility and flow behaviour. The biocompatibility of biomaterials of PLLA, PCL, collagen, chitosan, and gelatine is evaluated in terms of supporting neuron cells adhesion and outgrowth. Chitosan solution (2% w/v) in acetic acid is shown to be the most promising among the examined biomaterials for the fabrication of nerve tissue scaffolds. Its non-Newtonian flow behaviour is identified by using a commercial rheometer. <p>In the fabrication process, the flow rate of biomaterials dispensed, the profile of strand cross-sections, and the scaffold porosity are very important and must be precisely controlled. A model is developed to represent the flow rate of biomaterials dispensed under the assumptions that the flow is incompressible, steady, laminar, and axisymmetric. Also, the profile and size of line strands at different layers and portions are modeled based on the Young-Laplace equation. Thus the dispensing-based fabrication process can be predicted in terms of the flow rate and the scaffold porosity. <p>The effects of operation conditions on the fabrication result are identified theoretically and experimentally. Simulation result shows that a higher driving pressure, a higher temperature, and a larger needle diameter will result in a larger size of the strand cross-sections and lower scaffold porosity. The change pattern, however, is nonlinear, which is affected by the fluid surface tension and non-Newtonian flow behaviour of scaffold biomaterials. <p>To verify the effectiveness of the developed models, experiments were carried out on a commercial dispensing system (C-720, Asymtek, USA). To avoid the possible error derived from the temperature difference between the dispensing system and the rheometer, a new method is presented to characterize the fluid properties used for model predictions. Experimental results illustrate that the developed models, combined with the new identification method, are very promising to predict the dispensing-based fabrication process.

The depletion of oil reserves and higher oil prices has made chemical enhanced oil recovery (EOR) methods more attractive in recent years. Because of geological heterogeneity, unfavorable mobility ratio, and capillary forces, conventional oil recovery (including water flooding) leaves behind much oil in reservoir, often as much as 70% OOIP (original oil in place). Surfactant/polymer flooding targets these bypassed oil left after waterflood by reducing water mobility and oil/water interfacial tension. The complexity and uncertainty of reservoir characterization make the design and implementation of a robust and effective surfactant/polymer flooding to be quite challenging. Accurate numerical simulation prior to the field surfactant/polymer flooding is essential for a successful design and implementation of surfactant/polymer flooding. A recently developed unified polymer viscosity model was implemented into our existing polymer module within our in-house reservoir simulator, the Implicit Parallel Accurate Reservoir Simulator (IPARS). The new viscosity model is capable of simulating not only the Newtonian and shear-thinning rheology of polymer solution but also the shear-thickening behavior, which may occur near the wellbore with high injection rates when high molecular weight Partially Hydrolyzed Acrylamide (HPAM) polymers are injected. We have added a full capability of surfactant/polymer flooding to TRCHEM module of IPARS using a simplified but mechanistic and user-friendly approach for modeling surfactant/water/oil phase behavior. The features of surfactant module include: 1) surfactant component transport in porous media; 2) surfactant adsorption on the rock; 3) surfactant/oil/water phase behavior transitioned with salinity of Type II(-), Type III, and Type II(+) phase behaviors; 4) compositional microemulsion phase viscosity correlation and 5) relative permeabilities based on the trapping number. With the parallel capability of IPARS, commercial scale simulation of surfactant/polymer flooding becomes practical and affordable. Several numerical examples are presented in this dissertation. The results of surfactant/polymer flood are verified by comparing with the results obtained from UTCHEM, a three-dimensional chemical flood simulator developed at the University of Texas at Austin. The parallel capability and scalability are also demonstrated. / text

Yes / A novel Euler-Lagrangian (EL) computational uid dynamics (CFD) nite volume-based model to simulate the gas mixing of sludge for anaerobic digestion is developed and described. Fluid motion is driven by momentum transfer from bubbles to liquid. Model validation is undertaken by assessing the ow eld in a labscale model with particle image velocimetry (PIV). Conclusions are drawn about the upscaling and applicability of the model to full-scale problems, and recommendations are given for optimum application.