No description available.

L'oxydation sélective conduit à un appauvrissement du substrat qui contribue aux phénomènes de corrosion sous contrainte. Afin de prédire cet appauvrissement en élément oxydable, un modèle basé sur la diffusion dans le substrat et l'équilibre des flux au niveau de l'interface alliage/oxyde a été développé. La cinétique d'oxydation du matériau est utilisée comme donnée d'entrée du modèle et permet de piloter le déplacement de l'interface alliage/oxyde.Le modèle d'appauvrissement développé a été validé dans la cas de l'oxydation de l'alliage binaire Pt-Ni à 850 °C, en reproduisant la solution analytique du modèle de Wagner. L'utilisation de la méthode des différences finies a permis de lever les hypothèses d'une concentration à l'interface alliage/oxyde et d'un coefficient d'interdiffusion constants.L'appauvrissement initial en élément oxydable d'un alliage binaire (passivation, pré-oxydation) a peu d'influence sur l'évolution de l'appauvrissement de l'alliage lors de son oxydation à haute température. En revanche, lorsque le coefficient d'interdiffusion de l'alliage est plus grand en subsurface qu'en volume de cet alliage (écrouissage), la concentration en élément oxydable au niveau de l'interface alliage/oxyde décroit au cours de l'oxydation. Les profils d'appauvrissement présentent alors un point d'inflexion dû à une diffusion plus rapide au voisinage de l'interface.Enfin, le modèle a été appliqué à l'oxydation en milieu aqueux à une température de 290 °C pour étudier la déchromisation de l'alliage 690 en milieu primaire des REP. Ceci a permis de conforter la cinétique d'oxydation logarithmique directe proposée dans la littérature pour cet alliage. / High temperature selective oxidation of alloys usually results in substrate depletion. This phenomenon can be at the origin of stress corrosion cracking (SCC) initiation. To predict the depletion of the oxidizable element, a model based on diffusion in the substrate and fluxes balance at the oxide/alloy interface has been developed. The oxidation kinetics is used as input data of the model in order to drive the oxide/alloy interface displacement when the oxidizable element is "consumed" from the alloy.The depletion model presented in this work has been validated in the case of the oxidation at 850 °C of the Pt-Ni binary alloy by reproducing the Wagner's analytical solution. The two assumptions of a fixed concentration at the oxide/alloy interface, and of a diffusion coefficient homogeneous in the substrate and constant during oxidation, have been relaxed using finite difference method.An initial depletion of a binary alloy in oxidizable element (passivation, pre-oxidation) has a low impact on the substrate depletion evolution during high temperature oxidation. However, when the interdiffusion coefficient in the alloy is greater in subsurface than in the bulk (work-hardening), the oxidizable element concentration at the oxide/alloy interface decreases during the oxidation. Depletion profiles show an inflexion point due to faster diffusion close to the interface.Finally, the model has been applied to oxidation in aquaous media at a temperature of 290 °C to study Cr-depletion of alloy 690 in PWR primary water. The direct logarithmic oxidation kinetics proposed in the literature has been validated

