This thesis deals with flow measuring techniques applied on steady and pulsatingflows. Specifically, it is focused on gas flows where density changes canbe significant, i.e. compressible flows. In such flows only the mass flow ratehas a significance and not the volume flow rate since the latter depends onthe pressure. The motivation for the present study is found in the use of flowmeters for various purposes in the gas exchange system for internal combustionengines. Applications can be found for instance regarding measurements of airflow to the engine, or measurements of the amount of exhaust gas recirculation.However the scope of thesis is wider than this, since the thesis aims toinvestigate the response of flow meters to pulsating flows. The study is mainlyexperimental, but it also includes an introduction and discussion of several inindustry, common flow measuring techniques.The flow meters were studied using a newly developed flow rig, designedfor measurement of steady and pulsating air flow of mass flow rates and pulsefrequencies typically found in the gas exchange system of cars and smallertrucks. Flow rates are up to about 200 g/s and pulsation frequencies from 0 Hz(i.e. steady flow) up to 80 Hz. The study included the following flow meters:hot-film mass flow meter, venturi flowmeter, Pitot tube, vortex flowmeter andturbine flowmeter. The performance of these meters were evaluated at bothsteady and pulsating conditions. Furthermore, the flow under both steady andpulsating conditions were characterized by means of a resistance-wire basedmass flow meter, with the ability to perform time resolved measurements ofboth the mass flux ρu, and the stagnation temperature T0.Experiments shows that, for certain flow meters, a quasi-steady assumptionis fairly well justified at pulsating flow conditions. This means that thefundamental equations describing the steady flow, for each instant of time,is applicable also in the pulsating flow. In the set-up, back-flow occurred atcertain pulse frequencies, which can result in highly inaccurate output fromcertain flow meters, depending on the measurement principle. For the purposeof finding means to determine when back flow prevails, LDV measurementswere also carried out. These measurements were compared with measurementsusing a vortex flow meter together with a new signal processing technique basedon wavelet analysis. The comparison showed that this technique may have apotential to measure pulsating flow rates accurately.Descriptors: Flow measuring, compressible flow, steady flow, pulsating flow,hot-wire anemometry, cold-wire anemometry. / QC 20101208

The compressible turbulent boundary layer developing over atwo-dimensional bump which leads to a supersonic pocket with aterminating shock wave has been studied. The measurements havebeen made with hot-wire anemometry and Particle ImageVelocimetry (PIV). A method to calibrate hot-wire probes in compressible ow hasbeen developed which take into account not only the ow velocitybut also the inuence of the Mach number, stagnation temperatureand uid density. The calibration unit consists of a small jetow facility, where the temperature can be varied. The hot wiresare calibrated in the potential core of the free jet. The jetemanates in a container where the static pressure can becontrolled, and thereby the gas density. The calibration methodwas verfied in the at plate zero pressure gradient turbulentboundary layer in front of the bump at three different Machnumbers, namely 0.3, 0.5 and 0.7. The profiles were alsomeasured at different static pressures in order to see theinuence of varying density. Good agreement between the profilesmeasured at different pressures, as well as with the standardlogarithmic profile was obtained. The PIV measurements of the boundary layer ow in front ofthe 2D bump showed good agreement with the velocity profilesmeasured with hotwire anemometry. The shock wave boundary layerinteraction was investigated for an inlet Mach number of 0.69.A lambda shock wave was seen on the downstream side of thebump. The velocity on both sides of the shock wave as measuredwith the PIV was in good agreement with theory. The shock wavewas found to cause boundary layer separation, which was seen asa rapid growth of the boundary layer thickness downstream theshock. However, no back ow was seen in the PIV-data, probablybecause the seeding did not give enough particles in theseparated region. The PIV data also showed that the shock wavewas oscillating, i.e. it was moving approximately 5 mm back andforth. This distance corresponds to about five boundary layerthicknesses in terms of the boundary layer upstream theshock. <b>Descriptors:</b>Fluid mechanics, compressible ow,turbulence, boundary layer, hot-wire anemometry, PIV, shockwave boundary layer interaction, shape factor. / NR 20140805

No description available.

No description available.

