Magnesium and magnesium alloys are used in a wide variety of structural applications including automotive, aerospace, hand tools and electronic industries thanks to their light weight, high specific strength, adequate corrosion resistance and good castability. Al and Zn are the primary alloying elements in commercial Mg alloys and commonly used in automotive industries. AZ91 is one of the most popular Mg alloys containing 9% Al and 1% Zn. Hence, lots of research have been done during last decades on AZ91D. However, the existing data concerning mechanical properties and microstructural features showed large scatter and is even contradictory. This work focused on the correlation between the microstructure and the mechanical properties of as-cast AZ91 alloy. An exhaustive characterization of the grain size, secondary dendrite arm spacing (SDAS) distribution, and fraction of Mg17Al12 using optical and electron backscattered diffraction (EBSD) was performed. These microstructural parameters were correlated to offset yield point (Rp0.2), fracture strength and elongation to fracture. It was understood that the intermetallic phase, Mg17Al12, plays an important role in determining the mechanical and physical properties of the alloy at temperature range from room temperature up to 190oC. It was realized that by increasing the Mg17Al12 content above 11% a network of intermetallic may form. During deformation this rigid network should break before any plastic deformation happen. Hence, increase in Mg17Al12 content resulted in an increase in offset yield point. The presence of this network was supported by study of thermal expansion behaviour of the alloy containing different amount of Mg17Al12. A physically-based model was adapted and validated in order to predict the flow stress behaviour of as-cast AZ91D at room temperature up to 190ºC for various microstructures. The model was based on dislocation glide and climb in a single-phase (matrix) material containing reinforcing particles. The temperature dependant variables of the model were quite well correlated to the underlying physics of the material.

The sustainable development of the offshore wind and wave energy sectors requires optimising the exploitation of the resources, and it is in relation to this and the shared challenge for both industries to reduce their costs that the option of integrating offshore wind and wave energy arose during the past decade. The relevant aspects of this integration are addressed in this work, and in particular the evaluation of the Wave Energy Converter (WEC) sub-system of a hybrid wind-wave energy converter: the state of the art of combined technologies; the definition of a novel hybrid prototype, based on a preliminary feasibility analysis of a conceptual proposal; and the evaluation of a simplified version of this prototype by means of physical and numerical modelling as a mean to set the reference and define new tools and methods for future evaluation and optimisation of the prototype. Because of the novelty of combined wave and offshore wind systems, fundamental knowledge was lacking as, for example a comprehensive review and classification, which was published as a journal paper framed in the present work. In particular, the core of this PhD thesis deals with the WEC sub-system of a hybrid device that integrates an Oscillating Water Column (OWC) device into the typical monopile substructure of an offshore wind turbine. A new prototype of the hybrid energy converter has been proposed, and a patent application was filled. Furthermore, an experimental set-up was designed, built and tested at a wave flume. On the basis of this experimental campaign the performance of the device is analysed. Finally, a full 3D-numerical mirror of the experimental set-up, including the hybrid energy converter, is defined and validated, and the flume enclosure effects studied for regular waves.

Urbanisation will require significant expansion of underground infrastructure, which results in unavoidable ground displacements that affect the built environment. Predicting the interaction between a tunnel, the soil and existing structures remains an engineering challenge due to the highly non-linear behaviour of both the soil and the building. This thesis investigates the interaction between a surface structure and tunnelling-induced ground displacements. Specifically, novel three-dimensionally printed building models with brittle material behaviour, similar to masonry, were developed and tested in a geotechnical centrifuge. This enabled replication of building models with representative global stiffness values and realistic building features including strip footings, intermediate walls, a rough soil-structure interface, building layouts and façade openings. By varying building characteristics, the impact of structural features on both the soil and building response to tunnelling in dense sand was investigated. Results illustrate that the presence of surface structures considerably altered the tunnelling-induced soil response. The building-to-tunnel position notably influences the magnitude of soil displacements and causes localised phenomena such as embedment of building corners. An increase of the façade opening area and building length reduces the alteration of the theoretical greenfield settlements, in particular the trough width. Moreover, the impact of varying the building layout is discussed in detail. For several building-tunnel scenarios, building distortions are quantified and the crucial role of building features is demonstrated. Structures spanning the greenfield inflection point experienced more deformation than identical structures positioned in either sagging or hogging, and partitioning a structure either side of the greenfield inflection point is shown to lead to unconservative damage assessments. Results also quantify the significant extent to which structural distortions increase as façade openings and building length increases. Observed building damage and cracking patterns confirm the reported trends. The experimental results are used to evaluate the performance of available methods to assess the behaviour of buildings to tunnelling. Predictions ignoring soil-structure interaction are usually overly conservative, while approaches based on the relative stiffness of a structure and the soil result in inconsistent predictions, though some methods performed better than others. Practical improvements to consider structural details when assessing this tunnel-soil-structure system are finally proposed.

