No description available.

A computer program based on the discrete-element method has been developed to compute friction factors and Nusselt Numbers for fully-developed turbulent flows with randomly-rough surfaces. Formulations of the discrete-element model for fully-developed turbulent flows inside circular pipes and between infinite parallel plates with the necessary adaptations for randomly-rough surfaces are provided. Utilizing the output of a three-dimensional profilometer, proper description of the randomly-rough surface is necessary for use within the discrete-element model. Proper description of the randomly-rough surface is achieved by the McClain (2002) method of characterization. Predictions from the discrete-element model computer program are compared with the classical, laminar and turbulent, smooth-wall results. In addition to the smooth-wall evaluations, predictions are compared with experimental results for turbulent internal flows with deterministic surface roughness. Predictions from the model demonstrated excellent agreement in all cases. Friction factor and Nusselt Number predictions for fully-developed flows over randomly-rough surfaces are also presented. With the friction factor and Nusselt Number data, velocity profiles for flows over randomly-rough, deterministically-rough and smooth surfaces are provided for comparison.

Gas-solid flow is fundamental to many industrial processes. Extensive experimental and numerical studies have been devoted to understand the interphase momentum and heat transfer in these systems. Most of the studies have focused on spherical particle shapes, however, in most natural and industrial processes, the particle shape is seldom spherical. In fact, particle shape is one of the important parameters that can have a significant impact on momentum, heat and mass transfer, which are fundamental to all processes. In this study particle-resolved simulations are performed to study momentum and heat transfer in flow through a fixed random assembly of ellipsoidal particles with sphericity of 0.887. The incompressible Navier-Stokes equations are solved using the Immersed Boundary Method (IBM). A Framework for generating particle assembly is developed using physics engine PhysX. High-order boundary conditions are developed for immersed boundary method to resolve the heat transfer in the vicinity of fluid/particle boundary with better accuracy. A complete framework using particle-resolved simulation study assembly of particles with any shape is developed. The drag force of spherical particles and ellipsoid particles are investigated. Available correlations are evaluated based on simulation results and recommendations are made regarding the best combinations. The heat transfer in assembly of ellipsoidal particle is investigated, and a correlation is proposed for the particle shape studied. The lift force, lateral force and torque of ellipsoid particles in assembly and their variations are quantitatively presented and it is shown that under certain conditions these forces and torques cannot be neglected as is done in the larger literature. / Ph. D.

RANS models in CFD are used to predict the liner wall heat transfer characteristics of a gas turbine annular combustor with radial swirlers, over a Reynolds number range from 50,000 to 840,000. A three dimensional hybrid mesh of around twenty five million cells is created for a periodic section of an annular combustor with a single radial swirler. Different turbulence models are tested and it is found that the RNG k-e model with swirl correction gives the best comparisons with experiments. The Swirl number is shown to be an important factor in the behavior of the resulting flow field. The swirl flow entering the combustor expands and impinges on the combustor walls, resulting in a peak in heat transfer coefficient. The peak Nusselt number is found to be quite insensitive to the Reynolds number only increasing from 1850 at Re=50,000 to 2200 at Re=840,000, indicating a strong dependence on the Swirl number which remains constant at 0.8 on entry to the combustor. Thus the peak augmentation ratio calculated with respect to a turbulent pipe flow decreases with Reynolds number. As the Reynolds number increases from 50,000 to 840,000, not only does the peak augmentation ratio decrease but it also diffuses out, such that at Re=840,000, the augmentation profiles at the combustor walls are quite uniform once the swirl flow impinges on the walls. It is surmised with some evidence that as the Reynolds number increases, a high tangential velocity persists in the vicinity of the combustor walls downstream of impingement, maintaining a near constant value of the heat transfer coefficient. The computed and experimental heat transfer augmentation ratios at low Reynolds numbers are within 30-40% of each other. / Master of Science

