O estudo de métodos de modelagem de escoamentos aerodinâmicos em regime transônico é de grande importância para a engenharia aeronáutica. O maior desafio no tratamento desses escoamentos está na sua característica não linear devido aos efeitos de compressibilidade e formação de ondas de choque. Tais efeitos não lineares influenciam no desempenho de superfícies aerodinâmicas em geral, bem como são responsáveis pelo aparecimento de fenômenos danosos para a resposta aeroelástica de aeronaves. O equacionamento para esses tipos de escoamentos pode ser obtido via as equações básicas da mecânica dos fluidos. No entanto, apenas soluções numéricas de tais equações são possíveis de ser obtidas de forma prática no presente momento. Para o caso específico do tratamento de problemas transônicos, as equações de Euler formam um conjunto de equações diferenciais a derivadas parciais capazes de capturar os efeitos não lineares de escoamentos compressíveis, porém os efeitos da viscosidade não são levados em consideração. O objetivo desse trabalho é implementar uma rotina computacional capaz de resolver numericamente escoamentos em regime transônico em torno de aerofólios. Para isso as equações de Euler não lineares são utilizadas e o campo de fluido ao redor de um perfil aerodinâmico é discretizado pelo método das diferenças finitas. Uma malha estruturada do tipo C discretizando o fluido ao redor de um aerofólio NACA0012 é considerada. A metodologia para solução numérica é baseada no método explícito de MacCormack de segunda ordem de precisão no tempo e espaço. Baseados na aproximação upwind, termos de dissipação artificial com coeficientes não lineares também são adicionados ao método. A solução do escoamento transônico estacionário ao redor do aerofólio NACA0012 é obtida e as principais propriedades do escoamento são apresentadas. Observa-se a formação de ondas de choque através de contornos de número de Mach ao redor do aerofólio. Gráficos das distribuições de pressão no intra e extradorso do aerofólio são mostrados, onde se identificam aos efeitos da brusca variação de pressão devido as ondas de choque. Os resultados são validados com valores de distribuição de pressão para o mesmo aerofólio encontradas na literatura técnica. Os resultados obtidos combinam bem com os fornecidos em códigos computacionais para solução do mesmo problema aerodinâmico / The study of aerodynamic modeling methods for the transonic flow regime is of great importance in aeronautical engineering. Major challenge on the treatment of those flows is on their nonlinear features due to compressibility effects and shock waves (appearance). Such nonlinear effects present a strong influence on aerodynamic performance, as well as they are responsible for harmful aeroelastic response phenomena in aircraft. Equations for transonic flows can be obtained from the basic fluid mechanic equations. However, only numerical methods are able to attain practical solutions for those set of differential equations at the present moment. For the specific case of treating transonic flow problems, the nonlinear Euler equations provide a set of partial differential equations with features to capture nonlinear effects of typical compressible flows, despite of not accounting for viscous flows effects. The aim of this work is to implement a computational routine for the numerical solution of transonic flows around airfoils. The Euler equations are used and the flow field around a aerodynamic profile is discretized by finite difference method. A C-type structured mesh is used to discretize the flow around a NACA0012 airfoil. The methodology for numerical solution is based on the explicit MacCormack method which has second order accuracy in time and space. Based on the upwind approximation, artificial dissipation with nonlinear coefficients is incorporated to the method. The steady transonic flow around the NACA0012 airfoil numerical solution is assessed and the main flow properties are presented. Shock wave structure can also be observed by means of the Mach number contours around the airfoil. Pressure distributions on upper and lower surfaces for different flow conditions are also shown, thereby allowing the observation of the effects of the abrupt pressure change due to shock waves. The results are validated using data presented in the technical literature. The present code solutions agree well with the solution obtained in other computational codes used for the same problem

