During injection of a high-viscosity, non-Newtonian polymer into a long horizontal well, a significant pressure drop occurs along the well length. Computational Fluid Dynamics (CFD) modeling of the shear-thinning flow of polymer in the wellbore, coupled with the viscoelastic flow in composite gravel-pack/near-well formation zone, was carried out to develop convenient correlations for axial pressure values of both Newtonian and non-Newtonian fluids along the well length, for use in chemical EOR simulations. The detailed CFD modeling of the non-Newtonian flow behavior of polymer within the horizontal wellbore, completion zone and the near-well formation, not only allows accurate accounting of pressure distribution along the long horizontal well, but also can be employed for screening diagnosis for possible injectivity inefficiencies resulting from non-uniform pressure values. At both high and low injection rates, CFD modeling predicts non-uniform pressure distributions for highly viscous fluids. The inclusive pressure correlation was implemented into UTCHEM, a University of Texas at Austin research simulator, to determine the importance of including pressure drop in polymer injections. Early times (i.e., less than 100 days) yielded a significant oil recovery deviation from a uniform pressure wellbore. However, at later times the recovery loss generated by the pressure decrease was deemed negligible; therefore, the traditional assumption regarding uniform pressure in horizontal wellbores was still reasonable for highly viscous non-Newtonian flow. This CFD study is the first mechanistic investigation of the polymer injectivity with detailed description of the wellbore, completion zone and near-well formation, and with full accounting of the shear-thinning rheology for pipe flow and the viscoelastic rheology of polymer in porous media. With increased use of very high molecular-weight polymers for chemical EOR processes for mobility control, the latter mechanism is known to be critical. / text

This thesis develops methods to identify periodic solutions to the n-body problem by representing gravitational orbits with Fourier series. To find periodic orbits, a minimization function was developed that compares the second derivative of the Fourier series with Newtonian gravitation acceleration and modifies the Fourier coefficients until the orbits match. Software was developed to minimize the function and identify the orbits using gradient descent and quadratic curves. A Newtonian gravitational simulator was developed to read the initial orbit data and numerically simulate the orbits with accurate motion integration, allowing for comparison to the Fourier series orbits and investigation of their stability. The orbits found with the programs correlate with orbits from literature, and a number remain stable when simulated. / February 2016

The existence of Gravitational Waves from binary black holes is one of the most interesting predictions of General Relativity. These ripples in space-time should be visible to ground-based gravitational wave detectors worldwide in the next few years. One such detector, the Laser Interferometer Gravitational-wave Observatory (LIGO) is in the process of being upgraded to its Advanced sensitivity which should make gravitational wave detections routine. Even so, the signals that LIGO will detect will be faint compared to the detector noise, and so accurate waveform templates are crucial. In this thesis, we present a detailed analysis of the accuracy of hybrid gravitational waveforms. Hybrids are created by stitching a long post-Newtonian inspiral to the late inspiral, merger, and ringdown produced by numerical relativity simulations. We begin our investigation with a study of the systematic errors in the numerical waveform, and errors due to hybridization and choice of detector noise. For current NR waveforms, the largest source of error comes from the unknown high-order terms in the post-Newtonian waveform, which we first explore for equal-mass, non-spinning binaries, and also for unequal-mass, non-spinning binaries. We then consider the potential reduction in hybrid errors if these higher-order terms were known. Finally, we investigate the possibility of using hybrid waveforms as a detection template bank and integrating NR+PN hybrids into the LIGO detection pipeline.

This thesis reports on measurements of drag on an unbounded cylinder in a viscoelastic Stokes flow, and compares these values with a Newtonian equivalent. Cylinders of diameter 0.5 to 3.34 mm were submerged 10 to 36 mm into slowly rotating annular tanks with channel widths between 133 to 152 mm. Theoretical formulas and computer simulations were used to correct for the effects of ends and walls, yielding estimates of the unbounded drag. The methodology was verified by testing Newtonian fluids and comparing the results to Kaplun’s (1957) prediction for unbounded drag. The test fluids used were a silicone oil, a polybutene, and two Boger fluids. By comparing the Boger fluid results to equally viscous Newtonian results, the contributions of elasticity to the drag were determined. The Deborah number (De) was used to represent the magnitude of flow elasticity, and an onset of elastic effects was measured between 0.5 and 0.7.

