This thesis describes a method for using physically based techniques to model an explosion and the resulting side effects. Explosions are some of the most visually exciting phenomena known to humankind and have become nearly ubiquitous in action films. A realistic computer simulation of this powerful event would be cheaper, quicker, and much less complicated than safely creating the real thing. The immense energy released by a detonation creates a discontinuous localized increase in pressure and temperature. Physicists and engineers have shown that the dissipation of this concentration of energy, which creates all the visible effects, adheres closely to the compressible Navier-Stokes equation. This program models the most noticeable of these results. In order to simulate the pressure and temperature changes in the environment, a three dimensional grid is placed throughout the area around the detonation and a discretized version of the Navier-Stokes equation is applied to the resulting voxels. Objects in the scene are represented as rigid bodies that are animated by the forces created by varying pressure on their hulls. Fireballs, perhaps the most awe-inspiring side effects of an explosion, are simulated using massless particles that flow out from the center of the blast and follow the currents created by the dissipating pressure. The results can then be brought into Maya for evaluation and tweaking.

Conventional one-point turbulence closures have beenextended with an additional transported scalar for modeling ofmagnetohydrodynamic (MHD) turbulence. The new scalar, α ,captures the length scale anisotropy and tendency towardstwo-dimensionality, which is characteristic feature of MHDturbulence, and allows accurate modeling of the Jouledissipation of turbulence. The concept has been used for both afull Reynolds stress closure, and a three-equationK-ε -αmodel. An exact transport equation forαwas derived from the governing equations. All terms inthe equation require modeling, however. The proposed modeltransport equation for α includes terms for magneticdissipation, nonlinear energy transfer, and effects of meanshear and strain. Modeling of the magnetic and strain-relatedterms was based on rapid distortion analysis of the linearizedequations, while modeling of nonlinear effects isphenomenological in nature. For homogeneous turbulence, themodel was compared with linear theory, direct numericalsimulations and experiments. For turbulence subjected to astrong magnetic field, the model reproduces the energy andlength scale evolution predicted by linear theory. Whennonlinear effects are of importance, it predicts energy decayand length scale evolution in agreement with experiments. Theeddy viscosity and Reynolds stress versions of the modelcoincide with the respective conventional models in the absenceof a magnetic field. The objective of this project has been todevelop efficient MHD turbulence models for engineeringapplications, especially for modeling of continuous steelcasting. The novel MHD turbulence models appear to benumerically robust, and they have been implemented in acommercial flow solver, together with electromagnetic equationsfor the Lorentz forces in the mean momentum equations. <b>Keywords:</b>Turbulence model, magnetohydrodynamics, MHD,magnetohydrodynamic turbulence, computational fluid dynamics,continuous casting, dimensionality, Reynolds stresses, eddyviscosity

No description available.

In this thesis a set of boundary conditions for multiphase flow is suggested. Characteristic-based boundary conditions are reviewed for single-phase flow. The problem of open-boundary conditions is investigated, and to avoid drifting values, the use of control functions is proposed. The use of control functions is also verified with a new test which assesses the quality of the boundary conditions. Particularly, P- and PI-control functions are examined. PI-controllers have the ability to specify a given variable exactly at the outlet as well as at the inlet, without causing spurious reflections which are amplified. Averaged multiphase flow equations are reviewed, and a simplified model is established. This model is used for the boundary analysis and the computations. Due to the averaging procedure, signal speeds are reduced to the order of the flow speed. This leads to numerical challenges. For a horizontal channel flow, a splitting of the interface pressure model is suggested. This bypasses the numerical problems associated with separation by gravity, and a physical realistic model is used. In this case, the inviscid model is shown to possess complex eigenvalues, and still the characteristic boundary conditions give reasonable results. The governing equations are solved with a Runge-Kutta scheme for the time integration. For the spatial discretisation, a finite-volume and a finite-difference method are used. Both implementations give equivalent results. In single-phase flow, the results improve significantly when a numerical filter is applied. For two-dimensional two-phase flow, the computations are unstable without a numerical filter.

