Afin de concevoir des systèmes de propulsion innovants, l'amélioration des performances des prises d'air supersonique constitue un enjeu majeur. En particulier, les écoulements intervenant au sein des entrées d'air et/ou de diffuseurs supersoniques mettent en jeu des phénomènes complexes associés aux diverses échelles spatiales et temporelles: dynamique de la turbulence pariétale, interaction entre une onde de choc et une couche limite turbulente, décollements tridimensionnels et effets de coins. Malgré les contributions significatives et récentes des simulations numériques de haute fidélité sur les instationnarités de l'interaction onde de choc/couche limite sans paroi latérales, peu d'études numériques ont été menées sur l'influence des coins dans la dynamique de l'écoulement. En présence de parois latérales et à des nombres de Mach suffisamment élevés, l’interaction se modifie et un train de choc se forme dans le diffuseur. Dans le cadre de cette thèse, les équations de Navier-Stokes en régime compressible sont résolues à l'aide de schémas d'ordre élevé. Des simulations en régime supersonique de l'écoulement dans des diffuseurs rectangulaires de largeurs différentes sont effectuées. L'étude permet la mise en évidence de l'influence du confinement et des effets de coins. Une deuxième partie de l’étude est consacrée à la compréhension des instationnarités générées par un train de choc dans un diffuseur rectangulaire à l'aide d'outils de post-traitement avancés: décomposition modale dynamique et périodogramme. Les résultats montrent la présence d'un possible phénomène de résonance du diffuseur à des fréquences proches de celles émises par l'écoulement. / To design innovative propulsion systems, improving the performance of supersonic air intakes is a major issue. In particular, the flows through the air intakes and/or supersonic diffusers involve complex unsteady phenomena associated with various spatial and temporal scales such as: wall-bounded turbulence dynamics, interaction between a shock-wave and a turbulent boundary layer, three-dimensional separated flows and corners effects. Despite the significant contributions from recent high-fidelity simulations of unsteady shock-wave boundary layer interaction in the absence of side walls, few numerical studies were conducted with secondary flows due to corner effects. In the presence of side walls and at Mach numbers large enough, the topology of the interaction is modified and a shock-train forms in the diffuser. In this thesis, the Navier-Stokes compressible equations are solved using high-order schemes. Simulations of supersonic flows in rectangular diffusers of different widths are carried out. The study allows to highlight the influence of confinement and corners effects on the mean flow. A second part of the study is devoted to the understanding of the unsteadiness associated with a shock-train in a rectangular supersonic diffuser. For that purpose, advanced post-processing tools have been developed such as: dynamic mode decomposition and Fourier analysis. The results show the presence of a possible resonance phenomenon in the diffuser at frequencies close to those associated with the flow.

This master thesis is focused on aerodynamic design of transonic glider, which is assigned for following an airliner at cruising regime of flight. Main goal of the thesis is to determine basic geometrical design of airplane with respect to Whitcomb aera rule, mass analysis and drag polar. Weight analysis includes determining center of gravity according to longitudinal static stability margin. The drag polar of the airplane is determine for cruising regime of flight, that means the regime of following an airliner.

