L'objectif de ce travail de thèse est de prendre en compte à la fois l'hétérogénéité, l'anisotropie et des incertitudes dans la simulation 3D de la propagation d'ondes élastiques. Pour ce faire, dans un premier temps, on modélise le champ de propriétés mécaniques, ici le champ de tenseur d'élasticité, par un modèle de champ stochastique 3D des matrices définie-positives. La construction de ce modèle de champ est essentiellement fondée sur celle de Soize [2008]. Notre modèle conserve ainsi les propriétés principales du modèle de Soize comme le paramétrage minimal contrôlant l'amplitude de la fluctuation et la taille caractéristique de la variabilité spatiale, le comportement local a priori arbitrairement anisotrope (anisotropie triclinique) et les propriétés mathématiques fondamentales. De plus, un nouveau paramètre est introduit dans ce modèle pour imposer un niveau d'anisotropie moyen souhaité. Dans un deuxième temps, on effectue des adaptations du code de calcul d'éléments finis spectraux, à savoir le code parallèle SPEC3D, afin d'une part de générer les réalisations du champ stochastique du tenseur d'élasticité et d'autre part de prendre en compte l'anisotropie dans la résolution numérique du problème élastodynamique. Des études paramétriques utilisant SPEC3D sont ensuite réalisées mettant en évidence les influences de l'anisotropie et des paramètres d'hétérogénéité sur la propagation d'ondes sismiques. En particulier, elles démontrent une dépendance directe entre la longueur de corrélation du champ de propriétés et le temps caractéristique d'apparition de la diffusion. Ce régime se manifeste par l'équipartition d'énergie entre les mouvements irrotationnels et rotationnels.

L'objectif de ce travail de thèse est de prendre en compte à la fois l'hétérogénéité, l'anisotropie et des incertitudes dans la simulation 3D de la propagation d'ondes élastiques. Pour ce faire, dans un premier temps, on modélise le champ de propriétés mécaniques, ici le champ de tenseur d'élasticité, par un modèle de champ stochastique 3D des matrices définie-positives. La construction de ce modèle de champ est essentiellement fondée sur celle de Soize [2008]. Notre modèle conserve ainsi les propriétés principales du modèle de Soize comme le paramétrage minimal contrôlant l'amplitude de la fluctuation et la taille caractéristique de la variabilité spatiale, le comportement local a priori arbitrairement anisotrope (anisotropie triclinique) et les propriétés mathématiques fondamentales. De plus, un nouveau paramètre est introduit dans ce modèle pour imposer un niveau d'anisotropie moyen souhaité. Dans un deuxième temps, on effectue des adaptations du code de calcul d'éléments finis spectraux, à savoir le code parallèle SPEC3D, afin d'une part de générer les réalisations du champ stochastique du tenseur d'élasticité et d'autre part de prendre en compte l'anisotropie dans la résolution numérique du problème élastodynamique. Des études paramétriques utilisant SPEC3D sont ensuite réalisées mettant en évidence les influences de l'anisotropie et des paramètres d'hétérogénéité sur la propagation d'ondes sismiques. En particulier, elles démontrent une dépendance directe entre la longueur de corrélation du champ de propriétés et le temps caractéristique d'apparition de la diffusion. Ce régime se manifeste par l'équipartition d'énergie entre les mouvements irrotationnels et rotationnels.

L'effet des champs magnétiques sur les processus de transfert de matière peut se manifester essentiellement du fait de l'existence de deux champs de forces : les forces électromagnétiques qui résultent de l'interaction du champ magnétique B avec la densité de courant j selon : FL=jxB. L'autre champ de force a été mis en évidence beaucoup plus récemment. Il s'agit d'une force agissant sur les espèces paramagnétiques en solution : , avec χm la susceptibilité magnétique molaire, μ0 la perméabilité magnétique et C la concentration d'espèce paramagnétique). Cette expression, qui résulte d'une approche énergétique. Le travail présenté a pour objectif d'une part de mettre en évidence l'existence des forces d'origine magnétique, et d'autre part de quantifier leur influence sur les processus de transfert de matière. Le travail est basé sur une double approche expérimentale et théorique. Sur le plan expérimental une cellule électrochimique a été réalisée dans une configuration aussi simple que possible pour permettre une analyse qui s'écarte des géométries souvent adoptées utilisant des électrodes circulaires. Cette cellule est de forme parallélépipédique ce qui permet de mieux contrôler les distributions de courant électrique et d'imposer un champ magnétique parallèle à la densité de courant, éliminant au maximum les forces électromagnétiques. Les expériences réalisées permettent de mesurer les densités de courant locales et globales et de suivre leurs évolutions en fonction du champ magnétique. Ces mesures sont complétées par l'enregistrement du champ de vitesse effectué sous champ pour la première fois, à l'aide de sonde à ultrason. Sur le plan théorique, le travail se fonde sur la modélisation numérique des écoulements (CFD model). Cette approche a été réalisée en deux dimensions en supposant par ailleurs le régime permanent, ce qui a limité l'approche à des champs magnétiques modérés. Au delà de certaines valeurs, des instabilités se manifestent, ne permettant pas une exploitation du modèle. Des travaux ultérieurs seront très utiles pour une validation totale. Par ailleurs, une approche phénoménologique des processus a permis de retrouver les lois d'évolution des transferts de matière en B2/3 identifiées par la plupart des expérimentateurs. La confrontation des résultats théoriques et expérimentaux révèle un excellent accord, même si cette comparaison est limitée aux faibles valeurs du champ pour les raisons déjà énoncées, ce qui permet finalement de confirmer l'existence du champ de force magnétique comme élément susceptible de contrôler les transferts en électrochimie.