The original aim of the research project was to investigate the mechanism of power capture from sea waves and to optimise the performance of a vee-shaped floating Wave Energy Converter, the Floating Clam, patented by Francis Farley. His patent was based on the use of a pressurised bag (or ‘reservoir’) to hold the hinged Clam sides apart, so that, as they moved under the action of sea waves, air would be pumped into and out of a further air reservoir via a turbine/generator set, in order to extract power from the system. Such “Clam Action” would result in the lengthening of the resonant period in heave. The flexibility of the air bag supporting the Clam sides was an important design parameter. This was expected to lead to a reduction in the mass (and hence cost) of the Clam as compared with a rigid body. However, the present research has led to the conclusion that the Clam is most effective when constrained in heave and an alternative power take-off is proposed. The theoretical investigations made use of WAMIT, an industry-standard software tool that provides an analysis based on potential flow theory where fluid viscosity is ignored. The WAMIT option of Generalised Modes has been used to model the Clam action. The hydrodynamic coefficients, calculated by WAMIT, have been curve-fitted so that the correct values are available for any chosen wave period. Two bespoke mathematical models have been developed in this work: a frequency domain model, that uses the hydrodynamic coefficients calculated by WAMIT, and a time domain model, linked to the frequency domain model in such a way as to automatically use the same hydrodynamic and hydrostatic data. In addition to modelling regular waves, the time domain model contains an approximate, but most effective method to simulate the behaviour of the Clam in irregular waves, which could be of use in future control system studies. A comprehensive series of wave tank trials has been completed, and vital to their success has been the modification of the wave tank model to achieve very low values of power take-off stiffness through the use of constant force springs, with negligible mechanical friction in the hinge mechanism. Furthermore, the wave tank model has demonstrated its robustness and thus its suitability for use in further test programmes. The thesis concludes with design suggestions for a full-scale device that employs a pulley/counterbalance arrangement to provide a direct connection to turbine/generator sets, giving an efficient drive with low stiffness and inherently very low friction losses. At the current stage of research, the mean annual power capture is estimated as 157.5 kW, wave to wire in a far from energetic 18 kw/m mean annual wave climate, but with scope for improvement, including through control system development.