Au cours des dernières décennies, les métamatériaux ont montré qu'ils pouvaient avoir des propriétés étonnantes, permettant d'imaginer et de réaliser des dispositifs tels que des lentilles à indice négatif, des dispositifs de camouflage ou des lentilles parfaitement plates. Aussi, font-ils aujourd'hui l'objet d'une intense recherche. Cette thèse, inscrite dans le projet ANR METAPHOTONIQUE, avait pour objectif principal d'étudier la possibilité d'utiliser des métamatériaux pour l'optique guidée à la longueur d'onde télécom. Ceci a requis une meilleure compréhension du comportement des métamatériaux, ce qui nous a orientés vers l'étude de modèles effectifs. La récupération des paramètres optiques effectifs a nécessité de pouvoir simuler la réponse des métamatériaux sous une illumination en incidence oblique, ce qui a requis l'implémentation de telles méthodes en FDTD. Aux cours de nos travaux, nous avons aussi mis en évidence un phénomène de rotation de la polarisation de modes guidés suite à leur interaction avec des métamatériaux asymétriques. Nous avons alors pu concevoir et caractériser un dispositif pour l'optique guidée destiné faire tourner la polarisation de la lumière sur une longueur deux fois plus courte que la longueur d'onde. Enfin, un objectif secondaire a été de déterminer si les pertes dues à la présence de nanostructures métalliques dans nos métamatériaux pouvaient être compensées grâce à l'utilisation de matériaux à gain. Un modèle permettant de décrire de tels matériaux en FDTD a donc été étudié et implémenté / Over the last decades, metamaterials have been shown to exhibit extraordinary properties. These properties could allow the design of new devices such as negative refractive index lenses, cloaking devices and perfectly flat lenses. This PhD is part of the ANR METAPHOTONIQUE project and its main objective was to study the potential of metamaterials for integrated photonics at telecommunication wavelengths. This required to first better understand the behavior of metamaterials, which led us to study effective index models. The retrieval of effective optical parameters required to be able to simulate the optical response of metamaterials under oblique incidence excitation. Due to this, we had to study and implement such methods for FDTD computations. Our work also showed that asymmetrical metamaterials possessed the ability to make the polarization of guided modes rotate. This allowed us to design and characterize a potential device for integrated photonics, which effectively converts the polarization of guided modes over lengths shorter than half the wavelength. A secondary goal of this PhD was to study the potential compensation of losses occuring within the metallic nanostructures of our metamaterials through the use of materials with gain. Thus, a model allowing us to simulate such materials in FDTD has been implemented

Dans ce travail de thèse, nous étudions expérimentalement la dynamique vibrationnelle de nano-objets uniques métalliques ou semi-conducteurs par des techniques pompe/sonde femtosecondes. Nous démontrons d’abord que l’observation de nanofils uniques, suspendus au-dessus de tranchées permet à la fois de s’affranchir de l’effet d’étalement inhomogène des propriétés acoustiques et d’augmenter le confinement acoustique. Grâce à l’enrichissement du paysage vibrationnel ainsi obtenu, nous explorons les propriétés élastiques de nombreux systèmes : métalliques, semi-conducteurs, poreux, alliages, cœur/coquille, etc. Ensuite, nous tirons parti de l’augmentation du confinement acoustique pour observer la propagation d’ondes acoustiques gigahertz guidées le long de nanofils ou de nanopoutres. Nous montrons que la propagation de ces ondes acoustiques dans ces guides d’ondes nanométriques permet d’obtenir des informations indépendantes sur les propriétés élastiques des objets. A contrario, nous mettons en évidence que lorsque le nanofil est en contact avec un substrat, il agit comme une source acoustique monochromatique d’ondes longitudinales qui rayonne dans le substrat. Nous réalisons en transmission l’imagerie spatio-temporelle de ce champ acoustique généré et détectons acoustiquement l’orientation du nanofil sous-jacent grâce à l’anisotropie de forme du champ acoustique. Enfin, nous envisageons une preuve de concept d’un système de microscopie acoustique de résolution spatiale nanométrique en utilisant une pointe de microscopie à force atomique. / Vibrational dynamics of individual nano-objects is studied experimentally using pump and probe time-resolved spectroscopy. First, suspended and individual nano-objects avoid the inhomogeneous broadening of the acoustic properties and increase the acoustics confinement inside the nano-object. Elastic properties of metallic, semiconducting, porous, alloys or core-shell nanowires are thus studied in this advantageous geometry. The increased acoustic confinement in the suspended geometry also lead us to the observation of gigahertz coherent guided acoustic phonons in single copper nanowires and gold nanobeams. We show that the observation of propagating acoustic waves in nanoscale waveguides provide additional elastic informations. Furthermore, it gives the opportunity to unambiguously discriminate which mode is excited and detected using pump and probe time-resolved spectroscopy. On the contrary, nanowires can be used as monochromatic acoustic sources of longitudinal waves when deposited on a substrate. As the acoustic source radiates longitudinal waves inside the substrate, the spatiotemporal imaging of the generated acoustic field is undergone and the nanowire’s orientation is detected in transmission thanks to the acoustic field’s anisotropy. Finally, as another step toward acoustic microscopy with nanoscale spatial resolutions, an atomic force microscopy tip is used as a waveguide and an acoustic transducer with nanometric spatial extension.

