Computational fluid dynamics (CFD) is a tool that is widely used in most industries today. It is important to have rigorous techniques to estimate the error produced when using CFD. This thesis develops techniques to estimate discretization error for unsteady flows using the unsteady error transport equation (ETE) as well as defect correction. A framework to obtain exact truncation error and estimated truncation error is also presented. The technique and results for the steady-state cases are given and the algorithm used for the steady case is extended for the unsteady case. Numerical results are presented for the steady viscous Burgers' equation, unsteady viscous Burgers' equation, steady quasi-1D Euler equations, and unsteady 1D Euler equations when applied to a shock tube. Cases using either defect correction or ETE are shown to give higher orders of accuracy for the corrected discretization error estimates when compared to the discretization error of the primal solution. / Master of Science / Computational fluid dynamics (CFD) is a tool that is widely used in most industries today. It is used to understand complex flows that are difficult to replicate using experimental techniques or by theoretical methods. It is important to have rigorous techniques to estimate the error produced when using CFD even when the exact solution is not available for comparison. This paper develops techniques to estimate discretization error for unsteady flows. Discretization error has one of the largest error magnitudes in CFD solutions. The exact physics dictates the use of continuous equations but to apply CFD techniques, the continuous equations have to be converted to discrete equations. Truncation error is, the error obtained when converting the continuous equations to discrete equations. This truncation error is in turn, the local source term for discretization error. To reduce the discretization error in the discrete equations, the exact or estimated truncation error is either added as a source term to the discrete equations or is used along with the error transport equation to get a better estimate of the solutions. A framework to obtain exact truncation error and estimated truncation error is also presented. The framework is first applied to the steady equations and is verified with results from previous studies and is then extended to the unsteady flows.

An investigation was carried out on the current state of the art in bubble modelling for computational fluid dynamics, and comparisons made between the different methods for both polydisperse and monodisperse multiphase flows. A multigroup method for polydisperse bubbly flows with the bubbles binned in terms of mass was selected from the various alternatives, which included other multigroup models and moment methods. The latter of these involve the integration of moments of the bubble number density function and transport of these quantities. The equations from this multigroup solver were then changed to more accurately and efficiently model cases involving extremely small bubbles over significant amounts of time, as the original model which was subsequently adapted had, as its primary purpose, simulation of larger bubbles over shorter periods of time. This was done by decoupling the gas and liquid momentum equations and adding an empirical rise velocity term for the bubbles. This new model was then partially implemented into OpenFOAM. The functioning of this new solver was confirmed by comparisons between the results and basic analytical solutions to the problems, as well as by means of comparison with another similar multiphase CFD solver (pbeTransportFoam). Following this confirmation of its functionality, the bubble model was implemented into another solver specifically designed for modelling wakes. Finally, the newly created solver was used to run some cases of interest involving a submerged wake. / Master of Science / Bubbles in the ocean are significant for a number of reasons, ranging from mixing of the upper layer of the ocean to scavenging of biological matter, by which means they can also impact the state of the ocean's surface where they are present. They serve as an important mechanism by which air is dissolved in the ocean, and their breaking at the surface can cause particles or droplets to be ejected into the atmosphere. They can be created by a variety of sources, ranging from the movement of ship propellers and hulls to natural processes, both abiotic and from microorganisms or other living things. They can have exceedingly variable sizes, meaning bubbles behave very differently from one another in the same area. For these reasons, their study is both interesting and sometimes challenging. In this research, methods were developed to simulate the movement over a significant amount of time of a wide size variety of very small bubbles within the ocean. First, study was undertaken of preexisting methods of bubble simulation and the different cases they were intended to represent. One of these existing methods was selected for use and then changed to more accurately represent smaller bubbles, as well as including simplifications to allow the simulations to run faster. Lastly, these methods were implemented into OpenFOAM, an open-source set of solvers for computational fluid dynamics (CFD). These new methods for simulation were finally applied to some cases involving submerged bubbles in the ocean and the movement of bubbles in these cases studied.

