A model has been developed using the general-purpose Navier-Stokes solver CFX4 to simulate Atmospheric Boundary Layer flow over complex terrain. This model has been validated against the measured data from the Askervein Hill experiment, and has been shown to perform well. The CFD model is also compared to the WAsP linear model of wind flow over topography, and a significant improvement is noted for flow over complex topography. Boundary conditions, gridding issues and sensitivity to other solver parameters have all been investigated. An advanced roughness model has been developed to simulate flow over forest canopies, using a resistive body force within the canopy volume. The model is validated against measured data for simple 2D cases, and for a complex 3D case over real topography. The model is shown to give a more physically realistic profile for the wind speed in and just above forest canopies than the standard roughness length model used in most CFD simulations. An automated methodology for setting up CFD simulations using the models described has been developed. A custom pre-processing package to implement this has been written, to enable the use of the CFD methodology in a commercial environment.

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

Interactions between the aerodynamics and dynamics of bluff bodies are important in many engineering applications, including suspension bridges, tall buildings, oil platforms, wind turbine towers, air drops, and construction with cranes. In the rotorcraft field, bluff bodies are commonly suspended underneath the vehicle by tethers. This approach is often the only practical way to deliver a payload in a reasonable amount of time in disaster relief efforts, search-and-rescue operations, and military operations. However, currently a fundamental understanding of the aerodynamics of these bluff bodies is lacking, and accurate dynamic simulation models for predicting the safe flying speed are not available. In order to address these shortcomings, two main advancements are presented in this thesis. The aerodynamics of several three-dimensional canonical bluff bodies are examined over a range of Reynolds numbers representative of wind-tunnel-scale to full-scale models. Numerical experiments are utilized, with a focus on uncertainty analysis and validation of the computations. Mean and unsteady forces and moments for these bluff bodies have been evaluated, and empirical models of the shear layer characteristics have been extracted to quantify the behaviors and provide predictive capability. In addition, a physics-based reduced-order simulation model has been developed for bluff bodies. The physics-based approach is necessary to ensure that the predicted behavior of new configurations is accurate, and it is made possible by the breakthroughs in three-dimensional bluff body aerodynamics presented in this thesis. The integrated aerodynamic forces and moments and dynamic behavior predicted by model are extensively validated with data from wind tunnels, flight tests, and high-fidelity computations. Furthermore, successful stability predictions for tethered loads are demonstrated. The model is applicable to the simulation of any generic bluff body configuration, is readily extensible, and has low computational cost.

The problem of real-time simulation of ocean surface waves, ship movement and the coupling in between is tackled, and a number of different methods are covered and discussed. Among these methods, the finite volume method has been implemented in an attempt to solve the problem, along with the compressible Euler equations, an octree based staggered grid which allows for easy adaptive mesh refinement, the volume of fluid method and a variant of the Hyper-C advection scheme for compressible flows for advection of the phase fraction field. The process of implementing the methods that were chosen proved to be tricky in many ways, as they involve a large number of advanced topics, and the implementation that was implemented in this thesis work suffered from numerous issues. There were for example problems with keeping the interface intact, as well as a harsh restriction on the time step size due to the CFL condition. Improvements required to make the method sustainable for real-time applications are discussed, and a few suggestions on alternative approaches that are already in use for similar purposes are also given and discussed. Furthermore, a method for compensating for gain/loss of mass when solving the incompressible flow equations with an inaccurately solved pressure Poisson equation is presented and discussed. A momentum conservative method for transporting the velocity field on staggered grids without introducing unnecessary smearing is also presented and implemented. A simple, physically based illumination model for sea surfaces is derived, discussed and compared to the Blinn–Phong shading model, although it is never implemented. Finally, a two-dimensional partial differential equation in the spatial domain for simulating water surface waves for mildly varying bottom topography is derived and discussed, although it is deemed to be too slow for real-time purposes and is therefore never implemented. / <p>This publication differs from the printed version of the report in the sense that links are blue in this version and black in the printed version.</p>

