Les travaux de cette thèse se sont intéressés à l’obtention du comportement d’un composite sergé 2x2 verre/PA66 via un outil numérique basé sur une double homogénéisation : la première concerne les torons et la seconde le volume élémentaire représentatif du composite. A partir d’une campagne de caractérisation expérimentale sur le PA66, nous avons dans un premier temps identifié l’influence à la fois de l’hygrométrie et de la température sur le comportement de la matrice en quasi-statique. Ensuite, nous avons élaboré, implémenté et validé une loi de comportement isotrope élasto-plastique endommageable. Cette loi a servi à identifier par homogénéisation le comportement des torons en supposant un collage parfait des fibres et de la matrice ainsi qu’une homogénéité de l’eau dans le toron. Grâce à ces comportements identifiés et à l’élaboration d’une loi de comportement anisotrope élasto-plastique endommageable, nous avons pu déterminer les caractéristiques élastiques dans un premier temps et les comportements longitudinaux et en cisaillement dans un second temps. Les comparaisons aux résultats expérimentaux menés sur le composite offrent des résultats satisfaisants et permettent de penser que l’outil numérique développé, permettrait, à termes, d’aider à une conception rapide incluant ce genre de matériau en diminuant le nombre d’essais expérimentaux à faire. Des perspectives enfin sont proposées notamment pour l’extension aux comportements en dynamique (crash). / The objective of this PhD thesis is to create a methodology to simulate the behavior of a 2x2 glass/PA66 twill composite using homogenization. To achieve this goal, we use two steps: first one is dedicated to represent the yarn behavior and the second one to represent the behavior of the composite RVE. An experimental characterization campaign about the PA66 has been done in order to quantify the influence of both hygrometry and temperature on the behavior of the matrix in quasi-static. Then, we formulated, implemented and validated an isotropic elasto-plastic damageable constitutive law. It was used to identify by homogenization the behavior of the yarns by assuming a perfect bonding of the fibers and the matrix as well as a water homogeneity inside the yarn. Thanks to the identified behavior of the yarns and by formulating a compatible constitutive law (anisotropic elasto-plastic damageable), we have determined for the woven composite on one hand the elastic properties and in the other hand the longitudinal/transversal and shear nonlinear behaviors. The comparisons of the numerical and experimental results offer satisfactory results. It provides good outlook in short-term in regard of structures design which include this type of woven composite materials. The main advantage of the methodology is to decrease the number of required experimental tests. Some outlook dedicated to crash studies are finally suggested.

