Le SAGD (Steam Assisted Gravity Drainage) est un procédé de récupération des huiles lourdes qui remporte énormément de succès, en particulier pour le bitume. Le SAGD génère des interactions complexes entre la géomécanique et les écoulements polyphasiques en milieux poreux. Dans ce procédé, l'injection de vapeur modifie la pression et de la température dans le réservoir, ce qui peut augmenter ou diminuer les contraintes effectives dans le réservoir. La quantification de l'état de contrainte et déformation dans le réservoir est essentielle pour effectuer un bon pronostic de la productivité du réservoir, pour vérifier l'intégrité de la couverture et les risques de fracturation hydraulique, et également pour interpréter correctement la sismique 4D en termes d'évolution de la chambre de vapeur. Dans le procédé SAGD, les effets géomécaniques de l'injection de vapeur dans le réservoir sont liés aux écoulements de fluide.. Le couplage réservoir-géomécanique est un sujet de recherche important. Pour effectuer ce type de simulation, une solution consiste à utiliser un simulateur en éléments finis pour décrire la géomécanique et un simulateur en volumes finis pour décrire les écoulements. Dans cette thèse, une simulation couplée thermo-hydro-mécanique du SAGD a été effectuée à l'aide du simulateur de réservoir PumaFlow et du simulateur de géomécanique Abaqus. Les principaux thèmes étudiés dans cette étude ont été (1) la stratégie de couplage, (2) la géométrie du système et (3) le type de maillage utilisé. Ce travail a été effectué sur des cas synthétiques.

Les nano commutateurs magnétiques sont en cours d'étude et ils sont susceptibles d'être présents dans plusieurs domaines d'application comme les interrupteurs d'alimentation, les convertisseurs DC/DC...etc. Compte tenu de leur très petite taille et de leur procédé de fabrication, leurs performances sont sujettes à de grandes dispersions. Partant du besoin de modélisation et d'optimisation avec fiabilité de ces dispositifs, ces travaux de thèse se décomposent en trois axes. Axe modélisation : Développement d'un modèle semi analytique pour calculer la déformation avec l'analyse de contact mécanique. Ce modèle a été introduit pour le couplage magnéto-mécanique dans le logiciel de modélisation de MEMS magnétiques MacMMems. Axe dérivation automatique de code : Développement de JAP (Java Jacobian Automatic Programming) qui est un outil générique de dérivation d'algorithmes. Il a été appliqué pour le nano commutateur et a été introduit dans les outils de dimensionnement développés au G2Elab pour l'analyse de sensibilité, pour l'optimisation exploitant le Jacobien et aussi pour orienter des systèmes d'équations algébro-différentiels. Axe optimisation : Mise en oeuvre du modèle et des outils développés pour faire l'étude de sensibilité et l'optimisation sous contraintes de fiabilité du nano commutateur magnétique.

Avec 70% de la surface de la planète recouverte par les mers et les océans, la présence de particules salées, produites majoritairement par le déferlement de vagues, représente une composante majeure dans le cycle géochimique de l'atmosphère et dans le bilan radiatif terrestre. Nos connaissances doivent progresser vers une description et une modélisation des flux actuels plus précises afin d'être capable d'anticiper les changements éventuels de la composition de l'atmosphère. Ce travail de thèse s'inscrit dans le cadre général d'une meilleure caractérisation des particules injectées dans l'atmosphère par le déferlement des vagues. Un des objectifs est d'etudier la validité des modèles paramétriques et numériques pour prédire la variation spatio-temporelle des particules d'aérosols à l'échelle locale et régionale. Une première partie du travail est consacrée à l'étude de la variation spatiale des concentrations d'aérosols dans une zone côtière méditerranéenne. Pour ce faire, un couplage entre un code météorologique méso-échelle (RAMS) et un modèle déterministe de distributions d'aérosols développé par le LSEET (Medex)a été réalisé. Les prédictions du couplage sont confrontées aux données enregistrées lors d'une campagne expérimentale menée en mai 2007 à la station de mesures du laboratoire sur l'île de Porquerolles. Une deuxième partie du travail s'est concentrée sur l'amélioration du modèle paramétrique Medex. A partir, de données enregistrées en mai 2007, les prédictions de Medex ont été affinées, en particulier par correction de l'influence saisonnière. Dans une troisième partie, on s'est intéressé aux paramètres qui influencent les variations de concentrations en aérosols lors d'un épisode de Mistral. Il est montré que des paramètres autres que la vitesse du vent interviennent, tel que le fetch et la hauteur de la couche limite marine. Une dernière partie est consacrée au développement d'un modèle numérique de transport des aérosols. Les premiers résultats permettent de montrer l'influence de la stabilité atmosphérique sur les profils horizontaux et verticaux d'aérosols.