Le domaine du Térahertz a suscité beaucoup d’intérêt de la part de la communauté scientifique ces dernières années. La diode Schottky constitue la pierre angulaire des circuits de détection, des multiplieurs de fréquence ou encore des mélangeurs dans cette bande de fréquence, notamment grâce à son comportement non-linéaire. Les travaux menés durant cette thèse visent à optimiser les caractéristiques de ce composant pour deux fonctions non-linéaires – la détection et la multiplication de fréquence – celles-ci ayant des facteurs de mérites bien spécifiques. Un même dispositif ne saurait les satisfaire conjointement La non-linéarité capacitive est généralement mise à profit pour la multiplication de fréquence alors que la détection s’appuie sur la non-linéarité résistive associée à la caractéristique statique. Pour réaliser cette optimisation, un code particulaire Monte-Carlo (MC) résolvant l’équation de Boltzmann couplée à l’équation de Poisson a été développé. La diode Schottky est un composant largement contrôlé par l’interface métal/semi-conducteur et la gestion des conditions aux limites constitue une étape clef dans la modélisation du dispositif. Le principe d’exclusion de Pauli doit être considéré pour un semi-conducteur très dopé et une distribution spécifique pour les porteurs injectés du côté du contact ohmique a été optimisée puis utilisée dans la modélisation de la diode. D’autres effets physiques à l’interface métal/semi-conducteur ont été implémentés tels que l’effet tunnel suivant différents degrés de raffinement, le phénomène de force image ou encore l’abaissement de la hauteur de barrière par le champ électrique dû aux états de surface. Cette modélisation MC a permis de déduire un schéma équivalent électrique petit-signal aux fréquences Térahertz dont les différents paramètres sont ajustés en prenant en compte la déplétion possible du substrat pour des diodes courtes. L’extraction du schéma équivalent peut s’effectuer suivant différentes stratégies : en excitant la diode avec un signal de faible amplitude ou encore à partir de l’étude des densités spectrales associées aux fluctuations de courant et de tension. Les phénomènes physiques pouvant mettre en défaut ce schéma électrique tels que la vitesse de saturation des porteurs ou le phénomène d’ionisation par choc en polarisation inverse sont discutés. Le recours à un schéma électrique est motivé par une volonté de l’intégrer facilement au cœur d’un circuit pour une fonction spécifique et de l’exploiter avec un logiciel commercial tel que ADS (Advanced Design System) dans une logique d’optimisation. Des simulations de type « Harmonic Balance » ont été menées afin d’étudier le rendement d’un circuit de détection et d’un multiplieur de fréquence pour dégager les caractéristiques optimales de la diode sur chacun de ces circuits. Le GaAs est souvent un semi-conducteur de choix pour la réalisation de circuits aux fréquences Térahertz grâce à sa maturité technologique et à sa haute mobilité électronique. D’autres semi-conducteurs tels que l’InGaAs, le GaSb ou encore le GaN sont également étudiés. Une diode avec un couple métal/semi-conducteur présentant une faible hauteur de barrière donne les meilleurs rendements de conversion pour la détection. Pour le multiplieur de fréquence, il existe un dopage optimal en fonction de la longueur de la couche active permettant de maximiser le rendement du circuit. / The terahertz field has generated significant interest from the scientific community in recent years. The Schottky diode is the cornerstone of detection circuits, frequency multipliers or mixers in this frequency band, thanks to its non-linear behavior. The work carried out during this thesis aims at optimizing the characteristics of this component for two non-linear functions - detection and frequency multiplication - these having very specific figures of merit. The same device can not satisfy them together. Capacitive non-linearity is generally used for frequency multiplication while the detection is based on the resistive non-linearity associated with the current-voltage characteristic. To achieve this optimization, a Monte Carlo (MC) particle code solving the Boltzmann transport equation coupled to the Poisson equation was developed. The Schottky diode is a component essentially controlled by the metal / semiconductor interface and the boundary conditions are a key step in the device modeling. The Pauli exclusion principle must be considered for a highly doped semiconductor and a specific distribution for the carriers injected on the ohmic contact side has been optimized and then used in the diode modeling. Other physical effects at the metal / semiconductor interface have been implemented, such as tunneling following different degrees of refinement, image force barrier lowering or the barrier height lowering by the electric field due to surface states. A small-signal equivalent circuit at Terahertz frequencies was deduced from the MC modeling. The parameters of this circuit are adjusted by taking into account the possible depletion of the substrate for short diodes. The extraction of the equivalent circuit can be carried out according to different strategies: exciting the diode with a low amplitude signal or studying the spectral densities associated with the current and voltage fluctuations. The physical phenomena that can make the equivalent circuit defective, such as the saturation velocity of carriers or the impact ionization at reverse bias, are discussed. An electrical circuit of the diode was used to be easily integrated into a global circuit for a specific function and exploit it in commercial software such as ADS (Advanced Design System). Harmonic balance simulations were conducted to study the performance of a detection circuit and a frequency multiplier to reach the optimal characteristics of the diode on each of these circuits. GaAs is often a semiconductor of choice for the design of circuits at Terahertz frequencies thanks to its technological maturity and its high electronic mobility. Other semiconductors such as InGaAs, GaSb or GaN are also investigated. A diode with a metal / semiconductor couple having a low barrier height gives the best conversion efficiencies for detection. For the frequency multiplier, there is an optimal doping as a function of the active layer length which maximizes the circuit efficiency.

