The resolution of the complex physics of rotating configurations is critical for any engineering analysis that requires multiple frames of reference. Two well-known applications are in the rotorcraft and wind energy industries. Rotor wake impingement from rotor-fuselage and wind turbine-tower interactions impact structural and acoustic characteristics. Additionally, parasite drag resulting from rotorcraft hubs may result in severe limitations on forward flight vehicle performance. Complex turbulent wakes from rotors and hubs impinging on downstream empennage can create adverse aeroelastic behavior and can affect handling qualities. Numerical simulations of these flows require state-of-the-art Navier Stokes methods using dynamic overset grids. However, many current methods typically used in industry result in wakes that dissipate essential features. In order to address these concerns, two advancements are introduced in this thesis. Feature-based grid adaptation on dynamic overset grids has been developed and demonstrated with an unstructured Navier Stokes solver. The unique feature of the adaptation technique is that it is applied globally on the overset grid system except within the boundary layer. In concert with grid adaptation, an efficient parallelized search algorithm for solution interpolation over massively distributed systems has been created. This results in cost-effective interpolation that retains the numerical order of accuracy and has been verified in both space and time. The improvements have been demonstrated for rotor-fuselage interaction and a generic rotating hub. Detailed analysis of convergence of the methodology and sensitivity of the results to relevant parameters have also been included.

Tidal bores are surge waves propagating upstream rivers as the tide rushes into estuaries. They induce large turbulences and mixing of the river and estuary flow of which effects remain scarcely studied. Herein, tidal bores are investigated experimentally and numerically with an idealised model of positive surges propagating upstream an initially steady flow. The experimental work estimated flow changes and typical turbulent length scale evolution induced by undular bores with and without breaking roller. The bore passage was associated with large free surface and flow velocity fluctuations, together with some variations of the integral turbulent scales. Coherent turbulent structures appeared in the wake of leading wave near the bed and moved upward into the water column during the bore propagation. The numerical simulations were based on previous experimental work on undular bores. Some test cases were realised to verify the accuracy of the numerical methods. The results gave access to the detailed flow evolution during the bore propagation. Large velocity reversals were observed close to the no-slip boundaries. In some configurations, coherent turbulent structures appeared against the walls in the wake of the bore front.

Large-eddy simulation (LES) and unsteady Reynolds-averaged Navier-Stokes (URANS) calculations have been performed to investigate the effects of different mathematical models for scalar variance and its dissipation rate as applied to both a non-reacting bluff-body turbulent flow and an extension to a reacting case. In the conserved scalar formalism, the mean value of a thermo-chemical variable is obtained through the PDF-weighted integration of the local description over the conserved scalar, the mixture fraction. The scalar variance, one of the key parameters for the determination of a presumed β-function PDF, is obtained by solving its own transport equation with the unclosed scalar dissipation rate modelled using either an algebraic expression or a transport equation. The proposed approach is first applied to URANS and then extended to LES. Velocity, length and time scales associated with the URANS modelling are determined using the standard two-equation k-ε transport model. In contrast, all three scales required by the LES modelling are based on the Smagorinsky subgrid scale (SGS) algebraic model. The present study proposes a new algebraic and a new transport LES model for the scalar dissipation rate required by the transport equation for scalar variance, with a time scale consistent with the Smagorinsky SGS model.

