O método de Lattice-Boltzmann (MLB) vem ganhando destaque nas últimas décadas pela sua capacidade de solução de escoamentos complexos como escoamentos multifásicos e multicomponentes, meios porosos e magnetohidrodinâmicos. Também existem extensões do método para a solução de problemas de advecção-difusão, que permitem acoplamento dessa metodologia para a solução numérica do conjunto de equação Navier-Stokes-Fourier. No trabalho presente, os objetivos iniciais podem ser resumidos no estudo dos conceitos básicos necessários para entender a derivação do método a partir da teoria cinética e estudo do método de Lattice Boltzmann, com operador de colisão mais simples conhecido como operador BGK, para problemas bidimensionais térmicos e fluidodinâmicos. A implementação numérica do método foi realizada em linguagem C, Matlab e CUDA C, com foco na solução de cinco problemas incompressíveis e laminares em regime permanente, conhecidos na literatura: equação viscosa de Burgers, escoamento de Poiseuille com transferência de calor, convecção natural em uma cavidade quadrada, regimes de convecção natural, forçada e mista em uma cavidade com uma das fronteiras móvel, e por fim, convecção forçada em uma cavidade ventilada, com uma entrada e uma saída. Além disso, um pequeno estudo no tempo computacional utilizando três implementações distintas foram testadas: implementação em série, uso de interpolações entre malhas grosseiras como condição inicial para malhas mais refinadas, e por fim, a adição da implementação do código em paralelo. Os ganhos de tempo entre a primeira e segunda estratégia ficaram entre 1.5 e 6, ao passo que o código paralelizado mostrou-se entre 20 e 25 vezes mais rápido que a segunda estratégia testada, comprovando o benefício de utilizar o processamento em paralelo em unidades gráficas. Os resultados obtidos para os problemas foram comparados com outros trabalhos da literatura, mostrando boa concordância para os primeiros quatro problemas estudados. Para a cavidade ventilada, diferenças relativas de até 15.7% no coeficiente de troca de calor e de até 28.38% para o coeficiente de queda de pressão adimensional foram observadas. Análises a respeito dos termos de erros do método apresentado, e a utilização de outras metodologias com método de Lattice Boltzmann, como por exemplo, o emprego de outros operadores de colisão, para trazerem mais estabilidade e precisão, podem elucidar melhor as divergências observados entre o trabalho presente e outros trabalhos da literatura. / Lattice-Boltzmann Method (LBM) has gained attention over last decades due to its performance in solving complex flows such as multiphase and multicomponent flow, porous media and magnetohydrodynamics. There are also further techniques which makes LBM able to solve advection-difussion problems, which allows coupling this methodology to solve Navier-Stokes-Fourier equations. In this work, initial goals are, in a few words, studying main concepts required to understand numerical method origins from kinectic thery, and studying the method implemetantion to solve 2D fluid dynamic and thermal problems. Numerical implemetation was written in C, CUDA C and Matlab languages, keeping in focus five different cases of laminar incompressible flows in steady state: viscous Burgers\' equation, thermal Poiseuille, natural convection in square cavities, natural, forced and mixed convection in a lid driven cavity, and, finally, forced convection in a ventilated cavity, with one inlet and outlet ports. Besides that, a quick study regarding impact in simulation time of three distinct approaches was conducted. First approach consists in implementing a serial code, the second in using interpolation for coarser grids results as initial condition for finer meshes, and the last one is adding parallelized code implementation. Time gains between first and second approach range from 1.5 to 6, while parallelized code was able to converge from 20 to 25 times faster than second approach, confirming great benefits in using graphics processing units. Results obtained from numerical solutions of problems were compared with other works from literature, and a good agreement among them was observed, specially in the first four problems studied. In ventilated cavity problem, relative differences up to 15.7% in heat transfer coefficient and up to 28.38% in dimensionless pressure drop coefficient were observed. A further error terms analysis of method used in this work, and use of distinct approaches, such as different collision operators, in order to increase numerical solution stability and accurarcy, may shed a light on results divergences observed between this work and others from literature.

