Cette thèse s'intéresse à la problématique du contact entre matériaux hétérogènes. L'industrie (automobile, aéronautique, spatiale, ...) s'intéresse de plus en plus à ces types de matériaux. Il s'agira par exemple des alliages métalliques, des matériaux poreux, matériaux composites (composites tissés, interlocks 3D, interlocks 2D), des billes céramiques contenant des impuretés (porosités/précipités),... Dans ce manuscrit, un modèle de contact basé sur les méthodes semi-analytiques a été développé. Un algorithme de gradient conjugué est utilisé afin de résoudre rapidement le problème de contact. Le modèle permet de prendre en compte la présence d'une ou de plusieurs hétérogénéités isotropes/anisotropes dans le problème de contact. Une approche inspirée de la méthode de l'inclusion équivalente proposée par Eshelby est utilisée dans le solveur de contact pour prendre en compte l'effet de ces hétérogénéités. Les méthodes de transformées de Fourier rapides (FFT) permettent d'accélérer les calculs. Une méthode numérique a été mise en œuvre afin de prendre en compte l'interaction entre plusieurs hétérogénéités. Le massif peut être élastique ou viscoélastique. L'approche développée dans la thèse peut résoudre à la fois les problèmes d'indentation, de roulement/glissement ou de fretting en présence de matériaux élastiques hétérogènes, viscoélastiques homogènes ou hétérogènes. Les solutions sont données en termes de champs de pressions, de cisaillements et de contraintes. Dans le cas des matériaux viscoélastiques le code de calcul est capable de fournir le coefficient de frottement apparent ainsi que toutes les variables de contact aussi bien en régime permanent que transitoire. Le modèle a été validé par comparaison avec la méthode des éléments finis classiques en utilisant le logiciel commercial Abaqus v6.11. Le temps de calcul ainsi que l'espace mémoire nécessaire sont considérablement réduits par rapport à la méthode éléments finis. La parallélisation a été introduite dans le code de contact afin de réduire toujours plus le temps de calcul. Il s'agit d'un code robuste, rapide et facilement utilisable en Bureau d'Etudes. Une approche expérimentale originale a été mise en place afin de mesurer les champs de déplacements à l'interface des corps en contact. De bonnes corrélations essais/calculs ont été obtenues. Enfin quelques applications industrielles ont été présentées. Un couplage entre un code éléments finis structurel et le code semi-analytique de résolution de contact a été également réalisé. / The present PhD thesis deals with contact problems between heterogeneous materials. Nowadays heterogeneous materials are extensively used in several industrial domains (automotive, aeronautics, aerospace, ...). Heterogeneous materials involve porous materials, aluminum alloys, composites materials (woven composites, interlocks 3D, interlocks 2D), metallic or ceramics materials containing impurities (porosities/precipitates). In this work, a contact model based on semi-analytical method is proposed. A conjugate gradient algorithm is used for a fast resolution of contact equations. The model can account for one or more isotropic/anisotropic inhomogeneities. An approach taking inspiration from the Eshelby equivalent inclusion method is used in the contact solver to account for the effect of inhomogeneities. 2D and 3D Fast Fourier Transforms (FFT) are used to speed up the computation. A numerical method is implemented in order to take into account interactions between many heterogeneities. The semi-infinite space/ matrix can be either elastic or visco-elastic. The model developed in the present PhD thesis can solve indentation, rolling/sliding or fretting contact problems between heterogeneous elastic materials, homogeneous or heterogeneous visco-elactic materials. In the case of visco-elastic materials, the model permits to get the solution in terms of contact pressure distribution, subsurface stresses, apparent friction coefficient, both in the transient and then steady-state regimes. The model has been validated by performing a comparison with the results of a finite element model. The CPU time and memory necessary are greatly reduced in comparison with the classical finite element method. The model developed is fast, robust and extremely easy to use. An original experimental approach was proposed in order to measure the displacement fields at interface of two contacting bodies. A good agreement between experimental results and numerical simulations is obtained. Finally, the model is applied on some industrial applications. A coupling between a finite element model and the semi-analytical code allow to take into account the effects of structure on contact problem.