Aeroelasticita je nezbytná vědní disciplína zahrnuta do návrhu letounů. Zaměřuje se na předpovídání jevů, které vznikají vlivem interakce aerodynamických, elastických a setrvačných sil. Tyto jevy často vedou ke katastrofickým následkům, proto musí být prokázáno, že nevzniknou v rozsahu rychlostí ohraničujících letovou obálku. Aplikace moderních materiálů při konstrukci draku, spolu se snahou navrhnout aerodynamicky efektivní tvar křídel, vede ke zvyšování poddajnosti letounů. To má za následek změnu aerodynamických vlastností a také k výraznějšímu vliv na aeroelastické jevy, které mohou být vyvolány snadněji vlivem pohybů tuhého tělesa než v případě tužších konstrukcí. Aeroelastické jevy mohou vznikat v širokém rozsahu rychlostí zahrnujícím i transsonickou oblast. V této oblasti je ovlivněna zejména rychlost, při níž dochází k třepetání, a to vlivem nelineárních jevů v proudu. Běžné nástroje, které jsou založeny na lineárních teoriích, nejsou schopny tyto nelineární jevy popsat. Cílem práce je proto navrhnout, implementovat a otestovat nástroj pro výpočetní (numerickou) simulaci aeroelasticity. Nástroj má využívat řešič proudového pole, který je schopen předpovědět nelineární jevy. V práci je kladen důraz na simulaci statické aeroelasticity. V práci jsou popsány metody, které je nutno zahrnout do numerické simulace statické aeroelasticity. Dále je popsán vlastní nástroj a je provedeno zhodnocení konvergence statických aeroelastických výpočtů. Funkčnost nástroje byla ověřena na příkladech, kdy byly použity různé aerodynamické a strukturální modely. Nástroj byl také aplikován při aerodynamické tvarové optimalizaci poddajného křídla. Výsledky optimalizace a její výpočetní náročnost byly porovnány s případem optimalizace tuhého křídla. Na závěr je v práci prezentován příspěvek autora do výzkumu zaměřeného na zhodnocení vlivu časové synchronizace mezi CFD a CSM řešiči. Použitý testovací případ je transsonické obtékání křídla na začátku třepetání (flutteru). Výsledky byly srovnány s experimentálními daty poskytnutými NASA.

Cette thèse concerne l’aérodynamique de pales d'extrémité transsonique. Ces pales sont conçues pour maximiser le rendement en croisière, tout en générant la traction requise au décollage. Elles ont des profils fins et peu cambrés, travaillant à forte incidence au décollage, ce qui peut entraîner l’apparition d’un tourbillon de bord d’attaque (TBA). Or ce TBA présente des similitudes avec les tourbillons d’apex d’aile Delta, connus pour leur capacité à générer de la portance tourbillonnaire.Cette étude consiste à examiner l’intérêt du TBA pour les performances aérodynamiques.La démarche a consisté dans un premier temps à caractériser la topologie du TBA sur une maquette représentative d’une pale d’ Open Rotor, à l'aide d'essais PIV résolus en temps et de calculs RANS k-omega SST, et à évaluer la capacité de la simulation RANS à reproduire les caractéristiques d’intérêt pour cette étude. Un algorithme a été développé afin d'estimer la contribution de ce TBA à la portance à partir du champ de pression pariétal RANS.Afin d'expliciter l'influence des paramètres géométriques et de fonctionnement de la pale sur la portance tourbillonnaire, un modèle 1D de la portance tourbillonnaire a été développé puis couplé à la méthode de l'élément de pale.Les premières comparaison de géométries à iso-traction ont montré que la portance tourbillonnaire permet de générer la traction requise au décollage avec une surface alaire plus faible. Ces résultats ouvrent de nouvelles perspectives pour la conception de géométries avec un meilleur rendement en croisière. / This thesis deals with the aerodynamic properties of propeller blades. Those blades are designed to maximise cruise efficiency, while achieving the target thrust at take-off. Their thin, low-cambered profiles must work at high incidence at take-off, which may give rise to a leading-edge vortex (LEV).The topology of this LEV looks similar to Delta wing LEVs, which are known to generate vortex lift.the aim of this study is to explore the probable impact of the LEV on lift at take-off in order to reconsider propeller blade designs. The approach first consisted in caracterising the LEV topology on a model blade representative of an Open Rotor front blade, using both Time-Resolved PIV and RANS k-omega SST calculations. The comparison between both methods demonstrated the ability of RANS calculations to reproduce the LEV characteristics of interest to this study.Then, the LEV contribution to lift was evaluated thanks to an algorithm developed to estimate vortex lift contribution from RANS wall pressure fields.In order to explicit the influence of the blade's geometrical and functioning parameters on vortex lift, a 1D vortex lift model was developed and coupled to the Blade Element Momentum Theory.The first blade geometry comparative studies at iso-thrust showed that vortex lift enables to generate target thrust at take-off with a lower blade surface. This opens new perspectives for the design of blade geometries with enhanced cruise efficiency.