En dépit d'une amélioration remarquable de la performance des composants électroniques organiques, il y a encore un manque de compréhension théorique rigoureux sur le fonctionnement du composant. Cette thèse est consacrée à la création de modèles pratiques pour composants électroniques organiques à base physique complet, à savoir un modèle compact à base physique. Un modèle compact à base physique d'un élément de circuit est une équation mathématique qui décrit le fonctionnement du composant, et est généralement évaluée par trois critères: si elle est suffisamment simple pour être incorporé dans des simulateurs de circuits, précise pour rendre le résultat des simulateurs utile les concepteurs de circuits et rigoureux pour capturer des phénomènes physiques se produisant dans le composant. Dans ce contexte, les caractéristiques distinctives de l'injection de porteurs de charge et de transport dans les semi-conducteurs organiques sont incorporés dans les modèles avec un effort particulier pour maintenir la simplicité mathématique. L'effet concomitant sur les caractéristiques courant-tension des diodes et des transistors organiques prototypiques sont étudiés. Les méthodes d'extraction des paramètres cohérents aux modèles sont présentés qui permettent la détermination univoque des paramètres de le composant utilisé pour le fonctionnement du composant de modélisation et l'évaluation des performances de le composant et les propriétés des couches minces et des interfaces organiques. Les approches englobent le developement analytique des équations physiques, la simulation numérique à deux dimensions basé sur la méthode des éléments finis et la validation expérimentale. Les modèles compacts originaux et entièrement analytiques et des méthodes d'extraction de paramètres fournissent une compréhension fondamentale sur la façon dont le désordre énergétique dans une couche mince de semi-conducteur organique, décrit par la densité d’etats Gaussienne, affecte les caractéristiques courant-tension observables des composants.Mots-clés : Electronique organique, physique des composants électroniques, modélisation analytique, diodes, transistors à effet de champ, densité d’etats Gaussienne / In spite of a remarkable improvement in the performance of organic electronic devices, there is still a lack of rigorous theoretical understanding on the device operation. This thesis is dedicated to establishing practical models of organic electronic devices with a full physical basis, namely a physically-based compact model. A physically-based compact model of a circuit element is a mathematical equation that describes the device operation, and is generally assessed by three criteria: whether it is sufficiently simple to be incorporated in circuit simulators, accurate to make the outcome of the simulators useful to circuit designers, and rigorous to capture physical phenomena occuring in the device. In this context, distinctive features of charge carrier injection and transport in organic semiconductors are incorporated in the models with a particular effort to maintain mathematical simplicity. The concomitant effect on the current-voltage characteristics of prototypical organic diodes and transistors are studied. Parameter extraction methods consistent to the models are presented which enable unambiguity determination of device parameters used for modeling device operation and assessing device performance and properties of organic thin-films and interfaces. The approaches encompass analytical developement of physical equations, two-dimensional numerical simulation based on finite-element method and experimental validation. The original and fully analytical compact models and parameter extraction methods provide fundamental understanding on how energetic disorder in an organic semiconductor thin-film, described by the Gaussian density of states, affects the observable current-voltage characteristics of the devices.Keywords : Organic electronics, device physics, analytical modeling, diodes, field-effect transistors, Gaussian density-of-states

