Modelos algébricos para prever o comprimento de uma chama turbulenta têm sido foco de estudo de diversos grupos de pesquisa por suas aplicações na área de engenharia. O método experimental para obter o modelo varia desde visualizações simples, até técnicas fotográficas, este último com parâmetros fotográficos variando entre os autores. Técnicas fotográficas são usadas para estimar a altura de levantamento da base da chama, (Lift-Off) e o comprimento médio visível de chama (Visible Flame Length, VFL). Duas técnicas comuns que podem ser encontradas na literatura: por imagens de chama com baixo tempo de exposição e longo tempo de exposição, são comparados com um terceiro que se baseia na intensidade luminosa e na frequência de imagens de chama que ocupam um pixel. O melhor método foi utilizado para caracterizar o comportamento das chamas turbulentas de gás natural para diferentes regimes de velocidade do escoamento. Modelos algébricos que preveem o comprimento de chama, altura de levantamento e a velocidade crítica de extinção de chama são avaliados com os novos resultados experimentais. Logo após, os coeficientes numéricos dos melhores modelos algébricos são reajustados Finalmente, foram obtidos mapas de estabilidade relacionados à altura de levantamento e à velocidade crítica de extinção de chama para cada diâmetro em função da diluição com CO2 e do número adimensional de Reynolds. A terceira parte deste trabalho está focada no estudo da distribuição de radiação térmica. Em particular, foram consideradas três distâncias radiais medidas em comprimentos de chama (0,5 Lf, 1 Lf, 2 Lf) visando obter a distribuição do fluxo radiante experimental ao longo de um eixo vertical adjacente às chamas. Finalmente, os dados experimentais foram utilizados como dados de entrada em uma análise inversa com o objetivo de calcular os fatores de ponderação do modelo das múltiplas fontes ponderadas (por suas siglas em inglês WMPS). Nesta última parte, são apresentados frações radiantes e distribuições de fluxo de calor radiante de chamas de gás natural diluídas para diversas diluições com dióxido de carbono e diâmetros do queimador. / Predicting models for turbulent diffusion flame lengths have several applications driven the attention of many research groups. Since several studies use photographs to measure the flame length, with photographic parameters varying among authors, in other cases simple visualizations were used. It is important to explore possible discrepancies among measurement technics that could affect the results. Optical visualizations of turbulent diffusion flames are used to estimate the visible average flame length (VFL) and the lift-off. The study presents a study of three different methods to measure the VFL using optical techniques. The effect on the image of the main optic parameters such as focus, exposure time and ISO sensibility are analyzed. The VFL obtained with images in low exposure time and long exposure time are compared with a third optical method that is based on the luminous intensity and the frequency of flame images occupying a pixel. One method was used to characterize the behavior of turbulent diffusion flames of natural gas for a range of flames in function of the flow velocity. Universal non-dimensional models that describe the VFL, lift-off and the blow-out stability limit of gaseous jet diffusion flames in the still air have been compared with new experimental data. The numerical coefficients of the best models are adjusted. Finally, maps of stability related to lift-off and blow-out were obtained for each diameter in function of the dilution with CO2 and flow exit velocity expressed in non-dimensional Reynolds number The third part of this work focuses on the estimation of the thermal distribution of radiative flux from turbulent diffusion flames in laboratory-scale. The experimental measurements were gotten from the previous stability study. In particular, was considered three radial distances measured in flame lengths (0,5 Lf, 1 Lf, 2 Lf) aiming at obtaining the experimental radiant flux along a vertical axis adjacent to the flames. Finally, the experimental data was used as input data in an inverse analysis with the purpose of computing weight coefficients of the weighted multi-point source (WMPS) model. Then, experimental data that include: radiant fractions and radiative heat flux are presents for several flames with different dilutions with carbon dioxide and burner´s diameters.