We examine the Rayleigh–Bénard convection as modelled by the Boussinesq equation. Our aim is at deriving bounds for the heat enhancement factor in the vertical direction, the Nusselt number, which reproduce physical scalings. In the first part of the dissertation, we examine the the simpler model when the acceleration of the fluid is neglected (Pr=∞) and prove the non-optimality of the temperature background field method by showing a lower bound for the Nusselt number associated to it. In the second part we consider the full model (Pr<∞) and we prove a new upper bound which improve the existing ones (for large Pr numbers) and catches a transition at Pr~Ra^(1/3).

Le présent travail apporte une contribution à la compréhension des mécanismes des transferts combinés de chaleur et de masse en convection naturelle lors de l’évaporation d’un film liquide d’eau d’épaisseur négligeable dans un canal vertical ondulé. L’écoulement est laminaire et bidimensionnel. Les équations régissant le phénomène sont résolues à l’aide d’une méthode aux volumes finis et le traitement du couplage vitesse-pression est réalisé par la méthode de projection. Les influences de la densité de flux de chaleur, de la température ainsi que l’humidité de l’air à l’entrée et la forme de la paroi du canal sur les transferts sont étudiées. Les résultats sont présentés sous la forme de ligne de courant, d’isothermes et d’iso-concentrations.Les simulations numériques effectuées ont permis l’étude détaillée de la structure de l’écoulement ainsi que des champs thermiques et massiques. Nous représentons également, les nombres de Nusselt et de Sherwood. / The present work is a contribution to the understanding of the mechanisms of combined heat and mass transfers in natural convection during the evaporation of a liquid film with negligible thickness in a wavy vertical channel. The flow is laminar and two-dimensional. The equations governing the phenomenon are resolved using the finite volumes method and the treatment of the coupling between velocity and pressure is carried out by the projection method. The influences of the heat flux density, the temperature and the humidity of the inlet air and the shape of the channel wall on the transfers are studied. The results are presented in the form of cstreamlines, isotherms and iso-concentrations.The numerical simulations carried out have allowed the detailed study of the flow structure as well as the thermal and mass fields. We also represent the Nusselt and Sherwood numbers.

In this thesis, we consider two models from physics, which are characterized by the interplay of thermodynamical and fluid mechanical phenomena: demixing (spinodal decomposition) and Rayleigh--Bénard convection. In both models, we investigate the dependencies of certain intrinsic quantities on the system parameters. The first model describes a thermodynamically driven demixing process of a binary viscous fluid. During the evolution, the two components of the mixture separate into two domains of the different equilibrium volume fractions. One observes a clear tendency: Larger domains grow at the expense of smaller ones, and thus, the average domain sizes increases --- a phenomenon called coarsening. It turns out that two mechanisms are relevant for the coarsening process. At an early stage of the evolution, material transport is essentially mediated by diffusion; at a later stage, when the typical domain size exceeds a certain value, due to the viscosity of the mixture, a fluid flow sets in and becomes the relevant transport mechanism. In both regimes, the growth rates of the typical domain size obey certain power laws. In this thesis, we rigorously establish one-sided bounds on these growth rates via a priori estimates. The second model, Rayleigh--Bénard convection, describes the behavior of a fluid between two rigid horizontal plates that is heated from below and cooled from above. There are two competing heat transfer mechanisms in the system: On the one hand, thermodynamics favors a state in which temperature variations are locally minimized. Thus, in our model, the thermodynamical equilibrium state is realized by a temperature with a linearly decreasing profile, corresponding to pure conduction. On the other hand, due to differences in the densities of hot and cold fluid parcels, buoyancy forces act on the fluid. This results in an upward motion of hot parcels and a downward motion of cold parcels. We study the dependence of the average upward heat flux, measured in the so-called Nusselt number, on the temperature forcing encoded by the container height. It turns out that the efficiency of the heat transport is independent of the height of the container, and thus, the Nusselt number is a constant function of height. Using a priori estimates, we prove an upper bound on the Nusselt number that displays this dependency --- up to logarithmic errors. Further investigations on the flow pattern in Rayleigh--Bénard convection show a clear separation of length scales: Along the horizontal top and bottom plates one observes thin boundary layers in which heat is essentially conducted, whereas the large bulk is characterized by a convective heat flow. We give first rigorous results in favor of linear temperature profiles in the boundary layers, which indicate that heat is indeed essentially conducted close to the boundaries.:1 Introduction 2 Coarsening rates in binary viscous fluids 2.1 Background from physics 2.2 Background from mathematics 2.3 The model 2.4 The gradient flow structure 2.5 Heuristics 2.6 Numerical simulations 2.7 Main results 2.8 Preliminaries 2.9 Proof of upper bounds on coarsening rates 2.10 Appendix: Well-posedness and regularity of solutions 3 Scaling of the Nusselt number 3.1 Background from physics 3.2 The model and the Nusselt number 3.3 Heuristics 3.4 Main results 3.5 Scaling law in the linear regime 3.6 Preliminaries and review 3.7 Upper bound using the background field method 3.8 Upper bound using the maximum principle 3.9 Appendix: Some elementary estimates 4 The laminar boundary layer 4.1 Background, model, and motivation 4.2 Main results 4.3 Preparation: Bounds on the velocity field 4.4 On the energy distribution 4.5 Bounds on the second order derivatives of the temperature field 4.6 Bounds on the third order derivatives of the temperature field

