Thesis (MScEng)--Stellenbosch University, 2011. / ENGLISH ABSTRACT: The hydrodynamic characteristics of three typical ship-like hull sections, in different motions, are numerically investigated using FLUENT, 2009. These simple shapes, namely a v-bottom (triangle) hull, a at-bottom (square) hull and a round-bottom (semi-circle) hull, are investigated in uncoupled heave, sway and roll. The problem is described in two dimensions. A combination of numerical methods and models, found in literature, are used to conduct this investigation. Hull characterisation is achieved through the use of hull mass and damping coe cients. These numerically determined coe cients are compared to experimental work conducted by Vugts (1968). A good correlation between the numerical and experimental results exists for the heave and sway cases. By normalising the coe cients, different hulls are comparable to one another. The numerical models used are validated and veri ed. Roll motion remains largely unsolved for very large angles of roll (in excess of 11°). Different uid ow phenomena occurring around the hull sections have varying degrees of in uence on the motions of a hull. It is found that not one of the turbulence models investigated can be employed to globally solve each type of hull-motion case. Also, forced oscillations in computational simulations require considerably more computational time than free-decay oscillating hull simulations. / AFRIKAANSE OPSOMMING: Die hidrodinamiese karakteristieke van verskillende skeepsrompvorms, in verskeie bewegingswieë, is numeries ondersoek met behulp van FLUENT, 2009. Drie eenvoudige vorms ('n v-bodem (driehoek), plat-bodem (reghoek) en rondebodem (semi-sirkel) romp) is onderskeidelik ondersoek in opwieg, dwarswieg en rol. Die probleem is twee-dimensioneel. Daar is gebruik gemaak van 'n kombinasie van numeriese metodes en modelle, uit die literatuur, om die ondersoek uit te voer. Die rompe is gekarakteriseer met behulp van massa- en dempingskoëffi siënte. Hierdie numeries bepaalde koë ffisiënte is vergelyk met die eksperimentele werk van Vugts (1968). Daar bestaan 'n goeie korrelasie tussen die numeriese en eksperimentele resultate vir die opwieg en dwarswieg gevalle. Die koë ffisiënte is genormaliseer om die verskeie rompvorms te vergelyk. Die numeriese modelle is geverifi eer en valideer. Rolbewegings is onopgelos vir groot rolhoeke (groter as 11°). Die mate waartoe die romp se beweging beïnvloed word deur die verskillende vloei verskynsels wat om die rompe ontstaan, verskil. Daar is bevind dat geen van die turbulensie modelle gebruik kan word om alle skeepsbeweging-gevalle op te los nie. Gedwonge-ossilasie numeriese simulasies benodig meer berekeningstyd as vrye-verval ossilasie gevalle.

La compréhension de l'aérothermique d'un véhicule durant sa phase de développement est une question essentielle afin d'assurer, d'une part, un bon refroidissement et une bonne efficacité de ses composants et d'autre part de réduire la force de traînée et évidement le rejet des gaz à effet de serre ou la consommation d'essence. Cette thèse porte sur la simulation numérique et la validation expérimentale de l'aérothermique d'un véhicule à trois roues dont deux, en avant et une roue motrice en arrière. La simulation numérique est basée sur la résolution des équations de conservation de la masse, de la quantité de mouvement et de l'énergie en utilisant l'approche RANS (Reynolds-Averaged Navier-Stokes). Le rayonnement thermique est modélisé grâce à la méthode S2S (Surface to Surface) qui suppose que le milieu séparant les deux surfaces rayonnantes, ici de l'air, ne participe pas au processus du rayonnement. Les radiateurs sont considérés comme des milieux poreux orthotropes où la perte de pression est calculée en fonction de leurs propriétés inertielle et visqueuse; leur dissipation thermique est modélisée par la méthode Dual flow. Une première validation de l'aérodynamique est faite grâce à des essais en soufflerie. Ensuite, une deuxième validation de la thermique est faite grâce à des essais routiers. Un deuxième objectif de la thèse est consacré à la simulation numérique de l'aérodynamique en régime transitoire du véhicule. La simulation est faite à l'aide de l'approche Detached eddy simulation (DES). Une validation expérimentale est faite à partir d'étude en soufflerie grâce à des mesures locales de vitesse à l'aide de sondes cobra.

Doctor of Philosophy / Department of Mechanical and Nuclear Engineering / M.H. Hosni / Z.C. Zheng / In order to study the capability of computational methods in investigating the mechanisms associated with disease and contaminants transmission in aircraft cabins, the Computational Fluid Dynamics (CFD) models are used for the simulation of turbulent airflow, tracer gas diffusion, and particle dispersion in a generic aircraft cabin mockup. The CFD models are validated through comparisons of the CFD predictions with the corresponding experimental measurements. It is found that using Large Eddy Simulation (LES) with the Werner-Wengle wall function, one can predict unsteady airflow velocity field with relatively high accuracy. However in the middle region of the cabin mockup, where the recirculation of airflow takes place, the accuracy is not as good as that in other locations. By examining different k-ε models, the current study recommends the use of the RNG k-ε model with the non-equilibrium wall function as a Reynolds Averaged Navier Stokes (RANS) model for predicting the steady-state airflow velocity data. It is also found that changing the cabin air-inlet nozzle height has a significant effect on the flow behavior in the middle and upper part of the cabin, while the flow pattern in the lower part is not affected as much. Through the use of LES and species transport model in simulating tracer gas diffusion, very good agreement between predicted and measured tracer gas concentration data is observed for some monitoring locations, but the agreement level is not uniform for all the sampling point locations. The reasons for the deviations between predictions and measurements for those locations are discussed. The Lagrange-Euler approach is invoked in the particle dispersion simulations. In this approach, the equation of motion for the discrete phase is coupled with the continuous phase governing equations through the calculation of drag and buoyancy forces acting on particles. The continuous phase flow is turbulent and RANS is employed in order to calculate the continuous phase velocity field. A complete study on grid dependence for RANS simulation is performed through a controllable regional mesh refinement scheme. The grid dependence study shows that using unstructured grid with tetrahedral and hybrid elements in the refinement region are more efficient than using structured grid with hexahedral elements. The effect of turbulence on the particle dispersion is taken into account by using a stochastic tracking method (Discrete Random Walk model). One of the significant features of this study is the investigation of the effect of the number of tries on the accuracy of particle concentration predictions when Discrete Random Walk is used to model turbulent distribution of particles. Subsequently, the optimum number of tries to obtain the most accurate predictions is determined. In accordance with the corresponding experimental data, the effect of particle size on particle distribution is also studied and discussed through the simulation of two different sizes of mono-disperse particles in the cabin with straight injection tube, i.e., 3µm and 10µm. Due to the low particle loading, neglecting the effect of particles motion on the continuous phase flow-field seems to be a reasonable, simplifying assumption in running the simulations. However, this assumption is verified through the comparison of the results from 1-way and 2-way coupling simulations. Eventually through the simulations for the particle injection using the cone diffuser, the effects of cabin pressure gradient as well as the particle density on particles dispersion behavior are studied and discussed. In the last part of this dissertation, the turbulent airflow in a full-scale Boeing 767 aircraft cabin mockup with eleven rows of seats and manikins is simulated using steady RANS method. The results of this simulation cannot only be used to study the airflow pattern, but also can be used as the initial condition for running the tracer gas diffusion and particle dispersion simulations in this cabin mockup.

Airfoils &deg / apping in pitch and plunge are studied, and the &deg / apping motion parameters are op- timized to maximize thrust generation and the e&plusmn / ciency of the thrust generation. Unsteady viscous &deg / ow&macr / elds over &deg / apping airfoils are computed on overset grids using a Navier-Stokes solver. Computations are performed in parallel using Parallel Virtual Machine library routines in a computer cluster. A single &deg / apping airfoil and dual airfoils &deg / apping in a biplane con- &macr / guration are considered. A gradient based optimization algorithm is employed. The thrust production and the e&plusmn / ciency of the thrust production are optimized with respect to &deg / apping parameters / the plunging and pitching amplitudes, the &deg / apping frequency, and the phase shift between the pitch and plunge motions. It is observed that thrust generation of &deg / apping airfoils strongly depends on the phase shift and high thrust values may be obtained at the expense of reduced e&plusmn / ciency. For a high e&plusmn / ciency in thrust generation, the e&reg / ective angle of attack of the airfoil is reduced and large scale vortex formations at the leading edge are prevented. At a &macr / xed reduced &deg / apping frequency of 1, a single &deg / apping airfoil in pitch and plunge motion produces the maximum average thrust coe&plusmn / cient of 1:41 at the plunge amplitude of 1:60, the pitch amplitude of 23:5o, and the phase shift of 103:4o whereas the maximum e&plusmn / ciency of 67:5% is obtained at the plunge amplitude of 0:83, the pitch amplitude of 35:5o and the phase shift of 86:5o.

A computational approach to analyze the hydrodynamic performance of a hydrofoil or a marine propeller undergoing unsteady motion has been developed. In order to simulate heave and pitch motion of a hydrofoil, an unsteady boundary element method based modeling is performed. The wake of the hydrofoil is modeled by a continuous dipole sheet and determined in time by applying a force-free condition on its surface. An explicit vortex core model is adapted in this model to capture the rolling up shape and to avoid instability due to roll-up deformation of the wake. The numerical results of the developed model are compared with analytical results and those from the commercial Reynolds-Averaged Navier-Stokes solver (ANSYS/FLUENT). The results show close level of agreement with each other. The problem of flow around a marine propeller performing surge, roll and heave motion in an unbounded fluid is formulated and solved using both a vortex-lattice method and a boundary element method. A fully unsteady wake alignment algorithm is implemented into the vortex-lattice method in order to satisfy the force-free condition on the propeller wake surface. Finally, a comparative study of transient propeller forces on a propeller blade obtained from BEM and VLM (with or without fully aligned wake) is carried out and results are presented. In some cases, results from the presented methods are compared with those from RANS or other numerical methods available in the literature. / text

L'écoulement de jeu est un phénomène très important dans les turbomachines. Il provient du mouvement relatif entre la pale et la paroi d'extrémité, et la différence de pression à travers la pale. L'écoulement de jeu est extrêmement complexe pour sa nature tridimensionnelle et instable, et son existence conduit à de nombreux effets défavorables, par exemple, les pertes de performance aérodynamique et les émissions de bruit. C'est pourquoi l'écoulement de jeu a motivé de nombreuses recherches expérimentales et numériques. Afin d'améliorer la compréhension du écoulement de jeu et le bruit de large bande associé, une campagne de recherche a été menée au LMFA. En ce qui concerne l'écoulement de jeu, cette campagne de recherche comprend une expérience avec des technologies de mesure avancées, un calcul zonal LES et une série de calculs RANS / URANS. L'expérience et les simulations considèrent une configuration simple de l'écoulement de jeu à un faible nombre de Mach. Les résultats expérimentaux et numériques sont analysés de façon systématique et approfondie dans la présente étude. Enfin, des efforts sont déployés pour la modélisation et la prédiction du bruit à large bande avec des résultats expérimentaux et numériques. On observe dans l'expérience un système à multiple-tourbillon, avec une tourbillon de jeu intense. Les différentes analyses sur les caractéristiques d'écoulement montrent un bon accord entre l'expérience et le ZLES dans la région du écoulement de jeu. L'approche zonale (RANS-LES) s'avère être un outil puissant pour fournir une description détaillée du écoulement de jeu, avec un coût de calcul limité. Cependant, les calculs RANS et URANS surestiment globalement la diffusion de la tourbillon. En outre, l'oscillation du tourbillon de jeu est étudiée en utilisant des champs instantanés de PIV et l'amplitude d'oscillation est évaluée. La réponse dynamique de la tourbillon de jeu est également étudiée avec URANS aux fréquences choisies. Deux modèles de prédiction du bruit en champ lointain, correspondant à deux sources acoustiques différentes, sont reformulés et mis en oeuvre avec les données de champ proche des simulations numériques. Ces prédictions sont comparées aux mesures à champ lointain. En utilisant les données ZLES, le modèle de l’écoulement de jeu sur-estime le bruit généré dans la région de jeu. Le modèle de bruit de bord de fuite est implémenté avec les données ZLES et les données RANS et fournit une très bonne prédiction dans une large bande de fréquence. / The tip-leakage flow is a common flow feature in turbomachines. It originates from the relative motion between the blade tip and the end-wall, and the pressure difference across the blade. The tip-leakage flow is extremely complex for its three-dimensional unsteady nature, and its existence leads to many unfavourable effects, such as aerodynamic performance losses and noise emissions. These issues have motivated extensive experimental and numerical researches from both aerodynamic and aeroacoustic points of view. In order to improve the understanding of the tip-leakage flow and its associated broadband noise, a research campaign has been carried out at LMFA. Regarding the tip-leakage flow, this research campaign includes an experiment with advanced measurement technologies, a zonal LES computation and a series of RANS/URANS computations. Both the experiment and the simulations consider a single-airfoil configuration at low Mach number. Experimental and numerical results are analysed systematically and thoroughly in the current study. Finally, efforts are put on the broadband noise modelling and prediction based on the experimental and numerical results. A multi-vortex system with an intense tip-leakage vortex is observed in the experiment. The various analyses of the flow characteristics show a good agreement between the experiment and the ZLES in the blade tip region. The zonal (RANS-LES) approach proves itself to be a powerful tool to provide a detailed description of the tip-leakage flow, with a limited computational cost. However, the RANS and URANS computations globally over-estimate the diffusion of the tip-leakage vortex. Furthermore, the random oscillation of the tip-leakage vortex is investigated using PIV instantaneous flow fields and the wandering amplitude is evaluated. The dynamic response of the tip-leakage vortex is also studied with URANS at selected frequencies. Two far-field noise prediction models, corresponding to two different acoustic sources, are reformulated and implemented with the near-field data from the numerical simulations. These predictions are compared to the far-field measurements. Using the ZLES data as input, the blade-tip self-noise model is found to over-estimate the noise generated in the blade-tip region. The trailing-edge noise model is implemented with the time-averaged ZLES and the RANS near-field data, and yields a very good prediction within a broad range of frequency.

This Thesis presents simulation of flow past laboratory-scale and full-scale tidal stream turbines (TST) using EDF's open-source CFD solver Code_Saturne. The work shows that detailed results may be obtained with confidence and that greater information on the loading and wake structure is available than other methods, such as blade element momentum theory.Results are obtained using a new sliding-mesh method that has been implemented in Code_Saturne as part of this work. The sliding-mesh method uses internal Dirichlet boundary conditions with values on the interface prescribed via a halo-point method. Parallel performance is optimised by a carefully-chosen method of exchanging information between specific processes. Validation is provided for flow past a rotating cylinder and a sphere.For the laboratory-scale TST, Reynolds-Averaged Navier-Stokes models are used to model turbulence. The k-omega-SST and Launder-Reece-Rodi (LRR) models yield good agreement with experimental values of power and thrust coefficients as a function of tip-speed ratio (TSR). The standard k-epsilon model is shown to perform poorly due to an overprediction of turbulent kinetic energy upstream of the rotor plane. The k-omega-SST model is then used to examine wake behaviour for parametric studies of turbulence intensity and TSR. Increased turbulence levels are shown to reduce the downstream propagation of the wake because of increased mixing. The near wake is influenced by the TSR, whilst the far wake is independent of TSR.The predicted effect of tidal conditions typical of the EMEC test site are considered for flow past Tidal Generation Limited's 1MW TST. The effect of sheared-velocity profiles leads to an increase in loading on an individual turbine blade at the point of a rotation where velocity shear is greatest. The effect of increased yaw angle leads to large fluctuations of the power coefficient, but smaller fluctuations of the thrust coefficient. Mean values of thrust and power decrease as a function of the cosine of the yaw angle and yaw angle squared respectively.

<p>This work utilized Pointwise and Fluent to generate a two-dimensional axisymmetric model a Mk 48 torpedo, with the intention of informing methods to reduce the turbulence, and therefore hydrodynamic noise, of the torpedo’s wake. However, this work was unable to gather data on the unsteady nature of the turbulence expected around the torpedo due to Fluent providing unrealistic results when run using a transient solver. This work shows that the transient solver computed boundary layers greater than one order of magnitude smaller than expected, and in some cases there was no change in boundary layer thickness over the torpedo’s body. The work does contain steady state solutions that were validated by first performing a grid convergence study for a flat plate. The steady state results for the flat plate and torpedo both showed the expected growth for a turbulent boundary layer. Additionally, there was a high level of convergence with the Log-Law showing that the steady state data is valid. Future work should use a transient solver to determine the characteristics of the turbulence to resolve unsteady flow from vortex shedding, wake characteristics, and any broadband or narrowband noise to develop solutions to reduce the noise made by the Mk 48.</p>

This thesis focused on investigating flow, coherent structures, and mechanisms leading to scour around a series of three submerged groynes in a high curvature (radius of curvature (R)/channel width (B)=1.5) channel bend using a Large Eddy Simulation Numerical (LES) model. Flow was investigated during both an initial and a later stage of scour. The results showed that the groynes appeared effective in keeping the main core of high streamwise velocity away from the outer bank wall in the region where they were installed, although high potential still existed for local scour around the groynes. During the initial stage of scour, horseshoe vortices (HVs) showed the greatest propensity to induce scour immediately upstream of the groyne tips. During the later stage of scour, the HV in front of the first upstream groyne (G1) induced very high mean pressure fluctuations on the outer bank wall. Scour was also of very great concern around the tip of G1 due to severe mean bed pressure fluctuations. Downstream of the groyne field, the presence of a counter-rotating outer bank cell was capable of endangering the stability of the outer bank. The second focus of this thesis was to investigate flow in a 135° channel bend using both Reynolds Averaged Navier Stokes (RANS) and LES numerical models. The RANS study examined the effects of curvature ratio (R/B), and aspect ratio (B/H, where H is the inlet flow depth), on secondary circulation strength, and bed shear stresses. The study revealed that a decrease in R/B was associated with an increase in secondary circulation strength and peak bed shear stress. A change in B/H also substantially affected cross stream circulation strength. The LES study was conducted in a 135° (R/B = 1.5) bend flume with a fixed bed corresponding to near equilibrium scour conditions, and the results were compared to a similar high curvature 193° bend numerical study. Inner bank vortices and shear layers were present in both cases although their characteristics were substantially different. Distributions of boundary friction velocities, and turbulence were also quite different for each case.

Hot wire anemometers have been used in several wind velocity sensors deployed in Mars. They are based in keeping the temperature of a surface at a constant value, above the ambient. This is done by means of a heater controlled with an electronic system. The cooling rate of each point at the sensor surface can be used to calculate the wind velocity and direction. However, due to turbulent fluctuations, the cooling rate is not constant even in the case of constant velocity. Moreover, RANS simulations cannot estimate such fluctuations as they only provide an estimation of the averaged flow field. The goal of this work has been to estimate such fluctuations and the e↵ect they might have on the sensor readings. To do so, the turbulent cooling rate (Nusselt number) of a sensor with a generic shape, under the typical conditions to be find in Mars, has been simulated using high performance LES (Large Eddy Simulation) simulations and compared with RANS and URANS simulations.