Submitted by Dilene Paulo (dilene.fatima@ufcg.edu.br) on 2018-02-07T13:41:46Z No. of bitstreams: 1 VERA SOLANGE DE OLIVEIRA FARIAS – TESE PPGEP 2011.pdf: 7179434 bytes, checksum: 4a30c9a95f4a089e00fa550fbf1b42b8 (MD5) / Made available in DSpace on 2018-02-07T13:41:46Z (GMT). No. of bitstreams: 1 VERA SOLANGE DE OLIVEIRA FARIAS – TESE PPGEP 2011.pdf: 7179434 bytes, checksum: 4a30c9a95f4a089e00fa550fbf1b42b8 (MD5) Previous issue date: 2011-04-29 / CNPq / Este trabalho apresenta a solução numérica da equação de difusão tridimensional em regime transiente, para um domínio arbitrário. Para atingir os objetivos, a equação de difusão foi discretizada usando coordenadas generalizadas via método dos volumes finitos com uma formulação totalmente implícita, para condições de contorno de equilíbrio e convectiva. Para cada passo no tempo, o sistema de equações obtido para uma dada malha estruturada foi resolvido pelo método de Gauss-Seidel. O código computacional foi desenvolvido em FORTRAN, usando o estúdio CVF 6.6.0, na plataforma Windows Vista. A solução proposta foi validada usando soluções analíticas e numéricas da equação de difusão para várias geometrias, permitindo validar malhas ortogonais e não-ortogonais. A análise e comparação dos resultados mostraram que a solução proposta forneceu resultados coerentes para todos os casos investigados. O código computacional desenvolvido foi aplicado na simulação, a partir de dados experimentais da secagem de telhas cerâmicas para as seguintes condições experimentais: temperaturas de 55,6 °C, 69,7 °C, 82,7 °C e 98,6 °C e teor de umidade inicial variando de 0,2345 até 0,2405 (b.s.). A simulação tornou possível determinar o coeficiente de difusão efetivo em função da razão de umidade e da temperatura do ar de secagem e também o valor do coeficiente de transferência convectivo de massa correspondente para cada temperatura. / This work presents a three-dimensional numerical solution for the diffusion equation in transient state, in an arbitrary domain. The diffusion equation was discretized using the finite volume method with a fully implicit formulation and generalized coordinates, for the equilibrium and convective boundary condition. For each time step, the system of equations obtained for a given structured mesh was solved by the Gauss-Seidel method. A computational code in FORTRAN, using the CFV 6.6.0 Studio, in a Windows Vista platform was developed. The proposed solution was validated by analytical and numerical solutions of the diffusion equation for several geometries. The geometries tested enabled to validate both orthogonal and non-orthogonal meshes. The analysis and comparison of the results showed that the proposed solution provides correct results for all cases investigated. The developed computational code was applied in the simulation, using experimental data of the drying of ceramic roof tiles, for the following experimental conditions: temperature from 55.6; 69.7; 82.7; 72.8 and 98.7 °C, initial moisture content from 0.2345 up to 0.2405 (d.b.). The simulation makes it possible to determine an expression for the diffusion coefficient as a function of the moisture content and temperature of the drying air, and also the value of the convective mass transfer coefficient corresponding to each temperature.

