The mass and thermal transport in porous media play an important role in many engineering and geological processes. The hydrodynamic and thermal effects are two interesting aspects arising in the research of porous media. This thesis is concerned with numerical investigations of double-diffusive natural convective heat and mass transfer in saturated porous cavities with Soret and Dufour effects. An in-house FORTRAN code, named ALFARHANY, was developed for this study. The Darcy-Brinkman-Forchheimer (generalized) model with the Boussinesq approximation is used to solve the governing equations. In general, for high porosity (more than 0.6), Darcy law is not valid and the effects of inertia and viscosity force should be taken into account. Therefore, the generalized model is extremely suitable in describing all kinds of fluid flow in a porous medium. The numerical model adopted is based on the finite volume approach and the pressure velocity coupling is treated using the SIMPLE/SIMPLER algorithm as well as the alternating direction implicit (ADI) method was employed to solve the energy and species equations. Firstly, the model validation is accomplished through a comparison of the numerical solution with the reliable experimental, analytical/computational studies available in the literature. Additionally, transient conjugate natural convective heat transfer in two-dimensional porous square domain with finite wall thickness is investigated numerically. After that the effect of variable thermal conductivity and porosity investigated numerically for steady conjugate double-diffusive natural convective heat and mass transfer in two-dimensional variable porosity layer sandwiched between two walls. Then the work is extended to include the geometric effects. The results presented for two different studies (square and rectangular cavities) with the effect of inclination angle. Finally, the work is extended to include the Soret and Dufour effects on double-diffusive natural convection heat and mass transfer in a square porous cavity. In general, the results are presented over wide range of non-dimensional parameters including: the modified Rayleigh number (100 ≤ Ra* ≤ 1000), the Darcy number (10-6 ≤ Da ≤ 10-2), the Lewis number (0.1 ≤ Le ≤ 20), the buoyancy ratio (-5 ≤ N ≤ 5), the thermal conductivity ratio (0.1 ≤ Kr ≤ 10), the ratio of wall thickness to its height (0.1 ≤ D ≤ 0.4), the Soret parameter (-5 ≤ Sr ≤ 5), and the Dufour parameter (-2 ≤ Df ≤ 2).

In the present thesis turbulent natural convection of air within different challenging test cases are investigated numerically by means of an unstructured finite volume code, Code_Saturne. First, flow within both two-dimensional vertical and inclined differentially heated rectangular cavities at 60° and 15° to the horizontal for an aspect ratio of H/L=28.6 and Rayleigh number of 0.86×10e6 is computed using several high and low-Re models. Here the effectiveness of the RANS models in Code_Saturne is assessed through comparisons with a range of available experimental data. After some tests of thermal field inside vertical cavity, the “two-velocity-scale wall function” is chosen to be used with high-Re models. In both vertical and inclined cases the overall flow pattern appears similar, with a single circulation cell, and a boundary layer at the wall. The levels of turbulence energy are generally slightly lower in the inclined case. Most models give a reasonable prediction of measured Nusselt number, with the two low-Re approaches generally being closer to the data than the schemes employing wall functions. For the 15° inclined cavity, a multi cellular motion is shown by the high-Re models. Nevertheless, all the model predictions disagree with experimental data due to the presence in real flow of 3-D unsteady structures as found in Benard convection problems. These cannot, definitely, be reproduced using a 2-D geometry. Both highly resolved LES and unsteady RANS computations are then conducted, for turbulent natural convection of air inside 15° unstably and stably stratified cavities. In accordance with recent experimental data, the LES computations for both enclosures returned three-dimensional time-averaged flow fields. In the case of the unstably stratified enclosure, the flow is highly unsteady with coherent turbulent structures in the core of the enclosure. Results of LES computations show close agreement with the measured data. Subsequent comparisons of different URANS schemes with the present LES are used in order to explore to what extent these models are able to reproduce the large-scale unsteady flow structures. All URANS schemes have been found to be able to reproduce the 3-D unsteady flow features present in the 15° unstable cavity. However, the low-Re model tested as well as requiring a high resolution near-wall grid, also needed a finer grid in the core region than the high-Re models, thus making it computationally very expensive. Flow within the 15° stable cavity also shows some 3-D features, although it is significantly less unsteady, and the URANS models tested here have been less successful in reproducing this flow pattern. Finally, natural convection of CO2 inside a horizontal annular penetration enclosure, which can be found in AGR's, has been performed using a highly resolved LES and a set of RANS models. The Rayleigh number is 1.5×10e9. RANS models agree with the present LES on the fact that the flow is unsteady and there are large-scale oscillations present which decrease in amplitude as one moves from the open towards the closed end of the annular enclosure. Overall heat transfer and thermal quantitative and dynamic results show that RANS schemes are in close agreement with the current LES data except some discrepancies shown by the high-Re model which can be returned to the limitation of the simple wall function used to predict such complex flow.

