Le modèle autorégressif dans un espace de Banach (ARB) permet<br />de représenter des processus à temps continu. Nous considérons<br />l'estimation de l'opérateur d'autocorrelation d'un ARB(1). Les<br />méthodes classiques d'estimation (maximum de vraisemblance et<br />moindres carrées) s'avèrent inadéquates quand l'espace<br />paramétrique est de dimension infinie, Grenander (1983} a proposé<br />d'estimer le paramètre sur un sous espace de dimension en général<br />finie m, puis d'étudier la consistance de cet estimateur lorsque<br />la dimension m tend vers l'infini avec le nombres d'observations<br />à vitesse convenable. Cette méthode est dite méthode des cribles.<br />Notons que plus généralement il serait possible d'utiliser la<br />méthode des f-divergences. Nous définissons la méthode des<br />moindres carrées comme problème d'optimisation dans un espace de<br />Banach dans le cas ou l'opérateur est p-sommable,<br />p>1. Nous montrons la convergence de l'estimateur<br />crible et sa normalité asymptotique dans le cas d'un opérateur est<br />strictement -intégral. Nous utilisons la représentation duale<br />de la f-divergence pour définir l'estimateur du minimum des<br />f-divergences. Nous nous limitons ici à l'étude de<br />l'estimateur dit du minimum de KL-divergence (divergence de<br />Kullback-Leibler). Cet estimateur est celui<br /> du maximum de vraisemblance. Nous montrons par la suite qu'il<br /> converge presque surement vers la vraie valeur du paramètre<br />pour la norme des opérateurs p-sommables. La démonstration est<br />basée sur les techniques de Geman et Hwang (1982), utilisées pour<br />des observations indépendantes et identiquement distribuées, qu'on<br />a adapté au cas autorégressif.

Cette thèse a pour objet l'analyse et la discrétisation numérique des équations tridimensionnelles de l'électromagnétisme. Ces travaux débutent par l'étude de ces équations dans un domaine b orné multiplement connexe. Un théorème d'existence général a été établi, en proposant une nouvelle approche du problème, en le reformulant à l'aide d'un opérateur approprié, tenant compte des omplexités géométriques du domaine. Dans la suite, après avoir donnée un résultat de régularité, on propose une approximation numérique de la solution par une méthode spectrale. La méthode est, d'une part, analysée numériquement dans le cas d'une décomposition conforme du domaine, et d'autre, implantée dans le cadre d'un code 3D. Des tests numériques illustrant les prévisions théoriques sont exposés et comparés à ceux obtenus par une méthode d'éléments finis de type P1 qu'on présentera sommairement. En outre, les quatre premières valeurs propres du problème discret sont calculées et comparées à celui du spectre exact. La dernière partie de cette thèse est consacrée à l'étude d'une décomposition de domaine par une méthode spectrale avec joints pour le problème de Maxwell. Il est utile de souligner que les paramètres physiques sont considérés dans cette partie comme pouvant être hétérogènes. On applique cette méthode à un problème type présenté. Ce dernier permet d'unifier deux approches qui habituellement sont distinguées: le problème d'évolution de Maxwell, et le problème de Maxwell en régime harmonique. Des estimations d'erreurs sont démontrées, elles reposent sur un lemme, qui est une variante du second lemme de Strang, permettant de décomposer l'erreur en la somme de trois erreurs principales: l'erreur sur la meilleure approximation, l'erreur de consistance et l'erreur d'intégration numérique. Cette dernière étant obtenue de ma ière classique, les deux autres erreurs ont nécessité une recherche plus approfondie, notamment, la définition d'opérateurs discrets et un Lemme d'augmentation de degré pour l'erreur sur la meilleure approximation. Enfin des courbes d'erreurs et des tests numériques sont exposés validant un code de calcul tridimensionnel développé pour l'approximation de la solution du problème type (pour des paramètres physiques homogènes et hétérogènes).

