Un entier naturel est dit $y$-friable lorsque son plus grand facteur premier n'excède pas $y$. Ce travail est consacré à l'étude des entiers friables dans le cadre de la théorie analytique et probabiliste des nombres. La première partie est dévolue à un problème posé par Davenport en 1937, qui consiste à déterminer les conditions de validité de diverses généralisations de son développement de la fonction sinus en série de parties fractionnaires. Ces généralisations peuvent être décrites par un couple de fonctions arithmétiques, liées par la relation de convolution $f=g*\1$. Nous traitons le cas où $g$ est la fonction de Piltz d'ordre $z\in\CC$. La deuxième partie est consacrée à l'étude du comportement asymptotique de la constante optimale dans une version friable de l'inégalité de Turán-Kubilius. Précisant des résultats récents de La Bretèche et Tenenbaum, nous généralisons au cas friable une formule asymptotique de la variance d'une fonction arithmétique additive, établie par Hildebrand en 1983.

La résolution analytique du problème de la stabilité d'une suspension solide (particules rigides de forme sphérique) dans le système de Taylor-Couette cylindrique a été menée. On s'est basé sur les travaux de Ali and Lueptow (2002) pour formuler les équations régissant la stabilité de l'écoulement dans le cadre d'une théorie linéaire. Ces équations sont valables dans l'approximation du faible espace annulaire et ont pour but la prévision de l'instabilité primaire. A cet effet, nous avons utilisé une méthode variationnelle telle que la méthode de Galerkin pour résoudre le problème aux valeurs propres conduisant à établir le diagramme de stabilité liée au nombre d'onde au voisinage de l'état critique du développement de la première instabilité. Tout d'abord, on a cherché à mettre au point les calculs dans le cas de l'écoulement de Taylor-Couette classique en se référant aux travaux de Ali and Lueptow (2002). Ensuite on a procédé à la résolution systématique des équations du mouvement et l'on cherche à déterminer le critère d'apparition des instabilités en présence de particules en suspension et l'on détermine simultanément les paramètres de couplage entre forces d'interaction liquide-solide. L'ensemble des travaux ainsi réalisés permettront de lever la contradiction fondamentale entre la théorie et l'expérience. L'étude expérimentale a permis d'analyser les effets de la concentration des particules en suspension (disques) et du rapport d'aspect radial ' sur l'apparition des instabilités dans le système de Taylor-Dean. Le dispositif expérimental utilisé est constitué d'un cylindre intérieur tournant et le cylindre extérieur est maintenu fixe. Le comportement rhéologique du fluide utilisé est viscoplastique obéissant au modèle de Herschel Bulkley. L'apparition des instabilités est examinée en utilisant une technique de visualisation. Pour une concentration donnée dans l'intervalle étudié, la nature des structures apparaissant dans le système d'écoulement dépend ', alors que pour une valeur donnée de ' dans l'intervalle étudié, la valeur du nombre de Taylor critique Tac dépend de la concentration des particules. Nous avons obtenu que le nombre de Taylor critique Tac correspondant au déclenchement de la première instabilité évolue non linéairement en fonction de '. De plus, nous avons examiné expérimentalement les effets de limitation axiale (effet de bords) sur le déclenchement des instabilités dans le système de Taylor-Dean. Les résultats obtenus montrent que les bords tournants n'affectent pas le type de structures qui apparaissent dans le système d'écoulement de Taylor-Dean. Cependant, ils influencent le seuil critique d'apparition des instabilités qui est marquée par des valeurs élevées du nombre de Taylor critique pour des bords tournants ce qui indique un effet stabilisant des bords mobiles.