Orientador: Amauri Garcia / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecânica / Made available in DSpace on 2018-08-18T19:08:04Z (GMT). No. of bitstreams: 1 Silva_MariaAdrinaPaixaodeSouzada_D.pdf: 5780858 bytes, checksum: b602bf47a710b2fb6f161a596020054f (MD5) Previous issue date: 2011 / Resumo: Ligas de alumínio dispersas com bismuto, chumbo e índio apresentam aplicações promissoras em componentes automotivos resistentes ao desgaste. Essas dispersões de elementos de baixa temperatura de fusão diminuem a dureza e escoam facilmente em condições de deslizamento, resultando em um comportamento tribológico favorável. Muitos estudos têm sido realizados a fim de melhor compreender as distintas morfologias obtidas pela reação monotética Algumas pesquisas assumem que a evolução do espaçamento interfásico (?) nas Iigas monotéticas ou obedece à clássica relação utilizada para eutéticos: ?2v = C, ou ainda a relação de crescimento dentrítico ?.va.Gb = C, onde v é a velocidade de solidificação e C uma constante para ambos os casos, G é o gradiente de temperatura e a e b são constantes. Tais estudos utilizaram fornos de aquecimento à resistência do tipo Bridgman para produzir a solidificação direcional de ligas monotéticas. Existe uma falta de estudos consistentes no desenvolvimento microestrutural de ligas monotéticas durante condições de fluxo de calor transitório, que são de importância primordial, uma vez que esse tipo de fluxo de calor engloba a maioria dos processes industriais e de solidificação. No presente estudo, foram realizados experimentos de solidificação unidirecional em regime não-estacionário com as Iigas hipomonotéticas Al-0,9%Pb e Al-5,5%In, com a liga monotética Al-1,2% Pb e com as Iigas hipermonotéticas Al-5,0%Bi, Al-7,0%Bi e Al- 2,1%Pb, além da análise da rnacrossegregação da liga monotética Al-3,2%Bi. Os parâmetros térmicos como velocidades de crescimento, taxas de resfriamento e gradientes térmicos foram determinados experimentalmente por curvas de resfriamento adquiridas ao longo do comprimento do lingote. A evolução microestrutural foi caracterizada por técnicas metalográficas e os espaçamentos foram correlacionados com os parâmetros térmicos de solidificação. Verificou-se que a lei de crescimento ?2v= C é válida apenas para os casos onde há ocorrência de partículas de soluto dispersas irregularmente na matriz e que esta lei não se aplicou nos casos onde foi observado uma frente monotética celular e morfologia de fibras e cordões de pérolas. Já a lei ?.va.Gb = C se mostrou aplicável em todos os casos ascendentes. O aumento do teor de soluto das ligas e o fluxo convectivo induzido provocaram mudanças nas morfologias das fases, na magnitude dos espaçamentos interfásicos e nos diâmetros das partículas / Abstract: Aluminum alloys dispersed with bismuth, lead and indium show promising applications in wear­resistant automotive components. Such dispersions of low melting temperature elements decrease hardness and flow easily under sliding conditions, resulting in favorable tribological behavior. Much research has been devoted in order to better comprehend the distinct morphologies obtained by monotectic reactions. Some researches assume that the interphase spacing evolution in monotectic alloys follows the classical relationship used for eutectics: ?2v = C, or the dendritic growth relationship ?.va.Gb = C, where v is the solidification velocity and C a constant value for both cases, G is the thermal gradient and a and b are constants. Such studies have used Bridgman type het resistance furnaces to produce the directionally solidified monotectic samples. There is a lack of consistent studies on the microstructural development of monotectic alloys during transient heat flowconditions, which are of prime importance since this class of heat flow encompasses the majority of solidification industrial processes. In the present study, directional unsteady-state solidification experiments were carried out with hypomonotectic Al-0.9wt%Pb and Al-5.5wt%In, monotectic Al-1.2wt%Pb and hypermonotectic Al-5.0wt%Bi, Al-7.0wt%Bi and Al-2.1wt%Pb alloys, besides the macrosegregation analysis of the monotectic Al-3.2wt%Bi alloy. Thermal parameters such as the growth rate, cooling rate and thermal gradient were experimentally determined by cooling curves recorded along the casting length. The microstructural evolution was characterized by metallography and the spacings were correlated with the thermal parameters. It is show that the ?2v = C growth law is valid only for the cases where there is a morphology of solute particles irregularly arranged in the matrix and it is not applied on the cases where a cellular monotectic front or fibers and strings of pearls morphologies were observed. On the other hand, ?.va.Gb = C law seems to be able to characterize all the upward cases. The increase in the alloy solute content and induced convective flow affected the morphologies of the resulting phases, the magnitude of the interphase spacing and the diameters of the solute-rich particles / Doutorado / Materiais e Processos de Fabricação / Doutor em Engenharia Mecânica

