Les chocs induits par laser de puissance permettent d'investir le comportement hautement dynamique des matériaux, d'un grand intérêt tant pour la recherche fondamentale que pour l'industrie. L'évolution des technologies laser ces dernières années a permis d'accéder à des régimes plus courts, en dessous de la picoseconde. L'objectif de ce travail, résultat d'une collaboration entre l'institut P', le PIMM et le CEA‐DAM est de caractériser le comportement sous choc de matériaux métalliques (Aluminium, Tantale,...) dans ce régime ultra‐bref, conduisant à des sollicitations dynamiques extrêmes (>107s‐1). L'étude repose sur la comparaison et la validation de modèles numériques à des résultats expérimentaux obtenus sur la chaîne 100TW du LULI. Cette caractérisation passe dans un premier temps par l'étude de l'interaction laser‐matière afin caractériser le chargement équivalent en pression sur la cible. Les processus en régime ultra‐bref sont différents de ce qui est connu en régime nanoseconde : en effet, l'échelle de temps, quelques picosecondes, est du même ordre que bon nombres de phénomènes moléculaires tel que le déséquilibre électrons‐ions. Ensuite, nous avons étudié l'évolution de l'onde de choc et son amortissement, très prononcé dans ce régime. L'écaillage dans une telle configuration se produit par couches très minces (quelques μm) et régulières dans ce régime. L'endommagement obtenu est caractérisé par la mesure VISAR. Les résultats obtenus par observations post‐mortem jusqu'à présent montrent que plus l'épaisseur de cible est faible, plus l'épaisseur d'écaille diminue, pouvant atteindre l'échelle du micron. Dans le cadre de la modélisation de l'endommagement et le dimensionnement des critères d'endommagement utilisés et éprouvés en régime nanoseconde (Kanel), des essais à différentes épaisseurs de cible ont été réalisés afin d'observer les conséquences d'une variation de vitesse de déformation sur l'endommagement, et généraliser le modèle de Kanel au régime ultra‐bref, et plus généralement en fonction de la vitesse de déformation. L'ensemble des résultats relatifs à l'endommagement est généralisé à des configurations 2D, permettant notamment de caractériser l'évolution du diamètre d'écaille. En parallèle, des simulations microscopiques par dynamique moléculaire de choc laser ultra‐bref sur des cibles monocristallines de Tantale à l'échelle du micron ont été menées au CEA‐DAM donnent un point de vue complémentaire des processus microscopiques liés à l'endommagement à des vitesses de déformation aux abords de la limite de cohésion théorique. Les résultats obtenus sur les chocs ultra‐brefs et 2D présentent un grand intérêt pour le développement du test d'adhérence de revêtements par choc laser (LASAT), offrant la possibilité de nouvelles extensions pour le procédé LASAT. Par exemple, de nombreux domaines industriels utilisent des revêtements micrométriques (optique, électronique, ...) mais il existe peu de méthodes pour caractériser leurs propriétés avec fiabilité. Des essais de transposition de LASAT en régime femtoseconde sur des cellules photovoltaïques ont démontré la possibilité d'éjecter des revêtements submicrométriques et caractériser leur seuil d'adhérence.