Certaines applications industrielles nécessitent des niveaux d’étanchéité exceptionnels pour permettre la réalisation d’un vide poussé ou pour répondre à des enjeux de sécurité radiologique par exemple. Ces niveaux de haute étanchéité statique sur des assemblages démontables sont obtenus à l’aide de joints entièrement métalliques. La fuite résultante de l’assemblage n’est due qu’à la persistance d’un champ des ouvertures à l’interface entre le joint d’étanchéité et la bride d’assemblage, conséquence d’un contact imparfait entre les deux surfaces rugueuses. Le champ des ouvertures à l’interface de contact est assimilable à une fracture rugueuse hétérogène, de nature multi-échelle, et peut en principe être obtenu par un calcul de déformations mécaniques préalable. Dans ce travail, on s’intéressera plus particulièrement à l’écoulement gazeux raréfié dans le régime glissant au sein de cette fracture. Pour les régimes modérément raréfiés,l’écoulement est modélisé par l’équation de Reynolds faiblement compressible avec correction de glissement de premier ordre aux parois que l’on développe. On effectue ensuite un changement d’échelle par la méthode de la prise de moyenne volumique, permettant d’établir un modèle macroscopique d’écoulement à l’échelle d’un élément représentatif, où le débit massique est relié au gradient de pression via le tenseur de transmissivité. Celui-ci, caractéristique de l’élément représentatif de fracture, est calculé par résolution d’un problème de fermeture et est dépendant de la microstructure ainsi que du libre parcours moyen représentatif sur l’élément. Pour remonter à l’écoulement dans l’ensemble de la fracture, hétérogène à cette échelle, celle-ci est subdiviséeen pavés sur chacun desquels est calculé un tenseur de transmissivité local par la méthode sus-citée. Ensuite, l’écoulement dans ce champ de tenseurs est résolu par une méthode des éléments finis de frontière, donnant la transmissivité apparente glissante du joint dans son ensemble. Cette approche à deux échelles, vue comme outil d’aide à la conception, permet une réduction de la complexité de calcul par rapport à une simulation directe, rendant possible une analyse plus efficace du comportement d’un système d’étanchéité. Pour valider l’utilisation du modèle de glissement d’un point de vue macroscopique et s’affranchir des incertitudes sur le calcul de déformation mécanique, des puces nanofluidiques de type réseau hétérogène de canaux droits sont fabriquées par photolithographie par niveaux de gris. Des essais expérimentaux de mesure de fuite sont réalisés sur ces géométries modèles, représentant des joints idéalisés. Ces essais sont effectués en appliquant une forte différence de pression d’hélium par utilisation d’un spectromètre de masse mesurant la fuite, produisant une condition de vide en sortie de puce.Selon les puces, les régimes de raréfaction atteints vont alors du régime glissant au régime moléculaire. Le débit de fuite mesuré est alors supérieur à celui prédit par le modèle de premier ordre, l’écart restant inférieur à un ordre de grandeur quel que soit le régime / Some industrial applications require exceptional sealing levels to maintain ultra-high vacuumconditions or for radiological safety concerns for example. Such high performance static sealingconditions on mechanical assemblies are reached using entirely metallic gaskets. The resultingleak-rate is only due to the persistence of an aperture field at the seal-flange interface,consequence of a non-ideal contact between the two rough surfaces. This aperture field can beviewed as a rough and heterogeneous fracture, of multi-scale nature, and can be obtained by aprior contact mechanics computation. In this work, we are interested on the rarefied flow of a gasin this fracture, drawing our attention to the slip regime. For such moderately rarefied regime, theflow is described by the slightly compressible Reynolds equation with a first-order slip-flowcorrection at the walls, which we develop. Using the method of volume averaging, an upscalingprocedure is performed to derive the macroscopic flow model at the scale of a representativeelement, and where the mass flow rate is related to the pressure gradient by the transmissivitytensor. This latter is characteristic of the representative fracture element and is obtained by solvingan auxiliary closure problem which depends on the micro-structure as well as the representativemean free path on the element. To compute the flow in the whole fracture, heterogeneous at thisscale, it is subdivided in tiles on which a transmissivity tensor is locally computed by theaforementioned method. Then, the flow problem in this tensor field is solved using a boundaryelement method, leading to the apparent slip-corrected transmissivity of the entire aperture field.This two-scale approach is a conception tool which reduces the overall complexity with respect toa direct numerical simulation, allowing a more efficient analysis of the behavior of a sealingassembly. To validate the use of slip models at the macroscopic level and to eliminate theuncertainties of the contact mechanics computation, nanofluidic chips composed ofheterogeneous network of straight channels are fabricated using a grayscale photolithographytechnique. Experimental measurements of the leak-rate are performed on these idealizedgeometries that mimic a seal assembly. They are realized by applying a strong helium pressuredifference on the chip using a mass spectrometer to measure the leak, which produces a nearvacuum condition at the outlet. Depending of the chip, the rarefaction regime ranges from slip tofree-molecular. The measured leak-rate is greater than predicted by the first order model, thoughbeing of the same order of magnitude whatever the regime