This work approaches the topic of modeling the atmospheric boundary layer in four research projects, which are summarized below. i) The diurnal cycles of near-surface meteorological parameters over Antarctic sea ice in six widely used atmospheric reanalyses were validated against observations from Ice Station Weddell. The station drifted from February through May 1992 and provided the most extensive set of meteorological observations ever collected in the Antarctic sea ice zone. For the radiative and turbulent surface fluxes, both the amplitude and shape of the diurnal cycles varied considerably among different reanalyses. Near-surface temperature, specific humidity, and wind speed in the reanalyses all featured small diurnal ranges, which, in most cases, fell within the uncertainties of the observed cycle. A skill score approach revealed the superiority of the ERA-Interim reanalysis in reproducing the observed diurnal cycles. An explanation for the shortcomings in the reanalyses is their failure to capture the diurnal cycle in cloud cover fraction, which leads to errors in other quantities as well. Apart from the diurnal cycles, NCEP-CFSR gave the best error statistics. ii) The accuracy of prediction of stable atmospheric boundary layers depends on the parameterization of the surface layer which is usually derived from the Monin-Obukhov similarity theory. In this study, several surface-layer models in the format of velocity and potential temperature Deacon numbers were compared to observations from CASES-99, Cardington, and Halley datasets. The comparisons were hindered by a large amount of scatter within and among datasets. Tests utilizing R2 demonstrated that the Quasi-Normal Scale Elimination (QNSE) theory exhibits the best overall performance. Further proof of this was provided by 1D simulations with the Weather Research and Forecasting (WRF) model. iii) The increasing number of physics parameterization schemes adopted in numerical weather forecasting models has resulted in a proliferation of inter-comparison studies in recent years. Many of these studies concentrated on determining which parameterization yields results closest to observations rather than analyzing the reasons underlying the differences. In this work, the performance of two 1.5-order boundary layer parameterizations was studied, the QNSE and Mellor-Yamada-Janjić (MYJ) schemes, in the Weather Research and Forecasting (WRF) model. The objectives were to isolate the effect of stability functions on the near-surface values and vertical profiles of virtual temperature, mixing ratio and wind speed. The results demonstrate that the QNSE stability functions yield better error statistics for 2-m virtual temperature but higher up the errors related to QNSE are slightly larger for virtual temperature and mixing ratio. A surprising finding is the sensitivity of the model results to the choice of the turbulent Prandtl number for neutral stratification (Prt0): in the Monin-Obukhov similarity function for heat, the choice of Prt0 is sometimes more important than the functional form of the similarity function itself. There is a stability-related dependence to this sensitivity: with increasing near-surface stability, the relative importance of the functional form increases. In near-neutral conditions, QNSE exhibits too strong vertical mixing attributed to the applied turbulent kinetic energy subroutine and the stability functions including the effect of Prt0. iv) In recent years, many eddy-diffusivity mass flux (EDMF) planetary boundary layer (PBL) parameterizations have been introduced. Yet, most validations are based on idealized setups and/or single column models. To address this gap, this study focused on the effect the mass flux part has on the performance in the QNSE-EDMF PBL scheme in the WRF model by comparing the results to observations from the CASES-97 field campaign. In addition, two refined versions, one introducing the parameterized clouds to the WRF radiation scheme, and the second adding a different entrainment formulation, were evaluated. The introduction of mass flux reduced errors in the average moisture profile but virtual temperature and wind speed profiles did not change as much. The turbulent flux profiles for modeled virtual potential temperature were little affected, with consistent reasonable agreement with observations, if one allows for biases in the observed data and modeled surface fluxes. However, the water vapor flux divergences from QNSE tend to be more negative than observed, while including the mass flux part tends to make the divergences more positive, the latter at least partially due to deeper model PBLs resulting from excessive model surface virtual temperature fluxes. Further, both virtual potential temperature and water vapor flux profiles display spurious spikes attributed to the way the non-local and local terms interact in the model. The influence of the mass flux schemes extends to 60 – 100-km scale circulation features, which were greatly modified by both the inclusion of mass flux and the new entrainment formulation. Adding mass flux based clouds to the radiation calculation improved the time and space averaged modeled incoming shortwave flux. The choice of the representation for entrainment/detrainment often affected the results to the same extent as adding mass flux did.