During the launch sequence of the United Launch Alliance Delta IV launch vehicle, large amounts of pure hydrogen are introduced into the launch table and ignited by Radial-Outward-Firing-Igniters (ROFIs). This ignition results in a significant flame, or plume, that rises upwards out of the launch table due to buoyancy. The presence of the plume causes increased and unwanted heat loads on the surface of the vehicle. A proposed solution is to add a series of fans and structures to the existing launch table configuration that are designed to inject ambient air in the immediate vicinity of the launch vehicle's nozzles to suppress the plume rise. In addition to the air injection, secondary fan systems can be added around the launch table openings to further suppress the hydrogen plume. The proposed air injection solution is validated by computational fluid dynamics simulations that capture the combustion and compressible flow observed during the Delta IV launch sequence. A solution to the hydrogen plume problem will have direct influence on the efficiency of the launch vehicle: lower heat loads result in thinner vehicle insulation and thus allow for a larger payload mass. Current results show that air injection around the launch vehicle nozzles and air suppression around the launch table openings significantly reduces the size of the plume around the launch vehicle prior to liftoff.

Le présent manuscrit de thèse présente l'étude de différents phénomènes dans un convertisseur d?acier, grâce au développement de deux modèles distincts. Le premier modèle décrit la cavité produite à la surface libre du bain de métal par l'impact du jet d'oxygène supersonique. Il est basé sur le découpage du domaine de calcul en deux régions. Les effets de compressibilité du gaz sont pris en compte uniquement dans la région du jet où la vitesse est élevée, alors que partout ailleurs, le gaz est considéré comme incompressible. La méthode Volume Of Fluid (VOF) est utilisée pour suivre la déformation de la surface libre du bain. Les résultats de simulations sont présentés pour des systèmes bi- et triphasés et comparés à des données expérimentales obtenues dans diverses maquettes froides. L'influence sur la taille et la forme de la cavité de différents paramètres (parmi lesquels les conditions aux limites en sortie de la lance d'injection, le schéma d'advection de la méthode VOF et le modèle de turbulence) est étudiée. Le modèle est ensuite utilisé pour simuler l'interaction entre un jet supersonique d'oxygène et la surface libre d'un bain d'acier liquide dans un convertisseur de taille pilote. Le second modèle se focalise sur l'écoulement du gaz, le transfert de chaleur et la réaction de postcombustion dans la phase gazeuse au-dessus du bain de métal. Il utilise l'algorithme Simple Chemical Reaction Scheme pour décrire le transport des espèces chimiques, et prend en compte l'absorption d'oxygène dans le bain et les transferts thermiques par rayonnement. Les prédictions numériques sont en assez bon accord avec les mesures recueillies dans une expérience de laboratoire et dans un four à l'échelle pilote / The present thesis treats different phenomena taking place in a steelmaking converter through the development of two separate models. The first model describes the cavity produced at the free surface of the metal bath by the high speed impinging oxygen jet. It is based on a zonal approach, where gas compressibility effects are taken into account only in the high velocity jet region while elsewhere the gas is treated as incompressible. The Volume Of Fluid (VOF) method is employed to follow the deformation of the bath free surface. Calculations are presented for two- and three-phase systems and compared against experimental data obtained in various cold model experiments. The influence on the size and shape of the cavity of various parameters (including the jet inlet boundary conditions, the VOF advection scheme and the turbulence modelling) is studied. Next, the model is used to simulate the interaction of a supersonic oxygen jet with the surface of a liquid steel bath in a pilot-scale converter. The second model concentrates on fluid flow, heat transfer and the post-combustion reaction in the gas phase above the metal bath. It uses the Simple Chemical Reaction Scheme approach to describe the transport of the chemical species and takes into account the consumption of oxygen by the bath and thermal radiative transfer. The numerical predictions are in reasonable agreement with measurements collected in a laboratory experiment and in a pilot-scale furnace

