Ce travail a consisté à développer un modèle de canopée (CIM), qui pourrait servir d’interface entre des modèles méso-échelles de calcul du climat urbain et des modèles micro-échelles de besoin énergétique du bâtiment. Le développement est présenté en conditions atmosphériques variées, avec et sans obstacles, en s’appuyant sur les théories précédemment proposées. Il a été, par exemple, montré que, pour être en cohérence avec la théorie de similitude de Monin-Obukhov, un terme correctif devait être rajouté au terme de flottabilité de la T.K.E. CIM a aussi été couplé au modèle méso-échelle WRF. Une méthodologie a été proposée pour profiter de leurs avantages respectifs (un plus résolu, l’autre intégrant des termes de transports horizontaux) et pour assurer la cohérence de leurs résultats. Ces derniers ont montré que ce système, en plus d’être plus précis que le modèle WRF à la même résolution, permettait, par l’intermédiaire de CIM, de fournir des profils plus résolus près de la surface. / This study consisted in the development of a canopy model (CIM), which could be use as an interface between meso-scale models used to simulate urban climate and micro-scale models used to evaluate building energy use. The development is based on previously proposed theories and is presented in different atmospheric conditions, with and without obstable. It has been shown, for example, that to be in coherence with the Monin-Obukhov Similarity Theory, that a correction term has to be added to the buoyancy term of the T.K.E. CIM has also been coupled with the meteorological meso-scale model WRF. A methodology was proposed to take advantage of both models (one being more resolved, the other one integrating horizontal transport terms) and to ensure a coherence of the results. Besides being more precise than the WRF model at the same resolution, this system allows, through CIM, to provide high resolved vertical profiles near the surface.

L’objectif de ce travail est d’améliorer la caractérisation des sources à l’origine d’anomalies magnétiques, par le biais de développements théoriques et expérimentaux. Pour des structures invariantes dans une direction horizontale (structures à 2D), l’idée de Nabighian (1972) a été généralisée, ce qui implique une nouvelle façon d’étudier les méthodes potentielles à 2D. Ces développements ont permis de proposer une nouvelle approche théorique et de nouvelles représentations dans le plan complexe des fonctions magnétiques. La deuxième approche de ce travail est expérimentale. Une étude y est faite par le biais de données et de cas synthétiques, afin d’estimer l’apport d’acquisitions magnétiques à différentes altitudes pour caractériser une source. Dans un troisième temps, les développements théoriques et expérimentaux sont mis en œuvre pour deux exemples : des anomalies magnétiques en domaine océanique de type Vine et Mathews et une cartographie pour l’archéologie. / Magnetic anomalies recorded outside bodies provide high quality information relative to buried structures. By using theoretical and experimental developments, this thesis aims to improve the characterization of the source inducing the magnetic anomaly. Firstly, some properties of three dimensional magnetic functions are presented, for example the presence of several maxima of the analytic signal operator (3D) is demonstrated. The Nabighian (1972) equation is generalized, which imply a new process to study potential method in two dimensions. These developments enable a new visualization of the anomaly in the complex field. The second approach is experimental. Synthetic cases are used to estimate the contribution of different configurations of magnetic data acquisitions at different altitudes to characterize the magnetic source. Theoretical and experimental developments are finally applied to two field examples: oceanic magnetic anomalies and archaeological magnetic prospection