Les réseaux de tuyauterie industriels sont le siège de niveaux importants de vibrations induites par l’écoulement qui peuvent mener à la rupture par fatigue des installations. La présente étude se concentre sur l’analyse et modélisation simplifiée de la source vibratoire associée au passage d’un écoulement liquide turbulent par un coude à 90°. Une approche combinant expériences et simulation a été conduite. Une boucle de circulation d’eau munie d’un coude transparent a été conçue afin de permettre des me-sures de vitesse à l’intérieur du coude. Pour cette finalité, les techniques de Vélocimétrie par Imagerie de Particules (PIV) plane et stéréoscopique ont été utilisées. La pression pariétale et les vibrations du banc d’essais ont été mesurées simultanément. Plusieurs configurations d’écoulement ont été testées afin d’obtenir une riche base de données couplées reliant l’écoulement fluide à l’excitation dynamique des parois et, finalement, à la réponse vibratoire de la structure. En parallèle, l’écoulement instationnaire d’eau dans le coude a été simulé au moyen d’une approche du type Simulation des Grandes Echelles (LES). La simulation fluide a permis d’étudier en détails la topologie de l’écoulement turbulent au passage du coude ainsi que le champ instationnaire de pression fluctuante induit sur la paroi. Finale-ment, un ensemble d’outils statistiques a été appliqué aux données expérimentales et numériques afin de proposer un modèle simplifié des transferts qui relient l’écoulement turbulent à la sollicitation dynamique de la structure contenant le coude. / The flow of dense fluids within thin-walled piping systems may lead to significant levels of Flow-Induced Vibration, mainly in the vicinities of singularities such as obstacles inserted into the flow, sudden changes of cross-sectional area or flow direction. This study focuses on the analysis and reduced-order modelling of the vibrational source associated with the turbulent flow of liquids through a 90° elbow.A mixed experimental-computational approach is undertaken. A closed water loop containing a transparent elbow was designed in order to allow for fluid velocity measurements inside the singularity. To this purpose, planar and stereoscopic Particle Image Velocimetry (PIV) were employed. Wall pressure fluctuations and structural vibrations were measured simultaneously. Several flow configurations were tested in order to obtain a large coupled database linking the flow to the dynamic excitation, and then to the vibration response of the structure.In parallel, the unsteady water flow through the elbow was computed using Large-Eddy Simulation (LES). The fluid simulation allowed for a detailed study of the turbulent flow through the singularity and the unsteady pressure field induced on the piping walls. Finally, a set of statistical tools was applied to both experimental and computational data in order to propose a reduced-order model of the transfer function that links the tur-bulent flow to the dynamic excitation of the elbowed piping structure.

Le présent travail de thèse se consacre au développement d'un nouveau modèle micromécanique pour les composites en thermoplastique renforcé par des fibres de verre courtes. L'objectif est notamment la modélisation du comportement visco-endommageable en fatigue du PA66-GF30. Ce matériau, particulièrement utilisé dans l'industrie automobile, est sujet à une microstructure spécifique issue du procédé de moulage par injection. L'approche multi-échelles développée consiste en une méthode de Mori-Tanaka modifiée, appliquée à des renforts avec enrobage et prenant en compte l'évolution de l'endommagement à l'échelle microscopique. La description des mécanismes d'endommagement se base sur une investigation expérimentale poussée préalablement menée au sein de l'équipe. Des scénarios d'endommagement ont été proposés et incluent trois processus locaux différents : la décohésion de l'interface, la microfissuration de la matrice et les ruptures de fibres. Ceux-ci sont spécialement affectés par la microstructure. L'approche développée intègre ces cinétiques d'endommagement ainsi que la viscoélasticité non-linéaire de la matrice et la distribution d'orientation des inclusions due au procédé de fabrication. Chaque mécanisme d'endommagement est modélisé par une loi d'évolution basée sur les contraintes locales calculées à l'échelle microscopique. La loi constitutive finale, à l'échelle du volume élémentaire représentatif, est implémentée dans une bibliothèque scientifique en C++, SMART+, et est conçue pour être compatible avec une analyse de structures par éléments finis. L'identification du modèle est réalisée par rétro-ingénierie, en tirant profit de résultats expérimentaux multi-échelles, dont notamment des tests in-situ au MEB ainsi qu'une analyse qualitative et quantitative par μCT. / The current work focuses on a new micromechanical high cycle fatigue visco-damage model for short glass fiber reinforced thermoplastic composites, namely: PA66/GF30. This material, extensively used for automotive applications, has a specific microstructure which is induced by the injection process. The multi-scale developed approach is a modified Mori-Tanaka method that includes coated reinforcements and the evolution of micro-scale damage processes. Their description is based on the experimental investigations of damage mechanisms previously performed by the team. Damage chronologies have been proposed involving three different local degradation processes: fiber-matrix interface debonding/coating degradation, matrix microcracking and fiber breakage. Their occurrence strongly depends on the microstructure. The developed model integrates these damage kinetics and accounts for the complex matrix viscoelasticity and the reinforcement orientation distributions induced by the process. Each damage mechanism is introduced through an evolution law involving local stress fields computed at the microscale. The developed constitutive law at the representative volume element scale is implemented into a C++ scientific library, SMART+, and is designed to work with Finite Element Methods. The model identification is performed via reverse engineering, taking advantage of the multiscale experimental results: in-situ SEM tests as well as quantitative and qualitative μCT investigations.