O estudo de métodos de modelagem de escoamentos aerodinâmicos em regime transônico é de grande importância para a engenharia aeronáutica. O maior desafio no tratamento desses escoamentos está na sua característica não linear devido aos efeitos de compressibilidade e formação de ondas de choque. Tais efeitos não lineares influenciam no desempenho de superfícies aerodinâmicas em geral, bem como são responsáveis pelo aparecimento de fenômenos danosos para a resposta aeroelástica de aeronaves. O equacionamento para esses tipos de escoamentos pode ser obtido via as equações básicas da mecânica dos fluidos. No entanto, apenas soluções numéricas de tais equações são possíveis de ser obtidas de forma prática no presente momento. Para o caso específico do tratamento de problemas transônicos, as equações de Euler formam um conjunto de equações diferenciais a derivadas parciais capazes de capturar os efeitos não lineares de escoamentos compressíveis, porém os efeitos da viscosidade não são levados em consideração. O objetivo desse trabalho é implementar uma rotina computacional capaz de resolver numericamente escoamentos em regime transônico em torno de aerofólios. Para isso as equações de Euler não lineares são utilizadas e o campo de fluido ao redor de um perfil aerodinâmico é discretizado pelo método das diferenças finitas. Uma malha estruturada do tipo C discretizando o fluido ao redor de um aerofólio NACA0012 é considerada. A metodologia para solução numérica é baseada no método explícito de MacCormack de segunda ordem de precisão no tempo e espaço. Baseados na aproximação upwind, termos de dissipação artificial com coeficientes não lineares também são adicionados ao método. A solução do escoamento transônico estacionário ao redor do aerofólio NACA0012 é obtida e as principais propriedades do escoamento são apresentadas. Observa-se a formação de ondas de choque através de contornos de número de Mach ao redor do aerofólio. Gráficos das distribuições de pressão no intra e extradorso do aerofólio são mostrados, onde se identificam aos efeitos da brusca variação de pressão devido as ondas de choque. Os resultados são validados com valores de distribuição de pressão para o mesmo aerofólio encontradas na literatura técnica. Os resultados obtidos combinam bem com os fornecidos em códigos computacionais para solução do mesmo problema aerodinâmico / The study of aerodynamic modeling methods for the transonic flow regime is of great importance in aeronautical engineering. Major challenge on the treatment of those flows is on their nonlinear features due to compressibility effects and shock waves (appearance). Such nonlinear effects present a strong influence on aerodynamic performance, as well as they are responsible for harmful aeroelastic response phenomena in aircraft. Equations for transonic flows can be obtained from the basic fluid mechanic equations. However, only numerical methods are able to attain practical solutions for those set of differential equations at the present moment. For the specific case of treating transonic flow problems, the nonlinear Euler equations provide a set of partial differential equations with features to capture nonlinear effects of typical compressible flows, despite of not accounting for viscous flows effects. The aim of this work is to implement a computational routine for the numerical solution of transonic flows around airfoils. The Euler equations are used and the flow field around a aerodynamic profile is discretized by finite difference method. A C-type structured mesh is used to discretize the flow around a NACA0012 airfoil. The methodology for numerical solution is based on the explicit MacCormack method which has second order accuracy in time and space. Based on the upwind approximation, artificial dissipation with nonlinear coefficients is incorporated to the method. The steady transonic flow around the NACA0012 airfoil numerical solution is assessed and the main flow properties are presented. Shock wave structure can also be observed by means of the Mach number contours around the airfoil. Pressure distributions on upper and lower surfaces for different flow conditions are also shown, thereby allowing the observation of the effects of the abrupt pressure change due to shock waves. The results are validated using data presented in the technical literature. The present code solutions agree well with the solution obtained in other computational codes used for the same problem