<div>Combustion instability due to thermo-acoustic interactions in high-speed propulsion devices such as gas turbines and rocket engines result from pressure waves with very large amplitudes propagating back and forth in the combustion chamber. Exposure to the pressure fluctuations over a long period of time can lead to a cataclysmic failure of engines. The underlying physics governing the generation of the thermo-acoustic instability is a complex interaction among heat release, turbulence, and acoustic waves. Currently, it is very difficult to accurately predict the expected level of oscillations in a combustor. Hence development of strategies and engineering solutions to mitigate thermo-acoustic instability is an active area of research in both academia and industry. In this work, we carry out numerical modeling of thermo-acoustic instability in a self-excited, laboratory scale, model rocket combustor developed at Purdue University. Two different turbulent combustion models to account for turbulence-chemistry interactions are considered in this study, the flamelet model and the transported probability density function (PDF) method. </div><div><br></div><div>In the flamelet modeling approach, detailed chemical kinetics can be easily incorporated at a relatively low cost in comparison to other turbulent combustion models and it also accounts for turbulence-chemistry interactions. The flamelet model study is divided into two parts. In first part, we examine the effect of different numerical approaches for implementing the flamelet model. In advanced modeling and simulations of turbulent combustion, the accuracy of model predictions is affected by physical model errors as well as errors that arise from the numerical implementation of models in simulation codes. Here we are mainly concerned with the effect of numerical implementation on model predictions of turbulent combustion. Particularly, we employ the flamelet/progress variable (FPV) model and examine the effect of various numerical approaches for the flamelet table integration, with presumed shapes of PDF, on the FPV modeling results. Three different presumed-PDF table integration approaches are examined in detail by employing different numerical integration strategies. The effect of the different presumed-PDF table integration approaches is examined on predictions of two real flames, a laboratory-scale turbulent free jet flame, Sandia Flame D and the self-excited resonance model rocket combustor. Significant difference is observed in the predictions both of the flames. The results in this study further support the claims made in previous studies that it is imperative to preserve the laminar flamelet structure during integration while using the flamelet model to achieve better predictions in simulations. In the second part of the flamelet modeling study, computational investigations of the coupling between the transient flame dynamics such as the ignition delay and local extinction and the thermo-acoustic instability developed in a self-excited resonance combustor to gain deep insights into the mechanisms of thermo-acoustic instability. A modeling framework that employs different flamelet models (the steady flamelet model and the flamelet/progress variable approach) is developed to enable the examination of the effect of the transient flame dynamics caused by the strong coupling of the turbulent mixing and finite-rate chemical kinetics on the occurrence of thermo-acoustic instability. The models are validated by using the available experimental data for the pressure signal. Parametric studies are performed to examine the effect of the occurrence of the transient flame dynamics, the effect of artificial amplification of the Damkohler number, and the effect of neglecting mixture fraction fluctuations on the predictions of the thermo-acoustic instability. The parametric studies reveal that the occurrence of transient flame dynamics has a strong influence on the onset of the thermo-acoustic instability. Further analysis is then conducted to localize the effect of a particular flame dynamic event, the ignition delay, on the thermo-acoustic instability. The reverse effect of the occurrence of the thermo-acoustic instability on the transient flame dynamics in the combustor is also investigated by examining the temporal evolution of the local flame events in conjunction with the pressure wave propagation. The above observed two-way coupling between the transient flame dynamics (the ignition delay) and the thermo-acoustic instability provides a plausible mechanism of the self-excited and sustained thermo-acoustic instability observed in the combustor.</div><div><br></div><div>The second turbulent combustion model considered in this study is the transported PDF method. The transported PDF method is one of the most attractive models because it treats the highly-nonlinear chemical reaction source term without a closure requirement and it is a generalized model for a wide range of turbulent combustion problems.</div><div>Traditionally, the transported PDF method has been used to model low-Mach number, incompressible flows where the pressure is assumed to be thermodynamically constant. Since there is significant pressure fluctuations in the model rocket combustor, the flow is highly compressible and it is necessary to account for this compressibility in the transported PDF method. In the past there has been very little work to model compressible reactive flows using the transported PDF and no effort has been made to model thermo-acoustic instability using the transported PDF method. There is a pressing need to further examine and develop the transported PDF method for compressible reactive flows to broaden our understanding of physical phenomenon like thermo-acoustic instability, interaction between combustion and strong shock and expansion waves, coupling between acoustic and heat release which are observed in high-speed turbulent combustion problems. To address this, a modeling framework for compressible turbulent reactive flows by the using the transported PDF method is developed. This framework is validated in a series of test cases ranging from pure mixing to a supersonic turbulent jet flame. The framework is then used to study the thermo-acoustic interactions in the self-excited model rocket combustor.</div>