Dans les structures en béton précontraint, les phénomènes de fluage et de retrait tendent à réduire les efforts de précontrainte initialement prévus pour maintenir le béton dans un état minimisant les forces de traction et donc la fissuration. La compréhension et la prédiction de ces phénomènes par le biais de modèles sont donc primordiales pour la conception et la maintenance à long terme des ouvrages du génie civil tels que les enceintes de confinement des centrales nucléaires.L’objectif de cette thèse est d’élaborer un cadre de modélisation micromécanique pour décrire de manière unifiée le retrait et le fluage dans les matériaux cimentaires. Pour cela, l’étude se base sur l’échelle de la microstructure poreuse du gel de C-S-H où les mécanismes intrinsèques de ces déformations différées du béton opèrent. Une approche d’homogénéisation numérique modélisant ces phénomènes dans des microstructures poreuses à morphologies quelconques est développée. Une description explicite du réseau poreux ainsi que de la phase liquide de l’eau pendant les processus de séchage/humidification est prise en compte. Les mécanismes concernant lesdéformations différées dans la phase solide sont modélisés par la théorie de la microprécontrainte-solidification (MPS). Les simulations à l’échelle microscopique sont réalisées par une approche originale couplant la méthode de Lattice Boltzmann (LBM) et la méthode des éléments finis (FEM). La LBM est utilisée pour décrire la distribution du liquide capillaire à l’échelle du pore,tandis que la FEM est employée pour simuler la déformation du squelette solide sous l’action combinée de l’eau dans l’espace poreux et d’un chargement macroscopique.La démarche proposée permet, au travers des simulations, de mieux comprendre les mécanismes liés à la non saturation et aux effets capillaires dans les milieux poreux. En particulier, la prise en compte de morphologies réalistes de microstructures et des ménisques formés conduit à différents régimes de retrait/gonflement. Ainsi les effets de l’intensité de la pression capillaire,de la tension de surface et des surfaces de chargement sur la réponse élastique du squelette solide sont évalués. Enfin, nous proposons une extension des approches précédentes au cas d’un squelette viscoélastique se déformant sous les effets de la pression capillaire et des tensions de surface. A partir des observations numériques réalisées, nous proposons un modèle pour décrire le fluage et le retrait du gel de C-S-H de manière unifiée / In pre-stressed concrete structures, creep and shrinkage tend to reduce the pre-stress forces which are initially produced so as to maintain concrete in a state minimizing traction forces and then cracks. Understanding and predicting these phenomena through models are thus highly important for the design and durability of civil engineering structures, such as containment buildings in nuclear power plants.The objective of this thesis is to develop a micromechanical modeling framework to describe shrinkage and creep in cementitious materials in a unified manner. For this purpose, the study focuses on the scale of the porous structure of the C-S-H gel where the intrinsic mechanisms of delayed strains are active. A computational homogenization approach is developed to model these phenomena in porous structures with arbitrary morphologies. An explicit description of the porous network and of the liquid phase of water during the drying/humidification process is taken into account. The mechanisms related to delayed strains in the solid phase are modeled by the microprestress-solidification theory (MPS). The simulations at the microscale are conductedbased on an original approach coupling the Lattice Boltzmann method (LBM) and the finite element method (FEM). The LBM is used to describe the distribution of capillary water in the porous structure, whereas the FEM serves as modeling the strain of the solid skeleton under the capillary water effets and a macroscopic load.The proposed method allows, by means of the simulations, to better understand the mechanisms related to the capillary effects in the porous structure. More specifically, taking into account realistic morphologies of microstructures and of the formed menisci lead to different regimes of shrinkage/swelling. Then, the effects of capillary pressure intensity, of surface tension and of morphologies of capillary surfaces on the elastic response of the solid skeleton are evaluated. Finally, the above approaches are extended to the case of a viscoelastic solid deformed under the action of the capillary water. From numerical observations, we propose a model is proposed to describe the creep and shrinkage of C-S-H gel in a unified way

Ce travail de thèse, réalisé dans le cadre des projets de recherche ROSIRIS (IRSN) et Geant4-DNA, porte sur la construction d’une simulation multi-échelle dédiée au calcul des dommages radio-induits précoces à l’ADN qui peuvent apparaître suite à l’irradiation d’un noyau cellulaire. L’outil développé s’appuie sur une version modifiée du code de Monte Carlo Geant4-DNA et est capable de simuler dans le détail le transport et les interactions physiques entre l’irradiation ionisante et la matière biologique (étape physique), la création d’espèces chimiques (étape physico-chimique) et les réactions et processus de diffusion de ces dernières (étape chimique). Durant la simulation de ces trois étapes, un modèle géométrique de l’ADN, décrivant l’ensemble du génome humain avec une précision moléculaire, est généré avec un nouveau logiciel développé dans le cadre de cette thèse : DnaFabric. Les premiers résultats obtenus pour des irradiations avec des protons et des ions alpha sont détaillés et comparés à des données de la littérature. Un bon accord est observés avec ces dernières illustrant ainsi la cohérence de l’ensemble de la simulation. L’influence très significative du critère de sélection utilisé pour identifier les dommages à l’ADN est également démontrée. / This work was performed in the frame of the ROSIRIS (IRSN) and Geant4-DNA research projects and describes the development of a simulation tool to compute radioinduced early DNA damages in a cell nucleus. The modeling tool is based on a modified version of the Monte Carlo code Geant4-DNA and is able to simulate the physical interactions between ionizing particles and the biological target (physical stage), the creation of chemical species within the cell nucleus (physico-chemical stage) as well as the reactions and diffusion processes of these chemical species (chemical stage). During all the simulation, a geometrical model that describes the DNA content of a human diploid cell nucleus is taken into account. This model was generated with a new software (DnaFabric) developed in the frame of this work and has a molecular level of detail.The first results (in term of DNA strand breaks) obtained with this tool are detailed and compared with experimental data from the literature. The good agreement between the simulation results and those data shows the coherence of our modeling. The significant influence of the selection criteria used to identify the DNA damages is also demonstrated.

