An Axisymmetric Rocket Ejector Simulation (ARES) was developed to numerically analyze various configurations of an air augmented rocket. Primary and secondary flow field visualizations are presented and performance predictions are tabulated. A parametric study on ejector geometry is obtained following a validation of the flow fields and performance values. The primary flow is calculated using a quasi-2D, irrotational Method of Characteristics and the secondary flow is found using isentropic relations. Primary calculations begin at the throat and extend through the nozzle to the location of the first Mach Disk. Combustion properties are tabulated before analysis to allow for propellant property selection. Secondary flow calculations employ the previously calculated plume boundary and ejector geometry to form an isentropic solution. Primary and secondary flow computations are iterated along the new pressure distributions established by the 1D analysis until a convergence tolerance is met. Thrust augmentation and Specific Impulse values are predicted using a control volume approach. For the validation test cases, the nozzle characteristic net is very similar to that of previous research. Plume characteristics are in good agreement but fluctuate in accuracy due to flow structure formulation. The individual unit processes utilized by the Method of Characteristics are found to vary their outputs by up to 0.025% when compared to existing sources. Rocket thrust and specific impulse are increased by up to 22% for a static system and 15% for an ejector flow at Mach 0.5. Evidence of Fabri conditions were observed in the flow visualization and graphically through the performance predictions. It was determined that the optimum ejector divergence angle for an air augmented rocket greatly depends on the stagnation pressure ratio between the primary and secondary flows.

Computational Fluid Dynamics is used to analyze the formation of under-expanded jets and to investigate the three-dimensional flow field associated with the impingement of free jets onto stationary deflectors. This investigation was performed to develop a verified modeling ability for such problems. Predictions were compared with the experimental results obtained by Donaldson and Snedeker [1]. Computational models for free and impinging jets were created according to the data provided in Ref. 1. Numerical results for each of the experiments performed in this benchmark report are presented. Three different turbulent free jets produced by a simple convergent nozzle were analyzed. These include a subsonic jet with p₁/pâ =1 and M₁=0.57, a moderately under-expanded jet with p₁/pâ =1.42 and M₁=1, and a highly under-expanded jet with p₁/pâ =3.57 and M₁=1. The reflecting shocks associated with the moderately under-expanded jet as well as the shock disk associated with the highly under-expanded jet were fully resolved. Velocity profile data predicted at locations downstream of the nozzle exit agreed very well with the experimental results. The impingement of a moderately under-expanded jet with p₁/pâ =1.42 and M₁=1 was also investigated. The interaction of the high speed jet with circular flat plates at angles of 60° and 45° relative to the center axis of the jet are presented. Wall jet velocity profiles on the surface of the flat plate are fully resolved and compare well with experimental results. The CFD solver controls and method used to obtain these results are summarized and justified. / Master of Science

Two-dimensional supersonic flows out of rocket nozzles are one of three flow types: over-, perfectly-, or under-expanded. In under-expanded flows, an expansion fan is centered at the top and bottom tip of the rocket nozzle. When the waves from the expansion wave cross through the centerline and intersect the free boundary, the waves are reflected as compression waves. For higher exit-to-ambient pressure ratios, the compression waves coalesce and eventually form a barrel shock. The purpose of this study was to use the Method of Characteristics (MOC), a mathematical procedure for solving hyperbolic partial differential equations, to simulate the formation of the barrel shock. A MOC code was developed in the Python programming language to accomplish task. Results of the MOC code compared favorably with CFD results.

