Thesis (MScEng (Mechanical and Mechatronic Engineering))--University of Stellenbosch, 2010. / ENGLISH ABSTRACT: This study is aimed at the development of an efficient and reliable method of evaluating the performance of an air-cooled steam condenser (ACSC) under windy conditions, using computational fluid dynamics (CFD). A two-step modelling approach is employed as a result of computational limitations. The numerical ACSC model developed in this study makes use of the pressure jump fan model, amongst other approximations, in an attempt to minimize the computational expense of the performance evaluation. The accuracy of the numerical model is verified through a comparison of the numerical results to test data collected during full scale tests carried out on an operational ACSC. Good correlation is achieved between the numerical results and test data. Further verification is carried out through a comparison to previous numerical work. Satisfactory convergence is achieved for the most part and the few discrepancies in the results are explained. The effect of wind on ACSC performance at El Dorado Power Plant (Nevada, USA) is investigated and it is found that reduced fan performance due to distorted flow at the inlet of the upstream fans is the primary contributor to the reduction in performance associated with increased wind speed in this case. An attempt is subsequently made to identify effective wind effect mitigation measures. To this end the effects of wind screens, solid walkways and increasing the fan power are investigated. It is found that the installation of an appropriate wind screen configuration provides a useful means of reducing the negative effects of wind on ACSC performance and an improved wind screen configuration is suggested for El Dorado. Solid walkways are also shown to be beneficial to ACSC performance under windy conditions. It is further found that ACSC performance increases with walkway width but that the installation of excessively wide walkways is not justifiable. Finally, increasing the fan power during periods of unfavourable ambient conditions is shown to have limited benefit in this case. The model developed in this study has the potential to allow for the evaluation of large ACSC installations and provides a reliable platform from which further investigations into improving ACSC performance under windy conditions can be carried out. / AFRIKAANSE OPSOMMING: Hierdie studie is daarop gemik om die ontwikkeling van 'n geskikte en betroubare metode van evaluering van die verrigting van ’n lugverkoelde stoom-kondensator (air-cooled steam condenser, ACSC) onder winderige toestande, met behulp van numeriese vloei-dinamika. ’n Twee-stap modelleringsbenadering is aangewend as gevolg van rekenaar beperkings. Die numeriese ACSC-model wat in hierdie studie ontwikkel is, maak gebruik van die druksprong waaier model, asook ander benaderings, in ’n poging om die berekeningskoste van die verrigting-evaluering te verminder. Die akkuraatheid van die numeriese model is bevestig deur middel van ’n vergelyking van die numeriese resultate met toetsdata ingesamel tydens die volskaal toetse uitgevoer op ’n operasionele ACSC. Goeie korrelasie is bereik tussen die numeriese resultate en toetsdata. Verdere bevestiging is uitgevoer deur middel van ’n vergelyking met vorige numeriese werk. Bevredigende konvergensie is in die algemeen bereik en die paar verskille in die resultate word verduidelik. Die effek van wind op ACSC verrigting by El Dorado Power Plant (Nevada, VSA) is ondersoek, en daar is bevind dat verlaagde waaierverrigting, as gevolg van vervormde vloei by die inlaat van die stroomop waaiers, die primêre bydraer is tot die afname in ACSC werkverrigting geassosieer met verhoogde windsnelheid in hierdie geval. ’n Poging word dan aangewend om effektiewe wind-effek velagingsmaatreëls te identifiseer. Windskerms, soliede wandelvlakke en die verhoging van die waaierkrag word gevolglik ondersoek. Daar is bevind dat die installasie van ’n toepaslike windskerm-opset ’n nuttige middel tot ’n vermindering van die negatiewe effekte van wind op ACSC verrigting bied, en ’n verbeterde windskerm opset is voorgestel vir El Dorado. Soliede wandelvlakke word ook aanbeveel as voordelig vir ACSC verrigting onder winderige toestande. Dit is verder bevind dat die ACSC prestasie verhoog met wandelvlak breedte, maar dat die installasie van ’n te ruim wandelvlak nie regverdigbaar is nie. Ten slotte, word bewys dat die verhoging van die waaierkrag tydens periodes van ongunstige omgewingsomstandighede ’n beperkte voordeel in hierdie geval het. Die model wat ontwikkel is in hierdie studie het die potensiaal om voorsiening te maak vir die evaluering van groot ACSC- installasies en bied ’n betroubare platform vanwaar verdere ondersoeke tot die verbetering van ACSC verrigting onder winderige toestande uitgevoer kan word.