Cette thèse porte sur l'étude du remplissage rapide de réservoirs d'hydrogène sous haute pression (700 bars). Durant cette phase, se produisent simultanément une augmentation de la température du gaz et des contraintes internes dues à la pression. Un banc de fatigue thermomécanique instrumenté a été développé afin de simuler ce remplissage rapide sur des éprouvettes simples composées d'un liner en polyuréthane et d'une coque composite carbone/époxy, matériaux constitutifs des réservoirs étudiés. La validation thermique de ce banc a permis de montrer que le polyuréthane est une bonne barrière thermique. Une première étude sur un drapage composite non optimisé a permis de montrer l'influence néfaste de la température et d'un palier de maintien à la charge maximale sur le comportement du multicouche polymère/composite ainsi que l'effet bénéfique de l'alternance des plis du stratifié composite ; des couches épaisses conduisent à une apparition plus précoce de l'endommagement et à des durées de vie plus faibles. Un drapage ‘‘représentatif'' des conditions de service du réservoir a été optimisé par calcul analytique et par éléments finis afin d'atteindre les mêmes niveaux de contraintes maximales dans l'éprouvette sollicitée sur le banc de fatigue thermomécanique que ceux atteint dans le réservoir sous pression. Ces éprouvettes ‘‘représentatives'' conduisent néanmoins à des contraintes de cisaillement plus élevées que dans le réservoir et présentent des effets de bord qui n'existent pas dans le réservoir bobiné. Comme pour les drapages non optimisés, la fatigue thermomécanique conduit à des durées de vie plus faibles qu'en fatigue purement mécanique à 1Hz. L'observation par microscopie optique et par radiographie X des échantillons sollicités montre, en fatigue thermomécanique, une localisation des endommagements dans la partie la plus chauffée du composite alors qu'en fatigue mécanique à 1Hz, les endommagements sont répartis sur toute la longueur des éprouvettes. De plus, la fatigue thermomécanique conduit préférentiellement à des délaminages par rapport aux résultats obtenus en fatigue mécanique à 1Hz. Des mesures de champs de déplacement par corrélation d'images ont permis de mettre en évidence des concentrations de cisaillement sur le drapage ‘‘représentatif'' ainsi qu'un comportement viscoélastique, accentué en présence de cyclage thermique, qui seraient à l'origine des endommagements constatés.

Dans une première partie, notre étude porte sur les milieux faiblement continus avec une approche basée sur la géométrie non riemannienne. Nous considérons un corps solide déformable modélisé par une variété riemannienne munie d'une connexion affine. Un tel modèle est une extension d'un autre qui considère une connexion euclidienne, laquelle dérive du tenseur métrique imposé par l'espace ambiant. La masse par unité de volume peut être supposée non constante et le corps peut contenir des défauts décrits par des champs de discontinuité de champs scalaires ou de champs vectoriels définis sur la variété. Dans ce cas, en plus du tenseur métrique, nous utilisons nécessairement la torsion introduite par Cartan ou la courbure, deux tenseurs associés à la connexion affine. Nous disposons ainsi d'un modèle de milieu continu du second gradient. Les investigations prennent en compte les effets de ces champs tensoriels dans l'analyse de la déformation du corps. Comme application, nous décrivons l'atténuation spatiale d'une onde propagée dans un milieu non homogène. Dans une seconde partie, notre étude porte sur une modélisation de l'indentation Vickers du verre. Nous considérons un modèle qui utilise le schéma d'inclusion d'Eshelby dans une matrice semi-infinie, pour analyser les champs de contrainte et de déplacement durant le processus de charge. L'objectif est de déterminer la densification du verre sous l'indenteur. Les résultats semi-analytiques obtenus sont confrontés de manière positive avec des données expérimentales fournies par le LARMAUR.