La simulation numérique est largement utilisée pour évaluer les performances aérodynamiques des aéronefs ainsi qu'en optimisation de forme. Ainsi l'objectif de ces simulations est souvent le calcul de fonctions aérodynamiques. L'objet de cette thèse est d'étudier des méthodes d'adaptation de maillages basées sur la dérivée totale de ces fonctions par rapport aux coordonnées du maillage (notée dJ/dX). Celle-ci pouvant être calculée par la méthode adjointe discrète. La première partie de cette étude concerne l'application de méthodes d'adaptation de maillages appliquées à des écoulements de fluides parfaits. Le senseur qui détecte les zones de maillage à raffiner s'appuie sur la norme de cette dérivée pour adapter des maillages pour le calcul d'une fonction J. La seconde partie du travail est la construction et l'étude de critères plus fiables basés sur dJ/dX pour d'une part adapter des maillages et d'autre part estimer si un maillage est bien adapté ou non pour le calcul de la fonction J. De plus une méthode de remaillage plus efficace basée sur une EDP elliptique est aussi présentée. Cette nouvelle méthode est appliquée pour des écoulements bidimensionnels de fluides parfaits ainsi que pour un écoulement décrit par les équations RANS. La dernière partie de l'étude est consacrée à l'application de la méthode proposée à des cas tridimensionnels d'écoulement RANS sur des géométries d'intérêt industriel. / Numerical simulation is widely used for the assessment of aircraft aerodynamic performances and shape optimizations. Hence the objective of these simulations is often to compute aerodynamic outputs. The purpose of this thesis is to study mesh adaptation methods based on the total derivative of the outputs with respect to mesh coordinates (denoted dJ/dX). This derivative can be computed using the discrete adjoint method. The first part of this study is about the application of mesh adaptation methods applied for Eulerian flows. The mesh locations to refine are detected using a sensor based on the norm of the derivative dJ/dX. This study confirmed that this derivative is relevant in order to adapt a mesh for the computation of the output J. The second part of this work is the construction and the study of reliable criteria based on dJ/dX for both mesh adaptation and the quality assessment of a given mesh for the computation of the output J. Moreover a more efficient remeshing method based on an elliptic PDE is presented too. This new method is applied for both two-dimensional Eulerian flows and a flow described by the RANS equations. The last part of the study is devoted to the application of the proposed method to three-dimensional RANS flows on geometries of industrial interest.

O som gerado por escoamentos, também conhecido como aeroacústica, tem se tornado cada vez mais importante em áreas industriais diversas desde aviação comercial até aparelhos eletrodomésticos, afetando diretamente os requisitos necessários para o desenvolvimento de novos produtos. Um caso particular é o ruído gerado por válvulas de compressores herméticos de refrigeração, sendo o compressor a principal fonte de ruído em refrigeradores domésticos. O presente trabalho tem por objetivo iniciar o desenvolvimento de uma ferramenta confiável de simulação capaz de auxiliar engenheiros na predição de problemas de aeroacústica, especialmente um que possa no futuro ser utilizado para estudar o ruído gerado pelo escoamento em válvulas de compressores herméticos. Para isso, foi desenvolvido um código para simulação numérica direta de aeroacústica. Utilizou-se processamento paralelo com decomposição de domínio para usar Simulação Numérica Direta em um tempo factível; esquemas de discretização espaciais e temporais de alta ordem para minimizar ao máximo os fenômenos de dissipação e dispersão do escoamento e das ondas acústicas e uma série de tratamentos no domínio como filtragem e estiramento da malha como também condições de contorno características com o intuito de obter uma solução adequada para estudo de aeroacústica. Assim, são apresentadas todas as etapas desenvolvidas no equacionamento, implementação e verificação. A verificação foi realizada segundo um processo matemático formal (Método das Soluções Manufaturadas) com o qual obteve-se que a ordem de precisão dos cálculos era a mesma da ordem formal dos esquemas de discretização utilizados para todas as variáveis. Também obteve-se a mesma concordância para análise do divergente da velocidade, verificando o código para simulação numérica direta de aeroacústica. Posteriormente, foram realizadas simulações de escoamentos compressíveis cisalhantes e seus resultados comparados com dados apresentados em literatura. Também foram calculadas as taxas de amplificação de perturbações e comparadas com a Teoria de Estabilidade Linear. Novamente, foram obtidos resultados satisfatórios nessas etapas, mostrando que a implementação do código DNS está verificada. / Sound generated by flow, known as aeroacoustics, is becoming more important in several industrial areas from commercial aircraft to household appliances, affecting directly the requirements to the development of new products. A particular case is the noise generated by valves of refrigeration hermetic compressors, being the compressor the main noise source in household refrigerators. This work has the goal of initiate the development of a reliable tool able to help engineers to predict aeroacoustics problems, specially one that can be used in the future to study the noise generated by the flow in valves of hermetic compressors. To do so, it was developed a numerical code to perform direct numerical simulation of aeroacoustics. It was used parallel processing with domain decomposition to use Direct Numerical Simulation in a feasible time; high order temporal and spatial discretization schemes to minimize the most the dispersion and dissipation phenomena of the flow field and of the acoustics waves and a series of treatments in the domain as filtering and mesh stretching as well as characteristics boundary conditions aiming a proper solution to study aeroacoustics. Thus, here all the steps developed in the formulation, implementation and verification are presented . The verification was done according to a formal mathematical procedure (Method of Manufactured Solutions) with which was found that the precision order of the calculations was the same of the formal order of the used discretization schemes for all variables. The same agreement was also obtained to the analysis of the divergence of the velocity, verifying the code to direct numerical simulation of aeroacoustics. Posteriorly, it were simulated shear compressible flows and the results were compared to literature data. Also, it were calculated the amplification rates of the disturbances and compared to Linear Stability Theory. Once more, it was obtained satisfactory results in these steps, showing that the implementation of the DNS code is verified.