Cette thèse se concentre sur le développement d'une méthode précise et efficace pour la simulation des grandes échelles (LES) des écoulements turbulents. Une approche de la LES basée sur la méthode variationnelle multi-échelles (VMS) est considérée. La VMS applique aux équations de la dynamique des fluides une séparation d'échelles a priori sans recours à des hypothèses sur les conditions aux limites ou sur l'uniformité du maillage. Afin d'assurer effectivement une séparation d'échelles dans l'espace des nombres d'onde associé, nous choisissons d'utiliser les ondelettes de deuxième génération (SGW), une base polynomiale qui présente des propriétés de localisation spatiale-fréquence optimales. A partir de la séparation d'échelles ainsi réalisée, l'action du modèle sous-maille est limitée à un intervalle de nombres d'onde proche de la coupure spectrale. Cette approche VMS-LES basée sur les ondelettes est désignée par WAVVMS-LES. Elle est incorporée dans un solveur d'ordre élevé pour la simulation des écoulements incompressibles sur la base d'une méthode de Galerkin discontinue (DG-FEM) stabilisée pour la pression. La méthode est évaluée par réalisation de LES sur des maillages fortement sous-résolus pour le cas test du tourbillon de Taylor-Green 3D à deux nombres de Reynolds différents. / This thesis focuses on the development of an accurate and efficient method for performing Large-Eddy Simulation (LES) of turbulent flows. An LES approach based upon the Variational Multiscale (VMS) method is considered. VMS produces an a priori scale-separation of the governing equations, in a manner which makes no assumptions on the boundary conditions and mesh uniformity. In order to ensure that scale-separation in wavenumber is achieved, we have chosen to make use of the Second Generation Wavelets (SGW), a polynomial basis which exhibits optimal space-frequency localisation properties. Once scale-separation has been achieved, the action of the subgrid model is restricted to the wavenumber band closest to the cutoff. We call this approach wavelet-based VMS-LES (WAV-VMS-LES). This approach has been incorporated within the framework of a high-order incompressible flow solver based upon pressure-stabilised discontinuous Galerkin FEM (DG-FEM). The method has been assessed by performing highly under-resolved LES upon the 3D Taylor-Green Vortex test case at two different Reynolds numbers.

Cette thèse se penche sur la propagation d’ondes au sein du système coupléTerre-océan-atmosphère. La compréhension de ces phénomènes a une importance majeure pour l’étude de perturbations sismiques et d’explosions atmosphériques notamment dans le cadre de missions spatiales planétaires. Les formes d’ondes issues du couplage fluide-solide permettent d’obtenir de précieuses informations sur la source du signal ou les propriétés des milieux de propagation. On développe donc deux outils numériques d’ordre élevé pour la propagation d’ondes acoustiques et de gravité. L'u en différences finies et se concentre sur le milieu atmosphérique et la propagation d’ondes linéaires dans un milieu stratifié visqueux et avec du vent. Cette méthode linéaire est validée par des solutions quasi-analytiques reposant sur les équations de dispersion dans une atmosphère stratifiée. Elle est aussi appliquée à deux cas d’études : la propagation d’ondes liée à l’impact d’une météorite à la surface de Mars (mission NASA INSIGHT), et la propagation d’ondes atmosphériques liées au tsunami de Sumatra en 2004. La seconde méthode résout la propagation non-linéaire d’ondes gravito-acoustiques dans une atmosphère couplée, avec topographie, à la propagation d’ondes élastiques dans un solide visco-élastique. Cette méthode repose sur sur le couplage d’une formulation en éléments finis discontinus, pour résoudre les équations de Navier-Stokes la partie fluide, par éléments finis continus pour résoudre les équations de l’élastodynamique dans la partie solide. Elle a été validée grâce à des solutions analytiques ainsi que par des comparaisons avec les résultats de la méthode par différences finies. / This thesis deals with the wave propagation problem within the Earth-ocean-atmosphere coupled system. A good understanding of the these phenomena has a major importance for seismic and atmospheric explosion studies, especially for planetary missions. Atmospheric wave-forms generated by explosions or surface oscillations can bring valuable information about the source mechanism or the properties of the various propagation media. We develop two new numerical full-wave high-order modeling tools to model the propagation of acoustic and gravity waves in realistic atmospheres. The first one relies on a high-order staggered finite difference method and focus only on the atmosphere. It enables the simultaneous propagation of linear acoustic and gravity waves in stratified viscous and windy atmosphere. This method is validated against quasi-analytical solutions based on the dispersion equations for a stratified atmosphere. It has also been employed to investigate two cases : the atmospheric propagation generated by a meteor impact on Mars for the INSIGHT NASA mission and for the study of tsunami-induced acoutic and gravity waves following the 2004 Sumatra tsunami. The second numerical method resolves the non-linear acoustic and gravity wave propagation in a realistic atmosphere coupled, with topography, to the elastic wave propagation in a visco-elastic solid. This numerical tool relies on a discontinuous Galerkin method to solve the full Navier-Stokes equations in the fluid domain and a continuous Galerkin method to solve the elastodynamics equations in the solid domain. It is validated against analytical solutions and numerical results provided by the finite-difference method.

