Orientador: Evandro Conforti / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Eletrica e de Computação / Made available in DSpace on 2018-08-03T17:27:13Z (GMT). No. of bitstreams: 1 AnguloButron_JoseBoris_M.pdf: 3718604 bytes, checksum: 3c919d32bad48bb481ab0a8daa87bc5e (MD5) Previous issue date: 2003 / Mestrado

Une caractéristique remarquable de la convection de Rayleigh-Bénard qui concerne une couche de fluide horizontale chauffée par le bas et refroidie par le haut, est l’établissement spontané de l’ordre spatial et la formation d’une circulation cohérente à grande échelle. Comment les différents facteurs, tels que la géométrie du domaine et les conditions limites, influencent l’écoulement à grande échelle, restent une question ouverte. Malgré sa simplicité apparente, la convection de Rayleigh-Bénard présente une dynamique à grande échelle incroyablement riche et complexe, tels que des modes de torsion et du battement, la rotation du plan et des cessations de la circulation, qui coexistent souvent et se concourent. Une approche couramment utilisée dans l’étude des cessations, consiste à contraindre la circulation à grande échelle à un plan en limitant le domaine fluide à une cellule carrée (2D) ou à une cellule rectangulaire mince (quasi-2D). Cependant, il n’est pas tout à fait clair si les retournements 2-D et (quasi-)2-D correspondent au même phénomène. Le présent document est consacré à l’étude des modes d’écoulement à grande échelle dans la convection turbulente de Rayleigh-Bénard et de l’influence exercée par différents facteurs sur les structures d’écoulement et sur leur évolution temporelle. La caractérisation proposée combine une analyse statistique avec une approche physique s’appuyant sur le moment angulaire ainsi que sur les énergies cinétiques et potentielles pour mettre en évidence les mécanismes physiques sous-jacents. Nous essayons ensuite de relier ces mécanismes à chacun des modes d’écoulement distinctifs observés et à leur évolution. / One outstanding feature of the Rayleigh-Bénard problem which concerns a horizontal fluid layer heated from below and cooled from above, is the spontaneous establishment of spatial ordering and the formation of a coherent large-scale circulation. How different factors, such as the domain geometry and boundary conditions, influence the sizes and shapes of the large-scale flow remains an open question. Despite its apparent simplicity, Rayleigh-Bénard convection exhibits some incredibly rich and complex large-scale dynamics such as torsional modes, rotation, sloshing, and cessations, which often coexist and compete to each other. One approach, commonly used in the study of cessations is to constrain the large scale circulation to a plane by restricting the fluid domain to a (2-D) square cell or to a slim rectangular cell of small aspect ratio in the transversal direction. However, it is not entirely clear whether the 2-D and (quasi-)2-D reversals correspond to the same phenomenon. The present document is dedicated to the study of the large-scale flow patterns in turbulent Rayleigh-Bénard convection, and of the influence exerted by different factors on the flow structures and on their temporal evolution. The proposed characterization combines a statistical analysis with a physical approach relying on the angular momentum as well as the kinetic and potential energies to highlight the underlying physical mechanisms. We subsequently attempt to tie these mechanisms together to each of the distinctive flow patterns observed and to their evolution.

Rayleigh-Taylor instability has been an area of active research in fluid dynamics for the last twenty years, but relatively little attention has been paid to the dynamics of problems where Rayleigh-Taylor instability plays a role, but is only one component of a more complex system. Here, Rayleigh-Taylor instability between miscible fluids is examined in situations where it is confined by various means: by geometric restriction, by penetration into a stable linear stratification, and by impingement on a stable density interface. Water-based experiments are modelled using a variety of techniques, ranging from simple hand calculation of energy exchange to full three-dimensional numerical simulation. Since there are well known difficulties in modelling unconfined Rayleigh-Taylor instability, the confined test cases have been sequenced to begin with dynamically simple benchmark systems on which existing modelling approaches perform well, then they progress to more complex systems and explore the limitations of the various models. Some work on the phenomenology of turbulent mixing is also presented, including a new experimental technique that allows mixed fluid to be visualised directly, and an analysis of energy transport and mixing efficiency in variable density flows dominated by mixing.

No description available.