O coletor solar de tubos evacuados une uma alta absortividade de radiação solar a um ótimo grau de isolamento térmico. Estas características, aliadas a um custo relativamente baixo, fazem deste tipo de coletor o mais utilizado no mundo. Por isso, diversos tipos de abordagens estão sendo adotadas para descrever seu comportamento térmico. Nesta linha, este trabalho foi desenvolvido através de abordagem experimental e teórica, sendo a última subdividida em numérica por CFD e analítica. A abordagem experimental contou com a construção de uma bancada para medida de temperatura e radiação, em um coletor de 24 tubos evacuados acoplados em um reservatório de 178 L. A abordagem por CFD utilizou um modelo tridimensional transiente. Com o modelo numérico validado, utilizando resultados experimentais, simularam-se diferentes condições de operação, em termos de ângulo de inclinação, fluxo de calor sobre o coletor, tamanho do reservatório e temperatura de entrada da água. Em relação à abordagem analítica, esta é dividida em: modelo de irradiância, modelo do tubo e modelo do reservatório. O modelo de irradiância determina a distribuição da radiação solar ao longo da circunferência do tubo. Parte desta radiação é absorvida pelo coletor e transferida para água. Esta quantidade é determinada com o modelo do tubo, que é baseado no método de resistências térmicas. O modelo do reservatório descreve o comportamento térmico da água em seu interior, tanto em aquecimento quanto em resfriamento, analisando a interação energética com o coletor e com o meio externo. O desenvolvimento do modelo do tubo passa pela avaliação da vazão mássica entre o tubo e o reservatório, além disso, o coeficiente de transferência de calor por convecção no interior do tubo deve ser determinado. Tais variáveis são determinadas a partir de uma correlação para o número de Reynolds, a qual foi obtida com resultados da avaliação por CFD e é função dos números de Nusselt, Grashof e Prandtl. A correlação proposta apresentou bom ajuste com os resultados numéricos. Com a bancada de ensaio experimental foram feitas medidas de temperatura da água no reservatório ao longo de alguns dias. Para as mesmas condições do experimento, a temperatura média da água no reservatório foi estimada com resultados da integração dos modelos de irradiância, do tubo e do reservatório. A diferença entre os resultados experimental e teórico foi de 4,8% para a energia acumulada. / The evacuated tube solar collector combines high solar radiation absorptivity to a great thermal insulation degree. These characteristics, combined with a relatively low cost, make this type of collector the most used in the world. Therefore, various types of approaches are being adopted to describe its thermal behavior. In this way, this work was developed through experimental and theoretical approaches, the latter being subdivided into numeric, by CFD, and analytical analysis. For the experimental approach a test bench was built. The tests was carried on a solar collector with 24 evacuated tubes coupled to a 178 L tank, measuring temperature and solar radiation. The CFD approach used a transient three-dimensional model. After the numerical model validation using experimental data, simulations was carried over different operating conditions in terms of angle, heat flux on the collector, tank size and water inlet temperature. The analytical approach is divided into: irradiance model, tube model and tank model. The irradiance model determines the irradiance distribution of solar radiation along the circumference of the tube. Part of this radiation is absorbed by the collector and transferred to water, this amount is determined with the tube model, using the thermal resistance method. The tank model describes the thermal behavior of inside water, both in heating and in cooling, analyzing energy interaction with the collector and the external environment. The development of the tube model involves the assessment of the mass flow rate between the tube and the tank, furthermore the convection heat transfer coefficient inside the tube must be determined. These variables are determined from a correlation for the Reynolds number, which was obtained with evaluation results by CFD. Proposed Reynolds number is a function of the Nusselt, Prandtl and Grashof numbers. The correlation presented a good agreement with the numerical results. Using the experimental test bench the water temperature was measured into the tank over a few days. For the same experiment conditions, the average temperature of the water into the tank was estimated by results of integration of irradiance, tube and tank models. The difference between the experimental and theoretical results was 4.8% for the accumulated energy.