La surveillance de la qualité de l'air est actuellement effectuée avec des mesures de concentration et à partir d'outils de modélisation de la dispersion atmosphérique. Ces modèles numériques évaluent les concentrations des polluants avec une résolution spatio-temporelle plus fine que les mesures. Néanmoins, les estimations fournies par ces modèles sont moins précises que les mesures. Dans ce projet de recherche, nous avons étudié les approches de couplage d'échelles et d'assimilation de données pour améliorer les estimations fournies par le modèle de dispersion atmosphérique SIRANE, dédié à l'échelle urbaine. L'approche de couplage d'échelles consiste à déterminer les conditions aux limites d'une simulation à partir d'une autre simulation à plus grande échelle. Au cours de ce travail de thèse, nous avons analysé trois méthodes afin de coupler le modèle urbain SIRANE et le modèle à méso-échelle CHIMERE. Cette étude montre que ces méthodes permettent potentiellement d'estimer la qualité de l'air à l'échelle urbaine de manière plus satisfaisante que les modèles à méso-échelle (utilisés seuls). Cependant, elles n'améliorent pas forcément la modélisation des conditions aux limites d'une simulation à l'échelle urbaine et les estimations fournies par celles-ci. Cela est a priori lié au fait que les estimations fournies par le modèle CHIMERE ne sont pas suffisamment satisfaisantes sur notre cas d'étude. Il est néanmoins possible que ces méthodes améliorent les résultats à l'échelle urbaine en utilisant une simulation à l'échelle régionale de meilleure qualité. L'approche d'assimilation de données consiste à combiner les mesures et les données modélisées afin de déterminer la meilleure estimation de l'état d'un système. Durant cette thèse, nous avons étudié trois méthodes d'assimilation de données : la méthode de débiaisement, la méthode que nous avons nommée modulation de la contribution des sources et la méthode Best Linear Unbiased Estimator. Cette étude indique que ces méthodes permettent globalement d'améliorer les estimations fournies par le modèle SIRANE. L'étude de sensibilité vis-à-vis du nombre de mesures utilisées lors de l'assimilation de données indique qu'en général, plus ce nombre est élevé plus les résultats sont satisfaisants. Enfin, les résultats montrent que les performances statistiques associées à ces trois méthodes d'assimilation de données sont globalement comparables entre elles sur notre cas d'étude. / Air quality monitoring is currently carried out with concentration measurements and with atmospheric dispersion modeling tools. These numerical models evaluate pollutant concentrations with a finer spatio-temporal resolution than measurements. Nevertheless, the estimates provided by these models are less accurate than measurements. In this research project, we studied multiscale coupling and data assimilation approaches to improve the estimates provided by the SIRANE atmospheric dispersion model, dedicated to the urban scale. The multiscale coupling approach consists in determining the boundary conditions of a simulation from another simulation on a larger scale. In this thesis work, we analyzed three methods for coupling the SIRANE model with the CHIMERE mesoscale model. This study shows that these methods can potentially estimate the air quality at the urban scale more satisfactorily than the mesoscale models (used alone). However, they do not necessarily improve the modeling of the boundary conditions of a simulation at the urban scale and the estimates provided by them. This is a priori due to the fact that the estimates provided by the CHIMERE model are not sufficiently good on our case study. It is possible, however, that these methods improve the results at the urban scale by using a better simulation at the regional scale. The data assimilation approach consists of combining the measurements and the modelled data to determine the best estimate of the system state. During this thesis, we studied three data assimilation methods : the unbiased method, the method that we called source apportionment modulation, and the Best Linear Unbiased Estimator method. This study indicates that these methods generally improve the estimates provided by the SIRANE model. The sensitivity study on the number of measurements used during the data assimilation indicates that, in general, higher is this number, more satisfactory are the results. Finally, the results show that the statistical performances associated with these three data assimilation methods are globally comparable on our case study.

