Le "climat" à l'intérieur d'une serre dépend de son aération. Le processus d'aération est complexe, il participe à l'essentiel des échanges de chaleur et de masse avec l'extérieur, et sa maîtrise permet donc de contrôler les paramètres physiques tels que la température, l'humidité, ou les concentrations de gaz comme le CO2 par exemple. Ce contrôle est essentiel pour maintenir les plantes dans des conditions métaboliques favorables (respiration, photosynthèse, transpiration) et dans un état sanitaire satisfaisant.<br />La ventilation naturelle est le système le plus économique pour réguler le microclimat interne de la serre. Néanmoins, elle n'offre qu'un contrôle limité sur l'écoulement d'air dans la serre et reste difficile à maîtriser.<br />Cette étude contribue à l'analyse et à la modélisation des phénomènes mis en jeu dans l'aération naturelle des serres en verre, de grande portée, habituellement utilisées en culture ornementale (plantes en pots), sous climat tempéré, tel qu'en Anjou. Deux approches complémentaires incluant expérimentation in situ et modélisation mathématique du climat distribué sont mises en œuvre.<br />Des campagnes de mesures ont été menées à l'intérieur d'une serre de production et dans son environnement immédiat sous conditions réelles de culture ornementale. Des données météorologiques : température de l'air, vitesse et direction du vent, rayonnement solaire et atmosphérique, ont été collectées. L'ensemble de ces mesures constitue un jeu de données conséquent destiné à fournir les entrées du modèle numérique. Parallèlement à ces mesures, nous avons systématiquement procédé à des mesures du taux de renouvellement d'air qui ont été utilisées pour valider le modèle.<br />Un modèle numérique a été mis en œuvre. Il s'appuie sur un code de mécanique des fluides numérique (Computational Fluid Dynamics). Ce code permet de prédire les champs de vitesses et de températures à l'intérieur de la serre après résolution numérique des équations de base qui régissent les mouvements d'air (équations de Navier-Stokes couplées à l'équation de l'énergie) dans le domaine de calcul considéré. La turbulence, dont l'effet est loin d'être négligeable sous serre, a été modélisée à l'aide d'une fermeture de type k-e. Le taux d'aération a pu être déduit ensuite par résolution d'une équation de transport d'un gaz traceur virtuel. Un module radiatif a été ajouté dans le modèle numérique afin de prendre en compte le rayonnement d'origine solaire et atmosphérique. Ce module résout l'équation des Transferts Radiatifs qui est couplée à l'équation de l'énergie.<br />Ce modèle « complet » a pu être vérifié et validé pour différentes conditions climatiques. Il a été ensuite utilisé pour analyser l'impact de la configuration des ouvrants sur le climat et sur les flux de chaleur au niveau de la toiture de la serre. Cette analyse a porté non seulement sur la ventilation mais aussi sur l'homogénéité de la distribution des vitesses et des températures dans la serre et notamment au niveau des cultures.<br />Enfin, des indicateurs d'efficacité de l'aération de la serre sous climat estival ont pu être dégagés pour différentes configurations d'aération (ouverture) et différentes conditions climatiques.

Ce travail de thèse est consacré à l'étude des conditions d'échauffement des éléments de structure situés à l'extérieur des bâtiments en feu. Une synthèse bibliographique sur les feux de compartiment et plus particulièrement sur la problématique de la sortie de flammes par les ouvertures des locaux en feu a été réalisée. <br /><br />Un outil de simulation numérique avancé a été mis en oeuvre pour étudier les paramètres principaux affectant les caractéristiques des flammes extérieures. Cette tâche a nécessité au préalable la prise en main et le développement d'un modèle CFD de simulation numérique du feu, ainsi que sa validation sur des cas expérimentaux de complexité croissante. A l'étude numérique de la sortie de flammes d'un local en feu est venue s'ajouter une partie expérimentale concernant l'effet du vent sur les flammes extérieures. <br /><br />L'étude des flammes extérieures a permis de procéder à la vérification et à l'extension d'une méthode de calcul simplifiée destinée à prédire les conditions d'échauffement des éléments de structure situés à l'extérieur des bâtiments. L'application pratique de ce travail a consisté à regrouper dans un nouveau modèle global les résultats obtenus à partir de la synthèse bibliographique et des études numériques et expérimentales.