Os recentes desenvolvimentos na tecnologia de sistema aeronáutico de geometria variável têm sido motivados principalmente pela necessidade de melhorar o desempenho de aeronaves. O conceito de Morphing Aircraft, por meio da variação da linha de arqueamento, representa uma alternativa para sistemas aeronáuticos mais eficientes. No entanto, para aeronaves de alto desempenho, projetos com estes novos conceitos podem gerar reações aeroelásticas adversas, o que representa uma questão importante e pode vir a limitar esses novos projetos. A compreensão adequada do comportamento aeroelástico devido à variação da linha de arqueamento, particularmente em regimes transônico, compreende uma questão importante. Este trabalho consiste num estudo preliminar das consequências aeroelásticas de um sistema aeronáutico de geometria variável. O objetivo desse trabalho é explorar as repostas aeroelásticas transônicas de um aerofólio com arqueamento variável no tempo. A metodologia para análise aeroelástica é baseada num modelo de seção típica. A integração no tempo do sistema aeroelástico é obtida pelo método de Runge-Kutta de quarta ordem. A representação do escoamento transônico não estacionário foi computada por um código CFD em um contexto de malhas não estruturadas com uma formulação dada pelas equações de Euler-2D. Esses resultados preliminares podem fornecer aos projetistas informações importantes sobre as respostas aeroelásticas de um sistema aeronáutico com variação da linha de arqueamento, permitindo estabelecer um quadro adequado para futuras investigações de controle aeroelástico de sistema aeronáutico de geometria variável. / Recent developments on aircraft variable geometry technologies have been mainly motivated by the need for improving the flight performance. The morphing wing concept, by means of variable camber, represents an alternative towards more efficient lifting surfaces. However, for higher performance aircraft, this technology may lead to designs that create unsteady loads, which may result in adverse aeroelastic responses, which represents an important and limiting issue. Proper understanding of the aeroelastic behavior, particularly in transonic flight regimes, due to variations in camber comprises an important matter. This work is a primary study of aeroelastic consequences of an real-time adaptive aircraft. The objective of this work is to investigate prescribed variations to airfoil camberline and their influence to the aeroelastic response in transonic flight regime. The methodology is based on computational simulations of typical section with unsteady transonic aerodynamics solved with a Computational Fluid Dynamics (CFD) code. The time integration of the aeroelastic system is obtained by Runge-Kutta fourth order. The unsteady transonic flow was computed by a CFD code based on the 2D-Euler equations with unstructured mesh. Prescribed camber variation of a symmetrical airfoil is transferred to the CFD mesh, and aeroelastic responses and loading is assessed. These preliminary results may provide the designers valuable information on the interaction between changes in camber during airfoil aeroelastic reactions, allowing establishing an adequate framework for further aeroelastic control investigations of morphing wings.

In this thesis, we present a collection of newly obtained results concerning the existence of vanishing viscosity solutions to the one-dimensional compressible Euler equations of gas dynamics, with and without geometric structure. We demonstrate the existence of such vanishing viscosity solutions, which we show to be entropy solutions, to the transonic nozzle problem and spherically symmetric Euler equations in Chapter 4, in both cases under the simple and natural assumption of relative finite-energy. In Chapter 5, we show that the viscous solutions of the one-dimensional compressible Navier-Stokes equations converge, as the viscosity tends to zero, to an entropy solution of the Euler equations, again under the assumption of relative finite-energy. In so doing, we develop a compactness framework for the solutions and approximate solutions to the Euler equations under the assumption of a physical pressure law. Finally, in Chapter 6, we consider the Euler equations in special relativity, and show the existence of bounded entropy solutions to these equations. In the process, we also construct fundamental solutions to the entropy equations and develop a compactness framework for the solutions and approximate solutions to the relativistic Euler equations.

An obvious reason for studying unsteady flows is the prediction of the effect of unsteady aerodynamic forces on a flight vehicle, since these effects tend to increase the likelihood of aeroelastic instabilities. This is a major concern in aerodynamic design of aircraft that operate in transonic regime, where the flows are characterised by the presence of adjacent regions of subsonic and supersonic flow, usually accompanied by weak shocks. It has been a common expectation that the numerical approach as an alternative to wind tunnel experiments would become more economical as computers became less expensive and more powerful. However even with all the expected future advances in computer technology, the cost of a numerical flutter analysis (computational aeroelasticity) for a transonic flight remains prohibitively high. Hence it is vitally important to develop an efficient, cheaper (in the sense of computational cost) and physically accurate flutter simulation tech nique which is capable of reproducing the data, which would otherwise be obtained from wind tunnel tests, at least to some acceptable engineering accuracy, and that it is essentially appropriate for industrial applications. This need motivated the present research work on exploring and developing efficient and physically accurate computational techniques for steady, unsteady and time-linearised calculations of transonic flows over an aircraft wing with moving shocks. This dissertation is subdivided into eight chapters, seven appendices and a bibliography listing all the reference materials used in the research work. The research work initially starts with a literature survey in unsteady transonic flow theory and calculations, in which emphasis is placed upon the developments in these areas in the last three decades. Chapter 3 presents the small disturbance theory for potential flows in the subsonic, transonic and supersonic regimes, including the required boundary conditions and shock jump conditions. The flow is assumed irrotational and inviscid, so that the equation of state, continuity equation and Bernoulli's equation formulated in Appendices A and B can be employed to formulate the governing fluid equation in terms of total velocity potential. Furthermore for transonic flow with free-stream Mach number close to unity, we show in Appendix C that the shocks that appear are weak enough to allow us to neglect the flow rotationality. The formulations are based on the main assumption that aerofoil slopes are everywhere small, and the flow quantities are small perturbations about their free-stream values. In Chapter 4, we developed an improved approximate factorisation algorithm that solves the two-dimensional steady subsonic small disturbance equation with nonreflecting far-field boundary conditions. The finite difference formulation for the improved algorithm is presented in Appendix D, with the description of the solver used for solving the system of difference equations described in Appendix E. The calculation of steady and unsteady nonlinear transonic flows over a realistic aerofoil are considered in Chapter 5. Numerical solution methods, based on the finite difference approach, for solving the two-dimensional steady and unsteady, general-frequency transonic small disturbance equations are presented, with the corresponding finite difference formulation described in Appendix F. The theories and solution methods for the time-linearised calculations, in the frequency and time domains, for the problem of unsteady transonic flow over a thin planar wing undergoing harmonic oscillation are presented in Chapters 6 and 7, respectively. The time-linearised calculations include the periodic shock motion via the shock jump correction procedure. This procedure corrects the solution values behind the shock, to accommodate the effect of shock motion, and consequently, the solution method will produce a more accurate time-linearised solution for supercritical flow. Appendix G presents the finite difference formulation of these time-linearised solution methods. The aim is to develop an efficient computational method for calculating oscillatory transonic aerodynamic quantities efficiently for use in flutter analyses of both two- and three-dimensional wings with lifting surfaces. Chapter 8 closes the dissertation with concluding remarks and future prospects on the current research work.