In this thesis, the effects of split canard on the aerodynamic characteristics of missiles at high angles of attack are numerically investigated. Moreover, an enhanced semi-empirical engineering-level method is developed for prediction of normal force and pitching moment of split canard mounted missiles. In order to analyze the effects of split canard, a generic test case model is created by mounting a split canard to a generic test case model, NASA Dual Control Missile (NDCM), which was previously modeled and analyzed for the validation of CFD modeling. After obtaining a generic missile model with split canard, the effects of split canard on the aerodynamic characteristics of this missile in case of no control, pitch control, yaw control, and roll control deflections are numerically investigated. It is seen that the split canard decreases the local angle of attack of existing canard, increases the normal force and the maneuverability of the missile, and reduces the induced rolling moment at high angles of attack. Five different aerodynamic design parameters are determined for split canard and the effects of each parameter on missile aerodynamics are numerically investigated. It is seen that the roll orientation, deflection angle, size of the split canards have strong effects on missile&rsquo / s aerodynamic performance whereas longitudinal position of the split canards only affects the pitching moment of the missile. Finally, an enhanced semi-empirical engineering-level method, CFD-CBU, is developed for split canard mounted missiles in order to predict the normal force and the pitching moment coefficients. The developed method is validated with NDCM test case model. After this validation, the method is applied to the split canard mounted generic missile in case of no control deflection and pitch control deflection. The results of this method are compared with CFD results and it is seen that the results are in good agreement with each other.

This work presents a number of the practical aspects of developing reduced- order models (ROMs) based on proper orthogonal decomposition (POD). ROMS are derived and implemented for multiphase flow, quasi-2D nozzle flow and 2D inviscid channel flow. Results are presented verifying the ROMs against existing full-order models (FOM). POD is a method for separating snapshots of a flow field that varies in both time and space into spatial basis functions and time coefficients. The partial differential equations that govern fluid flow can then be pro jected onto these basis functions, generating a system of ordinary differential equations where the unknowns are the time coefficients. This results in the reduction of the number of equations to be solved from hundreds of thousands or more to hundreds or less. A ROM is implemented for three-dimensional and non-isothermal multiphase flows. The derivation of the ROM is presented. Results are compared against the FOM and show that the ROM agrees with the FOM. While implementing the ROM for multiphase flow, moving discontinuities were found to be a ma jor challenge when they appeared in the void fraction around gas bubbles. A point-mode POD approach is proposed and shown to have promise. A simple test case for moving discontinuities, the first order wave equation, is used to test an augmentation method for capturing the discontinuity exactly. This approach is shown to remove the unphysical oscillations that appear around the discontinuityin traditional approaches. A ROM for quasi-2D inviscid nozzle flow is constructed and the results are com- pared to a FOM. This ROM is used to test two approaches, POD-Analytical and POD-Discretized. The stability of each approach is assessed and the results are used in the implementation of a ROM for the Navier-Stokes equations. A ROM for a Navier-Stokes solver is derived and implemented using the results of the nozzle flow case. Results are compared to the FOM for channel flow with a bump. The computational speed-up of the ROM is discussed. Two studies are presented with practical aspects of the implementation of POD- based ROMs. The first shows the effect of the snapshot sampling on the accuracy of the POD basis functions. The second shows that for multiphase flow, the cross- coupling between field variables should not be included when computing the POD basis functions.

QC 20110118

The topic of this thesis concerns the LHC, the next large particle accelerator at CERN which will start operating in 2007. Being based on superconductivity, the LHC needs to operate at very low temperatures, which makes great demands on the cryogenic system of the accelerator. To cope with the heat loads induced by the particle beam, a beam screen cooled with forced flow of supercritical helium is used. There is an interest in upgrading the energy and luminosity of the LHC in the future and this would require a higher heat load to be extracted by the beam screen cooling system. The objective of this thesis is to quantify different ways to upgrade this system by mainly studying the effects of different pressure and temperatures levels as well as a different cooling medium, neon. For this a numerical program which simulates one-dimensional pipe flow was constructed. The frictional forces were accounted for by the empirical concept of friction factor. For the fluid properties, software using empirically made correlations was used. To validate the numerical program, a comparison with previous experimental work was done. The agreement with experimental data was good for certain flow configurations, worse for others. From this it was concluded that further comparisons with experimental data must be made in order to tell the accuracy of the mathematical model and the correlations for fluid properties used. When using supercritical helium, thermo-hydraulic instabilities may arise in the cooling loop. It was of special interest to see how well a numerical program could simulate and predict this phenomenon. It was found that the numerical program did not function for such unstable conditions; in fact it was much more sensitive than what reality is. For the beam screen cooling system we conclude that to cope with the increased heat loads of future upgrades, an increase in pressure level is needed regardless if the coolant remains helium, or is changed to neon. Increasing the pressure level also makes that the problems with thermo-hydraulic instabilities can be avoided. Of the two coolants, helium gave the best heat extraction capacity. Unlike neon, it is also possible to keep the present temperature level when using helium.