Les instabilités aéroélastiques de bandes d’acier constituent aujourd’hui l’un des problèmes majeurs dans les sections de refroidissement par jets impactants des lignes de recuit continues.En effet, le traitement thermique des nouveaux aciers nécessite de très fortes pentes de température impliquant constamment des augmentations de vitesse de soufflage susceptibles de mettre en jeu des instabilités aéroélastiques. Des flottements ainsi que des divergences de bande ont déjà été constatées et identifiées. Ces deux instabilités impliquent dans la plupart des cas des chocs entre la bande et les buses de soufflage ce qui engendre des défauts de surface sur la bande.Un banc d’essai a été conçu et fabriqué dans le but d’analyser ces instabilités et d’anticiper leur apparition. A partir d’observations, la dynamique structurelle de la bande a été simplifiée à un mode de rotation rigide. Le banc comporte une plaque oscillante en mouvement forcé.Celle-ci est impactée par un dispositif de plusieurs jets axisymétriques turbulents ayant une disposition identique à celle des tours industrielles. Les efforts aérodynamiques stationnaires et instationnaires agissant sur la plaque sont mesurés au moyen de capteurs de pression.L’impact de plusieurs jets en interaction crée de très importants gradients de pression sur la plaque il est donc nécessaire que la grille de prises de pression soit très fine pour que l’estimation des efforts aérodynamiques soit correcte. La plaque est donc instrumentée de 91capteurs de pression sur une surface de 18 cm². Elle peut également être translatée dans les ois directions de l’espace, ce qui permet d’obtenir la distribution des efforts instationnaires ainsi que des coefficients aéroélastiques sur une grande surface de plaque et à différentes distances d’impact.Les mesures de pression stationnaires ont permis d’établir les courbes d’évolution des efforts d’impact des jets sur la plaque en fonction de la distance jet-plaque ainsi que de la géométrie des buses. Les résultats ont permis de déterminer la stabilité statique de la plaque en mouvement de pompage. Les mesures de vitesses des jets libres ont été effectuées paranémométrie à fil chaud et ont permis de déterminer leurs propriétés statistiques.Les mesures de coefficients aéroélastiques sur la plaque en rotation ont été effectuées surune seule géométrie de soufflage, pour différentes vitesses réduites. Les résultats mettent en évidence l’importance des effets de bords sur la stabilité de plaque. Des méthodes de post traitements ont proposées afin d’extrapoler les résultats à différentes largeurs de bande. Ils sont confrontés aux travaux de Regardin et al. (réf. [1]) et mettent en évidence des désaccords avec le cas réel. Des suggestions sont apportées afin d’améliorer la représentativité du banc vis-à-vis des bandes industrielles. / Aeroelastic instabilities of steel strips impinged by arrays of cooling gas jets have becomeone of the main issues in cooling sections of continuous annealing lines. Indeed, the new steeltreatments require very high temperature variation rates which involve increases in jetvelocities that are likely to onset some aeroelastic instabilities. Strip flutter and divergencehave already been observed and identified. These two aeroelastic instabilities imply a strongrisk of contact with the blowing boxes, which can seriously blemish the strip.An experimental test rig was designed and built in order to analyze and predict of theseinstabilities. From observations, the strip’s structural dynamics was simulated by a rigidrotation mode. The rig included a forced oscillating plate which is impinged by an array ofaxisymmetric jets having the exact industrial geometry. The plate was instrumented withpressure sensors to measure the steady and unsteady surface pressures. The impingement ofinteracting jets creates very large pressure gradients on the plate and therefore a tight mesh ofpressure taps (91 over an 18cm² jet impingement surface) was necessary to allow a goodestimation of the aerodynamic loads The plate could also be moved in the three coordinatedirections as to obtain surface mappings of the unsteady jet forces and aeroelastic coefficientscan be obtained over a wide area and different jet-to-plate distances.The variation of the impinging aerodynamic forces was established as a function of the jetto-plate distance for different nozzle geometries. These results were used to determine the jetstatic stability in plunging motion. Velocity and turbulence measurements in free jets werecarried out using hot wire anemometry in order to determine their statistical properties.Aeroelastic coefficient measurements were carried out on the oscillating plate with onlyone nozzle geometry and different reduced velocities. Results show that the plate’s stability ismainly dependent on the boundary effects. Post processing methods are suggested in order toapply the results to larger plates. Results are compared to the data of Regardin et al. (réf. [1])and emphasize some discrepancies with respect to the real case. Finally some improvementsto the test-rig are suggested for it to be more representative of the industrial situation.