A investigação dos mecanismos de instabilidade pode ter um papel importante no entendimento do processo laminar para turbulento de um escoamento. Análise de instabilidade de uma camada limite de uma linha de estagnação compressível foi realizada no contexto de teoria linear BiGlobal. O estudo dos mecanismos de instabilidade deste escoamento pode proporcionar uma visão útil no desenho aerodinâmico das asas. Um novo procedimento foi desenvolvido e implementado computacionalmente de maneira sequencial e paralela para o estudo de instabilidade BiGlobal. O mesmo baseia-se em formar a matriz esparsa associada ao problema discretizado por dois métodos: pontos de colocação de Chebyshev-Gauss-Lobatto e diferenças finitas, além das combinações destes métodos. Isto permitiu o uso de bibliotecas computacionais eficientes para resolver o sistema linear associado ao problema de autovalor utilizando o algoritmo de Arnoldi. O desempenho do método numérico e código computacional proposto são analisados do ponto de vista do uso de métodos de ordenação dos elementos da matriz, coeficientes de preenchimento, memória e tempo computacional a fim de determinar a solução mais eficiente para um problema físico geral com técnicas de matrizes esparsas. Um estudo paramétrico da instabilidade da camada limite de uma linha de estagnação foi realizado incluindo o estudo dos efeitos de compressibilidade. O excelente desempenho código computacional permitiu obter as curvas neutras e seus respectivos valores críticos para a faixa de número de Mach 0 \'< OU =\' Ma \'< OU =\' 1. Os resultados confirmam a teoria assintótica apresentada por (THEOFILIS; FEDOROV; COLLIS, 2004) e mostram que o incremento do número de Mach reduz o numero de Reynolds crítico e a faixa instável do número de ondas. / Investigation of linear instability mechanisms is essential for understanding the process of transition from laminar to turbulent flow. An algorithm for the numerical solution of the compressible BiGlobal eigenvalue problem is developed. This algorithm exploits the sparsity of the matrices resulting from the spatial discretization of the enigenvalue problem in order to improve the performance in terms of both memory and CPU time over previous dense algebra solutions. Spectral collocation and finite differences spatial discretization methods are implemented, and a performance study is carried out in order to determine the best practice for the efficient solution of a general physical problem with sparse matrix techniques. A combination of spectral collocation and finite differences can further improve the performance. The code developed is then applied in order to revisit and complete the parametric analyses on global instability of the compressible swept Hiemenz flow initiated in (THEOFILIS; FEDOROV; COLLIS, 2004) and obtain neutral curves of this flow as a function of the Mach number in the 0 \'< OU =\' Ma \'< OU =\' 1 range. The present numerical results fully confirm the asymptotic theory results presented in (THEOFILIS; FEDOROV; COLLIS, 2004). This work presents a complete parametric study of the instability properties of modal three dimensional disturbances in the subsonic range for the flow conguration at hand. Up to the subsonic maximum Mach number value studied, it is found that an increase in this parameter reduces the critical Reynolds number and the range of the unstable spanwise wavenumbers.

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.