Parmi les diverses entités que l'architecte est amené à manipuler au cours du processus de conception, les « ambiances » constituent un domaine éminemment complexe, constitué des interactions que l'usager perçoit entre les phénomènes physiques et les formes de l'environnement construit. La lumière naturelle s'inscrit précisément dans ce cadre, où la variabilité individuelle vient s'ajouter à l'aspect dynamique des phénomènes. Or s'il est manifeste que les architectes s'attachent de manière permanente au travail sur la forme, ils ne semblent pas toujours traiter avec la même rigueur la disponibilité en lumière naturelle et les interactions diffuses. La complexité des notions abordées et le manque d'outils de simulation en sont sans doute les raisons majeures. Nous proposons donc, en nous appuyant sur des travaux préalables du laboratoire CERMA de l'École d'Architecture de Nantes, le développement d'un nouvel outil logiciel susceptible d'évaluer de manière performante la lumière naturelle au sein du projet architectural et urbain, à partir d'une maquette numérique tridimensionnelle de l'espace étudié. Si certains travaux concernant l'évaluation de la lumière naturelle dans les espaces intérieurs existent déjà, en revanche les espaces extérieurs restent, en dehors d'aspects relatifs à l'ensoleillement, rarement envisagés ; notre premier objectif est donc de considérer indifféremment espaces intérieurs et espaces extérieurs, et de les soumettre de manière globale aux conditions d'éclairement qui résultent de l'ensemble constitué du soleil et de la voûte céleste. Celles-ci sont basées sur des modèles de distribution statistiques. Notre perception de la lumière provenant exclusivement de sa réflexion sur les matériaux, ceux-ci jouent un rôle essentiel dans les échanges entre les différentes surfaces de l'environnement. Le second objectif est donc de prendre en compte ces réflexions complexes, au travers de vitrages éventuellement. Il est atteint en adoptant un traitement informatique s'inspirant de la théorie des échanges radiatifs établie par les thermiciens il y a plus d'un demi-siècle, reprise et instrumentée dans des recherches plus récentes (années 80) sur l'image numérique. En dernier lieu sont proposés quelques exemples d'application dans le domaine lumineux, mettant en évidence les potentialités d'analyse de l'outil développé ; ils sont complétés par les perspectives qu'offrent les modalités de traitement utilisées, notamment dans le champ de la micro-climatologie urbaine.

Ce travail utilise le calcul numérique pour prendre en compte les déformations plastiques en mécanique de la rupture. On peut ainsi vérifier les résultats de Rice concernant l'étendue de la zone plastique, en tenant compte de l'écrouissage. On étudie ensuite des problèmes liés à la propagation des fissures sous chargement monotone d'une part, sous chargement cyclique d'autre part. Dans le troisième chapitre, on vérifie l'estimation théorique de la force due aux singularités en fond de fissures. Ceci permet de mettre en évidence également une réduction de l'ouverture par rapport au cas élastique, qui explique le phénomène de fermeture des bords libres des fissures en fatigue. Un dernier chapitre applique les résultats précédents à la fatigue.

La récupération in situ (ISR) est une méthode d'exploitation minière qui consiste à extraire le minerai in situ par le biais de solutions d'attaque acides ou alcalines. Comme l'ISR a lieu en sous-sol, la gestion de l'ISR reste globalement empirique. La modélisation numérique de cette technique minière permettrait d'optimiser la conduite de ce type d'exploitation. Trois types de phénomènes doivent être pris en compte pour la modélisation de la lixiviation de l'uranium par ISR: (1) les réactions géochimiques; (2) les cinétiques de ces réactions et (3) l'influence du transport hydrodynamique sur ces réactions. Des tests de lessivage ont été menés sur des échantillons du site minier de Tortkuduk (Kazakhstan), où l'ISR de l'uranium est réalisée par voie acide. Deux types de tests de lessivage ont été menés en laboratoire : (1) des essais en batchs; (2) des essais sur colonnes. Les hypothèses concernant les mécanismes de lessivage ont été testées et validées en modélisant ces tests avec les codes de calcul CHESS et HYTEC, développés à Mines ParisTech. Un modèle 1D couplant géochimie et transport du procédé ISR à l'échelle du laboratoire a été proposé. Il permet de traduire la séquence des mises en solutions en une séquence de réactions géochimiques. Les facteurs contrôlant la dissolution de l'uranium et les précipitations secondaires ont été mis en évidence ainsi que leur importance relative.