Dans les écoulements de convection naturelle en canal, de nombreuses structures cohérentes se développent. Leur rôle dans les transferts de masse et de chaleurs reste aujourd’hui mal connu. Des études numériques basées sur des Simulations aux Grandes Echelles d’un canal vertical dont un des côtés est chauffé uniformément à flux constant ont donc été entreprises. Différentes étapes de la transition laminaire-turbulente sont identifiées et deux bandes de fréquences sont observées en proche paroi chauffée. Une méthode, basée sur la Proper Orthogonal Decomposition (POD), a aussi été utilisée et a permis une très nette séparation spatio-temporelle des modes les plus énergétiques de l’écoulement en fonction des fréquences précédemment identifiées. Il devient donc possible d’étudier la contribution de ces modes à l’écoulement. Une analyse POD des champs de vitesses expérimentaux permet l’identification de similarités avec les structures observées numériquement.Les variations du gradient de température extérieur au canal, ont été identifiées comme l’une des sources de ces différences et ont fait l’objet d’études expérimentales et numériques. Il est montré qu’un augmentation du gradient positif de température diminue le débit mais aussi déplace la transition plus bas dans le canal. Les études numériques permettent également d’étudier le cas de gradients de température faibles ou négatifs, gradients qui sont difficile à obtenir expérimentalement. Un modèle théorique de l’influence de la stratification extérieure sur le débit massique induit est aussi développé. L’accord entre les résultats prédits et les résultats expérimentaux et numériques est excellent. Ce dernier point souligne que la stratification thermique est un facteur clé dans ce type d’écoulement et que son influence est précisément quantifiée dans ce travail / In natural convective flows, complex coherent structures develop whose the role in heat and mass transfer are not well understood. A numerical study, based on Large-Eddy-Simulations of a vertical channel with one side uniformly was therefore carried out. Different stages of transitional flow development were identified numerically with two characteristic frequency bands being observed in the flow, near the heated wall. Methods derived from the Proper Orthogonal Decomposition (POD) was also used and allows the most energetic modes to be separated accordingly to two characteristic frequency bands found numerically. As result, the contribution of the two families of modes to the near wall turbulent heat transfer and velocity-temperature correlation has been evaluated. POD was also performed on experimental measurements showing similarities with the numerically observed structures.In this work, variations of the external thermal stratification have been identified as one possible source of these differences and its influence was therefore investigated experimentally and numerically. It is shown that the increase in the positive gradient of the external stratification not only decreases the mass flow rate but also displaces the transition height to a lower location in the channel. Numerical simulations also allow the study of cases of weak and negative thermal stratifications which are difficult to achieve in laboratories. A theoretical model of the influence of the external thermal stratification on the mass flow rate was also developed. There is an excellent agreement between the theoretical predictions and the experimentally and numerically obtained mass flow rates. This clearly highlights that external temperature distributions are key driving factors and their influence is accurately quantified in this work