Des couches minces contenant des phases intermétalliques présent et des propriétés et de combinaisons de propriétés qui ne sont que partiellement explorées. Elles portent des solutions potentielles pour conférer des multifonctionnalités à des matériaux avancés requis par les secteurs industriels et sont source derupture et de l'innovation. Leur élaboration par dépôt chimique en phase vapeur à partir de précurseurs métallo-organiques (MOCVD) permet un dépôt conforme sur, et la fonctionnalisation de surfaces complexes, avec un temps de traitement court et à un coût modéré.Pour ceci, il est nécessaire de contrôler les réactions chimiques complexes et des mécanismes de transport impliqués. La modélisation informatique du procédé, alimentée avec des données obtenues par des expériences de dépôt ciblées, est un outil intégré pour l'étude et la compréhension des phénomènes qui se produisent à différentes échelles,de l’échelle macroscopique à celle nanométrique. La MOCVD de composés intermétalliques Al-Fe est étudiée en tant que paradigme de la mise en oeuvre d'une telle approche combinée, expérimentale et théorique. La phase approximanteAl13Fe4 est particulièrement ciblée,en raison de son intérêt comme alternative peu onéreuse aux catalyseurs à base de métaux noble dans l'industrie chimique. La mise au point du dépôt de la phase Al13Fe4est subordonnée à l'étude des proc /min à 185oC. La simulation du procédé prédit des vitesses de croissance en bon accord avec ces résultats, en particulier dans la gamme 139oC-227oC. La modélisation multi-échelle prédit la rugosité RMS avec précision, permettant ainsi le contrôle des propriétés telles que la résistivité électrique. La possibilité d'obtenir des films MOCVD de Fe à faible contamination en O et C est explorée dans la gamme 130oC-250oC à partir de fer pentacarbonyle, Fe(CO)5. La morphologie de la surface des films dépend fortement de la température de dépôt; elle devient plus lisse au-dessus de200oC, qui correspond aussi à la vitesse de croissance maximale, 60nm/min. La vitesse de dépôt diminue fortement lorsque la pression augmente. Les prédictions de la modélisation macroscopique reproduisent précisément ce comportement. Elles indiquent que la diminution de la vitesse de croissance à des températures et des pressions élevées est due à l’augmentation de la décomposition du précurseur en phase gazeuse et à l'inhibition de la réactivité de surface par le ligand CO. Le modèle multi échelle conduit à des valeurs RMS en bon accord avec les mesures expérimentales, en particulier à des températures plus élevées. Suite à l’étude des deux procédés, des co-dépôts d'Al-Fe effectués à 200oCrésultent en des films riches en Alavec une microstructure poreuse et rugueuse. Ceux-ci ne contiennent pas de phases intermétalliques et sont riches en oxygène dû à la réaction d'Al avec les ligandscarbonyles. Afin d’éliminer la contamination, des dépôts séquentiels d'Al et de Fe sont réalisés, ce dernier dans des conditions modifiées à 140oC, 40Torr et 10 min. Ces films sont exemptes d’hétéroélements et présentent un rapport atomique Al:Fe13:4. Diffraction des rayons X et microscopies électroniques révèlent qu’un recuit in situ à 575oC pendant 1 h conduit à des films à gradient de la composition sur l'épaisseur,composés de la phase approximantem-Al13Fe4 conjointement avec des phases intermétalliques Al-Fesecondaires. Il est ainsi démontré que des procédés MOCVD sont appropriés pour obtenir des films constitués d'alliages intermétalliques.Ces films multifonctionnels,appliqués de façon conforme sur des surfaces complexes sont utiles pour un grand nombre d'applications. / Films containing intermetallic compounds exhibit properties and combination of properties which are only partially explored. They carry potential solutions to confer multifunctionality to advanced materials required by industrial sectors and to become a source of breakthrough and innovation.Metalorganic chemical vapor deposition (MOCVD) potentially allows conformal deposition on, and functionalization of complex surfaces, with high throughput and moderate cost. For this reason, it is necessary to control the complex chemical reactions and the transport mechanisms involved in a MOCVD process. In this perspective, computational modeling of the process, fed with experimental information from targeted deposition experiments, provides an integrated tool for the investigation and the understanding of the phenomena occurring at different length scales, from the macro- to the nanoscale. The MOCVD of Al-Fe intermetallic compounds is investigated in the present thesis as a paradigm of implementation of such a combined, experimental and theoretical approach. Processing of the approximant phase Al13Fe4 is particularly targeted, due to its potential interest as low-cost and environmentally benign alternative to noble metal catalysts in the chemical industry. The attainment of the targeted Al13Fe4 intermetallic phase passes through the investigation of the MOCVD of unary Al and Fe films. The MOCVD of Al from dimethylethylamine alane (DMEAA) in the range 139oC-241oC results in pure films. Increase of the deposition temperature yields higher film density and decreased roughness. The Aldeposition rate increases to a maximum of 15.5 nm/min at 185oC and then decreases. Macroscopic simulations of the process predictdeposition rates in sufficient agreement with experimental measurements, especially in the range 139oC-227oC. At higher temperatures, competitive gas phase and surface phenomena cannot be captured by the applied model. Multiscale modeling of the process predicts the RMS roughness of the films accurately, thus allowing the control of properties such as electrical resistivity which depend on the microstructure. The MOCVD of Fe from iron pentacarbonyl, Fe(CO)5, is investigated in the range 130oC-250oC for the possibility toobtain fairly pure Fe films with low Oand C contamination. The surface morphology depends strongly on the temperature and changes are observed above 200oC. The Fe deposition rate increases up to 200oC, to a maximum of 60 nm/min, and then decreases. Moreover, the deposition rate decreases sharply with increasing pressure. Computational predictions capture accurately the experimental behavior and they reveal that the decrease athigher temperatures and pressures is attributed to the high gas phase decomposition rate of the precursor and to inhibition of the surface fromCO. The multiscale model calculates RMS roughness in good agreement with experimental data, especially at higher temperatures. Upon investigation of the two processes, aseries of Al-Fe co-depositions performed at 200oC results in Al-rich films with a loose microstructure. They contain no intermetallic phases and they are O-contaminated due to the reaction of the Al with the carbonyl ligands. Sequential deposition of Al and Fe followed by in situ annealing at 575oC for 1 h is applied to bypass the Ocontamination. The process conditions of Fe are modified to 140oC, 40 Torr and 10 min resulting in O-free films with Al:Fe atomic ratio close to the targeted 13:4 one. Characterization techniques including X-ray diffraction, TEM and