La vaporisation rapide du sodium liquide surchauffé est supposée être à l’origine des arrêts automatiques pour réactivité négative du réacteur Phénix.Un dispositif expérimental a été mis en œuvre pour reproduire la détente d'un gaz pressurisé, repoussant un liquide dans un canal de section rectangulaire très allongée.L’interface qui sépare les deux fluides, initialement plate, ondule du fait d'instabilités de Rayleigh-Taylor dont le caractère 2D est garanti par le rapport d'aspect de la section du canal. L’aire interfaciale augmente d'un facteur 50.L’expansion du gaz peut être divisée en deux phases principales : une phase dite « de Rayleigh-Taylor » (linéaire et non-linéaire) et une phase dite « à multi-structures » (transitionnelle et chaotique). La première est caractérisée par la dynamique de l'interface et l’aire interfaciale qui en résulte est proportionnelle à l’amplitude des ondulations. La deuxième est influencée par le comportement des structures liquides, dispersées dans la matrice gazeuse et l’aire interfaciale est alors proportionnelle au nombre de structures.La distribution de fraction volumique suggère un modèle d’écoulement composé de trois régions : une région où la frontière des bulles est clairement définie et régulière, une région compartimentée par des membranes liquides issues des frontières des bulles, une région diphasique formée de la queue de ces structures. L’analyse de sensibilité à la tension superficielle confirme que plus la tension est faible, plus les interfaces sont instables. Les ondes sont plus prononcées et plus de structures sont produites, ce qui conduit à une majoration du taux de production de l’aire interfaciale. / The sharp vaporization of superheated liquid sodium is investigated. It is suspected to be at the origin of the automatic shutdown for negative reactivity, occurred in the Phénix reactor at the end of the eighties.An experimental apparatus has been designed and operated to reproduce the expansion of overpressurized air, superposed to water in a narrow vertical rectangular section channel.When expansion begins, the initial flat interface separating the two fluids becomes corrugated under the development of two-dimensional Rayleigh-Taylor instabilities. The interface area increases significantly and becomes even 50 times larger than the initial value. Since the channel is very narrow, instabilities along the channel depth do not develop.The gas expansion in a narrow channel can be divided into two main phases: Rayleigh-Taylor (linear and non-linear) and multi-structures (transition and chaotic) phases. The former is characterized by the dynamic of corrugated profile and the interface area results proportional to the amplitude of corrugation The latter is influenced by the behavior of the liquid structures dispersed in gas matrix and the interface area is mainly proportional to the number of liquid structures.The distribution of volume fraction suggests a model of channel flow consisting of three regions: the regular profile of peaks, the spike region and the structures tails. The analysis of sensibility to surface tension confirms that, with a lower surface tension, the fluids configuration is more unstable. The interface corrugations are more pronounced and more structures are produced, leading to a higher increment of the interface area.

Une étude des matériaux à base de polymères et de cristaux liquides de type PDLC (pour Polymer Dispersed Liquid Crystals), élaborés par la méthode de séparation de phases induite par rayonnement ultraviolet (UV), a été effectuée pour des mélanges comportant le cristal liquide nématique E7 et le monomère tripropylèneglycoldiacrylate (TPGDA). Ces matériaux possèdent des fonctionnalités électro-optiques intéressantes notamment dans les vitrages à transparence contrôlée. Une analyse par spectroscopie infrarouge des réseaux de polymères élaborés a permis de déduire qu’une conversion quasi-totale des fonctions réactives de type acrylique du monomère est obtenue pour un mélange contenant 70% de cristal liquide et 30% de monomère. Ces matériaux ont été, ensuite, étudiés par différentes techniques de caractérisation, telles que la calorimétrie différentielle, la microscopie optique à lumière polarisée, la microscopie électronique à balayage, et la spectroscopie UV-visible, qui ont permis d’obtenir des renseignements sur les propriétés thermophysiques, morphologiques, et spectrales de ces systèmes.Une modélisation de la réponse électro-optique des films PDLC a été effectuée en utilisant un modèle simple, basé sur une hiérarchie de paramètres d’ordre. Ce modèle, dont les calculs sont effectués dans deux approximations théoriques (RGA et ADA), a donné une description convenable du comportement électro-optique de ces systèmes complexes. L’étude menée sur la confrontation de ce modèle à l’expérience a permis l’obtention d’un certain nombre de résultats intéressants qui sont utiles à la compréhension et à l’amélioration de la réponse électro-optique des films PDLC. / A study of materials based on polymers and liquid crystals of type PDLC (for Polymer Dispersed Liquid Crystals), elaborated by the method of phase separation induced by ultraviolet radiation (UV), was carried out for mixtures containing the nematic liquid crystal E7 and the monomer tripropyleneglycoldiacrylate (TPGDA). These materials possess interesting electro-optical features in particular for privacy windows with controlled transparency.An analysis by infrared spectroscopy of the elaborated polymer networks allowed to deduce that a quasi-total conversion of the reactive acrylic functions of the monomer was obtained for a mixture containing 70 % of liquid crystal and 30 % of monomer. These materials were, then, studied by various techniques of characterization, such as differential scanning calorimetry, polarized optical microscopy, scanning electron microscopy, and UV-visible spectroscopy, who allowed to obtain informations on the thermophysical, morphological, and spectroscopical properties of these systems.An attempt was made to rationalize the electro-optical response of PDLC films by using a simple model, based on a hierarchy of order parameters. Using the theoretical description from two theoretical approximations (RGA and ADA), this model gave a proper description of the electro-optical behaviour of these complex systems. The study led on the confrontation of this model to experimental results, allowing to obtain some interesting results which are useful in understanding and improvement of the electro-optical response of PDLC films.