Le tatouage d’images (ou watermarking) consiste à insérer de l’information de copyright, ou d’autres messages de vérification, invisible à l’oeil humain. Parmi toutes les approches possibles, le tatouage psychovisuel est le plus performant. Il consiste à utiliser les connaissances du Système Visuel Humain (SVH) pour en déduire la force du tatouage de manière optimale. Dans cette étude, nous avons proposé une nouvelle méthode de tatouage d'image psychovisuel qui combine les avantages des ondelettes à celle du domaine spatial par l’utilisation d'une représentation en ondelettes à échelle mixée (JNDEM). En considérant la densité des coefficients d'ondelettes dominants, notre système est capable de différencier les zones uniformes, celles incluant des contours, et enfin celles texturées. Nous avons sélectionné les méthodes efficaces de tatouage psychovisuel adaptée à chaque zone, à savoir celle de Chou pour les zones uniformes, de nouveau celle de Chou pour les contours, et finalement celle de Qi pour les textures. Nous avons aussi pris en compte la sensibilité du SVH qui varie avec l’orientation locale de l’image. Nous avons comparé notre JNDEM à quatre autres méthodes de tatouage de la littérature en termes de qualité visuelle et de robustesse. La méthode JNDEM conduit à une grande qualité visuelle mesurée par des tests objectifs et subjectifs. De plus elle est robuste aux attaques. Ce dernier point est d'une grande importance pour les applications réelles. En conclusion, la méthode proposée est le meilleur compromis entre la qualité visuelle et la résistance aux attaques comparée aux quatre autres méthodes testées. / Watermarking consists in inserting copyright information, or other verification messages, invisible to the human eye. Among all the possible approaches, perceptual watermarking is the most efficient. Perceptual image watermarking consists in using knowledge of the Human Visual System (HVS) to choose the strength of the watermark. In this study, we propose a new perceptual image watermarking method (JNDEM ) that combines the advantages of both the wavelet domain and the spatial domain since a mixed-scale wavelet representation is applied. By considering the density of the dominant wavelet coefficients, our scheme is able to differentiate uniform, edge and texture areas. We selected effective Just Noticeable Difference (JND) models from the literature, i.e. luminance and edge masking developed by Chou, and texture masking developed by Qi. We also took into account the HVS sensitivity which varies with the orientation of the image activity. We have compared our JNDEM to four other JND models from the literature in terms of visual quality and of robustness to attacks. Results show that the proposed method yields a high visual quality and is very robust against attacks. This last point is of great importance for real applications. In conclusion, the proposed method is the best trade-off between visual quality and resistance to attacks among the four other tested methods.