Cette thèse vise à une meilleure compréhension de la dynamique du voilier et à la validation des outils numériques de prédiction de performances et d’optimisation par l'étude expérimentale in situ du problème aéro-élastique d'un gréement. Une instrumentation est développée sur un voilier de 8m de type J80 pour la mesure dynamique des efforts dans le gréement, de la forme des voiles en navigation, du vent et des attitudes du bateau. Un effort particulier est apporté à la mesure des caractéristiques géométriques et mécaniques des éléments du gréement, la calibration des capteurs et au système d'acquisition des données. Les principaux résultats montrent que le voilier instrumenté est un outil adapté pour les mesures instationnaires et soulignent l'amplitude de variation d'effort rencontrée en mer (20 à 50% de l'effort moyen dans une houle modérée). En outre, les variations du signal d'effort sont déphasées avec l'angle d'assiette, créant un phénomène d'hystérésis. Le comportement dynamique d'un voilier en mouvement diffère ainsi de l'approche quasi-statique. Les simulations numériques proviennent du code ARAVANTI, couplage implicite d’un code structure éléments finis ARA et d’un code fluide parfait, limitant son domaine de validité aux allures de près Les résultats de simulation sont très proches des cas stationnaires et concordent bien avec les mesures en instationnaire dans une houle de face. L'expérimentation numérique d'un gréement soumis à des oscillations harmoniques en tangage souligne l'importance de l'approche Interaction Fluide Structure (IFS) et montre que l’énergie échangée par le système avec la houle est reliée à la fréquence réduite et l'amplitude du mouvement. Certaines informations n'étant pas disponibles sur le voilier instrumenté, une expérience contrôlée en laboratoire est développée. Elle consiste en un carré de tissu tenu par deux lattes en oscillation forcée. Les mesures sur cette « voile oscillante » permettent d'étudier les phénomènes IFS avec décollement et sont utilisées pour la validation du couplage ARA-ISIS entre un code fluide Navier-Stokes (RANS) et le même code structure. / This work presents a full scale experimental study on the aero-elastic wind/sails/rig interaction in real navigation condition with the aim to give a reliable database of unsteady measurement. This database is used for the investigation of the dynamic behavior and loads in the rigging and for an experimental validation of an unsteady Fluid Structure Interaction (FSI) model. An inboard instrumentation system has been developed on a 8 meter yacht (J80 class) to simultaneously and dynamically measure the navigation parameters, yacht's motion, sails flying shape, wind and loads in the rigging. A special effort is made on mechanical and geometrical characteristics measurement, sensors calibration and data acquisition system synchronization. Results show that the instrumented boat is a reliable tool to measure the unsteady phenomena in navigation. Dynamic measurements at sea underline the load variation encountered, which represent 20 to 50% of the mean value in a moderate sea state. Oscillations of loads exhibit phase shift with the trim angle, reason for an hysteresis phenomenon, which shows that the dynamic behavior of a sail plan subject to yacht motion clearly deviates from the quasi-steady theory. Simulations are made with ARAVANTI, an implicit coupling of a Finite Element Method structural model ARA and an inviscid fluid model which restricts the simulation domain to upwind conditions. The simulation results compare very well with the experimental data for steady sailing conditions and show a good agreement in unsteady conditions (head swell). Numerical investigation of a sail plan submitted to harmonic pitching motion underlines the importance of FSI modeling and shows that the energy exchanged by the system with the swell increases with the motion reduced frequency and amplitude. Some information is not accessible on the instrumented boat and requires developing a controlled test case in laboratory. The experiment consists of a spinnaker fabric square mounted on two carbon battens moved in forced oscillation. This test case is used to study FSI phenomena with a separated flow and gives experimental results for the validation of the coupling ARA-ISIS of a RANS fluid model with the same structure model.