Dans cette thèse, nous nous intéressons au problème de tester en grande dimension l'uniformité sur la sphère-unité $S^{p_n-1}$ (la dimension des observations, $p_n$, dépend de leur nombre, $n$, et être en grande dimension signifie que $p_n$ tend vers l'infini en même temps que $n$). Nous nous restreignons dans un premier temps à des contre-hypothèses ``monotones'' de densité croissante le long d'une direction ${\pmb \theta}_n\in S^{p_n-1}$ et dépendant d'un paramètre de concentration $\kappa_n>0$. Nous commençons par identifier le taux $\kappa_n$ auquel ces contre-hypothèses sont contiguës à l'uniformité ;nous montrons ensuite grâce à des résultats de normalité locale asymptotique, que le test d'uniformité le plus classique, le test de Rayleigh, n'est pas optimal quand ${\pmb \theta}_n$ est connu mais qu'il le devient à $p$ fixé et dans le cas FvML en grande dimension quand ${\pmb \theta}_n$ est inconnu.Dans un second temps, nous considérons des contre-hypothèses ``axiales'', attribuant la même probabilité à des points diamétralement opposés. Elles dépendent aussi d'un paramètre de position ${\pmb \theta}_n\in S^{p_n-1}$ et d'un paramètre de concentration $\kappa_n\in\R$. Le taux de contiguïté s'avère ici plus élevé et suggère un problème plus difficile que dans le cas monotone. En effet, le test de Bingham, le test classique dans le cas axial, n'est pas optimal à ${\pmb \theta}_n$ inconnu et $p$ fixé, et ne détecte pas les contre-hypothèses contiguës en grande dimension. C'est pourquoi nous nous tournons vers des tests basés sur les plus grande et plus petite valeurs propres de la matrice de variance-covariance et nous déterminons leurs distributions asymptotiques sous les contre-hypothèses contiguës à $p$ fixé.Enfin, à l'aide d'un théorème central limite pour martingales, nous montrons que sous certaines conditions et après standardisation, les statistiques de Rayleigh et de Bingham sont asymptotiquement normales sous l'hypothèse d'invariance par rotation des observations. Ce résultat permet non seulement d'identifier le taux auquel le test de Bingham détecte des contre-hypothèses axiales mais aussi celui auquel il détecte des contre-hypothèses monotones. / In this thesis we are interested in testing uniformity in high dimensions on the unit sphere $S^{p_n-1}$ (the dimension of the observations, $p_n$, depends on their number, and high-dimensional data are such that $p_n$ diverges to infinity with $n$).We consider first ``monotone'' alternatives whose density increases along an axis ${\pmb \theta}_n\in S^{p_n-1}$ and depends on a concentration parameter $\kappa_n>0$. We start by identifying the rate at which these alternatives are contiguous to uniformity; then we show thanks to local asymptotic normality results that the most classical test of uniformity, the Rayleigh test, is not optimal when ${\pmb \theta}_n$ is specified but becomes optimal when $p$ is fixed and in the high-dimensional FvML case when ${\pmb \theta}_n$ is unspecified.We consider next ``axial'' alternatives, assigning the same probability to antipodal points. They also depend on a location parameter ${\pmb \theta}_n\in S^{p_n-1}$ and a concentration parameter $\kappa_n\in\R$. The contiguity rate proves to be higher in that case and implies that the problem is more difficult than in the monotone case. Indeed, the Bingham test, the classical test when dealing with axial data, is not optimal when $p$ is fixed and ${\pmb \theta}_n$ is not specified, and is blind to the contiguous alternatives in high dimensions. This is why we turn to tests based on the extreme eigenvalues of the covariance matrix and establish their fixed-$p$ asymptotic distributions under contiguous alternatives.Finally, thanks to a martingale central limit theorem, we show that, under some assumptions and after standardisation, the Rayleigh and Bingham test statistics are asymptotically normal under general rotationally symmetric distributions. It enables us to identify the rate at which the Bingham test detects axial alternatives and also monotone alternatives. / Doctorat en Sciences / info:eu-repo/semantics/nonPublished