Cette thèse a pour but de faire le lien entre l’évolution des galaxies, leur morphologie et les processus physiques internes, notamment la formation stellaire comme le résultat du milieu interstellaire turbulent et multiphase, en utilisant les simulations cosmologiques zoom-in, les simulations des galaxies isolées et en interaction, et le modèle analytique de la formation stellaire. Dans le chapitre 1, j’explique la motivation pour cette thèse et je passe brièvement en revue le contexte nécessaire lié à la formation des galaxies et la modélisation en utilisant les simulations numériques. Tout d’abord, j’explore l’évolution de la morphologie des galaxies du type de la Voie Lactée dans la série des simulations cosmologiques zoom-in à travers l’analyse des barres. J’analyse l’évolution de la fraction des barres avec le redshift, sa dépendance en fonction de la masse stellaire et l’histoire d’accrétion de galaxies individuelles. Je montre en particulier, que la fraction de barres décroit avec le redshift croissant, en accord avec les observations. Ce travail montre également que les résultats obtenus suggèrent que l’époque de la formation des barres correspond à la transition entre une phase précoce “violente” de la formation de galaxies spirales à z > 1, pendant laquelle elles sont souvent perturbées par les fusions avec les galaxies de masse comparable ou par multiple fusions avec les galaxies de petite masse, mais aussi les instabilités violentes de disque, et une phase "séculaire" tardive à z < 1, quand la morphologie finale est généralement stabilisée vers une structure dominée par le disque. Cette analyse est présentée dans le chapitre 2. Étant donné que ces simulations cosmologiques forment trop d'étoiles trop tôt par rapport aux populations de galaxies observées, je me concentre dans le chapitre 3 sur la formation stellaire dans un échantillon de simulation de galaxies en isolation, à bas redshift, et à résolution du parsec et sous-parsec. J'étudie l'origine physique de leurs relations de formation stellaire avec les cassures, et montre que le seuil de densité surfacique pour une formation stellaire efficace peut être lié à la densité caractéristique d'apparition de turbulence supersonique. Ce résultat s'applique aussi bien aux galaxies qui fusionnent, dans lesquelles l'augmentation de la turbulence compressive déclenchée par les marées compressives les conduit au régime de sursaut de formation d'étoiles. Un modèle analytique idéalisé de formation stellaire liant la densité surfacique de gaz au taux de formation stellaire comme une fonction de la présence de turbulence supersonique et la structure associée du milieu interstellaire est ensuite présenté dans le chapitre 4. Ce modèle prédit une cassure à basse densité de surface qui est suivie par un régime de loi de puissance à haute densité dans différents systèmes en accord avec les relations de formation stellaire des galaxies observées et simulées. La dernière partie de cette thèse est dédiée à la technique alternative de zoom-in cosmologique (Martig et al. 2009) et son implémentation dans le code à raffinement de maillage adaptatif RAMSES. Dans le chapitre 5, je présente les caractéristiques de base de cette technique aussi bien que certains de nos tout premiers résultats dans le contexte de l'accrétion cosmologique diffuse. / This thesis aims at making the link between galaxy evolution, morphology and internal physical processes, namely star formation as the outcome of the turbulent multiphase interstellar medium, using the cosmological zoom-in simulations, simulations of isolated and merging galaxies, and the analytic model of star formation. In Chapter 1, I explain the motivation for this thesis and briefly review the necessary background related to galaxy formation and modeling with the use of numerical simulations. I first explore the evolution of the morphology of Milky-Way-mass galaxies in a suite of zoom-in cosmological simulations through the analysis of bars. I analyze the evolution of the fraction of bars with redshift, its dependence on the stellar mass and accretion history of individual galaxies. I show in particular, that the fraction of bars declines with increasing redshift, in agreement with the observations. This work also shows that the obtained results suggest that the bar formation epoch corresponds to the transition between an early "violent" phase of spiral galaxies formation at z > 1, during which they are often disturbed by major mergers or multiple minor mergers as well as violent disk instabilities, and a late "secular" phase at z < 1, when the final morphology is generally stabilized to a disk-dominated structure. This analysis is presented in Chapter 2. Because such cosmological simulations form too many stars too early compared to observed galaxy populations, I shift the focus in Chapter 3 to star formation in a sample of low-redshift galaxy simulations in isolation at parsec and sub-parsec resolution. I study the physical origin of their star formation relations and breaks and show that the surface density threshold for efficient star formation can be related to the typical density for the onset of supersonic turbulence. This result holds in merging galaxies as well, where increased compressive turbulence triggered by compressive tides during the interaction drives the merger to the regime of starbursts. An idealized analytic model for star formation relating the surface density of gas and star formation rate as a function of the presence of supersonic turbulence and the associated structure of the ISM is then presented in Chapter 4. This model predicts a break at low surface densities that is followed by a power-law regime at high densities in different systems in agreement with star formation relations of observed and simulated galaxies. The last part of this thesis is dedicated to the alternative cosmological zoom-in technique Martig et al. 2009 and its implementation in the Adaptive Mesh Refinement code RAMSES. In Chapter 5, I will present the basic features of this technique as well as some of our very first results in the context of smooth cosmological accretion.