Les préoccupations environnementales actuelles visent à réduire les consommations énergétiques. Dans une démarche d’amélioration des bâtiments existants, l'étude du comportement thermique d’une paroi n'est pas aisée du fait de la méconnaissance de ses propriétés thermophysiques réelles. Ces paramètres sont pourtant prépondérants pour la phase d'optimisation économique des opérations de réhabilitation ou pour vérifier ses performances in situ. Il apparaît donc important de pouvoir caractériser les parois de bâtiment en place. Notre travail vise à développer une méthode de caractérisation thermique d’une paroi adaptée aux applications in situ basée sur une approche active. Le principe d'identification consiste à solliciter thermiquement une face d’accès en imposant un flux de chaleur sous forme d’un créneau et à étudier la réponse en température enregistrée par thermographie infrarouge sur l’autre face. A partir de signaux de flux et de températures mesurés aux limites de la paroi, les propriétés thermophysiques de la paroi seront estimées par méthode inverse. Nous nous sommes dans un premier temps intéressés aux parois homogènes. Le schéma d’inversion est construit autour d’un modèle numérique décrivant la réponse de la paroi suivant la méthode des différences finies en 1D. L’identification de la conductivité thermique et de la chaleur volumique de la paroi est réalisée en optimisant le groupement de paramètres qui permet de minimiser l’écart entre la température normalisée mesurée et la température normalisée simulée. Le coefficient d’échange surfacique global est également identifié à partir du même essai. Dans ce travail, la méthode a été appliquée à une paroi homogène en carreaux de plâtre mise en place au laboratoire. Elle a une épaisseur de 6.5 cm. Cette technique a été utilisée pour les parois multicouches de bâtiments. Les résultats issus de cette procédure d’inversion ont été comparés à des valeurs de référence obtenues à partir d’une procédure classique (NF EN 12664-méthode fluxmétrique). Une bonne concordance des résultats est obtenue. Une autre partie représente les essais in situ. / Current environmental concerns are intended to reduce energy consumption. In a process of improving existing buildings, the study of the thermal behavior of a wall is not easy because of the ignorance of its real thermophysical properties. These parameters are yet to dominate the economic optimization phase of the rehabilitation or to check its performance in situ. It therefore appears important to characterize the walls of existing building. Our work aims to develop a method of thermal characterization of a wall suitable for in situ applications based on an active approach. The principle of identification is to apply a heat-face access by imposing a heat flux in the form of a pulse and to study the temperature response recorded by infrared thermography on the other side. From signal flow and temperature measured at the limits of the wall, the thermophysical properties of the wall will be estimated by inverse method. We are at present interested in homogeneous walls. The inversion scheme is built around a digital model describing the response of the wall following the finite difference method in 1D. The identification of the thermal conductivity and heat volume of the wall is achieved by optimizing the group of parameters which minimizes the normalized difference between the temperature measured and the temperature standard simulated. The overall Global exchange coefficient is also identified from the same test. In this work, the method was applied to a homogeneous wall tile plaster introduction to the laboratory. It has a thickness of 6.5 cm. This technique was used for multilayer walls of buildings. The results of this inversion procedure were compared with reference values obtained from a standard procedure (DIN EN 12664-flow meter methods). A good agreement is obtained. Another part is the in situ tests.