O tratamento de águas residuárias industriais, submetido constantemente a cargas de choque, é pouco estudado no Brasil, sendo limitado o grau de conhecimento que se tem no mundo sobre os mecanismos de remoção de disruptores endócrinos nos sistemas de tratamento físico-químicos e biológicos. O presente trabalho traz um estudo de caso, envolvendo uma indústria com uma composição de produtos muito diversificada, com despejos líquidos de unidades isoladas de qualidades diversas e que tem como uma das principais matérias-primas, nonil fenóis, conhecidos disruptores endócrinos. Os mecanismos de remoção de carga orgânica e a avaliação da toxicidade da água residuária ao processo biológico de tratamento foram estudados em escala de laboratório, enquanto que a sua tratabilidade por coagulação, floculação, flotação com ar dissolvido seguido de sistema de lodos ativados, em escala piloto. Avaliou-se, também, a aplicabilidade do tratamento biológico com carvão ativado em pó. A grande variação qualitativa e quantitativa da água residuária bruta requer uma unidade de equalização com tempo de detenção de pelo menos 30 horas. O maior responsável por esta variação de qualidade é o processo da unidade química, cujo efluente é proveniente de lavagens de tanques e reatores. Apesar da reduzida eficiência de remoção de DQO obtida com a coagulação, floculação e flotação com ar dissolvido (20 a 30%), constatou-se que este tratamento é essencial para a redução da toxicidade ao processo biológico. A dosagem de coagulante e o pH ótimos variam conforme a característica do despejo e devem ser determinadas diariamente. No sistema biológico, observou-se uma remoção significativa da DQO do efluente bruto gerado na indústria, não apenas por biodegradação, mas também por volatilização e adsorção. Os testes de bancada evidenciaram remoções de DQO por arraste com ar de 47 a 77 % e por adsorção no lodo biológico, de 42%. Apesar dos constantes choques de carga orgânica e de poluentes tóxicos, conseguiu-se atingir o estado estacionário com variações máximas da concentração de sólidos em suspensão voláteis no tanque de aeração entre 20 e 30%. Neste período, a idade do lodo situou-se em torno de 25 a 30 dias e o tempo de detenção hidráulico foi de 3,8 dias. A eficiência média de remoção de DQO neste período foi de 86%. Os testes com a unidade piloto demonstraram que a utilização de carvão ativado em pó (CAP) produz resultados satisfatórios, comprovando a redução de poluentes tóxicos da água residuária e refletindo em uma significativa melhora na biodegradabilidade do efluente, com aumento da concentração de sólidos em suspensão voláteis no tanque de aeração (de 1380 mg/L para 3820 mg/L) e redução da amplitude de variação da DQO do efluente tratado que passou de 600 a 3200 mgO2/L para o sistema sem CAP para 300 a 600 mgO2/L para o sistema com CAP. Notou-se, também, com a adição de CAP, melhoria na sedimentação do lodo. Pode-se concluir também que o sistema de tratamento atendeu a legislação atual do Estado de São Paulo, com uma remoção média de 80% da DBO5. / Wastewater treatment systems continuously receiving shock loads and the behavior of contaminants under unsteady state conditions are not very well documented in our country, with limited degree of knowledge in the world on the removal of specific pollutants like endocrine disrupting chemicals (EDC) in physical chemical and biological wastewater treatment systems. The present work brings a case study regarding an industry which wastewaters were generated from plants with multiproducts and campaign production with variable composition, and with nonylphenols as one of their raw materials which is known as an endocrine disrupting chemical. The organic load and toxicity removal mechanisms of the wastewater to the biological treatment were studied in bench scale whilst their treatability was evaluated through coagulation, flocculation, dissolved air flotation followed by activated sludge, in pilot scale. The feasibility of the addition of powdered activated carbon to the biological system was also studied. The large qualitative and quantitative variability of the wastewater requires an equalization time of at least 30 hours. The main responsible for this variability in quality is the chemical unit process which wastewaters are originated from reactors and tanks cleanings. Although the low efficiency in COD removal obtained with coagulation, flocculation and dissolved air flotation (20 to 30%), it was found that this process is essential to the toxicity reduction for the biological process. Optimum pH and coagulant dose vary with the wastewater characteristics, and must be determined on a daily basis. On the biological system, it was observed important raw wastewater COD elimination not only through biodegradation, but also through volatilization and adsorption. Bench tests revealed COD elimination of 47 to 77% by air stripping and 42% by adsorption onto the biological sludge. Despite of constant organic and toxic shock loads, it had been possible to reach the steady state during which the maximum variation on the volatile suspended solids concentration was 20 to 30%. During this period, sludge age was around 25 to 30 days, and the hydraulic detention time was 3,8 days. Average efficiency on COD removal in this period was 86%. Tests results with PAC dosage on the pilot plant showed satisfactory, proving the reduction of toxic compounds from the wastewater and resulting in biodegradability improvement. The increase of volatile suspended solids in the aeration tank was from 1380 mg/L to 3820 mg/L, and reduction in the range of variation of remaining treated water COD from 600 to 3200 mgO2/L (system without CAP) down to 300 to 600 mgO2/L (system with CAP). The improvement on the sludge sedimentation with PAC addition was also remarkable. One can also conclude the treatment system attained the current legislation of the State of São Paulo, with 80% BOD5 removal.

