Abstract An experimental study is conducted to determine the velocity fields, from development to steady state, in a square enclosure due to movement of a constant velocity lid using Particle Image Velocitmetry (PIV). Experiments were conducted with water, seeded with hollow glass sphere particles 10 microns in diameter, at three different lid velocities leading to Reynolds numbers in the high laminar to transitional range. Driven Cavity Flow is a classic fluid dynamics case often used for benchmarking of computational codes. Previous work has primarily focused on improving computational codes, experimental work is lacking and focused on obtaining steady state readings. The test cavity is 1 inch (25.4mm) high by 1 inch (25.4 mm) wide leading to an aspect ratio of 1.0. The depth is taken to be 5 (127mm) inches to reduce the three dimensional effects. Readings are taken from development to steady state allowing for a full spectrum of flow characteristics. PIV technique is successful in capturing the development of driven cavity flow. Circulation is shown to increase strength with time and Reynolds number. PIV capture and processing settings are determined. Keywords: Driven Cavity Flow, Particle Image Velocimeter (PIV)

Les écoulements instationnaires au-dessus de surfaces discontinues produisent d'important bruit aérodynamique. L'objectif de ce travail de thèse est l'étude aéroacoustique d'écoulement au-dessus de cavités bidimensionnelles rectangulaires, et de trouver des stratégies de réduction du bruit. Des simulations numériques directes des équations bidimensionnelles de Navier-Stokes compressibles ont été réalisées afin d'étudier l'influence des conditions initiales sur le mode d'oscillation de l'écoulement pour des cavités profonde et peu profonde. Les résultats montrent que dans le cas de cavités profondes, l'écoulement oscille selon un régime de couche de cisaillement suivant le second mode de Rossiter, et ce quelle que soit la condition initiale choisie. En revanche, dans le cas de cavités peu profondes, le régime d'oscillation observé peut être en couche de cisaillement ou bien en mode de sillage suivant la condition initiale choisie. Une analyse de sensibilité d'écoulement dans le cas de cavités profondes a été réalisé en utilisant une méthode adjointe. Les équations adjointes ont été forcées par une perturbation localisée sinusoïdale soit de la quantité de mouvement suivant x adjointe (au voisinage de la couche de cisaillement), soit de la densité adjointe (loin de la cavité). Les résultats désignent une région de l'écoulement très sensible à l'ajout de masse, et localisée au voisinage du coin supérieur amont de la cavité. Par conséquent, un actionneur de type soufflage/aspiration placé au bord d'attaque de la cavité agira sur les fluctuations de quantité de mouvement suivant x au voisinage de la couche de cisaillement et sur les fluctuations de pression au loin. / The unsteady flow over surface discontinuities produces high aerodynamic noise. The aim of this thesis is to study the aeroacoustics of two-dimensional rectangular cavities and to find strategies for noise reduction. Direct Numerical Simulation of the compressible Navier-Stokes equations is performed to investigate the influence of the initial condition on the oscillation modes in deep and shallow cavities. Results show that the deep cavity oscillates in shear layer regime at the second Rossiter mode regardless of the initial condition. On the other hand different initial conditions lead to a shear layer or wake mode in the shallow cavity case. A sensitivity analysis of the deep cavity is done by the use of adjoint methods. Local sinusoidal perturbations of x-momentum and density are applied to the adjoint equations. The results show a high sensitivity region to mass injection at the upstream corner. Therefore an actuator placed at the leading edge will modify the velocity fluctuations reaching the trailing edge and hence the pressure fluctuations in the far-field.

This thesis deals with the lattice-Boltzmann method (LBM) in combination with other methods to solve thermal flow problems. The three primary, current approaches for thermal lattice-Boltzmann method (TLBM) will be introduced. The three approaches are the multispeed approach by McNamara and Alder , the passive scalar approach by Shan, and the thermal distribution model proposed by He et al. Shi et al. simplified the thermal distribution model for incompressible thermal flows. The model proposed by Shi et al. was simulated and compared to a hybrid LBM and energy equation model proposed by Khiabani et al. The thermal lattice-Boltzmann method will be compared to the temperature fields generated by the energy equation of the hybrid method. To determine which method is better suited from computer simulations the two will be compared for computational demands, and the speed of both convergence and computation.