Cette thèse est consacrée au développement d’un noyau du logiciel MEMSALab de modélisation parcalcul symbolique qui sera utilisé pour la génération automatique de modèles asymptotiques pourdes matrices de micro et nano-systèmes. Contrairement à des logiciels traditionnels réalisant des simulationsnumériques utilisant des modèles prédéfinis, le principe de fonctionnement de MEMSALabest de construire des modèles asymptotiques qui transforment des équations aux dérivées partiellesen tenant compte de leurs caractéristiques. Une méthode appelée ”par extension-combinaison” pourla modélisation asymptotique, qui permet la construction de modèle de façon incrémentale de sorteque les caractéristiques désirées soient incluses étape par étape est tout d’abord proposé pour lemodèle d’homogénéisation dérivation. Il repose sur une combinaison de méthodes asymptotiquesissues de la théorie des équations aux dérivés partielles et de techniques de réécriture issues del’informatique. Cette méthode concentre sur la dérivation de modèle pour les familles de PDEs aulieu de chacune d’entre elles. Un modèle d’homogénéisation de l’électro thermoélastique équationdéfinie dans un domaine mince multicouche est dérivé par utiliser la méthode mathématique danscette approche. Pour finir, un outil d’optimisation a été développé en combinant SIMBAD, une boite `aoutils logicielle pour l’optimisation et développée en interne, et COMSOL-MATLAB. Il a ´ et ´e appliquépour étudier la conception optimale d’une classe de sondes de microscopie atomique thermique et apermis d’ établir des règles générale pour leurs conception / This thesis is dedicated to develop a kernel of a symbolic asymptotic modeling software packageMEMSALab which will be used for automatic generation of asymptotic models for arrays of micro andnanosystems. Unlike traditional software packages aimed at numerical simulations by using pre-builtmodels, the purpose of MEMSALab is to derive asymptotic models for input equations by taking intoaccount their own features. An approach called ”by extension-combination” for the asymptotic modelingwhich allows an incremental model construction is firstly proposed for the homogenization modelderivation. It relies on a combination of the asymptotic method used in the field of partial differentialequations with term rewriting techniques coming from computer science. This approach focuses onthe model derivation for family of PDEs instead of each of them. An homogenization model of theelectrothermoelastic equation defined in a multi-layered thin domain has been derived by applyingthe mathematical method used in this approach. At last, an optimization tool has been developed bycombining a house-made optimization software package SIMBAD and COMSOL-MATLAB simulationand it has been applied for optimization of a SThM probe.

L'étude du comportement du béton en température est un problème majeur qui vise in fine à évaluer le niveau de sécurité des structures sous des sollicitations thermiques sévères, lors d'incendies par exemple. Pour cela, de nombreux modèles sont développés dans un cadre de couplage thermo-hydro-mécanique (THM), pour prendre en compte les différents processus physico-chimiques et mécaniques mobilisés par ces sollicitations et conditionnant la stabilité du matériau béton.Cependant, ces modélisations ne prennent souvent pas en compte explicitement la nature hétérogène du matériau. En effet, le béton est un matériau constitué d'inclusions noyées dans une matrice cimentaire pouvant être schématisée comme un milieu poreux ouvert partiellement saturé en eau. Les inclusions sont caractérisées par leurs natures minéralogiques, leurs morphologies et leurs tailles classées en fuseaux granulaires. Cette hétérogénéité introduit une complexité supplémentaire : la nécessité de prendre en compte la microstructure pour quantifier l'effet de l'incompatibilité (thermique, hydrique et mécanique) inclusion-matrice sur le comportement THM du béton. Ce travail constitue une première avancée dans ce sens. A ce titre, une modélisation élément fini multi-échelle tridimensionnelle (3D) est développée permettant d'affecter des comportements spécifiques à la matrice et aux inclusions. Pour la matrice, siège des transports de masse dans son réseau poreux connecté, une approche thermo-hydrique à trois fluides (eau, vapeur et air sec) est adoptée et est couplée à une loi de comportement poro-mécanique endommageable. Les inclusions (granulats) sont considérées hydriquement inertes, une approche thermo-mécanique avec endommagement est alors retenue.Une analyse, par simulations numériques, des effets de la nature minéralogique des granulats (calcaires ou silico-calcaires) de leurs distributions et de leurs morphologies a été menée sur des configurations 2D et 3D. Les effets étudiés ont notamment concerné l'influence de ces paramètres sur les fluctuations locales des champs de température, de pression de gaz et d'endommagement au regard de la dispersion des mesures expérimentales. L'analyse est limitée à l'échelle mésoscopique, celle où les granulats de taille caractéristique supérieure à 5 mm sont considérés, le reste (stable et pâte de ciment) étant une matrice homogène. Enfin, cette analyse a mis en évidence le besoin de mettre en œuvre une approche expérimentale cohérente avec une analyse mutli-échelle, à la fois pour la caractérisation des propriétés (thermiques, hydriques et mécaniques) de chaque constituant et pour l'étude des évolutions des champs lors des changements d'échelles. Un protocole expérimental a été définit et des premiers résultats d'essais sont présentés et analysés au travers de résultats obtenus dans la partie modélisation / The investigation of the behavior of heated concrete is a major research topic which concerns the assessment of safety level of structures when exposed to high temperatures, for instance during a fire. For this purpose, several modeling approaches were developed within thermo-hydro-mechanical (THM) frameworks in order to take into account the involved physic-chemical and mechanical processes that affect stability of heated concrete. However, existing models often do note account explicitly for the heterogeneity of the material : concrete is composite material that may be schematized as an assembly of inclusions (aggregates) embedded in a cementitious matrix (cement paste). This latter may be described as a partially saturated open porous medium. The aggregates are characterized by their mineralogical nature together with their morphology and size distribution. The material heterogeneity bring an additional complexity : the need to take into account the microstructure in order to quantify the effect of matrix-inclusion thermal, hygral and mechanical incompatibilities on the THM behavior of concrete. This work is a first step in this direction. For this purpose, a three-dimensional (3D) multi-scale finite element model is developed. It allows affecting specific behaviors to matrix and inclusions. For the former, where mass transports occur within the connected porous network, a three-fluids approach (liquid water, vapor and dry air) is adopted and is coupled to a poro-mechanical damage based approach. For inclusions (aggregates) no hygral component arises a pure thermo-mechanical model is considered. The developed model is then used to investigate, either by 2D or 3D numerical simulations, effects of mineralogical nature, morphology and distribution of aggregates. Studied effects have mainly concerned the influence of these parameters on local fluctuations of simulated temperature, gas pressure and damage fields with regard to experimentally observed dispersion. The analysis is here limited to the mesoscale, at which only aggregates with a characteristic size above 5 mm are meshed while the remaining inclusions together with the cement paste are considered to be a homogeneous matrix. Finally, the numerical analysis carries out the need to perform an experimental campaign that is consistent with a multi-scale approach of the THM behavior of concrete : an experimental campaign that allows to identify thermal, hygral and mechanical properties of each concrete constituent and that permit to assess evolution of fields during upscalling. An experimental protocol is then elaborated for this purpose and some obtained results are presented and analyzed with regard to results obtained in the modeling part

Le plissement dans les films minces sur un substrat plus mou a été largement observé dans la nature. Ces phénomènes ont suscité un intérêt considérable au cours de la dernière décennie. L’évolution en post-flambage d’instabilités morphologiques implique souvent de forts effets de non-linéarité géométrique, de grandes rotations, de grands déplacements, de grandes déformations, une dépendance par rapport au chemin de chargement et de multiples brisures de symétrie. En raison de ces difficultés notoires, la plupart des analyses non-linéaires de flambement ont recouru à des approches numériques parce qu’on ne peut obtenir qu’un nombre limité de solutions exactes de manière analytique. Cette thèse propose un cadre général pour étudier le problème de flambage de systèmes film/substrat de manière numérique : de la modélisation 2D ou 3D, d’un point de vue classique ou multiéchelle. L’objectif principal est d’appliquer des méthodes numériques avancées pour des analyses de bifurcations multiples aux divers modèles de film/substrat, en particulier en se concentrant sur l’évolution en post-flambement et la transition du mode à la surface. Les modèles intègrent la Méthode Asymptotique Numérique (MAN) comme une technique robuste de pilotage et des indicateurs de bifurcation qui sont bien adaptés à la MAN pour détecter une séquence de bifurcations multiples ainsi que les modes d’instabilité associés sur leur chemin d’évolution de post-flambement. La MAN donne un accès interactif aux branches d’équilibre semi-analytique, qui offre un avantage considérable en termes de la fiabilité par rapport aux algorithmes itératifs classiques. En outre, une stratégie originale de couplage non-local est développée pour coupler les modèles classiques et les modèles multi-échelles concurremment, où les forces de chaque modèle sont pleinement exploitées, et leurs lacunes surmontées. Une discussion sur la transition entre les différentes échelles est fournie d’une manière générale, qui peut également être considéré comme un guide pour les techniques de couplage impliquant d’autres modèles réduits. A la fin, un cadre général de modélisation macroscopique est développé et deux modèles spécifiques de type Fourier sont dérivés de modèles classiques bien établis, qui permettent de prédire la formation des modes d’instabilités avec beaucoup moins d’éléments et donc de réduire le coût de calcul de manière significative / Surface wrinkles of stiff thin layers attached on soft materials have been widely observed in nature and these phenomena have raised considerable interests over the last decade. The post-buckling evolution of surface morphological instability often involves strong effects of geometrical nonlinearity, large rotation, large displacement, large deformation, loading path dependence and multiple symmetry-breakings. Due to its notorious difficulty, most nonlinear buckling analyses have resorted to numerical approaches since only a limited number of exact analytical solutions can be obtained. This thesis proposes a whole framework to study the film/substrate buckling problem in a numerical way: from 2D to 3D modeling, from classical to multi-scale perspective. The main aim is to apply advanced numerical methods for multiple-bifurcation analyses to various film/substrate models, especially focusing on post-buckling evolution and surface mode transition. The models incorporate Asymptotic Numerical Method (ANM) as a robust path-following technique and bifurcation indicators well adapted to the ANM to detect a sequence of multiple bifurcations and the associated instability modes on their post-buckling evolution path. The ANM gives interactive access to semi-analytical equilibrium branches, which offers considerable advantage of reliability compared with classical iterative algorithms. Besides, an original nonlocal coupling strategy is developed to bridge classical models and multi-scale models concurrently, where the strengths of each model are fully exploited while their shortcomings are accordingly overcome. Discussion on the transition between different scales is provided in a general way, which can also be seen as a guide for coupling techniques involving other reduced-order models. Lastly, a general macroscopic modeling framework is developed and two specific Fourier-related models are derived from the well-established classical models, which can predict the pattern formation with much fewer elements so as to significantly reduce the computational cost

Les voies sèches de valorisation du bois énergie sont de plus en plus exigeantes vis-à-vis de la qualité de la matière première et notamment quant à sa teneur en eau. Ainsi, une étape intermédiaire de préconditionnement, sur site industriel, tend à se développer. La mise en œuvre de séchoirs en continu est alors séduisante de par son faible coût et sa parfaite intégration dans la ligne de production. Cependant, elle n'est pas sans écueil : la variabilité de la biomasse et l'hétérogénéité des conditions climatiques au sein du séchoir rendent fastidieux le dimensionnement du séchoir et de la source de chaleur associée d'une part et la maximisation du flux matière d'autre part.Notre étude a pour ambition de développer un outil informatique d'aide à la conception et à l'optimisation de séchoirs en continu traversés par un lit condensé de particules de bois. Pour y parvenir, nous associons une approche expérimentale et une approche numérique du procédé. La démarche scientifique s'articule autour de deux échelles représentatives :- A l'échelle de la particule de bois, le modèle de Van Meel, reposant sur le concept de courbe caractéristique de séchage, est étendu : le couplage entre transferts de chaleur et de masse est rendu explicite. Nous disposons dès lors d'un modèle réactif et prédictif en conditions climatiques variables. Parallèlement, un premier dispositif expérimental original est conçu et construit pour l'étude du séchage de particules de bois isolées. Les données recueillies sont alors analysées en vue d'alimenter en paramètres et de valider le modèle semi-analytique ;- A l'échelle du séchoir, une modélisation double-échelle est adoptée pour tenir compte de l'évolution des conditions climatiques. Un soin particulier est apporté à la modélisation, se voulant être la plus proche possible des installations existantes : à ce titre, les transferts au niveau de la paroi du séchoir sont pris en compte, ainsi que la variabilité des particules de bois. Parallèlement, un second dispositif expérimental original est conçu et construit pour l'étude du séchage de particules disposées en lit. Les données recueillies sont alors confrontées aux simulations du modèle double-échelle.La validation du modèle s'étant révélée probante tant à l'échelle de la particule qu'à celle du séchoir, une utilisation du code pour l'aide à la conception et à l'optimisation d'installations industrielles a pu être envisagée. Ainsi, ce travail s'achève-t-il par deux études de cas à travers lesquelles nous explorons le potentiel du modèle pour maximiser le flux matière tout en garantissant la qualité du produit et l'efficacité énergétique du séchoir. / The dry conversion routes of wood to energy require more and more improvement in the quality of raw material, particularly regarding to its moisture content. That is why a preliminary step of preconditioning tends to be developed on the industrial sites. The use of continuous dryers is then tempting because of their low cost and their perfect integration in the production line. However it is not without pitfall: the climatic conditions heterogeneity inside the dryer, on the one hand, and the biomass variability, on the other hand, make the utilization of these industrial plants tedious.Our study aims to develop a computing tool to help engineers optimizing and to designing continuous dryers passed through by a condensed bed of wood particles. For this result, an experimental approach of the process is associated to a numerical one. The scientific method revolves around two representative scales:• At the wood particle scale, the Van Meel model, based on the concept of characteristic drying curve, is enhanced: the coupling between heat and mass transfer is made explicit. So we obtain a reactive model which is predictive when the climatic conditions are variable. Meanwhile, a first original experimental device is designed and built in order to study the drying of single wood particles. Thanks to the collected data, the semi-analytical model is provided with parameters and is validated;• At the dryer scale, a dual-scale approach is chosen in order to take into account the evolution of the climatic conditions. A particular attention is paid to the modeling in order to make it as similar as possible to the existing installations: as such, the transfers at the wall surface and the variability of wood particles are included. Meanwhile, a second original experimental device is designed and built in order to study the drying of a wood particles stack. The data collected are then confronted with simulations of the double-scale model.The validation of the model has proven to be convincing at the particle scale and at the dryer one. Therefore a utilization of the code to help in the design and optimization of industrial plants can be considered. Thus, the work ends by two case studies through which we explore the potential of the model to maximize material flow while ensuring product quality and efficiency of the dryer.

Ce travail de thèse porte sur la modélisation mathématique et la simulation numérique de la migration par gravité et capillarité du CO₂ supercritique injecté dans un site de séquestration géologique hétérogène. Les simulations sont réalisées à l'aide du code DuMux. Particulièrement, on s'intéresse à la mise à l'échelle, de l'échelle de la cellule à l'échelle du réservoir, d'un modèle d'écoulement diphasique CO₂ -saumure, au sein d'un milieu stratifié périodique constitué d'un réseau de barrières peu perméables horizontales, continues ou discontinues. La mise à l'échelle est effectuée par la méthode asymptotique à double échelle. Dans un premier temps, on considère le cas d'une colonne verticale parfaitement stratifiée. Un modèle homogénéisé est développé puis validé par simulation numérique pour différentes valeurs du nombre capillaire et du flux incident de CO₂ . La méthode d'homogénéisation est appliquée au cas d'un écoulement dans un milieu bidimensionnel constitué de strates discontinues. Par l'effet de gravité, le CO₂ s'accumule sous les strates peu perméables, ce qui conduit à un problème mathématique local non standard. Cette stratification est modélisée à l'aide de l'approche des courants de gravité. L'approche est étendue au cas des strates semi-perméables et en prenant en compte la capillarité. Le modèle mis à l'échelle est comparé à des simulations numériques effectuées pour différents types de strates, avec ou sans pression capillaire, et sa limite de validité est discutée pour chacun de ces cas. La dernière partie de la thèse est dédiée à l'étude des performances du code DuMux pour simuler par calcul parallèle l'injection et la migration de CO₂ dans des milieux hétérogènes tridimensionnels (milieu périodique stratifié, milieu fluviatile et milieu réservoir SPE10). / This work deals with the mathematical modeling and the numerical simulation of the migration under gravity and capillarity effects of the supercritical CO₂ injected into a geological heterogeneous sequestration site. The simulations are performed with the code DuMux. Particularly, we consider the upscaling, from the cell scale to the reservoir scale, of a two-phase (CO₂ -brine) flow model within a periodic stratified medium made up of horizontal low permeability barriers, continuous or discontinuous. The upscaling is done by the two-scale asymptotic method. First, we consider perfectly layered media. An homogenized model is developed and validated by numerical simulation for different values of capillary number and the incident flux of CO₂ . The homogenization method is then applied to the case of a two-dimensional medium made up of discontinuous layers. Due to the gravity effect, the CO₂ accumulates under the low permeability layers, which leads to a non-standard local mathematical problem. This stratification is modeled using the gravity current approach. This approach is then extended to the case of semi-permeable stratas taking into account the capillarity. The upscaled model is compared with numerical simulations for different types of layers, with or without capillary pressure, and its limit of validity is discussed in each of these cases. The final part of this thesis is devoted to the study of the parallel computing performances of the code DuMux to simulate the injection and migration of CO₂ in three-dimensional heterogeneous media (layered periodic media, fluvial media and reservoir model SPE 10).