The Richtmyer-Meshkov instability occurs when density-stratified interfaces are impulsively accelerated, typically by a shock wave. We present a numerical method to simulate the Richtmyer-Meshkov instability in cylindrical geometry. The ideal MHD equations are linearized about a time-dependent base state to yield linear partial differential equations governing the perturbed quantities. Convergence tests demonstrate that second order accuracy is achieved for smooth flows, and the order of accuracy is between first and second order for flows with discontinuities. Numerical results are presented for cases of interfaces with positive Atwood number and purely azimuthal perturbations. In hydrodynamics, the Richtmyer-Meshkov instability growth of perturbations is followed by a Rayleigh-Taylor growth phase. In MHD, numerical results indicate that the perturbations can be suppressed for sufficiently large perturbation wavenumbers and magnetic fields.

We investigate the Richtmyer-Meshkov instability (RMI) that occurs when an incident shock impulsively accelerates the interface between two different fluids. RMI is important in many technological applications such as Inertial Confinement Fusion (ICF) and astrophysical phenomena such as supernovae. We consider RMI in the presence of the magnetic field in converging geometry through both simulations and analytical means in the framework of ideal magnetohydrodynamics (MHD). In this thesis, we perform linear stability analyses via simulations in the cylindrical geometry, which is of relevance to ICF. In converging geometry, RMI is usually followed by the Rayleigh-Taylor instability (RTI). We show that the presence of a magnetic field suppresses the instabilities. We study the influence of the strength of the magnetic field, perturbation wavenumbers and other relevant parameters on the evolution of the RM and RT instabilities. First, we perform linear stability simulations for a single interface between two different fluids in which the magnetic field is normal to the direction of the average motion of the density interface. The suppression of the instabilities is most evident for large wavenumbers and relatively strong magnetic fields strengths. The mechanism of suppression is the transport of vorticity away from the density interface by two Alfv ́en fronts. Second, we examine the case of an azimuthal magnetic field at the density interface. The most evident suppression of the instability at the interface is for large wavenumbers and relatively strong magnetic fields strengths. After the shock interacts with the interface, the emerging vorticity breaks up into waves traveling parallel and anti-parallel to the magnetic field. The interference as these waves propagate with alternating phase causing the perturbation growth rate of the interface to oscillate in time. Finally, we propose incompressible models for MHD RMI in the presence of normal or azimuthal magnetic field. The linearized equations are solved numerically using inverse Laplace transform. The incompressible models show that the magnetic field suppresses the RMI, and the mechanism of this suppression depends on the orientation of the initially applied magnetic field. The incompressible model agrees reasonably well with compressible linear simulations.

In this thesis, the dynamics of an isolated cavitation bubble submerged in a steady flow is studied numerically. A Lagrangian-Eulerian approach is considered, in which properties of the fluid are computed first by means of Eulerian methods (in this study the commercial CFD software Ansys Fluent 19 was used) and the trajectory of the bubble is then computed in a Lagrangian fashion, i.e. the bubble is considered as a small particle moving relative to the fluid, due to the effect of several forces depending on fluid's pressure field, fluid's velocity field and bubble's radius. Bubble's radius dynamics, modeled by Rayleigh-Plesset equation, has a big influence on its kinetics, so a special attention is given to it. Two study cases are considered. The first one, motivated by acoustic cavitation is concerned with the response of the bubble's radius in a static flow under the influence of an oscillatory pressure field, the second one studies the trajectory of the bubble submerged in a fluid passing by a Venturi tube and a sharp-edged orifice plate.