Le but de ce travail était de fournir un cadre cohérent pour le traitement des problèmes de contact en utilisant une discrétisation de type noeud à segment. Trois aspects principaux de la mécanique numérique du contact ont été particulièrement considérés : la description de la géométrie, le problème de détection de contact et les techniques de résolution. Le manuscrit contient cependant une présentation complète de la mécanique du contact et des algorithmes numériques qui lui sont attachés. Un nouveau formalisme mathématique -- les s-structures -- est employé dans l'ensemble de la thèse. Il fournit un cadre de formulation intrinsèque qui permet d'exprimer de façon compacte un grand nombre de problèmes de mécanique et de physique. La thèse propose plusieurs idées originales et des extensions des techniques classiques, qui ont toutes été mises en œuvre dans le code de calcul par éléments finis ZéBuLoN (ZSeT). Plusieurs études de cas, présentées dans la thèse, viennent démontrer les performances et la robustesse des méthodes numériques utilisées pour la détection et la résolution.

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.

Les systèmes de santé accordent une attention croissante au point de vue des usagers dans l'organisation de l'offre de soins. L'instauration d'une offre de soins sensible aux besoins et préférences des patients constitue un enjeu majeur de qualité et d'efficacité des soins. Ce travail analyse le point de vue des patients pour l'offre de soins de médecine générale en utilisant différentes méthodes permettant d'obtenir des informations complémentaires en termes d'expérience de soins, de satisfaction, d'importance ou encore de préférences. Il s'agit des méthodes Delphi, de classement du meilleur au pire et de révélation des préférences par les choix discrets. Ces méthodes sont appliquées sur deux échantillons : en population générale pour la première et chez des patients poly-pathologiques pour les 2 autres afin d'identifier les principaux enjeux actuels et à venir de la réorganisation de l'offre de soins de médecine générale du point de vue des patients. Les résultats montrent le rôle central de la relation médecin-patient et plus particulièrement de l'échange d'informations entre le médecin et le patient. Cependant une relation médecin-patient de qualité ne doit pas pour autant être réalisée au détriment de la qualité technique du soin et de la coordination de la prise en charge du patient. Ce travail montre également l'importance de prendre en compte l'expérience de soins des patients lors de l'analyse de leur point de vue, et plus particulièrement de leur disposition au changement. L'évaluation systématique et régulière des préférences des patients en pratique quotidienne peut permettre d'améliorer la communication médecin-patient ainsi que le contenu de l'offre de soins du point de vue des patients.

L'objet de cette thèse a été de développer un code de calcul pour la simulation d'écoulements diphasiques de fluides non miscibles, incompressibles et isothermes afin de l'appliquer au phénomène de fragmentation d'un jet liquide, et plus particulièrement à la désintégration assistée d'une nappe de liquide par deux écoulements d'air portés à haute vitesse.<br /><br />Les choix des hypothèses physiques et des schémas numériques associés ont été fait afin de respecter au mieux la physique complexe de brisure d'un jet liquide. La résolution des équations de Navier-Stokes incompressibles est faite de façon directe par le biais d'une méthode de projection. La méthode naissante et prometteuse Level-Set assure quant à elle le suivi numérique de la surface de séparation entre deux fluides non miscibles. Enfin, la méthode Ghost Fluid permet un traitement correct des discontinuités à la traversée de l'interface en préservant au niveau discret les conditions de saut entre équilibre des forces de capillarité, de pression et de viscosité. Les bonnes aptitudes de tels schémas numériques sont montrées à travers une batterie de cas tests académiques.<br /><br />Les mécanismes physiques mis en jeu lors de l'atomisation primaire d'un jet liquide ont fait l'objet de nombreuses études théoriques et expérimentales. Les instabilités se développant à la surface du liquide sont multiples et clairement tridimensionnelles. Actuellement et malgré beaucoup d'efforts de recherche, aucune théorie ou modèle ne rend compte rigoureusement de ce phénomène. Or, il est montré que l'outil numérique peut apporter une lumière nouvelle. Par exemple, l'impact de la couche limite gazeuse sur la fréquence d'oscillation est étudié. Nous avons également pu retrouver par la simulation la dynamique ligamentaire intrinsèque à la désintégration assistée d'un jet liquide.