Computational Fluid Dynamics (CFD) is one of the fields that has strongly developed since the recent development of faster computers and numerical modeling. CFD is also finding its way into chemical engineering on several levels. We have used CFD for detailed modeling of heat and mass transfer in a packed bed. One of the major questions in CFD modeling is whether the computer model describes reality well enough to consider it a reasonable alternative to data collection. For this assumption a validation of CFD data against experimental data is desired. We have developed a low tube to particle, structured model for this purpose. Data was gathered both with an experimental setup and with an identical CFD model. These data sets were then compared to validate the CFD results. Several aspects in creating the model and acquiring the data were emphasized. The final result in the simulation is dependent on mesh density (model detail) and iteration parameters. The iteration parameters were kept constant so they would not influence the method of solution. The model detail was investigated and optimized, too much detail delays the simulation unnecessarily and too little detail will distort the solution. The amount of data produced by the CFD simulations is enormous and needs to be reduced for interpretation. The method of data reduction was largely influenced by the experimental method. Data from the CFD simulations was compared to experimental data through radial temperature profiles in the gas phase collected directly above the packed bed. It was found that the CFD data and the experimental data show quantitatively as well as qualitatively comparable temperature profiles, with the used model detail. With several systematic variances explained CFD has shown to be an ample modeling tool for heat and mass transfer in low tube to particle (N) packed beds.

The DigIndy tunnel is an extension of the Indianapolis combined sewer system that stores the combined sewer overflow during heavy rain conditions. The tunnel system has several openings in and around the city of Indianapolis. Gasses emitted from the tunnel may create health concerns and affect the quality of life for nearby residents. Understanding the air circulation patterns provides valuable insight into where gases are likely to emerge from the tunnel and what steps may be taken to mitigate gas emissions in undesirable locations. The objective of the present work is to develop a computational fluid dynamics (CFD) model capable of predicting the air circulation patterns in the DigIndy tunnel under dry weather conditions. In order to inform and validate the CFD model, an experimental campaign was designed and executed to measure weather data and air flow rates within the DigIndy tunnel. Obtaining accurate results requires careful consideration of key physical phenomena to include in the model, geometric simplification strategies, mesh generation strategies, and numerical modeling strategies. Results showed that the seasonal effect, manifest by thermally-driven flow, plays a significant role in the air circulation patterns within the tunnel. Furthermore, results show that tunnel alignment affects the natural air circulation within the tunnel. Large diameter shafts, as the working and retrieval shafts, lead to significant circulation rates in the new tunnel alignments.