Modelos algébricos para prever o comprimento de uma chama turbulenta têm sido foco de estudo de diversos grupos de pesquisa por suas aplicações na área de engenharia. O método experimental para obter o modelo varia desde visualizações simples, até técnicas fotográficas, este último com parâmetros fotográficos variando entre os autores. Técnicas fotográficas são usadas para estimar a altura de levantamento da base da chama, (Lift-Off) e o comprimento médio visível de chama (Visible Flame Length, VFL). Duas técnicas comuns que podem ser encontradas na literatura: por imagens de chama com baixo tempo de exposição e longo tempo de exposição, são comparados com um terceiro que se baseia na intensidade luminosa e na frequência de imagens de chama que ocupam um pixel. O melhor método foi utilizado para caracterizar o comportamento das chamas turbulentas de gás natural para diferentes regimes de velocidade do escoamento. Modelos algébricos que preveem o comprimento de chama, altura de levantamento e a velocidade crítica de extinção de chama são avaliados com os novos resultados experimentais. Logo após, os coeficientes numéricos dos melhores modelos algébricos são reajustados Finalmente, foram obtidos mapas de estabilidade relacionados à altura de levantamento e à velocidade crítica de extinção de chama para cada diâmetro em função da diluição com CO2 e do número adimensional de Reynolds. A terceira parte deste trabalho está focada no estudo da distribuição de radiação térmica. Em particular, foram consideradas três distâncias radiais medidas em comprimentos de chama (0,5 Lf, 1 Lf, 2 Lf) visando obter a distribuição do fluxo radiante experimental ao longo de um eixo vertical adjacente às chamas. Finalmente, os dados experimentais foram utilizados como dados de entrada em uma análise inversa com o objetivo de calcular os fatores de ponderação do modelo das múltiplas fontes ponderadas (por suas siglas em inglês WMPS). Nesta última parte, são apresentados frações radiantes e distribuições de fluxo de calor radiante de chamas de gás natural diluídas para diversas diluições com dióxido de carbono e diâmetros do queimador. / Predicting models for turbulent diffusion flame lengths have several applications driven the attention of many research groups. Since several studies use photographs to measure the flame length, with photographic parameters varying among authors, in other cases simple visualizations were used. It is important to explore possible discrepancies among measurement technics that could affect the results. Optical visualizations of turbulent diffusion flames are used to estimate the visible average flame length (VFL) and the lift-off. The study presents a study of three different methods to measure the VFL using optical techniques. The effect on the image of the main optic parameters such as focus, exposure time and ISO sensibility are analyzed. The VFL obtained with images in low exposure time and long exposure time are compared with a third optical method that is based on the luminous intensity and the frequency of flame images occupying a pixel. One method was used to characterize the behavior of turbulent diffusion flames of natural gas for a range of flames in function of the flow velocity. Universal non-dimensional models that describe the VFL, lift-off and the blow-out stability limit of gaseous jet diffusion flames in the still air have been compared with new experimental data. The numerical coefficients of the best models are adjusted. Finally, maps of stability related to lift-off and blow-out were obtained for each diameter in function of the dilution with CO2 and flow exit velocity expressed in non-dimensional Reynolds number The third part of this work focuses on the estimation of the thermal distribution of radiative flux from turbulent diffusion flames in laboratory-scale. The experimental measurements were gotten from the previous stability study. In particular, was considered three radial distances measured in flame lengths (0,5 Lf, 1 Lf, 2 Lf) aiming at obtaining the experimental radiant flux along a vertical axis adjacent to the flames. Finally, the experimental data was used as input data in an inverse analysis with the purpose of computing weight coefficients of the weighted multi-point source (WMPS) model. Then, experimental data that include: radiant fractions and radiative heat flux are presents for several flames with different dilutions with carbon dioxide and burner´s diameters.

