La réduction des nuisances sonores est un enjeu permanent pour les acteurs de l’aéronautique. L’optimisation de la réduction de bruit apportée par les traitements acoustiques tapissant la nacelle des réacteurs turbofan passe par une caractérisation précise des matériaux employés dans l’environnement aéroacoustique d’utilisation, qui met en jeu un écoulement rasant de vitesse importante combiné à de forts niveaux sonores. L’objectif de cette thèse est de développer une méthode inverse pour la détermination de l’impédance acoustique de liners soumis à un écoulement rasant, basée sur des mesures non intrusives du champ de vitesse acoustique au-dessus du matériau par Vélocimétrie Laser Doppler (LDV). L’impédance de liner est obtenue par minimisation de l’écart entre le champ de vitesse acoustique mesuré et le champ simulé numériquement en résolvant les équations d’Euler linéarisées bidimensionnelles harmoniques, discrétisées par un schéma Galerkin discontinu. Le gradient de la fonction objectif minimisée est calculé via la résolution, à chaque itération, des équations directes et adjointes. Une première étape de validation du solveur est effectuée sur des cas-tests académiques, puis sur des cas expérimentaux impliquant des mesures de pression acoustique en paroi rigide opposée au liner. Dans un second temps, la méthode est appliquée à des mesures de vitesse acoustique obtenues par LDV dans le banc B2A de l’ONERA en l’absence d’écoulement. La dernière étape consiste à prendre en compte l’effet d’un écoulement rasant de profil cisaillé. Les impédances identifiées à partir de mesures LDV en présence d’écoulement ont notamment permis de gagner en compréhension sur les phénomènes d’absorption intervenant dans le banc B2A. / While aircraft noise constraints become increasingly stringent, efficient duct treatment of turbofan engines requires an accurate knowledge of liner impedance with grazing flow at high acoustic levels. This thesis aims at developing an impedance eduction method in the presence of grazing flow. The inverse process is based on acoustic velocity fields acquired by Laser Doppler Velocimetry (LDV) above the liner. The liner acoustic impedance is obtained by minimization of the distance between the measured acoustic velocity field and the simulated one. Computations rely on the resolution of the 2D linearized Euler equations in the harmonic domain, spatially discretized by a discontinuous Galerkin scheme. The gradient of the objective function is achieved by the resolution, at each iteration on the liner impedance, of the direct and adjoint equations. The solver is first validated on academic test cases, then on experimental results of acoustic pressure measurements at the rigid wall opposite the liner. Secondly the method is applied to acoustic velocity measurements obtained by LDV above the liner without flow, in the ONERA B2A test bench. The last step consists in taking into account the effects of a sheared grazing flow. The impedances educed from LDV measurements in the presence of flow namely allowed to gain insight into the absorption phenomena occuring in the B2A test bench.

On s'est intéressé dans ce travail à la résolution numérique des équations d'Euler linéarisées autour d'un écoulement stationnaire, subsonique et assez régulier. Dans le but d'obtenir des matrices symétriques dans ces équations, et in fine, une équation d'équilibre énergétique, nous considérons la linéarisation d'une forme symétrique des équations d'Euler tridimensionnelles. Nous proposons un schéma non-diffusif de type Galerkin discontinu en domaine temporel (GDDT) s'appuyant sur une formulation centrée-élément avec des ux numériques totalement centrés et un schéma en temps explicite de type saute-mouton, ce qui permet d'obtenir une approximation sans dissipation et fournit une estimation précise des variations de l'énergie aéroacoustique. En effet, dans le cas général de la linéarisation autour d'un écoulement non-uniforme, il existe une équation d'équilibre énergétique au niveau continu que nous vérifions au niveau discret. Nous montrons qu'il existe un terme source discret permettant de conserver exactement l'énergie ce qui permet de prouver la stabilité de notre schéma. Ainsi, notre schéma non-diffusif de type GDDT fournit un outil précis pour contrôler des phénomènes tel que les instabilités de Kelvin-Helmholtz. Nous illustrons la capacité de notre méthode aussi bien sur plusieurs cas tests académiques que sur différentes configurations complexes grâce à une implémentation parallèle.