A técnica de transformada integral generalizada é utilizada para a obtenção de soluções analíticas do problema de convecção forçada em regime laminar, na região de entrada de dutos de geometria simples, tais como canais de placas paralelas e dutos circulares. Nessa região as camadas limites hidrodinâmica e térmica estão se desenvolvendo simultaneamente. O fluido é considerado como newtoniano e suas propriedades físicas como constantes. São utilizadas distribuições de velocidades longitudinais na forma analítica, as quais estão disponíveis na literatura e foram obtidas por métodos de linearização da equação de qualidade de movimento na direção axial. A partir da distribuição de velocidade axial obtém-se a distribuição de velocidade normal e analisa-se a influência desta velocidade sobre os parâmetros de transferência de calor: temperatura média de mistura e números de Nusselt local e médio. Analisa-se também a influência de dois perfis de velocidades diferentes sobre os resultados. Obtém-se resultados numéricos para a temperatura média de mistura e números de Nusselt local e médio considerando-se a geometria de um canal de placas planas paralelas, resolvendo-se um sistema completo de equações diferenciais ordinárias acopladas, para vários números de Prandtl. Implementam-se também, soluções aproximadas para um cálculo mais rápido dos resultados, verificando-se a precisão de tais soluções. Quando possível os resultados são comparados com resultados existentes na literatura.

Resumo: O objetivo deste trabalho é a análise numérica e experimental de escoamento viscoso, incompressível, permanente, com transferência de calor, em um canal estreito contendo um arranjo de obstáculos retangulares. A análise experimental envolveu determinação de coeficiente de transferência de calor médio bem como o número de Nusselt médio e medidas de temperatura em esteira térmica para comparação com os resultados obtidos por simulação numérica. Para a análise numérica usamos o programa comercial de mecânica dos fluidos e transferência de calor computacional ICEPAK®. Verificamos que quanto mais adentro o obstáculo estiver no arranjo maior é a transferência de calor por convecção forçada. Determinamos coeficientes de transferência de calor médio e número de Nusselt médio (com incerteza entre 6 e 15%) e verificamos que o efeito da posição diminui à medida que a velocidade aumenta. Concluímos também que ambos os modelos de turbulência utilizados, k-ε padrão e k-ε RNG, foram incapazes de predizer o efeito da posição apropriadamente. Entretanto, o modelo k-ε RNG apresentou melhor comportamento, pois o seu uso resultou em soluções com valores de temperatura intermediários aos experimentais / Abstract: The purpose of this work is the study of the numerical and experimental viscous incompressible steady flow with heat transfer into a narrow channel containing a rectangular array of obstacles. The experimental approach involves determining the coefficient of heat transfer and temperature measurements in thermal wake for comparison with the results obtained in numerical simulations. For the numerical analysis we use the commercial program of fluid mechanics and heat transfer computational ICEPAK™. We confirmed that in the last lines of the array the biggest is the heat transfer by forced convection. We determined the average heat transfer coefficients (with uncertainty between 6 and 15%) and found that the effect of the position decreases as flow speed increases. We use in the simulations the k-ε turbulence model and the k-ε RNG turbulence model. We conclude that both turbulence models used were unable to predict the effect of the position properly. However, the k-ε RNG model showed better behavior. The numerical temperatures with this model were consistent to the experimental temperature / Orientador: João Batista Campos Silva / Coorientador: Amarildo Tabone Paschoalini / Coorientador: Marcio Antonio Bazani / Banca: Ricardo Alan Verdú Ramos / Banca: Marcio Higa / Mestre

