Cette thèse est dédiée à la résolution numérique tridimensionnelle des équations de Maxwell harmoniques, par des méthodes de décomposition de domaine couplant des résolutions par équations intégrales entre elles. Pour traiter les problèmes de diffraction d'ondes, la méthode des équations intégrales est un outil précieux. Elle consiste à paramétrer le champ électromagnétique solution par une source définie sur la surface de l'objet diffractant, solution d'une nouvelle équation linéaire (l'équation intégrale). Pour des applications à haute fréquence, le grand nombre d'inconnues (de l'ordre du million) nous oblige à utiliser un solveur itératif pour résoudre l'équation intégrale. Le problème du conditionnement des systèmes linéaires est alors crucial. De récents développements ont permis de construire une équation intégrale performante (la GCSIE) et de conditionnement stable avec la montée en fréquence. Cependant, la présence d'une cavité large et résonnante dans l'objet diffractant (telle que la cavité moteur d'un avion) dégrade le conditionnement de cette équation. Nous proposons deux méthodes de décomposition de domaine (DDM) afin de découpler le problème de la cavité du problème extérieur. La première (DDM en Y) s'exprime en fonction des opérateurs Dirichlet-to-Neumann Y, qui sont synthétisés via la résolution de problèmes métalliques par équations intégrales dans chaque sous-domaine. La seconde (DDM en S) s'exprime en fonction des opérateurs de scattering S, synthétisés par résolution de problèmes de type métal-impédant, donc bien posés à toute fréquence. La DDM en S permet ainsi de se débarrasser des phénomènes de résonance dans les cavités. Nous proposons dans un premier temps un préconditionneur analytique pour la DDM en Y, basé sur l'opérateur électromagnétique de simple couche. Nous calculons ensuite les modes guidés le long d'un cylindre infini tangent à la cavité près de l'interface, et nous diagonalisons les opérateurs Dirichlet-to-Neumann et scattering dans la base des traces de modes guidés sur l'interface. On extrait de cette étude deux préconditionneurs spectraux respectivement pour la DDM en Y et la DDM en S. Les résultats numériques confirment l'efficacité des préconditionneurs proposés / This work is dedicated to the numerical solution of the tridimensional harmonic Maxwell equations, using domain decomposition methods coupling integral equations between them. To deal with scattering problems, integral equations methods are a precious tool. They allow to look for the electromagnetic field by parameterizing it with a source only defined on the boundary of the scattering object, solution of a new linear equation (the integral equation). For applications at high frequency, the great number of unknowns forces the use of iterative methods. To accelerate the solution of integral equations, one moreover has to ensure the good condition number of the linear systems, or to propose well-suited preconditioners. An efficient method, the GCSIE, was developed in Onera. It is an intrinsically well-conditioned integral equation whose condition number remains stable whith the frequency increase. However, the existence of large and resonant cavities (such as air intakes) deteriorates the condition number. In order to circumvent this problem, we propose two domain decomposition methods (DDM) allowing to decouple the exterior problem from the problem of the cavities. The first one (Y-DDM) is based on Dirichlet-to-Neumann operators Y, which are built through the solution of metallic problems using integral equations in each subdomain. The second one (S-DDM) is based on scattering operators S, built through the solution of problems of metallic-impedant type, which are well-posed at any frequency. The S-DDM allows to avoid the resonance phenomena inside the cavities. First, we propose an analytic preconditioner for the Y-DDM, based on the electromagnetic single layer operator. We then calculate the modes guided along an artificial infinite cylinder, that is tangent to the cavity near the interface. We diagonalize the Dirichlet-to-Neumann and scattering operators in the basis of the traces of the guided modes on the interface. We deduce from this study two spectral preconditioners for the Y-DDM and the S-DDM. The numerical results confirm the efficiency of the employed preconditioners.