Les temps de trajet, amplitudes et angles à la source des ondes sismiques sont utilisés dans de nombreuses applications telles que la migration, la tomographie, l'estimation de la sensibilité de détection et la localisation des microséismes. Dans le contexte de la microsismicité, il est nécessaire de calculer en quasi temps réel ces attributs avec précision. Nous avons développé ici un ensemble d'algorithmes rapides et précis en 3D pour des modèles à fort contraste de vitesse.Nous présentons une nouvelle méthode pour calculer les temps de trajet, les amplitudes et les angles à la source des ondes correspondant aux premières arrivées. Plus précisément, nous résolvons l'équation Eikonal, l'équation de transport et l'équation des angles en nous basant sur une approche par différences finies pour des modèles de vitesse en 3D. Nous proposons une nouvelle méthode hybride qui bénéficie des avantages respectifs de plusieurs approches existantes de résolution de l'équation Eikonal. En particulier, les approches classiques proposent généralement de résoudre directement les équations et font l'approximation localement d'une onde plane. Cette approximation n'est pas bien adaptée au voisinage de la source car la courbure du front d'onde est importante. Des erreurs de temps de trajet sont alors générées près de la position de la source, puis propagées à travers tout le modèle de vitesse. Ceci empêche de calculer correctement les amplitudes et les angles à la source puisqu'ils reposent sur les gradients des temps. Nous surmontons cette difficulté en introduisant les opérateurs sphériques ; plus précisément nous reformulons les temps de trajet, amplitudes et angles à la source par la méthode des perturbations.Nous validons nos nouvelles méthodes pour différents modèles à fort contraste de vitesse en 2D et 3D et montrons notre contribution par rapport aux approches existantes. Nos résultats sont similaires à ceux calculés en utilisant la modélisation de la forme d'onde totale alors qu'ils sont bien moins coûteux en temps de calcul. Ces résultats ouvrent donc de nouvelles perspectives pour de nombreuses applications telles que la migration, l'estimation de la sensibilité de détection et l'inversion des mécanismes au foyer. / Traveltimes, amplitudes and take-off angles of seismic waves are used in many applications such as migration, tomography, detection sensitivity estimation and microseism location. In the microseismicty context it is necessary to compute in near real time accurately these attributes. Here we developed a set of fast and accurate algorithms in 3D for highly contrasted velocity models.We present a new accurate method for computing first arrival traveltimes, amplitudes and take-off angles; more precisely we solve the Eikonal, transport and take-off angle equations based on a finite difference approach for 3D velocity models. We propose a new hybrid method that benefits from the advantages of several existing Eikonal solvers. Common approaches that solve directly these equations assume that we are locally propagating a plane wave. This approximation is not well adapted in the neighborhood of the source since the wavefront curvature is important. Travel times errors are generated near the source position and then propagated through the whole velocity model. This prevents from properly calculating the amplitudes and the take-off angles since they rely on the travel time gradients that are not accurate. We overcome this difficulty by introducing spherical operators. Indeed we reformulate the traveltimes, amplitudes and take-off angles with the perturbation method.We validate our new methods on various highly contrasted velocity models in 2D and 3D and show our contribution compared to other existing approaches. Our results are similar to those computed using full waveform modeling while they are obtained in a much shorter CPU time. These results open thus new perspectives for several applications such as migration, detection sensitivity estimation and focal mechanism inversion.