Modelos algébricos para prever o comprimento de uma chama turbulenta têm sido foco de estudo de diversos grupos de pesquisa por suas aplicações na área de engenharia. O método experimental para obter o modelo varia desde visualizações simples, até técnicas fotográficas, este último com parâmetros fotográficos variando entre os autores. Técnicas fotográficas são usadas para estimar a altura de levantamento da base da chama, (Lift-Off) e o comprimento médio visível de chama (Visible Flame Length, VFL). Duas técnicas comuns que podem ser encontradas na literatura: por imagens de chama com baixo tempo de exposição e longo tempo de exposição, são comparados com um terceiro que se baseia na intensidade luminosa e na frequência de imagens de chama que ocupam um pixel. O melhor método foi utilizado para caracterizar o comportamento das chamas turbulentas de gás natural para diferentes regimes de velocidade do escoamento. Modelos algébricos que preveem o comprimento de chama, altura de levantamento e a velocidade crítica de extinção de chama são avaliados com os novos resultados experimentais. Logo após, os coeficientes numéricos dos melhores modelos algébricos são reajustados Finalmente, foram obtidos mapas de estabilidade relacionados à altura de levantamento e à velocidade crítica de extinção de chama para cada diâmetro em função da diluição com CO2 e do número adimensional de Reynolds. A terceira parte deste trabalho está focada no estudo da distribuição de radiação térmica. Em particular, foram consideradas três distâncias radiais medidas em comprimentos de chama (0,5 Lf, 1 Lf, 2 Lf) visando obter a distribuição do fluxo radiante experimental ao longo de um eixo vertical adjacente às chamas. Finalmente, os dados experimentais foram utilizados como dados de entrada em uma análise inversa com o objetivo de calcular os fatores de ponderação do modelo das múltiplas fontes ponderadas (por suas siglas em inglês WMPS). Nesta última parte, são apresentados frações radiantes e distribuições de fluxo de calor radiante de chamas de gás natural diluídas para diversas diluições com dióxido de carbono e diâmetros do queimador. / Predicting models for turbulent diffusion flame lengths have several applications driven the attention of many research groups. Since several studies use photographs to measure the flame length, with photographic parameters varying among authors, in other cases simple visualizations were used. It is important to explore possible discrepancies among measurement technics that could affect the results. Optical visualizations of turbulent diffusion flames are used to estimate the visible average flame length (VFL) and the lift-off. The study presents a study of three different methods to measure the VFL using optical techniques. The effect on the image of the main optic parameters such as focus, exposure time and ISO sensibility are analyzed. The VFL obtained with images in low exposure time and long exposure time are compared with a third optical method that is based on the luminous intensity and the frequency of flame images occupying a pixel. One method was used to characterize the behavior of turbulent diffusion flames of natural gas for a range of flames in function of the flow velocity. Universal non-dimensional models that describe the VFL, lift-off and the blow-out stability limit of gaseous jet diffusion flames in the still air have been compared with new experimental data. The numerical coefficients of the best models are adjusted. Finally, maps of stability related to lift-off and blow-out were obtained for each diameter in function of the dilution with CO2 and flow exit velocity expressed in non-dimensional Reynolds number The third part of this work focuses on the estimation of the thermal distribution of radiative flux from turbulent diffusion flames in laboratory-scale. The experimental measurements were gotten from the previous stability study. In particular, was considered three radial distances measured in flame lengths (0,5 Lf, 1 Lf, 2 Lf) aiming at obtaining the experimental radiant flux along a vertical axis adjacent to the flames. Finally, the experimental data was used as input data in an inverse analysis with the purpose of computing weight coefficients of the weighted multi-point source (WMPS) model. Then, experimental data that include: radiant fractions and radiative heat flux are presents for several flames with different dilutions with carbon dioxide and burner´s diameters.

Ce travail de thèse est dédié à la modélisation du thermoformage du verre. Le procédé consiste à déformer une plaque de verre sous l'effet de son propre poids. Posée sur un support et placée dans un four, la température de la pièce augmente et sa viscosité diminue, ce qui permet d'obtenir la forme désirée. Les simulations numériques, qui se basent sur un modèle thermomécanique, doivent permettre de mieux comprendre l'influence, sur le produit final, des différents paramètres d'essai, comme le chargement thermique, la géométrie et le matériau du moule ou encore la forme initiale de la pièce. Pour ce faire, le logiciel commercial Abaqus®, qui utilise une méthode de résolution des calculs par éléments finis, prend en charge les aspects mécaniques et conductifs. En revanche, comme le verre est un milieu semi-transparent, la modélisation du transfert radiatif est complexe et nécessite le développement d'un code se basant sur une méthode de Monte Carlo dite réciproque. La méthode a été validée en deux dimensions sur des cas-tests de la littérature scientifique. Le code a ensuite été implémenté dans le logiciel Abaqus® afin de réaliser des simulations de thermoformage sur moule et en suspension. Le verre est considéré comme un matériau élasto-visco-plastique obéissant à un modèle de Maxwell simple et la thermodépendance de la viscosité est prise en compte par une loi WLF. Une attention particulière a été accordée au modèle radiatif. Différentes hypothèses, issues de la littérature scientifique, sont testées afin de vérifier leur validité dans notre cas d'étude / This dissertation is dedicated to the modelling of the glass sagging process. This operation consists in forming a sheet or a plate of glass by heating it in a furnace. Glass temperature rises and reaches a work temperature at which viscosity is low enough to allow glass to sag under its own weight due to gravity. Numerical simulation, based on a thermomechanical model, can help to better understand the influence of the different parameters on the final product, such as the thermal loading, the shape and material of the mould or even the initial geometry of the glass plate. Thus, the commercial software Abaqus® is used to solve the problem by a finite elements method. However, it cannot render the complexity of the radiative heat transfer in glass. So, a Monte Carlo code based on a reciprocal method was developed and validated using benchmarks from the scientific literature. Then, the code was implemented into Abaqus® in order to simulate glass sagging on a mould or glass forming by the draping process. Glass is considered as an elasto-viscoplastic material which obeys a Maxwell model. Viscosity is dependant to temperature according to a WLF law. Special care was devoted to the radiative heat transfer. Different hypothesis are reviewed and performed to check their validity when applied to our numerical set-up