Les matériaux nanostructurés suscitent un intérêt qui va croissant en raison de leurs propriétés chimiques et physiquesexceptionnelles. A cause de la complexité et du coût des développements expérimentaux à l’échelle nano, la simulationnumérique devient une alternative de plus en plus populaire aux études expérimentales. Dans ce travail de thèse, nous avons essayé de combiner des simulations à l’échelle atomique avec de la modélisation en milieu continu pour évaluer la conductivité thermique et la réponse élastique de matériaux nanostructurés. Nous avons utilisé des simulations de dynamique moléculaire pour calculer la réponse mécanique et thermique des matériaux sur des volumes à l’échelle nano. Des méthodes de micromécanique et la méthode des éléments finis, qui utilisent la mécanique des milieux continus, ont permis d’évaluer les propriétés mécaniques des matériaux à l'échelle macroscopique. Les résultats obtenus par ces simulations numériques ont été ensuite comparés avec ceux issus de l’expérience. / Nanostructured materials are gaining an ongoing demand because of their exceptional chemical and physical properties. Due to complexities and costs of experimental studies at nanoscale, computer simulations are getting more attractive asexperimental alternatives. In this PhD work, we tried to use combination of atomistic simulations and continuum modeling for the evaluation of thermal conductivity and elastic stiffness of nanostructured materials. We used molecular dynamics simulations to probe and investigate the thermal and mechanical response of materials at nanoscale. The finite element and micromechanics methods that are on the basis of continuum mechanics theories were used to evaluate the bulk properties of materials. The predicted properties are then compared with existing experimental results.

Cette thèse présente une étude de la dynamique des plages sableuses ouvertes dominées par la houle, au travers d’une base de données originale, couvrant une période de 29 mois, et composée de 150 levés DGPS couvrant 750m de linéaire côtier, donnant accès à la morphodynamique du site de Biscarrosse à différentes échelles de temps. Dans un premier temps, l’analyse du jeu de données à l’échelle des événements (tempêtes et successions de tempêtes appelées clusters) nous a permis de montrer que la réponse des plages sableuses aux clusters ne résulte pas de la somme des impacts induits par chaque tempête d’un cluster. Ainsi, l’effet cumulé des clusters, rapporté sur d'autres sites dans la littérature, n’est ici pas vérifié. L'impact de l’enchainement des tempêtes a également été étudié et il en résulte que lors d’un cluster, un changement des conditions hydrodynamique, à savoir, une augmentation des hauteurs de vagues et/ou du niveau d’eau, est nécessaire pour que la tempête suivante ait un impact érosif significatif sur le système. Dans un second temps, nous avons étudié la dynamique saisonnière du système plage/dune, que ce soit la saison hivernale ou estivale, dans le but de mettre en relief les principaux processus impliqués à cette échelle. Nos travaux montrent que la réponse hivernale de la plage ne dépend pas uniquement des conditions énergétiques et du profil pré-hivernal de la plage, mais également du séquençage des événements, comme lors d'un cluster. Mes travaux confirment également la nécessité de prendre en compte de nombreux paramètres dans l’étude de la dynamique hivernale des littoraux sableux : les conditions hydrodynamiques, le séquençage des évènements érosifs mais également reconstructifs, en particulier le ré-engraissement post-évènement, les transports sédimentaires cross-shore et longshore, ainsi que la position de la barre interne et des courants d’arrachements. La saison estivale est, quant-à-elle, marquée par la reconstruction de berme. Elle semble être liée à la fois aux conditions hydrodynamiques et aux caractéristiques des barres sableuses. L’étude de deux étés et deux hivers successifs a ainsi permis d’identifier les interactions entre les saisons et l’impact de la saison hivernale sur l’estivale, et l’influence de la dynamique événementielle sur la dynamique saisonnière. Elle a aussi permis de mettre en relief l’impact de l’urbanisme et des stratégies d’aménagement dans la réponse du système, à différentes échelles de temps. / This thesis presents a study of an open sandy beach wave-dominated, based on an original dataset, covering 29 months and composed by 150 DGPS surveys recorded along 750m of sandy shore, giving an access to the morphodynamic of Biscarrosse beach at different timescales. In a first time, event scale analysis showed that sandy beach response to clusters is not the result of the sum of the impact generated by each storm of a cluster on the system. Thus, the cumulated effect of clusters, described in the literature is not verified here. The storm sequencing has also been studied: during a cluster, changes in hydrodynamics conditions (rising of the water level and/or wave height) are necessary to provoke a significant erosion of the system by the second storm. In a second time, we studied the seasonal scale dynamic of the beach/dune system (winter and summer seasons) with the purpose to highlight dominant processes involved at this timescale. Beach response to winter seasons not only depends on hydrodynamic conditions and previous beach profile, but also on erosion/recovery event sequencing, post-storm recovery, cross-shore and longshore sediment transport, the barline characteristics and RIP current positions. Summer seasons are here defined by the berm reconstruction. Recovery periods are both linked to hydrodynamic conditions and barline characteristics (e.g. position and shape).The study of successive winters and summers allowed us to identify interactions between seasons, and the influence of short-scale dynamics on the seasonal one. It also emphasizes the impact of urbanism and coastal management strategies on the system’s response, at different timescales.