This study investigates the relationship between wave-induced hydrodynamics and the resulting loading on Core-Loc concrete armour units below the still water level in a breakwater structure. Physical modelling experiments were performed at the National Research Council in Ottawa in which a 3D-printed 12 cm Core-Loc armour unit was instrumented and fixed in place within a rubble mound structure. Testing featured simultaneous measurement of force on this instrumented unit, pressure head at the base of the unit, and flow velocities below the SWL. Two main scenarios were tested, the isolated unit and fully armoured scenarios, under a range of regular waves and irregular sea states. Analysis of force development on the instrumented unit indicates that maximum slope-normal forces (both into and away from the structure) are associated with extremes in pressure head above the instrumented unit, while slope-parallel force extremes (both upslope and downslope) occur at times of the fastest change in water level. These loadings are consistent with Morison’s equation and imply drag dominance in the slope-parallel direction and inertia dominance in the slope-normal direction. Significant differences in forces were observed between isolated (no neighbouring units) and embedded (with neighbouring units) armour unit test cases. The presence of the armour layer significantly increased the normal force exerted on the unit and reduced the parallel force. Irregular sea state testing shows force peaks following normal distribution. Analysis of flow above the armour layer showed that force, flow velocity and flow acceleration are symmetrical in the slope-parallel direction, but largely asymmetrical in the slope-normal direction, with the flow velocity and force on the unit in particular experiencing large asymmetries. Wave height analysis indicated that each wave height follows a similar force development pattern with a magnitude proportional to wave height. Wave period analysis showed the formation of small secondary waves as the period increases. Wave steepness affected the peak force loading of the instrumented unit in a mostly linear fashion up to the critical Iribarren number.

This thesis work aimed at increasing the knowledge regarding phenomena occurring when gas is injected using a top-blown lance on to a bath. All results are based on physical modeling studies carried out both using low and high gas flow rates and nozzle diameters ranging from 0.8 mm to 3.0 mm. At the low gas flow rates, the penetration depth in the bath was studied. The experiments focused on studying the effect of nozzle diameter, lance height and gas flow rate on the penetration depth. It was found that the penetration depth increases with decreasing nozzle diameter, decreasing lance height and with increasing gas flow rate. The results were also compared with previous work. More specifically, it was studied how the previous published empirical relationships fitted the current experimental data. It was found that the relationships of Banks [1], Davenport [2], Chatterjee [3] and Qian [4] agreed well with the experimental data of this investigation for nozzle diameters of 2.0 mm and 3.0 mm. However, for smaller nozzle diameters there were considerable deviations. Therefore, a new correlation heuristically derived from energy conservation consideration was suggested and showed better agreement for small nozzle diameters. The experiments carried out at higher gas flow rates focused on the study of swirl motion. The effects of nozzle diameter, lance height, gas flow rate and aspect ratio on the swirl motion were investigated. The amplitude and period of the swirl as well as the starting time and the damping time of the swirl were determined. The amplitude was found to increase with an increased nozzle diameter and gas flow rate, while the period had a constant value of about 0.5 s for all nozzle diameters, gas flow rates and lance heights. The starting time for the swirl motion was found to decrease with an increased gas flow, while the damping time was found to be independent of gas flow rate, nozzle diameter, lance height and ratio of depth to diameter. / QC 20101217