Aquesta tesi centra els seus esforços en l'àmbit de la microfluídica, un camp relativament recent dins de la Mecànica de Fluids, amb un futur prometedor i amb un ritme d'investigació intens en les seves diferents especialitzacions. En aquest sentit, la tesi presenta dos aportacions científiques principals. Primer, aporta una eina numèrica d'elaboració pròpia per realitzar simulacions de fluxos reactius en microcanals. Eina que s'aplica satisfactòriament a la identificació dels principals processos de transport involucrats en la oxidació parcial del metà per a produir gas de síntesi, i a l'estudi de l'efecte que tenen alguns paràmetres d'operació en aquest procés reactiu. Segon, estén el coneixement dels fluxos multifàsics en microunions en T, estudiant experimentalment fluxos de dues fases amb fluids principalment miscibles i en condicions supercrítiques, que son portats al seu equilibri vapor-líquid. Durant aquest estudi, a més, reporta un succés inesperat que presenta futurs reptes en l'aplicació d'aquest tipus de fluxos multifàsics. / The present thesis focuses on microfluidics, a relatively recent field of Fluid Mechanics with promising expectations and with an intense scientific interest on its different areas. In this regard, the thesis aims to provide two main scientific contributions. First, it presents an in-house numerical tool to carry out simulations of reactive flows within microchannels. The tool is successfully applied to the identification of the main transport phenomena involved on the partial oxidation of methane to produce synthesis gas, and to the analysis of the effect of several operating parameters on this reactive process. Second, it extends the knowledge on multiphase flows in microfluidic T-junctions with an experimental study of two-phase flows of mixtures of potentially miscible fluids, in supercritical conditions and in vapour-liquid equilibrium. In this study it is also reported an unexpected phenomenon, which brings new challenges to the application of these kind of multiphase flows.

Three models for the generation and propagation of hydrodynamic noise near the ocean surface are presented, and are compared for their ability to predict hull noise generated by CFAV Quest. The simulated fluctuating pressure field on the hull is also validated against experimental results. The near field flow is first solved using the NWT CFD package, and the hydrodynamic noise is calculated using the Lighthill-Curle acoustic analogy. The far field sound is obtained using three methods: a method of images solution to the Lighthill-Curle equations, a simple source model of the transmission loss, or a normal mode model of the transmission loss. Both the simple source and method of images models improve the SPL predictions of the Lighthil-Curle equations. Best performance is obtained from the method of images, improving predictions by approximately 40 dB. The normal mode model is shown to give poor results, due to assumed sea-floor boundary conditions.

A computational and experimental study of the Richtmyer-Meshkov instability is presented here for an inclined interface perturbation. The computational work is composed of simulation studies of the inclined interface RMI performed using the Arbitrary Lagrange Eulerian (ALE) code called ARES. These simulations covered a wide range of Mach numbers (1.2 to 3.5), gas pairs (Atwood numbers 0.23to 0.95), inclination angles (30° to 85°), and explored various perturbation types (both inclined interface and sinusoidal). The computational work included the first parametric study of the inclined interface RMI. This study yielded the first scaling method for the inclined interface RMI mixing width growth rates. It was extended to explore the effect of perturbation linearity and identified that a sharp transition in growth regimes occurs for an initial perturbation inclination angle of 75° with angles below (above) this growing faster (slower). Finally a study of the effects of incident shock strength on the refracted shock wave perturbation decay rate is presented. This study examined how the perturbations induced on the transmitted shock front by the RMI decay with time and found that the decay rates follow a power law model, Alpha=Beta∗S^(Epsilon). When the coefficients from the power law decay model were plotted versus Mach number, a distinct transition region was found which is likely a result of the post-shock heavy gas velocity transitioning from the subsonic to supersonic range. The experimental portion of this work was conducted using the TAMUFMSTF, completed in May of 2012. This facility uses a variable inclination shock tube, with a modular construction design for incident shock strengths of up to Mach 3.0. It employs optical systems for measuring density and velocity fields simultaneously using the planar laser induced fluorescence and particle imaging velocimetry techniques. The design and construction of this facility is reviewed in detail in chapter 4 of this work. The initial experiments performed in the TAMUFMSTF provided the first known extensive experimental data for an inclined interface RMI. Planar laser Mie scattering images and velocity vectors were obtained for a N_(2)/CO_(2) interface at a 60° inclination angle and an incident shock strength of Mach 1.55. These images have been compared with simulations made using the ARES codes and have been shown to have some distinct differences. Some of these differences indicate that the initial conditions in the experiments deviate from the ideal planar interface. Other differences have revealed features which have not been resolved by the simulations due to resolution limitations.

No description available.

No description available.

<p>Nuclear reactors release small amounts of radioactivity during their normal operations. The most common method of calculating the dose to the public that results from such releases uses Gaussian Plume models. We are investigating these methods using CAP88-PC, a computer code developed for the Environmental Protection Agency (EPA) in the USA that calculates the concentration of radionuclides released from a stack using Pasquill stability classification. A buoyant or momentum driven part is also included. The uptake of the released radionuclide by plants, animals and humans, directly and indirectly, is then calculated to obtain the doses to the public. This method is well established but is known to suffer from many approximations and does not give answers that are accurate to be better than 50% in many cases. More accurate, though much more computer-intensive methods have been developed to calculate the movement of gases&nbsp / using fluid dynamic models. Such a model, using the code FLUENT can model complex terrains and will also be investigated in this work. This work is a preliminary study to compare the results of the traditional Gaussian plume model and a fluid dynamic model for a simplified case. The results indicate that Computational Fluid Dynamics calculations give qualitatively similar results with the possibility of including much more effects than the simple Gaussian plume model.</p>