This thesis reports on measurements of drag on an unbounded cylinder in a viscoelastic Stokes flow, and compares these values with a Newtonian equivalent. Cylinders of diameter 0.5 to 3.34 mm were submerged 10 to 36 mm into slowly rotating annular tanks with channel widths between 133 to 152 mm. Theoretical formulas and computer simulations were used to correct for the effects of ends and walls, yielding estimates of the unbounded drag. The methodology was verified by testing Newtonian fluids and comparing the results to Kaplun’s (1957) prediction for unbounded drag. The test fluids used were a silicone oil, a polybutene, and two Boger fluids. By comparing the Boger fluid results to equally viscous Newtonian results, the contributions of elasticity to the drag were determined. The Deborah number (De) was used to represent the magnitude of flow elasticity, and an onset of elastic effects was measured between 0.5 and 0.7.

The existence of Gravitational Waves from binary black holes is one of the most interesting predictions of General Relativity. These ripples in space-time should be visible to ground-based gravitational wave detectors worldwide in the next few years. One such detector, the Laser Interferometer Gravitational-wave Observatory (LIGO) is in the process of being upgraded to its Advanced sensitivity which should make gravitational wave detections routine. Even so, the signals that LIGO will detect will be faint compared to the detector noise, and so accurate waveform templates are crucial. In this thesis, we present a detailed analysis of the accuracy of hybrid gravitational waveforms. Hybrids are created by stitching a long post-Newtonian inspiral to the late inspiral, merger, and ringdown produced by numerical relativity simulations. We begin our investigation with a study of the systematic errors in the numerical waveform, and errors due to hybridization and choice of detector noise. For current NR waveforms, the largest source of error comes from the unknown high-order terms in the post-Newtonian waveform, which we first explore for equal-mass, non-spinning binaries, and also for unequal-mass, non-spinning binaries. We then consider the potential reduction in hybrid errors if these higher-order terms were known. Finally, we investigate the possibility of using hybrid waveforms as a detection template bank and integrating NR+PN hybrids into the LIGO detection pipeline.

This thesis examines how the non-Newtonian characteristics of liquid hog manure affect the flow conditions within a steady-flow anaerobic digester. There are three main parts to this thesis. In the first part of this thesis, the physical properties of liquid hog manure and their variation with temperature and solids concentration are experimentally determined. Naturally¬¬-settled manure sampled from an outdoor storage lagoon is studied, and density, viscosity, and particle size distribution are measured. Hog manure with total solids concentrations of less than 3.6% exhibits Newtonian behaviour; manure between 3.6% and 6.5% total solids is pseudoplastic, and fits the power law; manure with more than 6.5% total solids exhibits non-Newtonian and time-dependent characteristics. The second part of this thesis investigates the flow of Newtonian and non-Newtonian fluids—represented by tap water and xanthan gum solution, respectively—within four lab-scale reactor geometries, using residence time distribution (RTD) experiments. The effect of reactor geometry, flow rate, and fluid viscosity are evaluated. In the third part of this thesis, flow conditions within lab-scale and pilot-scale anaerobic digester reactors are simulated using three-dimensional modeling techniques. The RTDs of lab-scale reactors as predicted by the 3D numerical models compare well to the experimental results. The 3D models are also validated using data from particle image velocimetry (PIV) experiments. Finally, the viscous properties of liquid hog manure at 3% and 8% total solids are incorporated into the models, and the results are evaluated.