Ce travail porte sur le problème des partitions minimales, à l'interface entre théorie spectrale et optimisation de forme. Une introduction générale précise le problème et présente des résultats, principalement dûs à B. Helffer, T. Hoffmann-Ostenhof et S. Terracini, qui sont utilisés dans le reste de la thèse.Le premier chapitre est une étude spectrale asymptotique du laplacien de Dirichlet sur une famille de domaines en dimension deux qui tend vers un segment. L'objectif est d'obtenir une localisation des lignes nodales dans la limite des domaines minces. En appliquant les résultats de Helffer, Hoffmann-Ostenhof et Terracini, on montre ainsi que les domaines nodaux des premières fonctions propres forment des partitions minimales.Le deuxième chapitre étudie les valeurs propres de certains opérateurs de Schrödinger sur un domaine plan avec condition au bord de Dirichlet. On considère des opérateurs qui ont un potentiel électrique nul et un potentiel magnétique d'un type particulier, dit d'Aharonov-Bohm, avec des singularités en un nombre fini de points appelés pôles. On démontre que les valeurs propres dépendent continuement des pôles. Dans le cas de pôles distincts et éloignés du bord, on prouve que cette dépendance est analytique lorsque la valeur propre est simple. On exprime de plus une condition suffisante pour que la fonction qui aux pôles associe une valeur propre présente un point critique. On utilise alors la caractérisation magnétique des partitions minimales pour montrer que l'énergie minimale est une valeur critique d'une de ces fonctions.Le troisième chapitre est un article écrit en collaboration avec Virginie Bonnaillie-Noël. Il porte sur une famille d'exemples, les secteurs angulaires de rayon unité et d'ouverture variable, dont on tente de déterminer les partitions minimales. On applique pour cela les théorèmes généraux rappelés dans l'introduction afin de déterminer les partitions nodales qui sont minimales. On s'intéresse plus particulièrement aux partitions minimales en trois domaines. En appliquant les idées du deuxième chapitre, on montre que pour certaines valeur de l'angle, il n'existe aucune partition minimale qui soit symétrique par rapport à la bissectrice du domaine. D'un point de vue quantitatif, on obtient des encadrements précis de l'énergie minimale.Le quatrième chapitre consiste en l'étude des partitions minimales de tores plats dont on fait varier le rapport entre longueur et largeur. On utilise une méthode numérique très différente de celle du troisième chapitre, basée sur un article de B. Bourdin, D. Bucur et É. Oudet. Elle consiste en une relaxation suivie d'une optimisation par un algorithme de gradient projeté. On peut ainsi tester des résultats théoriques antérieurs. Les résultats présentés suggèrent de plus la construction explicite de familles de partitions (en liaison avec des pavages du tore) qui donnent une nouvelle majoration de l'énergie minimale.Un dernier chapitre de perspectives présente plusieurs applications possibles des méthodes décrites dans la thèse. / This work is concerned with the problem of minimal partitions, at the interface between spectral theory and shape optimization. A general introduction gives a precise statement of the problem and recall results, mainly due to B. Helffer, T. Hoffmann-Ostenhof and S.Terracini, that are used in the rest of the thesis.The first chapter is an asymptotic spectral study of the Dirichlet Laplacian on a familly of two-dimensional domains converging to a line segment. The aim is to localize the nodal lines when the domains become very thin. With the help of the results of Helffer, Hoffmann-Ostenhof, and Terracini, we then show that the nodal domains of the first eigenfunctions give minimal partitions.The second chapter studies the eigenvalues of some Schrödinger operators on a domain with Dirichlet boundary conditions. We consider operators that have no electric potential and a so-called Aharonov-Bohm magnetic potential, which has singularities at a finite number of points called poles. We prove that the eigenvalues are continuous functions of the poles. When the poles are distinct and far from the boundary, we prove that this function is analytic, assuming the eigenvalue is simple. We also give a sufficient condition for the function to have a critical point. Using the magnetic characterization of minimal partitions, we show that the minimal enery is a critical value for one of these functions.The third chapter in an article written in collaboration with Virginie Bonnaillie-Noël. It studies minimal partitions for sectors of unit radius with a variable angular opening. We apply the general results presented in the introduction, together with numerical computations, to determine nodal partitions that are minimal. We focus on partitions into three domains. Using ideas from the second chapter, we show that, for some values of the angle, there is no minimal partition that is symmetric with respect to the bisector. Form a quantitative point of view, we obtain precise bounds on the minimal energy.The fourth chapter studies the minimal partitions of flat tori in function of the ratio between width and length. We use a numerical method that is quite different from chapter three, and is based on an article by B. Bourdin, D. Bucur, and É. Oudet. It consists in a relaxation of the problem, followed by optimization with the help of a projected gradient algorithm. The results shown here additionally suggest explicit families of partitions, which consist in tilings of tori by polygons, that give upper bounds on the minimal energy. In the last chapter we consider several possible applications of the methods described in the thesis.