Cette thèse porte sur la résolution théorique et numérique des équations de Maxwell dans le domaine temporel ou fréquentiel. Dans une première partie, on démontre l'existence et l'unicité mathématique de la solution du problème d'évolution. On s'intéresse également au comportement asymptotique en temps de cette solution lorsque le second membre des équations est sinusoïdal en temps. L'approche utilisée fait appel à la théorie des systèmes hyperboliques linéaires du premier ordre, au théorème de Hille-Yosida, aux principes d'amplitude-limite et d'absorption-limite, ainsi qu'à des théorèmes de traces (dans le cas du problème aux limites). Dans un second temps, on développe une approximation par éléments finis discontinus du problème fréquentiel, basée sur une décomposition de la matrice des flux en partie positive et négative (méthode de flux-splitting). Cette approche autorise l'utilisation de maillages totalement déstructurés. Une étude d'erreur lorsque le pas h du maillage tend vers zéro est proposée. Un algorithme itératif de résolution du problème discret, basé sur une décomposition de domaine sans recouvrement, est ensuite décrit. On démontre sa convergence vers l'unique solution discrète. L'implémentation sur un ordinateur à architecture massivement parallèle (IPSC 860) a été réalisée. Enfin, on construit une équation intégrale adaptée à la méthode, pour la résolution des problèmes en domaine non borné. Des expériences numériques sont décrites dans le cas d'éléments finis de type P0 (approximation constante par élément).

On s'intéresse dans ce travail à la simulation et l'approximation d'écoulements diphasiques avec transition de phase. Il s'agit de modéliser la dynamique d'une bulle de cavitation. Le modèle repose sur les équations d'Euler en coordonnées sphériques, l'interface gaz-liquide étant indiquée par une fonction de couleur. Dans une première partie, aucun transfert de masse n'est supposé. Un schéma Lagrange-projection est d'abord proposé : seule l'interface est déplacée à la vitesse du fluide. La projection repose sur un algorithme de suppression-découpage qui assure que les volumes ne deviennent pas négatifs. Le second chapitre traite du terme source géométrique. On construit un schéma équilibre articulé autour du schéma VFRoe-ncv pour lequel on propose une correction entropique non paramétrique. Un méthode d'ordre élevé de type Galerkin discontinu est ensuite étudiée dans le cadre de la magnétohydrodynamique. L'intégration en temps est réalisée par une méthode de type Adams-Bashforth, qui s'avère bien adaptée aux algorithmes de pas de temps local. La deuxième partie de la thèse traite de la modélisation du changement de phase liquide-vapeur. L'inf-convolution et la transformée de Legendre définissent une structure naturelle pour la construction de loi de pression de mélange. En particulier on montre que la construction de Maxwell est équivalente à la construction de l'enveloppe convexe de l'énergie de van der Waals. Un algorithme de transformée de Legendre rapide, adapté à la prise en compte correcte des bords du domaine de calcul de loi d'état, est développé. La méthode est appliquée à la construction de lois de pression tabulées de mélanges binaires miscibles ou immiscibles de gaz raides.

Pour la conception des cœurs de réacteurs de 4ème génération, une précision accrue est requise pour les calculs des différents paramètres neutroniques. Les ressources mémoire et le temps de calcul étant limités, une solution consiste à utiliser des méthodes de raffinement de maillage afin de résoudre l'équation de transport des neutrons. Le flux neutronique, solution de cette équation, dépend de l'énergie, l'angle et l'espace. Les différentes variables sont discrétisées de manière successive. L'énergie avec une approche multigroupe, considérant les différentes grandeurs constantes sur chaque groupe, l'angle par une méthode de collocation, dite approximation Sn. Après discrétisation énergétique et angulaire, un système d'équations hyperboliques couplées ne dépendant plus que de la variable d'espace doit être résolu. Des éléments finis discontinus sont alors utilisés afin de permettre la mise en place de méthodes de raffinement dite hp. La précision de la solution peut alors être améliorée via un raffinement en espace (h-raffinement), consistant à subdiviser une cellule en sous-cellules, ou en ordre (p-raffinement) en augmentant l'ordre de la base de polynômes utilisée.Dans cette thèse, les propriétés de ces méthodes sont analysées et montrent l'importance de la régularité de la solution dans le choix du type de raffinement. Ainsi deux estimateurs d'erreurs permettant de mener le raffinement ont été utilisés. Le premier, suppose des hypothèses de régularité très fortes (solution analytique) alors que le second utilise seulement le fait que la solution est à variations bornées. La comparaison de ces deux estimateurs est faite sur des benchmarks dont on connaît la solution exacte grâce à des méthodes de solutions manufacturées. On peut ainsi analyser le comportement des estimateurs au regard de la régularité de la solution. Grâce à cette étude, une stratégie de raffinement hp utilisant ces deux estimateurs est proposée et comparée à d'autres méthodes rencontrées dans la littérature. L'ensemble des comparaisons est réalisé tant sur des cas simplifiés où l'on connaît la solution exacte que sur des cas réalistes issus de la physique des réacteurs.Ces méthodes adaptatives permettent de réduire considérablement l'empreinte mémoire et le temps de calcul. Afin d'essayer d'améliorer encore ces deux aspects, on propose d'utiliser des maillages différents par groupe d'énergie. En effet, l'allure spatiale du flux étant très dépendante du domaine énergétique, il n'y a a priori aucune raison d'utiliser la même décomposition spatiale. Une telle approche nous oblige à modifier les estimateurs initiaux afin de prendre en compte le couplage entre les différentes énergies. L'étude de ce couplage est réalisé de manière théorique et des solutions numériques sont proposées puis testées.