Master’s thesis deals with numerical finding the global minimum. A theoretical part of project presents a problem of global optimization, her basic division and shows basic examples the specific problems of global optimization. Further, the reader is familiar with the basics of some numerical methods and with description of their strengths and weaknesses. A practical part shows concrete algorithm of global optimization, its use on data of specified function. Data refer to the issue of radiative heat transfer.

A three-dimensional computational tool was developed that simulates the heat and mass transfer interaction in a soil-root-stem system (SRSS) for a tree in a seasonally varying deciduous forest. The development of the SRSS model involved the modification and coupling of existing heat and mass transport tools to reproduce the three-dimensional diurnal internal and external temperatures, internal fluid distribution, and heat flow in the soil, roots, and stems. The model also required the development of a parallel Monte-Carlo algorithm to simulate the solar and environmental radiation regime consisting of sky and forest radiative effects surrounding the tree. The SRSS was tested, component-wise verified, and quantitatively compared with published observations. The SRSS was applied to simulate a tree in a dense temperate hardwood forest that included the calculations of surface heat flux and comparisons between cases with fluid flow transport and periods of zero flow. Results from the winter simulations indicate that the primary influence of temperature in the trunk is solar radiation and radiative energy from the soil and surrounding trees. Results from the summer simulation differed with previous results, indicating that sap flow in the trunk altered the internal temperature change with secondary effects attributed to the radiative energy from the soil and surrounding trees. Summer simulation results also showed that with sap flow, as the soil around the roots become unsaturated, the flow path for the roots will be changed to areas where the soil is still saturated with a corresponding increase in fluid velocity.

Multi-core, multi-processor, and Graphics Processing Unit (GPU) computer architectures pose significant challenges with respect to the efficient exploitation of parallelism for large-scale, scientific computing simulations. For example, a simulation of the human tonsil at the cellular level involves the computation of the motion and interaction of millions of cells over extended periods of time. Also, the simulation of Radiative Heat Transfer (RHT) effects by the Photon Monte Carlo (PMC) method is an extremely computationally demanding problem. The PMC method is example of the Monte Carlo simulation method—an approach extensively used in wide of application areas. Although the basic algorithmic framework of these Monte Carlo methods is simple, they can be extremely computationally intensive. Therefore, an efficient parallel realization of these simulations depends on a careful analysis of the nature these problems and the development of an appropriate software framework. The overarching goal of this dissertation is develop and understand what the appropriate parallel programming model should be to exploit these disparate architectures, both from the metric of efficiency, as well as from a software engineering perspective. In this dissertation we examine these issues through a performance study of PathSim2, a software framework for the simulation of large-scale biological systems, using two different parallel architectures’ distributed and shared memory. First, a message-passing implementation of a multiple germinal center simulation by PathSim2 is developed and analyzed for distributed memory architectures. Second, a germinal center simulation is implemented on shared memory architecture with two parallelization strategies based on Pthreads and OpenMP. Finally, we present work targeting a complete hybrid, parallel computing architecture. With this work we develop and analyze a software framework for generic Monte Carlo simulations implemented on multiple, distributed memory nodes consisting of a multi-core architecture with attached GPUs. This simulation framework is divided into two asynchronous parts: (a) a threaded, GPU-accelerated pseudo-random number generator (or producer), and (b) a multi-threaded Monte Carlo application (or consumer). The advantage of this approach is that this software framework can be directly used within any Monte Carlo application code, without requiring application-specific programming of the GPU. We examine this approach through a performance study of the simulation of RHT effects by the PMC method on a hybrid computing architecture. We present a theoretical analysis of our proposed approach, discuss methods to optimize performance based on this analysis, and compare this analysis to experimental results obtained from simulations run on two different hybrid, parallel computing architectures. / Ph. D.