Thesis (MScEng (Mechanical and Mechatronic Engineering))--University of Stellenbosch, 2010. / ENGLISH ABSTRACT: Cooling tower inlet losses and effective flow diameter under no crosswind conditions and the pressure distribution around a circular cylinder subjected to a crosswind are modelled using CFD. The CFD model used to evaluate the inlet losses is validated with data measured in an experimental cooling tower sector model and data obtained from literature. The effect of different inlet geometries on the inlet loss coefficient and the effective diameter are investigated in order to improve cooling tower inlet designs. CFD models are developed to investigate the pressure distribution around infinite and finite circular cylinders. The infinite cylinder is modelled with a smooth surface and a rough surface so that the results can be compared to experimental data from literature. Ultimately a finite cylinder model with a rough surface is developed and the results are compared to experimental data from literature. / AFRIKAANSE OPSOMMING: Koeltoring inlaatverlies en effektiewe vloei deursnit onder geen teenwind toestande en die drukverdeling rondom ‘n sirkelvormige silinder, onderworpe aan ‘n teenwind, word gemodelleer deur gebruik te maak van “CFD”. Die “CFD” model wat gebruik word om die inlaatverlies te evalueer is gevalideer met data verkry vanaf ‘n eksperimentele koeltoring sektor model. Verder word die “CFD” model gebruik in ‘n ondersoek om te bebaal wat die effek is van verskillende inlaat geometrieë op die inlaat verlies koeffisiënt en die effektiewe diameter sodat die inlaat geometrie van koeltorings verbeter kan word. ‘n “CFD” model word dan ontwikkel om die druk verdeling rondom ‘n sirkelvormige silinder te ondersoek. Die silinder word as oneindig gesimuleer met ‘n glade en ruwe wand sodat die resultate vergelyk kan word met eksperimentele data verkry vanaf literatuur. Die afdeling word afgesluit deur die silinder as eindig met ‘n ruwe wand te simuleer en dan word die resultate vergelyk met eksperimentele data verkry vanaf literatuur.

Thesis (MScEng)--Stellenbosch University, 2012. / ENGLISH ABSTRACT: The use of micro gas turbines (MGTs) for the propulsion of unmanned aerial vehicles (UAVs) has become an industry standard. MGTs offer better performance vs. weight than similar sized, internal combustion engines. The front component of an MGT serves the purpose of compressing air, which is subsequently mixed with a fuel and ignited to both power the turbine which drives the compressor, and to produce thrust. Centrifugal compressors are typically used because of the high pressure ratios they deliver per stage. The purpose of this project was to design a centrifugal compressor impeller, and to devise a methodology and the tools with which to perform the aforementioned. A compressor impeller adhering to specific performance and dimensional requirements was designed. The new compressor was designed using a mean-line performance calculation code. The use of the code was vindicated through comparison with the results from a benchmark study. This comparison included mean-line, Computational Fluid Dynamic (CFD), and experimental results: the new design mean-line results were compared to the results of CFD simulations performed on the same design. The new design was optimised using an Artificial Neural Network (ANN) and Genetic Algorithm. Prior to and during optimisation, the ANN was trained using a database of sample CFD calculations. A Finite Element Analysis (FEA) was done on the optimised impeller geometry to ensure that failure would not occur during operation. According to CFD results, the final design delivered good performance at the design speed with regards to pressure ratio, efficiency, and stall margin. The mechanical stresses experienced during operation were also within limits. Experimental results showed good agreement with CFD results of the optimised impeller. Keywords: micro gas turbine, centrifugal compressor, impeller, CFD, experimental, optimisation, FEA. / AFRIKAANSE OPSOMMING: Die gebruik van mikrogasturbines vir die aandrywing van onbemande vliegtuie het ‟n standaard geword in die industrie. Mikrogasturbines bied beter werkverrigting teen gewig as binnebrandenjins van soortgelyke grote. Hierdie eienskap verseker dat mikrogasturbines as aandryfmotors vir onbemande vliegtuie uiters voordelig is. Die voorste komponent van ‟n mikrogasturbine dien om lug saam te pers, wat dan met brandstof gemeng en daarna aan die brand gesteek word om krag aan die kompressor en stukrag te voorsien. Sentrifugaalkompressors word tipies gebruik as gevolg van die hoë drukverhoudings wat hierdie komponente per stadium kan lewer. Die doel van hierdie projek was om ‟n sentrifugaalkompressor te ontwerp, en ‟n metode en die hulpmiddels te ontwikkel om laasgenoemde uit te voer. ‟n Kompressor rotor wat voldoen het aan sekere werkverrigtings en dimensionele vereistes is ontwerp. Die nuwe kompressor rotor is met behulp van 1-dimensionele werkverrigting-berekeningskode ontwerp. Die berekeningsakkuraatheid van die kode en díé van ‟n kommersiële Berekenings Vloeidinamika pakket is bevestig deur die berekende resultate te vergelyk met die van eksperimente. Die nuwe rotor is gevolglik deur middel van ‟n Kunsmatige Neurale Netwerk en Genetiese Algoritme geoptimeer. Die Kunsmatige Neurale Netwerk is voor en gedurende optimering deur Berekenings Vloeidinamika simulasies opgelei. Die meganiese sterkte van die geoptimeerde rotor is nagegaan met behulp van ‟n Eindige Element Analise. Dit is gedoen om te verseker dat die rotor nie sal faal by die bedryfspunt nie. Berekenings Vloeidinamika resultate het getoon dat die finale rotor ontwerp ‟n goeie werkverrigting lewer by die ontwerpspoed, met betrekking tot drukverhouding, bennutingsgraad, en stakingsmarge. Eksperimentele resultate het goeie ooreenstemming met die Berekenings Vloeidinamika resultate van die geoptimeerde rotor getoon. Sleutelwoorde: mikrogasturbine, sentrifigaalkompressor, rotor, Berekenings Vloeidinamika, eksperimenteel, optimering, Eindige Element Analise.