La simulation d'éclairage peut être utilisée pour l'étude et l'analyse du confort visuel ou de la performance de dispositifs d'éclairage. Pour répondre à de tels objectifs, les méthodes utilisées doivent résoudre de manière précise et réaliste la problématique de l'illumination globale. De plus, les logiciels de simulation d'éclairage doivent souvent manipuler des scènes géométriquement complexes mais aussi exploiter les propriétés photométriques réalistes des sources artificielles étendues, des sources naturelles et des matériaux. L'objectif du travail de thèse est d'étendre les possibilités de ces outils à la prise en compte d'environnements enfumés dans lesquels la densité et la répartition des fumées évoluent avec le temps, tout en considérant le déplacement d'un observateur virtuel dans la scène. De telles possibilités ouvriraient le champ des possibilités de la simulation d'éclairage à des cas d'étude de la vision dans les fumées pour la sécurité incendie par exemple. Suite à une analyse globale du problème (interaction lumière/matériaux, lumière/fumée, évolution de la fumée dans le temps), le travail de recherche est décomposé en trois parties. Nous présentons dans un premier temps une nouvelle méthode de résolution de l'illumination globale pour les objets surfaciques basée sur la méthode de cache d'éclairement avec des enregistrements dont les zones d'influence s'adaptent à la géométrie et aux variations d'éclairage. Nous appelons ces enregistrements des \textit{enregistrements adaptatifs} Cette technique permet de contrôler plus finement la densité du cache. Par la suite, les travaux s'intéressent en détail à la problématique des milieux participatifs statiques et de leur interaction avec la lumière. Une méthode de résolution, s'appuyant sur les travaux de la première partie, est alors proposée. Les enregistrements adaptatifs sont créés dans l'espace en fonction des caractéristiques de la fumée (coefficients de diffusion et d'absorption) et de son influence sur l'éclairage global. Enfin, l'aspect dynamique est étudié et une extension temporelle de la méthode de simulation d'éclairage en présence de milieux participatifs est alors proposée. Nous introduisons le concept d'enregistrements adaptatifs spatio-temporels (pour les surfaces et les volumes) pour interpoler les variations d'éclairement à la fois dans l'espace et dans le temps.

This work is dedicated to the control of the accuracy of computational simulations of sound propagation and scattering. Assuming time-harmonic behaviour, the mathematical models are given as boundary value problems for the Helmholtz equation <i>Delta u+k2u=0 </i> in <i>Oméga</i>. A distinction is made between interior, exterior and coupled problems and this work focuses mainly on interior uncoupled problems for which the Helmholtz equation becomes singular at eigenfrequencies. As in other application fields, error control is an important issue in acoustic computations. It is clear that the numerical parameters (mesh size h and degree of approximation p) must be adapted to the physical parameter k. The well known ‘rule of the thumb’ for the h version with linear elements is to resolve the wavelength <i>lambda=2 pi k-1</i> by six elements characterising the approximability of the finite element mesh. If the numerical model is stable, the quality of the numerical solution is entirely controlled by the approximability of the finite element mesh. The situation is quite different in the presence of singularities. In that case, <i>stability</i> (or the lack thereof) is equally (sometimes more) important. In our application, the solutions are ‘rough’, i.e., highly oscillatory if the wavenumber is large. This is a singularity inherent to the differential operator rather than to the domain or the boundary conditions. This effect is called the <i>k-singularity</i>. Similarly, the discrete operator (“stiffness” matrix) becomes singular at eigenvalues of the discretised interior problem (or nearly singular at damped eigenvalues in solid-fluid interaction). This type of singularities is called the <i>lambda-singularities</i>. Both singularities are of global character. Without adaptive correction, their destabilizing effect generally leads to large error of the finite element results, even if the finite element mesh satisfies the ‘rule of the thumb’. The k- and lambda-singularities are first extensively demonstrated by numerical examples. Then, two <i>a posteriori</i> error estimators are developed and the numerical tests show that, due to these specific phenomena of dynamo-acoustic computations, <i>error control cannot, in general, be accomplished by just ‘transplanting’ methods that worked well in static computations</i>. However, for low wavenumbers, it is necessary to also control the influence of the geometric (reentrants corners) or physical (discontinuities of the boundary conditions) singularities. An <i>h</i>-adaptive version with refinements has been implemented. These tools have been applied to two industrial examples : the GLT, a bi-mode bus from Bombardier Eurorail, and the Vertigo, a sport car from Gillet Automobiles. As a conclusion, it is recommanded to replace the rule of the thumb by a criterion based on the control of the influence of the specific singularities of the Helmholtz operator. As this aim cannot be achieved by the <i>a posteriori</i> error estimators, it is suggested to minimize the influence of the singularities by modifying the formulation of the finite element method or by formulating a “meshless” method.

L'étude expérimentale présentée vise à étudier l'effet de la température sur le comportement hydraulique et poromécanique d'un calcaire. de nombreux essais (mesures de porosité et de perméabilité, essais de compressions uniaxiales et hydrostatiques) ont été menés afin d'évaluer l'influence des cycles thermiques et donc de la microfissuration thermique sur le comportement de la roche. En parrallèle, un dispositif de mesure de la perméabilité sous températures élevées (jusqu'à 200°C) a été conçu, permettant ainsi d'évaluer directement la variation de la perméabilité du calcaire sous sollicitations thermiques. Par la suite, notre étude a été étendue au cas des matériaux cimentaires, l'objectif étant d'analyser l'effet de la température sur la perméabilité d'un mortier normalisé.