Les méthodes de décomposition de domaine en espace ont prouvé leur utilité dans le cadre des architectures parallèles. Pour les problèmes d'évolution en temps, il est nécessaire d'introduire une dimension supplémentaire de parallélisme dans la direction du temps. Ceci peut alors être couplé avec des méthodes de type optimisé Schwarz waveform relaxation. Nous nous intéressons dans cette thèse aux méthodes directes de décomposition en temps. Nous en étudions particulièrement deux. Dans une première partie nous étudions la méthode de produit tensoriel, introduite par R. E. Lynch, J. R. Rice, et D. H. Thomas in 1963. Nous proposons une méthode d'optimisation des pas de temps, basée sur une étude d'erreur en variable de Fourier en temps. Nous menons cette étude sur les schémas d'Euler et de Newmark pour la discrétisation en temps de l'équation de la chaleur. Nous présentons ensuite des tests numériques établissant la validité de cette approche. Dans la seconde partie, nous étudions les méthodes dites de Bloc, introduites par Amodio et Brugnano en 1997. Nous comparons diverses implémentations de la méthode, basées sur différentes approximations de l'exponentielle de matrice. Nous traitons l'équation de la chaleur et l'équation des ondes, et montrons par une étude numérique bidimensionnelle la puissance de la méthode.