Abrasive flow machining (AFM) is a non-traditional manufacturing technology used to expose a substrate to pressurised multiphase slurry, comprised of superabrasive grit suspended in a viscous, typically polymeric carrier. Extended exposure to the slurry causes material removal, where the quantity of removal is subject to complex interactions within over 40 variables. Flow is contained within boundary walls, complex in form, causing physical phenomena to alter the behaviour of the media. In setting factors and levels prior to this research, engineers had two options; embark upon a wasteful, inefficient and poor-capability trial and error process or they could attempt to relate the findings they achieve in simple geometry to complex geometry through a series of transformations, providing information that could be applied over and over. By condensing process variables into appropriate study groups, it becomes possible to quantify output while manipulating only a handful of variables. Those that remain un-manipulated are integral to the factors identified. Through factorial and response surface methodology experiment designs, data is obtained and interrogated, before feeding into a simulated replica of a simple system. Correlation with physical phenomena is sought, to identify flow conditions that drive material removal location and magnitude. This correlation is then applied to complex geometry with relative success. It is found that prediction of viscosity through computational fluid dynamics can be used to estimate as much as 94% of the edge-rounding effect on final complex geometry. Surface finish prediction is lower (~75%), but provides significant relationship to warrant further investigation. Original contributions made in this doctoral thesis include; 1) A method of utilising computational fluid dynamics (CFD) to derive a suitable process model for the productive and reproducible control of the AFM process, including identification of core physical phenomena responsible for driving erosion, 2) Comprehensive understanding of effects of B4C-loaded polydimethylsiloxane variants used to process Ti6Al4V in the AFM process, including prediction equations containing numerically-verified second order interactions (factors for grit size, grain fraction and modifier concentration), 3) Equivalent understanding of machine factors providing energy input, studying velocity, temperature and quantity. Verified predictions are made from data collected in Ti6Al4V substrate material using response surface methodology, 4) Holistic method to translating process data in control-geometry to an arbitrary geometry for industrial gain, extending to a framework for collecting new data and integrating into current knowledge, and 5) Application of methodology using research-derived CFD, applied to complex geometry proven by measured process output. As a result of this project, four publications have been made to-date – two peer-reviewed journal papers and two peer-reviewed international conference papers. Further publications will be made from June 2014 onwards.