L’objectif de ces travaux est d’étudier expérimentalement la dispersion d’hydrogène pour évaluer l’impact des fuites chroniques ou accidentelles qui peuvent intervenir en milieu libre ou encombré, à faible ou à forte pression, sur un véhicule fonctionnant avec une pile à combustible. Les fuites étudiées sont assimilées à des jets verticaux, turbulents, axisymétriques, à densité variable, et issus d’orifices cylindriques de 1 à 3 mm de diamètre. Un banc expérimental a été conçu pour étudier ces fuites : l’hydrogène a été remplacé par de l’hélium pour des raisons de sécurité. Il résiste à une pression de 200 bars et permet de positionner un obstacle dans le jet. La technique BOS (Background Oriented Schlieren) a été adaptée aux jets millimétriques et à la présence d’un obstacle. Un soin particulier à été apporté à la mise en place de cette technique. Les résultats obtenus sont en accord avec ceux de la littérature quand il en existe. De nouvelles lois de similitude sont proposées, plus représentatives de la physique des jets : libres subsoniques, libres sousdétendus, et subsoniques en présence d’un obstacle (sphère). À partir de l’analyse de la structure compressible des jets sous-détendus, de nouvelles lois sont aussi proposées pour estimer la position et le diamètre du disque de Mach puis la longueur du cône potentiel. Enfin, deux lois d’estimation du volume et de la masse inflammables de jets libres d’hydrogène sont proposées : elles sont fonction du débit massique de la fuite. / The aim of this work is to experimentally investigate the hydrogen dispersion to evaluate the impact of chronic or accidental leaks that may occur in a free or in a congested environment, from a low or a high pressure tank, on a fuel cell vehicle. The leaks are assimilated to vertical turbulent and axisymmetric jets with variable density. They are issued from cylindrical orifices from 1 to 3 mm diameter. An experimental set-up was designed to investigate the leaks : hydrogen has been replaced by helium for safety reasons. It supports a 200 bar pressure and allows to position an obstacle in the jet flow. The BOS (Background Oriented Schlieren) technique has been adapted to millimeter jets and to the presence of an obstacle. Particular attention has been given to the development of this technique. The results show good agrement with the available litterature data.New similarity laws are proposed, more representative of the flows of : subsonic free jets, under-expanded free jets, and subsonic impinging (a sphere) jets. From the structure analysis of compressible under-expanded free jets, new relations are proposed to estimate the Mach disk position and diameter as well as the potential core length. Finally, two new laws are found from the analysis of flammable volume and flammable mass : they are based on the leak mass flow rate.

Cette thèse est une contribution expérimentale à l’étude des résonances aéroacoustiques des jets sous-détendus : le Screech. Diverses méthodes expérimentales sont utilisées à ces fins, telles que la mesure de pression acoustique, la strioscopie et Vélocimétrie par Image des Particules, et associées à des techniques classiques de post-traitement comme les décompositions en mode de Fourier et aux valeurs propres. Ces Techniques permettent d’évaluer les effets d’épaisseur de la lèvre de la buse sur l’écoulement, et fournissent des informations sur les différences de comportement d’un même jet montrant des modes oscillatoires différents. Enfin, on entreprend d’étudier la présence de divers mécanismes de fermeture de la boucle de résonance pour divers modes de screech. La présence d’ondes intrinsèques du jet, se propageant vers l’aval pour les modes axisymétrique (A2) et hélicoïdal (C) suggèrent que ces ondes puissent jouer un rôle dans la résonance. La signature de ces ondes n’est en revanche pas attestée pour les modes battants (B). Ces résultats semblent donc indiquer que plusieurs mécanismes de rétroaction différents puissent être à l’oeuvre dans la résonance du jet sous-détendus. / This work provides an experimental contribution to the study of the Screech phenomenon. Various experimental techniques such as microphones array, Schlieren and Particle Image Velocimetry (PIV) together with advanced post-processing techniques like azimuthal Fourier decomposition and Proper Orthogonal Decomposition (POD) are employed. These techniques enable the evaluation of the lip thickness effects on the jets generated by two different round nozzles. The differences on the flow aerodynamics and acoustics are discussed. Then, we carry out experiments to analyse the effects of the different dominant Screech modes (B and C) on the flow characteristics. No noticeable differences are found in the mean fields. However, the fluctuation fields shows the contrary: B mode has larger fluctuation. In the last part, we investigate the Screech closure mechanism. The signature of upstream jet waves is revealed in the axisymmetric (A2) and helical (C) mode. However, the mode B does not present evidence of this instability in the flow, indicating that its closure mechanism may be bonded to another kind of waves. The conclusion from the results is that the Screech phenomenon seems be driven by different closure mechanisms.