Ce travail porte sur la simulation numérique des structures magnétiques dans les matériaux ferromagnétiques. Nous nous sommes intéressés au mouvement de parois et de lignes de Bloch en utilisant une approche bidimensionnelle. Le problème initial, équations de Landau-Lifshitz-Gilbert (LLG), est tridimensionnel. Il fallait donc simplifier le problème et le formuler en dimension deux avant de faire la simulation. Dans un premier temps, nous avons résolu numériquement le modèle simplifié des équations de (LLG) donné par Slonczewski. Nous avons calculé le champ démagnétisant dans le cas d'une paroi sans structure. Des études locales des ingularités du champ au niveau de la paroi ont été menées. Une formule de quadrature pour approcher le champ est aussi proposée. L'étude est ensuite généralisée au cas d'une paroi avec structure. La discrétisation du modèle de Slonczewski nous a conduit à étudier entre autres l'évolution d'une courbe sous l'influence de sa courbure. Un schéma semi-implicite est proposé et une méthode d'approximation de la courbure pour une courbe définie par un nombre fini de points est détaillée et analysée. Nous avons également développé un algorithme pour la résolution d'une équation aux dérivées partielles posée sur une courbe en mouvement. Dans un second temps, nous avons proposé un nouveau modèle simplifié des équations de (LLG). Il est donné par deux équations couplant la vitesse normale d'une courbe (centre de la paroi) et l'angle azimutal. Contrairement à Slonczewski, nous avons tenu compte des interactions entre lignes de Bloch. Nous avons terminé l'exposé par la validation des différentes discrétisations et approximations, des résultats de la physique sont obtenus: domaines en labyrinthes, translation d'une bulle, annihilation de deux lignes déroulantes, déplacement de deux lignes enroulantes. La différence entre les deux modèles est aussi dégagée

Ce mémoire porte sur la simulation numérique de la réponse des capteurs de pression capacitifs micro-électroniques. Les courbes de réponse du capteur capacitif classique, en fonction d'une pression constante, sont établies. L'auteur examine ensuite les modifications apportées à ces courbes lorsqu'une pression électrostatique est superposée à la pression constante. La suite de ce travail est consacrée à l'étude de capteurs capacitifs de structures plus complexes. Les capteurs à armature déformable d'épaisseur variable sont traités. L'auteur montre l'influence de la variation de l'épaisseur sur la réponse du capteur. Il présente ensuite une étude des capteurs de géométrie polygonale et propose deux applications permettant d'élargir la gamme de pressions. Une analyse du comportement des capteurs en régime de contact est effectuée : la sensibilité des courbes de réponse est nettement améliorée par ce mode de fonctionnement. L'auteur propose finalement une modélisation des capteurs capacitifs en utilisant la théorie non linéaire des fortes déflexions de von Karman. Il est indispensable de tenir compte de cette théorie dès que la déflexion de la plaque sensible est de l'ordre de grandeur de son épaisseur. La simulation des problèmes développés dans ce mémoire fait intervenir des équations et des sytèmes d'équations aux dérivées partielles, linéaires et non linéaires, à deux variables d'espace. La discrétisation de ces équations est effectuée par la méthode générale des différences finies. La résolution des équations résultantes est obtenue par les méthodes itératives du type gradient conjugué et Newton-Raphson. Un logiciel convivial de discrétisation des équations aux dérivées partielles stationnaires linéaires ou non linéaires a été développé. L'auteur présente finalement une étude de la méthode des différences finies curvilignes et de son application à la deflexion des plaques.