Les chocs thermiques provoquent, en général, l’endommagement et la fissuration des matériaux. Ces phénomènes sont observés, par exemple, dans le revêtement de barrière thermique pour les moteurs des turbines, le traitement des surfaces ou la soudure par laser etc. Plusieurs travaux de recherche ont été réalisés au cours des dernières décennies dans l’objectif d’améliorer les performances thermiques et/ou mécaniques des matériaux sous chargement thermique. L’étude des dommages et de la fissuration des matériaux provoqués par les chocs thermiques, tels que le décollement des interfaces et de décohésion de revêtements, a reçu également une attention considérable par les chercheurs. La majorité de ces travaux utilisent les théories classiques, tels que la loi de Fourier de conduction thermique et l'hypothèse de quasi-statique. Malheureusement ces théories ne sont pas adaptées dans le cas de charges extrêmes provoqués par le choc thermique et dans le cas des matériaux micro-fissurés. En conséquence, les théories conventionnelles doivent être enrichies.L'objectif de la thèse est de montrer le rôle crucial des termes non Fourier et les termes inertiels dans le cas de choc thermique sous conditions sévères et dans le cas où les fissures sont petites. Pour cela nous avons mené des études sur deux structures particulières soumises à des chocs thermiques. Chaque structure contient une fissure parallèle au bord libre de la structure située au voisinage de ce dernier. L’influence de la présence de fissure sur la conductivité thermique est prise en compte. Nous avons utilisé la théorie Hyperbolique de transfert de chaleur par conduction pour les champs thermique et mécanique à la place de la théorie traditionnelle classique de Fourier. Pour mener cette étude, nous avons utilisé les Transformées de Laplace et de Fourier aux équations de mouvement et à l’équation de transfert de chaleur. En s’intéressant en particulier aux champs de contrainte au voisinage de la pointe de fissure et aux facteurs d'intensité de contrainte dynamiques. Le problème se ramène à la résolution d’un système d'équations intégrales singulières dans l'espace de Laplace-Fourier. On utilise une méthode d'intégration numérique pour obtenir les différents champs. Nous résolvons ensuite un système d'équations algébriques linéaires. En effectuant des inversions numériques des transformées, nous obtenons les champs de contrainte de température et les facteurs d'intensité de contrainte dynamiques dans le domaine temporel.Les résultats numériques montrent que la conductivité thermique du milieu est affectée par l’ouverture de la fissure ce qui perturberait fortement le champ de température ainsi que l'amplitude des facteurs d'intensité de contrainte dynamiques. Les amplitudes sont supérieures à celles obtenues à partir de la théorie classique de Fourier ainsi que dans le cadre de l'hypothèse quasi-statique. On constate également qu’elles oscillent au cours du temps. La prise en compte simultanément de l’influence de la fissure sur la conductivité thermique, de l'effet non-Fourier ainsi que les effetsIVd'inertie induit un couplage entre les trois phénomènes qui rendrait le problème de choc thermique très complexe. L'effet de barrière thermique induit par la fissure affecte d’une manière significative les champs de température et des contraintes. Les effets d’inertie, et des termes non-Fourier joueraient également un rôle non négligeable lorsque la longueur de la fissure est petite. Comme dans de nombreux problèmes d'ingénierie, l'initiation et la propagation des micro-fissures sont des mécanismes dont il faut tenir compte dans les prévisions de la rupture des structures. Ces effets non conventionnels ne sont plus négligeables et doivent être inclus dans l'analyse de la fracture des structures soumises à des chocs thermiques. / Thermal shock problems occur in many engineering materials and elements, which are used in high temperature applications such as thermal barrier coatings (TBCs), solid propellant of rocket-engine, pulsed-laser processing of materials, and so on. The thermal shock resistance performances and the thermal shock damages of materials, especially the interface debonding and spallation of coatings, have received considerable attention in both analysis and design. Some conventional theories, such as the Fourier’s law of thermal conduction and the quasi-static assumption of the thermoelastic body, may no longer be appropriate because of the extreme loads provoked by the thermal shock. Therefore, these conventional theories need to be enriched or revised.The objective of this thesis is to develop the solutions of the transient temperature field and thermal stresses around a partially insulated crack in a thermoelastic strip under thermal shock loading. The crack lies parallel to the heated traction free surface. The thermal conductivity of the crack gap is taken into account. Hyperbolic heat conduction theory is used in solving the temperature field instead of the traditional Fourier thermal conduction theory. Equations of motion are applied to obtain the stress fields and the dynamic stress intensity factors of the crack. The Laplace and Fourier transforms are applied to solve the thermal-elastic governing equations such that the mixed boundary value problems are reduced to solving a singular integral equations system in Laplace-Fourier space. The numerical integration method is applied to get the temperature field and stress fields, respectively. The problems are then solved numerically by converting the singular integral equations to a linear algebraic equations system. Finally, numerical inversions of the Laplace transform are performed to obtain the temperature field and dynamic stress intensity factors in the time domain.Numerical results show that the thermal conductivity of the crack gap strongly affects the uniformity of the temperature field and consequently, the magnitude of the dynamic stress intensity factors of the crack. The stress intensity factors would have higher amplitude and oscillating feature comparing to those obtained under the conventional Fourier thermal conduction and quasi-static hypotheses. It is also observed that the interactions of the thermal conductivity of the crack gap, the non-Fourier effect and the inertia effects would make the dynamic thermal shock problem more complex. The magnitude of the thermal barrier, non-Fourier and inertia effects is estimated for some practical cases.