La réduction des nuisances sonores est un enjeu permanent pour les acteurs de l'aéronautique. L'optimisation de la réduction de bruit apportée par les traitements acoustiques tapissant la nacelle des réacteurs turbofan passe par une caractérisation précise des matériaux employés dans l'environnement aéroacoustique d'utilisation, qui met en jeu un écoulement rasant de vitesse importante combiné à de forts niveaux sonores. L'objectif de cette thèse est de développer une méthode inverse pour la détermination de l'impédance acoustique de liners soumis à un écoulement rasant, basée sur des mesures non intrusives du champ de vitesse acoustique au-dessus du matériau par Vélocimétrie Laser Doppler (LDV). L'impédance de liner est obtenue par minimisation de l'écart entre le champ de vitesse acoustique mesuré et le champ simulé numériquement en résolvant les équations d'Euler linéarisées bidimensionnelles harmoniques, discrétisées par un schéma Galerkin discontinu. Le gradient de la fonction objectif minimisée est calculé via la résolution, à chaque itération, des équations directes et adjointes. Une première étape de validation du solveur est effectuée sur des cas-tests académiques, puis sur des cas expérimentaux impliquant des mesures de pression acoustique en paroi rigide opposée au liner. Dans un second temps, la méthode est appliquée à des mesures de vitesse acoustique obtenues par LDV dans le banc B2A de l'ONERA en l'absence d'écoulement. La dernière étape consiste à prendre en compte l'effet d'un écoulement rasant de profil cisaillé. Les impédances identifiées à partir de mesures LDV en présence d'écoulement ont notamment permis de gagner en compréhension sur les phénomènes d'absorption intervenant dans le banc B2A.

Ce travail porte sur l'étude d'une méthode d'éléments finis discontinus (ou méthode de type Galerkin Discontinu, DGTD) basée sur l'utilisation d'un maillage héxaédrique régulier, proposée pour la résolution des équations de Maxwell dans le domaine temporel, afin de l'adapter à la caractérisation de structures rayonnantes et de l'associer à des techniques de domaines fictifs.<br />On présente tout d'abord une méthode Galerkin Discontinu s'appuyant sur une formulation centrée pour approcher les flux numériques aux interfaces du maillage et sur un schéma en temps explicite de type saute-mouton. Ainsi, le schéma obtenu est non-diffusif, stable, peu dispersif, parfaitement adapté à l'utilisation de maillages localement raffinés de manière non-conforme. La méthode a été dotée de parois absorbantes performantes (modèle Unsplit-PML), permettant de prendre en compte facilement des objets à cheval entre le domaine de calcul et la couche absorbante. Nous avons ensuite utilisé la méthode pour effectuer des calculs d'impédances, de paramètres S et de T.O.S. sur des structures rayonnantes planaires. La comparaison entre la simulation et la mesure de ces structures montre le bon fonctionnement de la méthode.<br />Nous avons alors couplé une méthode de domaines fictifs avec la méthode DGTD afin de prendre en compte la présence d'objets métalliques à géométries complexes. La méthode des domaines fictifs utilise deux maillages de manière indépendante: un maillage cartésien, pour faire évoluer le champ électromagnétique dans l'espace libre, et un maillage surfacique qui permet de prendre en compte l'objet métallique. La convergence de la méthode (pour la méthode FDTD) est liée à une relation de compatibilité entre le maillage volumique et le maillage surfacique: le plus petit élément du maillage surfacique impose la taille des éléments du maillage volumique. Ainsi, pour des objets présentant de tout petits détails, cette condition n'est assurée que si le maillage volumique est de l'ordre du plus petit élément du maillage surfacique, ce qui peut devenir extrêmement contraignant sans le recours à des techniques de raffinement local, telle que celle rendue possible par la méthode Galerkin Discontinu et utilisée ici.