Pour des raisons de sécurité, les moteurs aéronautiques doivent pouvoir redémarrer en vol sur toute leur plage d'opération. Mais les contraintes sur les émissions polluantes nécessitent le développement de nouvelles chambres de combustion dont la conception peut détériorer les capacités d'allumage du moteur. Afin d'améliorer la compréhension du processus d'allumage et d'aider à l'optimisation de la conception, les recherches actuelles combinent les études expérimentales de plus en plus complexes et les simulation numériques hautes fidélités. Dans ce travail, l'étude numérique du processus d'allumage des chambres de combustion aéronautiques, de l'étincelle à la propagation azimutale de la flamme, est conduite avec plusieurs objectifs: améliorer la robustesse et la confiance de l'outil LES pour l'étude de l'allumage, étudier les mécanismes qui affectent l'allumage dans des conditions représentatives des conditions réelles et enfin améliorer les méthodes bas-ordre pour la prédiction des performances d'allumage. Dans une première partie, la SGE d'un monobruleur installé au CORIA permet de mettre en évidence les bons résultats de la LES et de construire une base de données pour l'analyses des mécanismes d'extinction. Ces données sont aussi utilisées pour développer une méthodologie permettant de prédire les performances d'allumage à bas coût en utilisant les résultats d'une SGE non-réactive. Dans une seconde partie, la propagation inter-secteur est investiguée par l'étude de deux cas expérimentaux et la SGE est capable de reproduire les modes de propagation mais aussi les temps d'allumage avec précision. Sur la bases de ces bons résultats, une analyse plus fine de la simulation permet d'identifier les mécanismes qui contrôlent la propagation de la flamme. / For safety reasons, in-flight relight of the engine must be guaranteed over a wide range of operating conditions but the increasing stringency of pollutant emission constraints requires the development of new aero-engine combustor whose design might be detrimental to the ignition capability. To improve the knowledge of the ignition process in aeronautical gas turbines and better combine conflicting technological solutions, current research relies on both complex experimental investigation and high fidelity numerical simulations. In this work, numerical study of the ignition process in gas turbines from the energy deposit to the light-around is performed with several objectives: increase the level of confidence of Large Eddy Simulations tool for the analysis of the ignition process, investigate the mechanisms controlling ignition in conditions representative of realistic aeronautical gas turbine flows and improve the low-order methodologies for the prediction of ignition performance. In a first part, LES of the single burner installed at CORIA (France) is carried out and allows to highlight the LES accuracy and to build a database on which the main mechanisms controlling the ignition success are identified. Based on these results, a methodology is developed to predict the ignition performance at a low computational cost using the non-reacting flow statistics only. In a second part, the light-around process is studied on two experimental set-ups and the very good agreement of the LES results with experiments is the starting point from an analysis of the mechanisms driving the flame propagation process.