<p>We study the effect of minute perturbations by using blowing jets at upstream and bio-inspired micro denticles on turbulence large-scale motions which are observed to be crucial in controlling heat transfer, noise and drag reduction. This work is divided into two phases. In first phase, we studied the effect of blowing perturbations at upstream on large-scale motions and associated co?herent vortical structures which are crucial in enhancing heat transfer by promoting mixing. The second phase is focused on impact of flow dynamics in preventing the biofouling using micro bioinspired structures and the importance of flow regime in designing the antifouling coating us?ing bioinspired structures is demonstrated, and subsequently, separation bubble dynamics and its characterization is carried out for a transonic channel imposed with pressure gradient to further expand our thesis outcomes to utilize micro bioinspired structures in aerospace applications, noise reduction, and to delay separation.</p> <p><br></p> <p>Extensive studies were focused on the importance of large-scale motions (LSM) and their con?tribution to TKE and turbulence mixing. Although there are studies focusing on the λ2 coherent vortical structures and large-scale motions separately, there are no studies addressing the control?ling using upstream perturbations on the large-scale motions and their associated λ2 vortices. In the first phase of our studies, we used the DNS data of channel flow for Reτ = 394 generated using in-house code. In these simulations, we created blowing perturbations using spanwise jets of low blowing ratio, 0.2, placed at upstream. The spatial large-scale motions are extracted using a a novel 3D adaptive Gaussian filtering technique developed based on Lee and Sung [1] for turbulent pipe flows. POD is used to extract the energetic large-scale motions and coherent vortical structures are extracted using λ2-criterion for its efficiency in educing coherent structures in cross flow jets. The results show that the upstream perturbations enhance streamwise heat flux via energetic LSM and also create a secondary peak of scalar production in the log-layer showing that the perturbations alter LSMs to enhance the heat transfer. Filtered large-scale field from Gaussian filtering technique have an integral length scale greater than 2h (where h is channel half-height) are used to obtain λ2 vortices. The resulted λ2 vortices are of ring-type and have higher signature of temperature than their counterpart. The pre-multiplied spectra shows that the upstream perturbations can excite the large-scale wave-numbers which are in the same order as the jet diameter and spacing between them. Simulations show the presence of secondary peak in the log-layer and increased turbulence production which are eminent of large-scales. Furthermore, our results suggest that jet spacing and diameter are crucial in exciting large-scale field to control turbulent flows.</p> <p><br></p> <p>Evans, Hamed, Gorumlu, et al. [2] modeled the denticles present on Mako shark skin into a diverging micro-pillars. They conducted experimental studies in a water tunnel using these on the back of airfoil exposed to an adverse pressure gradient flow. They observed that presence of these pillars reduced the re-circulation bubble (form drag) by 50%. They proposed a blowing and suction type mechanism by which the micro pillars interact with the boundary layer. However, the details of underlying interfacial mechanism is not completely understood. The unique impact of flow conditions on anti-biofouling and the corresponding mechanisms for the first time is illustrated. We employed commercially available bioinspired structures as micro-diverging pillars making it feasible to apply in real life. We demonstrated the underlying mechanism by which bio?inspired structures are responsible for anti-biofouling. To study the pressure gradient effects on the separation under transonic conditions, we performed direct numerical simulations (DNS) in a non?equilibrium flow created by a sinsuoidal contraction and also, we quantified the separation length,</p> <p>detachment, and attachment points of separation bubble imposed with various pressure gradients and their variation in the transonic and subsonic regimes. We noticed that the resultant shear at the attachement led to the enhancement of coherent structures which are extended into the outer layer under transonic flow which is quite different than the subsonic flow.</p>