Este trabalho propõe um modelo de transferência de calor e massa para prever a formação de gelo em um container preenchido com legumes congelados. O problema físico é modelado como um meio poroso composto pelo próprio produto e o ar em seu entorno. O regime de convecção natural é assumido dentro do container, o qual promove o transporte de massa. Como uma primeira validação, o modelo é simulado considerando diferentes temperaturas de ar externo, causadas por flutuações da vizinhança. Resultados para quatro ciclos de temperaturas foram comparados, variando separadamente a temperatura média do ar, amplitude e frequência de oscilação. De modo geral, é observado que a temperatura do produto se comporta assim como era esperado e este resultado é diretamente associado à formação de gelo dentro do container. A formação de gelo cresce com uma maior amplitude de oscilação, porém decresce com um aumento na frequência e na temperatura média. Os parâmetros do modelo foram obtidos para dois diferentes produtos: fatias de cenouras congeladas e vagens congeladas, ambos em meio ao ar. As definições de parâmetros são oriundas de revisão bibliográfica, medições experimentais e simulações numéricas. Os parâmetros encontrados para a caracterização desses meios porosos foram similares para ambos os produtos, mesmo eles possuindo diferentes geometrias. A validação experimental foi feita para as fatias de cenoura considerando dois ciclos de temperatura O modelo numérico é capaz de prever o campo de velocidades do ar, as temperaturas do produto e a formação de gelo local. Os resultados foram validados em relação a um grupo independente de resultados numéricos, tal comparação apresentou uma boa concordância. A circulação de ar encontrada é, de fato, devido à convecção natural. O comportamento da temperatura dos produtos simulados concorda com os valores medidos e os valores de temperaturas diferem por menos de 12%. Com respeito à formação de gelo, o modelo é capaz de prevê-la corretamente nas regiões mais suscetíveis a este fenômeno. Porém, a quantidade de gelo formado prevista pelo modelo (1,56 g/semana) é menor do que a experimental (4,67 g/semana), apesar de serem de mesma ordem de magnitude. O efeito de cada parâmetro no modelo é estudado visando detectar maneiras de aprimorar o modelo. Foi encontrado que os parâmetros mais importantes para a formação de gelo total são a difusividade de massa efetiva e o coeficiente de transferência de calor convectivo dentro do container. Ajustando estes parâmetros duas vezes foi possível encontrar resultados melhores com respeito à formação de gelo (3,09 g/semana). / A model of heat and mass transfer is proposed in order to predict frost formation into a closed container filled with frozen vegetables. The physical problem is modeled as a macroporous media composed by the product itself and the surrounding air. Natural convection air flow is assumed into the container, who promotes water mass transport. As a first validation, the model is simulated for several exterior air temperatures, under environmental fluctuations (boundary conditions). Results of four temperature cycles were compared, varying average air temperature, amplitude and frequency of oscillation, one by one. As a general result, it is observed that the product temperature behavior is as expected, and it is directly associated with frost formation into the container. Frost formation increases with large amplitude of oscillation, but decreases with higher frequencies and higher mean temperatures. Model parameters were obtained for two assembling: frozen slices of carrots and air, and frozen extra thin green beans and air. Parameter definition and evaluation combines literature review, measurements and numerical simulation. In general, parameters which characterize these porous media were similar for both products, even though they display different geometries. The experimental validation is performed for carrot slices with two temperature cycles The numerical model is able to predict air velocity field, air and product temperatures, and local frost formation. Results are validated in respect to a set of independent experimental results that shown a good agreement. Air flow circulation is as expected due to natural convection. Product temperature simulated behavior agrees with measurements, and temperature values differ by less than 12%. Respect to frost formation predictions, the model predicts correctly the most susceptible regions to frost formation. However, the quantity of frost formed predicted by the model (1.56 g/ week)is lower than the experimental one (4.67g/week), despite being of the same order of magnitude. The effect of each parameter in the model is study in order to detect how to improve the model. The most important parameters affecting total frost formation are effective mass diffusivity and convective heat coefficient into the storage container. Adjusting these parameters to twice, better results in terms of frost formation could be obtained (3.09 g/ week).