Dans ce manuscrit, nous nous inéressons à l'étude du mouvement des vagues soumises uniquement à leur poids par le biais d'équations asymptotiques. Nous commençons par rappeler la dérivation des principaux modèles généralement utilisés (Boussinesq, Green-Naghdi,...). Nous introduisons également un nouveau modèle exprimé en amplitude-flux qui correspond à une variante des équations de Nwogu. Dans le second chapitre, nous démontrons un résultat d'existence en temps long pour ces nouvelles équations et nous étudions l'existence d'ondes solitaires pour les équations de Boussinesq. Ce travail permet notamment de calculer avec une grande précision ces solutions exactes. Le troisième chapitre détaille les différences non linéaires que l'on retrouve entre les différentes équations de Boussinesq (modèles en flux-amplitude comparés aux modèles en vitesse-amplitude). Enfin, les deux derniers chapitres introduisent un nouveau paradigme pour trouver des schémas numériques adaptés aux modèles asymptotiques. L'idée est d'appliquer une analyse asymptotique aux équations d'Euler discrétisées. Ce nouveau paradigme est appliqué aux équations de Peregrine, de Nwogu et de Green-Naghdi. Plusieurs cas tests sont proposés dans ces deux chapitres. / In this work, we are interested in the evolution of water waves under the gravity force using asymptotics models. We start by recalling the derivation of most used models (Boussinesq, Green-Naghdi,...) and we introduce a new model expressed amplitude-flux, which is an alternative version of the Nwogu equations. In the second chapter, we prove a long time existence result for the new model and we investigate the existence of solitary waves for the Boussinesq models. This work allow us to compute these solutions with a good precision. The third chapter highlights the nonlinear differences between the Boussinesq equations (amplitude-flux models versus amplitude-velocity models). Finally, the two last chapter introduce a new paradigm in order to find numerical schemes adapted to asymptotics models. The idea is to apply an asymptotic analysis to a discretized Euler system. This new paradigm is applied to Peregrine equations, Nwogu equations and Green-Naghdi equations. Test cases are presented in these two chapters

La résolution analytique du problème de la stabilité d’une suspension solide (particules rigides de forme sphérique) dans le système de Taylor-Couette cylindrique a été menée. On s’est basé sur les travaux de Ali and Lueptow (2002) pour formuler les équations régissant la stabilité de l’écoulement dans le cadre d’une théorie linéaire. Ces équations sont valables dans l’approximation du faible espace annulaire et ont pour but la prévision de l’instabilité primaire. A cet effet, nous avons utilisé une méthode variationnelle telle que la méthode de Galerkin pour résoudre le problème aux valeurs propres conduisant à établir le diagramme de stabilité liée au nombre d’onde au voisinage de l’état critique du développement de la première instabilité. Tout d’abord, on a cherché à mettre au point les calculs dans le cas de l’écoulement de Taylor-Couette classique en se référant aux travaux de Ali and Lueptow (2002). Ensuite on a procédé à la résolution systématique des équations du mouvement et l’on cherche à déterminer le critère d’apparition des instabilités en présence de particules en suspension et l’on détermine simultanément les paramètres de couplage entre forces d’interaction liquide-solide. L’ensemble des travaux ainsi réalisés permettront de lever la contradiction fondamentale entre la théorie et l’expérience. L’étude expérimentale a permis d’analyser les effets de la concentration des particules en suspension (disques) et du rapport d’aspect radial ’ sur l’apparition des instabilités dans le système de Taylor-Dean. Le dispositif expérimental utilisé est constitué d’un cylindre intérieur tournant et le cylindre extérieur est maintenu fixe. Le comportement rhéologique du fluide utilisé est viscoplastique obéissant au modèle de Herschel Bulkley. L’apparition des instabilités est examinée en utilisant une technique de visualisation. Pour une concentration donnée dans l’intervalle étudié, la nature des structures apparaissant dans le système d’écoulement dépend ’, alors que pour une valeur donnée de ’ dans l’intervalle étudié, la valeur du nombre de Taylor critique Tac dépend de la concentration des particules. Nous avons obtenu que le nombre de Taylor critique Tac correspondant au déclenchement de la première instabilité évolue non linéairement en fonction de ’. De plus, nous avons examiné expérimentalement les effets de limitation axiale (effet de bords) sur le déclenchement des instabilités dans le système de Taylor-Dean. Les résultats obtenus montrent que les bords tournants n’affectent pas le type de structures qui apparaissent dans le système d’écoulement de Taylor-Dean. Cependant, ils influencent le seuil critique d’apparition des instabilités qui est marquée par des valeurs élevées du nombre de Taylor critique pour des bords tournants ce qui indique un effet stabilisant des bords mobiles. / The analytical solution of the stability problem of a solid suspension (rigid spherical particles) in the system of cylindrical Taylor-Couette was conducted. We are based on the work of Ali and Lueptow (2002) to formulate the equations governing the stability of the flow in a linear theory. These equations are valid in the approximation of small gap configuration and aim to predict the primary instability. For this purpose, we used a variational method such as the Galerkin method to solve the eigenvalue problem leading to establish the stability diagram related to the wave number in the vicinity of the critical state of development of the first instability. First, we develop the calculations in the case of Taylor-Couette flow with reference to classic work of Ali and Lueptow (2002). Then, we carried out a systematic solution of the equations of motion and we search to determine the criterion of onset of instabilities in the presence of suspended particles and coupling parameters are simultaneously determined from liquid-solid interaction force. All work carried out and will remove the fundamental contradiction between theory and experiment. The experimental study has analyzed the effect of the concentration of suspended particles (disks) and radial aspect ratio ' on the occurrence of instabilities in the Taylor-Dean flow system. The experimental device used consists of a rotating inner cylinder and the outer cylinder is stationary. The rheological behavior of the fluid is viscoplastic obeying to Herschel Bulkley model. The onset of instability is examined using a visualization technique. For a given concentration in the range studied, the nature of the structure appearing in the flow system depends on ', while for a given value of ' in the range studied, the value of the critical Taylor number Tac depends on the particle concentration. We obtain that the critical Taylor number Tac corresponding to the onset of the first instability evolves nonlinearly versus '. In addition, we examined experimentally the effect of axial limitation (endwall effects) on the onset of instabilities in the Taylor-Dean flow system. The results show that the rotating endwalls do not affect the type of structures that appear in the Taylor-Dean flow system. However, they influence the threshold of appearance of instabilities which is characterized by high values of the critical Taylor number for rotating endwalls indicating a stabilizing effect of the rotating endwalls.