This work examines the detailed flow characteristics of direct measuring skin friction gages with computational methods. This type of device uses a small movable head mounted flush to a wall such that the head is assumed to be exposed to the same shear stress from the flow as the surrounding wall. The force caused by the action of the shear stress on the head deflects a flexure system monitored by instruments such as strain gages mounted at the base of a beam. The goal of the study was to develop an understanding of the effects that the geometric design and installation parameters of the sensor have on the surrounding flow and the ability of the sensor to reflect the undisturbed shear stress value. Disruption of the external flow due to poor design and/or improper installation of the sensor can take the form of intrusion into the flow, recession into the wall, and/or tilted alignment of the sensor such that the head is not flat in the plane of the wall, as well as flow into or out of the small gap surrounding the sensing head. Further, the performance of a direct measuring skin friction sensor in the presence of a pressure gradient has always been a concern. These effects are studied here with a three-dimensional, Navier-Stokes code based on a finite element method technique. Numerical solutions for cases in which one or more design parameters were varied are shown for a variety of flow situations. These situations include: (a) a laminar fully-developed channel flow at a low Reynolds number, (b) a turbulent flat plate boundary layer flow at a high Reynolds number, and (c) strong favorable and adverse pressure gradient turbulent boundary layer flows created by converging and diverging channels at high Reynolds number. Reported results for all cases include detailed flow visualization and stress field imagery, and total surface forces on the sensing head and gage flexure. Under ideal circumstances, these total forces should reflect as accurately as possible the average value of undisturbed shear stress times the exposed sensing head area (the friction force). Any deviation from this value was considered an "error" in the simulated measurement. The laminar channel flow case with a strong favorable pressure gradient showed the importance of proper alignment of the sensor. Protrusion or recession of the sensing head proved to be the dominant effect on resulting forces seen by the gage, changing the output by up to 15% for head protrusion and 10% for head recession for misalignments up to +/-1% of the head diameter. The thickness of the lip on the edge of the head also proved to have a significant effect on the output, with a smaller lip thickness generally showing better performance than a large one. Zero lip thickness indicated accuracy to within 1% of the desired wall shear result, since the pressure differences had little influence on the sensing head. Finally, the assumption of a linear pressure variation from the surface to the cavity along the lip as has been suggested in the past was investigated. The results indicate that the linear assumption works well only for large ratios of lip thickness to gap size, a fact which is correlated with previous experimental results. For the turbulent external flat plate case, misalignment remained the dominant effect on the sensor response. Results indicated that, in general, protrusion is more costly than the same level of recession, and a protrusion of +1% of the head diameter was shown to cause in excess of 100% error in indicated wall shear output. Both protrusion and recession produced large variations in both force and moment on the sensing flexure, but the outcome was that for protrusion the errors caused by these two effects tended to sum together, while for recession they tended to partially cancel out. The gap size played an increased role in the high Reynolds number boundary layer cases. Gap sizes of 1.67% up to 6.67% of the head diameter were studied and were shown to produce output errors between 4% and 22% (with larger errors corresponding to larger gap sizes), thus showing the importance of minimizing the gap for high Reynolds number flows. The lip was shown to have no significant effect for a flow without a pressure gradient. Finally, the favorable and adverse pressure gradient flows showed reasonable performance of the skin friction gage. Errors in output were shown to be -6% for the favorable pressure gradient case and 17% for the adverse pressure gradient case. Only the baseline gage design was studied for these situations, but the results from the two cases indicate that further reducing the lip thickness may not improve the performance of the gage. The error in output was caused almost entirely by applied moment for the adverse pressure gradient, while the applied force and applied moment had a cancellation effect in the favorable pressure gradient case. As a general result, the use of computational fluid dynamics has been shown to be an effective tool in the design and analysis of skin friction gages. Using a computational approach has the advantage of being able to resolve the small, confined gap regions of the gage, providing information that has been shown to be unavailable from previous experimental studies. This work has contributed to a much better understanding of the detailed flow over, in, and around skin friction gages. This will lead to improved gage design and reduced uncertainty in these important measurements. / Ph. D.

Master of Science / Department of Mechanical and Nuclear Engineering / Steven Eckels and Mohammad Hosni / This research presents results from experimental and numerical investigations of two-phase flow pattern analysis in a staggered tube bundle. Shell-side boiling tube bundles are used in a variety of industries from nuclear power plants to industrial evaporators. Fluid flow patterns in tube bundles affect pressure drop, boiling characteristics, and tube vibration. R-134a was the working fluid in both the experimental and computational fluid dynamics (CFD) analysis for this research. Smooth and enhanced staggered tube bundles were studied experimentally using a 1.167 pitch to diameter ratio. The experimental tube bundles and CFD geometry consist of 20 tubes with five tubes per pass. High speed video was recorded during the experimental bundle boiling. Bundle conditions ranged in mass fluxes from 10-35 kg/m[superscript]2.s and inlet qualities from 0-70% with a fixed heat flux. Classification of the flow patterns from these videos was performed using flow pattern definitions from literature. Examples of smooth and enhanced bundle boiling high speed videos are given through still images. The flow patterns are plotted and compared with an existing flow pattern map. Good agreement was found for the enhanced tube bundle while large discrepancies exist for the smooth tube bundle. The CFD simulations were performed without heat transfer with non-symmetrical boundary conditions at the side walls, simulating rectangular bundles used in this and other research. The two-phase volume of fluid method was used to construct vapor interfaces and measure vapor volume fraction. A probability density function technique was applied to the results to determine flow patterns from the simulations using statistical parameters. Flow patterns were plotted on an adiabatic flow pattern map from literature and excellent agreement is found between the two. The agreement between simulation results and experimental data from literature emphasizes the use of numerical techniques for tube bundle design.

Computational fluid dynamics and heat transfer (CFD/HT) models have been employed as the dominant technique for the design and optimization of both new and existing data centers. Inviscid modeling has shown great speed advantages over the full Navier-Stokes CFD/HT models (over 20 times faster), but is incapable of capturing the physics in the viscous regions of the domain. A coupled inviscid-viscous solution method (CIVSM) for bounded domains has been developed in order to increase both the solution speed and accuracy of CFD/HT models. The methodology consists of an iterative solution technique that divides the full domain into multiple regions consisting of at least one set of viscous, inviscid, and interface regions. The full steady, Reynolds-Averaged Navier-Stokes (RANS) equations with turbulence modeling are used to solve the viscous domain, while the inviscid domain is solved using the Euler equations. By combining the increased speed of the inviscid solver in the inviscid regions, along with the viscous solver’s ability to capture the turbulent flow physics in the viscous regions, a faster and potentially more accurate solution can be obtained for bounded domains that contain inviscid regions which encompass more than half of the domain, such as data centers.