This thesis reviews the development and current methods of numerical convection heat transfer from available literature, encompassing an analysis of the various finite element formulations available for investigating convection. It further describes the finite element formulation for the primitive variable convection heat transfer equations via a Galerkin weighted residual scheme and using mixed interpolation, and it demonstrates the capability of this method by means of five practical examples, namely natural convection in a thermally driven square cavity, a thermally driven vertical slot, a thermally driven triangular cavity, and a liquid convective diode, and forced convection in a cooling pond. This study also provides the background and framework for the problem of transient convection heat transfer, and for further steady-state studies using parameters outside those considered herein.

This dissertation focuses on the effects of superhydrophobic (SHPo) surfaces on thermal transport. The work is divided into two main categories: thermal transport without phase change and thermal transport with phase change. Thermal transport without phase change is the topic of four stand-alone chapters. Three address jet impingement at SHPo surfaces and the fourth considers natural convection at a vertical, SHPo wall. Thermal transport with phase change is the topic of a single stand-alone chapter exploring pool boiling at SHPo surfaces.Two chapters examining jet impingement present analytical models for thermal transport; one considered an isothermal wall and the other considered an isoflux wall. The chapter considering the isothermal scenario has been archivally published. Conclusions are presented for both models. The models indicated that the Nusselt number decreased dramatically as the temperature jump length increased. Further, the influence of radial position, jet Reynolds number, Prandtl number and isoflux versus isothermal heating become negligible as temperature jump length increased. The final chapter concerning jet impingement reports an experimental exploration of jet impingement at post patterned SHPo surfaces with varying microfeature pitch and cavity fraction. The empirical results show a decrease in Nusselt number relative to smooth hydrophobic surfaces for small pitch and cavity fraction and the isoflux model agrees well with this data when the ratio of temperature jump length to slip length is 3.1. At larger pitch and cavity fractions, the empirical results have higher Nusselt numbers than the SHPo surfaces with small pitch and cavity fraction but remain smaller than the smooth hydrophobic surface. We attribute this to the influence of small wetting regions. The chapter addressing natural convection presents an analytical model for buoyant flow at a vertical SHPo surface. The Nusselt number decreased dramatically as temperature jump length increased, with greater decrease occurring near the lower edge and at higher Rayleigh number. Thermal transport with phase change is the topic of the final stand-alone chapter concerning pool boiling, which has been archivally published. Surface heat flux as a function of surface superheat was reported for SHPo surfaces with rib and post patterning at varying microfeature pitch, cavity fraction, and microfeature height. Nucleate boiling is more suppressed on post patterned surfaces than rib patterned surfaces. At rib patterned surfaces, transition superheat decreases as cavity fraction increases. Increasing microfeature height modestly increases the transition superheat. Once stable film boiling is achieved, changes in surface microstructure negligibly influence thermal transport.