Ce travail de thèse est consacré à la corrélation entre propriétés structurales à l'échelle atomique et propriétés magnétiques de multicouches (Pt/Co)3/Pt /IrMn à anisotropie d'échange perpendiculaire. Ces multicouches sont utilisées dans le domaine de l'électronique de spin. Une attention particulière a été portée à la caractérisation nanostructurale de l'interface Co/IrMn, cette interface jouant un rôle prédominant dans le phénomène d'anisotropie d'échange, caractérisée par le champ d'échange HE. L'étude structurale, réalisée principalement en sonde atomique tomographique a permis de révéler la structure des multicouches étudiées dont certaines couches sont d'épaisseur sub-nanométrique. Les résultats des mesures magnétiques ont montré que l'insertion d'une couche mince (spacer) de Pt à l'interface Co/IrMn augmente les valeurs de HE pour la multicouche dont l'épaisseur des couches de Co est tCo = 0,4 nm, mais diminue HE lorsque tCo = 0,6 nm. Les analyses nanostructurales montrent qu'en l'absence du spacer et pour tCo =0,4 nm, les atomes d'Ir et de Mn diffusent dans toute l'épaisseur de la couche de Co sur laquelle repose la couche IrMn alors que pour tCo = 0,6 nm la diffusion de ces atomes est beaucoup plus réduite. Quelle que soit l'épaisseur de la couche de Co, la présence du spacer réduit très fortement la diffusion d'Ir et de Mn, montrant ainsi que le spacer agit comme une barrière de diffusion. Pour tCo =0,4 nm, l'ajout du spacer, en limitant très fortement l'interdiffusion à l'interface entraîne l'augmentation de HE . Pour tCo =0,6 nm, la diffusion étant limitée, l'ajout d'un spacer entraîne une diminution de HE en raison de l'éloignement des spins de Co et Mn.

Le cadre général de ce travail concerne la modélisation du comportement mécanique sous impact de merlons de protection contre les chutes de blocs rocheux. Les ouvrages ciblés sont à technologie cellulaire, et sont composés de gabions à enveloppe grillagée remplis de matériaux granulaires (pierres concassées ou mélange de sable et de déchiquetas de pneus). L'absence de recommandations relatives au dimensionnement de ce type de structure de protection a conduit à la mise en place du projet ANR REMPARe couplant expérimentations in situ et modélisations numériques. Dans le cadre de ce projet, une méthode numérique innovante couplant les approches discrète et continue a été développée. Localement, la méthode des éléments discrets est employée pour conserver la précision de modélisation dans les zones fortement sollicitées proches de l'impact. Le remblai technique aval, peu sollicité, a été discrétisé par une méthode continue aux différences finies. La jonction entre ces deux zones fait l'objet d'un couplage basé sur des critères de conservation d'énergie et de condition cinématique en déplacement. Un modèle de comportement pour les gabions a été proposé et calibré à partir d'essais expérimentaux réalisés à l'échelle de la cellule. Des simulations d'impacts sur un merlon de protection prototype ont été menées puis validées par comparaison avec les mesures issues des expérimentations en vraie grandeur. Afin d'améliorer la compréhension des mécanismes engendrés, l'étude paramétrique portant sur les propriétés des matériaux et les caractéristiques d'impact a permis d'évaluer leur influence sur le mécanisme d'impact.