Orientador: Eugênio Spanó Rosa / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecânica / Made available in DSpace on 2018-08-23T04:15:03Z (GMT). No. of bitstreams: 1 Souza_MarcoAntonioSampaioFerrazde_D.pdf: 6867490 bytes, checksum: 80a84e4d200eeba5b304454e706824a5 (MD5) Previous issue date: 2013 / Resumo: A fração de vazio é um parâmetro essencial para a descrição da estrutura do escoamento intermitente gás-líquido. Na produção de petróleo, por exemplo, em sistemas de elevação artificial, a técnica de gas-lift é muito utilizada em poços verticais para aumentar a produção e requer a previsão dos fenômenos que influenciam o escoamento intermitente para sua otimização. Estes efeitos incluem principalmente as distribuições radial e axial da fração de vazio. Neste trabalho medidas da fração de vazio no pistão de líquido e na bolha alongada em escoamento bifásico vertical ascendente de ar-água no padrão intermitente em golfadas foram realizadas em diversas razões de vazões de gás e líquido utilizando uma sonda elétrica e um sensor de impedância de anéis. Com um conjunto formado por dois sensores de impedância de anéis gêmeos espaçados entre si foi monitorada a passagem dos pistões de líquido e das bolhas alongadas e foram determinados os comprimentos dos pistões e das bolhas assim como a velocidade de translação do nariz da bolha alongada em uma tubulação de 26 mm de diâmetro. As técnicas de medidas foram validadas resolvendo as equações resultantes do balanço volumétrico na célula unitária com os valores das frações de vazio do pistão de líquido, da bolha alongada e da unidade e das velocidades, comprimentos e frequências obtidos experimentalmente. Esta validação tem por objetivo mostrar a consistência que as medidas experimentais apresentam com o modelo físico conferindo às medidas um alto grau de confiabilidade. Os resultados das frações de vazio do pistão de líquido, da bolha alongada e da unidade pistão-bolha foram comparados com correlações disponíveis na literatura assim como os resultados das frações de líquido do pistão de líquido foram comparados e analisados com os modelos fenomenológicos de fração de líquido disponíveis. Com os testes utilizando a sonda elétrica foi possível levantar além da distribuição axial a distribuição radial da fração de vazio local na seção transversal da tubulação. Os resultados obtidos demonstram a consistência das técnicas de medidas e permitem discutir vários fenômenos que influenciam o escoamento bifásico vertical ascendente de ar-água no padrão intermitente em golfadas / Abstract: The void fraction is an essential parameter for describing the structure of the gas-liquid slug flow. In oil production, for example, in artificial lift systems, the gas-lift technique is widely used in vertical wells to increase the production and requires the prediction of the phenomena which influence the slug flow optimization. These effects primarily include axial and radial distributions of void fraction. In this work measures the void fraction of the liquid piston and the elongated bubble in two-phase flow vertical upward air-water in slug flow pattern were performed at various ratios of gas and liquid flow rates using an electric probe and a rings impedance sensor. With a set of two twin rings impedance sensors spaced was monitored passage of liquid slugs and elongated bubbles and were determined and the lengths of the pistons and bubbles as well as the translation velocity of the nose elongated bubble in a pipe 26 mm id. The measurement techniques were validated by solving the equations resulting from the volume balance in unit cell with the values of the void fractions of the piston liquid, the bubble elongated and the unit and velocities, lengths and frequencies obtained experimentally. This validation is to show the consistency with the experimental measurements show the physical model measures conferring a high degree of reliability. The results of the void fractions of the liquid piston, the elongated bubble and bubble-piston unit were compared with correlations available in the literature and the results of the liquid fractions of liquid piston were compared and analyzed with phenomenological models of liquid fraction available. With the tests using the electric probe was possible to rise beyond of the axial distribution the radial distribution of local void fraction in the cross section of the pipe. The results demonstrate the consistency of measurement techniques and allow discuss various phenomena that influence the two-phase flow in vertical upward air-water slug flow pattern / Doutorado / Termica e Fluidos / Doutor em Engenharia Mecânica