<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>

A spur dike is a hydraulic structure built on the bank of a river at some angle to the main flow direction. A series of spur dikes in a row may also be placed on one side or both sides of a river to form a spur dike field. Spur dikes are used for two main purposes, namely river training and bank protection. For river training, spur dikes may be used to provide a desirable path for navigation purposes or to direct the flow to a desirable point such as a water intake. For bank protection, spur dikes may be used to deflect flow away from a riverbank and thus protect it from erosion. It has also been observed that spur dikes provide a desirable environment for aquatic habitat. Despite the fact that spur dikes are useful hydraulic structures, they have been found to increase the flow resistance in rivers and hence increase the flow stage. The present study focuses on the quantification of the flow resistance and associated flow stage increase due to a single spur dike and also that of a spur dike field. Increased flow stage is referred to herein as a backwater effect.<p> In the first stage of the study, the flow resistance due to a single spur dike, expressed as a drag force exerted on the flow in an open channel, was studied and quantified. The work was carried out in a rigid bed flume, with the model spur dike being simulated using various sizes of a two-dimensional (2-D) rectangular plate. Several discharge conditions were studied. The drag force exerted by the spur dike for both submerged and unsubmerged flow conditions was determined directly from measurements made using a specially designed apparatus and also by application of the momentum equation to a control volume that included the spur dike. It was found that the unit drag force (i.e., drag force per unit area of dike) of an unsubmerged spur dike increases more rapidly with an increase in the discharge when compared with that of a submerged spur dike. The results also showed that an increase in the blockage of the open channel cross-section due to the spur dike is the main parameter responsible for an increase in the spur dike drag coefficient, hence the associated flow resistance and backwater effect. Based on these findings, relationships were developed for estimating the backwater effect due to a single spur dike in an open channel.<p> In the second stage of the study, the flow resistance due to a spur dike field expressed as a drag force exerted on the flow was quantified and subsequently related to the backwater effect. The work was carried out in a rigid bed flume, with the model spur dikes simulated using 2-D, rectangular plates placed along one side of the flume. For various discharges, the drag force of each individual spur dike in the spur dike field was measured directly using a specially-designed apparatus. For these tests, both submerged and unsubmerged conditions were evaluated along with various numbers of spur dikes and various relative spacings between the spur dikes throughout the field. It was concluded that the configuration of a spur dike field in terms of the number of spur dikes and relative spacing between the spur dikes has a substantial impact on the drag force and hence the flow resistance and backwater effect of a spur dike field. The most upstream spur dike had the highest drag force amongst the spur dikes in the field, and it acted as a shield to decrease the drag force exerted by the downstream spur dikes. From the experiments on the submerged spur dikes, it was observed that the jet flow over the spur dikes has an important effect on the flow structure and hence the flow resistance.<p> In the third stage of the study, the flow field within the vicinity of a single submerged spur dike was modeled using the three-dimensional (3-D) computational fluid dynamic (CFD) software FLUENT. Application of the software required solution of the 3-D Reynolds-averaged Navier-Stokes equations wherein the Reynolds stresses were resolved using the RNG ê-å turbulence model. One discharge condition was evaluated in a smooth, rectangular channel for two conditions, including uniform flow conditions without the spur dike in place and one with the spur dike in place. The CFD model was evaluated based on some experimental data acquired from the physical model. It was found that the CFD model could satisfactorily predict the flow resistance and water surface profile adjacent to the spur dike, including the resulting backwater effect. Furthermore, the CFD model gave a good prediction of the velocity field except for the area behind the spur dike where the effects of diving jet flow over the spur dike was not properly modeled.