O presente trabalho objetiva um estudo de aplicação do Método dos Elementos de Contorno à análise de problemas de propagação de ondas eletromagnéticas. O Método baseia-se numa formulação integral que elimina todas as operações de integração em domínio, restando apenas as de contorno. Inicialmente faz-se um estudo dos fundamentos teóricos do método, apresentando-o de forma genérica e encaminhando sua aplicação à equação de Helmholtz. Os procedimentos computacionais desenvolvidos para a implementação do método viabilizam a solução eficiente de problemas de interesse, envolvendo diferentes meios com ou sem perdas. São apresentados resultados de simulações realizadas que confirmam a aplicabilidade do método, permitindo também uma análise de seu desempenho através da variação de parâmetros, como, por exemplo, número de elementos na discretização e função de interpolação. / This work deals with a study of application of the Boundary Element Method (BEM) directed to electromagnetic guided wave propagation. This method relies on an integral formulation that does not need any domain integration. Only boundary integrations have to be performed. The work begins with a study of the theoretical foundations of the method, presenting its general formulation and then directing it to Helmholtz’s equation solution. Developed computational procedures allow efficient application of the method to real problems with more than one medium, with or without losses. Simulations results are presented which confirm the applicability of the method and allow the analysis of its performance through parameters variation as, for example, the number of discretized elements and interpolation function.

Le principe du contrôle par flux de fuite magnétique (Magnetic Flux Leakage MFL) consiste à aimanter la pièce à contrôler par un champ magnétique et à détecter à l'aide d'un capteur magnétique les fuites des lignes du champ qui résultent de la présence d'un défaut dans la pièce. Dans le but d'améliorer les performances d'un dispositif de détection, le CEA et la société Vallourec collaborent pour développer un modèle numérique dédié au contrôle virtuel des défauts longitudinaux dans les tubes ferromagnétiques. Le dispositif expérimental comprend un circuit magnétique tournant à une vitesse constante autour du tube qui défile. Dans le cadre de cette thèse, on débute le problème de la modélisation sans tenir compte des effets de la vitesse de rotation, il s'agit donc de résoudre un problème d'électromagnétisme en régime magnétostatique.Pour résoudre ce problème, on propose de comparer une approche semi-analytique basée sur le formalisme des équations intégrales (EI) et une approche purement numérique utilisant les éléments finis (EF).Dans la première partie de cette thèse, après avoir établi le formalisme théorique par EI, un premier modèle considérant des matériaux ferromagnétiques à perméabilité magnétique constante (régime linéaire) a été mis en œuvre en 2D. Ce modèle a été appliqué pour un exemple de système extrait de la littérature et validé numériquement par une comparaison des résultats EI/EF. Pour une meilleure détection, il est opportun de saturer magnétiquement la pièce. Le matériau ferromagnétique est alors caractérisé par une courbe B(H) non-linéaire. Par conséquent, la deuxième partie de la thèse a été consacrée à la mise en œuvre du modèle en régime non linéaire qui tient compte de cette caractéristique.Différentes méthodes de discrétisation ont été étudiées afin de réduire le nombre d'inconnues et le temps de calcul. L'originalité de la thèse réside dans l'utilisation des fonctions d'interpolation d'ordre élevé (polynôme de Legendre) pour une discrétisation des équations intégrales par une approche de type Galerkin. Les premiers essais de validation numérique de ce modèle ont été effectués sur un système MFL simplifié. Des premiers essais de validation expérimentale pour des données obtenues par EF ont été effectués en deux phases : La première a consisté à vérifier le distribution du champ magnétique pour un tube sain et en régime magnétostatique. La deuxième phase a consisté à calculer la réponse d'un défaut dans le tube ferromagnétique en tenant en compte les effets éventuels de la rotation du circuit magnétique par rapport au tube. / The principle of the non destructive testing by magnetic flux leakage (MFL) is to magnetize the part to be inspected by a magnetic field and to detect a flaw thanks to magnetic leakage field lines due to the strong decreasing of the magnetic permeability in the flawed region. In order to improve the performance of detection, the CEA and the Vallourec society collaborate to develop a numerical model dedicated to the virtual NDT of longitudinal defects in ferromagnetic tubes. The experimental system includes a magnetic circuit rotating at a constant speed around the tube to be inspected. The modeling task is started without considering the effects of the rotational speed, so the magnetostatic regime is considered to solve the modeling problem. In the framework of this thesis, we propose to compare a semi-analytical approach based on the formalism of integral equations method (IEM) and a purely numerical approach using finite element method (FEM).In the first part of this thesis, the theoretical formalism was established. A first simple discretization scheme is been implemented in the linear regime considering a constant magnetic permeability. This first numerical model has been validated for a simplified MFL configuration extracted and modified from the literature.For better detection, it is wishable to magnetically saturate the piece under-test. The ferromagnetic material is then characterized by a B(H) curve. Therefore, the second part of the thesis was devoted to the implementation of the model in the non-linear regime that takes into account this non-linear characteristic. Different discretization schemes have been studied in order to reduce the number of unknowns and the computational time. The originality of the thesis lies in the use of basis function of high order (Legendre polynomials) associated to a Galerkin approach for the discretization of integral equations. The first numerical result has been validated on a simplified MFL system. The first results of the experimental validation based on simulated data obtain by FEM have been performed in two steps. The first one consists to verify the distribution of the magnetic field for a ferromagnetic tube without any defect and in the magnetostatic regime. The objective of the second one was to compute the response of the flaw and to evaluate the effects of the rotational speed of the magnetic circuit around the tube.