Le contrôle 3D de la lumière est réalisé, à l'échelle de la longueur d'onde, dans des circuits photoniques intégrés. La brique élémentaire choisie dans cette étude est le cristal photonique (CP) membranaire qui, par ses propriétés de dispersion, permet un contrôle de la lumière à la fois dans le plan (optique guidée) et hors du plan (dispositifs adressables par la surface). En particulier, l'exploitation de modes de Bloch situés au point Gamma de la courbe de dispersion (k//=0) permet l'émission de la lumière dans la direction verticale. Dans cette étude, nous nous sommes focalisés sur des cristaux photoniques à réseaux de micropiliers, comme alternative aux réseaux de trous, mais également en envisageant la possibilité d'intégrer ces structures dans des systèmes microfluidiques, les fluides ayant la capacité de circuler au travers des réseaux de piliers. Nous décrirons dans une première partie, des dispositifs à CP à émission par la surface. Nous démontrerons, pour la première fois, que l'utilisation de modes fortement résonants permet de réaliser des microlasers à réseaux de piliers en InP. Les modes faiblement résonants peuvent être utilisés pour la réalisation de miroirs à CP et de microcavités Fabry-Pérot constituées uniquement de tels miroirs. Les facteurs de qualité obtenus (>10000) rendent possible la fabrication de nouveaux types de VCSEL. Dans une seconde partie, nous nous intéresserons à la problématique de l'intégration de ces dispositifs dans un circuit photonique 3D. Tout d'abord, nous expliquerons comment il est possible d'optimisation le diagramme de rayonnement des composants. Ensuite, nous étudierons le couplage de dispositifs à cristaux photoniques avec un ou deux guides d'onde ruban en silicium. De fortes efficacités de couplage sont obtenues en simulation FDTD 3D (95%). Ces dispositifs en cours de fabrication en collaboration avec le LETI-CEA de Grenoble permettront de démontrer expérimentalement ce couplage.

L'objectif de l'étude consiste en l'établissement d'une modélisation prédictive du procédé avec des temps de calcul acceptables en vue d'optimiser le procédé vis-à-vis de la tenue mécanique de l'assemblage soudé. Le modèle mis en place est constitué de deux étapes sucessives.<br /><br />Lors de la première étape, un calcul stationnaire thermohydrodynamique couplé permet d'obtenir le champ de température et l'écoulement de matière autour de l'outil. Pour cela, un comportement purement visqueux est adopté et la zone de contact outil-matière est modélisée par des conditions de glissement. Des données expérimentales obtenues par thermographie infrarouge ont permis de fixer deux paramètres de ce modèle.<br /><br />La seconde étape du modèle s'appuie sur un comportement élastoviscoplastique pour déterminer l'état mécanique résiduel d'un assemblage. Pour cela, l'écoulement de matière permet de suivre l'histoire de chaque particule, en particulier lors de leur passage autour de l'outil. Le champ de température permet de quantifier la dissolution des précipités durcissants au sein de l'alliage et de tenir compte de la variation avec la température des paramètres de la loi de comportement adoptée. En outre, la fraction de précipités dissous est reliée à la limite d'élasticité du matériau. Une méthode de calcul stationnaire est utilisée pour obtenir directement l'état stationnaire des variables mécaniques.<br /><br />Les résultats obtenus sont en accord avec des profils de contraintes résiduelles publiés dans la littérature. De plus, l'influence de la vitesse de soudage et de la vitesse de rotation de l'outil a été étudiée au regard des contraintes résiduelles simulées.

