Orientador: Vicente Luiz Scalon / Banca: Helio Aparecido Navarro / Banca: Sergio Rodrigues Fontes / Resumo: A fluidodinâmica computacinal (CFD) tem sido utilizada, estudadda e implementada ao longa das duas últimas décadas na solução dos mais diversos problemas de engenharia. O princípio básico desta ciência é a aplicação de métodos numéricos em problemas que envolvam mecânica dos fluidos. Nesse contexto, este trabalho utiliza essa técnica para analisar o comportamento de um fluido incompreensível, que se encontra numa cavidade cilíndrica fechada onde as faces inferior e superior são adiabáticas e as superfícies laterais se encontram em diferentes temperaturas. Os perfis de velocidade e temperatura resultantes - ocasionados pela convecção natural - serão avaliados em todo o domínio do problema. Existe uma série de técnicas para a solução de problemas envolvendo escoamentos, sendo as mais comuns as que se utilizam do "Esquema de Passo Fracionado" proposto por Chorin no final da década de 60. Dentre as diversas soluções que se utilizam desta técnica, este trabalho optou pelo uso do método das características e do algoritmo CBS de solução proposto por Zienkiewicz e Codina (1995). Para a implementação do algoritmo de solução do problema proposto foi realizada uma discretização geral através do método dos elementos finitos usando-se de uma malha formada por elementos bilineares. A solução foi obtida a partir de um ambiente matemático adequado, a GNU-Octave (2008). Os resultados foram analisados para diferentes razões de curvatura, números de Rayleigh e métodos de solução, sendo plotados para as suas diversas variáveis buscando descrever o comportamento do fenônemo / Abstract: The Computational Fluid Dynamics (CFD) has been used, studied and performed through the last two decades to solve the series of problems in Engineering. The most basic aim of this science is the appliance of numerical methods in cases that envolve fluid mechanics. In this context, this work uses this technic to analyze the behaviour of an incompressible fluid, which is found in a closed cylindrical cavity, a place where the inferior and superior surfaces are adiabatic and the lateral surfaces are shown in different temperatures. The resultant profiles of speed and temperature - induced by the free convection - are going to be appraised in all the dominion of the problem. There is a set of technics to solve the problems which involve the drainage, but the most usual are those which use the techic "Fractional Step Method" offered by Chorin in the final of 60s. Among the several solutions that are solved through this technic, this research used the characteristics method and of the CBS algorthm, offered by Zienkiewicz e Codina (1995). For the implementation of the algorithm, it was realized a general discretization through the finite elements method, making use of a loop formed by bilinear elements. The resolution was obtained from an adequated mathematical ambient, the GNU-Octave (2008). The results were analysed for different curvature ratios, Rayleigh numbers and methods of solution, being plotted for its different variables searching to describe the behavior of the phenomenon / Mestre