L'étoile beta Pictoris est à l'heure actuelle la seule étoile de la séquence principale autour de laquelle un disque de poussières a été visualise (1984). L'observation de son spectre visible et ultra-violet a révélé la présence de raies en absorption dues à la partie gazeuse du disque. L'étude de cette partie gazeuse a révélé que les raies observées variaient dans le temps: parfois, des composantes additionnelles en absorption, décalées par effet doppler vers les grandes longueurs d'onde y apparaissent. Ces phénomènes sporadiques ont été reliées à la chute vers l'étoile de petit corps de type comètaire vaporisant des éléments métalliques. La dynamique d'ions métalliques évaporés d'un corps parent a été ensuite etudiée, et une modélisation numerique du comportement d'un nuage de ces ions a été entreprise. La simulation a montré que les observations étaient reproduites par le modèle. De plus, la différence de comportement entre les raies visibles et uv a été expliquée. Des contraintes déduites de la simulation ont alors permis de formuler une hypothèse selon laquelle une planète située dans le disque pourrait par perturbations être responsable de la chute d'objets vers l'étoile. Un modèle fondé sur cette hypothèse a été ensuite développé. La simulation qui en a résulté a permis de montrer que si l'orbite de la planète est suffisamment elliptique, des perturbations sont capables d'envoyer un grand nombre de petits objets vers l'étoile. Il est ressorti de cette étude que l'évolution globale, constatée sur quelques années, des raies de beta Pictoris apparaissait comme une conséquence naturelle du modèle. En poursuivant l'étude théorique de ce disque, nous saurons peut-être s'il constitue le premier exemple de système planétaire découvert en dehors du système solaire

Dans une optique de maîtrise du procédé d'incinération des ordures ménagères et de ses possibles émissions polluantes, nous avons développé un modèle mathématique qui simule un lit d'ordures ménagères en combustion sur une grille mobile. Ce modèle décrit la plupart des phénomènes physicochimiques et thermiques intervenant lors de l'incinération : séchage et pyrolyse de la charge, combustion et gazéification du carbone résiduel, transferts thermiques, effondrement du lit, brassage... Il intègre également une description des mécanismes de volatilisation des métaux lourds et de formation des NOx. La cinétique de départ des métaux lourds est modélisée en tenant compte des différentes étapes de transport (transfert externe, diffusion intraparticulaire, volatilisation) au moyen de l'approche des temps caractéristiques additifs. Dans le cas simulé du cadmium, la prédiction d'une volatilisation quasi-complète est conforme aux résultats de la littérature. Le sous-modèle NOx prend en compte les mécanismes de formation thermique, prompt, combustible, par l'intermédiaire de N2O, ainsi que les mécanismes de réduction homogène par recombustion et hétérogène par le carbone résiduel. Les calculs révèlent que prédominent la formation par le mécanisme combustible et la destruction par la réduction hétérogène. Enfin, le modèle de lit a été utilisé pour étudier l'influence des divers paramètres opératoires : température, débit et distribution d'air primaire, taille des particules de déchets, brassage et schéma de brassage. Les résultats sont présentés et discutés en détail. L'influence des conditions opératoires sur l'efficacité du procédé et sur les émissions de Cd et NOx est analysée.

Cette contribution à l'étude du comportement dynamique du procédé de rectification plane passe profonde est axée sur l'influence des paramètres d'usinage (formes géométriques, états de surface

La distillation est un procédé industriel de séparation de phases qui fait typiquement intervenir un écoulement diphasique caractérisé par un film liquide laminaire ou faiblement turbulent s'écoulant par gravité et cisaillé à contre-courant par un courant gazeux turbulent. Afin de comprendre la dynamique de ce genre d'écoulements, nous avons développé un modèle numérique de simulation d'écoulements diphasiques prenant en compte la présence éventuelle des structures turbulentes. Ce modèle s'appuie sur un couplage entre les méthodologies Volume of Fluid sans étape de reconstruction pour le suivi d'interface et la simulation des grandes échelles pour le traitement de la turbulence. Les contraintes de sous-maille sont évaluées par une approche dynamique mixte, ce qui permet au modèle de s'adapter aux caractéristiques locales de la turbulence et de fonctionner même dans des zones laminaires. Le modèle développé est ensuite testé en simulant différentes configuration d'écoulements de films liquides cisaillés ou non par un courant gazeux. / Distillation is an industrial process of phase separation which involves a two-phase flow characterized by a laminar or weakly turbulent gravity- riven liquid film sheared by a countercurrent turbulent gas stream. To understand the dynamics of such flows, we developed a numerical technique aimed at computing incompressible turbulent two-phase flows. A large eddy simulation (LES) approach based on a dynamic mixed model is used to compute turbulence while the two-phase nature of the flow is described through a Volume of Fluid (VOF) approach with no interface reconstruction step. The use of a dynamic mixed approach for modelling the subgrid stresses allows the developed model to self-adapt to local characteristics of turbulence, so that it also works in laminar flows. The whole methodology is then applied to the computation of different configurations of liquid films sheared or not by a gas stream.