Heat exchangers are used in a wide variety of industrial applications. Augmentation of heat transfer can realize a reduction in heat transfer size and increase the effectiveness and efficiency of heat exchangers. Heat transfer can be enhanced with various methods where the turbulence of the fluid flow is enhanced: by adding ribs, grooves or steps to the channel wall, using helical inserts, or by adding bluff bodies in the channel flow. By using these methods, there is also an increase in pressure drop penalty and larger pumping power is required to achieve the same flow rate. Circular cylindrical bluff bodies have been found to have smaller drag coefficients than square, rectangular or triangular cylindrical bluff bodies in the channel flow. Heat transfer and pressure drop experimental tests were done for eight different circular cylindrical cross-bar arrays at 15 different Reynolds numbers, in the range of 640 to 12 500. Eight different cross-bar configurations were tested: the cylinder diameter to pitch ratios were, d/p = 0.025, d/p = 0.05, d/pi=i0.1 and d/p = 0.2, and the angle to the flow direction, was θ = 90° and θ = 45° for each of the four different diameter-to-pitch ratios. Transient CFD simulations were done using Ansys fluent for d/p = 0.05 and d/p = 0.2, for θ = 90°, at Reynolds numbers 920 and 9 700, to analyze the secondary flow structures in the wake of the cylinders, partly responsible for the heat transfer and pressure drop increase in the channel flow in comparison to the smooth channel. The k-Ω shear stress transport (SST) model was used for the simulations. A mesh dependence study was done for spatial discretization, temporal discretization and validated against the experimental setup. The pressure drop gradient was found from the test data for the hydraulically developed part of the test section to calculate the friction factors. With an increase in Reynolds number, the friction factors decreased until reaching an asymptotic value for all the cross-bar configurations. For θi=i90° the friction factors were larger than for θ = 45° for the same d/p ratio and Reynolds number. With an increase in d/p, the friction factors increased. The largest measured friction factor was f = 0.3, for configuration d/p = 0.2, θ = 90°, at Re = 640 and the smallest measured friction factor f = 0.02, for d/pi= 0.025, θ = 45°, at Re = 12 500. The friction factor ratio was then used to quantify the pressure penalty for using cylindrical cross-bars in the channel flow to enhance heat transfer. The maximum friction factor ratio, f/f0 = 16.7 occurred at Re = 9 700, for d/pi=i0.2, θ = 90° and the minimum friction factor ratio, f/f0 = 2.1, at Re = 640, for d/pi=i0.025, θ = 45°. The average Nusselt numbers were then calculated using the spatial integral average of the local Nusselt numbers. With an increase in Reynolds number, there was an increase in the average Nusselt number for all the cylindrical cross-bar configurations. For larger d/p ratios and θ = 90° cases, the average Nusselt numbers were larger than for smaller d/p ratios and θ = 45°. The largest average Nusselt number was Nuavg = 66.3, at Re = 9 700 for d/p = 0.2, θ = 90° and the smallest average Nusselt number, Nuavg = 8.7, at Re = 640 for d/p = 0.025, θ = 45°. The Nusselt number ratio could then be used to quantify the heat transfer enhancement of the cylindrical cross-bar channel to that of the smooth channel, where the largest Nusselt number ratio was, Nuavg /Nu0,avg = 3.3, for d/p = 0.2, θ = 90°, at Rei=i3 000 and the smallest Nuavg /Nu0,avg = 1.1, for d/p = 0.025, θ = 45°, at Re = 640. The CFD results concluded that the pressure drop increase and heat transfer enhancement were caused by the flow acceleration, flow separation, eddy formation, vorticity increase, and boundary layer deformation next to and behind the cylinders. The Strouhal number for the larger d/p ratios suggested that the unsteadiness in the flow is higher for the cylinder arrays with a larger diameter, increasing both the heat transfer enhancement and friction factor in comparison with the smaller diameter cylinder arrays. Finally, the thermal performance coefficients could be calculated by using the friction factor ratios and Nusselt number ratios. The thermal performance coefficient combines the effects of the heat transfer and pressure penalty increase. The thermal performance coefficients increased from Re = 640 until Rei=i3 000 after which it decreased with an increase in Reynolds number. This is because the pressure penalty starts to outweigh the heat transfer increase caused by the turbulators. The largest thermal performance coefficient was η = 1.6, for d/p = 0.025, θ = 45°, at Re = 3 000, and the lowest, η = 0.79, for d/p = 0.05, θ = 90°, at Re = 640. / Dissertation (MEng (Mechanical Engineering))--University of Pretoria, 2021. / Mechanical and Aeronautical Engineering / MEng (Mechanical Engineering) / Unrestricted