A lateral jet in a supersonic crossflow creates a highly complex three-dimensional flow field which is not easily predicted. The aim of this research was to assess the use of a RANS based CFD method to simulate a lateral jet in supersonic crossflow interaction by comparing the performance of available RANS turbulence models. Four turbulence models were trialled in increasingly complex configurations; a flat plate, a body of revolution and a body of revolution at incidence. The results of this numerical campaign were compared to existing experimental and numerical data. Overall the Spalart-Allmaras turbulence model provided the best fit to experimental data. The performance of the lateral jet as a reaction control system was assed by calculating the force and moment amplification factors. The predicted flowfield surrounding the interaction was analysed in detail and was shown to predict the accepted shock and vortical structures. The lateral jet interaction flowfield over a body of revolution was shown to be qualitatively the same as that over a flat plate. An experimental facility was designed and manufactured allowing the study of the lateral jet interaction in Cranfield University’s 2 ½” x 2 ½” supersonic windtunnel. The interaction was studied with a freestream Mach number of 1.8, 2.4 & 3.1 and over a range of pressure ratios (50≤PR≤200). Levels of unsteadiness in the interaction were measured using high bandwidth pressure transducers. The level of unsteadiness was quantified by calculating the OASPL of the pressure signal. OASPL was found to increase with increasing levels of PR or MPR and to decrease with increases of Mach number. The levels of unsteadiness found were low with the highest levels found downstream of the jet.

No description available.

Un nouveau procédé propre et performant de préparation de surface reposant sur l'utilisation d'un jet d'azote sous très haute pression et basse température est développé. Comparativement aux techniques classiques de préparation de surface, ce nouveau procédé apporte des avantages environnementaux (pas de génération de déchets additionnels), techniques (modes d'action sur le revêtement à traiter) et énergétiques (consommation d'énergie). Cependant, le fonctionnement du procédé de traitement de surface en question, qui existe sous forme de prototype à l'heure actuelle, repose sur des bases empiriques et les phénomènes physiques qu'il met en jeu sont encore bien mal appréhendés. Pour comprendre ces phénomènes et les modéliser, une analyse numérique est présentée. Cette analyse s'inclut dans une stratégie d'étude qui vise à augmenter progressivement le degré de complexité de la modélisation. Un premier modèle vise à décrire l'évolution du fluide en temps et en espace en supposant un écoulement compressible visqueux et axisymétrique. Dans ce modèle, le fluide est considéré comme parfait. Une extension directe de ce modèle est ensuite présentée, elle propose de tenir compte des effets de fluide réel dans l'écoulement compressible monophasique. Ces modèles sont implémentés dans un code de calculs par volumes finis. Des cas tests sont étudiés afin de valider les modèles numériques. Une étude sur des configurations de type industrielles, représentatives des conditions d'utilisation du procédé de traitement de surface par jet d'azote, est ensuite menée / A new and efficient process of surface treatment is developed and exists in prototype form at present. The process aims at injecting an inert gas, such as nitrogen, at supercritical conditions through a nozzle. The jet resulting from the expansion of the fluid at high pressure and supercritical temperature impinges normally a flat surface. Alternatively to water jet technologies, which need expansive purification of water after use, and to other classical surface treatment process, this process provides environmental (no generation of additional waste), technical (action on the coating) and energetical (efficiency) benefits. However, the physical phenomena involved in the jet are still poorly understood at present. To understand and model these phenomena, numerical analysis is presented. This analysis is part of a study strategy that aims at gradually increasing the complexity of modeling. A first model aims at describing the evolution of the fluid in time and space, assuming a compressible axisymmetric viscous flow. In this model, the fluid is assumed perfect. A direct extension of this model is then presented where real fluid effects are taken into account in the compressible flow. These models are implemented in finite volume CFD code. Test cases are studied to validate the numerical models. A study of industrial-type configurations, representative of the conditions of use of the process of surface treatment by nitrogen is then conducted