Les matériaux vitreux mous (émulsion concentrée, mousse, suspension concentrée...) présentent un comportement rhéologique entre solide et liquide. Aux petites contraintes le système reste élastique, mais au-dessus d’une contrainte seuil, le système s’écoule comme un liquide visqueux. Ce comportement se trouve dans de nombreux fluides industriels comme dans les cosmétiques, l’agro-alimentaire ou encore le béton. Une contrainte seuil n’apparait que au dessus d’une certaine fractions volumique. Au dessus de cette fraction les particules se bloquent entre elle, la relaxation n’est plus possible et l’écoulement devient fortement coopératif.Cette coopérativité influe sur la rhéologie à petite échelle, Quand le confinement devient de l’ordre de quelques particules, la viscosité ne dépend plus uniquement, comme habituellement, de la contrainte locale mais aussi de la contrainte au voisinage. C’est ce qu’on appelle rhéologie non-localeJ’ai étudié expérimentalement ce comportement en utilisant les outils de micro fluidiques. J’ai étudié une micro-émulsion concentrée s’écoulant dans un microcanal en observant directement l’écoulement des gouttes avec un microscope confocal. Les résultats sont comparés au model “Kinetic-Elasto-Plastic” de Bocquet et al. 2009 et à des simulations de dynamique moléculaire. / Soft glassy materials (concentrated emulsion, foams, concentrated suspension…) present rheological properties between solids and liquid. Under small stress they stay elastic but at stress higher than a yield stress they begin to flow as a liquid. Those fluids are used in cosmetics, food industry or building materials as concrete. The yield stress behavior only appears when the volume fraction is high enough, where the particles are blocked by their neighbors. So the systems cannot relax and the flow become highly cooperative.This cooperativity impacts the rheology at small scale. When the confinement is of the order of few particles, the viscosity does not only depend on the local stress as usually but also on the stress in the neighborhood. This is called non-local rheology.I studied experimentally this behavior by flowing concentrated emulsion in a microchannel and observing directly the flow of the droplet with a confocal microscopy. The results from these microfluidics experiments are compared to predictions of the Kinetic Elasto Plastic model of Bocquet et al. 2009 and molecular dynamic simulation of jammed soft particles.