Les instabilités paramétriques non linéaires (NL) ont été observées sur les ondes magnéto-élastiques dans le cas d'un couplage de trois quasi-phonons sous pompage électromagnétique. La théorie en prédit une dynamique supercritique explosive, mais limitée expérimentalement par le décalage de fréquence dû aux fortes nonlinéarités. La dynamique supercritique des instabilités paramétriques NL est étudiée dans deux matériaux antiferromagnétiques "plan facile" (AFEP): l'hématite α-Fe2O3 et le borate de fer FeBO3. Ces matériaux possèdent une très grande NL acoustique effective en raison du couplage magnéto-élastique élevé. Les mécanismes de limitation de la dynamique explosive ont été analysés à l'aide de l'approximation anharmonique. La compensation du décalage fréquentiel NL par une modulation de phase singulière du pompage a été proposée et théoriquement vérifiée, puis utilisée pour l'observation expérimentale de la dynamique supercritique explosive des excitations de trois quasi-phonons dans les résonateurs magnéto-élastiques. Les études sur FeBO3 ont été réalisées dans la gamme de température 77 K - 293 K où les paramètres magnéto-élastiques du cristal varient de façon significative. Un modèle fortement non linéaire des excitations de trois quasi-phonons dans les AFEPs a été développé. Les simulations numériques sont en accord avec les résultats expérimentaux. Les études théoriques de couplage de trois ondes magnéto-élastiques progressives ont été effectuées sur la base de modèles théoriques prenant en compte la non-linéarité cubique des cristaux AFEP réels. Les simulations numériques prévoient un comportement explosif et une localisation spatiale des triades générées

Ce travail constitue une première étape vers la compréhension des couplages involontaires des composant bobinés. Il se rapporte plus précisément aux inductances entrant dans un cube de 10 cm d'arête et travaillant entre 100 kHz et 10 MHz. Deux études sont particulièrement développées : l'outillage mathématique indispensable pour formuler, dans l'approximation du champ proche, le couplage magnétique de composants voisins et la métrologie nécessaire à la caractérisation de l'induction rayonnée par un composant de ce type. Quelques caractérisations de composants industriels sont présentées pour illustrer l'intérêt de ces travaux. L'induction magnétique entourant une inductance se décrit par la donnée de ses trois composantes en une multitude de points. Cette présentation est difficile à illustrer et elle ne permet pas de comparer facilement deux champs. Grâce au développement multipolaire, ces données se synthétisent, avec une précision suffisante, par 15 constantes au maximum. Habituellement, c'est le potentiel pseudo scalaire qui est développé. Malheureusement, ce potentiel n'est facile à relier ni aux courants qui le créent ni aux tensions qu'il induit. Le développement multipolaire du potentiel vecteur permet d'atteindre ce but et de prouver que les 15 composantes de l'émissivité d'un composant caractérisent aussi sa susceptibilité. Les harmoniques sphériques réelles et les harmoniques sphériques vectorielles sont des notions originales introduites par ce travail. Le banc de caractérisation déplace une antenne boucle selon les coordonnées angulaires. Un ordinateur commande ces déplacements et enregistre tes données issues d'une détection synchrone. De multiples précautions ramènent le bruit aux alentours de 30 nV, ce qui est indispensable pour mesurer des tensions de l'ordre du ^V. Grâce à ce banc, les rayonnements d'inductances électriquement équivalentes sont comparés, l'efficacité du blindage de certains composant est chiffrée et nous savons maintenant que le moment dipolaire d'un tore n'est pas axial.

Ce travail porte sur la simulation numérique tridimensionnelle par éléments finis du refroidissement primaire lors de la coulée continue d'aciers via le logiciel THERCAST. L'objectif de cette étude est de modéliser avec précision les régions solides, liquides (et pâteuses) présentes dans une lingotière de coulée continue. Dans ce contexte, la stratégie globale instationnaire est utilisée : le domaine de calcul croît à la vitesse de coulée. Afin de s'affranchir des difficultés liées aux approches purement eulériennes ou langrangiennes, nous proposons une méthode de type Eulérien-Lagrangien-ALE (Arbitrairement Lagrangienne Eulérienne). La peau solidifiée est modélisée par une loi de comportement élasto-visco-plastique. Une méthode de remaillage adaptatif automatique permet de la représenter finement grâce à un maillage anisotrope très fin dans le sens des gradients de température et de contrainte. Afin de modéliser les écoulements liquides, nous avons implémenter une méthode de transport sur la densité de bulles d'argon, couplée aux équations de Navier-Stokes par le terme de force gravitationnelle. Les effets turbulents sont pris en compte par le modèle k − " standard. Tous ces développements ont été confrontés avec succès à des résultats académiques ou à des benchmarks, bien que la simulation d'une configuration industrielle prenant en compte l'ensemble de ces méthodes et phénoménologies ne soit pas encore à notre portée.