Cette thèse s'articule autour de deux thématiques : les méthodes numériques pour la propagation d'ondes acoustiques sous écoulement et les méthodes de réduction de modèles. Dans la première thématique, nous développons une méthode de couplage d'éléments finis et d'éléments de frontière pour résoudre l'équation d'Helmholtz convectée, lorsque l'écoulement est uniforme à l'extérieur d'un domaine borné. En particulier, nous proposons une formulation bien posée à toutes les fréquences de la source. Dans la deuxième thématique, nous proposons une solution au problème classique d'accumulation d'arrondis machine qui survient en calculant l'estimateur d'erreur a posteriori dans la méthode des bases réduites. Par ailleurs, nous proposons une méthode non intrusive pour calculer une approximation sous forme séparée des systèmes linéaires résultant de l'approximation en dimension finie de problèmes aux limites dépendant d'un ou plusieurs paramètres / This thesis has two topics : numerical methods for acoustic wave propagation in a flow and reduced order models. In the first topic, we develop a coupled finite element and boundary element method to solve the convected Helmholtz equation, when the flow is uniform outside a bounded domain. In particular, we propose a formulation that is well-posed at all the frequencies of the source. In the second topic, we propose a solution to the classical problem of round-off error accumulation that occurs when computing the a posteriori error bound in the reduced basis method. Furthermore, we propose a non intrusive method for the approximation, in a separated representation form, of linear systems resulting from the finite-dimensional approximation of boundary-value problems depending on one or several parameters

O presente trabalho objetiva um estudo de aplicação do Método dos Elementos de Contorno à análise de problemas de propagação de ondas eletromagnéticas. O Método baseia-se numa formulação integral que elimina todas as operações de integração em domínio, restando apenas as de contorno. Inicialmente faz-se um estudo dos fundamentos teóricos do método, apresentando-o de forma genérica e encaminhando sua aplicação à equação de Helmholtz. Os procedimentos computacionais desenvolvidos para a implementação do método viabilizam a solução eficiente de problemas de interesse, envolvendo diferentes meios com ou sem perdas. São apresentados resultados de simulações realizadas que confirmam a aplicabilidade do método, permitindo também uma análise de seu desempenho através da variação de parâmetros, como, por exemplo, número de elementos na discretização e função de interpolação. / This work deals with a study of application of the Boundary Element Method (BEM) directed to electromagnetic guided wave propagation. This method relies on an integral formulation that does not need any domain integration. Only boundary integrations have to be performed. The work begins with a study of the theoretical foundations of the method, presenting its general formulation and then directing it to Helmholtz’s equation solution. Developed computational procedures allow efficient application of the method to real problems with more than one medium, with or without losses. Simulations results are presented which confirm the applicability of the method and allow the analysis of its performance through parameters variation as, for example, the number of discretized elements and interpolation function.

In this paper, we study ill-posedness concepts of nonlinear and linear inverse problems in a Hilbert space setting. We define local ill-posedness of a nonlinear operator equation $F(x) = y_0$ in a solution point $x_0$ and the interplay between the nonlinear problem and its linearization using the Frechet derivative $F\acent(x_0)$ . To find an appropriate ill-posedness concept for the linarized equation we define intrinsic ill-posedness for linear operator equations $Ax = y$ and compare this approach with the ill-posedness definitions due to Hadamard and Nashed.

No description available.

No description available.

No description available.