Multiple high-fidelity, time-accurate computational fluid dynamics simulations were performed to investigate the effects of upstream stator loading and rotor shock strength on vortex shedding characteristics in a single stage transonic compressor. Various configurations of a transonic axial compressor stage, including three stator/rotor axial spacings of close, mid, and far in conjunction with three stator loadings of decreased, nominal, and increased were simulated in order to understand the flow physics of transonic blade-row interactions. Low-speed compressors typically have reduced stator/rotor axial spacing in order to decrease engine weight, and also because there is an increase in efficiency with reduced axial spacing. The presence of a rotor bow shock in high-speed compressors causes additional losses as the shock interacts with the upstream stator trailing edge. This research analyzes the strength of shock-induced vortices due to these unsteady blade-row interactions. The time-accurate URANS code, TURBO, was used to generate periodic, quarter annulus simulations of the Blade Row Interaction compressor rig. Both time-averaged and time-accurate results compare well with experimentally-observed trends. It was observed that vortex shedding was synchronized to the passing of a rotor bow shock. Normal and large shock-induced vortices formed on the stator trailing edge immediately after the shock passing, but the large vortices were strengthened at the trailing edge due to a low-velocity region on the suction surface. This low velocity region was generated upstream of mid-chord on the suction surface from a shock-induced thickening of the boundary layer or separation bubble, due to the rotor bow shock reflecting off the stator trailing edge and propagating upstream. The circulation of the shock-induced vortices increased with shock strength (decreased axial spacing) and stator loading. Most design tools do not directly account for unsteady effects such as blade-row interactions, so a model is developed to help designers account for shock-induced vortex strength with varying shock strength and stator loading. An understanding of the unsteady interactions associated with blade loading and rotor shock strength in transonic stages will help compressor designers account for unsteady flow physics early in the design process.

No description available.

Multidisciplinary design optimization (MDO) requires many designs to be evaluated while searching for an optimum. As a result, the calculations done to evaluate the designs must be quick and simple to have a reasonable turn-around time. This makes aerodynamic calculations in the transonic regime difficult. Running computational fluid dynamics (CFD) calculations within the MDO code would be too computationally expensive. Instead, CFD is used outside the MDO to find two-dimensional aerodynamic properties of a chosen airfoil shape, BACJ, at a number of points over a range of thickness-to-chord ratios, free-stream Mach numbers, and lift coefficients. These points are used to generate surrogate models which can be used for the two-dimensional aerodynamic calculations required by the MDO computational design environment. Strip theory is used to relate these two-dimensional results to the three-dimensional wing. Models are developed for the center of pressure location, the lift curve slope, the wave drag, and the maximum allowable lift coefficient before buffet. These models have good agreement with the original CFD results for the airfoil. The models are integrated into the aerodynamic and aeroelastic sections of the MDO code. / Master of Science