This paper will analyze the effects of using hybrid optimization methods for optimizing objective functions that are determined by computational fluid dynamics solvers for compressible viscous flow for optimal design of airfoils. Previous studies on this topic by the authors had examined the application of deterministic optimization methods and stochastic optimization methods such as Simulated Annealing and Simultaneous Perturbation Stochastic Analysis (SPSA). The studies indicated that SPSA method has a greater or equal efficiency as compared with SA method in reaching optimal airfoil designs for the design problem in question. However, in some situations SPSA method has a tendency to demonstrate an oscillatory behavior in the vicinity of a local optima. To overcome this tendency, a hybrid method designed to take full advantage of SPSA’s high rate of reduction of the objective function at the inception of the design process to drive the design cycles towards the optimal zone at first, and then combining with other methods to perform the final stages of the convergence towards the optimal solutions is considered. SPSA method has been combined with the gradient-based Broydon-Fletcher-Goldfarb-Shanno (BFGS) method as well as Simulated Annealing method for the transonic inverse airfoil design problem that is concerned with the specification of a target airfoil surface pressure distribution and starting from an initial guess of an airfoil shape, the target airfoil shape is reached by way of minimization of a quantity that depends on the difference between the target and current airfoil surface pressure distribution. For a typical transonic flow test case, the effects of using hybrid optimization techniques such as SPSA+BFGS and SPSA+SA as opposed to using SPSA alone can be seen in Figure 1. After 800 design cycles using SPSA, the hybrid SPSA+SA method took 2521 function evaluations of SA while the SPSA+BFGS method took 271 function evaluations to reach similar values which are much better than that reached by using SPSA alone in the entire minimization process. Results indicate that both of the two hybrid methods have capability to find a global optimum more efficiently than the SPSA method. The paper will address issues related to hybridization and its impact on the optimal airfoil shape designs in various contexts. / Singapore-MIT Alliance (SMA)

The current limits of supersonic parachute technology may constrain the ability to safely land future robotic assets on the surface of Mars. This constraint has led to a renewed interest in supersonic inflatable aerodynamic decelerator (IAD) technology, which offers performance advantages over the DGB parachute. Two supersonic IAD designs of interest include the isotensoid and tension cone, named for their respective formative structural theories. Although these concepts have been the subject of various tests and analyses in the 1960s, 1970s, and 2000s, significant work remains to advance supersonic IADs to a technology readiness level that will enable their use on future flight missions. In particular, a review of the literature revealed a deficiency in adequate aerodynamic and aeroelastic data for these two IAD configurations at transonic and subsonic speeds. The first portion of this research amended this deficiency by testing flexible IAD articles at relevant transonic and subsonic conditions. The data obtained from these tests showed that the tension cone has superior drag performance with respect to the isotensoid, but that the isotensoid may demonstrate more favorable aeroelastic qualities than the tension cone. Additionally, despite the best efforts in test article design, there remains ambiguity regarding the accuracy of the observed subscale behavior for flight scale IADs. Due to the expense and complexity of large-scale testing, computational fluid-structure interaction (FSI) analyses will play an increasingly significant role in qualifying flight scale IADs for mission readiness. The second portion of this research involved the verification and validation of finite element analysis (FEA) and computational fluid dynamic (CFD) codes for use within an FSI framework. These verification and validation exercises lend credence to subsequent coupled FSI analyses involving more complex geometries and models. The third portion of this research used this FSI framework to predict the static aeroelastic response of a tension cone IAD in supersonic flow. Computational models were constructed to mimic the wind tunnel test articles and flow conditions. Converged FSI responses computed for the tension cone agreed reasonably well with wind tunnel data when orthotropic material models were used and indicated that current material models may require unrealistic input parameters in order to recover realistic deformations. These FSI analyses are among the first results published that present an extensive comparison between FSI computational models and wind tunnel data for a supersonic IAD.