La simulation numérique adaptative d'écoulements présentant des phénomènes multi-échelles s'effectue généralement au moyen d'une hiérarchie de modèles différents selon l'échelle mise en jeu et le niveau de précision requis. Ce type de modélisation numérique entraîne des problèmes de couplage multi-échelles complexes. Cette thèse est ainsi dédiée au développement, à l'analyse et à la mise en œuvre de méthodes performantes permettant de résoudre des problèmes de couplage en espace de modèles décrits par des systèmes d'équations aux dérivées partielles hyperboliques.Dans une première partie, nous développons et analysons une méthode numérique dédiée au couplage interfacial des équations d'Euler mono-dimensionnelles. Chacun des systèmes de lois de conservation est muni d'une loi de pression distincte et l'interface de couplage séparant ces modèles est supposée fixe et infiniment mince. Les conditions de transmission sont modélisées par un terme source mesure localisé à l'interface de couplage. Le poids associé à cette mesure modélise les pertes de conservation à l'interface (typiquement des pertes de charge) et sa définition permet l'application de plusieurs stratégies de couplage. Notre méthode d'approximation repose sur les techniques d'approximation par relaxation de type Suliciu. La résolution exacte du problème de Riemann pour le système relaxé nous permet de définir un schéma numérique équilibre pour le modèle de couplage. Ce schéma préserve certaines solutions stationnaires du modèle de couplage et est applicable pour des lois de pression générales. L'implémentation de notre méthode permet de mener des expériences numériques illustrant les propriétés de notre schéma. Par exemple, nous montrons qu'il est possible de contrôler l'écoulement à l'interface de couplage en calculant des poids solutions de problèmes d'optimisation sous contraintes.La deuxième partie de cette thèse est dédiée au développement de deux schémas numériques d'ordre arbitrairement élevé en espace pour l'approximation des solutions stationnaires du problème mixte associé au modèle de Jin et Xin. Nos schémas d'approximation reposent sur la méthode de Galerkin discontinue. L’approximation des solutions du problème mixte par notre premier schéma fait intervenir uniquement des erreurs de discrétisation tandis que notre deuxième schéma est constitué à la fois d'erreurs de modélisation et de discrétisation. L'erreur de modélisation provient du remplacement, dans certaines régions spatiales, de la résolution du modèle de relaxation par celle de l'équation scalaire équilibre associée. Sous l'hypothèse d'une interface de couplage éventuellement caractéristique, la résolution du problème de Riemann associé au modèle couplé nous permet de construire un schéma numérique d'ordre arbitrairement élevé prenant en compte l'éventuelle existence de couches limites à l'interface de couplage. Enfin, la mise en œuvre de ces méthodes nous permet d'analyser quantitativement et qualitativement les erreurs de modélisation et de discrétisation commises lors de l'utilisation du schéma couplé. Ces erreurs sont fonction du niveau de raffinement de maillage utilisé, du degré de polynôme choisi et de la position de l'interface de couplage. / The adaptive numerical simulation of multiscale flows is generally carried out by means of a hierarchy of different models according to the specific scale into play and the level of precision required. This kind of numerical modeling involves complex multiscale coupling problems. This thesis is thus devoted to the development, analysis and implementation of efficient methods for solving coupling problems involving hyperbolic models.In a first part, we develop and analyze a coupling algorithm for one-dimensional Euler systems. Each system of conservation laws is closed with a different pressure law and the coupling interface separating these models is assumed fix and thin. The transmission conditions linking the systems are modelled thanks to a measure source term concentrated at the coupling interface. The weight associated to this measure models the losses of conservation and its definition allows the application of several coupling strategies. Our method is based on Suliciu's relaxation approach. The exact resolution of the Riemann problem associated to the relaxed system allows us to design an extremely accurate scheme for the coupling model. This scheme preserves equilibrium solutions of the coupled problem and can be used for general pressure laws. Several numerical experiments assess the performances of our scheme. For instance, we show that it is possible to control the flow at the coupling interface when solving constrained optimization problems for the weights.In the second part of this manuscript we design two high order numerical schemes based on the discontinuous Galerkin method for the approximation of the initial-boundary value problem associated to Jin and Xin's model. Our first scheme involves only discretization errors whereas the second approximation involves both modeling and discretization errors. Indeed in the second approximation, we replace in some regions the resolution of the relaxation model by the resolution of its associated scalar equilibrium equation. Under the assumption of a possible characteristic coupling interface, we exactly solve the Riemann problem associated to the coupled model. This resolution allows us to design a high order numerical scheme which captures the possible boundary layers at the coupling interface. Finally, the implementation of our methods enables us to analyze quantitatively and qualitatively the modeling and discretization errors involved in the coupled scheme. These errors are functions of the mesh size, the degree of the polynomial approximation and the position of the coupling interface.

