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>

The purpose of clarification in the sugar industry is to remove soluble, insoluble and colloidal matter from cane juice. Efficient clarification is required to produce high quality sugar and to prevent entrainment of solids in downstream equipment. The objective of this study is to produce a Computational Fluid Dynamics (CFD) model of the Magra Ultrasep clarifier. This was accomplished by: • Modelling the hydrodynamics of a laboratory scale clarifier In the Fluent CFD program. • Incorporating the flocculation process into the CFD model. • Performing experiments on a pilot scale clarifier to obtain parameter values for the flocculation model. The hydrodynamic model of the clarifier showed the presence of a recirculation zone above the baffle plate. Particle injections using Fluent's discrete phase modelling option determined that particles within the size range of IOOj..Lm to 4mm would circulate in this region, forming the bed of floc particles required for the Magra Ultrasep to work efficiently. The flocculation process in Fluent was represented using three different solid phases of different particle sizes. Small and medium sized particles were allowed to combine to form larger particles by changing the volume fractions according to three rate equations. A fibre glass laboratory scale model was set up at Maidstone Sugar Mill and fed the same sugar cane juice that enters the Rapi-Dorr clarifiers. The experimental results were then fed into a simplified flocculation model in MATLAB. An overall rate constant (k) of 5kg.m-3.s-) for the flocculation kinetic equation satisfied the experimental result. / Thesis (M.Sc.Eng.)-University of KwaZulu-Natal, Durban, 2008.

Anaerobic digestion is a process that simultaneously treats waste and produces renewable energy in the form of biogas. Applications include swine and cattle waste management, which is still dominated by aerobic digestion, a less environmental alternative. The low adoption rates of anaerobic digestion is partly caused by the lack of modelling basis for the technology. This is due to the complexity of the process, as it involves dozens of interrelated biochemical reactions driven by hundreds of species of micro-organisms, immersed in a three-phase, non-Newtonian fluid. As a consequence, no practical computer models exist, and therefore, unlike most other engineering fields, the design process for anaerobic digesters still relies heavily on traditional methods such as trial and error. The current state-of-the-art model is Anaerobic Digestion Model No. 1 (ADM1), published by the International Water Association in 2001. ADM1 is a bulk model, therefore it does not account for the effects of concentration gradients, stagnation regions, and particle settling. To address this, this thesis works toward the creation of the first three-dimensional spatially resolved anaerobic digestion model, called Anaerobic Digestion Model with Multi-Dimensional Architecture (ADM-MDA), by developing a framework. The framework, called Coupled Reaction-Advection Flow Transient Solver (CRAFTS), is a general reaction solver for single-phase, incompressible fluid flows. It is a novel partial differential and algebraic equation (PDAE) solver that also employs a novel programmable logic controller (PLC) emulator, allowing users to define their own control logic. All aspects of the framework are verified for proper function, but still need validation against experimental results. The biochemistry from ADM1 is input into CRAFTS, resulting in a manifestation of ADM-MDA; however the numerical stiffness of ADM1 is found to conflict with the second order accuracy of CRAFTS, and the resulting model can only operate under restricted conditions. Preliminary results show spatial effects predicted by the CRAFTS model, and non-observable in the bulk model, impact the digester in a non-trivial manner and lead to measurable differences in their respective outputs. A detailed discussion of suggested work to arrive at a practical spatially resolved anaerobic digestion model is also provided.