Ce travail est dédié à l´étude et la simulation numérique de la rupture dynamique des séismes en deux et trois dimensions d´espace par une méthode d´éléments finis discontinus. Après avoir transformé le système de l´élastodynamique en un système hyperbolique symétrique du premier ordre, nous proposons un schéma numérique basé sur des flux centrés et un schéma explicite en temps de type saute-mouton. A travers l´étude d´une énergie discrète du système, nous spécifions les conditions aux limites sur la faille afin de prendre en compte de manière faible la rupture en mode cisaillant que nous traitons. Nous montrons, qu´en l´absence de tractions tangentielles sur la faille, cette énergie est parfaitement conservée. Nous illustrons la capacité de notre méthode à travers divers cas tests sur des configurations complexes grâce à une implémentation parallèle.

L'étude de la susceptibilité électromagnétique de circuits électroniques nécessite l'utilisation d'un outil de simulation rapide, précis et suffisamment flexible pour intégrer les dernières innovations technologiques. La méthode itérative basée sur le concept d'onde (notée WCIP pour Wave Concept Iterative Procedure) initialement proposée par H. Baudrand est particulièrement adaptée pour la modélisation numérique de circuits multicouches à plusieurs niveaux de métallisation. Pour ce type de circuits, elle se révèle être l'une des méthodes qui utilise le plus petit nombre d'inconnues pour atteindre une précision donnée. Néanmoins, la WCIP n'est pas adaptée à la prise en compte des diélectriques inhomogènes et des trous d'interconnexion. L'objectif de la thèse est de s'affranchir de ces limitations par un couplage avec des méthodes numériques volumiques. En premier lieu, l'hybridation a été mise en œuvre avec une méthode basée sur la théorie des lignes de transmission pour des raisons de correspondance de maillages. Par la suite, le couplage avec une technique d'éléments finis de type Galerkin Discontinu (notée GD) Hybridée permet d'atteindre des objectifs de précision et de rapidité car GD apporte une flexibilité dans la discrétisation. En effet, c'est une méthode d'éléments finis non conforme qui permet notamment de faire varier d'un élément à l'autre l'ordre polynomial d'approximation. On a ainsi développé une nouvelle méthode numérique hybride couplant la WCIP avec des méthodes volumiques qui offrent plus de souplesse pour la prise en compte des milieux complexes. Enfin, une stratégie de résolution par décomposition de domaines est également abordée à la fin du manuscrit. / Electromagnetic susceptibility study of electronic circuits requires the use of a simulation tool which is fast, accurate and flexible enough to incorporate the latest technological innovations. The Wave Concept Iterative Procedure (WCIP) initially proposed by H. Baudrand is particularly adapted for numerical modeling of multilayered circuits with multilevel metallization. For this kind of circuits, it turns out to be one of the methods that uses the smallest number of unknowns to reach a given accuracy. However, the WCIP is not appropriate for inhomogeneous dielectric substrates and metallized via holes. The aim of this PhD thesis is to overcome these limitations coupling the WCIP with volume numerical methods. First, hybridization is carried out with the Frequency Domain Transmission Line Matrix (denoted FDTLM) assuming matching meshes at the interface between computational domains of both methods. Subsequently, the coupling with a finite element technique like a Hybridized Discontinuous Galerkin (denoted DG) method is considered to achieve the objectives of accuracy and speed because DG brings flexibility in the discretization. Indeed, it is a nonconforming finite element method which allows in particular changing the polynomial approximation order from one element to another. Therefore, a new hybrid method is developed coupling the WCIP with volume numerical methods which offer more flexibility for dealing with complex environments. Finally, a domain decomposition solution strategy is also discussed at the end of the manuscript.