Le travail actuel traite le problème de l'écoulement gaz-liquide compositionnel pour la représentation d'un puits dans des simulations de réservoir. L'objectif est de développer les rapports analytiques qui pourraient lier la pression de puits, la saturation et les concentrations de composant à leurs valeurs moyennes dans chaque zone de l'influence de puits. Nous avons montre que des N-2 équations décrivant le transport des concentrations de phase peuvent être transformées en équations ordinaires (différentiation en ce qui concerne la pression) indépendantes du temps et de l'espace examinant le long des lignes de courant. Ces équations transformées représentent des relations thermodynamiques additionnelles réduisant le degré de liberté thermodynamique. En raison de ceci le variance thermodynamique du modèle compositionnel limite s'avère égal à 1 pour tout nombre de composants chimiques. Cette transformation assurent se découplage total du modèle compositionnel limite dans le nouveau modèle thermodynamique et le modèle hydrodynamique, qui peut être resoved inedpendently d'un un autre. Le modèle thermodynamique décompose est totalement indépendant sur l'hydrodynamique, et décrit le comportement d'équilibre d'un système gaz-liquide ouvert. Ce modèle contient les équations d'équilibre et la EOS classiques, aussi bien que les N-2 nouvelles équations appelées la "Delta-loi", qui déterminent la variation de composition d'un système ouvert dans lequel la masse de chaque composant n'est pas conservée. Le modèle hydrodynamique décompose a été utiliser pour développer les solutions asymptotiques des problèmes d'écoulement de gaz-condensat. Le problème a été montré perturbé singulièrement avec la formation d'une couche limite à voisinage du puits. Dans cette couche la propriété basique de contraste des mobilities de gaz et de liquide est perturbée. On développe une technique spéciale qui permet de construire des expansions asymptotiques sous forme de deux diverses séries: le primer est valide loin du puits (l'expansion extérieure), alors que le second dans valide à voisinage du puits (la couche limite ou l'expansion intérieure). En appliquant la méthode asymptotique suggérée, nous avons développé les solutions asymptotiques pour le problème de l'écoulement multicompositionnel de gaz-condensat àu puits dans un domaine borné à un débit variable. En plusieurs cas la solution peut être obtenue sous la forme analytique, alors qu'en cas général de l'écoulement la méthode mène à une solution de semi-analytical, présentée comme problème initial pour une équation. Cette solution, même étant présenté en forme non-analytique, est beaucoup plus simple que le modèle compositionnel original, car l'équation pour la saturation ne dépend pas de la pression locale, mais dépend de la pression de bord seulement. Dans le dernier chapitre nous avons prolongé cette approche au cas quand la pression capillaire n'est pas négligée. Nous avons supposé cependant que les forces capillaires sont inférieures à la différence de pression entre le puits et la bord de réservoir, dus à ce que nous avons appliqué la méthode de perturbation pour petit nombre capillaire inverse. On obtient les solutions asymptotiques améliorées qui tiennent compte de l'effet capillaire. Simulations numériques montrées que ces effets sont maximaux àu voisinage du puits. Le cas d'une exploitation à long terme du réservoir. Tout d'abord, la simulation traditionnelle du comportement de réservoir peut être effectuée avec l'ECLIPSE en ajoutant la méthode de Peaceman de représentation bonne, qui est une relation analytique pour la pression de puits par l'intermédiaire du débit de production. Cette relation inclut une saturation condensat qui peut être évaluée comme saturation moyenne de réservoir. Une telle simulation fournit un bon résultat pour la pression de puits (ou le débit de production), et un bon résultat pour la saturation de bord, mais des données faibles pour la saturation de puits. Cette valeur peut être calculée en utilisant les solutions asymptotiques suggérées dans le présent projet. Le cas d'un puits de production à court terme (un essai de puits). Il est suffisant de simuler le comportement de réservoir dans le domaine de l'influence de puits, en supposant que l'état de frontière demeure invariable (et connu a priori). Dans ce cas-ci les solutions asymptotiques suggérées dans le travail de presnet peuvent être directement employées pour simuler le problème (sans employer l'ECLIPSE). Le problème de l'écoulement de gaz-condensat à une fracture. Nous avons construit un champ plutôt arbitraire avec des lignes de courant orientées à la fracture, en supposant que la fracture joue le rôle d'une décharge, et les lignes de courant sont stationnaire. Pour une ligne de courant arbitraire nous avons reformulé le modèle d'écoulement de gaz-condensat dans des coordonnées cartésiennes. Pour ce problème nous avons développé les expansions asymptotiques / The present work deals with the problem of the compositional gas-liquid flow for the well representation in reservoir simulations. The objective is to develop analytical relationships which would be able to link the wellbore pressure, saturation and component concentrations to their mean values within each zone of the well influence. It is shown that N-2 equations describing the transport of phase concentrations can be transformed into the space- and time-independent ordinary differential equations (differentiation with respect to pressure) when examined along flow streamlines. These transformed equations represent additional thermodynamic relations reducing the thermodynamic degree of freedom. Due to this the thermodynamic variance of the limit compositional model is shown to be equal to 1 for any number of chemical components. This transformation ensure a total splitting of the limit compositional model into the new thermodynamic model and a hydrodynamic model, which may be resoved inedpendently of one another. The split thermodynamic model is totally independent on the hydrodynamic one, and describes the equilibrium behaviour of an open gas-liquid system. This model contains the classic equilibrium equations and EOS, as well as N-2 new differential equations called the "delta-law" which determine the composition variation in an open system, in which the mass of each component is not conserved. The split hydrodynamic model consists of two equations for pressure and saturation. The split hydrodynamic model was used to develop asymptotic solutions of gas-condensate flow problems. The problem was shown to be singularly perturbed with formation of a boundary layer in the vicinity of the well. In this layer the basic contrast property of gas and liquid mobilities is perturbed. A special technique is developed which enables to construct asymptotic expansions in the form of two various series, one of them is valid far from the well (the exterior expansion), while the second one in valid in the vicinity of the well (the boundary-layer or interior expansion). By applying the suggested asymptotic method, we have developed the asymptotic solutions for the problem of multi-component gas-condensate flow to a well in a bounded domain at a variable flow rate. In several cases the solution may be obtained in the analytical form, while in general case of flow the method leads to a semi-analytical solution presented as an initial problem for a differential equation. This solution, even being presented in non-analytical form, is much simpler than the original compositional model, as the equation for saturation does not depend on the local pressure, but on the boundary pressure only. In the last chapter we extended this approach to the case when the capillary pressure is not neglected. We assumed however that the capillary forces are lower than the pressure difference between the wellbore and reservoir boundary, due to which we applied the perturbation method over the small inverse capillary number. The improved asymptotic solutions are obtained which take into account the capillary effect. Numerical simulations shown that these effects are maximal in the vicinity of the well. For the practice, the obtained asymptotic solutions may be used in the following way to resolve the problem of gas-condensate well representation. The case of a long-term exploitation of the reservoir}. First of all, the traditional simulation of the reservoir behaviour can be performed with ECLIPSE by adding the Peaceman method of well representation, which is an analytical relation for the wellbore pressure via the production rate. This relation includes a condensate saturation which can be evaluated as a mean reservoir saturation. Such a simulation provides a good result for the wellbore pressure (or the production rate), and a good result for the boundary saturation, but poor data for the wellbore saturation. This value can be calculated next by using the asymptotic solutions suggested in the presented project. The case of a short-term well production (a well test). It is sufficient to simulate the reservoir behaviour in the domain of the well influence, by assuming that the boundary state remains invariable (and known a priori). In this case the asymptotic solutions suggested in the presnet work can be directly used to simulate the problem (without using ECLIPSE)

Les films minces sont abondamment utilisés pour protéger ou fonctionnaliser les surfaces. Les contraintes mécaniques peuvent cependant entraîner la ruine de ces structures par la propagation de fissures ou leur décollement. Nous présentons deux systèmes expérimentaux qui amènent à s'interroger sur les limites de validité des critères de stabilité communément admis. Pour une couche de sol- gel présentant une énergie d'adhésion modérée avec son substrat de silicium, nous observons la propagation de fissures dans des films réputés stables. Les motifs de propagation sont étonnement réguliers (en forme de spirales ou d'allées de croissants). Nous caractérisons expérimentalement cette propagation. Nous mettons en évidence un nouveau mode de propagation de fissure où la délamination et la fissure transverse collaborent. Nous montrons l'existence d'une taille d'interaction robuste entre fissures indépendantes du chargement et des énergies de rupture. Nous proposons un modèle simple permettant de prédire les conditions de propagation et la sélection des motifs. Nous caractérisons enfin la dynamique de propagation de telles fissures. Dans un second système expérimental motivé par les applications en électronique flexible, nous nous intéressons à des films nanométriques d'oxyde de zinc déposés sur des substrats de polymère. La rigidité du substrat est alors faible devant celle de la couche. Nous décrivons la fragmentation du film et proposons une extension d'un modèle de couche élastique sur fondation élasto-plastique. Enfin, nous nous intéressons à la compression et à la fragmentation d'un multicouche comportant une fine épaisseur d'argent. / Nanometer to micrometer thin film coatings are extensively used in material science to protect and functionalize surfaces. Coating processes however generally result into residual strains and may induce the propagation of cracks with dramatic consequences on film reliability. We present two experimental models which challenges the usual stability criteria for thin films. For a thin sol-gel layer moderatly adherent to a silicon substrate, we observe fracture propagation in coatings considered stable. Moreover, intriguing cracks of regular geometry are found to self-replicate an initial triggering pattern independently from the local material properties. We present a novel mode of fracture which involves the simultaneous progression of both, a crack and a front of delamination. We describe the fracture mechanism with universal energy arguments and provide a simple model describing the crack path. We eventually characterize the dynamics of the propagation. In a second model experiment, we study the fragmentation of coating deposited on a compliant substrate (motivated by flexible electronics application). We describe the fragmentation and propose an extension to an elastic layer deposited on an elastic foundation to an elasto-plastic case. Eventually, we study the compression and fragmentation of a zinc oxyde and silver multilayer

Lorsqu'un fluide complexe s'écoule à travers un milieu poreux, à la non-linéarité de l'écoulement s'ajoute la spécificité de la structure poreuse, qui est souvent multi-échelle. Il émerge alors un grand nombre de problématiques fondamentales liées à l'interaction entre le fluide et la structure poreuse. L'interprétation et la modélisation de la grande variété des phénomènes physiques à petite échelle ainsi que leurs répercussions à grande échelle soulèvent de nombreuses questions. Dans cette thèse, les fluides étudiés sont des solutions de polymères, et les milieux poreux sont, entre autre, des roches issues de réservoirs pétroliers. Dans le contexte des méthodes de récupération améliorée pour les gisements pétroliers, l’injection d’eau avec polymères fait en effet partie des méthodes couramment utilisées, permettant d’augmenter l’efficacité du balayage et donc la production d’huile sur différents types de réservoirs. De part la rhéologie non-Newtonienne ainsi que les phénomènes particuliers proches de la paroi que développent les molécules de polymères, les fluides modélisés dans cette thèse peuvent être qualifiés de complexes. L'objectif de cette thèse est d'étudier la rhéologie non-Newtonienne ainsi que le comportement des molécules de polymère proches de la paroi. On relie alors ces phénomènes aux propriétés effectives macroscopiques. Pour cela, on simule numériquement les écoulements à travers desimages tomographiques de milieux poreux. Ainsi, on souhaite mieux cerner la physique qui est en jeu et également proposer des pistes d'amélioration des modèles actuellement implémentés dans les simulateurs de réservoirs. Dans une première partie de ce travail de thèse, on s'intéresse à la transition du régime d'écoulement macroscopique, de Newtonien à non-Newtonien, induite parune solution de polymère. Par des simulations numériques de l'écoulement à travers un large panel de milieux poreux, on étudie la transition entre ces deux régimes. Une analyse de la mécanique de l'écoulement permet de proposer un modèle simple et d'évaluer en ordre de grandeur la vitesse moyenne de transition. Ensuite, on étudie le glissement apparent induit par un mécanisme de répulsion des chaînes de polymère à la paroi liquide/solide. On propose un modèle à l'échelle des pores de ce phénomène et, par comparaison avec des données expérimentales,on montre que ce modèle permet de retrouver avec une précision acceptable les comportements macroscopiques observés. Enfin, avec des simulations directes sur des milieux périodiques, on relie les phénoménologies micro- et macroscopiques d'écoulements non-Newtoniens. D'un point de vue fondamental, on étudie notamment la compétition entre la non-linéarité induite par la rhéologie non-Newtonienne et le désordre inhérent à la structure poreuse. Les modèles actuellement utilisés dans les simulateurs de réservoirs sont reconsidérés au vu des résultats. / When a complex fluid is flowing through a porous medium, in addition to the fluid intrinsic physics, the multi-scale properties of the porous structure play a significant role. From the interaction between these features arise a great number of complex physical phenomena. The understandingand the modelling of the variety of these phenomena involved at the small scales and their impact on the large scales is the subject of intense work. In this thesis, the fluids we consider are polymer solutions and the porous media are typical sandstones met in petroleum applications. In petroleumengineering, the injection of a polymer slug into the oil-bearing reservoir is indeed a method commonly used in enhanced oil recovery. This method allows to increase the sweep efficiency, hence to improve the oil production of the reservoir. Due to the non-Newtonian rheology induced by the polymer molecules as well as specific mechanisms occurring at the liquid/solid interface, a polymer solution may be qualified as a complex fluid. The goal of this thesis is to investigate the non-Newtonian rheology and the behaviour of the polymer molecules near the liquid/solid interface. These phenomena are related to macro-scale effective properties. To thoroughly address this goal, we perform numerical simulations of flows through porous media. The goal is to obtain a better understanding of the underlying physics and, furthermore, we wish to propose possible improvements of the models that are currently used in reservoir simulators. Primarily, we are interested in the macro-scale transition from a Newtonian to a non-Newtonian flow. This transition is induced by the non-Newtonian rheology. By simulating flows through a wide panel of porous media, we study the macro-scale transition. An analysis of the fluid mechanics involved allows us to propose a simple model for the critical average velocity at which the transition occurs. In addition, we study the apparent slip induced by a repulsive mechanism of the polymer chains from the liquid/solid interface. We propose a pore-scale model to this mechanism. By performing comparisons with experimental datasets, we show that the model allows for a good description of the observed macro-scale behaviour. Finally, with direct simulations over periodic media, we link the phenomenology at the micro- and macro-scale for the flow of non-Newtonian fluids. On a fundamental level, we study the competition that emerges between the nonlinearity induced by the rheology and the disorder inherent to the porous structure. The results are related to modelscommonly used in reservoir simulators.