Dans les grands systèmes de réfrigération, l’utilisation de réfrigérants comme lubrifiants, à la place des huiles, simplifie la conception du système, l’allège et réduit son impact environnemental. La très faible viscosité du réfrigérant conduit à une épaisseur de film séparant les deux surfaces comparable à leur rugosité. Néanmoins, des travaux expérimentaux avec le réfrigérant R1233zd suggèrent que la lubrification est possible dans ces conditions grâce à la formation d’une couche adsorbée sur la surface d’oxyde de fer. Les analyses expérimentales in situ dans le contact sont très difficiles à cause des conditions de fort confinement et haute pression. C’est pourquoi une approche numérique multi-échelles est développée, afin d’étudier l’impact des réactions physico-chimiques à l’interface réfrigérant-surface sur les propriétés de lubrification. La théorie de la fonctionnelle de la densité est utilisée pour quantifier au niveau quantique, l’adsorption d’une molécule de réfrigérant sur une surface d’oxyde de fer. Des énergies de liaison allant de -0.92 eV à -0.22 eV sont observées et reliées à différents cas d’adsorption. Ces résultats sont exploités pour paramétrer un champ de forces interfacial, qui prédit des structures moléculaires à l’interface, différentes de celles obtenues avec des potentiels basés sur les règles de mélange classiques. Des simulations de dynamique moléculaire utilisant ce champ de forces paramétré confirment l’existence d’une couche fortement adsorbée de R1233zd sur une surface d’oxyde de fer. Avec des surfaces atomiques lisses et seulement 2 nm d’épaisseur de film de réfrigérant, les couches adsorbées résistent à des pressions allant jusqu’à 4 GPa et des vitesses de cisaillement atteignant 100 m/s. Une valeur minimale de 5 molécules de réfrigérant par nm² est nécessaire à la formation de deux couches adsorbées à 0.5 GPa. De plus, des simulations en cisaillement avec une surface rugueuse prédisent une rupture totale du film à environ 13 GPa. / In large refrigeration systems, using the refrigerant as lubricant instead of oil can help to simplify the design, lighten the systems, and reduce their environmental impact. However, the very low viscosity of refrigerants leads to ultra-thin films separating the surfaces, with a thickness comparable to surface roughness. Nevertheless, experiments with the R1233zd refrigerant suggests that lubrication is still possible in that situation thanks to an adsorbed layer formed on iron oxide surfaces. Experimental in situ analysis area is very difficult because of high confinement and high pressure. That is why a multiscale numerical approach is developed here, to explore the impact of chemical reactions and physical processes at the refrigerant-surface interface on large-scale lubrication properties. Density functional theory is used to quantify the adsorption of a refrigerant molecule on an iron oxide surface at the quantum level. Binding energies ranging from -0.92 eV to -0.22 eV are measured and related to different adsorption cases. These results are then used to parametrize an interfacial force field, whose predictions of interfacial molecular structure differs from those obtained using potentials based on standard mixing rules. Large-scale molecular dynamics simulations involving this parametrized force field confirm the existence of a strongly adsorbed layer of R1233zd molecules on iron oxide surface. With atomically smooth surfaces, and a refrigerant film thickness as small as 2 nm, the adsorbed refrigerant layers resists pressures as high as 4 GPa and sliding velocities as high as 100 m/s. A minimum value of 5 refrigerant molecules per nm² is necessary to the formation of two adsorbed layers at 0.5 GPa. Moreover, sliding simulations with a rough surface reveal total film breakdown for ca. 13 GPa.