Cette thèse présente les résultats de deux études : la première concerne la convection naturelle turbulente dans une cavité rectangulaire chauffée uniformément par le bas et remplie d’un nanofluide et la seconde concerne l’investigation du transfert de chaleur conjugué dans un échangeur de chaleur à tubes ailetés.L’enceinte de la première étude a un faible rapport d’aspect. Ses parois gauche, droite et supérieure sont maintenues à une température relativement basse. Le fluide de travail est un nanofluide constitué d’eau et de nanoparticules, soient d’alumine (Al2O3), ou de cuivre (Cu) ou d’oxyde de cuivre (CuO). L’influence des paramètres tels que le nombre de Rayleigh (basé sur la hauteur H de la cavité et la densité de flux de chaleur), le type de nanofluide et la fraction volumique des nanoparticules sur la performance de refroidissement est présentée. Les équations de Navier-Stokes et les équations de conservation de la masse et de l'énergie sont résolues pour une géométrie bidimensionnelle par la méthode numérique des volumes finis. L'algorithme SIMPLE est utilisé pour le couplage pression-vitesse. La discrétisation des termes convectifs est faite avec le schéma QUICK. Le modèle de turbulence k-epsilon standard est utilisé. Le maillage du domaine simulé est généré par le code Gambit. Les résultats montrent que pour toutes les valeurs de Ra, le nombre de Nusselt moyen augmente d’une façon linéaire et monotone avec l’augmentation de la concentration des nanoparticules dans le fluide de base. Le flux de chaleur moyen prend des valeurs qui diminuent en fonction de l’ordre suivant : Cu, CuO et Al2O3.La deuxième étude est une investigation numérique des caractéristiques dynamique et thermique d'un échangeur de chaleur. Les calculs supposent un transfert de chaleur et un écoulement en régime permanent. Le nombre de Nusselt et le coefficient de frottement qui caractérisent l'échangeur de chaleur sont déterminés pour différentes valeurs du nombre de Reynolds. L’équation de conservation de l'énergie dans le fluide et l’équation de conduction de la chaleur dans le solide en trois dimensions ont été résolues avec les équations de la conservation de la masse et de la quantité de mouvement afin de déterminer ces caractéristiques. Les deux régimes d'écoulement laminaire et turbulent sont considérés. L'effet de la modélisation de la turbulence a été étudié en utilisant trois modèles différents (Spalart-Allmaras modèle de turbulence à une équation, le modèle k-epsilon ; standard et le modèle RSM). La validation du modèle a été effectuée en comparant les facteurs de frottement, f, et le facteur, j, de Colburn avec les données expérimentales trouvées dans la littérature. Les résultats tracés ont montré un bon accord qualitatif entre les résultats numériques et les données expérimentales. Les résultats montrent également que le plus simple des trois modèles de turbulence testés (à savoir. Spalart-Allmaras) donne les valeurs les plus proches des données expérimentales. / This thesis presents the results of two studies: the first concerns natural turbulent convection in a rectangular cavity heated from the bottom wall and filled with a nanofluid and the second relates to the investigation of conjugate heat transfer in a fin-and-tube heat exchanger.The cavity of the first study is tall and has a heat source embedded on its bottom wall, while its left, right and top walls are maintained at a relatively low temperature. The working fluid is a water based nanofluid having three nanoparticle types: alumina, copper and copper oxide. The influence of pertinent parameters such as the Rayleigh number, the type of nanofluid and solid volume fraction of nanoparticles on the cooling performance is studied. Steady forms of twodimensional Reynolds-Averaged-Navier-Stokes equations and conservation equations of mass and energy, coupled with the Boussinesq approximation, are solved by the volume control based on the discretisation method employing the SIMPLE algorithm for pressure-velocity coupling. Turbulence is modeled by using the standard k-epsilon model. The Rayleigh number, Ra, is varied from 2.49xE09 to 2.49xE11. The volume fractions of nanoparticles where varied too. Stream lines, isotherms, velocity profiles and temperature profiles are presented for various combinations of Ra, the type of nanofluid and solid volume fraction of nanoparticles. The results are reported in the form of average Nusselt number on the heated wall. It is shown that for all values of Ra, the average heat transfer rate from the heat source increases almost linearly and monotonically as the solid volume fraction increases. Finally the average heat transfer rate takes on values that decrease according to the ordering Cu, CuO and Al2O3.In the second study We determined the heat transfer and friction characteristics of a realistic fin-and-tube heat exchanger. The computations assume steady-state heat transfer and fluid flow. Nusselt number and friction factor characteristics of the heat exchanger are presented for various values of Reynolds numbers. The energy conservation and the heat conduction equations in 3 dimensions have been solved in the fluid and the solid respectivelyalong with the mass and momentum conservation equations in order to determine these characteristics. Both laminar and turbulent flow regimes are considered. The effect of turbulence modeling was investigated using three different models (the one equation Spalart-Allmaras turbulence model, the standard k-epsilon; model and the RSM model). The computations allowed the determination of the dynamic and thermal fields. Model validation was carried out by comparing the calculated friction factor f and Colburn j-factor to experimental results found in the literature. The plotted results showed a qualitatively good agreement between numerical results and experimental data. The results obtained also showed that the simplest of the three turbulence models tested(i.e. Spalart-Allmaras) gives the closest values to the experimental data.