Afin d'estimer les forces statiques équivalentes du vent, qui produisent les réponses quasi-statiques et dynamiques extrêmes dans les structures soumises au champ de pression instationnaire induit par les effets du vent, une nouvelle méthode probabiliste est proposée. Cette méthode permet de calculer les forces statiques équivalentes du vent pour les structures avec des écoulements aérodynamiques complexes telles que les toitures de stade, pour lesquelles le champ de pression n'est pas gaussien et pour lesquelles la réponse dynamique de la structure ne peut être simplement décrite en utilisant uniquement les premiers modes élastiques (mais nécessitent une bonne représentation des réponses quasi-statiques). Généralement, les mesures en soufflerie du champ de pression instationnaire appliqué à une structure dont la géométrie est complexe ne suffisent pas pour construire une estimation statistiquement convergée des valeurs extrêmes des réponses dynamiques de la structure. Une telle convergence est nécessaire pour l'estimation des forces statiques équivalentes afin de reproduire les réponses dynamiques extrêmes induites par les effets du vent en tenant compte de la non-gaussianité du champ de pression aléatoire instationnaire. Dans ce travail, (1) un générateur de réalisation du champ de pression instationnaire non gaussien est construit en utilisant les réalisations qui sont mesurées dans la soufflerie à couche limite turbulente; ce générateur basé sur une représentation en chaos polynomiaux permet de construire un grand nombre de réalisations indépendantes afin d'obtenir la convergence des statistiques des valeurs extrêmes des réponses dynamiques, (2) un modèle d'ordre réduit avec des termes d'accélération quasi-statique est construit et permet d'accélérer la convergence des réponses dynamiques de la structure en n'utilisant qu'un petit nombre de modes élastiques, (3) une nouvelle méthode probabiliste est proposée pour estimer les forces statiques équivalentes induites par les effets du vent sur des structures complexes décrites par des modèles éléments finis, en préservant le caractère non gaussien et sans introduire le concept d'enveloppes des réponses. L'approche proposée est validée expérimentalement avec une application relativement simple et elle est ensuite appliquée à une structure de toiture de stade pour laquelle des mesures expérimentales de pressions instationnaires ont été effectuées dans la soufflerie à couche limite turbulente / In order to estimate the equivalent static wind loads, which produce the extreme quasi-static and dynamical responses of structures submitted to random unsteady pressure field induced by the wind effects, a new probabilistic method is proposed. This method allows for computing the equivalent static wind loads for structures with complex aerodynamic flows such as stadium roofs, for which the pressure field is non-Gaussian, and for which the dynamical response of the structure cannot simply be described by using only the first elastic modes (but require a good representation of the quasi-static responses). Usually, the wind tunnel measurements of the unsteady pressure field applied to a structure with complex geometry are not sufficient for constructing a statistically converged estimation of the extreme values of the dynamical responses. Such a convergence is necessary for the estimation of the equivalent static loads in order to reproduce the extreme dynamical responses induced by the wind effects taking into account the non-Gaussianity of the random unsteady pressure field. In this work, (1) a generator of realizations of the non-Gaussian unsteady pressure field is constructed by using the realizations that are measured in the boundary layer wind tunnel; this generator based on a polynomial chaos representation allows for generating a large number of independent realizations in order to obtain the convergence of the extreme value statistics of the dynamical responses, (2) a reduced-order model with quasi-static acceleration terms is constructed, which allows for accelerating the convergence of the structural dynamical responses by using only a small number of elastic modes of the structure, (3) a novel probabilistic method is proposed for estimating the equivalent static wind loads induced by the wind effects on complex structures that are described by finite element models, preserving the non-Gaussian property and without introducing the concept of responses envelopes. The proposed approach is experimentally validated with a relatively simple application and is then applied to a stadium roof structure for which experimental measurements of unsteady pressures have been performed in boundary layer wind tunnel