Pour la conception des cœurs de réacteurs de 4ème génération, une précision accrue est requise pour les calculs des différents paramètres neutroniques. Les ressources mémoire et le temps de calcul étant limités, une solution consiste à utiliser des méthodes de raffinement de maillage afin de résoudre l'équation de transport des neutrons. Le flux neutronique, solution de cette équation, dépend de l'énergie, l'angle et l'espace. Les différentes variables sont discrétisées de manière successive. L'énergie avec une approche multigroupe, considérant les différentes grandeurs constantes sur chaque groupe, l'angle par une méthode de collocation, dite approximation Sn. Après discrétisation énergétique et angulaire, un système d'équations hyperboliques couplées ne dépendant plus que de la variable d'espace doit être résolu. Des éléments finis discontinus sont alors utilisés afin de permettre la mise en place de méthodes de raffinement dite hp. La précision de la solution peut alors être améliorée via un raffinement en espace (h-raffinement), consistant à subdiviser une cellule en sous-cellules, ou en ordre (p-raffinement) en augmentant l'ordre de la base de polynômes utilisée.Dans cette thèse, les propriétés de ces méthodes sont analysées et montrent l'importance de la régularité de la solution dans le choix du type de raffinement. Ainsi deux estimateurs d'erreurs permettant de mener le raffinement ont été utilisés. Le premier, suppose des hypothèses de régularité très fortes (solution analytique) alors que le second utilise seulement le fait que la solution est à variations bornées. La comparaison de ces deux estimateurs est faite sur des benchmarks dont on connaît la solution exacte grâce à des méthodes de solutions manufacturées. On peut ainsi analyser le comportement des estimateurs au regard de la régularité de la solution. Grâce à cette étude, une stratégie de raffinement hp utilisant ces deux estimateurs est proposée et comparée à d'autres méthodes rencontrées dans la littérature. L'ensemble des comparaisons est réalisé tant sur des cas simplifiés où l'on connaît la solution exacte que sur des cas réalistes issus de la physique des réacteurs.Ces méthodes adaptatives permettent de réduire considérablement l'empreinte mémoire et le temps de calcul. Afin d'essayer d'améliorer encore ces deux aspects, on propose d'utiliser des maillages différents par groupe d'énergie. En effet, l'allure spatiale du flux étant très dépendante du domaine énergétique, il n'y a a priori aucune raison d'utiliser la même décomposition spatiale. Une telle approche nous oblige à modifier les estimateurs initiaux afin de prendre en compte le couplage entre les différentes énergies. L'étude de ce couplage est réalisé de manière théorique et des solutions numériques sont proposées puis testées. / The different neutronic parameters have to be calculated with a higher accuracy in order to design the 4th generation reactor cores. As memory storage and computation time are limited, adaptive methods are a solution to solve the neutron transport equation. The neutronic flux, solution of this equation, depends on the energy, angle and space. The different variables are successively discretized. The energy with a multigroup approach, considering the different quantities to be constant on each group, the angle by a collocation method called Sn approximation. Once the energy and angle variable are discretized, a system of spatially-dependent hyperbolic equations has to be solved. Discontinuous finite elements are used to make possible the development of $hp-$refinement methods. Thus, the accuracy of the solution can be improved by spatial refinement (h-refinement), consisting into subdividing a cell into subcells, or by order refinement (p-refinement), by increasing the order of the polynomial basis.In this thesis, the properties of this methods are analyzed showing the importance of the regularity of the solution to choose the type of refinement. Thus, two error estimators are used to lead the refinement process. Whereas the first one requires high regularity hypothesis (analytical solution), the second one supposes only the minimal hypothesis required for the solution to exist. The comparison of both estimators is done on benchmarks where the analytic solution is known by the method of manufactured solutions. Thus, the behaviour of the solution as a regard of the regularity can be studied. It leads to a hp-refinement method using the two estimators. Then, a comparison is done with other existing methods on simplified but also realistic benchmarks coming from nuclear cores.These adaptive methods considerably reduces the computational cost and memory footprint. To further improve these two points, an approach with energy-dependent meshes is proposed. Actually, as the flux behaviour is very different depending on the energy, there is no reason to use the same spatial discretization. Such an approach implies to modify the initial estimators in order to take into account the coupling between groups. This study is done from a theoretical as well as from a numerical point of view.