Os recentes desenvolvimentos na tecnologia de sistema aeronáutico de geometria variável têm sido motivados principalmente pela necessidade de melhorar o desempenho de aeronaves. O conceito de Morphing Aircraft, por meio da variação da linha de arqueamento, representa uma alternativa para sistemas aeronáuticos mais eficientes. No entanto, para aeronaves de alto desempenho, projetos com estes novos conceitos podem gerar reações aeroelásticas adversas, o que representa uma questão importante e pode vir a limitar esses novos projetos. A compreensão adequada do comportamento aeroelástico devido à variação da linha de arqueamento, particularmente em regimes transônico, compreende uma questão importante. Este trabalho consiste num estudo preliminar das consequências aeroelásticas de um sistema aeronáutico de geometria variável. O objetivo desse trabalho é explorar as repostas aeroelásticas transônicas de um aerofólio com arqueamento variável no tempo. A metodologia para análise aeroelástica é baseada num modelo de seção típica. A integração no tempo do sistema aeroelástico é obtida pelo método de Runge-Kutta de quarta ordem. A representação do escoamento transônico não estacionário foi computada por um código CFD em um contexto de malhas não estruturadas com uma formulação dada pelas equações de Euler-2D. Esses resultados preliminares podem fornecer aos projetistas informações importantes sobre as respostas aeroelásticas de um sistema aeronáutico com variação da linha de arqueamento, permitindo estabelecer um quadro adequado para futuras investigações de controle aeroelástico de sistema aeronáutico de geometria variável. / Recent developments on aircraft variable geometry technologies have been mainly motivated by the need for improving the flight performance. The morphing wing concept, by means of variable camber, represents an alternative towards more efficient lifting surfaces. However, for higher performance aircraft, this technology may lead to designs that create unsteady loads, which may result in adverse aeroelastic responses, which represents an important and limiting issue. Proper understanding of the aeroelastic behavior, particularly in transonic flight regimes, due to variations in camber comprises an important matter. This work is a primary study of aeroelastic consequences of an real-time adaptive aircraft. The objective of this work is to investigate prescribed variations to airfoil camberline and their influence to the aeroelastic response in transonic flight regime. The methodology is based on computational simulations of typical section with unsteady transonic aerodynamics solved with a Computational Fluid Dynamics (CFD) code. The time integration of the aeroelastic system is obtained by Runge-Kutta fourth order. The unsteady transonic flow was computed by a CFD code based on the 2D-Euler equations with unstructured mesh. Prescribed camber variation of a symmetrical airfoil is transferred to the CFD mesh, and aeroelastic responses and loading is assessed. These preliminary results may provide the designers valuable information on the interaction between changes in camber during airfoil aeroelastic reactions, allowing establishing an adequate framework for further aeroelastic control investigations of morphing wings.

