Dans de nombreux problèmes physiques, on s'intéresse à l'interaction de phénomènes ayant des longueurs caractéristiques très différentes. On parle alors de phénomènes multi-échelles. Les méthodes numériques classiques, comme les éléments finis ou les différences finies, nécessitent alors un pas de maillage de l'ordre de la plus petite longueur caractéristique. Ceci a pour effet de faire exploser le nombre de degrés de liberté et les coûts de calcul. Afin de pallier à cette difficulté, on trouve dans la littérature différentes méthodes qui consistent soit à 1. développer des méthodes purement numériques de raffinement local. Ainsi, on limitera le nombre de degrés de liberté. 2. dériver à l'aide de l'analyse asymptotique des modèles approchés dont la solution peut être approchée numériquement sans raffinement local. De nombreux auteurs se sont intéressés à ces problématiques. 3. combiner une analyse asymptotique avec une méthodes de type éléments finis en augmentant l'espace de Galerkin par des fonctions reproduisant les propriétés locales de la solution du modèle exacte. Cette problématique est particulièrement présente en propagation d'ondes. En effet, des détails géométriques de petites tailles (fil mince, fente mince, petit trou, couche mince) devant la longueur d'ondes peuvent avoir une influence significative. Lors de ces dernières années, j'ai étudié ces phénomènes dans le cadre des problèmes de propagation d'ondes linéaires scalaires en régime fréquentiel. Dans cet exposé, je vous présenterai quelques résultats que j'ai obtenu à l'aide de la technique des développements asymptotiques raccordés.

Cette thèse est consacrée à l'étude de l'influence d'une fine couche hétérogène sur le comportement élastique linéaire d'une structure tridimensionnelle.Deux types d'hétérogénéités sont pris en compte : des cavités et des inclusions élastiques. Une étude complémentaire, dans le cas d'inclusions de grande rigidité, a été réalisée en considérant un problème de conduction thermique.Une analyse formelle par la méthode des développements asymptotiques raccordés conduit à un problème d'interface qui caractérise le comportement macroscopique de la structure. Le comportement microscopique de la couche est lui déterminé sur une cellule de base. Le modèle asymptotique obtenu est ensuite implémenté dans un code éléments finis. Une étude numérique permet de valider les résultats de l'analyse asymptotique. / This thesis is devoted to the study of the influence of a thin heterogeneous layeron the linear elastic behavior of a three-dimensional structure. Two types of heterogeneties are considered : cavities and elastic inclusions. For inclusions of high rigidty a further study was performed in the case of a heat conduction problem.A formal analysis using the matched asymptotic expansions method leads to an interface problem which characterizes the macroscopic behavior of the structure. The microscopic behavior of the layer is determined in a basic cell.The asymptotic model obtained is then implemented in a finite element software.A numerical study is used to validate the results of the asymptotic analysis.