Numerous efforts have contributed to the study of phase-change problems for over a century|both analytical and numerical. Among those numerical approximations applied to solve phase-transition problems, phase-field models attract more and more attention because they not only capture two important effects, surface tension and supercooling, but also enable explicitly labeling the solid and liquid phases and the position of the interface. In the research of this dissertation, a phase-field model has been employed to simulate 2-D dendrite growth of pure nickel without a flow, and 2-D ice crystal growth in a high-Reynolds-number lid-driven-cavity flow. In order to obtain the details of ice crystal structures as well as the flow field behavior during freezing for the latter simulation, it is necessary to solve the phase-field model without convection and the equations of motion on two different scales. To accomplish this, a heterogeneous multiscale method is implemented for the phase-field model with convection such that the phase-field model is simulated on a microscopic scale and the equations of motion are solved on a macroscopic scale. Simulations of 2-D dendrite growth of pure nickel provide the validation of the phase-field model and the study of dendrite growth under different conditions, e.g., degree of supercooling, interface thickness, kinetic coefficient, and shape of the initial seed. In addition, simulations of freezing in a lid-driven-cavity flow indicate that the flow field has great effect on the small-scale dendrite structure and the flow eld behavior on the large scale is altered by freezing inside it.

Dans ce travail, nous étudions l'écoulement au sein d'un canal symétrique avec une expansion et une contraction soudaines. Cette configuration peut être considérée comme constituée de deux cavités face à face, deux cavités en miroir, que nous dénommons "double cavité". Le sujet est traité expérimentalement, numériquement et analytiquement, en faisant varier la vitesse d'entrée et de la distance entre cavités, mais en restant à des nombres de Reynolds modérés. L'accent est mis sur l'interaction entre les deux couches de cisaillement et sur le mécanisme de rétroaction intracavitaire dans la limite des écoulements incompressibles. Expérimentalement, on mesure la vitesse par Vélocimétrie par Images de Particules non résolue en temps (PIV 2D2C) dans un plan longitudinal permettent de quanti er le champ de vitesse en moyenne temporelle. Par ailleurs, des mesures par Vélocimétrie Laser à effet Doppler (LDV) et des mesures résolues en temps par PIV 2D2C permettent d'accéder à la composition spectrale des fluctuations de vitesse dans la direction de l'écoulement. L'écoulement est caractérisé à partir des séries temporelles, enregistrées dans les couches de cisaillement d'une des deux cavités, pour une large gamme de vitesses d'entrée et des distances entre cavités. Des simulations numériques directes 2D et 3D, permettent d'étudier le mécanisme hydrodynamique de rétroaction intracavitaire, à partir des champs de vitesse complet. Le champ de vorticité issu des simulations numériques 2D montre l'importance de la rotation d'ensemble au sein de la cavité qui transporte les injections de vorticité induites par les oscillations de la couche de mélange conduisant à une structure de type "carrousel" elle-même à l'origine du mécanisme de rétroaction responsable des oscillations auto-entretenues de la couche de cisaillement. La quanti cation des temps caractéristiques de cette rotation permet d'identifier le régime dans lequel se trouve l'écoulement. Une analyse de stabilité en temps seul, ainsi qu'en temps et espace est réalisée pour des écoulements non visqueux, en prenant un écoulement de base unidimensionnel pour chacun des cas: cavité simple ou double. Pour prendre en compte l'extension finie du système, dans le cas de l'analyse de stabilité linéaire spatio-temporel, on ajoute la condition dite de Kulikowskii, qui permet de prendre en compte la réflexion des ondes d'instabilité hydrodynamique aux bornes du domaine de la cavité. Ce mécanisme de rétroaction produit un ensemble discret de fréquences non-harmoniques, dont certaines correspondent effectivement aux données expérimentales. / This work investigates the flow in a symmetric channel with a sudden expansion and contraction, creating two facing cavities, a so called double cavity. Double cavity flow at moderate Reynolds numbers is studied experimentally, numerically and analytically, as the inflow velocity and the distance between the cavities are varied. The focus is put on the interaction of the two shear layers and the intracavitary hydrodynamic feedback in the incompressible limit.Experimentally, standard 2D2C particle image velocimetry (PIV) measurements in a given spanwise plane provide information on the instantaneous and mean velocity flow fields. Laser Doppler velocimetry and time resolved 2D2C PIV measurements reveal the richness of the streamwise fluctuating velocity spectra. The flow is characterized based on times series, recorded in one of the cavity's shear layers, for a wide range of inflow velocities and cavity distances.Two dimensional and three dimensional direct numerical simulations, which give easy access to the entire flow field, are used to study the intracavitary hydrodynamic feedback mechanism. Vorticity fields, obtained from 2D numerical simulations, show the importance of the recirculating intracavitary back flow. Vorticity packages, injected by the oscillating and impinging shear layer at the downstream cavity edge, are advected upstream in the recirculation region, creating a ``carousel-like'' pattern. The interaction of this vortex carousel with the oscillating shear layer is found to be responsible for the self-sustained oscillations observed experimentally in single and double cavity flow. The quantification of three characteristic time scales of the rotation allows to identify in which regime the flow resides.Temporal and spatio-temporal inviscid linear stability analyses are applied to a one dimensional base flow of single and double cavity flows. To account for the finite extent of the system, the spatio-temporal linear stability analysis is conditioned by a so called Kulikowskii condition, which allows the reflection of hydrodynamic instability waves within the cavity domain. This feedback mechanism yields a set of discrete, non-harmonic frequencies, some of which compare well with experimental results.