Made available in DSpace on 2014-06-11T19:25:27Z (GMT). No. of bitstreams: 0 Previous issue date: 2010-08-20Bitstream added on 2014-06-13T19:32:32Z : No. of bitstreams: 1 jordam_a_me_bauru.pdf: 3535049 bytes, checksum: a22a97d375d5d91449be8bd46442ee9e (MD5) / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) / A fluidodinâmica computacinal (CFD) tem sido utilizada, estudadda e implementada ao longa das duas últimas décadas na solução dos mais diversos problemas de engenharia. O princípio básico desta ciência é a aplicação de métodos numéricos em problemas que envolvam mecânica dos fluidos. Nesse contexto, este trabalho utiliza essa técnica para analisar o comportamento de um fluido incompreensível, que se encontra numa cavidade cilíndrica fechada onde as faces inferior e superior são adiabáticas e as superfícies laterais se encontram em diferentes temperaturas. Os perfis de velocidade e temperatura resultantes - ocasionados pela convecção natural - serão avaliados em todo o domínio do problema. Existe uma série de técnicas para a solução de problemas envolvendo escoamentos, sendo as mais comuns as que se utilizam do Esquema de Passo Fracionado proposto por Chorin no final da década de 60. Dentre as diversas soluções que se utilizam desta técnica, este trabalho optou pelo uso do método das características e do algoritmo CBS de solução proposto por Zienkiewicz e Codina (1995). Para a implementação do algoritmo de solução do problema proposto foi realizada uma discretização geral através do método dos elementos finitos usando-se de uma malha formada por elementos bilineares. A solução foi obtida a partir de um ambiente matemático adequado, a GNU-Octave (2008). Os resultados foram analisados para diferentes razões de curvatura, números de Rayleigh e métodos de solução, sendo plotados para as suas diversas variáveis buscando descrever o comportamento do fenônemo / The Computational Fluid Dynamics (CFD) has been used, studied and performed through the last two decades to solve the series of problems in Engineering. The most basic aim of this science is the appliance of numerical methods in cases that envolve fluid mechanics. In this context, this work uses this technic to analyze the behaviour of an incompressible fluid, which is found in a closed cylindrical cavity, a place where the inferior and superior surfaces are adiabatic and the lateral surfaces are shown in different temperatures. The resultant profiles of speed and temperature - induced by the free convection - are going to be appraised in all the dominion of the problem. There is a set of technics to solve the problems which involve the drainage, but the most usual are those which use the techic Fractional Step Method offered by Chorin in the final of 60s. Among the several solutions that are solved through this technic, this research used the characteristics method and of the CBS algorthm, offered by Zienkiewicz e Codina (1995). For the implementation of the algorithm, it was realized a general discretization through the finite elements method, making use of a loop formed by bilinear elements. The resolution was obtained from an adequated mathematical ambient, the GNU-Octave (2008). The results were analysed for different curvature ratios, Rayleigh numbers and methods of solution, being plotted for its different variables searching to describe the behavior of the phenomenon

Nous présentons dans cette étude le développement de la convection à partir de différentes perturbations de l'état conductif d'une couche fluide confinée dans une cavité cylindrique, chauffée par le bas et avec une surface supérieure libre. La discrétisation spatiale du domaine repose sur la méthode des éléments spectraux et les itérations temporelles sont assurées par une méthode splitting.Au déclenchement de la convection, les structures convectives correspondent à des modes de Fourier, et les seuils critiques dépendent du rapport de forme de la cavité, et des nombres de Biotet de Marangoni qui caractérisent la surface libre. Les transitions d'écoulements au-delà du seuil primaire sont caractérisées quantitativement en fonction du nombre de Rayleigh pour différentes valeurs du nombre de Biot et Ma = 0. Les résultats présentés sont obtenus en résolvant l'ensemble des équations non-linéaires de conservation à travers une méthode de continuation. Lorsque la convection se déclenche sous la forme d'un mode axisymétrique m = 0, l'évolution non-linéaire montre la coexistence de différentes structures convectives, des structures axisymétriques avec écoulement montant ou descendant au centre de la cavité et des structures correspondant à des combinaisons de modes qui apparaissent sur des branches secondaires sous-critiques.L'action d'un champ magnétique constant est ensuite étudiée pour des fluides conducteurs dans une même configuration comprenant une surface supérieure libre. Nous montrons l'effet stabilisateur du champ magnétique sur les seuils primaires ainsi que son action sélective sur les différents modes de convection. Nous analysons l'évolution des structures convectives au delà de ces seuils et montrons comment le champ magnétique modifie les transitions entre ces structures.En soumettant le bain fondu à un champ magnétique tournant, le mouvement de rotation du fluide se superpose aux mouvements de convection thermique et on observe une diminution des fluctuations de température et un retard du déclenchement de l'instabilité de Rayleigh-Bénard(lorsque les deux parois haut/bas du bain sont rigides). La rotation influe sur ce déclenchement qui de stationnaire devient oscillatoire, à l'exception du mode m = 0 de Fourier, pour qui la transition reste stationnaire jusqu'à une certaine valeur critique du nombre de Taylor magnétique.La dynamique de l'écoulement axisymétrique de part et d'autre de cette valeur critique sera étudiée en détail. / The growth of thermal convection out of different perturbations of the conductive base state is investigated using a spectral element time-stepping code. The fluid is subject to a vertical heat transfer in a cylindrical cavity with an upper free surface corresponding to the so-called Rayleigh-Bénard-Marangoni situation and the heat exchange through the free surface is evaluated via the Biot number. The results of the stability diagrams show that the evolution of the primary thresholds are largely influenced by the Biot number, the Marangoni number, and the aspect ratio of the cavity. Flow transitions are elucidated in quantitative detail as a function of the Rayleigh number for different Biot numbers in the tension free limit Ma = 0. The results shown are obtained by solving the full nonlinear field equations numerically among a continuation method. When an axisymmetric m = 0 Fourier mode is obtained at onset, the non-linear evolution shows the coexistence of different convective structures, the axisymmetric structures with up-ow or down-ow at the center and mixed-mode structures which appear on secondary subcritical branches. The action of a constant magnetic field is then considered for melts in the same type of configuration with an upper free surface. We show the global stabilizing effect of the magnetic field on the primary bifurcation thresholds and the selective effect on the different instability modes. We also analyze the evolution of the convective structures above the thresholds and show how the magnetic field modifies the transitions between these structures. When applying a magnetic body forcing in the azimuthal direction (RMF), one can damp the unavoidable thermal fluctuations inside the melt and delay the transition to the Rayleigh-Bénard instability (for rigid-rigid circular plates at top and bottom). The rotation effect also changes the transitions from steady to oscillatory, except for the m = 0 Fourier mode where the transitionis first steady until a critical Taylor number and then becomes oscillatory. The dynamics of the transitions to the axisymmetric flow, below and above this value of critical magnetic Taylor number, is particularly interesting and will be described.