The study presented herein describes and analyzes a detailed experimental investigation of the effects of density ratio on endwall thermal performance at varying blowing rates for a typical nozzle guide vane platform purge jet cooling scheme. An axisymmetric converging endwall with an upstream doublet staggered cylindrical hole purge jet cooling scheme was employed. Nominal exit flow conditions were engine representative and as follows: {rm Ma}_{Exit} = 0.85, {rm Re}_{Exit,C_{ax}} = 1.5 times {10}^6, and large-scale freestream Tu = 16%. Two blowing ratios were investigated corresponding to the upper and lower engine extrema. Each blowing ratio was investigated amid two density ratios; one representing typical experimental neglect of density ratio, at DR = 1.2, and another engine representative density ratio achieved by mixing foreign gases, DR = 1.95. All tests were conducted on a linear cascade in the Virginia Tech Transonic Blowdown Wind Tunnel using IR thermography and transient data reduction techniques. Oil paint flow visualization techniques were used to gather quantitative information regarding the alteration of endwall flow physics due two different blowing rates of high-density coolant. High resolution endwall adiabatic film cooling effectiveness, Nusselt number, and Net Heat Flux Reduction contour plots were used to analyze the thermal effects. The effect of density is dependent on the coolant blowing rate and varies greatly from the high to low blowing condition. At the low blowing condition better near-hole film cooling performance and heat transfer reduction is facilitated with increasing density. However, high density coolant at low blowing rates isn't adequately equipped to penetrate and suppress secondary flows, leaving the SS and PS largely exposed to high velocity and temperature mainstream gases. Conversely, it is observed that density ratio only marginally affects the high blowing condition, as the momentum effects become increasingly dominant. Overall it is concluded density ratio has a first order impact on the secondary flow alterations and subsequent heat transfer distributions that occur as a result of coolant injection and should be accounted for in purge jet cooling scheme design and analysis. Additionally, the effect of increasing high density coolant blowing rate was analyzed. Oil paint flow visualization indicated that significant secondary flow suppression occurs as a result of increasing the blowing rate of high-density coolant. Endwall adiabatic film cooling effectiveness, Nusselt number, and NHFR comparisons confirm this. Low blowing rate coolant has a more favorable thermal impact in the upstream region of the passage, especially near injection. The low momentum of the coolant is eventually dominated and entrained by secondary flows, providing less effectiveness near PS, near SS, and into the throat of the passage. The high momentum present for the high blowing rate, high-density coolant suppresses these secondary flows and provides enhanced cooling in the throat and in high secondary flow regions. However, the increased turbulence impartation due to lift off has an adverse effect on the heat load in the upstream region of the passage. It is concluded that only marginal gains near the throat of the passage are observed with an increase in high density coolant blowing rate, but severe thermal penalty is observed near the passage onset. / Master of Science / Gas turbine technology is used frequently in the burning of natural gas for power production. Increases in engine efficiency are observed with increasing firing temperatures, however this leads to the potential of overheating in the stages following. To prevent failure or melting of components, cooler air is extracted from the upstream compressor section and used to cool these components through various highly complex cooling schemes. The design and operational adequacy of these schemes is highly subject to the mainstream and coolant flow conditions, which are hard to represent in a laboratory setting. This experimental study explores the effects of various coolant conditions, and their respective response, for a purge jet cooling scheme commonly found in engine. This scheme utilizes two rows of staggered cylindrical holes to inject air into the mainstream from platform, upstream of the nozzle guide vane. It is the hope that this air forms a protective layer, effectively shielding the platform from the hostile mainstream conditions. Currently, little research has been done to quantify these effects of purge flow cooling scheme while mimicking engine geometry, mainstream and coolant conditions. For this study, an endwall geometry like that found in engine with a purge jet cooling scheme is studied. Commonly, an upstream gap is formed between the combustor lining and first stage vane platform, which is accounted for in this testing. Mainstream and coolant flow conditions can have large impacts on the results gathered, so both were matched to engine conditions. Varying of coolant density and injection rate is studied and quantitative results are gathered. Results indicate coolant fluid density plays a large role in purge jet cooling, and with neglection of this, potential thermal failure points could be overlooked This is exacerbated with less coolant injection. Interestingly, increasing the amount of coolant injected decreases performance across much of the passage, with only marginal gains in regions of complex flow. These results help to better explain the impacts of experimental neglect of coolant density, and aid in the understanding of purge jet coolant injection.