A research program was designed to investigate the characteristics of three-dimensional laminar wall jet flow of both Newtonian and two shear-thinning non-Newtonian fluids. The non-Newtonian fluids were prepared from xanthan gum solutions of various concentrations. Both experimental and numerical methodologies were employed in this study. The wall jet was created using a circular pipe of diameter 7 mm and flows into an open fluid tank. The initial Reynolds numbers based on the pipe diameter and jet exit velocity ranged from 250 to 800. The velocity measurements were conducted using a particle image velocimetry technique. The measurements were conducted at several streamwise locations to cover both the developing and self-similar regions. For the numerical study, the complete nonlinear Navier-Stokes equation was solved using an in-house colocated finite volume based CFD code. A Carreau model was employed for the non-Newtonian fluids. The viscosity in the governing equations was obtained explicitly. From the PIV measurements and CFD results, velocity profiles and jet half-widths were extracted at selected downstream locations to study the effects of Reynolds number and specific fluid type on the jet characteristics. It was observed that the numerical results are in reasonable agreement with the experimental data. The decay of maximum velocity, jet spread rates, skin friction coefficient, streamwise velocity profiles, and secondary flows depend strongly on the initial Reynolds number irrespective of the fluid. The results also show that the jet spreads more in the spanwise direction than in the transverse direction in the early flow development whereas the reverse is true in the downstream region. Important differences were observed when the results for the non-Newtonian fluids were compared with those for Newtonian fluid.

The blood is a complex fluid that contains, in addition to water, cells, macro-molecules and a large number of smaller molecules. The physical properties of the blood are therefore the result of non-linear interactions of its constituents, which are influenced by the local flow field conditions. Hence, the local blood viscosity is a function of the local concentration of the blood constituents and the local flow field itself. This study considers the flow of blood-like fluids in generalised 90-degree bifurcating pipes and patient-specific arterial bifurcations relevant to the large aortic branches in humans. It is shown that the Red Blood Cell (RBC) distribution in the region of bifurcations may lead to large changes in the viscosity, with implications on the concentrations of the various cells in the blood plasma. This in turn implies that the flow in the near wall regions is more difficult to estimate and predict than that under the assumption of a homogeneous fluid. The rheological properties of blood are complex and are difficult to measure, since the results depend on the measuring equipment and the inherent flow conditions. We attempt to model the viscosity of water containing different volume fractions of non-deforming RBC-like particles in tubes. The apparent viscosities of the mixtures obtained from these model experiments have been compared to the predictions of the different rheological models found in the literature. The same rheological models have also been used in the different simulations, where the local RBC concentration and local shear rate are used in the viscosity models. The flow simulations account for the non-linearity due to coupling between the flow and fluid rheology. Furthermore, from a physiological perspective, it is shown that oscillatory wall shear stresses are affected by changes in RBC concentration in the regions of the bifurcation associated with atherogenesis. The intrinsic shear thinning rheological property of the blood, in conjunction with stagnation in separated flows, may be responsible for elevated temporal wall shear stress gradients (TWSSG) influencing endothelial cell behaviour, which has been postulated to play a role in the development of atherosclerosis. The blood-like fluid properties along with variations in the RBC concentration could also lead to variations in the developing flow structures in the larger arteries that could influence the work the heart has to bear. / <p>QC 20131206</p>

The goal of this thesis is to study the ground or metastable state structure of nematic liquid crystal systems confined inside handled shapes such as a torus or double torus. We begin our work by introducing a new method to generate a toroidal droplet from a Newtonian liquid inside another, immiscible, Newtonian liquid. In this situation, a toroidal droplet is unstable and follows one of two routes in transforming into a spherical droplet: (i) its tube breaks in a way reminiscent to the breakup of a cylindrical jet, or (ii) its tube grows until it finally coalesces onto itself. However, to be able to probe the nematic structure, we need to address the issue of instabilities. This is done by replacing the outer liquid with a yield stress material, which ultimately leads to the stabilization of the toroidal droplet. Through the experimental investigation, we are able to establish the stabilization conditions. Finally, we generate and stabilize toroidal droplets with a nematic liquid crystal as the inner liquid and a yield stress material as the outer medium. Here we observe that in the ground state, the nematic liquid crystal exhibits an intriguing twisted structure irrespective of the aspect ratio of the torus. While there are no defects observed in a toroidal droplet case, two defects with -1 topological charge each emerge each time we increase the number of handles.