Continued growth of the aviation industry and increasingly strict noise requirements set by international bodies and airport authorities alike means that novel methods of reducing aircraft noise must be found. Engine noise represents a majority contribution to total aircraft noise during take-off and turbulent mixing of the exhaust gases is the dominant noise source of the engine at take-off. While bypass ratio has been the historical, and rather convenient means, of reducing jet noise, an upper limit to bypass ratio is now being approached and additional means of reducing jet noise must be found. One method that has shown potential for reducing aeroacoustic jet noise is the application of small, high pressure jets to the circumference of the jet nozzle. These jets, termed microjets, have the advantage over static devices that the microjets can be activated only when the noise benefit is required and deactivated when emitted noise is not an issue, such as in cruise, thereby reducing the thrust penalty associated with the devices over the majority of the flight. Large eddy simulations have been performed to investigate the impact that the addition of microjets has on the aerodynamic flowfield and radiated far-field noise of a high Reynolds number, Mach 0.9, propulsive, laboratory scale jet. Far-field noise was predicted through a new implementation of the permeable Ffowcs Williams Hawkings surface method in the solver. In addition to single-point flowfield statistics and far-field noise, spatio-temporal second- and fourth-order correlations are investigated. Two pairs of simulations were conducted, a coarse mesh containing 100 million elements and a fine mesh with 200 million elements. The coarse mesh included an azimuthal clustering of the cells in the near-microjet region. The non-uniformity of the azimuthal cell size was shown to adversely affect the development of the initial shear layer, yielding a delay in transition to a fully turbulent state and larger coherent structures in regions with larger cells. Radial velocity and turbulent kinetic energy profiles show good agreement with experimental results. A previously unidentified periodic interaction between the main jet and microjets was found. The dynamic interaction gives rise to velocity and pressure fluctuations in the near microjet region that match a tonal frequency found in the microjet far-field spectra that is absent from the clean jet case. Second- and fourth-order correlation distributions show large periodic regions of high correlation amplitude in the near microjet region. The evidence demonstrates that the main-microjet interaction is a clear high-frequency noise source. Despite the high-frequency noise associated with the main-microjet interaction, the addition of microjets yields a 1-2 dB reduction in overall sound pressure level. Additionally, over a significant portion of the length of the potential core the microjets reduce the amplitude of the majority of the six main correlation amplitudes that can be used in far-field noise prediction. Finally, the generation of the counter-rotating vortex pair downstream of the microjets was investigated. It is commonly presumed that this vortex pair is similar in origin to the counter-rotating vortex pair present in a jet in a crossflow. Vortex identification methods, velocity vectors and streamlines in the near microjet region demonstrate that the horseshoe-like vortex is the source of the counter rotating vortex pair that is present downstream of the microjets. The horseshoe-like vortex in the microjet case has the same sense as the vortices in the microjet shear layer and appears to be generated by the development of a recirculation region of microjet fluid during the main-microjet interaction.

Abdominal aortic aneurysm (AAA) disease involves a dilation of the aorta below the renal arteries. If the aneurysm becomes sufficiently dilated and tissue strength is less than vascular pressure, rupture of the aorta occurs entailing a high mortality rate. Despite improvements in surgical technique, the mortality rate for emergency repair remains high and so an accurate predictor of rupture risk is required. Inflammation and the associated recruitment of monocytes into the aortic wall are critical in the pathology of AAA disease, stimulating the degradation and remodeling of the vessel wall. Areas with high concentrations of macrophages may experience an increase in tissue degradation and therefore an increased risk of rupture. Determining the magnitude and distribution of monocyte recruitment can help us understand the pathology of AAA disease and add spatial accuracy to the existing rupture risk prediction models. In this study finite element computational fluid dynamics simulations of AAA haemodynamics are seeded with monocytes to elucidate patterns of cell deposition and probability of recruitment. Haemodynamics are first simulated in simplified AAA geometries of varying diameters with a patient averaged flow waveform inlet boundary condition. This allows a comparison with previous experimental investigations as well as determining trends in monocyte adhesion with aneurysm progression. Previous experimental investigations show a transition to turbulent flow occurring during the deceleration phase of the cardiac cycle. There has thus far been no investigation into the accuracy of turbulence models in simulating AAA haemodynamics and so simulations are compared using RNG κ − ε, κ − ω and LES turbulence models. The RNG κ − ε model is insufficient to model secondary flows in AAA and LES models are sensitive to inlet turbulence intensity. The probability of monocyte adhesion and recruitment depends on cell residence time and local wall shear stress. A near wall particle residence time (NWPRT)model is created incorporating a wall shear stress-limiter based on in vitro experimental data. Simulated haemodynamics show qualitative agreement with experimental results. Peaks of maximum NWPRT move downstream in successively larger geometries, correlating with vortex behaviour. Average NWPRT rises sharply in models above a critical maximum diameter. These techniques are then applied to patient-specific AAAs. Geometries are created from CT slices and velocity boundary conditions taken from Phase Contrast-MRI (PC-MRI) data for 3 patients. There is no gold standard for inlet boundary conditions and so simulations using 3 velocity components, 1 velocity component and parabolic flow profiles at the inlet are compared with each other and with PC-MRI data at the AAA midsection. The general trends in flow and wall shear stress are similar between simulations with 3 and 1 components of inlet velocity despite differences in the nature and complexity of secondary flow. Applying parabolic velocity profiles, however, can cause significant deviations in haemodynamics. Axial velocities show average to good correlation with PC-MRI data though the lower magnitude radial velocities produce high levels of noise in the raw data making comparisons difficult. Patient specific NWPRT models show monocyte infiltration is most likely at or around the iliac bifurcation.