On propose une méthode générale de modélisation numérique des équations aux dérivées partielles à l'aide de la méthode des éléments finis. On expose son application aux équations de la physique des milieux continus. On construit un système informatique expert mettant en œuvre ces résultats, on le teste et on l'applique à plusieurs études : phénomènes magnétiques, phénomènes thermodynamiques, équations de Navier-Stokes

.

L'hydrodynamique radiative est le domaine dans lequel gaz et rayonnement interagissent dynamiquement. Ses champs d'application sont très vastes, de l'astrophysique à la fusion par confinement inertiel.<br /><br />Lors de cette thèse, un code numérique parallèle tridimensionnel d'hydrodynamique radiative baptisé HERACLES a été développé. Il s'appuie sur le modèle M1 qui permet de traiter un rayonnement à forte anisotropie. De nombreux tests ont validé HERACLES sur une grande variété de conditions physiques, dont le régime semi-transparent et la diffusion des photons, ses résultats étant comparables aux codes Monte-Carlo. Il a ensuite été utilisé dans deux thématiques.<br /> <br />La première concerne les chocs radiatifs, phénomènes astrophysiques reproduits sur Terre grâce aux lasers de puissance. HERACLES a mis en évidence l'influence de différents paramètres sur l'évolution d'un tel choc : le rapport de la largeur du canal de propagation sur le libre parcours moyen de photons, l'albédo des parois... Il a ensuite contribué à l'analyse d'une expérience réalisée avec le laser PALS. Il a permis de reproduire la courbe de décélération du précurseur observée dans l'expérience ainsi que la transmission du diagnostic transverse.<br /><br />La seconde thématique concerne les jets générés par les étoiles en formation et interagissant avec le nuage moléculaire environnant. Les opacités du milieu interstellaire montrant qu'une partie significative du rayonnement est absorbée, nous avons mené les premières simulations d'un jet tenant compte du transfert radiatif. Elles ont montré que le jet pouvait être fortement comprimé et que le transfert radiatif semble donc pouvoir jouer un rôle important dans sa propagation.

Le massif Alpin est régulièrement soumis à l'automne à des systèmes dépressionnaires amenant un flux d'air méditerranéen chaud et humide sur son versant sud . De par la forme et la topographie du massif, ce phénomène peut conduire localement à des précipitations intenses souvent à l'origine d'inondations éclairs, dommageables pour l'homme. L'un des objectifs de l'expérience internationale MAP (Mesoscale Alpine Program) est l'étude des mécanismes dynamiques et microphysiques de ces précipitations, dans le but d'en améliorer la prévision. La campagne s'est déroulée durant l'automne 1999 et a donné lieu à 17 Périodes d'Observations Intensives (POI).<br />Les travaux de recherche effectués durant cette thèse sont basés sur l'analyse de plusieurs d'entre elles, ceci au moyen de données de radars Doppler sole et aéroportés, et du modèle numérique méso-échelle non-hydrostatique Méso-NH. Ces travaux visent à établir le scénario et les processus fondamentaux qui ont mené à leur occurrence. Différents types d'évènements pluvieux ont été étudiés : 2 cas de précipitations convectives, le premier à fort cumul de pluie dans la région du Lac Majeur (ouest des Alpes, POI 2B), et un autre sur la région du Frioul (est des Alpes, POI 5) ; puis des conditions pour lesquelles les précipitations se sont révélées plus modérées sur la région du Lac Majeur (POIs 4, 8, 9 et 15) permettant d'étudier un autre mécanisme lié au relief ou non : le processus de blocage.<br />Sur ces 6 cas d'étude, le facteur dominant responsable du cumul de pluie au sol est le facteur durée : la persistance des précipitations associée au ralentissement du système dépressionnaire au niveau des Alpes mène aux plus forts cumuls. Ceux-ci ont été enregistrés pour les POIs 2B et 8, soit des évènements convectif et stratiforme respectivement. Dans le cas de la POI 2B, c'est un système quasi-stationnaire d'altitude qui est responsable du ralentissement, alors que pour la POI 8 il s'agit de la nature stable du flux incident sur les flancs alpins, et par la suite de son blocage au niveau du relief. Les POIs 4 et 15, caractérisées également par la stabilité du flux d'alimentation des précipitations, illustrent des acs pour lesquels l'arrivée du front froid associé au système dépressionnaire, et par conséquent d'un flux post-frontal antagoniste au flux d'alimentation, est responsable du déplacement des cellules précipitantes. La POI5 est une situation caractérisant le passage rapide du front froid, et l'évolution d'une ligne de cellules convectives associée.<br />L'effet du relief a été mis en évidence, dans les épisodes pluvieux stratiforme et convectif. Dans tous les cas, ce dernier contribue à renforcer/déclencher les précipitations. Il est responsable de la déviation vers l'ouest des flux incidents stables. L'impact de l'orographie de petite échelle a également été démontré à l'aide de tests numériques, en particulier les irrégularités du relief de l'ordre du kilomètre ont joué un rôle important dans la progression du système précipitant de la POI 5.