Dans le contexte économique et environnemental actuel, la prochaine génération de moteurs d’avion devra offrir opérabilité, compacité et hauts rendements. Les compresseurs demeurent une des pièces critiques de ces moteurs, et leur conception un challenge. À débit réduit, leur plage de fonctionnement est contrainte par la limite de pompage, phénomène hautement instable et dangereux. À ce jour, peu d’études expérimentales sur un compresseur en situation de pompage ont été réalisées, étant donné le danger inhérent pour les installations. Dans ce cadre, la simulation numérique peut apporter des informations sur le développement des instabilités aérodynamiques et aider à la prévision de la limite de pompage. L’objectif du travail présenté dans cette thèse est de mettre en place une méthode afin de simuler numériquement l’entrée en pompage et un cycle complet de l’instabilité avec le code elsA. Le cas test retenu est le compresseur de recherche axial multi-étage CREATE dessiné par Snecma, et étudié expérimentalement par le LMFA. Des études antérieures ont montré le rôle joué par les volumes entourant le compresseur ; l’originalité de cette étude réside donc dans l’inclusion des volumes du banc d’essai dans la simulation du compresseur. Une des difficultés inhérentes à la simulation de ces instabilités est leur temps caractéristique, qui représente plus d’une centaine de rotations de la machine. Le calcul a donc nécessité le recours à une approche massivement parallèle ; environ un million d’heures CPU ont été utilisées pour décrire le cycle. Enfin, compte tenu du retournement de l’écoulement dans le compresseur, les conditions aux limites ont été modifiées pour pouvoir s’adapter aux changements de sens de l’écoulement. La simulation a permis de décrire l’entrée en pompage et un cycle complet de l’instabilité. La comparaison avec les données expérimentales montre que les caractéristiques du cycle sont correctement prédites (phénomènes physiques précurseurs de l’instabilité, durée du cycle..). En parallèle, une étude acoustique a été menée afin de mettre en évidence les modes propres du banc d’essai. L’analyse de ces résultats a notamment montré le rôle de l’acoustique dans le déclenchement du pompage. Les différentes phases du cycle de pompage sont ensuite étudiées, et caractérisées (déclenchement, débit inversé, récupération et recompression). Ce travail a généré une base de données qui permet de mieux comprendre les instabilités qui se développent dans ce type de machine. À terme, ces résultats pourront être utilisés pour élaborer et valider des modélisations du phénomène de pompage moins coûteuses, pouvant intervenir dans un cycle de conception. / In order to deal with the current economical and environmental context, the next engine generation will need to offer great operability, compactness and high efficiency. In aircraft engines, the compressor remains one of the critical components, and its design is still a challenging task. At low massflow rate, their operability is bounded by the surge limit, surge being a highly unstable and dangerous phenomenon. Today, few experimental studies on compressor surge are available because of the inherent threat to the facility. In that context, numerical simulation can bring about information on the onset of aerodynamic instabilities and help to predict the surge limit. The work presented in this PhD thesis aims at setting up a method to perform the numerical simulation of surge inception and of an entire cycle of the instability with the CFD code elsA. The chosen test case is the axial multistage research compressor CREATE designed and built by Snecma, and experimentally studied at LMFA. Previous studies have pointed out the role of the volumes adjacent to the compressor ; the originality of this work is thus the inclusion of the volumes of the test-rig in the simulation of the compressor. One of the difficulties inherent to the simulation of those instabilities is their characteristic time of at least one hundred revolutions of the machine. Hence the computation has required a massively parallel approach and about one million CPU hours. Finally, given that the flow reverses during a surge cycle, the boundary conditions have been modified to be able to cope with the flow inversions. The simulation was able to capture surge inception and the entire cycle of the instability. The comparison with the experimental data showed that the main patterns of the cycle are correctly predicted (precursor phenomena of surge, duration of the cycle...). In the meantime, an acoustic study has been performed in order to isolate the eigenmodes of the test-rig. The analysis of the results pointed out the role of acoustic phenomena in surge inception. The different phases of the cycle are then studied and characterized (surge inception, reversed-flow phase, recovery and repressurization). This work has incremented a database that allows a better understanding of the instabilities that develop in this kind of machine. From now on, those results may help to elaborate and validate cheaper models of the surge phenomenon to be used in the design process.