Aujourd'hui le monde des télécommunications est en plein essor et le nombre de services proposés aux consommateurs augmente d'année en année. Les technologies employées font appel à l'optique. Le simple constat du nombre d'applications "sans fils" qui se développent permet de se rendre compte de l'importance des microondes dans l'ensemble des technologies modernes de communication. Cette thèse se situe dans ce cadre. Elle présente des techniques de modélisation de base pour concevoir et optimiser un guide d'onde de surface fonctionnant dans le domaine des microondes. Dans le premier chapitre, des outils de simulation permettant de calculer la réponse de diffraction par les réseaux de strips métalliques ont été développés. Les méthodes rigoureuses qui ont été retenues sont la méthode MMFE, CBCM et la méthode C qui prennent en compte la nature vectorielle de la lumière. La difficulté de convergence des calculs numériques et le problème de discontinuité des champs sur les strips métalliques est mise en évidence. Pour traiter ces problèmes, un changement de coordonnées est proposé, c'est le système de " Coordonnées adaptative " qui permet de resserrer les lignes de coordonnées au voisinage des points des discontinuités. Il en résulte une diminution du saut de discontinuité en ces points et une amélioration de la convergence des calculs numériques. Dans le second chapitre, ces méthodes sont étendues au problème de la diffraction en incidence oblique encore appelé diffraction conique. Dans le troisième chapitre, nous avons appliqué ces méthodes pour étudier et caractériser les réseaux de strips métalliques déposés sur une couche diélectrique. Nous avons mis en évidence le phénomène de résonance sur les facteurs de réflexion et nous avons pu montré que ces effets des résonances sont dus au couplage de l'onde de surface et de l'onde plane incidente. Afin d'analyser les couplages résultants qui existent entre le mode, nous avons étudié les pôles de la matrice S et déterminé la sensibilité et l'influence des paramètres optoélectroniques sur le pic de résonance. Cette étude a permis de déterminer le triplet (hauteur, facteur de forme, l'angle de polarisation ) relatif à la structure pour que le guide soit optimisé.

Lorsqu'une onde acoustique se propage dans un milieu granulaire, elle est susceptible de provoquer la mobilité des grains, aussi infime soit-elle. Inversement, la mobilité d'un grain dans une matrice fluide peut induire un champ acoustique et dans les deux cas, l'énergie acoustique peut être transférée à la fois au travers des pores et des contacts entre grains. Nous avons mis au point un modèle original permettant de considérer ces deux modes de transfert d'énergie pour simuler la propagation d'ondes acoustiques dans les milieux granulaires immergés. Dans le cas des milieux granulaires secs, l'inertie du fluide est telle que l'énergie transférée dans l'air peut être négligée et le milieu modélisé avec des algorithmes de type "dynamique moléculaire". Au contraire, dans le cas de milieux immergés, l'énergie portée par le fluide ne peut pas être négligée et nous montrons que la méthode des domaines fictifs basée sur les multiplicateurs de Lagrange distribués permet de coupler les équations de la dynamique et l'équation d'onde. Nous utilisons la méthode des éléments finis pour propager l'onde dans le fluide, les grains étant modélisés en 2D par des sphères rigides et incompressibles afin de satisfaire les hypothèses de l'algorithme de dynamique moléculaire. Les résultats du modèle sur des expériences numériques simples mais pour lesquelles existent des solutions analytiques de l'acoustique mettent en évidence la validité du nouveau modèle. Nous en donnons une illustration pour l'étude des interactions subies par un empilement réaliste de multiples grains mobiles soumis à un signal acoustique.

Décrit comme l'un des algorithmes les plus prometteurs du 20ème siècle, le formalisme multipolaire appliqué à la méthode des éléments de frontière, permet de nos jours de traiter de larges problèmes encore inconcevables il y a quelques années. La motivation de ce travail de thèse est d'évaluer la capacité, ainsi que les avantages concernant les ressources numériques, de ce formalisme pour apporter une solution de référence aux problèmes de propagation sonore tri-dimensionnels en environnement urbain, dans l'objectif d'améliorer les algorithmes plus rapides déjà existants. Nous présentons la théorie nécessaire à l'obtention de l'équation intégrale de frontière pour la résolution de problèmes non bornés. Nous discutons également de l'équation intégrale de frontière conventionnelle et hyper-singulière pour traiter les artefacts numériques liés aux fréquences fictives, lorsque l'on résout des problèmes extérieurs. Nous présentons par la suite un bref aperçu historique et technique du formalisme multipolaire rapide et des outils mathématiques requis pour représenter la solution élémentaire de l'équation de Helmholtz. Nous décrivons les principales étapes, d'un point de vue numérique, du calcul multipolaire. Un problème de propagation sonore dans un quartier, composé de 5 bâtiments, nous a permis de mettre en évidence des problèmes d'instabilités dans le calcul par récursion des matrices de translations, se traduisant par des discontinuités sur le champs de pression de surface et une non convergence du solveur. Ceci nous a conduits à considérer le travail très récent de Gumerov et Duraiswamy en lien avec un processus récursif stable pour le calcul des coefficients des matrices de rotation. Cette version améliorée a ensuite été testée avec succès sur un cas de multi diffraction jusqu'à une taille dimensionnelle de problème de 207 longueur d'ondes. Nous effectuons finalement une comparaison entre un algorithme d'élément de frontière, Micado3D, un algorithme multipolaire et un algorithme basé sur le tir de rayons, Icare, pour le calcul de niveaux de pression moyennés dans une cour ouverte et fermée. L'algorithme multipolaire permet de valider les résultats obtenus par tir de rayons dans la cour ouverte jusqu'à 300 Hz (i.e. 100 longueur d'ondes), tandis que concernant la cour fermée, zone très sensible par l'absence de contribution directes ou réfléchies, des études complémentaires sur le préconditionnement de la matrice semblent requises afin de s'assurer de la pertinence des résultats obtenus à l'aide de solveurs itératifs / Described as one of the best ten algorithms of the 20th century, the fast multipole formalism applied to the boundary element method allows to handle large problems which were inconceivable only a few years ago. Thus, the motivation of the present work is to assess the ability, as well as the benefits in term of computational resources provided by the application of this formalism to the boundary element method, for solving sound propagation problems and providing reference solutions, in three dimensional dense urban environments, in the aim of assessing or improving fast engineering tools. We first introduce the mathematical background required for the derivation of the boundary integral equation, for solving sound propagation problems in unbounded domains. We discuss the conventional and hyper-singular boundary integral equation to overcome the numerical artifact of fictitious eigen-frequencies, when solving exterior problems. We then make a brief historical and technical overview of the fast multipole principle and introduce the mathematical tools required to expand the elementary solution of the Helmholtz equation and describe the main steps, from a numerical viewpoint, of fast multipole calculations. A sound propagation problem in a city block made of 5 buildings allows us to highlight instabilities in the recursive computation of translation matrices, resulting in discontinuities of the surface pressure and a no convergence of the iterative solver. This observation leads us to consider the very recent work of Gumerov & Duraiswamy, related to a ``stable'' recursive computation of rotation matrices coefficients in the RCR decomposition. This new improved algorithm has been subsequently assessed successfully on a multi scattering problem up to a dimensionless domain size equal to 207 wavelengths. We finally performed comparisons between a BEM algorithm, extit{Micado3D}, the FMBEM algorithm and a ray tracing algorithm, Icare, for the calculation of averaged pressure levels in an opened and closed court yards. The fast multipole algorithm allowed to validate the results computed with Icare in the opened court yard up to 300 Hz corresponding, (i.e. 100 wavelengths), while in the closed court yard, a very sensitive area without direct or reflective fields, further investigations related to the preconditioning seem required to ensure reliable solutions provided by iterative solver based algorithms