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

La thèse porte sur l'apport de la méthode de Laplace pour l'approximation du filtre bayésien dans des modèles de Markov cachés généraux, c'est-à-dire dans un cadre séquentiel, avec comme domaine d'application privilégié la poursuite de cibles mobiles. A la base, la méthode de Laplace est une méthode asymptotique pour le calcul d'intégrales, c'est-à-dire dans un cadre statique, valide en théorie dès que la fonction à intégrer présente un maximum de plus en plus significatif, lequel apporte la contribution essentielle au résultat. En pratique, cette méthode donne des résultats souvent très précis même en dehors de ce cadre de validité théorique. Les deux contributions principales de la thèse sont les suivantes. Premièrement, nous avons utilisé la méthode de Laplace en complément du filtrage particulaire : on sait en effet que les méthodes de Monte Carlo séquentielles basées sur l'échantillonnage pondéré sont mises en difficulté quand la fonction de pondération (ici la fonction de vraisemblance) est trop localisée, par exemple quand la variance du bruit d'observation est trop faible, or c'est précisément là le domaine où la méthode de Laplace est efficace et justifiée théoriquement, d'où l'idée naturelle de combiner les deux points de vue. Nous proposons ainsi un algorithme associant la méthode de Laplace et le filtrage particulaire, appelé le Laplace particle filter. Deuxièmement, nous avons analysé l'approximation du filtre bayésien grâce à la méthode de Laplace seulement (c'est-à-dire sans génération d'échantillons aléatoires) : il s'agit ici de contrôler la propagation de l'erreur d'approximation d'un pas de temps au pas de temps suivant, dans un cadre asymptotique approprié, par exemple quand le bruit d'observation tend vers zéro, ou quand le bruit d'état et le bruit d'observation tendent conjointement (et à la même vitesse) vers zéro, ou plus généralement quand l'information contenue dans le système tend vers l'infini, avec une interprétation en terme d'identifiabilité. / The thesis deals with the contribution of the Laplace method to the approximation of the Bayesian filter in hidden Markov models with continuous state--space, i.e. in a sequential framework, with target tracking as the main application domain. Originally, the Laplace method is an asymptotic method used to compute integrals, i.e. in a static framework, valid in theory as soon as the function to be integrated exhibits an increasingly dominating maximum point, which brings the essential contribution to the integral. The two main contributions of the thesis are the following. Firstly, we have combined the Laplace method and particle filters: indeed, it is well-known that sequential Monte Carlo methods based on importance sampling are inefficient when the weighting function (here, the likelihood function) is too much spatially localized, e.g. when the variance of the observation noise is too small, whereas this is precisely the situation where the Laplace method is efficient and theoretically justified, hence the natural idea of combining the two approaches. We thus propose an algorithm associating the Laplace method and particle filtering, called the Laplace particle filter. Secondly, we have analyzed the approximation of the Bayesian filter based on the Laplace method only (i.e. without any generation of random samples): the objective has been to control the propagation of the approximation error from one time step to the next time step, in an appropriate asymptotic framework, e.g. when the variance of the observation noise goes to zero, or when the variances of the model noise and of the observation noise jointly go (with the same rate) to zero, or more generally when the information contained in the system goes to infinity, with an interpretation in terms of identifiability.