This thesis focuses on design of device to generate hot air. This generator is supposed to be used to perform thermal tests for a company. Introduction summarizes general knowledge about ventilators, heat transport in air flow, and pressure losses. Main part of the thesis describes development of the generator. Initial design was followed by selection of suitable solution variants, design of major parts of the generator supported by calculations and overall conception of the generator. Final part summarizes the achievements and economic analysis.

Directional solidification has a wide interest due to its importance to the iron and steel industry. Examples of further application can be found in the aerospace industry regarding the manufacture of turbine blades and the semiconductor industry regarding single-crystal growth applications. Solute convection in the solidification process results in channel formation, which has a freckle-like appearance in cross-section and has a critical effect on the mechanical strength of a casting. For a solidification process that occurs via planar solidification from a solid boundary, one may consider the presence of three distinct regions often identified as horizontal layers, i.e. a fluid binary mixture (the melt), the solid layer and a two-phase (fluid-solid) mushy layer, separating the other two. The mushy layer is practically a porous medium consisting of an interconnected solid phase having its voids filled with the melt binary fluid. Channelling in the mushy layer and the creating of freckles are being considered the main reasons for non-homogeneous solidification and production of defects in the resulting solid product. The production of defects adversely affects the mechanical properties of the solid product leading to undesirable constraints on its industrial use. The purpose of this study is to evaluate the effect the vibrations have on the heat transfer during the solidification process as well as on the average density of the solid product and void formation. Experimental as well as theoretical investigations related to the solidification process were undertaken. Two effects that have been observed in previous experimental studies when metals and metal alloys are vibrated during solidification are a decrease in dendritic spacing, which directly affects density, and faster cooling rates and associated solidification times. Because these two effects happen simultaneously during solidification it is challenging to determine the one effect independently from the other. Most previous studies were on metals and metal alloys. In these studies, the one effect, i.e. the decrease in dendritic spacing, might influence the other, i.e. the faster cooling rates, and vice versa. The direct link between vibration and heat transfer has not yet been studied independently. The purpose of this study was to experimentally investigate the effect of vibration only on heat transfer and thus solidification rate. Experiments were conducted on paraffin wax, because it had a clearly defined macroscopic crystal structure consisting of mostly large straight-chain hydrocarbons. The advantage of the large straight-chain hydrocarbons was that the dendritic spacing was not affected by the cooling rate. Experiments were done with paraffin wax inside hollow plastic spheres of 40 mm diameter with 1 mm wall thickness. The paraffin wax was initially in a liquid state at a uniform temperature of 60°C and then submerged into a thermal bath at a uniform constant temperature of 15°C, which was approximately 20°C below the mean solidification temperature of the wax. Experiments were conducted in approximately 300 samples, with and without vibration at frequencies varying from 10 – 300 Hz. The first set of experiments were conducted to determine the solidification times. In the second set of experiments, the mass of wax solidified was determined at discrete time steps, with and without vibration. The results showed that paraffin wax had vibration independent of solid density contrary to other materials, eg. metals and metal alloys. Enhancement of heat transfer resulted in quicker solidification times and possible control over the heat transfer rate. The increase in heat transfer leading to faster solidifcation times was observed to first occur, as frequency increased and then to decrease. Experimental results showed that paraffin wax had vibration independent of solid density contrary to other materials, eg. metals and metal alloys. Enhancement of heat transfer resulted in quicker solidification times and possible control over the heat transfer rate. The increase in heat transfer leading to faster solidifcation times was observed to first occur, as frequency increased and then to decrease. Theoretical results of heat convection in a porous layer heated from below and subject to vibrations are presented by using a truncated spectral method in space. The partial differential equations governing the mass, momentum, heat, and solute transport were tranformed into a set of ordinary differential equations via a truncated modal expansion. Then the resutling equations were solved to identify the variety of regimes, and transitionbetween them, i.e. from steady convection, via periodic and quasi-periodic convection, towards chaotic or weak turbulent convection. The theoretcial results show that the heat convection subject to vibration is generally reduced when compared with the corresponding convection without vibrations. The exception for a certain frequency range shows about a 10% enhancement in the weak turbulent regime of convection, however, a 10% enhancement is still lower than the heat transfer prior to the transition to weak turbulence. Therefore, the heat transfer mechanism can be excluded as the main reason behind the improvement in solidification when vibrations are applied. Both experimental and theoretical results show an enhancement in heat transfer which correlate qualitativally. / Thesis (PhD)--University of Pretoria, 2014. / tm2015 / Mechanical and Aeronautical Engineering / PhD / Unrestricted