Le bruit de cavité est un phénomène très fréquent dans le domaine des transports aériens.Il survient notamment lors de l’approche à l’atterrissage, où des interactions entre la cellule de l’aéronef et l’écoulement sont à l’origine de fortes émissions tonales. Il devient dès lors une source de pollution acoustique non-négligeable pour les populations résidant à proximité de zones aéroportuaires. Les études numériques et expérimentales décrites jusqu’à présent dans la littérature abordent essentiellement le cas des cavités rectangulaires. Pourtant, les cavités rencontrées en pratique dans l’industrie aéronautique impliquent des géométries souvent plus complexes. Lorsque ces cavités sont soumises à une excitation de nature aérodynamique, leur spécificité géométrique conduit le plus souvent à des réponses acoustiques assez éloignées des estimations issues de modèles académiques construits sur l’observation de cavités rectangulaires. Quelques travaux seulement abordent le cas des cavités cylindriques.Ce travail est consacré à l’étude aéroacoustique des cavités cylindriques, à l’initiative d’Airbus. Il s’inscrit dans le cadre du projet AEROCAV soutenu par la Fondation de Recherche pour l’Aéronautique & l’espace (FRAE). Son objectif est de déterminer les mécanismes impliqués dans les émissions acoustiques intenses et tonales pour les configurations étudiées.Une première partie présente les résultats expérimentaux issus des campagnes de mesures menées dans la soufflerie anéchoïque du Centre Acoustique du LMFA et de l’école Centrale de Lyon. Un modèle semi-empirique, reposant sur l’hypothèse d’une résonance acoustique pilotée par les instabilités présentes dans la couche de cisaillement à l’ouverture de la cavité,est construit à partir du modèle d’Elder (1978). Le modèle permet d’estimer les fréquences susceptibles de dominer l’acoustique rayonnée en champ lointain à partir de la donnée du champ moyen de vitesse longitudinale, que l’on mesure dans le plan de l’écoulement par Vélocimétrie par Imagerie des Particules (PIV).Une seconde partie est destinée au calcul direct du bruit rayonné par un écoulement laminaire ou turbulent affleurant une cavité cylindrique de référence. Il consiste à calculer le champ acoustique directement à partir de la résolution des équations tridimensionnelles de la mécanique des fluides. Le solver Alesia est présenté dans une version modifiée et adaptée à la mise en oeuvre d’une approche multidomaine d’ordre élevé faisant intervenir plusieurs maillages se recouvrant. Des techniques d’interpolation sont spécifiquement développées en vue d’assurer une communication bidirectionnelle entre les différents maillages, malgré des contraintes géométriques fortes. Un modèle d’excitation de l’écoulement est aussi développé afin de disposer de fluctuations dans l’écoulement incident, pour le cas turbulent. Ces deux points font l’originalité des calculs réalisés.Les simulations, menées sur une cavité de rapport d’aspect géométrique égal à 1 et soumise à un écoulement incident à Mach 0.2, montrent que le rayonnement acoustique peut être fidèlement reproduit numériquement. La couche de cisaillement est caractérisée par la présence de deux larges structures tourbillonnaires s’amplifiant lors de leur convection. Leur présence s’accompagne de fortes fluctuations de vitesse à l’origine d’un débit aérodynamique de fluide à l’ouverture qui excite la cavité acoustiquement. Une résonance forcée s’établit dans celle-ci, excitant la couche de mélange au voisinage du point de séparation. Ce couplage auto-entretenu est à l’origine du rayonnement acoustique intense et fortement tonal de la cavité. Il s’établit à une fréquence proche de la fréquence prédite par le modèle semi-empirique développé. / Cavity noise is a very frequent phenomenon in air transport. It occurs in particular during landing approaches, where airframe–flow interactions are responsible for strong tonal emissions. Accordingly, it turns to be a non negligible source of acoustic pollution for populations living near airport areas. Numerical and experimental studies reported in the literature tackle essentially the case of rectangular cavities. Nevertheless, cavities may often exhibit more complex shapes in practice. When subject to aerodynamic excitations, and because of their geometrical specificity, these cavities may have acoustic responses which can be rather far from estimations resulting from academic models designed for rectangular cavities. Only asmall number of studies tackle the case of cylindrical cavities.The present work requested by Airbus is dedicated to the study of aeroacoustics in cylindrical cavities. This work was been supported by the Fondation de Recherche pour l’Aéronautique& l’Espace (FRAE) under contract reference AEROCAV. It aims at discerning the mechanisms responsible for strong and tonal acoustic emissions for the studied configurations.Experimental data resulting from measurements performed in the anechoic wind-tunnel of the Centre Acoustique at ´Ecole Centrale de Lyon are presented in a first part. A semi-empirical model based on the hypothesis of a shear-layer driven acoustic resonance is constructed from the Elder model (1978). The model provides an estimation of the frequences which are likely to be predominant in the far-field acoustics, given the mean streamwise velocity field, currently measured in the flow plane by Particle Image Velocimetry (PIV).A second part deals with the direct computation of the noise radiated by a laminar or turbulent grazing flow over a standard cylindrical cavity. The method consists in the calculationof the acoustic field directly from the resolution of the tridimensional Navier–Stokes equations. The Alesia solver is presented in a modified form, adapted to the implementationof a high-order chimera method involving several overlapping grids. Interpolation techniques have been specifically developed to achieve a bidirectional communication between the meshes in spite of strong geometrical constraints. A flow excitation model has also been constructed in order to obtain fluctuations into the incoming flow in the turbulent case. These two last points make the present computations original. The simulations, which are performed on a cavity of geometric ratio taken as 1 and subject to a grazing flow of Mach 0.2, reveal that it is possible to retrieve the radiated noise numerically with high fidelity. They indicate the presence of two large amplifying vortices in the shearlayer. These vortices go with strong velocity fluctuations giving rise to an inflow of fluid at the cavity mouth which excites the cavity acoustically. A forced acoustic resonance occurs into the cavity, then destabilises the shear layer near the separation point. This self-sustained coupling is responsible for strong tonal radiations from the cavity. The frequency of the radiated noise is close to the one predicted by the semi-empirical model.