Cette thèse présente une étude numérique d’un séchoir solaire indirect à convection naturelle destiné à sécher les produits agricoles (les figues). La première partie analyse un séchoir solaire contenant une chambre de séchage couplée à un absorbeur inversé et une cheminée solaire. Les simulations ont été faites afin de déterminer les champs dynamique et thermique sous l’influence de la variation de la configuration de la cheminée solaire et la taille de l’ouverture d’admission. Les équations de conservations basées sur le modèle de turbulence k-ε standard sont résolues par la méthode des volumes finis à l’aide du code commercial ANSYS-Fluent. La deuxième partie analyse l’effet d’ajouter un stockage thermique sous forme d’un lit en gravier dans la chambre de séchage. Le lit est modélisé comme un milieu poreux. En plus, cette étude propose l’utilisation d’une deuxième entrée d'air dans la chambre de séchage afin d’assurer une distribution thermique plus homogène au niveau des claies et de garantir un séchage plus uniforme. / This thesis presents a numerical study of an indirect natural convection solar dryer for drying agricultural products (Figs). The first part analyzes a solar dryer containing a drying chamber coupled to a reversed absorber and a solar chimney. Simulations were made to determine the dynamic and thermal fields under the influence of the variation of the solar chimney configuration and the size of the inlet opening. The governing equations based on the standard k-ε turbulence model are solved by the finite volume method using the ANSYS-Fluent commercial code. The second part analyzes the effect of adding a thermal storage in the form of a gravel bed in the drying chamber. The bed is modeled as a porous medium. Furthermore, this study proposes the use of a second air inlet in the drying chamber in order to ensure a more homogeneous thermal distribution at the level of the racks and to guarantee a more uniform drying.