La communauté scientifique, en relation avec les gestionnaires forestiers, travaille depuis trente ans à l'élaboration d'outils d'aide à la décision. Cependant, aucune approche de modélisation ne permet une évaluation simultanée de l'impact sur la forêt des changements globaux et de la gestion des services écosystémiques. L'élaboration d'une nouvelle génération de modèles dédiés au couplage sol-plante est indispensable pour aider les gestionnaires forestiers à adapter leurs pratiques sylvicoles face aux changements globaux. Les modèles phénoménologiques issus de la dendrométrie proposent des estimations de croissance et de production tenant compte des pratiques sylvicoles. Leur utilisation est limitée par leur dépendance à un indice de fertilité combinant de manière indifférenciée l'influence du climat et des cycles biogéochimiques. La remise en jeu de la notion d'indice de fertilité nécessite de tirer avantage des concepts issus des modèles à base écophysiologiques et biogéochimiques. La modélisation de la dynamique mensuelle du renouvellement foliaire par le modèle Stand Leaf Canopy Dynamics a permis d'estimer: i) la dynamique de l'indice foliaire, ii) la production de litière foliaire. Ce modèle probabiliste tient compte de l'influence du climat, de la disponibilité en eau de l'écosystème et repose sur des hypothèses écologiques fortes sur les mécanismes impliqués dans la longévité des feuilles. Ce modèle pose les bases de l'intégration du cycle de l'eau aux modèles dendrométriques par couplage avec un modèle écophysiologique. Le modèle de décomposition de la matière organique le long d'un profil de sol vertical proposé par Bosatta et Ågren (1996) a été analysé mathématiquement afin d'étudier le devenir des chutes de litières dans le sol forestier. Nous avons proposé: i) une preuve d'existence et l'unicité de solutions à l'équation de transport intégro-différentielle du modèle, ii) des schémas de différences finies implicites-explicites convergents estimant les solutions du modèle et iii) une discussion autour des hypothèses sous-jacentes à une troncature du modèle proposée par Bosatta et Ågren. Ces deux étapes sont des leviers indispensables à l'étude des cycles biogéochimiques dans une approche à base dendrométrique. L'influence du climat et de la phénologie sur le cycle de l'eau, la décomposition des litières foliaires et la chimie du sol constituera une prochaine étape de modélisation. / The scientific community, in collaboration with forestry managers, developed decision support tools since thirty years. However, current modelling approaches do not allow simultaneous estimations of global changes impact management policies on forests. The design of a new generation of models dedicated to soil-plant coupling is necessary to help forest managers to adapt forestry practices to face global changes. Phenomenological models arising from dendrometry estimate forest growth \& yield and take into account of silvicultural practices. Their use is limited by their dependence to a fertility index, which combines implicitly the influence of climate and biogeochemical cycles. A re-examination of fertility index concept is necessary to take advantage of concepts developed in ecophysiological and biogeochemical models. Modelling the monthly dynamics of foliar renewal with the Stand Leaf Canopy Dynamics model permited: i) to estimate leaf area index dynamics, ii) to model foliar litter production. This probabilistic model takes into account climate influence, ecosystem water availability and lies on strong ecological hypothesis on mechanisms involved in leaf longevity. This model, associated with an ecophysiological model, permits to integrate water cycle to dendrometric models. The soil organic matter decomposition along a vertical soil profile model introduced by Bosatta et Ågren (1996) was analyzed to study the future of litterfall in forest soil. We proposed: i) a proof of existence and uniqueness of solutions to the integro-differential transport equation of the model, ii) convergent implicit-explicit finite differences schemes estimating model solutions and iii) a discussion on hypothesis underlying a model truncation introduced by Bosatta et Ågren. These two modelling steps are essential to consider biogeochemical cycle with a dendrometric approach. Climate and phenology influence on water cycle, litterfall decomposition and soil chemistry are the elements of a future modelling phase.

