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.

Les matelas fibreux à base de fibres de verre produits par voie humide ont des applications industrielles en tant que séparateur de batterie et matériaux d’isolation (cœur de panneaux d’isolation sous vide). Ces matériaux utilisent principalement des fibres submicroniques, relativement chères et pouvant présenter des risques pour la santé. Ce projet est une contribution à l'élaboration d'un procédé de fabrication par voie humide de matelas de fibres de verre visant à valoriser des fibres plus grossières, le produit final devant respecter un cahier des charges précis. Nous avons été amenés à étudier le comportement des fibres de verre dans les différentes étapes du procédé et à caractériser les matelas résultants. Nous nous sommes intéressés en particulier au comportement physico-chimique des suspensions aqueuses de fibres de verre. Nous avons caractérisé les propriétés de contexture, la résistance mécanique en traction, la compressibilité et la conductivité thermique des matelas fibreux. Les travaux expérimentaux nous ont permis de proposer une formulation de la composition optimisée et des conditions opératoires du procédé afin que le matelas final soit conforme au cahier des charges. Cette nouvelle composition intègre des fibres de renfort en faible quantité. Elle permet d’améliorer les caractéristiques mécaniques sans affecter les autres propriétés.Enfin, nous avons quantifié les coûts de production et les avons comparés à ceux du procédé actuellement utilisé avec des fibres grossières (production par voie sèche).Mots-clés : Génie papetier, physico-chimie, milieu poreux, fibres de verre, caractérisation. / Glass fibre - based mats produced by a wetlaid process have industrial applications as a battery separator and insulation materials (core of vacuum insulation panels). These materials are mainly made with sub-micron fibres which relatively expensive and can present a risk to health. This project is a contribution to the production of glass fibre-based mats by a wet-laid process to add value to coarser fibres, the final product should respect precise specifications. We have been led to study the behavior of glass fibers in different stages of the process and to characterize the resulting mats. We investigated especially the physico-chemical behavior of aqueous suspensions of glass fiber. We have characterized structure properties, the mechanical resistance to traction, the compressibility and the thermal conductivity of fibrous mats. The experimental work has allowed us to give a formulation of the optimized composition and operational conditions of the process so that the final mattress conforms to the specifications. This new composition includes reinforcement fibres in small quantities. It allows improving the mechanical characteristics without affecting the other properties. Finally, we quantified the production costs and compared them to those of the process currently used with coarse fiber (dry-laid production).Keywords: Paper engineering, physical chemistry, porous media, glass fibre, characterization.

Dans les études expérimentales de la micromécanique des matériaux granulaires, la mesure des forces de contact entre particules est de nos jours toujours un challenge en comparaison avec les outils et techniques bien mieux établis pour la caractérisation cinématique à l'échelle des particules. Cette thèse de doctorat s'attaque à cet ambitieux problème. L'approche proposée implique deux aspects : (i) la caractérisation expérimentale du réseau de contact et de la cinématique à l'échelle des particules, qui peut être réalisée avec des techniques d'imagerie standards ; (ii) une approche numérique capable d'exploiter ces mesures afin de déduire les forces de contact.L'une des contraintes qu'on s'était imposée était de ne s'appuyer que sur la connaissance de la géométrie des particules ainsi que du réseau de contacts pour réaliser la déduction des forces de contact. Trois techniques numériques différentes ont été proposées à cet effet : une méthode basée sur l'élasticité des contacts (CEM), une méthode basée sur la dynamique de contact (CDM) et une méthode basée sur l'équilibre élasto-plastique de l'assemblage granulaire (QSM). Chacune de ces techniques repose sur une approche de la famille des méthodes en éléments discrets ; il s'agit respectivement de le DEM de type Cundall, la dynamique des contacts non régulière, et une approche de calcul statique élastoplastique. La non-unicité de la solution est le principal problème avec les techniques choisies, et elles sont étroitement liées à l’indétermination des forces dans le système.Les trois méthodes sont d'abord présentées et validées en les appliquant à l'estimation des forces dans les systèmes granulaires 2D générés au moyen de simulations DEM explicites. Nous prenons ces simulations comme des expériences "idéales" dans le sens où elles fournissent des données similaires à celles extraites des expériences, mais dépourvues d'erreurs de mesure. Un avantage évident de cette stratégie est d’obtenir des ensembles de forces faisant office de référence faisant foi. Sur cette base, les principaux aspects affectant la détermination des forces peuvent être étudiés. En particulier, le rôle crucial de l'histoire du chargement est mis en évidence et certaines solutions pour les prendre en compte dans la détermination des forces ont été prospectées. Une évaluation de l'influence de l'erreur de mesure a également été réalisée pour prédire l'applicabilité de chaque méthode à des expériences réelles. Une brève analyse de la variabilité des solutions est également fournie.Finalement, des tentatives ont été faites pour déduire des forces issues d’expériences effectuées dans le dispositif 1gamma2epsilon. La cinématique des particules et la connectivité ont été évaluées au moyen de la technique de corrélation d'image numérique. Les avantages et inconvénients des trois méthodes ont été éclaircis. Ils nous conduisent à envisager une utilisation combinée des trois méthodes pour tirer parti de leurs atouts respectifs. À l'avenir, il conviendra de réfléchir à la prise en compte de la stabilité de la solution -- dans l'algorithme de convergence vers une solution -- avec l'espoir de limiter la variabilité des solutions. / In the experimental study of the micro-mechanics of granular materials, measuring inter-particle contact forces is still a challenging task, if compared to the well-established tools and techniques for the kinematic characterisation at particle scale. This doctoral thesis addresses this problem. The proposed approach consists of two parts: an experimental characterisation of the granular network geometry and of particle-scale kinematics, which can be carried out with common imaging techniques such as Digital Image Correlation; a numerical approach aiming to exploit these measurements for the estimation of forces.One imposed constraint was to only make use of the rigid motions of particles, together with the knowledge of the contact network, to infer contact forces. Three different numerical techniques have been proposed to this purpose, referred to as Contact Elasticity Method (CEM), Contact Dynamics-based Method (CDM) and Quasi-Static Method (QSM). Each of these techniques is based on the formulation of common approaches in the family of Discrete Element Methods, respectively the classical Cundall-like DEM, the Non Smooth Contact Dynamics and a quasi-static approach accounting for both contact elasticity and plasticity. It is shown that memory of the history of the packing is the main concern with all the chosen techniques.The three methods are first presented and validated by applying them to the estimation of forces in 2D granular systems generated by means of explicit-time DEM simulations. We refer to these simulations as "ideal" experiments since they are meant to provide the same information that can be extracted from experiments, but without any measurement error. An obvious benefit of this strategy is to get reference force sets that are taken as ground truth. Based on this, the main aspects that affect the determination of forces can be investigated. In particular, the crucial role of history is emphasised here, and some solutions to take it into account in the force inference have been investigated.An assessment of the influence of measurement error has also been carried out, to predict the applicability of each method to real experiments. A short analysis of the variability of the solutions is also provided.Finally, some attempts have been made to infer forces from experiments carried out in the 1gamma2epsilon device. Particle kinematics and connectivity have been assessed by means of the Digital Image Correlation technique.The benefits and drawbacks of the three methods have been demonstrated. They conduct us to envision a combined usage of the three methods. In the future, studying the stability of equilibrium might help reducing the variability of the solutions.