Cette thèse concerne la fissuration dynamique des matériaux à gradient de propriétés. Plus particulièrement, elle traite de la tenue des gaines de combustible, tubes minces et élancées, à fort taux de combustion lors de transitoires accidentels. La fissuration est étudiée à l'échelle locale et basée sur la notion de zone cohésive, laquelle est mise en oeuvre à l'aide d'une approche multicorps. Afin de prendre en compte une fissuration en mode mixte, les modèles de zone cohésive considérés sont couplés à du contact frottant. L'approche Non-Smooth Contact Dynamics est mise en oeuvre pour le traitement dynamique des problèmes de contact frottant non réguliers. Compte tenu des dimensions mises en jeu, une démarche multiéchelle permet de prendre en compte la microstructure dans un calcul de structure. A l'échelle microscopique, les propriétés effectives surfaciques, entre les corps, sont obtenues par homogénéisation numérique périodique. Une formulation à deux champs du problème Eléments Finis multicorps est ainsi obtenu. La méthode NSCD est alors étendue à cette formulation. L'outil de simulation numérique développé permet des simulations dynamiques en grandes transformations depuis l'amorçage de fissures jusqu'à la ruine du matériau. La mise en oeuvre à l'échelle microscopique est effectuée par calcul sur des cellules de bases pour plusieurs tirages de microstructures aléatoires. A l'échelle macroscopique, des calculs dynamiques de structures sur des portions de gaine en fonction du taux moyen et du gradient d'hydrogène mettent en évidence le rôle prépondérant de ces deux paramètres dans la tenue de la gaine sous chargement dynamique rapide.

La méthode des éléments finis, bien qu'employée dans la grande majorité des codes de simulation industriels, présente un certain nombre de limitations pour la simulation des procédés de mise en forme, notamment lors de grandes déformations, de la localisation des phénomènes comme c'est le cas pour la formation de bandes de cisaillement ou de fissures. La nécessité de reconstruire un maillage vérifiant un certain nombre de critères de qualité entraîne des coûts de calcul additionnels, ainsi que des problèmes de robustesse, particulièrement pour les problèmes tridimensionnels complexes. Depuis une dizaine d'années, de nouvelles méthodes numériques, alternatives à la méthode des éléments finis, ont été développées. Ces méthodes, appelées méthodes sans maillage, construisent une partie ou la totalité de l'approximation entièrement à partir du nuage de nœuds. L'objectif de cette thèse est de développer un certain nombre d'outils numériques basés sur les concepts précédemment cités affin de mettre au point un logiciel de simulation des procédés permettant de surmonter certaines des difficultés liées à la méthode des éléments finis. La technique proposée, baptisée méthode C-NEM, est une extension de la méthode des éléments naturels (NEM) dont les améliorations permettent un traitement plus simple pour l'analyse des procédés. La première partie présente la méthode et son intérêt pour les problèmes présentant des discontinuités fixes ou mobiles. La deuxième partie présente des approches permettant de traiter les problèmes en grandes déformation avec éventuellement une localisation des phénomènes. La dernière partie illustre la méthode au travers d'exemples numériques.