Dans cette thèse on a développé des algorithmes semi-implicites procédures partagées et monolithiques pour l'interaction entre un fluide gouverné par le modèle de Navier Stokes et une structure. Dans le premier chapitre, on présente un algorithme semi-implicite procédures partagées pour l'interaction entre un fluide et une structure gouvernée soit par les équations d'élasticité linéaire ou soit par le modèle de Saint-Venant Kirchhoff non linéaire. Dans le second chapitre, on propose un algorithme semi-implicite procédures partagées pour l'interaction entre un fluide et une structure de modèle linéaire et on montre un résultat de stabilité inconditionnelle en temps de l'algorithme. Un problème d'optimisation est résolu dans les deux algorithmes précédents, afin de satisfaire les conditions de continuité des vitesses et d'égalité des contraintes à l'interface. Durant les itérations de BFGS pour résoudre le problème d'optimisation, le maillage fluide reste fixe et la matrice fluide n'est factorisée qu'une seule fois, ce qui réduit l'effort de calcul. Dans le troisième chapitre, un algorithme semi-implicite monolithique pour l'interaction entre un fluide et une structure de modèle linéaire est proposé. L'algorithme utilise un maillage global pour le domaine fluide structure. La condition de continuité des vitesses à l'interface est automatiquement satisfaite et celle de l'égalité des contraintes n'apparaît pas explicitement dans la formulation faible. A chaque pas de temps on résout un système monolithique d'inconnues vitesse et pression définies sur le domaine global. Le temps CPU est réduit quand l'approche monolithique est utilisée à la place des procédures partagées. / Our aim was to develop some partitioned procedures and monolithic semi-implicit algorithms for solving the interaction between a fluid governed by Navier Stokes equations and a structure. In the first chapter, we propose a partitioned procedures semi-implicit algorithm for solving fluid-structure interaction problems, with a structure governed either by linear elasticity equations or by the non-linear Saint-Venant Kirchhoff model. In the second chapter, we present a partitioned procedures semi-implicit algorithm for solving fluid-structure interaction problem with a linear model for the structure and we prove an unconditional stability result of the algorithm. In the above algorithms, an optimization problem must be solved in order to get the continuity of the velocity as well as the continuity of the stress at the interface. During the iterations of BFGS for solving the optimization problem, the fluid mesh does not move and the fluid matrix is only factorized once, which reduces the computational effort. In the fast chapter, we present a monolithic semi-implicit algorithm for solving fluid-structure interaction problem with linear model for the structure. The algorithm uses one global mesh for the fluid-structure domain. The continuity of velocity at the interface is automatically satisfied and the continuity of stress does not appear explicitly in the global weak form due to the action and reaction principle. At each time step, we have to solve a monolithic system of unknowns velocity and pressure defined on the global fluid-structure domain. When the monolithic approach is used the CPU time is reduced compared to a particular partitioned procedures strategy.

Dans la présente étude, la conductivité thermique et le module d'élasticité de revêtementsd’YPSZ élaborés par projection plasma ont été prédits par modélisations numériques 2D et3D de type différences finies et éléments finis.L'influence de la résolution d'image, de la taille et de la valeur du seuil sur les propriétésprédites du revêtement a été étudiée. En outre, les effets de la méthode numérique et du typede condition aux limites ont été étudiés. En particulier, la quantification de l'effet Knudsen(effet de raréfaction) sur le transfert de chaleur à travers une structure poreuse a été réaliséepar modélisation numérique en combinaison avec l'analyse d'image. Les conductivitéseffectives obtenues par modélisation 3D s'avèrent plus élevées que celles obtenues en 2D, etaussi en meilleur accord avec les résultats mesurés. Une corrélation 2D/3D a été trouvéepour la modélisation de la conductivité thermique : cette corrélation permet de prédire lesvaleurs 3D à partir des valeurs calculées en 2D. / In the present study, the thermal conductivity and elastic modulus of thermal spray YPSZcoatings were predicted by 2D and 3D finite differences and finite elements numericalmodeling based on cross-sectional images.The influence of the image resolution, size and threshold on the predicted properties of thecoating was studied. Moreover, the effects of the numerical method and of the boundarycondition were investigated. In particular, the quantification of the Knudsen effect(rarefaction effect) on the heat transfer through a porous structure was realized by numericalmodeling in combination with image analysis. The predicted thermal conductivities obtainedby 3D modeling were found to be higher than those obtained by 2D modeling, and in betteragreement with the measured results. A 2D/3D correlation was sucessfully found for themodeling of thermal conductivity: this correlation allows predicting 3D computed valuesfrom 2D ones.