This thesis is devoted to the investigation of passivesafety and inherent features of subcritical nucleartransmutation systems - accelerator-driven systems. The generalobjective of this research has been to improve the safetyperformance and avoid elevated coolant temperatures inworst-case scenarios like unprotected loss-of-ow accidents,loss-of-heat-sink accidents, and a combination of both theseaccident initiators. The specific topics covered are emergencydecay heat removal by reactor vessel auxiliary cooling systems,beam shut-off by a melt-rupture disc, safety aspects fromlocating heat-exchangers in the riser of a pool-type reactorsystem, and reduction of pressure resistance in the primarycircuit by employing bypass routes. The initial part of the research was focused on reactorvessel auxiliary cooling systems. It was shown that an 80 MWthPb/Bi-cooled accelerator-driven system of 8 m height and 6 mdiameter vessel can be well cooled in the case of loss-of-owaccidents in which the accelerator proton beam is not switchedoff. After a loss-of-heat-sink accident the proton beam has tobe interrupted within 40 minutes in order to avoid fast creepof the vessel. If a melt-rupture disc is included in the wallof the beam pipe, which breaks at 150 K above the normal coreoutlet temperature, the grace period until the beam has to beshut off is increased to 6 hours. For the same vessel geometry,but an operating power of 250 MWth the structural materials canstill avoid fast creep in case the proton beam is shut offimmediately. If beam shut-off is delayed, additional coolingmethods are needed to increase the heat removal. Investigationswere made on the filling of the gap between the guard and thereactor vessel with liquid metal coolant and using water spraycooling on the guard vessel surface. The second part of the thesis presents examinationsregarding an accelerator-driven system also cooled with Pb/Bibut with heat-exchangers located in the risers of the reactorvessel. For a pool type design, this approach has advantages inthe case of heat-exchanger tube failures, particularly if wateris used as the secondary uid. This is because a leakage ofwater from the secondary circuit into the Pb/Bi-cooled primarycircuit leads to upward sweeping of steam bubbles, which wouldcollect in the gas plenum. In the case of heatexchangers in thedowncomer steam bubbles may be dragged into the ADS core andadd reactivity. Bypass routes are employed to increase the owspeed in loss-of-ow events for this design. It is shown thatthe 200 MWth accelerator-driven system with heat-exchangers inthe riser copes reasonably well with both a loss-of-ow accidentwith the beam on and an unprotected loss-of-heat-sink accident.For a total-loss-of-power (station blackout) and an immediatebeam-stop the core outlet temperature peaks at 680 K. After acombined loss-of-ow and loss-of-heat-sink accident the beamshould be shut off within 4 minutes to avoid exceeding the ASMElevel D of 977 K, and within 8 minutes to avoid fast creep.Assuming the same core inlet temperature, both the reactordesign with heat-exchanger in the risers and the downcomershave similar temperature evolutions after a total-loss-ofpoweraccident. A large accelerator-driven system of 800 MWth with a 17 mtall vessel may eventually become a standard size. For thishigher power ADS, the location of the heat-exchangers hasgreater impact on the natural convection capability. This isdue to that larger heatexchangers have more inuence on thedistance between the thermal centers during a lossof- owaccident. The design with heat-exchangers in the downcomers,the long-term vessel temperature peaks at 996 K during aloss-of-ow accident with the beam on. This does not pose athreat of creep rupture for the vessel. However, the locationof the heat-exchangers in the downcomers will probably requiresecondary coolant other than water, like for example oil (fortemperatures not higher than 673 K) or Pb/Bi coolant.

In this thesis, we consider a projection-based stabilization method for solving buoyancy driven flows (natural convection problems). The method consists of adding global stabilization for all scales and then anti-diffusing these effects on the large scales defined by projections into appropriate function spaces. In this way, stabilization acts only on the small scales. We consider two different variations of buoyancy driven flows based on the projection-based stabilization. First, we focus on the steady-state natural convection problem of heat transport through combined solid and fluid media in a classical enclosure. We present the mathematical analysis of the projection-based method and prove existence, uniqueness and convergence of the approximate solutions of the velocity, temperature and pressure. We also present some numerical tests to support theoretical findings. Second, we consider a system of combined heat and mass transfer in a porous medium due to the natural convection. For the semi-discrete problem, a stability analysis of the projectionbased method and a priori error estimate are given for the Darcy-Brinkman equations in double-diffusive convection. Then we provide numerical assessments and a comparison with some benchmark data for the Darcy-Brinkman equations. In the last part of the thesis, we present a fully discrete scheme with the linear extrapolation of convecting velocity terms for the Darcy-Brinkman equations.

Finned heat sink use for electronics cooling via natural convection is numerically investigated. An experimental study from the literature that is for vertical surfaces is taken as the base case and the experimental setup is numerically modeled using commercial CFD software. The flow and temperature fields are resolved. A scale analysis is applied to produce an order-of-magnitude estimate for maximum convection heat transfer corresponding to the optimum fin spacing. By showing a good agreement of the results with the experimental data, the model is verified. Then the model is used for heat transfer from inclined surfaces. After a large number of simulations for various forward and backward angles between 0-90 degrees, the dependence of heat transfer to the angle and Rayleigh number is investigated. It is observed that the contributions of radiation and natural convection changes with the angle considerably. Results are also verified by comparing them with experimental results available in literature.