Dans cette thèse, nous nous intéressons à la résolution des équations de Maxwell dans une structure périodique constituée d'un anneau mince de matériau diélectrique de rayon moyen r à l'intérieur duquel s'enroulent deux nappes de fils hélicoïdaux. L'épaisseur de l'anneau et la distance entre deux fils consécutifs sont du même ordre de grandeur d et nous supposons que d est bien inférieur à la longueur d'onde de l'onde incidente ainsi qu'au rayon moyen r. La présence des deux échelles rend les simulations numériques directes difficiles (il est alors nécessaire de mailler la structure à l'échelle du fil). C'est pourquoi nous construisons des modèles approchés dans lesquels l'anneau périodique est remplacé par une condition de transmission posée sur l'interface médiane S. La résolution du modèle approché par une méthode d'éléments finis est bien moins coûteuse que celle du problème exact car il n'y a plus besoin de mailler les fils. La construction des modèles approchés repose sur un développement asymptotique de la solution en fonction du petit paramètre d. Nous utilisons une méthode couplant les techniques d'homogénéisation et des développements asymptotiques raccordés. Les conditions de transmission approchées se construisent alors à l'aide du développement asymptotique tronqué. Nous accordons une attention particulière à la stabilisation des modèles approchés ainsi qu'à leur justification théorique. Enfin, nous validons nos modèles par des simulations numériques.

Dans cette thèse, nous souhaitons décrire la dynamique du front d'avancement d'un film mince s'écoulant sur un plan incliné non rugueux. Nous nous intéressons surtout au problème de point triple situé à l'interface entre la paroi solide, le fluide en mouvement et l'air, par exemple lors de l'écoulement d'une goutte sur une surface inclinée. Dans une première partie, nous expliquons pourquoi on peut se ramener aux équations de Stokes et pourquoi le problème résultant est mal posé. Pour y remédier, la condition de non-glissement à la paroi est remplacée par une condition de glissement lorsqu'on est proche du front. Ainsi on réussit à trouver une solution dans H1. Puis nous développons la dynamique de l'écoulement à l'amont du front : un film mince. Cet écoulement peut se modéliser sous la forme d'équations de type Saint-Venant sur la hauteur et le débit. Nous justifions cette construction à partir des équations de Navier-Stokes en utilisant un développement asymptotique en fonction du paramètre onde longue. Dans la zone du front nous résolvons le système de Stokes stationnaire avec glissement au fond par un développement asymptotique en fonction du nombre capillaire. Le front est divisé en une zone interne près du front et une zone externe loin du front, puis les solutions de chaque zone sont soit raccordées directement (angles dynamique et statique égaux), soit raccordées au moyen d'une zone intermédiaire (angles dynamique et statique différents). Cela nous conduit à deux familles de modèles. En réunissant les modèles type Saint-Venant et les différents modèles de front, nous obtenons un modèle de Saint-Venant tenant compte de la dynamique du front. À partir de ce modèle à deux équations nous pouvons écrire un modèle plus simple à une équation sur la hauteur. Ce modèle permet d'étendre les modèles existants avec adhérence à des modèles avec glissement. On peut alors réaliser des simulations numériques combinant un front d'avancement et un film mince / In the present work, we describe the dynamics of a moving contact line for thin films flowing down an inclined plane. Our focus is the problem of triple point located at the interface between the solid wall, the moving fluid and air, for example the spreading of a drop on a plane dry wall (horizontal or inclined) due to gravity and capillarity. In the first part, we explain how we can reduce to the Stokes equations and why the resulting problem is ill-posed. This singularity is removed by permitting the fluid to slip along the wall close to the contact line. Thus we manage to find a solution in H1 constructed by asymptotic expansions. Then we focus on the upstream dynamic of the flow, which is set to a thin film flow. We develop the classical system of Shallow-water equations (Saint-Venant equations) from the full Navier-Stokes system using the classical long-wavelength expansion. We obtain a set coupled equations for the flow depth and the flow-rate. In the neightboorhood of the contact line, we develop an asymptotic expansion of the steady Stokes system with slip at bottom in function of the capillary number. The solution in the vicinity of the contact line is developped in the inner region and the outer region. Then, a direct matching can be done (assuming dynamic and static angles are equals) or using an intermediate region (with different angles). This leads to two different families of models. Bringing together the upstream Shallow-water equations and the contact line models, we write a new Shallow-water model taking into account the dynamic of the moving contact line. Then, we deduce a simplier one-equation model for the film thickness. This model extends existing models with no slip at bottom to models with slip. Direct numerical simulations of the last models are performed, combining a moving contact line and a thin liquid film

Cette thèse est consacrée à l'étude de l'influence d'une fine couche hétérogène sur le comportement élastique linéaire d'une structure tridimensionnelle. Deux types d'hétérogénéités sont pris en compte : des cavités et des inclusions élastiques. Une étude complémentaire, dans le cas d'inclusions de grande rigidité, a été réalisée en considérant un problème de conduction thermique. Une analyse formelle par la méthode des développements asymptotiques raccordés conduit à un problème d'interface qui caractérise le comportement macroscopique de la structure. Le comportement microscopique de la couche est lui déterminé sur une cellule de base. Le modèle asymptotique obtenu est ensuite implémenté dans un code éléments finis. Une étude numérique permet de valider les résultats de l'analyse asymptotique.

Nous nous intéressons à des problèmes de transmission entre diélectriques et métamatériaux, milieux présentant des propriétés électromagnétiques inhabituelles comme des caractéristiques effectives négatives à certaines fréquences. Par exemple, ces milieux peuvent être construits comme des assemblages périodiques de microstructures résonantes et dans ce cas la théorie de l'homogénéisation permet de justifier mathématiquement ces propriétés effectives. En régime harmonique et dans des géométries à variables séparables, des calculs analytiques peuvent être menés. Ils révèlent dans des cas dits critiques des difficultés mathématiques: les solutions n'ont pas la régularité standard, voire le problème peut être mal posé.La première partie étudie ces problèmes de transmission en régime temporel pour lequel les métamatériaux sont modélisés par des modèles dispersifs (modèle de Drude ou de Lorentz). Les difficultés résident dans le choix d'un schéma de discrétisation mais surtout dans la construction de conditions absorbantes. La méthode retenue ici est celle des Perfectly Matched Layers (PMLs). Comme les PMLs classiques sont instables pour ces modèles du fait de la présence d'ondes inverses, nous proposons une nouvelle classe de PMLs pour lesquelles nous menons une analyse de stabilité. Cette dernière permet de construire des PMLs stables. Elles sont ensuite utilisées pour simuler le comportement en temps long d'un problème de transmission; nous illustrons alors le fait que le principe d'amplitude limite peut être mis en défaut en raison de résonances d'interface.La deuxième partie vise à pallier ces phénomènes d'interface en régime harmonique en revenant sur le processus d'homogénéisation classique, pour un milieu dissipatif. Pour des problèmes de transmission, il est connu que les modèles issus de cette méthode perdent en précision du fait de la présence de couches limites à l'interface. Nous proposons un modèle enrichi au niveau de l'interface. En combinant la méthode d'homogénéisation double-échelle et celle des développements asymptotiques raccordés, nous construisons des conditions de transmission non standards faisant intervenir des opérateurs différentiels le long de l'interface. Le calcul de ces conditions nécessite la résolution de problèmes de cellule et de problèmes non standards posés dans des bandes périodiques infinies. Une analyse d'erreur confirme l'amélioration de la précision du modèle. Des simulations numériques illustrent l'efficacité de ces nouvelles conditions. Enfin, cette démarche est reproduite formellement dans le cas des matériaux à fort contraste se comportant comme des métamatériaux. Nous montrons alors que ces nouvelles conditions permettent de régulariser le problème de transmission dans les cas critiques. / We are interested in transmission problems between dielectrics and metamaterials, that is to say media with unusual electromagnetic properties such as negative constants at some frequencies. These media are often made of periodic assemblies of resonant micro-structures and in this case the homogenization theory can justify mathematically these effective properties. A preliminary part deals with these problems in the harmonic domain and in geometry with separation of variables.Analytical computations are done and reveal in the so-called critical cases some mathematical diffculties: the solutions do not have the standard regularity and the problem can even be ill-posed.The first part examines these transmission problems in the time domain for which metamaterials are modelled by dispersive models (Drude model or Lorentz model for instance). The diffculties reside in the choice of a discretization scheme but especially in the construction of absorbing conditions. The method used here is the use of Perfectly Matched Layers (PMLs). Since classical PMLs are unstable for these models due to the presence of backward waves, we propose a new class of PMLs for which we conduct a stability analysis. The latter allows us to build stable PMLs. They are then used to simulate the long-time behaviour of a transmission problem; we illustrate the fact that the limit amplitude principle can be faulted because of interface resonances.The second part aims to overcome these phenomena by coming back to the classical homogenization in the harmonic domain, for dissipative media. For transmission problems, it is known that models resulting from this method lose accuracy due to the presence of boundary layers at the interface. We propose an enriched model at the interface: by combining the method of two-scale homogenization and that of matched asymptotic expansions, we build non-standard transmission conditions involving tangential derivatives along the interface (Laplace-Beltrami operators). This requires to solve cell problems and non-standardproblems in infinite periodic bands. An error analysis confirms the improvement of the accuracy of the model and numerical simulations show the effectiveness of these new conditions. Finally, this approach is formally reproduced in the case of high contrast materials which behave like metamaterials. We show that these new conditions regularise the transmission problem in the critical cases.

Cette thèse concerne la modélisation mathématique et l'approximation numérique de modèles homogénéisés de métamatériaux. Dans la première partie on étudie des problèmes de propagation d'ondes en présence de métamatériaux homogénéisés tels que les équations de Maxwell, le système de l'acoustique ou de l'élasticité linéaire. Nous établissons des résultats d'existence et d'unicité pour ces systèmes sous des hypothèses phénoménologiques sur le métamatériau en accord avec certains modèles de la littérature. Nous abordons ensuite leurs approximations numériques. Nous présentons des résultats concernant les éléments finis pour l'approximation de l'équation de Helmholtz qui montrent que ce schéma peut ne pas converger en présence de métamatériaux. On propose alors un schéma adapté aux métamatériaux, le schéma EF-AL, qui converge dès que le problème est bien-posé. On termine par l'étude du schéma Galerkin Discontinu dont on montre numériquement sa convergence sur des exemples de métamatériaux. La seconde partie présente l'homogénéisation non-périodique formelle de métamatériaux acoustiques. Les travaux d'A.G. Ramm sur la création de milieux à partir d'assemblages d'obstacles sont repris afin de préciser l'asymptotique fine du comportement du champ diffracté par un nombre fini de petites boules de rayon \delta. On utilise pour cela la méthode des développements asymptotiques raccordés. On établit l'existence et l'unicité de ce dernier et des estimations d'erreurs qui valident l'approche formelle. On suppose ensuite que le nombre de petits objets tend vers l'infini lorsque \delta tend vers 0 et passons à la limite dans le développement. Une approximation de Born permet d'obtenir l'indice du milieu contenant tous les objets qui, dans certains cas, est celui d'un métamatériau.

Dans cette thèse, nous souhaitons décrire la dynamique du front d'avancement d'un film mince s'écoulant sur un plan incliné non rugueux. Nous nous intéressons surtout au problème de point triple situé à l'interface entre la paroi solide, le fluide en mouvement et l'air, par exemple lors de l'écoulement d'une goutte sur une surface inclinée. Dans une première partie, nous expliquons pourquoi on peut se ramener aux équations de Stokes et pourquoi le problème résultant est mal posé. Pour y remédier, la condition de non-glissement à la paroi est remplacée par une condition de glissement lorsqu'on est proche du front. Ainsi on réussit à trouver une solution dans H1. Puis nous développons la dynamique de l'écoulement à l'amont du front : un film mince. Cet écoulement peut se modéliser sous la forme d'équations de type Saint-Venant sur la hauteur et le débit. Nous justifions cette construction à partir des équations de Navier-Stokes en utilisant un développement asymptotique en fonction du paramètre onde longue. Dans la zone du front nous résolvons le système de Stokes stationnaire avec glissement au fond par un développement asymptotique en fonction du nombre capillaire. Le front est divisé en une zone interne près du front et une zone externe loin du front, puis les solutions de chaque zone sont soit raccordées directement (angles dynamique et statique égaux), soit raccordées au moyen d'une zone intermédiaire (angles dynamique et statique différents). Cela nous conduit à deux familles de modèles. En réunissant les modèles type Saint-Venant et les différents modèles de front, nous obtenons un modèle de Saint-Venant tenant compte de la dynamique du front. À partir de ce modèle à deux équations nous pouvons écrire un modèle plus simple à une équation sur la hauteur. Ce modèle permet d'étendre les modèles existants avec adhérence à des modèles avec glissement. On peut alors réaliser des simulations numériques combinant un front d'avancement et un film mince.

De cette thèse, une nouvelle méthode éléments finis hybride MAX-FEM dédiée à la modélisation thermomécanique des structures avec couches minces a été développée. Cette nouvelle approche se base sur un couplage analytique-numérique de deux méthodes : les Développements Asymptotiques Raccordés (MAE) et la Partition de l'Unité (PUM). Ce couplage consiste à construire l'enrichissement de la PUM par MAE est mène à une forme corrigée de la méthode des éléments finis classique (FEM). Cette correction est obtenue à travers des matrices de correction contenant les informations géométriques et caractéristiques du matériau de la couche mince. Les matrices introduites par l'approche MAX-FEM simplifient son implémentation numérique sous différents codes de calculs (MATLAB, ABAQUS, ...) et permettent l'obtention de la solution globale en un seul calcul. Les résultats obtenus par la MAX-FEM pour des applications 1D et 2D thermomécaniques montrent une très bonne précision avec un temps de calcul minimal et sans raffinement de maillage. De plus, la MAX-FEM surmonte les limitations de la MAE ainsi que celle de la PUM en termes de nombre de calculs, de la sensibilité aux propriétés des matériaux, des conditions aux limites ainsi que l'intégration numérique. Finalement, l'approche MAX-FEM est exploitée pour le développement d'un nouveau protocole expérimental dédié à la caractérisation thermique des couches minces. Ce protocole vise l'identification, de manière simple, de la conductivité thermique de la couche mince après son élaboration et sous les deux régimes transitoire et permanent. L'approche consiste à confronter la nature du transfert thermique d'une éprouvette homogène à une contenant une couche mince. La différence relevée est directement liée à la conductivité thermique de la couche mince. Les résultats obtenus, après réalisation du banc d'essais, montrent une bonne précision de l'approche avec une méthodologie de mesure simple à mettre en oeuvre / This work introduces a new simplified finite elements method MAX-FEM based on hybrid analytical-numerical coupling. This method is intended to the multi-scales analysis of transient thermomechanical behavior of mediums containing thin layers such as bounded and coated structures. The MAX-FEM consists in correcting the classical Finite Elements Method (FEM) by correction matrices taking into account the presence of thin layers without any mesh refinement. The proposed correction is based on the analytical approach of Matched Asymptotic Expansions (MAE) and the numerical method of Partition of Unity Method (PUM). The developed approach can easily implemented under different numerical codes (MATLAB, ABAQUS, ...) and can be used to perform mechanical, thermal and thermomechanical analyses of 1D and 2D bounded and coated structures. The obtained results show a good accuracy with short computation time, and without any required mesh refinement. Also, the developed method overcomes the limitation of the MAE and PUM methods by exploiting the advantages of their coupling. Finally, the MAX-FEM approach was also used to develop an experimental test bench intended to the thermal characterization of thin layers. Indeed, a simple confrontation between the heat transfer in an homogeneous structure and a second structure with thin layer allows identifying the thermal conductivity in both transient and stationary regimes. The test bench is simple to release and the obtained results for brazed structure show a good accuracy of the developed approach.