No description available.

No description available.

In this work, we develop a method for solving viscoelastic fluid flows using the Navier-Stokes equations coupled with the Oldroyd-B model. We solve the Navier-Stokes equations in skew-symmetric form using the mixed finite element method, and we solve the Oldroyd-B model using the discontinuous Galerkin method. The Crank-Nicolson scheme is used for the temporal discretization of the Navier-Stokes equations in order to achieve a second-order accuracy in time, while the optimal third-order total-variation diminishing Runge-Kutta scheme is used for the temporal discretization of the Oldroyd-B equation. The overall accuracy in time is therefore limited to second-order due to the Crank-Nicolson scheme; however, a third-order Runge-Kutta scheme is implemented for greater stability over lower order Runge-Kutta schemes. We test our numerical method using the 2D cavity flow benchmark problem and compare results generated with those found in literature while discussing the influence of mesh refinement on suppressing oscillations in the polymer stress. / Master of Science / Viscoelastic fluids are a type of non-Newtonian fluid of great importance to the study of fluid flows. Such fluids exhibit both viscous and elastic behaviors. We develop a numerical method to solve the partial differential equations governing viscoelastic fluid flows using various finite element methods. Our method is then validated using previous numerical results in literature.

A two-dimensional (2-D), mathematical model is adopted to investigate the development of circulation patterns for compressible, laminar, and shear driven flow inside a rectangular cavity. The bottom of the cavity is free to move at a specified speed and the aspect ratio of the cavity is changed from 1.0 to 1.5. The vertical sides and the bottom of the cavity are assumed insulated. The cavity is filled with a compressible fluid with Prandtl number, Pr =1. The governing equations are solved numerically using the commercial Computational Fluid Dynamics (CFD) package ANSYS FLUENT 2015 and compared with the results for the primitive variables of the problem obtained using in house CFD code based on Coupled Modified Strongly Implicit Procedure (CMSIP). The simulations are carried out for the unsteady, lid driven cavity flow problem with moving boundary (bottom) for different Reynolds number, Mach numbers, bottom velocities and high initial pressure and temperature.

Os mecanismos de instabilidade hidrodinâmica têm um papel importante no processo da transição do escoamento de laminar para turbulento. A análise da instabilidade hidrodinâmica em uma cavidade com tampa deslizante foi realizada através da decomposição em modos globais (biglobal) para avaliar o efeito da compressibilidade neste fenômeno. O escoamento base foi obtido através de simulação numérica direta (DNS). Para tal, foi desenvolvido um código DNS compressível com discretização espacial por diferenças finitas compactas de alta resolução espectral e capacidade de processamento paralelo, com um método de decomposição de domínio que mantém a precisão das diferenças finitas compactas. O escoamento base é usado para montar o problema de autovalor oriundo das equações de Navier-Stokes linearizadas para a perturbação, discretizadas por diferenças finitas explícitas. O uso de diferenças finitas em conjunto com a implementação em matrizes esparsas reduz sensivelmente o uso de memória. Através do algoritmo de Arnoldi, a ordem do problema de autovalor é reduzida e os autovalores de interesse são recuperados. Os resultados indicam o efeito estabilizante da compressibilidade nos modos dominantes da cavidade e revelam modos inerentes ao escoamento compressível, para os quais a compressibilidade tem efeito desestabilizante. Dentre estes modos compressíveis, estão presentes modos de propagação sonora em dutos e modos relacionados à geração de som na cavidade. / Hydrodynamic instability mechanisms play an important role in laminar to turbulent transition. Hydrodynamic instability analysis of a lid-driven cavity flow was performed by global mode decomposition (biglobal) to evaluate compressibility effects on this phenomenon. The basic flow was calculated by direct numerical simulation (DNS). A compressible DNS code was developed with spectral-like compact finite difference spatial discretization. The code allows parallel processing with a domain decomposition method that preserves the compact finite difference accuracy. The basic flow is used to form the eigenvalue problem associated to the linear Navier- Stokes equations for the perturbation, which were discretized by an explicit finite difference scheme. The combination of sparse matrix techniques and finite difference discretization leads to a significant memory reduction. The order of the eigenvalue problem was reduced using the Arnoldi algorithm and the eigenvalues of interest were calculated. Results show the stabilizing effect of compressibility on the leading modes and reveal some modes intrinsic to compressible flow, for which compressibility has a destabilizing effect. Among these compressible modes, there are some related to sound propagation in ducts and to sound generation inside the cavity.