No description available.

Master of Science / Department of Mechanical and Nuclear Engineering / B. Terry Beck / Mohammad H. Hosni / Cavitation is the change of a liquid to a two-phase mixture of liquid and vapor, similar to boiling. However, boiling generates a vapor by increasing the liquid temperature while cavitation generates vapor through a decrease in pressure. Both processes are endothermic, removing heat from the surroundings. Both the phase change and heat absorption associated with cavitation provide many engineering applications, including contributing to a new type of refrigeration cycle under development. Cavitation can occur at or below the vapor pressure; conditions that delay cavitation and allow for a metastable liquid are not well understood. A converging-diverging nozzle was designed and fabricated to create a low pressure region at the nozzle throat. The converging section of the nozzle increased the water velocity and decreased the pressure, according to Bernoulli’s principle. A cavitation front was formed slightly past the nozzle throat. The cavitation location suggested that the water was metastable near the nozzle throat. Flow through the system was controlled by changing the nozzle inlet and outlet pressures. The flowrate of water was measured while the outlet pressure was lowered. The flowrate increased as the outlet pressure dropped until cavitation occurred. Once cavitation initiated, the flow became choked and remained constant and independent of the nozzle outlet pressure. High-speed imagery was used to visualize the flow throughout the nozzle and the formation and collapse of cavitation in the nozzle’s diverging section. High-speed video taken from 1,000 to 35,000 frames per second captured the formation of the cavitation front and revealed regions of recirculating flow near the nozzle wall in the diverging section. Particle Image Velocimetry (PIV) was used to measure the velocity vector field throughout the nozzle to characterize flow patterns within the nozzle. PIV showed that the velocity profile in the converging section and throat region were nearly uniform at each axial position in the nozzle. In the diverging section, PIV showed a transient, high-velocity central jet surrounded by large areas of recirculation and eddy formation. The single-phase experimental results, prior to cavitation onset, were supplemented by Computational Fluid Dynamics (CFD) simulations of the velocity distribution using Fluent software.

Automotive catalytic converters are used extensively in the automotive industry to reduce toxic pollutants from vehicle exhausts. The flow across automotive exhaust catalysts is distributed by a sudden expansion and has a significant effect on their conversion efficiency. The exhaust gas is pulsating and flow distribution is a function of engine operating condition, namely speed (frequency), load (flow rate) and pressure loss across the monolith. The aims of this study are to provide insight into the development of the pulsating flow field within the diffuser under isothermal conditions and to assess the steady-state computational fluid dynamics (CFD) predictions of flow maldistribution at high Reynolds numbers. Flow measurements were made across an automotive catalyst monolith situated downstream of a planar wide-angled diffuser in the presence of pulsating flow. Cycle-resolved Particle Image Velocimetry (PIV) measurements were made in the diffuser and hot wire anemometry (HWA) downstream of the monoliths. The ratio of pulse period to residence time within the diffuser (J factor) characterises the flow distribution. During acceleration the flow remained attached to the diffuser walls for some distance before separating near the diffuser inlet later in the cycle. Two cases with J ~ 3.5 resulted in very similar flow fields with the flow able to reattach downstream of the separation bubbles. With J = 6.8 separation occurred earlier with the flow field resembling, at the time of deceleration, the steady flow field. Increasing J from 3.5 to 6.8 resulted in greater flow maldistribution within the monoliths; steady flow producing the highest maldistribution in all cases for the same Re. The oblique entry pressure loss of monoliths were measured using a one-dimensional steady flow rig over a range of approach Reynolds number (200 < Rea < 4090) and angles of incidence (0o < α < 70o). Losses increased with α and Re at low mass flow rates but were independent of Re at high flow rates being 20% higher than the transverse dynamic pressure. The flow distribution across axisymmetric ceramic 400 cpsi and perforated 600 cpsi monoliths were modelled using CFD and the porous medium approach. This requires knowledge of the axial and transverse monolith resistances; the latter being only applicable to the radially open structure. The axial resistances were measured by presenting uniform flow to the front face of the monolith. The transverse resistances were deduced by best matching CFD predictions to measurements of the radial flow profiles obtained downstream of the monolith when presented with non-uniform flow at its front face. CFD predictions of the flow maldistibution were performed by adding the oblique entry pressure loss to the axial resistance to simulate the monolith losses. The critical angle approach was used to improve the predictions, i.e. the oblique entry loss was limited such that the losses were assumed constant above a fixed critical angle, αc. The result showed that the perforated 600 cpsi monolith requires the entrance effect to be restricted above αc = 81o, while the losses were assumed constant above αc = 85o for the ceramic 400 cpsi monolith. This might be due to the separation bubble at the monolith entrance being restricted by the smaller hydraulic diameter of the perforated monolith thus limiting the oblique entry loss at the lower incidence angle.