Avec l’avènement de la prochaine génération de télescope de plus de 30m de diamètre, il devient primordial de réduire le coût des observations et d’améliorer leur rendement scientifique. De plus il est essentiel de construire ces instruments sur des sites disposant d’une qualité optique maximale. J’ai donc essayé, au cours de ma thèse, de développer un outil fiable, facile d’utilisation et économique permettant de satisfaire ces exigences. J’ai donc utilisé le modèle de prévision météorologique Weather Research and Forecasting et le modèle de calcul de la turbulence optique Trinquet-Vernin pour prédire, plusieurs heures à l’avance, les conditions optiques du ciel tout au long de la nuit. Cette information permettrait d’améliorer la gestion du programme d’observation, appelée "flexible scheduling", et ainsi de réduire les pertes dues à la variation des conditions atmosphériques. Les résultats obtenus et les améliorations apportées au modèle WRF-TV lui permettent de présenter un bon accord entre les mesures et les prévisions ce qui est prometteur pour une utilisation réelle. Au delà de cette gestion, nous avons voulu créer un moyen d’améliorer la recherche et le test de sites astronomiquement intéressants. Nous avons donc définit un paramètre de qualité qui prend en compte les conditions météorologiques et optiques. Ce paramètre a été testé au-dessus de l’île de La Palma aux Canaries et a montré que l’Observatorio del Roque de los Muchachos est situé au meilleur emplacement de l’île. Enfin nous avons créé une routine d’automatisation du modèle WRF-TV afin d’avoir un outil opérationnel fonctionnant de manière autonome. / With the next generation of extremely large telescope having mirror with a diameter larger than 30m, it becomes essential to reduce the cost of observations and to improve their scientific efficiency. Moreover it is fundamental to build these huge infrastructures in location having the best possible optical quality. The purpose of my thesis is to bring a solution easier and more economical than before. I used the Weather Research and Forecasting (WRF) model and the Trinquet-Vernin parametrization, which computes the values of the optical turbulence, to forecast a couple of hours in advance the evolution of the sky optical quality along the coming night. This information would improve the management of observation program, called "flexible scheduling", and thereby reduce losses due to the atmospheric variations. Our results and improvements allow the model us WRF-TV to have a good agreement between previsions and in-situ measurements in different sites, which is promising for a real use in an observatory. Beyond the flexible scheduling, we wanted to create a tool to improve the search for new sites or site testing for already existing sites. Therefore we defined a quality parameter which takes into account meteorological conditions (wind, humidity, precipitable water vapor) and optical conditions (seeing, coherence time, isoplanatic angle). This parameter has been tested above La Palma in Canary island showing that the Observatorio del Roque de los Muchachos is located close to the best possible location of the island. Finally we created an automated program to use WRF-TV model in order to have an operational tool working routinely.