Ce présent ouvrage s'intéresse à la quantification de l'impact de différents paramètres sur les performances en isolation vibratoire d'un actionneur piézoélectrique THUNDER. Une récapitulation des travaux de recherche dans le domaine est tout d'abord proposée. Plus spécifiquement, la modélisation et le procédé de fabrication de l'actionneur sont décrits. Un modèle semi-analytique Rayleigh-Ritz à une base de décomposition par champ de déplacements, résolu à l'aide de MATLAB, est alors conçu. Pour ce faire, l'énergie potentielle et l'énergie cinétique totale du système sont déterminées en posant l'hypothèse de Kirchhoff pour un composite. Le modèle est par la suite validé à l'aide d'un modèle éléments finis Abaqus. Il est alors constaté qu'il est difficile de représenter les discontinuités associées à l'extension métallique de l'actionneur à l'aide d'une unique base de décomposition. Un modèle à trois bases de décomposition par champ de déplacements, résolu à l'aide du logiciel MAPLE, est donc conçu et validé à l'aide du modèle éléments finis précédemment considéré. Finalement, une quantification de l'impact de différents paramètres tels que la longueur de l'extension métallique, les conditions limites et le rayon de courbure sur le comportement fréquentiel d'un THUNDER est proposée. Il est observé qu'à l'aide du modèle Rayleigh-Ritz, il est possible d'optimiser les performances en isolation vibratoire de l'actionneur.

Dans cette thèse, nous nous intéressons à la résolution des <br />équations de Maxwell en régime fréquentiel, afin de calculer<br />précisément la signature radar de cibles diverses. Pour avoir<br />une grande précision nécessaire pour des expérience de grande taille,<br /> nous utilisons des méthodes d'ordre élevé.<br /><br />Dans le cas scalaire, les éléments finis spectraux hexaédriques<br />avec condensation de masse, permettent d'obtenir un produit matrice vecteur <br />rapide et peux coûteux en stockage. Dans le cas vectoriel, les hexaèdres<br />de la première famille ne réalisent pas la condensation de masse, mais on peut<br />écrire un algorithme rapide de produit matrice-vecteur. Des résultats<br />numériques 3-D montrent la performance de l'algorithme proposé.<br /><br />Nous traitons également le cas où la géométrie présente<br />une symétrie de révolution. On est alors ramenés à une succession<br />de problèmes 2-D indépendants.<br />Nous proposons une méthode éléments finis d'ordre élevé <br />couplée à des équations intégrales d'ordre élevé.

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.