This dissertation focuses on aerosol transport modeling in occupational environments and mining sites in Arizona using computational fluid dynamics (CFD). The impacts of human exposure in both environments are explored with the emphasis on turbulence, wind speed, wind direction and particle sizes. Final emissions simulations involved the digitalization process of available elevation contour plots of one of the mining sites to account for realistic topographical features. The digital elevation map (DEM) of one of the sites was imported to COMSOL MULTIPHYSICS® for subsequent turbulence and particle simulations. Simulation results that include realistic topography show considerable deviations of wind direction. Inter-element correlation results using metal and metalloid size resolved concentration data using a Micro-Orifice Uniform Deposit Impactor (MOUDI) under given wind speeds and directions provided guidance on groups of metals that coexist throughout mining activities. Groups between Fe-Mg, Cr-Fe, Al-Sc, Sc-Fe, and Mg-Al are strongly correlated for unrestricted wind directions and speeds, suggesting that the source may be of soil origin (e.g. ore and tailings); also, groups of elements where Cu is present, in the coarse fraction range, may come from mechanical action mining activities and saltation phenomenon. Besides, MOUDI data under low wind speeds (<2 m/s) and at night showed a strong correlation for particles 1-micrometer in diameter between the groups: Sc-Be-Mg, Cr-Al, Cu-Mn, Cd-Pb-Be, Cd-Cr, Cu-Pb, Pb-Cd, As-Cd-Pb. The As-Cd-Pb group correlates strongly in almost all ranges of particle sizes. When restricted low wind speeds were imposed more groups of elements are evident and this may be justified with the fact that at lower speeds particles are more likely to settle. When linking these results with CFD simulations and Pb-isotope results it is concluded that the source of elements found in association with Pb in the fine fraction come from the ore that is subsequently processed in the smelter site, whereas the source of elements associated to Pb in the coarse fraction is of different origin. CFD simulation results will not only provide realistic and quantifiable information in terms of potential deleterious effects, but also that the application of CFD represents an important contribution to actual dispersion modeling studies; therefore, Computational Fluid Dynamics can be used as a source apportionment tool to identify areas that have an effect over specific sampling points and susceptible regions under certain meteorological conditions, and these conclusions can be supported with inter-element correlation matrices and lead isotope analysis, especially since there is limited access to the mining sites. Additional results concluded that grid adaption is a powerful tool that allows to refine specific regions that require lots of detail and therefore better resolve flow detail, provides higher number of locations with monotonic convergence than the manual grids, and requires the least computational effort. CFD simulations were approached using the k-epsilon model, with the aid of computer aided engineering software: ANSYS® and COMSOL MULTIPHYSICS®. The success of aerosol transport simulations depends on a good simulation of the turbulent flow. A lot of attention was placed on investigating and choosing the best models in terms of convergence, independence and computational effort. This dissertation also includes preliminary studies of transient discrete phase, eulerian and species transport modeling, importance of saltation of particles, information on CFD methods, and strategies for future directions that should be taken.