Thesis (MEng)--Stellenbosch University, 2014. / ENGLISH ABSTRACT: In modern steelmaking a tundish serves as an important metallurgical reactor to remove inclusions and maintain thermal and chemical homogeneity in the product. In this study the flow behaviour in a four strand tundish was investigated by means of a 1/2 scale water model, based on Froude number similarity, as well as by using numerical modelling. Both the numerical study and physical model were used to characterise residence time distribution (RTD) in the vessel and to calculate properties pertaining to the tundish flow regime. The three different tundish configurations investigated in this study are: a bare tundish with no flow control devices, a tundish with a turbulence inhibitor and a tundish using a turbulence inhibitor with holes in combination with dams. Preliminary investigations focussed on the framework for obtaining an accurate numerical solution within reasonable computational times. The effect of assuming symmetry and dynamically steady flow in the numerical model was shown to be small relative to the effect of grid size and justifiable by the savings in computational time. The grid independence study indicated the importance of using a finer mesh in areas of high velocity gradients to obtain realistic results and also to limit the number of computational cells. A procedure using gradient adaptation was used to refine the mesh automatically in the required regions for different tundish geometries. Results also showed that the inlet boundary of the numerical model should be selected at the ladle outlet, since assuming a flat velocity profile at the nozzle port resulted in significant changes in the RTD response. Comparison of the results obtained using the numerical model with those from physical experiments yielded an average error of less than 10%. This was assumed to be a good prediction, considering the assumptions employed in the numerical model. Both the physical and numerical models showed that a tundish without flow control devices was prone to significant short circuiting. The addition of a turbulence inhibitor was shown to be successful in preventing short circuiting and provided surface directed flow, which is thought to aid inclusion removal in the slag. Additionally, the minimum, peak and mean residence times and plug flow volume fraction were increased significantly, while the dead volume fraction decreased. However, using a turbulence inhibitor with holes in combination with dams showed that this configuration may cause increased refractory wear together with increased risk of slag entrainment due to flow patterns with increased surface turbulence. It also showed that the short-circuiting might not be eliminated completely. This indicates that certain design changes to tundish flow control systems can introduce problems that outweigh the benefits of the altered flow patterns. Furthermore, the numerical method, which was based on the water model, was modified to simulate the high temperature steel process. A very good match was obtained between the results using the two different numerical models. This serves as additional evidence that tundish water modelling based on Froude number similarity provides a good representation of the actual industrial process. Using the numerical model based on the high temperature steel process the effect of turbulence inhibitor shape was studied for four different turbulence inhibitor designs. Results showed the best performance, based on flow characteristic properties and surface turbulence values, was achieved for the design using a rectangular box-like shape with flanges at the top. However, the comparison emphasized the effect of the turbulence inhibitor shape on the flow behaviour, as each design yielded completely different flow patterns. It was also observed that a good turbulence inhibitor provided an optimum amount of turbulent suppression. Insufficient suppression would cause fast flows, which will result in insufficient residence time for inclusion flotation and high surface turbulence values, which may cause slag entrainment. On the other hand, too much suppression may increase the variation between strands. / AFRIKAANSE OPSOMMING: Die verdeeltrog speel ‘n belangrike rol in die moderne staalvervaardigingsproses deur inklusies te verwyder en termiese en chemiese homogeniteit in die produk te verseker. In hierdie studie is die vloeigedrag in ‘n verdeeltrog met vier uitlate bestudeer deur middel van ‘n 1:2-skaal watermodel, gebaseer op Froude-getal gelykheid, sowel as deur numeriese modellering. Beide die numeriese en watermodel is gebruik om die verblyftyd-distribusie in die trog te karakteriseer en om waardes te bereken wat die vloeigedrag in die verdeeltrog beskryf. Drie verskillende verdeeltrog-opstellings is in hierdie studie bestudeer, naamlik: ‘n leë verdeeltrog met geen vloeibeheertoestelle nie, ‘n verdeeltrog met ‘n turbulensie- inhibeerder en ‘n verdeeltrog wat gebruik maak van ‘n turbulensie-inhibeerder met gate, gekombineer met lae damwalle. Vroeë ondersoeke het gefokus op ‘n metode om akkurate numeriese resultate binne aanvaarbare tye te verkry. Die invloed van die aannames van simmetrie en dinamiese gestadigde vloei op die resultate is bepaal om klein te wees teenoor dié van die roostergrootte. Die gebruik van die aannames is dus geregverdig deur die afname in berekeningstyd wat dit meebring. Die roosteronafhanklikheidstudie het getoon dat dit belangrik is om die vloeivergelykings oor klein volume-eenhede op te los in areas van hoë snelheidsgradiënte, beide om realistiese resultate te verseker en om die aantal volume-eenhede te beperk. ‘n Prosedure wat gebruik maak van gradiëntaanpassing is gebruik om die roostergrootte outomaties te verklein in die areas met hoë snelheidsgradiënte vir verskillende verdeeltrog-opstellings. Resultate het ook getoon dat dit belangrik is om die inlaatgrens van die numeriese model by die smeltpot se uitlaat te kies, aangesien noemenswaardige verskille in die verblyftyd-distribusie waargeneem is wanneer ‘n uniforme snelheidsprofiel aanvaar is waar die metaal die verdeeltrog binnegaan. Daar is gevind dat die resultate by die numeriese en watermodelle verskil met ‘n gemiddelde fout van ongeveer 10%. Dit word beskou as ‘n goeie ooreenstemming, as die aannames wat in die numeriese model gebruik word in ag geneem word. Beide die numeriese en watermodelle het getoon dat die metaalvloei in die leë verdeeltrog geneig is om die kortste roete na die uitlate te kies. Deur van ‘n turbulensie inhibeerder gebruik te maak word hierdie probleem opgelos deur die metaalvloei opwaarts te forseer. In die proses neem die minimum-, piek- en gemiddelde verblyftye, sowel as die propvloei volumefraksie, toe, terwyl die dooievloei fraksie afneem. Die turbulensie inhibeerder met die gate en die lae damwalle is egter meer geneig tot erosie van die vuurvaste materiaal en kan ook vloeipatrone veroorsaak wat metaalskuimdruppels in die metaal kan opneem. Verder volg ‘n klein fraksie van die vloei steeds ‘n kortpad na die uitlaat. Hierdie resultate toon dat sekere veranderings aan verdeeltrog vloeibeerheerstelsels nadelige gevolge kan hê wat enige positiewe gevolge uitkanselleer. Verder is die numeriese metode wat tot dusver op die watersisteem gebaseer is, aangepas om die vloei in die hoëtemperatuur staalproses te simuleer. ‘n Baie goeie vergelyking is verkry tussen die resultate wat die twee numeriese modelle voorspel. Hierdie dien as ‘n verdere bewys dat ‘n watermodel, wat op Froude-getal gelykheid gebaseer is, die industriële proses akkuraat kan simuleer. Deur gebruik te maak van die numeriese model van die hoë temperatuur staalproses, is die effek van die turbulensie-inhibeerder se vorm vir vier verskillende ontwerpe bestudeer. Die studie het getoon dat die beste resultate vir vloeikarakteriserende eienskappe en oppervlak- turbulensie verkry is met ‘n reghoekige ontwerp. Die resultate beklemtoon egter die belangrikheid van die invloed van die turbulensie-inhibeerder se vorm op die vloeigedrag in die verdeeltrog, aangesien elke vorm noemenswaardige verskillende vloeipatrone opgelewer het. Daar is opgemerk dat ‘n goeie turbulensie-inhibeerder ‘n optimale hoeveelheid turbulensie onderdrukking veroorsaak. In die geval van te min onderdrukking is die verblyftye te kort en die oppervlak turbulensie te hoog. Te veel onderdrukking kan egter lei tot groot variasie in die eienskappe van die metaal by die verskillende uitlate.