Cette thèse a pour objectif d'améliorer la compréhension de la dynamique d'une interaction entre une onde de choc droite ou oblique et une couche limite laminaire ou turbulente. En particulier, nous nous sommes intéressés aux mécanismes responsables de l'apparition d'oscillations auto-entretenues basses fréquences. Ce phénomène survient dans de nombreux cas applicatifs comme dans des entrées d'air d'avions supersoniques, autour d'un profil d'aile en régime transsonique et au sein de tuyère en régime de sur-détente. La première partie de ce mémoire traite des différentes études réalisées pour déterminer la phénoménologie de ce type de dynamique. Dans un deuxième temps, nous expliquons la stratégie retenue pour effectuer notre étude qui consiste à développer un outil d'étude des instabilités, adapté à des écoulements turbulents présentant une interaction entre une onde de choc et une couche limite. Le développement d'un outil CFD linéarisé couplé à une méthode de résolution d'un problème aux valeurs propres par une approche dite sans matrice ou de « time-stepping », a permis la réalisation d'une telle étude. Après une étape de validation de notre outil, nous avons étudié des cas d'écoulements présentant une interaction entre une onde de choc et une couche limite. Trois cas en particulier ont été traités. Le premier cas correspond à une interaction entre une onde de choc oblique impactant une couche limite laminaire se développant sur une plaque plane. Ce cas est généralement qualifié dans la littérature de cas de « réflexion de choc ». Nous montrons qu'un tel écoulement est globalement stable et que sa dynamique peut être caractérisée par des mécanismes de réceptivité et par la réponse de l'écoulement vis-à-vis de perturbations extérieures. Les deux autres cas abordés dans ce travail ont été le cas d'un écoulement transsonique autour d'un profil d'aile de type NACA0012 en régime d'entrée en tremblement aérodynamique et un cas de tuyère transsonique plane de type Sajben en régime de sur-détente. Dans le premier cas, l'analyse de stabilité nous permet de mettre en évidence le phénomène de « buffet » sur le profil NACA0012, ce qui montre que le phénomène est lié à une instabilité globale linéaire . Dans le second cas, l'analyse de stabilité ne permet pas d'expliquer le phénomène auto-entretenues basses fréquences, et montre que l'écoulement est linéairement globalement stable. Dans ce cas, la dynamique est convective, transitoire et pilotée par des mécanismes de réceptivité. / The general purpose of this study is to provide a better understanding of the dynamics of an interaction between a shock wave and a laminar or turbulent boundary layer. In particular, we were interested in mechanisms responsible for the emergence of low-frequency self-sustained oscillations. This phenomenon arises in numerous industrial cases as in air inlets of supersonic aircrafts, around a profile of wing in transonic regime and within over-extended nozzle. The first part of this report handles various studies carried out to determine the phenomenology of this kind of dynamics. Secondly, we explain the strategy adopted to make our study which consists in developing a tool of study of the instabilities adapted to turbulent flows including an interaction between a shock wave and a boundary layer. The development of a linearized CFD tool coupled with a method of resolution of a eigenvalue problem by a free-matrix approach ( " time-stepping " approach), allowed the realization of such a study. After a stage of validation of our tool, we studied cases of flows including an interaction between a shock wave and a boundary layer. Three cases in particular were handled. The first case corresponds to an interaction enter an oblique shock wave impacting on a laminar boundary layer developing on a flat plate. This case is generally qualified in the literature of case as " reflected shock wave". We show that such a flow is globally stable and that the dynamics of such a flow behaves as a selective noise amplifier, the dynamic is mainly driven by receptivity mechanisms and by the response of upstream white nose disturbance. Two other cases have been studied on this work, the case of a transonic flow around a profile wing of NACA0012 type around the onset of buffet phenomenon and the case of transonic nozzle of Sajben type on over-extended regime. In the first case, the global stability analysis allows us to highlight the buffet phenomenon of on the profile NACA0012, what shows that the phenomenon is linked to a linear global instability. In the second case, the analysis of stability does not allow to explain the self-sustained low frequencies phenomenon, and shows that the flow is linearly globally stable. In this case, the dynamics is convective, passing and piloted by receptivity mechanisms.