Cette thèse est consacrée à la modélisation, à l'analyse mathématique et à la simulation numérique d'écoulements de divers fluides complexes dans des domaines de faible épaisseur. En effet, les modèles de fluides newtoniens ne sont pas toujours suffisants pour décrire de manière réaliste les écoulements considérés. Plusieurs phénomènes peuvent être pris en compte :<br /> * le caractère complexe des fluides eux-mêmes, comme pour des fluides non-newtoniens ;<br /> * l'hétérogénéité de l'écoulement, dans le cas de mélanges de fluides par exemple.<br />Il est important d'analyser comment ces modèles peuvent être simplifiés dans le cas de domaines minces, et d'étudier rigoureusement les modèles approchés.<br />Dans une première partie, des écoulements de fluides non newtoniens visco-élastiques représentés par une loi de comportement de type Oldroyd-B couplée aux équations de Navier-Stokes sont étudiés. Dans le cas de géométries minces, un modèle approché a été proposé. On justifie la validité de cette approximation ; la démonstration repose sur des estimations et des résultats de régularité fins.<br />Dans une deuxième partie, on considère un modèle d'écoulement piezovisqueux utilisé en lubrification hydrodynamique. Ce modèle fait aussi intervenir la déformation élastohydrodynamique du domaine (déformation du type Hertz), et l'aspect diphasique de la cavitation, qui est décrit par le modèle d'Elrod-Adams (en pression-saturation). On montre l'existence d'une solution à ce problème pour des lois pression-viscosité réalistes.<br />Dans une troisième partie, on introduit un modèle diphasique à interface diffuse, permettant de rendre compte de phénomènes plus fins tels que les gouttes. Pour cela, un paramètre d'ordre est introduit (fraction volumique d'une phase dans le mélange), gouverné par le modèle de Cahn-Hilliard. Un système approché est obtenu de manière heuristique pour un domaine de faible épaisseur. On étudie les propriétés mathématiques de ce système, et on montre un résultat d'existence, avec prise en compte ou non de la tension de surface.<br />Dans la dernière partie, un schéma numérique est mis en place pour simuler le modèle décrit précédemment d'écoulements diphasiques en domaines minces. Il permet de prendre en compte différents phénomènes physiques, comme de grandes variations de la viscosité ou la présence de recirculations à l'intérieur d'une goutte, ainsi que de simuler des mélanges dans le cadre d'écoulements lubrifiés.

Dans un premier temps, nous décrivons le complexe des solutions du transformé de Mellin algébrique d'un D-module M en fonction des solutions de M. Pour cela, nous définissons un foncteur de transformation de Mellin faisceautique. Nous montrons alors que le transformé de Mellin du complexe des solutions à décroissance rapide en 0 et à l'infini d'un D-module holonome régulier M est quasi-isomorphe au complexe des solutions du transformé de Mellin algébrique de M, l'hypothèse de régularité n'étant plus nécessaire à une variable.<br />Dans un second temps, nous faisons un travail analogue avec la transformation de Mellin inverse : les résultats sont plus partiels. Nous définissons une transformation de Mellin inverse faisceautique. Nous démontrons alors qu'il existe des morphismes naturels reliant le complexe des solutions du transformé de Mellin inverse algébrique d'un module aux différences avec le transformé de Mellin inverse faisceautique du complexe des solutions à croissance au plus exponentielle d'ordre 1 à l'infini dans des bandes verticales. Nous montrons ensuite que dans le cas d'un module aux différences à une variable et à une seule pente strictement positive, ces morphismes sont des isomorphismes.

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.

Dans une première partie, nous présentons des équations de Saint-Venant. Sur le modèle proprement dit, nous remarquons tout d'abord que, suivant le lien entre la viscosité et le rapport des échelles caractéristiques, il est indispensable de conserver l'expression complète de la force de Coriolis : nous obtenons ainsi un nouveau modèle, avec un "effet cosinus". Nous montrons alors que les preuves d'existence de solutions faibles peuvent être adaptées à ce nouveau système. Des simulations numériques de certaines ondes soulignent l'importance de ce terme. Nous étudions ensuite l'influence des conditions limites (surface, fond) sur des modèles de type Saint-Venant. Nous présentons également des modèles obtenus en utilisant des échelles multiples en espace et en temps. Enfin, nous analysons théoriquement et numériquement un nouveau modèle de sédimentation puis nous donnons certains résultats pour les fluides visco-plastiques.<br />Dans une deuxième partie, nous nous intéressons aux équations limites que sont les équations quasi-géostrophiques (QG) et les équations des lacs. L'étude numérique des équations QG 2d nous permet de voir le rôle de l'effet cosinus de la force de Coriolis. En fonction de la topographie considérée, nous montrons que celui-ci peut être non négligeable. Toujours sur les équations QG, nous donnons un schéma, basé sur des développements asymptotiques, qui permet de bien capter la couche limite mais aussi d'ajouter le terme de topographie à la solution obtenue avec fond plat, sans refaire tous les calculs. Enfin, nous expliquons l'obtention des équations des lacs avec effet cosinus, et nous prouvons que les propriétés d'existence de solutions restent valables.

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.

La première réalisation de ce travail est la construction unifiée et automatique d’un milieu continu équivalent à un treillis de poutres, dans le domaine élastique, en adoptant un modèle de poutres de Bernoulli. Une extension a été réalisée au domaine plastique, selon un algorithme de suivi de la loi de comportement après écrouissage. Suivant l’ordre des développements asymptotiques choisi, on obtient pour le comportement élastique un milieu continu classique ou micropolaire. On se restreint dans ce dernier cas aux treillis à cellules élémentaires centro-symétriques. Les codes de calculs obtenus fournissent de façon automatique les lois de comportement effectives et les modules mécaniques homogénéisés. Une grande variété de treillis, existants ou originaux, a été étudiée. Les résultats ont été systématiquement comparés aux données de la littérature et vérifiés par des simulations éléments finis avec une bonne concordance. La méthode utilisée montre également une capacité à prédire et comprendre le comportement atypique de certains treillis dits auxétiques présentant des coefficients de contraction négatifs. L’homogénéisation dans le domaine plastique a été limitée aux treillis à dominante extensionnelle. Le domaine de résistance élastique a été construit pour différents treillis, et un algorithme d’évolution du comportement avec écrouissage, de type retour-radial a été conçu et implémenté dans un code dédié. Un modèle de poutre élastoplastique à écrouissage isotrope est utilisé. L’application de l’algorithme à une simulation de charge-décharge montre une bonne concordance entre le treillis homogénéisé et les simulations éléments finis / The first achievement of this work is to construct a unified and effective continuum equivalent to a lattice of beams, in the elastic domain, using a Bernoulli beam model. An extension has been done to calculate the elastic domain resistance of such lattices and to build an algorithm for monitoring the constitutive law taking into account work hardening. The choice of the asymptotic expansions leads to a classical continuous or to a micropolar elastic continuum. We restrict in this last case our study to lattices with centro-symmetric unit cells. The numerical codes developed provide the stress-strain relationship and the effective mechanical moduli. A wide variety of trusses has been studied, either existing or original, including typical geometries of foams and various auxetic lattices, exhibiting negative contraction coefficients. The results were systematically compared with data from literature and verified by finite element simulations with a good agreement. The homogenization in the plastic range has been limited to stretching dominated lattices. The equilibrium equations of the discrete asymptotic homogenization have been used to automatically obtain the elastic resistance domain for several trusses, and a return-mapping algorithm for the follow up of the stress-strain relationship including hardening has been conceived and implemented in a dedicated code. An isotropic hardening elastoplastic model of the beam has been used. The application of the algorithm to the simulation of a loading-unloading cycle shows a good agreement between the homogenized lattice and finite element simulations

Dans cette thèse, nous considérons le problème d'Euler surface libre sur un domaine à fond non plat, dans le cadre du régime d'ondes longues de faible amplitude. L'objectif est de construire, justifier et comparer de nouveaux modèles asymptotiques pour ce problème, permettant de prendre en compte les effets liés aux variations bathymétriques. En premier lieu, nous construisons rigoureusement deux classes de modèles de Boussinesq symétriques dans le cadre de deux régimes topographiques distincts, celui de faible variations bathymétriques et celui de fortes variations. Dans un second temps, nous retrouvons et discutons dans le cas de faibles variations topographiques l'approximation classique de Korteweg-de Vries, et proposons une nouvelle approximation via l'ajout de termes bathymétriques. Dans une troisième partie, ces deux modèles, ainsi que les modèles de Boussinesq construits dans la première partie, sont simulés numériquement et comparés sur des cas tests de topographie. Enfin, il est présenté une étude numérique des équations de Green-Naghdi, dont le domaine de validité physique est plus étendu, ainsi qu'une comparaison numérique de ce modèle avec les modèles précédents sur des bathymétries spécifiques.

La prédiction de la propagation des feux de forêts consiste à trouver l'évolution du front du feu. Cette thèse reprend le modèle monodimensionnel de propagation du feu sous forme d'ellipses de Richards afin de trouver l'expression intrinsèque de la vitesse de ce front. Une formulation équivalente variationnelle de ce modèle en utilisant les principes de l'optique géométrique est aussi dérivée. Un modèle tridimensionnel de la propagation du feu grâce à une homogénéisation par prise de moyennes et à l'utilisation de la thermodynamique des processus irréversibles est alors obtenu. Une simplification de ce modèle, suivie d'une réduction bidimensionnelle sur la surface gauche du sol, nous conduisent alors à un modèle bidimensionnelle, qui tient compte des principaux paramµetres de la propagation des feux de forêts. Celui-ci permet d'obtenir des formes de corrélations utilisées par le modèle des ellipses. Le caractère non local de l'intervention du flux radiatif en provenance de la zone en feu qui se trouve au-dessus de la végétation est alors introduit. Des simulations numériques, de ce modèle bidimensionnel de propagation du feu avec terme de rayonnement non local, ont enfin été réalisées.