L'objectif de ce travail est d’évaluer le positionnement et le potentiel des récepteurs à lit fluidisé à changement de section par rapport aux autres méthodes de chauffage de gaz à haute température par voie solaire. A cette fin, une connaissance approfondie des phénomènes thermiques et hydrodynamiques du récepteur est nécessaire. Pour acquérir cette connaissance, nous avons modélisé les transferts thermiques dans le récepteur en portant une attention particulière sur les transferts radiatifs en prenant en compte les diffusions multiples de la lumière dans le milieu particulaire, les effets de parois sur les transferts radiatifs et la directionnalité du rayonnement solaire concentré. La détermination précise de la distribution de particules dans le ciel du lit fluidisé s'est avérée un paramètre critique pour le calcul des transferts thermiques. Ces modèles, plus tard affinés par une confrontation avec des références expérimentales, nous ont permis d'explorer l'effet de la géométrie sur les transferts thermiques dans le récepteur. Il ont permis entre autres de mettre en évidence l'intérêt d'utiliser une colonne de fluidisation à changement de section et l'importance de l'optimisation du couple concentrateur solaire / récepteur afin d'éviter d'éventuelles surchauffes au niveau des parois du récepteur. De même, il semble que l'homogénéisation de la température dans le lit fluidisé contenu dans le récepteur soit favorable à son rendement. / The aim of this work is to evaluate the position and the potential of solar fluidized bed receivers compared to other methods for the solar heating of gases at high temperature. To this end, a thorough knowledge of the heat transfer and hydrodynamic of the receiver is necessary. To acquire this knowledge, we modeled the heat transfer in the receiver with a focus on the radiative transfer by taking into account the multiple scattering of light in the particle medium, the effect of walls on radiative heat transfer and the directionality of the concentrated solar radiation. The accurate determination of the distribution of particles within the fluidized bed has been a critical parameter for the calculation of heat transfer. With these models, later refined by a confrontation with experimental references, we have studied the effect of geometry on heat transfer in the receiver. This study highlighted the necessity to use a switching section fluidization column and the importance to optimize the pair : solar concentrator / receiver to avoid any overheating at the walls of the receiver. Moreover, it appears that the homogenization of the temperature in the fluidized bed of the receiver increase its performance.

Recently, the supercritical carbon dioxide (s-CO2) Brayton cycle has been identified as a promising candidate for solar-thermal energy conversion due to its potentially high thermal efficiency (50%, for turbine inlet temperatures of ~ 1000 K). Realization of such a system requires development of solar receivers which can raise the temperature of s-CO2 by over 200 K, to a receiver outlet temperature of 1000 K. Volumetric receivers are an attractive alternative to tubular receivers due to their geometry, functionality and reduced thermal losses. A concept of a ceramic pressurized volumetric receiver for s-CO2 has been developed in this work. Computational Fluid Dynamics (CFD) analysis along with a Discrete Ordinate method (DOM) radiation heat transfer model has been carried out, and the results for temperature distribution in the receiver and the resulting thermal efficiency are presented. Issues regarding material selection for the absorber structure, window, coating, receiver body and insulation are also addressed. A modular small scale prototype with 0.5 kWth solar heat input has been designed. The design of a small scale s-CO2 loop for testing this receiver module is also presented in this work. There is a lot of ongoing investigation for design and simulation of different configurations of heat exchangers and solar receivers using s-CO2 as the working fluid, in which wall temperatures up to 1000 K are encountered. While CO2 is considered to be transparent as far as solar radiation spectrum is concerned, there may be considerable absorption of radiation in the longer wavelength range associated with radiation emission from the heated cavity walls and tubes inside the receivers. An attempt has been made, in this study, to include radiation modelling to capture the effect of absorption bands of s-CO2 and the radiative heat transfer among the equipment surfaces. As a case study, a numerical study has been performed to evaluate the contribution of radiative heat transfer as compared to convection and conduction, for s-CO2 flow through a circular pipe. The intent is to provide a guideline for future research to determine the conditions for which radiation heat transfer modelling inside the pipe can be significant, and what errors can be expected otherwise. The effect of parameters such as Reynolds number, pipe diameter, length to diameter ratio, wall emissivity and total wall heat flux has been studied. The effect of radiation modelling on wall temperatures attained for certain amount of heat flux to be transferred to s-CO2 is also studied. The resulting temperature distribution, in turn, affects the estimation of heat loss to the environment

The present thesis is concerned with methods for determination and modeling of characteristic parameters of combustion of gaseous fuels. The focus is stressed on formation of nitrogen oxides and heat transfer from hot flue gases into combustion chamber’s walls. Experimental work, which is focused on testing of two burners with suppressed formation of nitrogen oxides, is an important part of the thesis. Its aim is to obtain data that is necessary for further processing and modeling. The work presents two methods that may be used in modeling of characteristic combustion parameters, namely the method based on statistical processing of data and the method based on computational fluid dynamics. The approaches are applied to two devices (burner with two-staged fuel supply, burner with two-staged air supply) with the objective to analyze their parameters. First approach covers detailed planning of burner test prior to its own carrying out (definition of the goal of experiment, choice of input factors and response, experimental plan) and subsequent statistical processing of experimental data. On the contrary, CFD approach offers simulations as an alternative option to traditional experimental methods. The simulation of combustion includes building of computational grid, setup of boundary conditions, turbulence model, heat transfer model and chemical kinetics. Results of simulations are compared with experimental measured data.