O trabalho desenvolvido nesta tese propõe a aplicação do método da expansão em funções hierárquicas elaborado por Zienkiewics e Morgan (1983), para a solução das equações de conservação da massa (continuidade), conservação da quantidade de movimento (Navier-Stokes) e conservação da energia, para fluidos compressíveis em duas dimensões e em alta velocidade. Esse método consiste no emprego do método de elementos finitos utilizando a formulação Petrov-Galerkin, conhecido como SUPG (“Streamline Upwind Petrov-Galerkin"), desenvolvido por Brooks e Hughes (1982), aplicado em conjunto com uma expansão das variáveis em funções hierárquicas. A fim de testar e validar o método numérico proposto, assim como o programa computacional elaborado, foram simulados alguns casos conhecidos da literatura. Os casos estudados foram os seguintes: teste de Continuidade; teste de convergência e estabilidade; problema do degrau de temperatura e problema do choque oblíquo, onde o objetivo desse último caso era, basicamente, verificar a captura da onda de choque pelo método numérico desenvolvido. Através dos casos estudados e em função dos resultados obtidos nas simulações realizadas, conclui-se que o objetivo desse trabalho foi alcançado de maneira satisfatória, pois os resultados obtidos com o método desenvolvido nesse trabalho foram qualitativamente e quantitativamente bons, quando comparados com os resultados teóricos. / The Thesis develops a new application for the Hierarchical Function Expansion Method, proposed by Zienkiewics and Morgan (1983), for the solution of the Navier-Stokes equations for compressible fluids in two dimensions and in high velocity. This method is based on the finite elements method using the Petrov-Galerkin formulation, know as, SUPG (Streamline Upwind Petrov-Galerkin) developed by Brooks and Hughes (1982), and applied in conjunction with the expansion of the variables into hierarchical functions. To test and validate the numerical method proposed as well as the computational program developed some cases whose theoretical solution are known simulated. These cases are the following: continuity test; stability and convergence test; temperature step problem; and several oblique shocks. The objective of the last cases is basically to verify the capture of the shock wave by the method developed. The results obtained in the simulations of the cases performed with the proposed method were good both qualitatively and quantitatively when compared with the teorethical solutions. This allows us to conclude that the objective of this Thesis was satisfactorily reached.

No description available.

O trabalho desenvolvido nesta tese propõe a aplicação do método da expansão em funções hierárquicas elaborado por Zienkiewics e Morgan (1983), para a solução das equações de conservação da massa (continuidade), conservação da quantidade de movimento (Navier-Stokes) e conservação da energia, para fluidos compressíveis em duas dimensões e em alta velocidade. Esse método consiste no emprego do método de elementos finitos utilizando a formulação Petrov-Galerkin, conhecido como SUPG (“Streamline Upwind Petrov-Galerkin”), desenvolvido por Brooks e Hughes (1982), aplicado em conjunto com uma expansão das variáveis em funções hierárquicas. A fim de testar e validar o método numérico proposto, assim como o programa computacional elaborado, foram simulados alguns casos conhecidos da literatura. Os casos estudados foram os seguintes: teste de Continuidade; teste de convergência e estabilidade; problema do degrau de temperatura e problema do choque oblíquo, onde o objetivo desse último caso era, basicamente, verificar a captura da onda de choque pelo método numérico desenvolvido. Através dos casos estudados e em função dos resultados obtidos nas simulações realizadas, conclui-se que o objetivo desse trabalho foi alcançado de maneira satisfatória, pois os resultados obtidos com o método desenvolvido nesse trabalho foram qualitativamente e quantitativamente bons, quando comparados com os resultados teóricos. / The Thesis develops a new application for the Hierarchical Function Expansion Method, proposed by Zienkiewics and Morgan (1983), for the solution of the Navier-Stokes equations for compressible fluids in two dimensions and in high velocity. This method is based on the finite elements method using the Petrov-Galerkin formulation, know as, SUPG (Streamline Upwind Petrov-Galerkin) developed by Brooks and Hughes (1982), and applied in conjunction with the expansion of the variables into hierarchical functions. To test and validate the numerical method proposed as well as the computational program developed some cases whose theoretical solution are known simulated. These cases are the following: continuity test; stability and convergence test; temperature step problem; and several oblique shocks. The objective of the last cases is basically to verify the capture of the shock wave by the method developed. The results obtained in the simulations of the cases performed with the proposed method were good both qualitatively and quantitatively when compared with the teorethical solutions. This allows us to conclude that the objective of this Thesis was satisfactorily reached.

Ce travail de thèse porte sur l'étude d'un système d'interaction fluide-structure. Nous en traitons de nombreux aspects allant de sa modélisation jusqu'à l'étude de sa stabilisation et de sa simulation numérique. Le premier chapitre du manuscrit aborde la modélisation du système ainsi que l'existence de solutions fortes en temps petits. Le fluide est représenté par les équations de Navier-Stokes incompressibles. La structure est déformable et dépend d'un nombre fini de paramètres. Nous obtenons ses équations en appliquant un principe des travaux virtuels. Le système d'équations final est non linéaire. Nous prouvons l'existence locale d'une solution à ce système, dans un premier temps sur le système linéarisé autour de l'état nul. Puis, nous prouvons l'existence de solutions en temps petits au système non linéaire grâce à un argument de point fixe. Le deuxième chapitre traite de la stabilisation par feedback autour d'un état stationnaire non nul du système présenté dans le Chapitre 1. L'opérateur de feedback est déterminé à partir de l'analyse du problème linéarisé autour de l'état stationnaire et de la résolution d'une équation de Riccati. Le résultat de stabilisation portant sur le système non linéaire requiert des données petites et est obtenu par un argument de point fixe. Le troisième chapitre se concentre sur les aspects numériques de ce problème. La construction de l'opérateur de feedback correspond à la version discrétisée de celle proposée dans le Chapitre 2. Le système fluide-structure est simulé en utilisant une méthode de domaines fictifs. / This PhD thesis deals with the study of a fluid-structure interaction system. We are interested in several aspects such as modelling, stabilization and numerical simulation. In the first chapter of the manuscript, we show the modelling of the system and prove the existence of strong solutions in small times. The fluid is modelled by the incompressible Navier- Stokes equations. The structure is deformable and depends on a finite number of parameters. The equations are obtained with a virtual work principle. The final system of equations is nonlinear. We prove local existence of a solution to this system, first on the linearized system. Then, existence of solutions in small times to the full nonlinear system is obtained with a fixed point argument. In the second chapter, we prove feedback stabilization of the problem around a non-null stationary state. The feedback operator is computed with the solution to a Riccati equation obtained by the analysis of the linearized problem around the stationary state. The stabilization result holds on the full nonlinear system and requires small data. It is proven by a fixed point argument. In the third chapter, we focus on the numerical aspects of the problem. The feedback operator used corresponds to a discretization of the feedback operator of Chapter 2. The solution to the full nonlinear system is computed by the use of a fictitious domain method.

Lors de sa rentrée dans l’atmosphère d’une planète, un engin spatial subit un échauffement important dû aux frottements des gaz atmosphériques sur la paroi. Cette élévation de température conduit à une dégradation physico-chimique du bouclier thermique de l’objet constitué de matériaux composites. Un composite est constitué de divers matériaux qui s’ablatent différemment. Dans cette thèse, nous nous intéressons essentiellement à la fusion d’un matériau durant sa phase de rentrée atmosphérique. Nous sommes donc en présence de trois phases : solide, liquide et gaz. Pour simuler ce phénomène, des méthodes numériques robustes ont été mises au point pour calculer l’écoulement diphasique compressible autour de l’objet. Le couplage entre le solide et l’écoulement fluide a aussi été étudié. Les méthodes numériques développées durant cette thèse sont basées sur une approche volumes finis. Une stratégie de décomposition d’opérateurs est utilisée pour résoudre le modèle diphasique à cinq équations avec les termes de dissipation modélisant l’écoulement fluide. L’idée principale de cette décomposition d’opérateurs est de séparer les phénomènes acoustiques et dissipatifs des phénomènes de transport. Un traitement implicite de l’étape acoustique est réalisé tandis que l’étape de transport est résolue explicitement. Le schéma semi-implicite global est alors très robuste, conservatif et préserve les discontinuités de contact. Les conditions d’interface entre les domaines fluide et solide sont déduites des bilans de masse et d’énergie à la paroi. Le front de fusion est suivi explicitement grâce à une formulation ALE des équations. La robustesse de l’approche et l’apport de la formulation semi-implicite sont finalement démontrés grâce à des expériences numériques mono et bidimensionnelles sur maillages curvilignes mobiles. / During atmospheric re-entry phase, a spacecraft undergoes a sudden increase of the temperature due to the friction of atmospheric gases. This rise drives to a physical-chemical degradation of the thermal protective system of the object made of composite material. A composite is made of several materials with ablates differently. In this thesis, we mainly focus on the melting of an object during its re-entry phase. Therefore there are three phases: solid, liquid and gas phases. In order to simulate this phenomenon, robust numerical methods have been developed to compute a compressible multiphase flow. The coupling strategy between the solid and the fluid have also been studied. Solvers developed in the present work are based on Finite Volume Method. A splitting strategy is used to compute compressible two-phase flows using the five-equation model with viscous and heat conduction effects. The main idea of the splitting is to separate the acoustic and dissipative phenomena from the transport one. An implicit treatment of the acoustic step is performed while the transport step is solved explicitly. The overall scheme resulting from this splitting operator strategy is very robust, conservative, and preserves contact discontinuities. The boundary interface condition between the solid and the multiphase flow is enforced by mass and energy balances at the wall. The melting front is tracked explicitly using an ALE formulation of the equations. The robustness of the approach and the interest of the semi-implicit formulation are demonstrated through numerical simulations in one and two dimensions on moving curvilinear grids.

Dans cette thèse, nous étudions les propriétés des écoulements de fluides qui interagissent avec un corps rigide ou avec une source et un puits. Dans le cas d'un fluide visqueux incompressible qui satisfait les équations de Navier Stokes dans un domaine borné 2D, les solutions faibles de Leray-Hopf sont bien comprises. L'existence et l'unicité sont prouvées. De plus, les solutions sont continues en temps `a valeurs dans L 2 (Omega) et satisfont l’égalité d'énergie classique. Plus récemment, le problème d'un corps rigide en mouvement dans un fluide visqueux incompressible modélisé par les équations de Navier-Stokes couplées aux lois de Newton qui décrivent le mouvement du solide a également été abordé dans le cas où des conditions aux limites sans glissement ont été prescrites. Des résultats analogues concernant les solutions de Leray-Hopf ont également été démontrés dans ce contexte. Dans ce manuscrit, nous étudions le cas de conditions aux limites de Navier-Slip. Dans ce cadre, le résultat d'existence pour le système couplé a été prouvée par G'erard-Varet et Hillairet en 2014. Ici, nous montrons que les solutions sont continues en temps, qu'elles satisfont l’égalité d'énergie et qu’elles sont uniques. De plus, nous montrons un résultat d'existence des solutions faibles dans le cas d'un fluide incompressible visqueux auquel s'ajoute un corps rigide dans le cas où la vitesse du fluide a une partie orthonormale d'énergie infinie.Pour un fluide incompressible non visqueux modélisé par les équations d'Euler dans un domaine borné 2D, le cas où le fluide est autorisé à entrer et à sortir de la frontière a été traité par Judovic qui a introduit certaines conditions limites consistant à prescrire la composante normale de la vitesse et de la vorticité entrante. Dans ce manuscrit, nous considérons un domaine borné qui possède deux trous. L'un d'eux est une source, ce qui signifie que le fluide est autorisé à entrer dans le domaine et l'autre est un puits où le fluide peut sortir. En particulier, nous établissons les équations limites vérifiées par le fluide lorsque la source et le puits se contractent en deux points différents. Le système limite est caractérisé par un point source/puits et un point vortex en chacun des deux points où les trous se sont contractés. / In this thesis, we investigate properties of incompressible flows that interact with a rigid body or a source and a sink. In the case of an incompressible viscous fluid that satisfies the Navier Stokes equations in a 2D bounded domain well-posedness of Leray-Hopf weak solutions is well-understood. Existence and uniqueness are proved. Moreover solutions are continuous in time with values in L 2 (Omega) and they satisfy the energy equality. Recently the problem of a rigid body moving in a viscous incompressible fluid modeled by the Navier-Stokes equations coupled with the Newton laws that prescribe the motion of the solid, was also tackled in the case where the no-slip boundary conditions were imposed. And the correspondent well-posedness result for Leray-Hopf type weak solutions was proved. In this manuscript we consider the case of the Navier-slip boundary conditions. In this setting, the existence result for the coupled system was proved by G'erard-Varet and Hillairet in 2014. Here, we prove that solutions are continuous in time, that they satisfy the energy equality and that they are unique. Moreover we show an existence result for weak solutions of a viscous incompressible fluid plus rigid body system in the case where the fluid velocity has an orthoradial part of infinite energy.For an inviscid incompressible fluid modelled by the Euler equations in a 2D bounded domain, the case where the fluid is allowed to enter and to exit from the boundary was tackled by Judovic who introduced some conditions which consist in prescribing the normal component of the velocity and the entering vorticity. In this manuscript we consider a bounded domain with two holes, one of them is a source which means that the fluid is allowed to enter in the domain and the other is a sink from where the fluid can exit. In particular we find the limiting equations satisfied by the fluid when the source and the sink shrink to two different points. The limiting system is characterized by a point source/sink and a point vortex in each of the two points where the holes shrunk.

La résolution des équations de Navier-Stokes linéarisées compressibles est utilisée pour 2 types de problèmes : 1. Pour la résolution de problèmes d'aéroélastcité et d'aéroacoutstique 2. Les exercices d'optimisation par méthode du gradient Les algorithmes proposés sont généralement basés sur une approche dite time-marching simplifiant le développement numérique. Dans ce doctorat nous avons développer une méthode sans intégration temporelle pour stabiliser la résolution des équations de Navier-Stokes linéarisées. Les systèmes linéaires obtenus sont résolus par une méthode itérative de type GMRes-ILU0. Les résultats numériques sont comparés à une approche AF-ADI et GMRes-time-marching pour des calculs 2D relatif à une perturbation harmonique de pression. La méthode de résolution est aussi validée pour 2 exercices d'optimisation. Une méthode de résolution pseudo-Newton-GMRes des équations de Navier-Stokes non linéaires faiblement couplée a aussi été validée dans le cas d'écoulements 2D.

Helicopter performance is greatly influenced by its drag. Pylons, fuselage, landing gear, and especially the rotor hub of a helicopter experience large separated flow regions, even under steady level flight conditions the vehicle has been designed for, contributing to the helicopter drag. Several passive and active flow control concepts have been studied for reducing helicopter drag. While passive flow control methods reduce drag, they do so at one optimized design condition. Therefore, passive drag reduction methods may not work for helicopters that operate under widely varying flight conditions. Active flow control (AFC) methods overcome this disadvantage and consequently are widely being pursued. The present investigator has studied some of these AFC methods using computational fluid dynamics (CFD) techniques and has found synthetic (or pulsed) jets as one of the more effective drag reduction devices. Two bluff bodies, representative of helicopter components, have been studied and the mechanism behind drag reduction has been analyzed. It was found that the increase in momentum due to the jet, and a resultant reduction in the separated flow region, is the main reason for drag reduction in these configurations. In comparison with steady jets, synthetic jets were found to use less power for a greater drag reduction. The flow inside these synthetic jet devices is incompressible. It is computationally inefficient to use compressible flow solvers in incompressible regions. In such regions, using Lattice Boltzmann equations (LBE) is more suitable compared to solving the incompressible Navier-Stokes equations. The length scales close to the synthetic jet devices are very small. LBE may be used to better resolve these small length scale regions. However, using LBE throughout the whole domain would be computationally expensive since the grid spacing in the LBE solver has to be of the order of the mean free path. To address this need, a coupled Lattice Boltzmann-Navier-Stokes (LB-NS) methodology has been developed. The LBE solver has been successfully validated in a standalone manner for several benchmark cases. The solver has also been shown to be of second order accuracy. This LBE solver has been subsequently coupled with an existing Navier-Stokes (NS) solver. Validation of the coupled methodology has been done for analytical problems with known closed form solution. This LB-NS methodology is further used to simulate the flow past a cylinder where synthetic jet devices have been used to reduce drag. The LBE solver is used in the cavity of the synthetic jet nozzle while the NS solver is employed in the rest of the domain. The cylinder configuration was chosen to demonstrate drag reduction on helicopter hub shape geometries. Significant drag reduction is observed when synthetic jets are used, compared to the baseline no flow control case.

Hydrodynamic behavior in reactors used for water treatment, particularly in ozone contactors with serpentine flow, is known to strongly affect the process efficiency. However, exact flow characteristics inside these reactors are not well understood, as traditional approach either considers these reactors as black box or relies on less accurate Reynolds-Averaged Navier-Stokes (RANS) simulation. In order to provide a deep understanding of the hydrodynamics and solute transport phenomena in these reactors, high resolution numerical studies using the Large Eddy Simulation (LES) method are performed. The reactor geometries investigated in this research are Constant Baffle Spacing Multi-Chamber (CBSMC) ozone contactors and a Variable Baffle Spacing ozone contactor Model (VBSM). The LES results in two multi-chamber ozone contactors (CBSMC -Normal-Width and -Half-Width) suggest that the flow through these reactors is characterized by the presence of extensive short-circuiting and large internal recirculation. The results also suggest that the flow is highly three dimensional with a pair of symmetric counter-rotating secondary vortices. LES studies based on VBSM, the baffle spacing of which varies between 0.5 times to 5 times the size of the base chamber; suggest that the width of the recirculation zone grows at about the same rate as the baffle spacing. Instantaneous turbulent eddies are prevalent in the chamber and increase turbulent mixing. The elevated levels of turbulence are found in the short-circuiting flow path. The tracer is dispersed along the short-circuiting path and strongly into the recirculation zone due to turbulent diffusion. Baffle spacing greater than the entrance gate height, but also smaller baffle spacing, worsens the disinfection efficiency. Finally, the turbulent Schmidt number of RANS simulation was investigated by employing the previously validated LES simulation. Due to the presence of very strong turbulent diffusion in the reactors, the turbulent Schmidt number is found out to be much less than the values commonly used, and is also specific to the baffle spacing.

In this thesis we describe the development of Reynolds- averaged Navier Stokes code for the flow past two- dimensional configuration. Particularly, emphasis has been laid on the study of low Reynolds number airfoil aerodynamics. The thesis consists of five chapters covering the back ground history, problem formulation, method of solution and discussion of the results and conclusion. Chapter I deals with a detailed background history of low Reynolds number aerodynamics, problem associated with it, state of the art, its importance in practical applications in aircraft industries. Chapter II describes the mathematical model of the flow physics and various levels of approximations. Also it gives an account of complexity of the equations at low Reynolds number regarding flow separation, transition and reattachment. Chapter III describes method of solution, numerical algorithm developed, description of various upwind schemes, grid system, finite volume discrieti-zation of the governing equations described in Chapter II. Chapter IV describes the application of the newly developed Navier Stokes code for the test cases from GAMM Workshop proceedings. Also it describes validation of the code for Euler solutions, Blasius solution for the flow past flat plate and compressible Navier Stokes solution for the flow past NACA 0012 Airfoil at low Reynolds number. Chapter V describes the application of the Navier Stokes code for the more test cases of current practical interest . In this chapter laminar separation bubble characteristics are investigated in detail regarding formation, growth and shedding in an unsteady environment. Finally the conclusion is drawn regarding the robustness of the newly developed code in predicting the airfoil aerodynamic characteristics at low Reynolds number both in steady and unsteady environment. Lastly, suggestion for future work has been highlighted.