Thesis (MEng)--Stellenbosch University, 2015. / ENGLISH ABSTRACT: The capability of computational fluid dynamics (CFD), in particular the FLUENT ™ commercial software suite, to predict wind loadings on heliostats has been investigated. If CFD proves useful in this area then the overall development costs of heliostats and concentrating solar thermal power plants could be reduced. Due to the largest loading on the heliostat originating from wind loads, by using CFD to determine these loads it could be possible to ensure heliostats are not overdesigned. This thesis contains a first study within the Solar Thermal Energy Research Group (STERG) at Stellenbosch University into the use of CFD for determining heliostat wind loads. The relevant theoretical background concerning the turbulence models used in this study, namely, the RNG k-ε, Realisable k-ε and SST k-ω turbulence models is reiterated. The „standard‟ k-ε model and the large eddy simulation (LES) approach, due to their relevance to bluff body flows, are also revisited. Some analysis is also provided around each model to gain insight as to the role of respective modelling sensitivities and their advantages. Previous work done in the area of heliostat wind studies is reviewed. The geometric considerations when dealing with heliostats leads onto the discussion concerning the requirement of modelling boundary layer profiles. Hence some background is provided on boundary layer modelling techniques. Further insight is drawn from more general previous bluff body CFD reported in the literature, from which observations and recommendations regarding the use of variations of the k-ε turbulence model can be inferred. The simulation procedure from geometry creation to results obtained for the flow over a vertical flat plate is reported. This investigation led to the conclusion that the Realisable k-ε should be used for the heliostat simulations on account of its accurate drag prediction under steady state flow conditions. It was also found that for transient simulations for heliostat like geometries, the SST k-ω model appears most suitable. The Realisable k-ε model is then used to model the flow about a heliostat using the same procedures as for the flat plate; both with flat and boundary layer inlet profiles. The overall conclusions drawn from this work are that the Realisable k-ε would not be suitable for predicting wind loads used in the final design of heliostats although it may be used with flat velocity and turbulence profiles to compare differences between early heliostat designs. The conclusion that the Realisable k-ε model should not be used to predict the flow field in the vicinity of a heliostat is also reached. It is recommended that further work should be carried out by using more advanced modelling techniques, such as the LES, to determine wind loads on heliostats. Furthermore, additional studies focused on accurately reproducing the velocity and turbulence profiles should be done. Lastly a larger set of data containing the orientations mentioned in literature should be generated using the methods contained within this study. / AFRIKAANSE OPSOMMING: Die vermoë van Numeriese Vloei Meganika (NVM), spesifiek die van die FLUENT ™ kommersiële sagtewarepakket, om die windlaste op heliostate te voorspel was ondersoek. As daar gevind word dat NVM wel betekinsvolle resultate kan lewer, kan dit die totale ontwikkelingskoste van heliostate en gekonsentreerdesonkragstasies verlaag. Wind plaas die grootste las op heliostate, dus deur gebruik te maak van NVM om die windlaste op heliostate te voorspel, kan dit gebruik word om te verseker dat heliostate nie oorontwerp word nie. Hierdie tesis bevat „n eerste studie binne die Sontermiese Energie Navorsings Groep aan die Universiteit van Stellenbosch, wat die gebruik van NVM om windlaste op heliostate te voorspel ondersoek. Alle relevante teoretiese agtergrond wat turbulensiemodelle aanbetref, naamlik die RNG k-ε, Realiseerbare k-ε en SST k-ω turbulensiemodelle, word bespreek. Hulle relevansie tot stompligaamvloei toegestaan, word die „standaard‟ k-ε model en die groot werwel simulasie (GWS) benaderings ook bespreek. Elke model word bespreek om die leser insig te gee in dié model se sensitiwiteite en voordele. Vorige studies wat betrekking het tot die studie van heliostate en wind word bespreek. Die geometrie van heliostate lei tot „n bespreking oor die noodsaklikheid vir „n model vir die grenslaagprofiel, dus word grenslaagmodelleringstegnieke bespreek. Verdere insig word verkry van vorige NVM studies uit die literatuur met meer algemene stomp liggame, wat waarnemings en voorstelle vir die gebruik van die k-ε turbulensiemodel en variante verskaf. Die simulasieproses, vanaf geometrieskepping tot die resultate vir die vloei oor 'n vertikale vlak, word bespreek. Hierdie ondersoek het tot die gevolgtrekking gelei dat die realiseerbare k-ε model gebruik moet word vir die heliostaat simulasies, as gevolg van die akkurate sleurvoorspellings onder bestendigetoestande. Daar was ook gevind dat vir heliostaatagtige liggame onder oorgangskondisies, die SST k-ω model mees geskik sal wees. Die Realiseerbare k-ε model word dan gebruik om die vloei om 'n heliostaat te modelleer deur gebruik te maak van dieselfde proses wat gebruik word om vloei oor 'n plat plaat te analiseer: albei met plat en grenslaaginlaatprofiele. Die gevolgtrekkings van hierdie studie is dat die Realiseerbare k-ε model nie gebruik kan word tydens die finale ontwerpfase om die windlaste op 'n heliostaat te voorspel nie. Dit kan wel gebruik word met plat snelheids- en turbulensieprofile om die versikille tussen vroeë heliostaatkonsepte te vergelyk. Daar was ook bepaal dat die Realiseerbare k-ε model nie gebruik moet word om die vloeiveld om 'n heliostaat te voorspel nie. Daar word voorgestel dat verdere studies in hierdie vakgebied met meer gevorderde modelleringstegnieke aangepak word. Dit word aanbeveel dat verdere werk uitgevoer moet word deur die gebruik van meer gevorderde modellering tegnieke, soos GWS, om die wind kragte op heliostats te bepaal. Verder, studies wat akkurate snelheid en turbulensieprofiele produseer sal nog bygelas moet word. Laastens 'n groter stel data met oriëntasies soos wat in die literatuur beskryf word, moet deur middel van die metodes van dié studie gegenereer word.