<p>This thesis is devoted to the investigation of passivesafety and inherent features of subcritical nucleartransmutation systems - accelerator-driven systems. The generalobjective of this research has been to improve the safetyperformance and avoid elevated coolant temperatures inworst-case scenarios like unprotected loss-of-ow accidents,loss-of-heat-sink accidents, and a combination of both theseaccident initiators. The specific topics covered are emergencydecay heat removal by reactor vessel auxiliary cooling systems,beam shut-off by a melt-rupture disc, safety aspects fromlocating heat-exchangers in the riser of a pool-type reactorsystem, and reduction of pressure resistance in the primarycircuit by employing bypass routes.</p><p>The initial part of the research was focused on reactorvessel auxiliary cooling systems. It was shown that an 80 MWthPb/Bi-cooled accelerator-driven system of 8 m height and 6 mdiameter vessel can be well cooled in the case of loss-of-owaccidents in which the accelerator proton beam is not switchedoff. After a loss-of-heat-sink accident the proton beam has tobe interrupted within 40 minutes in order to avoid fast creepof the vessel. If a melt-rupture disc is included in the wallof the beam pipe, which breaks at 150 K above the normal coreoutlet temperature, the grace period until the beam has to beshut off is increased to 6 hours. For the same vessel geometry,but an operating power of 250 MWth the structural materials canstill avoid fast creep in case the proton beam is shut offimmediately. If beam shut-off is delayed, additional coolingmethods are needed to increase the heat removal. Investigationswere made on the filling of the gap between the guard and thereactor vessel with liquid metal coolant and using water spraycooling on the guard vessel surface.</p><p>The second part of the thesis presents examinationsregarding an accelerator-driven system also cooled with Pb/Bibut with heat-exchangers located in the risers of the reactorvessel. For a pool type design, this approach has advantages inthe case of heat-exchanger tube failures, particularly if wateris used as the secondary uid. This is because a leakage ofwater from the secondary circuit into the Pb/Bi-cooled primarycircuit leads to upward sweeping of steam bubbles, which wouldcollect in the gas plenum. In the case of heatexchangers in thedowncomer steam bubbles may be dragged into the ADS core andadd reactivity. Bypass routes are employed to increase the owspeed in loss-of-ow events for this design. It is shown thatthe 200 MWth accelerator-driven system with heat-exchangers inthe riser copes reasonably well with both a loss-of-ow accidentwith the beam on and an unprotected loss-of-heat-sink accident.For a total-loss-of-power (station blackout) and an immediatebeam-stop the core outlet temperature peaks at 680 K. After acombined loss-of-ow and loss-of-heat-sink accident the beamshould be shut off within 4 minutes to avoid exceeding the ASMElevel D of 977 K, and within 8 minutes to avoid fast creep.Assuming the same core inlet temperature, both the reactordesign with heat-exchanger in the risers and the downcomershave similar temperature evolutions after a total-loss-ofpoweraccident.</p><p>A large accelerator-driven system of 800 MWth with a 17 mtall vessel may eventually become a standard size. For thishigher power ADS, the location of the heat-exchangers hasgreater impact on the natural convection capability. This isdue to that larger heatexchangers have more inuence on thedistance between the thermal centers during a lossof- owaccident. The design with heat-exchangers in the downcomers,the long-term vessel temperature peaks at 996 K during aloss-of-ow accident with the beam on. This does not pose athreat of creep rupture for the vessel. However, the locationof the heat-exchangers in the downcomers will probably requiresecondary coolant other than water, like for example oil (fortemperatures not higher than 673 K) or Pb/Bi coolant.</p>