La modélisation des transferts thermiques en écoulements diphasiques est l'une des pierres angulaires de l'étude de la sûreté des réacteurs nucléaires. À l'échelle du réacteur, elle repose sur des corrélations expérimentales. L'utilisation croissante de la mécanique des fluides numérique pour les études de sûreté renforce la demande d'expertise dans les outils de simulation, en particulier du point de vue de la modélisation. En soutien aux modèles moyennés à deux fluides, nous souhaitons apporter des informations de fermetures locales pour considérer la physique des transferts interfaciaux et les effets 3D. Pour cela, comme la résolution directe des équations de bilan locales par SND est trop coûteuse, nous souhaitons développer un outil de SGE diphasique pour modéliser les petites échelles turbulentes et les petites déformations interfaciales. Comme le changement de phase est à l'origine de l'écoulement diphasique pour les applications visées, nous étendons dans ce mémoire le modèle Interfaces and Subgrid-Scales (ISS, Toutant et al., 2009a) aux interfaces avec changement de phase, pour lesquelles l'hypothèse de continuité de la vitesse à l'interface n'est plus valable. Le suivi explicite des interfaces permet d'évaluer précisément les transferts comme le taux de transfert de masse. Dans un premier temps, nous établissons une description mésoscopique du problème où l'interface est diffuse en filtrant les équations locales instantanées et en modélisant les transferts sous-filtres aux interfaces. Les principales difficultés de modélisations proviennent (i) de la détermination de la vitesse de l'interface, (ii) de l'effet de la discontinuité des vitesses sur les modèles sous-maille, (iii) de la discontinuité du flux et (iv) de la condition de saturation de l'interface. Les modèles proposés sont qualifiés a priori en observant leur prédiction par filtrage explicite de solutions de SND. Dans un deuxième temps, nous établissons un système macroscopique discontinu équivalent au problème diffus pour bénéficier de l'expertise acquise pour les méthodes numériques de SND. Aux interfaces, les modèles sous-maille sont concentrés pour modifier les conditions de raccord entre les phases. Les conditions de saut ainsi déterminées montrent que la vitesse de l'interface est affectée par la courbure et par le saut de vitesse. Un saut de vitesse tangentielle est introduit pour modéliser la couche limite dynamique. Sur le plan thermique, nous retrouvons la condition de saturation caractéristique du changement de phase ; le taux de changement de phase ne dépend plus uniquement du saut de flux conductif mais, pour pallier la sous-résolution de la couche limite thermique au voisinage de l'interface, nous proposons de lui ajouter la contribution sous-maille des corrélations vitesse/température. Comme en SGE monophasique, le gain apporté par la modélisation ISS permet d'envisager l'utilisation de simulations fines pour des problèmes appliqués. C'est la première étape d'une démarche multi-échelle pour fournir des fermetures aux modèles moyennés à deux fluides. Nous illustrons son potentiel sur une SND multi-bulles complexe.