On s'intéresse dans cette thèse à la modélisation et à la simulation multidimensionnelle de l'hydrodynamique fluviale, notamment près des estuaires. Le modèle physique de référence est le modèle 3D, mais au vu de son important coût de calcul, il est intéressant de disposer de modèles plus simples en 1D, 2D ou 2.5D, que l'on peut utiliser dans des zones adéquates du fleuve, en fonction de sa bathymétrie. Ainsi, à partir du modèle 3D basé sur les équations instationnaires et incompressibles de Navier-Stokes, des modèles plus simples sont dérivés par projection par formulations faibles du problème 3D. On obtient ainsi un modèle en 1D, écrit sur la courbe médiane de la surface libre du fleuve, ainsi que deux modèles en 2D, le 2D-vertical écrit sur la surface longitudinale médiane du fleuve et le 2D-horizontal écrit sur la surface libre. Enfin on définit un modèle en quasi-3D, le modèle 2.5D, écrit dans la somme des espaces 2D-vertical et 2D-horizontal. Tous ces modèles prennent en compte la géométrie du fleuve et fournissent une vitesse tridimensionnelle ainsi que la pression, qui n'est pas supposée hydrostatique mais qui est une inconnue entière du problème. En outre, on définit et justifie un estimateur de modèles entre le modèle 3D et chacune de ses approximations en 1D, 2D et 2.5D. Cet estimateur calcule l'erreur entre le modèle 3D et son approximation, et donne ainsi une indication sur la qualité des résultats obtenus à partir des modèles 1D, 2D ou 2.5D, dans leurs zones respectives de calcul. Tous ces modèles hydrodynamiques sont implémentés dans des codes d'éléments finis, écrits en C++. Enfin, ils sont couplés numériquement à l'aide de l'estimateur de modèles.

Les travaux présentés dans cette thèse traitent de la modélisation par la méthode des éléments finis de configurations de contrôle non destructif par courants de Foucault. Un code de calcul programmé en langage C++, s'appuyant sur la discrétisation en trois dimensions des équations de Maxwell en régime harmonique a été développé. Deux formulations magnétodynamiques duales en potentiels combinés "a − psi" et "t −phi " ont été mises en oeuvre. Des éléments de Whitney étraédriques du premier ordre ont été utilisés. Les potentiels scalaires sont discrétisés aux noeuds des éléments et les potentiels vecteurs à leurs arêtes. Le travail présenté dans ce mémoire aborde un cas précis mais de grande importance parmi les applications du CND par CF. Il s'agit des défauts de faible ouverture (fissures), qui sont fréquemment rencontrés en CND. Leur détection permet de prévenir la destruction des pièces en fonctionnement et d'augmenter la fiabilité des produits industriels. Leur prise en compte par la MEF peut être délicate car leur maillage conduit à une forte densité d'éléments ou à des éléments déformés. C'est dans ce contexte qu'un modèle a été développé pour faciliter la modélisation de ce type de défauts. L'idée consiste à considérer la fissure comme une surface non conductrice imperméable au courant. Pour assurer le bon comportement des différentes grandeurs électromagnétiques, des conditions aux limites appropriées sont appliquées sur la surface de la fissure. Ces conditions sont prises en compte par la formulation électrique a− psi en dédoublant les degrés de liberté (potentiel scalaire électrique) attachés aux noeuds de la surface de la fissure de part et d'autre de celle-ci. Par la formulation magnétique t− phi en annulant la circulation du potentiel vecteur électrique t sur les arêtes appartenant à la surface de celle-ci. La résolution simultanée des deux formulations a−psi et t−phi par la MEF permet de vérifier au sens fort toutes les équations de Maxwell. Les résultats obtenus se complètent et les erreurs numériques dues aux discrétisations se traduisent sous forme d'une non-vérification des lois de comportement. Cela a permis, en utilisant cette propriété de complémentarité de deux formulations, de définir des indicateurs d'erreur afin de développer une procédure d'adaptation automatique de maillage. Plusieurs types d'estimateurs d'erreur ont été définis et étudiés sur différents cas test.

Dans cette thèse nous nous intéressons à l'analyse d'erreur a priori et a posteriori de méthodes d'éléments finis mixtes et non-conformes. Nous considérons en particulier les équations de Darcy à perméabilité fortement variable et les équations de convection-diffusion-réaction en régime de convection dominante. Nous discrétisons les équations de Darcy par une méthode d'éléments finis mixtes non-conformes de type Petrov-Galerkin appelée schéma boîte. Les techniques d'estimations d'erreur a posteriori par résidu et hiérarchique conduisent à des estimateurs d'erreur a posteriori fiables et optimaux indépendamment des fluctuations de la perméabilité. Les résultats théoriques sont validés numériquement sur différents cas tests présentant de forts contrastes de perméabilité. Enfin, nous montrons comment les indicateurs d'erreur obtenus permettent de générer des maillages adaptatifs. Nous discrétisons les équations de convection-diffusion-réaction par des éléments finis nonconformes. Deux méthodes de stabilisation sont étudiées: la stabilisation par viscosité de sous-maille, conduisant à un schéma boîte et la méthode de pénalisation sur les faces. Nous montrons que les deux schémas ainsi obtenus ont les mêmes propriétés de convergence que les approximations par éléments finis conformes. Grâce aux techniques d'estimations d'erreur par résidu nous obtenons des estimateurs d'erreur a posteriori fiables et optimaux. Certains des indicateurs d'erreur sont robustes au sens de Verfürth, c'est à dire que le rapport des constantes intervenant dans les inégalités de fiabilité et d'optimalité explose en au plus l'inverse du nombre de Péclet. Les résultats théoriques sont validés numériquement et les indicateurs d'erreur a posteriori obtenus permettent de générer des maillages adaptatifs sur des problèmes présentant des couches intérieures.