L’étude effectuée au cours de cette thèse a permis de caractériser numériquement les instabilités d’origine aérodynamique rencontrées dans un compresseur centrifuge dessiné par Turbomeca. Ce compresseur est composé d’une roue directrice d’entrée, d’un rouet centrifuge, d’un diffuseur radial et de redresseurs axiaux. Le module expérimental, dénommé Turbocel, sera accueilli au LMFA courant 2016. Le contenu de cette étude repose donc exclusivement sur des résultats numériques dont certains sont cependant comparés à des résultats expérimentaux partiels obtenus par Turbomeca sur une configuration proche. _ Le fonctionnement du compresseur est analysé à différentes vitesses de rotation, à partir de simulations RANS et URANS menées avec le code elsA. Du point de vue de la méthodologie, deux points importants sont à retenir :- Du fait du caractère transsonique de l’écoulement dans le rouet et le diffuseur radial à haut régime de rotation, les simulations RANS stationnaires ne permettent pas d’accéder à une description satisfaisante des phénomènes physiques. Cela est dû à l’utilisation d’un plan de mélange aux différentes interfaces rotor-stator qui a pour effet d’empêcher les ondes de choc de remonter à l’amont, et qui affecte tant la physique de l’écoulement que l’étendue de la plage de fonctionnement stable.- En-dessous d’un certain débit, les calculs URANS sur période machine révèlent que le comportement de l’étage n’obéit plus à la périodicité spatio-temporelle mono-canal. Une plage instable est alors obtenue à toutes les iso-vitesses simulées. A bas régime de rotation, une autre plage stable existe lorsque le compresseur est suffisamment vanné. L’étage retrouve alors une périodicité spatio-temporelle, à condition d’étendre le domaine de calcul dans le stator à deux canaux inter-aubes. En ce qui concerne les limites de stabilité de Turbocel, différentes évolutions sont décrites selon la vitesse de rotation considérée :- A haut régime de rotation, une basse fréquence commence à émerger près du point de rendement maximal et son intensité ne fait qu’augmenter jusqu.au pompage.- A bas régime, une signature basse fréquence comparable se manifeste près du point de rendement maximal mais disparaît passé un certain vannage, et n’est donc présente que sur une plage de débit délimitée. La seconde zone stable peut alors être numériquement parcourue jusqu.au pompage proprement dit. La signature basse fréquence est imputée à l’instauration d’une recirculation dans l’inducteur qui une fois établie est quasi-stationnaire. Les résultats numériques mettent en évidence que la source d’instabilité sévère sur Turbocel provient du diffuseur aubé. En fonction du point de fonctionnement, ce composant adopte des comportements différents, entre lesquels une certaine continuité existe, et ses performances chutent progressivement lorsque le débit diminue. Au final, les domaines de stabilité de l’étage de compression peuvent être reliés au type d’écoulement qui se développe dans le diffuseur radial, et apparaissent dictés par le diffuseur semi-lisse à haut régime de rotation. Enfin, afin d’étendre les plages de fonctionnement stable, une stratégie de contrôle basée sur l’aspiration de couche limite dans le diffuseur aubé a également été déterminée dans le cadre de cette thèse. Son évaluation fera l’objet d’études ultérieures sur Turbocel. / The present study aims at characterizing the aerodynamic instabilities involved in a centrifugal compressor designed by Turbomeca, by means of numerical simulation. This compressor is composed of inlet guide vanes, a centrifugal impeller, a radial vaned diffuser and axial outlet guide vanes. The test module, named Turbocel, will be delivered to the LMFA in 2016. Thus, the results presented in this manuscript are only based on CFD, although some of them are compared to experimental results obtained by Turbomeca on a close configuration.RANS and URANS simulations are performed for several rotational speeds, using the elsA software.Two methodological key points are to be emphasized:- As the flow in both the impeller and the radial diffuser is transonic at high rotational speed, steady RANS simulations cannot provide a satisfactory description of the physical phenomena taking place. This can be explained by the use of the mixing plane approach which prevents shock waves to extend upstream the rotor-stator interfaces, and which impacts the flow field predicted as well as the prediction of the stable operating range.- Below a given massflow rate, URANS simulations covering the spatial period of the compressor prove that the stage behavior does not obey to the single passage spatio-temporal periodicity anymore. An unstable operating range then appears at all the simulated rotational speeds. At low rotational speed, another stable range is however obtained if the compressor is further throttled’ A new periodicity arises on this massflow range, provided that the stator domain is extended to two neighboring blade passages. Concerning the stability domains of Turbocel, different evolutions are obtained depending on the rotational speed:- At high rotational speed, a low frequency phenomenon starts to develop near the peak efficiency point and its intensity keeps increasing until surge happens.- At low rotational speed, a low frequency signature also appears near the peak efficiency point, but it then vanishes when the compressor is further throttled, so that only a restricted operating range exhibits this instability. It then gives rise to a second stable operating range which can be described numerically, ending with surge itself. The low frequency signature is attributed to the enhancement of a flow recirculation in the inducer which, once fully established, is quasi-steady. The numerical results underline that the source of severe instability in the compressor comes from the vaned diffuser. Depending on the operating point, this component can adopt different behaviors, between which a relative continuity exists, and its performances decrease when the massflow rate decresases. The overall stage performances prove that at high rotational speed, the global stability is driven by the semi-vaneless diffuser and depends on the flow developing in the radial diffuser. Finally, in order to extend the stable operating range of the compressor, a flow control strategy based on boundary layer suction has also been determined in the diffuser. Its impact on the performances of Turbocel will be deeply studied later on.