En phase d’approche, le bruit rayonné par l’entrée d’air des turboréacteurs est principalement dû aux interactions entre le rotor et le stator. Cependant les ondes de choc (ou ondes en N) générées par le rotor en régime transsonique peuvent devenir une source de bruit dominante durant le décollage et la montée de l’appareil. L’étude des ondes en N nécessite de se concentrer sur deux processus majeurs : 1) la génération des chocs par un rotor parfait (dont toutes les aubes sont identiques) et par un rotor réel (en tenant compte des irrégularités géométriques des aubes), et 2) la propagation de ces ondes en N à travers la nacelle, produisant du bruit dont le spectre se compose des harmoniques de la fréquence de passage des aubes pour un rotor régulier, et des harmoniques aux fréquences multiples de la rotation du rotor (FMR) pour un rotor irrégulier. Plusieurs approches analytiques et numériques ont été développées durant les 40 dernières années.Cette thèse relate dans un tout premier temps les principales théories de la propagation des ondes de choc ainsi que les modèles majeurs de génération de FMR. Une attention particulière est portée sur les liens entre les équations générales de la mécanique des fluides et ces modèles de propagation non linéaire afin de mettre en évidence les différentes hypothèses formulées dans ces modèles. Dans un deuxième temps, les principales méthodes semi-analytiques de génération et de propagation des chocs seront évaluées et comparées en les appliquant à des configurations de turboréacteurs. En outre, un nouveau modèle de génération de FMR basé sur des considérations géométriques est élaboré par l’intermédiaire d’une campagne d’essais comportant d’une part des mesures de signaux de pression dans la nacelle et d’autre part les mesures des angles de calage des aubes pendant le fonctionnement du moteur. Le deuxième volet de la thèse concerne le développement d’une méthodologie de simulation numérique basée sur l’utilisation du code elsA de l’ONERA en résolvant les équations d’Euler (approche CAA). L’objectif de cette approche est de s’affranchir des limitations des modèles de propagation semi-analytiques et de tenir compte de la géométrie réelle de la nacelle ainsi que d’un écoulement réaliste. Des ondes de choc régulières et irrégulières sont directement injectées dans un plan proche de la soufflante et se propagent en remontant l’écoulement. Ces ondes de choc sont injectées par l’intermédiaire d’une condition limite de non-réflexion qui nécessite d’imposer le champ conservatif. La signature des chocs peut provenir d’un RANS, de mesures ou d’un signal analytique. Étant donné que les mesures ou le signal théorique ne permettent d’obtenir que la pression, une méthode de reconstruction du champ conservatif à partir des variations de pression induites par le choc a été élaborée. Cette méthode d’injection est tout d’abord appliquée à un conduit annulaire infiniment mince et validée par la méthode de propagation semi-analytique de McAlpine & Fisher. Ensuite, les effets de propagation 3D sont étudiés en augmentant l’épaisseur du conduit. Enfin, la méthode CAA est appliquée à des configurations de turboréacteurs modernes et des ondes de choc régulières et irrégulières sont propagées numériquement. Les résultats sont comparés aux solutions RANS ainsi qu’aux mesures disponibles. / Whereas the sound radiated from the inlet of turbofans is mainly due to rotor–stator interactions in approach flight, the shock waves (or N-waves) emitted by the rotor at transonic rotation speeds can be a dominant noise source during takeoff and climb. The study of N-waves needs to take account of two main processes: 1) the generation of N-waves for a perfect rotor (in which all blades are identical) and for a real rotor (considering small geometrical blade dispersion), and 2) the N-wave propagation through the inlet duct producing the blade passing harmonics for a perfect rotor, and the multiple pure tones (harmonics of the rotation frequency) for a real rotor. Several analytical and numerical approaches have been investigated for the past 40 years.This thesis first intends to relate the main propagation theories and to address the foremost MPT generation method hypotheses. The links between fluid dynamics equations and practical non-linear theories currently adopted are emphasized and discussed. In a second step, the main relevant semi-analytical methods are cross-checked by applying them to representative turbofan configurations. Moreover, a novel model of irregular N-wave generation based on geometrical considerations is investigated thanks to test data related to in-duct pressure signatures and blade stagger angle measurements during the engine operation. Then, a second part of the work investigates a numerical strategy based on elsA ONERA code, solving the full Euler’s equations (CAA approach). The objective is to prevent from the limitations of 2D analytical models and to take into account actual inlet geometry and realistic convection flow. Regular and non-regular shock waves are directly injected in a plane close to the fan and propagated through the inlet. These shock waves are injected through a non-reflective boundary condition which requires the conservative field. The initial shock description near the fan is provided either by a RANS computation or by experiment, or else from analytical model. As experiment or analytical signals only provide pressure signatures, a theory is set up to re-built the whole conservative field from the basis of a pressure shockwave. This injection method is firstly applied on an infinitely narrow annular duct and validated through the comparison with the McAlpine & Fisher analytical method. Then, the 3D propagation effects are pointed out by increasing the duct height. Finally, the CAA method is applied on actual intake geometry of modern turbofan demonstrators, and propagation of regular and irregular shock-waves are simulated. The numerical results are compared to RANS solutions and to available measurements.