<p>This licentiate thesis describes development of modellingtools, experimental physical modelling and numerical modellingto simulate real combustion processes for advanced industrialutility boiler before and after retrofit.</p><p>The work presents extended study about formation,destruction and control of pollutants, especially NOx, whichoccur during combustion process.</p><p>The main aim of this work is to improve mixing process incombustion chamber. To do this, the optimization of placementand direction of additional air and fuel nozzles, the physicalmodelling technique is used. By using that method, it ispossible to obtain qualitative information about processes,which occur in the real boiler. The numerical simulationsverify the results from physical modelling, because duringmathematical modelling quantitative informations about flow andmixing patterns, temperature field, species concentration areobtained.</p><p>Two 3D cases, before and after retrofit, of pulverized fuelfired boiler at 125 MW output thermal power are simulated. Theunstructured mesh technique is also used to discretize theboiler. The number of grid was 427 656 before retrofit and 513362 after retrofit. The comparisons of results of numericalsimulation before and after retrofit are presented. The resultsfrom physical modelling and numerical simulation are alsoshown.</p><p>Results present that nozzles of additional air and fuel givea considerably better mixing process, uniform temperature fieldand CO2 mass fraction. The whole combustion chamber worksalmost as a "well stirred reactor", while upper part of boilerworks as a "plug flow reactor".</p><p>Differences between from measured of temperatures andpredicted temperatures are not too big, the maximum differenceis about 100 K. It seems, that calculated temperatures showgood agreement with measurement data.</p><p>The results illuminate the potential of physical andnumerical modelling methods as promising tools to deal with thecomplicated combustion processes, even for practicalapplication in the industry.</p><p><b>Keywords:</b>air staging, fuel staging, boiler, furnace,computational fluid dynamics, numerical simulation, pollutants,physical modeling, pulverized fuel combustion.</p>