Thesis (MEng)--Stellenbosch University, 2015. / ENGLISH ABSTRACT: This thesis investigates the flow field and pressure distributions around cylinders at trans-critical Reynolds numbers using the k-ε Realizable turbulence model. A steady state 2-D and 3-D Fluent® model is successfully developed to evaluate the effects of changing various modelling parameters on the static pressure distribution around an infinite and finite cylinder. These parameters include surface roughness, cylinder rotation and air viscosity at the cylinder surface. The subsequent results obtained are compared to each other and to data trends from literature as well as measured experimental results and are found to be in good agreement. In addition a method for calibrating all developed methods based on their shear stress curves over a flat plate model is also successfully developed. The main objective is to find an appropriate single parameter which can be used for the rigorous adjustment of the pressure distribution around a cooling tower, which will allow for improved sensitivity analysis and modelling of cooling tower performance under wind conditions with and without meridional ribs located on the outer shell surface. / AFRIKAANSE OPSOMMING: Hierdie tesis ondersoek die vloeiveld en druk verdelings rondom silinders by trans-kritiese Reynolds getalle deur gebruik te maak van die k-ε Realizable turbulensie model. ‘n Bestendige toestand 2-D en 3-D Fluent® model is suksesvol ontwikkel om die uitwerking van die verandering van verskeie model parameters op die statiese druk verdeling rondom ‘n oneindige en eindige silinder te evalueer. Die laasgenoemde parameters sluit in oppervlak grofheid, silinder rotasie en lug viskositeit by die silinder wand. Die daaropeenvolgende resultate wat verkry word, word met data tendense uit die literatuur asook gemete data vanuit eksperimente vergelyk en goeie ooreenkoms i.t.v die data tendense is gevind. Verder is ‘n metode vir die suksesvolle kalibrasie van die ontwikkelde numeriese tegnieke ontwikkel. Die laasgenoemde kalibrasie metode is gebaseer op die vergelyking van skuifspanning kurwes vir vloei oor ‘n plat plaat model. Die hoofdoel van die navorsing is om ‘n geskikte enkele parameter te vind wat gebruik kan word vir die effektiewe aanpassing van die druk verdeling rondom ‘n koeltoring wat sal lei tot verbeterde sensitiwiteits analise en modellering van koeltoring verrigting onder wind toestande met en sonder meridionale ribbes geleë op die buitenste dop oppervlak.