Résumé:<p><p>Cette dissertation s’intéresse aux méthodes numériques dans le domaine de l’acoustique. Les propriétés acoustiques d’un produit sont devenues une part intégrante de la conception. En effet, de nos jours le bruit est perçu comme une nuisance par le consommateur et constitue un critère de vente. Il y a de plus des normes à respecter. Les méthodes numériques permettent de prédire la propagation sonore et constitue dès lors un outil de conception incontournable pour réduire le temps et les coûts de développement d’un produit.<p><p>Cette dissertation considère la propagation d’ondes acoustiques dans le domaine fréquentiel en tenant compte de la présence d’un écoulement. Nous pouvons citer comme application industrielle, le rayonnement d’une nacelle de réacteur d’avion. Le but de la thèse est de proposer une nouvelle méthode et démontrer ses performances par rapport aux méthodes actuellement utilisées (i.e. la méthode des éléments finis).<p><p>L’originalité du travail consiste à étendre la méthode de partition de l’unité polynomiale dans le cadre de la propagation acoustique convectée, pour des domaines extérieurs. La simulation acoustique dans des domaines de dimensions infinies est réalisée dans ce travail à l’aide d’un couplage entre éléments finis et éléments infinis.<p><p>La dissertation présente la formulation de la méthode pour des applications axisymétriques et tridimensionnelles et vérifie la méthode en comparant les résultats numériques obtenus avec des solutions analytiques pour des applications académiques (i.e. propagation dans un conduit, rayonnement d’un multipole, bruit émis par la vibration d’un piston rigide, etc.). Les performances de la méthode sont ensuite analysées. Des courbes de convergences illustrent à une fréquence donnée, la précision de la méthode en fonction du nombre d’inconnues. Tandis que des courbes de performances présentent le temps de calcul nécessaire pour obtenir une solution d’une précision donnée en fonction de la fréquence d’excitation. Ces études de performances montrent l’intérêt de la méthode présentée.<p><p>Le rayonnement d’un réacteur d’avion a été abordé dans le but de vérifier la méthode sur une application de type industriel. Les résultats illustrent la propagation pour une nacelle axisymétrique en tenant compte de l’écoulement et la présence de matériau absorbant dans la nacelle et compare les résultats obtenus avec la méthode proposée et ceux obtenus avec la méthode des éléments finis.<p><p>Les performances de la méthode de la partition de l’unité dans le cadre de la propagation convectée en domaines infinis sont présentées pour des applications académiques et de type industriel. Le travail effectué illustre l’intérêt d’utiliser des fonctions polynomiales d’ordre élevé ainsi que les avantages à enrichir l’approximation localement afin d’améliorer la solution sans devoir créer un maillage plus fin.<p><p><p>Summary:<p><p>Environmental considerations are important in the design of many engineering systems and components. In particular, the environmental impact of noise is important over a very broad range of engineering applications and is increasingly perceived and regulated as an issue of occupational safety or health, or more simply as a public nuisance. The acoustic quality is then considered as a criterion in the product design process. Numerical prediction techniques allow to simulate vibro-acoustic responses. The use of such techniques reduces the development time and cost.<p><p>This dissertation focuses on acoustic convected radiation in outer domains such as it is the case for turbofan radiation. In the current thesis the mapped infinite partition of unity method is implemented within a coupled finite and infinite element model. This method allows to enrich the approximation with polynomial functions. <p><p>We present axisymmetric and three-dimensional formulations, verify and analyse the performance of the method. The verification compares computed results with the proposed method and analytical solutions for academic applications (i.e. duct propagation, multipole radiation, noise radiated by a vibrating rigid piston, etc.) .Performance analyses are performed with convergence curves plotting, for a given frequency, the accuracy of the computed solution with respect to the number of degrees of freedom or with performance curves, plotting the CPU time required to solve the application within a given accuracy, with respect to the excitation frequency. These performance analyses illustrate the interest of the mapped infinite partition of unity method.<p><p>We compute the radiation of an axisymmetric turbofan (convected radiation and acoustic treatments). The aim is to verify the method on an industrial application. We illustrate the radiation and compare the mapped infinite partition of unity results with finite element computations.<p><p>The dissertation presents the mapped partition of unity method as a computationally efficient method and illustrates its performances for academic as well as industrial applications. We suggest to use the method with high order polynomials and take the advantage of the method which allows to locally enrich the approximation. This last point improves the accuracy of the solution and prevent from creating a finer mesh.<p> / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished

Le développement des moyens de communications par l'intermédiaire des ondes<br /> électromagnétiques connaît une croissance sans précédent depuis quelques années, grâce<br /> notamment au développement de la téléphonie mobile. La chambre réverbérante est un<br /> moyen d'essais qui permet d'étudier l'influence de ces ondes électromagnétiques sur un<br /> appareil électronique particulier. Cependant, le fonctionnement d'une chambre<br /> réverbérante étant complexe, il est primordial de procéder à des simulations afin de<br /> déterminer quels sont les paramètres cruciaux entrant en jeu.<br /> <br /> Le travail de cette thèse consiste à modéliser et à simuler le fonctionnement d'une<br /> chambre réverbérante à l'aide de la méthode des différences finies dans le domaine<br /> temporel. Après une brève étude portant sur quelques résultats de mesures de champ et<br /> de puissances effectuées dans une chambre réverbérante, le chapitre~2 aborde les<br /> différents problèmes liés à la modélisation de la chambre. La notion de pertes étant<br /> déterminante pour évaluer le fonctionnement d'une chambre réverbérante, deux méthodes<br /> implémentant ces pertes sont aussi exposées dans ce chapitre. L'étude menée dans le<br /> chapitre~3 consiste à analyser l'influence du brasseur sur les premiers modes propres<br /> de la chambre, ceux-ci pouvant être décalés de plusieurs MHz. Le chapitre~4 présente<br /> des résultats de simulations en hautes fréquences comparés à des résultats<br /> statistiques théoriques. Le cas de la présence d'un objet au sein de la chambre<br /> pouvant perturber le champ est aussi abordé. Enfin, le chapitre~5 montre une<br /> comparaison des résultats statistiques dans le cas où l'on considère plusieurs formes<br /> de brasseurs.