The thermophysical properties of nanofluids have attracted the attention of researchers to a far greater extent than the heat transfer characteristics of nanofluids have. Contradictory results on the thermal-fluid behaviour of nanofluids have been numerically and experimentally reported on in the open literature. Natural convection has not been investigated experimentally as much as the other properties of nanofluids. In this study, the characteristics and stability of Al2O3-water nanofluids (d = 20 30 nm) were analysed using a Malvern zetasizer, zeta potential and UV-visible spectroscopy. The natural convection of Al2O3- water nanofluids (formulated with a single-step method) was experimentally studied in detail for the volume fractions 0, 0.05, 0.1, 0.2, 0.4 and 0.6% in a rectangular cavity with an aspect ratio of 1, heated differentially on two opposite vertical walls for the Rayleigh number (Ra) range 3.49 x 10⁸ to 1.05 x 10⁹. The viscosity of Al2O3-water nanofluids measured between 15 and 50 °C. The effect of temperature and volume fraction on viscosity was also investigated. A detailed study of the nanoparticle concentration effect on the natural convection heat transfer coefficient was performed. It was found that increasing the concentration of nanoparticles improves the heat transfer coefficient by up to 15% at a 0.1% volume fraction. Further increasing the concentration of nanoparticles causes the natural convection heat transfer coefficient to deteriorate. This research also supports the idea that "for nanofluids with thermal conductivity more than the base fluids an optimum concentration may exist that maximises heat transfer in an exact condition as natural convection, laminar force convection or turbulence force convection". / Dissertation (MEng)--University of Pretoria, 2016. / Mechanical and Aeronautical Engineering / MEng / Unrestricted