Search

A space-time description of the finite transformations of thermo-mechanical continua is developed: the use of such a four-dimensional approach guarantees the general covariance of the proposed models. The conservation equations are written in this context and a constitutive model is derived for reversible transformations. We use projection operators to obtain the space and time components of the 4D governing equations and to interpret the results. We next propose a weak formulation of the problem along with its finite-element discretization, to be solved for the finite transformations of a solid. The advantage of this description is that the integration on space and time is performed in one step. We discuss why the 4D convective coordinate system is of interest to solve the problem. Finally, we illustrate the approach with analytical examples and solve thermo-mechanical problems numerically with an implementation on FEniCS software.; Une description spatio-temporelle des grandes déformations thermomécaniques des milieux continus est développée: l'utilisation d'une telle approche quadridimensionnelle garantit la covariance générale des modèles proposés. Les équations de conservation sont écrites dans ce contexte et un modèle constitutif est dérivé pour les transformations réversibles. Nous utilisons des opérateurs de projection pour obtenir les composantes spatiales et temporelles des équations régissant 4D et pour interpréter les résultats. Nous proposons ensuite une formulation faible du problème ainsi que sa discrétisation par éléments finis, à résoudre pour les grandes déformations d'un solide. L'avantage de cette description est que l'intégration sur l'espace et le temps se fait en une seule étape. Nous discutons pourquoi le système de coordonnées convectives 4D est intéressant pour résoudre le problème. Enfin, nous illustrons la démarche par des exemples analytiques et résolvons numériquement des problèmes thermomécaniques avec une implémentation sur le logiciel FEniCS.

Rockburst is a violent explosion of rock that can occur in underground mines. In the current research, the main objective is to demonstrate the causes that may influence the rockburst proneness by using the numerical modeling tool. However, firstly, the pre-mining stress state and the induced stresses due to surrounding excavations have to be studied precisely. The Provence coal mine, where a rockburst took place in its shaft station that is surrounded by many longwall caving panels, has been chosen as a case study. A large-scale 3D numerical model has been constructed to include the shaft station area with its small-scale pillars and galleries, and the large-scale longwall panels with their accompanying goaf area. Many problems appeared while developing such large-scale numerical model, the first problem was the initialization of stress state at a large-scale, where the measured vertical stresses are in disagreement with the overburden weight, and the in-situ stresses are highly anisotropic. The second problem was the simulation of the goaf area accompanying longwall panels. The third problem was the assessment of pillars instability in terms of its strength/average stress ratio, and its volume. The Fourth problem was the assessment of rockburst proneness in the shaft station based on different rockburst criteria. Five methods were developed to initialize the heterogeneous pre-mining stress in the large-scale numerical model. These methods are based on the Simplex Method, which is mainly based on optimizing the difference between the in-situ measured stress values and the numerical stress values to develop stress gradients able to express the stress heterogeneity and compatible with the in-situ measurements. The method that is based on initiating the stress state with 3D stress gradients was found to be more efficient than the traditional method that is based on the horizontal-to-vertical stress ratios. Regarding the goaf simulation, three models were developed and implemented in the numerical model to express the mechanical behavior within the goaf area above longwall panels. Two of these models are based on an elastic behavior, and the third one is based on the strain-hardening elasto-plastic behavior that takes the consolidation phenomenon into consideration. It was found that the goaf area above longwall panels could reach up to 32 times the seam thickness, and the elastic modulus of caved area (the first few meters in the goaf area) did not exceed 220 MPa to fulfill the roof-floor convergence. But, with advance of the exploitation, this soft material consolidated under the pressure of the overlying strata. In case of critical and super-critical width, the vertical stress in the goaf area exceeded the overburden weight, and it increased up to 4 times the overburden weight on the rib-sides. The vertical stress increased in the shaft station pillars as a result of exploiting the nearby longwall panels. It was found that the pillar volume plays an important role in its stability. And, the strength/stress ratio was found to be insufficient to quantify the rockburst proneness in underground mines. Many rockburst criteria were implemented in the numerical model to assess the rockburst proneness. It was found that the criteria that are based on stress and strain changes were able to assess the rockburst proneness; Le phénomène de coup de terrain est une explosion violente de roche qui peut se produire dans les mines souterraines. Dans la présente recherche, nous avons essayé de démontrer les causes qui peuvent influer sur la prédisposition aux coups de terrain en utilisant la modélisation numérique. Cependant, avant tout, l'état de contrainte avant l'exploitation minière et les contraintes induites par les excavations environnantes doivent être étudiés avec précision. La mine de charbon de Provence, qui a subi un phénomène de coup de terrain au niveau de son puits vertical entouré de nombreux panneaux de longue taille, a été choisie comme cas d’étude. Un modèle numérique 3D à grande échelle a été construit pour inclure la zone du puits vertical avec ses piliers et galeries à petite échelle et les panneaux de longue taille à grande échelle avec leurs zones de foudroyage associées. Plusieurs problèmes ont été rencontrés lors du développement de ce modèle numérique à grande échelle. Le premier porte sur l'initialisation de l'état de contrainte à grande échelle, où les contraintes verticales mesurées divergent avec le poids des déblais et les contraintes in situ sont très anisotropes. Le deuxième porte sur la simulation de la zone de foudroyage associée aux panneaux de longue taille. Le troisième concerne l'évaluation de l’instabilité du pilier en fonction de son ratio résistance/contrainte moyenne et de son volume. Le quatrième concerne l'évaluation de la prédisposition aux coups de terrain au niveau du puits vertical en fonction de différents critères. Cinq méthodes ont été développées pour initialiser l’état de contrainte hétérogène dans le modèle numérique à grande échelle avant l’exploitation minière. Elles sont basées sur la méthode de corrélation Simplex, qui consiste à optimiser la différence entre les valeurs de contrainte mesurées in-situ et les valeurs numériques. Le but est de développer des gradients qui soient capables d'exprimer l'hétérogénéité de la contrainte et qui soient compatibles avec les mesures in-situ. La méthode basée sur l’initialisation de l'état de contrainte avec des gradients 3D s’est avérée plus efficace que celle traditionnelle basée sur les ratios de contrainte horizontale à verticale. Concernant la simulation du foudroyage, trois modèles ont été développés et intégrés dans le modèle numérique pour exprimer le comportement mécanique dans la zone de foudroyage au-dessus des panneaux de longue taille. Deux d’entre eux sont basés sur un comportement élastique alors que le troisième est basé sur un comportement elasto-plastique avec écrouissage un phénomène de consolidation. Il a été constaté que la zone de foudroyage au-dessus des panneaux de longue taille peut atteindre 32 fois l'épaisseur de la couche exploitée et que le module d'élasticité de la partie la plus endommagée de la zone foudroyée ne doit pas excéder 220 MPa pour satisfaire la convergence toit-mur. Mais, avec l'avancée de l'exploitation, ce matériau souple se compacte sous la pression des couches supérieures. Dans le cas d'une largeur critique et super-critique, la contrainte verticale dans la zone de foudroyage pourrait dépasser le poids des déblais et pourrait augmenter jusqu'à 4 fois ce poids sur les bords. La contrainte verticale a augmenté dans les piliers au niveau du puit vertical suite à l'exploitation des panneaux de longue taille à proximité. Il a été constaté que le volume du pilier joue un rôle important dans sa stabilité. Le rapport contrainte/résistance a été jugé insuffisant pour expliquer un coup de terrain. Plusieurs critères ont été intégrés au modèle numérique pour évaluer la prédisposition aux coups de terrain. Il a été constaté que les critères basés sur les contraintes et les déformations sont capables d'évaluer la prédisposition aux coups de terrain

The current evolution of the social and economic context requires from the different actors of the public works sector to adapt their practices to the challenges of sustainable development. In the field of earthworks, companies must offer technical solutions to reuse the materials located in the vicinity of the projects in order to limit the borrowing materials and unused soils deposits. Soil treatment may allow the reuse of these materials. This study was focused on four treatment products (kaolinite, bentonite, lime and cement) as well as their effects on the hydro-mechanical behavior and internal erosion resistance of a fine silt. An enhanced HET device was designed in the framework of this study to determine the internal erosion characteristics of treated soils especially with lime and cement. Beyond the characterization of treatment effects on short-term hydro-mechanical behavior of soils, the work of this study extends to the durability of treatment and the evolution of long-term hydro-mechanical behavior of treated soils subjected to drying-wetting cycles. This study of sustainability was carried out through a multi-scale approach based on laboratory study data on soil samples and field investigations on hydraulic structures made of treated soil. The study showed that hydraulic conditions variations can decrease the performance of treated soils. These degradations result on an increase in hydraulic conductivity, a decrease of the mechanical strength and also a reduction in the critical shear stress which expresses a decrease of the internal erosion resistance. The kinetic of performance loss depends to the nature of the treatment product and percentage used as well as the exposure level, the number and amplitude of the hydraulic variations. However, the field study showed that it is possible to reduce the kinetic degradation of the soil performance through an appropriate construction design., L’évolution actuelle du contexte socio-économique oblige les différents acteurs du secteur des travaux publics à s’adapter aux problématiques du développement durable. Dans le domaine des ouvrages en terre, les entreprises doivent proposer des solutions techniques de valorisation des matériaux situés dans l’emprise des projets afin de limiter les emprunts extérieurs et la mise en dépôt des sols non utilisés. Les techniques de traitement des sols constituent une possibilité de valoriser ces matériaux. Cette étude se focalise sur quatre produits de traitement (kaolinite, bentonite, chaux et ciment) ainsi que sur leurs effets sur le comportement hydromécanique et la résistance à l’érosion interne d’un limon fin. Un dispositif d’érosion interne HET optimisé a été spécialement conçu dans le cadre de ce travail pour déterminer les caractéristiques d’érosion des sols traités notamment à la chaux et au ciment. Au-delà de la caractérisation des effets du traitement sur le comportement hydromécanique à court terme des sols, le travail de cette étude s’étend sur la durabilité des traitements et l’évolution du comportent hydromécanique à long terme des sols traités soumis à des sollicitations hydriques de séchage-humidification. Cette étude de la durabilité est effectuée à travers une approche multi-échelle fondée sur des données d’études en laboratoire sur des éprouvettes de sol traité et des investigations in situ sur des ouvrages hydrauliques en sol traité. L’étude montre que les sollicitations hydriques peuvent dégrader les performances du sol traité. Cette dégradation se traduit par une augmentation de la conductivité hydraulique, une diminution de la résistance mécanique et aussi une diminution de la contrainte critique qui exprime une diminution de la résistance à l’érosion interne. La cinétique de perte de performances est conditionnée par la nature du produit de traitement et son dosage utilisé, et aussi par le niveau d’exposition, le nombre et l’intensité des cycles hydriques. Cependant, l’étude in situ montre qu’il est possible de diminuer la cinétique de dégradation des performances de sols à travers des dispositions constructives adaptées.

The current evolution of the social and economic context requires from the different actors of the public works sector to adapt their practices to the challenges of sustainable development. In the field of earthworks, companies must offer technical solutions to reuse the materials located in the vicinity of the projects in order to limit the borrowing materials and unused soils deposits. Soil treatment may allow the reuse of these materials. This study was focused on four treatment products (kaolinite, bentonite, lime and cement) as well as their effects on the hydro-mechanical behavior and internal erosion resistance of a fine silt. An enhanced HET device was designed in the framework of this study to determine the internal erosion characteristics of treated soils especially with lime and cement. Beyond the characterization of treatment effects on short-term hydro-mechanical behavior of soils, the work of this study extends to the durability of treatment and the evolution of long-term hydro-mechanical behavior of treated soils subjected to drying-wetting cycles. This study of sustainability was carried out through a multi-scale approach based on laboratory study data on soil samples and field investigations on hydraulic structures made of treated soil. The study showed that hydraulic conditions variations can decrease the performance of treated soils. These degradations result on an increase in hydraulic conductivity, a decrease of the mechanical strength and also a reduction in the critical shear stress which expresses a decrease of the internal erosion resistance. The kinetic of performance loss depends to the nature of the treatment product and percentage used as well as the exposure level, the number and amplitude of the hydraulic variations. However, the field study showed that it is possible to reduce the kinetic degradation of the soil performance through an appropriate construction design.; L’évolution actuelle du contexte socio-économique oblige les différents acteurs du secteur des travaux publics à s’adapter aux problématiques du développement durable. Dans le domaine des ouvrages en terre, les entreprises doivent proposer des solutions techniques de valorisation des matériaux situés dans l’emprise des projets afin de limiter les emprunts extérieurs et la mise en dépôt des sols non utilisés. Les techniques de traitement des sols constituent une possibilité de valoriser ces matériaux. Cette étude se focalise sur quatre produits de traitement (kaolinite, bentonite, chaux et ciment) ainsi que sur leurs effets sur le comportement hydromécanique et la résistance à l’érosion interne d’un limon fin. Un dispositif d’érosion interne HET optimisé a été spécialement conçu dans le cadre de ce travail pour déterminer les caractéristiques d’érosion des sols traités notamment à la chaux et au ciment. Au-delà de la caractérisation des effets du traitement sur le comportement hydromécanique à court terme des sols, le travail de cette étude s’étend sur la durabilité des traitements et l’évolution du comportent hydromécanique à long terme des sols traités soumis à des sollicitations hydriques de séchage-humidification. Cette étude de la durabilité est effectuée à travers une approche multi-échelle fondée sur des données d’études en laboratoire sur des éprouvettes de sol traité et des investigations in situ sur des ouvrages hydrauliques en sol traité. L’étude montre que les sollicitations hydriques peuvent dégrader les performances du sol traité. Cette dégradation se traduit par une augmentation de la conductivité hydraulique, une diminution de la résistance mécanique et aussi une diminution de la contrainte critique qui exprime une diminution de la résistance à l’érosion interne. La cinétique de perte de performances est conditionnée par la nature du produit de traitement et son dosage utilisé, et aussi par le niveau d’exposition, le nombre et l’intensité des cycles hydriques. Cependant, l’étude in situ montre qu’il est possible de diminuer la cinétique de dégradation des performances de sols à travers des dispositions constructives adaptées.

The stability of open-cast operations (mines or quarries), excavated in rock mass depends on their geometry, the geometrical characteristics of the rock mass fractures (orientation and spacing) and the mechanical characteristics of the rock mass and the discontinuities. The assessment of the rock instability risk depends on the quantity and quality of the available information on the rock mass and the analysis methods used for the evaluation of the mechanical behavior of the rock block collection. Different sites of natural rock slopes, of road cuts, and of open-cast mines and quarries, described in the present document, show that the known information can vary a lot from one site to the other. The steps of data gathering and analysis, then the modelling step used to analyze the stability of the rock slope, lead to uncertainties. The present work deals in particular with the influence of the geometrical uncertainties, and the different clustering methods, to define families of fractures, are examined. The combination of works suggested by several authors leads to a new method called PSMY. This method, together with the spectral method, was coded in the Mathematica platform, and the obtained results are compared with "hand-made" clusters, done with the DIPS software. The clustering methods are presented according to the ratio of classified orientations. The orientation and spacing of fracture families are fitted by statistical law. The statistical parameters are compared according to the clustering method used. These parameters have an influence on the construction of the rock mass geometrical model, called DFN (Discrete Fracture Network). The influence of the clusters on the stability analysis is studied from the stochastic models based on the limit equilibrium analysis, in the SWEDGE and RESOBLOK software. These programs do not take into account the same uncertainties, and as a result, they give different stability indicators. A sensitivity analysis of these indicators (number of unstable blocks, average volume of unstable blocks, and total volume of unstable blocks) versus the clustering methods used, and the orientation of the slopes, is carried out. A variance analysis allows an evaluation of the influence of these factors. The assessment of the global stability condition of the rock mass, depending on the cohesion and the friction angle, is proposed. The influence of the modeling method is analyzed by comparing 3D calculations using the limit equilibrium calculations, and 2D and 3D models using discrete rigid and deformable blocks. A coupling between RESOBLOK (limit equilibrium) and LMGC90 (discrete elements) allows the comparison of results on the same original geometry. For the case of Ax-les-Thermes road-cut, and for various unstable geometries, the influence of the model parameters is tested. Several cases are compared. The 3D simulation of an excavation, at different steps, is performed, and the mobilization index is studied, in order to compare several types of contacts, within the LMGC model, in relation to the possible sliding of blocks, at those different steps; La stabilité des exploitations à ciel ouvert (mines ou carrières) excavées dans des massifs rocheux dépend de leur géométrie, des caractéristiques géométriques de la fracturation (orientation et espacement) du massif et des caractéristiques mécaniques du massif et des discontinuités. L'évaluation du risque d'instabilité rocheuse dépend de la quantité et la qualité des informations disponibles sur le massif rocheux et des méthodes d'analyse mises en œuvre pour évaluer le comportement mécanique de l'assemblage de blocs rocheux. Différents sites de talus naturels, d'excavation routière et de carrières et mines à ciel ouvert présentés dans le mémoire montrent que les informations connues peuvent beaucoup varier d'un site à l'autre. Les étapes de collecte et de synthèse des données, puis de modélisation afin d'analyser la stabilité d'une pente rocheuse conduisent à des incertitudes. Le mémoire traite plus particulièrement de l'influence des incertitudes géométriques et différentes méthodes de regroupement des fractures en familles sont explorées. La combinaison du travail proposé par différents auteurs conduit à une nouvelle méthode nommée PSMY. Cette méthode ainsi que la méthode spectrale ont été programmées sous l'environnement Mathematica et les résultats sont comparés à des regroupements "manuels" réalisé avec le logiciel DIPS. Les méthodes de regroupement sont présentées selon le pourcentage d'orientations classé par la méthode. L’orientation et l’espacement des familles de fractures sont ajustées à des lois statistiques. Les paramètres statistiques sont comparés en fonction de la méthode regroupement. Ces paramètres influencent la construction d'un modèle géométrique du massif rocheux appelé DFN (Discrete Fracture Network). L'influence des regroupements sur l'analyse de stabilité est étudiée à partir de modélisation stochastiques à l'équilibre limite utilisant les logiciels SWEDGE et RESOBLOK. Ces logiciels ne prennent pas en compte les mêmes incertitudes et en sortie ils fournissent des indicateurs de stabilité différents. Une analyse de sensibilité des indicateurs (nombre de blocs instables, volume moyen du bloc instable, volume total instable) aux méthodes de regroupement et aux orientations de talus est réalisée. Une analyse de variance permet de préciser l'influence de ces deux facteurs. Une évaluation de l’état de stabilité global du massif, en fonction de la cohésion et de l'angle de frottement est proposée. L'influence de la méthode de modélisation est évoquée en comparant des calculs tridimensionnels à l'équilibre limite et des modélisations bidimensionnel et tridimensionnels en éléments discrets rigides ou déformables. Un couplage entre le logiciel RESOBLOK (équilibre limite) et LMGC90 (éléments discrets) permet de comparer les analyses pour une même géométrie. Pour le cas de la déviation d'Ax-les-Thermes et certaines configurations géométriques instables l’influence des paramètres de modèles sur stabilité du talus est testée. Plusieurs variantes sont comparées. La simulation d’une excavation en 3D multi phases est réalisée et l’indice de mobilisation est étudié pour comparer les divers types de contacts dans LMGC90 par rapport à l’éventuel glissement des blocs dans les étapes successives

This thesis systematically investigates the breakup and coalescence processes of the involved droplet (bubble) interface under magnetic field or not in two-phase microfluidic flow, by using a high speed digital camera. The whole breakup processes of ferrofluid interface under different magnetic fields were investigated and compared. It was found that the morphological structure and necking velocity of the interface can be actively controlled by the magnetic force. Thus the volumes and the formation frequencies of ferrofluid droplets can be actively adjusted. The breakup of Liquid-Liquid interface usually leads to the formation of satellite droplet with its size proportional to the capillary number of the continuous phase. The coalescences of droplets were investigated. The evolution of the neck connecting two droplets was analyzed. It was found that the formation of liquid bridge or neck could occurs with a visible gap in the order of tens of micrometers between the leading edges under magnetic field and the inertia of the ferrofluid originating from the magnetic attraction fields becomes the driving force at the initial stage of coalescence instead of capillary force; Cette thèse a étudié de façon systématique les processus de la rupture et de la coalescence impliquant une interface liquide-liquide (gouttelettes) ou gaz-liquide (bulles) en présence ou pas d’un champ magnétique dans des dispositifs microfluidiques, à l’aide d’une caméra rapide. Les mécanismes de la rupture d’une interface ferrofluide sous différents champs magnétiques ont été étudiés et comparés. On a constaté que la structure morphologique et la vitesse d’amincissement du cou peuvent être contrôlées activement par la force magnétique. Ainsi, le volume et la fréquence de formation de gouttelettes de ferrofluide peuvent être pilotés aisément. La rupture de l'interface liquide-liquide a généralement conduit à la formation de gouttelettes satellites avec sa taille proportionnelle au nombre capillaire de la phase continue. La coalescence des gouttelettes ont aussi été étudiée avec l’analyse de l'évolution du cou reliant deux gouttelettes. Il a été constaté que la formation du pont ou du cou liquide pouvait se produire dans la gamme de l'ordre de dizaines de micromètres entre les bords d'attaque sous champ magnétique. L'inertie d’origine d'attraction magnétique sur des gouttelettes de ferrofluide devenait la force motrice pour la coalescence lors de la première étape au détriment de la force capillaire

This thesis systematically investigates the breakup and coalescence processes of the involved droplet (bubble) interface under magnetic field or not in two-phase microfluidic flow, by using a high speed digital camera. The whole breakup processes of ferrofluid interface under different magnetic fields were investigated and compared. It was found that the morphological structure and necking velocity of the interface can be actively controlled by the magnetic force. Thus the volumes and the formation frequencies of ferrofluid droplets can be actively adjusted. The breakup of Liquid-Liquid interface usually leads to the formation of satellite droplet with its size proportional to the capillary number of the continuous phase. The coalescences of droplets were investigated. The evolution of the neck connecting two droplets was analyzed. It was found that the formation of liquid bridge or neck could occurs with a visible gap in the order of tens of micrometers between the leading edges under magnetic field and the inertia of the ferrofluid originating from the magnetic attraction fields becomes the driving force at the initial stage of coalescence instead of capillary force, Cette thèse a étudié de façon systématique les processus de la rupture et de la coalescence impliquant une interface liquide-liquide (gouttelettes) ou gaz-liquide (bulles) en présence ou pas d’un champ magnétique dans des dispositifs microfluidiques, à l’aide d’une caméra rapide. Les mécanismes de la rupture d’une interface ferrofluide sous différents champs magnétiques ont été étudiés et comparés. On a constaté que la structure morphologique et la vitesse d’amincissement du cou peuvent être contrôlées activement par la force magnétique. Ainsi, le volume et la fréquence de formation de gouttelettes de ferrofluide peuvent être pilotés aisément. La rupture de l'interface liquide-liquide a généralement conduit à la formation de gouttelettes satellites avec sa taille proportionnelle au nombre capillaire de la phase continue. La coalescence des gouttelettes ont aussi été étudiée avec l’analyse de l'évolution du cou reliant deux gouttelettes. Il a été constaté que la formation du pont ou du cou liquide pouvait se produire dans la gamme de l'ordre de dizaines de micromètres entre les bords d'attaque sous champ magnétique. L'inertie d’origine d'attraction magnétique sur des gouttelettes de ferrofluide devenait la force motrice pour la coalescence lors de la première étape au détriment de la force capillaire

The effects of lime treatment on the mechanical properties of soils are usually not accounted for in the design of geotechnical structures. As a result the potential of lime treatment has not been fully exploited. In this thesis, a comprehensive experimental program has been carried out to identity the key features of the mechanical behaviour of structured materials. The chemical modifications arising from lime treatment were quantified using thermal analysis methods. From these results a non-linear chemo-mechanical coupling was established between the concentration of cementitious compounds and the yield stress. Using these results, a new formulation to model the degradation of the structure at yield has been developed and implemented in a constitutive model for structured materials. This new model, developed in the framework of the Modified Cam Clay model, requires a limited number of additional parameters that all have a physical meaning and can all be determined from a single isotropic compression test. The model has proven to be successful in reproducing the key features of structured materials and for the modelling of the mechanical behaviour of lime treated specimens under various stress paths. Due to similarities in behaviour, it is shown that the formulation is also suitable for naturally structured soils. To account for a structured material in the design of geotechnical structures, a fully functional finite element program for elasto-plastic problems was developed including the pre- and post-processing of the results. A thorough validation has confirmed the good implementation of the finite element method and its suitability for the modelling of complex geometries involving structured materials, Le traitement des sols à la chaux est une méthode couramment utilisée pour améliorer les propriétés mécaniques de sols aux performances insuffisantes. Cependant, ces améliorations mécaniques ne sont pas prises en compte dans les calculs de dimensionnement. Cette thèse propose une méthodologie pour pallier à ce problème. Un programme expérimental approfondi est réalisé afin de décrire avec précision le comportement mécanique d’un sol traité et les processus associés à la structuration introduite par le traitement. La composition chimique, et notamment la quantité de composés cimentaires hydratés, est déterminée par analyses thermogravimétriques et thermodifférentielles. Un couplage non-linéaire entre la quantité d’hydrates et la limite élastique se basant sur ces résultats expérimentaux est proposé. À partir de ces résultats, une nouvelle formulation a été développée afin de modéliser la dégradation de la structure en plasticité, et a servi au développement d’une nouvelle loi de comportement élasto-plastique basée sur le modèle de Cam Clay Modifié. Cette dernière a montré reproduire correctement les principaux traits de comportement spécifiques aux sols traités. Il est démontré que ce modèle est également adapté aux sols naturellement structurés. Afin de prendre en compte les effets du traitement et la présence de structure dans le dimensionnement, un programme basé sur la méthode des éléments finis, comprenant le pré- et le post-processing de la géométrie et des résultats du problème, a été développé. Une validation rigoureuse a confirmé l'implémentation correcte de la méthode et son potentiel pour l’optimisation du dimensionnement des ouvrages

Pleated filters are widely used in air treatments because of the advantageous effective surface to overall dimension ratio they offer. Their major drawback though resides in their reduced lifetime which still needs to be controlled. Indeed, when clogging, the pressure drop considerably increases, the filtration flow is then no longer maintained which might lead to the deterioration of the media. It is then crucial to characterize the evolution of the pressure drop under operating conditions in order to best design these equipments. Part of our work consisted in studying how the operating conditions influence the geometry of the deposit. To do so, we used Single-Photon Emission Computed Tomography (SPECT), a non-destructive imaging technique that keeps intact the particle structuring. The visualization of aerosol deposit at the beginning of the filtration process allows observing preferential particle deposition on the whole height of the pleat. A numerical approach was used to study the permeability of bimodal fibrous media and we experimentally studied the local velocity as well as the biphasic flow inside pleated filter media. Comparison between experiments and simulations allowed us to validate the Geodict code for a wide range of media properties and velocities. Regarding bimodal fibrous media, the fast data acquisition has allowed testing several existing models which resulted in classifying them in a unique way. If the experimental results on the initial deposition in pleated filters are encouraging, those related to beforehand clogging point to several improvements regarding the technique we used; L'utilisation de filtres plissés dans le domaine du traitement de l'air est extrêmement fréquente en raison de leur simplicité d'utilisation et de maintenance. Néanmoins, au cours du colmatage, la perte de charge augmente considérablement avec le dépôt de particules. Par conséquent le débit de filtration n'est plus maintenu et une détérioration des médias peut se produire. Il est essentiel de caractériser cette évolution en fonction des conditions opératoires pour mieux concevoir ces équipements. Les travaux ont consisté en une étude expérimentale de l'influence des conditions opératoires sur la géométrie du dépôt par des observations en Tomographie d'Émission MonoPhotonique. Une approche numérique a consisté, dans un premier temps, en une étude de la perméabilité de médias bien caractérisés afin de valider le code GeoDict. Dans un deuxième temps, une étude de la vitesse de l'air a été entreprise sur les médias plissés. La confrontation des résultats numériques avec les données expérimentales nous a permis une validation du code dans les gammes de propriétés des médias et de vitesses typiques. Concernant les médias bimodaux, l'acquisition rapide de données a permis de tester différentes corrélations analytiques afin de classifier les modèles existants d'une manière qui facilite leur utilisation. Si les résultats expérimentaux obtenus sur la localisation du dépôt préférentiel initial dans les filtres plissés sont encourageants, notamment vis-à-vis des simulations d'écoulement réalisées avec GeoDict, ceux issus des acquisitions sur les filtres précolmatés sont à prendre avec prudence et suggèrent d'apporter un certain nombre d'amélioration à notre technique

Pleated filters are widely used in air treatments because of the advantageous effective surface to overall dimension ratio they offer. Their major drawback though resides in their reduced lifetime which still needs to be controlled. Indeed, when clogging, the pressure drop considerably increases, the filtration flow is then no longer maintained which might lead to the deterioration of the media. It is then crucial to characterize the evolution of the pressure drop under operating conditions in order to best design these equipments. Part of our work consisted in studying how the operating conditions influence the geometry of the deposit. To do so, we used Single-Photon Emission Computed Tomography (SPECT), a non-destructive imaging technique that keeps intact the particle structuring. The visualization of aerosol deposit at the beginning of the filtration process allows observing preferential particle deposition on the whole height of the pleat. A numerical approach was used to study the permeability of bimodal fibrous media and we experimentally studied the local velocity as well as the biphasic flow inside pleated filter media. Comparison between experiments and simulations allowed us to validate the Geodict code for a wide range of media properties and velocities. Regarding bimodal fibrous media, the fast data acquisition has allowed testing several existing models which resulted in classifying them in a unique way. If the experimental results on the initial deposition in pleated filters are encouraging, those related to beforehand clogging point to several improvements regarding the technique we used, L'utilisation de filtres plissés dans le domaine du traitement de l'air est extrêmement fréquente en raison de leur simplicité d'utilisation et de maintenance. Néanmoins, au cours du colmatage, la perte de charge augmente considérablement avec le dépôt de particules. Par conséquent le débit de filtration n'est plus maintenu et une détérioration des médias peut se produire. Il est essentiel de caractériser cette évolution en fonction des conditions opératoires pour mieux concevoir ces équipements. Les travaux ont consisté en une étude expérimentale de l'influence des conditions opératoires sur la géométrie du dépôt par des observations en Tomographie d'Émission MonoPhotonique. Une approche numérique a consisté, dans un premier temps, en une étude de la perméabilité de médias bien caractérisés afin de valider le code GeoDict. Dans un deuxième temps, une étude de la vitesse de l'air a été entreprise sur les médias plissés. La confrontation des résultats numériques avec les données expérimentales nous a permis une validation du code dans les gammes de propriétés des médias et de vitesses typiques. Concernant les médias bimodaux, l'acquisition rapide de données a permis de tester différentes corrélations analytiques afin de classifier les modèles existants d'une manière qui facilite leur utilisation. Si les résultats expérimentaux obtenus sur la localisation du dépôt préférentiel initial dans les filtres plissés sont encourageants, notamment vis-à-vis des simulations d'écoulement réalisées avec GeoDict, ceux issus des acquisitions sur les filtres précolmatés sont à prendre avec prudence et suggèrent d'apporter un certain nombre d'amélioration à notre technique

Industrial photocatalysis applications requires further research to control and design photocatalytic reactors. The micro structuring of reactors is an innovative solution. It increases the catalytic efficiency by improving the contact of the pollutant with the catalyst. Thesis aim concerned the destruction of a model pollutant, salicylic acid, with deposited catalyst (TiO2) at the bottom of a microchannel. The research is based on experiments and simulations to predict the photocatalytic degradation for a given reactor configuration. In our experimental conditions of low flow, an axial dispersion and a concentration gradient between the catalyst surface and the reaction medium have been demonstrated. We have shown that the gradient of velocity is more pronounced in low cross section microreactors. This gradient allows the formation of recirculation zones near the catalytic surface thus creating favorable conditions for contact between the pollutant and the catalytic surface. Using numerical soft, we have simulated the reaction of photocatalytic degradation to represent our experimental results. On a one hand, the role of diffusion to the catalytic surface has been demonstrated. On the other hand, numerical simulation has allowed to establish a relationship between the Sherwood number as a function of Reynolds number, Schmidt and the hydraulic diameter of the microchannel. This relationship is used to estimate the mass transfer constant to represent correctly the photocatalytic degradation, L'implantation dans l'industrie de la photocatalyse nécessite de nouvelles recherches afin de contrôler et dimensionner les réacteurs photocatalytiques. Une solution innovante concerne la microstructuration des réacteurs. En effet, elle permet d'augmenter l'efficacité catalytique en améliorant le contact du polluant avec le catalyseur. L'objectif de la thèse concernait la destruction d'un polluant modèle, l'acide salicylique, par un catalyseur (TiO2) déposé au fond d'un microcanal. La recherche repose sur des expérimentations mais aussi sur des simulations permettant de prévoir la dégradation photocatalytique pour une configuration de réacteur donnée. Dans nos conditions expérimentales de faible débit, nous avons mis en évidence une dispersion axiale et un gradient de concentration entre la surface catalytique et le milieu réactionnel. En tenant compte des conditions d'écoulement, nous avons montré que le gradient de vitesse est plus marqué dans les microréacteurs de faible section de passage. Ce gradient favorise la formation de zones de recirculation près de la surface catalytique créant ainsi des conditions favorables au contact entre le polluant et la surface catalytique. A l'aide d'un logiciel numérique, nous avons simulé la réaction de dégradation photocatalytique qui permet de rendre compte de nos résultats expérimentaux. D'une part, le rôle primordial de la diffusion vers la surface catalytique a été mis en évidence. D'autre part, la simulation numérique nous a permis d'établir la relation du nombre de Sherwood en fonction des nombres de Reynolds, de Schmidt et du diamètre hydraulique du microcanal. Cette relation permet d'estimer la constante de transfert de matière pour représenter correctement la dégradation quel que soit le type de microréacteur envisagé

Restoration is a structural geology tool which usefulness has been shown in various fields, such as the understanding of deformation mechanisms, fracture forecasting, or structural models validation. In 3D, the geomechanical approach is particularly promising as it allows the deformation at any point in the model to be infered, but also to account for material property contrasts. However, it requires to generate a tetrahedral mesh, a step that may be, so far, extremely difficult, particularly when the geometry of the model is complex (interconnected faults, faults with small throws, unconformities or layer pinch-out).Implicit modeling consists in representing stratigraphic horizons by property isovalues, allowing most meshing issues to be overcame. In a first part, a geomechanical restoration method suitable for implicit models is presented. The second part of this work focuses on the application of this method to complex cases. New behavioral laws and boundary conditions have been developed, and sensitivity studies to these parameters have been performed. This work is applied to a petroleum basin located in the Niger delta toe where the restoration allows the structural evolution of the system to be constrained. The last part presents the integration of a 3D isostatic decompaction method into the 3D restoration process. Applied to the Annot Sandstones syncline (SE France), this method allows to characterize the deformation and burial history of the basin, La restauration est un outil de géologie structurale qui a prouvé son utilité dans divers domaines, tels que la compréhension des mécanismes de déformation, la prédiction de fractures, la validation de modèles structuraux, etc. En 3D, l'approche géomécanique est particulièrement prometteuse, puisqu'elle permet non seulement de connaître la déformation en tout point du modèle, mais aussi de prendre en compte les contrastes de propriété des matériaux. Cependant, elle nécessite un maillage tétraédrique du modèle dont la génération est à ce jour très difficile dans le cas de modèles à géométrie complexe (réseaux de failles interconnectés ou à faible rejet, discordances ou amincissements de couche). La modélisation implicite consiste à représenter les horizons stratigraphiques par des isovaleurs de propriété, ce qui permet de s'affranchir de la plupart des contraintes de maillage. Dans une première partie, une méthode de restauration géomécanique 3D applicable à des modèles implicites est présentée. La deuxième partie de ce travail s'intéresse à l'application de celle-ci à des cas complexes. Dans ce cadre, de nouvelles lois de comportement géomécaniques et conditions aux limites ont été développées, et des études de sensibilité à ces paramètres ont été menées. Ces travaux sont illustrés sur un bassin pétrolier complexe situé dans le delta du Niger où la restauration permet de définir une chronologie de déformation. Enfin, la troisième partie présente l'intégration d'une méthode de décompaction isostatique en volume au sein du processus de restauration 3D. Appliquée au synclinal des grès d'Annot (SE de la France), cette méthode permet de caractériser l'histoire de déformation et d'enfouissement du bassin

Gas-liquid two-phase flow is an important research project in microfluidics. The gas-liquid two-phase flow, the bubble formation and moving behaviours in microchannels were investigated, by using a high speed digital camera and a micro Particle Image Velocimetry (micro-PIV). The gas-liquid two-phase flow in vertical rectangular microchannels was investigated and a flow pattern map was constructed; the bubble formation in both Newtonian and non-Newtonian fluids in cross-flowing microfluidic T-junctions and flow-focusing devices was investigated; the bubble formation process could be divided into expansion, collapse and pinch-off stages; the collapse speed of the gaseous thread in the second stage is controlled by the squeezing pressure, and is proportional to the liquid flow rates; while the minimum width of the neck of the gaseous thread in the third stage for bubble formation in flow-focusing devices could be scaled with the remaining time to the ultimate pinch-off as a power law relationship with an exponent of 1/3; the PAAm solutions prolong the gaseous thread in the tangential direction of the neck; bubble coalescence in a microchannel with an expansion section was studied; the bubble behavior in a microchannel with a loop was also investigated, Les écoulements gaz-liquide constituent un axe de recherche très actif en microfluidique. Le rapport des débits entre les deux phases, la formation de bulles et les champs de vitesse des microcanaux ont été étudiés dans cette thèse, en utilisant une caméra numérique rapide et un microsystème de Particule Image Velocimetry (micro-PIV). En particulier, le diagramme des phases gaz-liquide ont été établi dans des microcanaux carrés ; la formation des bulles en fluides tant newtoniens que non newtoniens a été étudiée en détail dans plusieurs configurations géométriques telles que T-injonction et flow-focusing. Les mécanismes régissant la formation d'une bulle ont été modélisés pour chaque étape : expansion, amincissement et rupture. L'étape amincissement de la traînée d'une bulle est notamment contrôlée par une pression orthogonale qui dépend du débit du liquide. Dans le cadre de flow-focusing, le mécanisme de la rupture du film gazeux peut être décrit par une loi d'échelle reliant l'épaisseur minimale du film au temps restant juste avant la rupture avec un exposant 1/3. Le caractère non newtonien de fluides PAAm allonge la traînée d'une bulle par rapport aux fluides newtoniens. Enfin, l'étude de la coalescence entre bulles a été entreprise à l'échelle microscopique ainsi que le comportement complexe des trains de bulles dans des réseaux de microcannaux

Gas-liquid two-phase flow is an important research project in microfluidics. The gas-liquid two-phase flow, the bubble formation and moving behaviours in microchannels were investigated, by using a high speed digital camera and a micro Particle Image Velocimetry (micro-PIV). The gas-liquid two-phase flow in vertical rectangular microchannels was investigated and a flow pattern map was constructed; the bubble formation in both Newtonian and non-Newtonian fluids in cross-flowing microfluidic T-junctions and flow-focusing devices was investigated; the bubble formation process could be divided into expansion, collapse and pinch-off stages; the collapse speed of the gaseous thread in the second stage is controlled by the squeezing pressure, and is proportional to the liquid flow rates; while the minimum width of the neck of the gaseous thread in the third stage for bubble formation in flow-focusing devices could be scaled with the remaining time to the ultimate pinch-off as a power law relationship with an exponent of 1/3; the PAAm solutions prolong the gaseous thread in the tangential direction of the neck; bubble coalescence in a microchannel with an expansion section was studied; the bubble behavior in a microchannel with a loop was also investigated; Les écoulements gaz-liquide constituent un axe de recherche très actif en microfluidique. Le rapport des débits entre les deux phases, la formation de bulles et les champs de vitesse des microcanaux ont été étudiés dans cette thèse, en utilisant une caméra numérique rapide et un microsystème de Particule Image Velocimetry (micro-PIV). En particulier, le diagramme des phases gaz-liquide ont été établi dans des microcanaux carrés ; la formation des bulles en fluides tant newtoniens que non newtoniens a été étudiée en détail dans plusieurs configurations géométriques telles que T-injonction et flow-focusing. Les mécanismes régissant la formation d'une bulle ont été modélisés pour chaque étape : expansion, amincissement et rupture. L'étape amincissement de la traînée d'une bulle est notamment contrôlée par une pression orthogonale qui dépend du débit du liquide. Dans le cadre de flow-focusing, le mécanisme de la rupture du film gazeux peut être décrit par une loi d'échelle reliant l'épaisseur minimale du film au temps restant juste avant la rupture avec un exposant 1/3. Le caractère non newtonien de fluides PAAm allonge la traînée d'une bulle par rapport aux fluides newtoniens. Enfin, l'étude de la coalescence entre bulles a été entreprise à l'échelle microscopique ainsi que le comportement complexe des trains de bulles dans des réseaux de microcannaux

The numerical work described in this thesis deals with model reduction in the field of heat transfer and fluid mechanics. The reduced models being solved rapidly, they can be coupled with some optimal algorithms. The construction of reduced models is performed based on the Modal Identification Method (MIM) that has been developed in our laboratory for many years. This method is based on the state representation under a modal form. The reduced model is identified through the minimization of a function containing the parameters of the model and based on the difference between the outputs of a detailed model and those of a reduced one. An application is proposed dealing on a flow over a backward-facing step, with a time-varying heat flux density applied upstream of the step. Reduction modelling leads to important reduction of simulation times, still while maintaining a very good accuracy. We then propose an approach combining the theory of optimal control with the model obtained by the MIM. This approach is used to find, on-line, the fluxes applied upstream the step such that the obtained temperature stays as close as possible to the set point., Les travaux numériques décrits dans ce mémoire portent sur la réduction de modèle dans le domaine de la thermique et de la mécanique des fluides. Ces modèles réduits étant rapides en temps d'exécution, ils permettent la mise en place d'algorithmes de contrôle optimal. La construction des modèles réduits est réalisée à partir de la Méthode d'Identification Modale (MIM) développée au laboratoire depuis de nombreuses années. Cette méthode s'appuie sur la représentation d'état sous forme modale : le modèle réduit est identifié à travers la minimisation d'une fonctionnelle comprenant les paramètres du modèle et basée sur l'écart entre les réponses du modèle détaillé et celles du modèle réduit. Une application est proposée sur un écoulement le long d'une marche descendante, avec un flux de chauffage intervenant en amont de la marche. On montre comment on peut réduire de façon importante les temps de simulation, tout en gardant une très bonne précision. Sur cette application, nous proposons également une approche combinant la théorie de la commande optimale avec les modèles réduits obtenus par la MIM. Cette approche est utilisée pour trouver une loi de commande en flux permettant de maintenir un profil de température proche d'une consigne dans l'écoulement en aval de la marche.

The steel quality is an important focus for ArcelorMittal. The presence of inclusions in the form of precipitates of different composition can produce defects in the slab (exfoliation) which can affect rolling operations and giving a poor surface sheet of steel. To provide actuators for efficient removal of these inclusions, we must first know the phenomena that govern their birth and their lives. The objective of this thesis is to understand and to predict the stages of germination / growth of solid inclusions present in liquid steel. In this context, a numerical tool has been implemented. It takes into account not only the mechanics of fluid through the use of Fluent® but also thermodynamics with a homemade tool called MIPPHASOLACIDO and an inclusions population size tracking with a new method based on the work of Marchisio et al. (2004). Multi-QMOM is an elegant tool built with a finite number of QMOM (Quadrature Method Of Moments) used to monitor changes of population size. It represents the major innovation of this thesis. The numerical tool whose results were compared to analytical and industrial cases can model the phenomena of precipitation, growth, aggregation and extraction of various solid inclusions. It allows the study of different parameters and their impacts on inclusions population.; La qualité des aciers est un axe important pour ArcelorMittal. La présence d'inclusions, sous forme de précipités de compositions diverses, peut produire des défauts dans les brames (exfoliations) qui pourront se répercuter lors des opérations de laminage, en donnant un mauvais état de surface à la tôle d'acier. Afin de proposer des actionneurs permettant une élimination efficace de ces inclusions, nous devons avant tout connaître les phénomènes qui régissent leur naissance et leur vie. L'objectif de cette thèse est donc de comprendre et de pouvoir prédire les étapes de germination/croissance des inclusions solides présentes dans l'acier liquide. Dans cette optique, un outil numérique a été mis en place. Il prend en compte non seulement la mécanique des fluides à travers l'utilisation de Fluent® mais aussi la thermodynamique avec un outil fait maison et le suivi en taille des populations d'inclusions avec une nouvelle méthode basée sur les travaux de Marchisio et al. (2004). Le module Multi-QMOM est un outil élégant, construit avec un nombre fini de méthode QMOM (Quadrature Method Of Moments), utilisé pour suivre l'évolution en taille d'une population multiple. Il représente l'innovation majeure de cette thèse. L'outil numérique dont les résultats ont été confrontés à des cas analytiques et industriels permet de modéliser les phénomènes de précipitation, de croissance, d'agrégation et d'extraction de diverses inclusions solides et il permet l'étude de différents paramètres et leurs impacts sur la population inclusionnaire.

The two-phase flows through fractures occur in many industrial and environmental processes.Modelling these flows needs understanding the role of geometric characteristics of the fractures surfaces. The objective of this thesis is to study experimentally the two-phase flow and the dispersion of tracers in fractures which the aperture map is known and propose a suitable model.The inertial two-phase flow model is based on the generalization of the single-phase full cubiclaw that accounts for non-Darcian effects by using the F function approach. On this approach,the presence of a second fluid is taken into account through a multiplier function introduced into the superficial velocity of each fluid. This model let us predict the relative permeability and the saturation as a function of the two fluid flows with two adjustment parameters.The concentration profiles obtained during dispersion experiments were interpreted with different models to study their ability to describe the dispersion in fractures and hence the heterogeneity of them. The results show that the classical convection-dispersion is not appropriate due to the heterogeneity of the fractures. On the other hand, the stratified model,with only one estimated parameter, gives satisfactory results and allows quantifying heterogeneity. Finally, with five estimated parameters, from the concentration curve at the outlet of the fractures, the CTRW model can describe the concentration profiles of all abscissas with good accuracy, Les écoulements diphasiques en fracture se produisent dans de nombreuses applications industrielles et environnementales. Pour modéliser ces écoulements, il est important de connaître le rôle des caractéristiques géométriques des surfaces constituant les fractures.L'objectif de ce travail de thèse est d'étudier expérimentalement les écoulements diphasiques et la dispersion de traceurs dans des fractures réelles dont on connaît la cartographie des surfaces et d'en proposer une modélisation adéquate.L'écoulement diphasique inertiel a été modélisé à l'aide de la généralisation de la loi de Darcy couplée à l'approche de la fonction F. Le modèle proposé permet de prédire les perméabilités relatives et les saturations en fonction des débits des deux fluides avec deux paramètres d'ajustement.Les profils de concentration obtenus lors des expériences de dispersion ont été interprétés avec différents modèles afin d'évaluer leur capacité à décrire la dispersion dans les fractures et donc l'hétérogénéité de celles-ci. Les résultats montrent que le modèle stratifié, avec un seul paramètre estimé, donne des résultats satisfaisants et permet de quantifier l'hétérogénéité.Enfin, avec cinq paramètres estimés, à partir de la courbe de la concentration de traceur à la sortie des fractures, le modèle CTRW permet de décrire avec une bonne précision les profils de concentration de toutes les abscisses

Lattice structures are widely used, especially in civil engineering. The finite element analysis of such structures might require a consequent amount of computational time when the periodical mesh of this lattice is complex. Defining an equivalent continuous medium in order to reduce the size of the problem appears to be interesting. The aim of the first part of this document is to apply a homogenization method in order to find a Timoshenko beam model macroscopically equivalent to a slender structure which is periodical in the longitudinal direction. One of the unusual aspects tackled reside in the study of structures with periodical cells having a longitudinal asymmetry. In addition, the case of periodical structures with free deformation (e.g. thermal dilatation) at microscopic scale is dealt. The second part is consecrated to the study of the axonemal structure of the biological cell flagella and Cilia. A shorten version of the axonemal structure is studied at first and homogenized afterward, Les structures treillis constituées d’un nombre important de barres sont largement utilisées, notamment en génie civil. L’étude par éléments finis de telles structures se révèle très coûteuse dès que la maille répétitive du treillis est complexe. Il s’avère intéressant de réduire la taille du problème en définissant un milieu continu équivalent. L’objectif de la première partie de ce travail est de proposer, en se plaçant dans le cadre des méthodes d’homogénéisation des milieux périodiques, une poutre de Timoshenko équivalente à une structure périodique dont l’une des dimension est grande par rapport aux deux autres. Une des originalités réside dans l’étude de cellules de base non symétriques. Par ailleurs, on s’intéresse à la prise en compte de déformations libres (par exemple, d’origine thermique) apparaissant à l’échelle microscopique. La seconde partie est consacrée à l’étude de la structure axonémale du flagelle et des cils vibratiles. Il s’agit de proposer et valider un modèle pour cette structure biomécanique complexe et d’appliquer ensuite la méthode d’homogénéisation proposée

Computational Fluid Dynamics (CFD) techniques are used to model the precipitation of Barium Carbonate in a solid-liquid fluidized bed reactor. The study is divided in two sections: The hydrodynamic behavior and the barium carbonate precipitation. The CFD model is validated with experimental results for both cases. In the first part, a neutralization reaction in the fluidized bed column with and without solids is carried out. In order to simulate the liquid phase reaction, the Eddy Dissipation Model (EDM), the Eddy Dissipation - Multiple Time Scale turbulent mixer model (EDM-MTS) and the Eddy Dissipation Concept (EDC) reaction models are tested. The model EDM-MTS is chosen for giving the best approach and for being the less computationally expensive. In the second part, the population balance equation is added to the model in order to calculate the Particle Size Distribution (PSD) in the fluidized bed reactor. The discrete method is chosen to solve the population balance equation coupled with the multi-phase Eulerian-Eulerian model, the k-e turbulence model and the EDM-MST model. The nucleation, growth and aggregation kinetics of BaCO3 are included in the precipitation model. The experimental BaCO3 precipitations realized in a fluidized bed reactor allowed the CFD precipitation model validation. The PSD obtained by CFD are in good agreement with the experimental PSD, La mécanique des fluides numérique (CFD) est utilisée pour modéliser la précipitation du carbonate de baryum dans un réacteur à lit fluidisé. L'étude est divisée en deux parties : la modélisation de l'hydrodynamique du réacteur et la modélisation de la précipitation du carbonate de baryum. Pour ces deux parties, les modèles sont validés par l'expérience. Dans la première partie de l'étude, des réactions instantanées de neutralisation en absence et en présence de particules solides inertes sont mises en oeuvre dans le réacteur à lit fluidisé. Pour représenter la réaction chimique dans la phase liquide, plusieurs modèles son testés : le modèle Eddy Dissipation (EDM), le modèle Eddy Dissipation Concept (EDC) et le modèle modifié Eddy Dissipation- Multiple Time Scale turbulent mixer (EDM-MTS). Le modèle qui donne la meilleure prédiction de la réaction chimique est choisi : il s'agit du modèle EDM-MTS, qui demande aussi le moins de temps de calcul. Dans la deuxième partie, l?équation de bilan de population est incorporée au code de calcul pour prédire la distribution de taille des particules (PSD). La méthode des classes est implémentée couplée avec le modèle des écoulements multiphasiques Eulérien-Eulérien, le modèle de turbulence k-e et le modèle EDM-MTS. Les cinétiques de précipitation de BaCO3 sont incluses dans le modèle. Des expériences de précipitation du BaCO3 en réacteur à lit fluidisé permettent de valider le modèle de CFD développé. La PSD donnée par le modèle de CFD est en bonne concordance avec les résultats expérimentaux

VOCs are a source of discomfort for occupants and can have harmful effects on health. Photocatalysis is a promising process to remove these compounds. ArcelorMittal Liège Research has indeed developed TiOB2 Bthin film coated on steel which can be easily integrated into purification air systems. Objectives of this research are firstly, identifying parameters which influence the photocatalytic degradation of these VOCs using of such catalyst. Secondly, validate the possibility to predict performances of the photocatalytic reactor by combining fluid flow and photocatalytic reaction. An experimental set-up including an annular reactor has been designed and realized in order to conduct kinetic studies. Acetaldehyde has been chosen as a pollutant model and studied mainly with concentrations of few ppm ranges. Its conversion yield depends on pollutant concentration, relative humidity, temperature, UV light intensity and presence of others VOCs. Based on a theoretical surface model of the photocatalytic reaction, and using experimental data, this reaction is modelled. In order to predict performances of the photocatalytic reactor packed thin film of TiOB2B, a methodology which consists to combine fluid flow and photocatalytic reaction has been developed and validated. Space-time distribution of pollutant concentration is then determined with a computational fluid dynamics approach. Comparison between experimental data and numerical simulations obtained with a photocatalytic reactor at the laboratory scale and another one in semi-pilot scale are well satisfactory, Les COV présents en air intérieur sont une source d'inconfort pour les occupants et peuvent être nocifs pour la santé. La photocatalyse apparaît comme un procédé prometteur pour l'abattement de ces polluants. ArcelorMittal Liège Research développe à cet effet des films minces de TiOB2B déposé sur acier qui peuvent être intégrés dans des systèmes de purification d'air. Les objectifs de cette recherche sont, d'une part, identifier les paramètres influençant la dégradation photocatalytique de ces COV avec ce type de catalyseur et d'autre part, valider la possibilité de prédire les performances d'un réacteur photocatalytique par un modèle couplant l'écoulement et la réaction photocatalytique. Pour réaliser l'étude cinétique à des teneurs de quelques ppm, un dispositif expérimental adapté a été conçu ; il utilise un réacteur annulaire continu. L'acétaldéhyde a été choisi comme polluant modèle. La conversion dépend de la concentration en polluant, du taux d'humidité, de la température, de l'intensité lumineuse en UV et de la présence d'autres COV. La réaction est modélisée avec un modèle théorique surfacique de la réaction photocatalytique. Afin de prédire les performances d'un réacteur photocatalytique utilisant des films minces de TiOB2B, une méthodologie consistant à coupler l'écoulement et la réaction photocatalytique a été développée et validée. La distribution spatio-temporelle de la concentration en polluant est alors déterminée par une approche de mécanique des fluides numériques. La comparaison des résultats expérimentaux et des simulations obtenus avec un réacteur photocatalytique de laboratoire et un autre à l'échelle semi-pilote est très satisfaisante

Dans une optique d'amélioration ou de construction de moteurs fusée dérivés du modèle Vulcain, ce projet propose une nouvelle géométrie d'injecteur permettant l'atomisation efficace d'une nappe annulaire liquide en un brouillard de gouttelettes, nécessaire à un rendement élevé de combustion de propergols. Les images de tomographie et de visualisations rapides ont permis de comprendre les mécanismes de rupture de la nappe liquide par un courant interne et/ou externe tandis que la phase diluée a été caractérisée par les techniques de Vélocimétrie par Images de Particules et d'Interférométrie par Phase Doppler. Plusieurs modes de rupture ont été identifiés dépendant principalement du rapport du flux de quantité de mouvement gaz/liquide et une étude fréquentielle a été réalisée sur les battements de la nappe. Le mode qui nous intéresse plus particulièrement, nommé « arbre de Noël », permet une atomisation primaire beaucoup plus efficace que celui d'un jet liquide assisté par une couronne annulaire de gaz mais aussi d'obtenir une granulométrie plus faible des gouttes produites. L'étude des différents paramètres : rapport des densités de flux de quantité de mouvement gaz/liquide, pression ambiante et rotation du gaz a montré que les deux premiers cités ont une influence très importante sur la granulométrie tandis que la longueur de rupture de la nappe dépend principalement du premier et du dernier. La rotation du gaz modifie fortement la forme du profil de vitesse en sortie d'injecteur et donc du spray généré par la suite. Elle permet également une fluctuation plus faible de la longueur de rupture de la nappe et surtout une bonne homogénéisation de la vitesse des gouttes. Bien que ce type d'atomisation soit déjà particulièrement efficace, l'ajout d'une couronne de gaz externe améliore nettement l'atomisation.

The development of new materials with improved properties seems to rely nowadays more on blending and compounding than on the synthesis of chemically new polymers. Mixing may have a great effect on the morphology and structure of multi-component polymer materials. Twin-screw extruders (TSE) are widely used as mixers/reactors for blending, compounding, and reactive processing. This work aimed at developing a new instrument to measure in real time the residence time distribution (RTD) which characterizes the axial mixing and transport abilities of different screw elements based on the analysis of the transient flow pattern and systematic evaluation of mixing theory in TSE. Distributive mixing of polymer melts is characterized by the generation of interfacial area, which is experimentally much more difficult to measure. This 3D numerical simulation based on CFD is adopted, Aujourd'hui le développement de nouveaux matériaux polymères ayant de bonnes propriétés repose de plus en plus sur des procédés de mélange ou de compoundage de polymères au lieu de recourir à la synthèse de nouvelles molécules. L'action du mélange peut fortement influer sur la morphologie des matériaux polymères multi-constituants. Les extrudeuses bi-vis (TSE) sont souvent utilisées comme mélangeurs/réacteurs pour des procédés de mélange, de compoundage et d'extrusion réactive. Cependant, l'étude sur la qualité du mélange dans les TSE demeure un grand défi en raison de la complexité géométrique et du caractère transitoire de l'écoulement. Cette thèse a pour objet de développer un nouvel instrument en line pour mesurer en temps réel la distribution des temps de séjour (DTS) qui caractérise la performance du mélange axial et la capacité de convoyage de différents types d'éléments de vis basées sur l'analyse de l'écoulement transitoire et l'évaluation systématique de la théorie de mélange dans les TSE. Le mélange distributif des polymères fondus est caractérisé par la génération de l'aire des interfaces, un paramètre difficile à mesurer expérimentalement. Alors on fait appel à des simulations numériques de type CFD

One of the difficulties faced in turbulent flows comes from the fact that the statistics are heavily dependent on phenomena occurring in the vicinity of the walls. The presence of walls is usually accounted for by introducing damping functions in the model equations. This approach make's the model specific to a particular situation and does not favour it's universality. The present thesis aims at developing explicit algebraic models accounting for the wall effects by the elliptic blending approach developed by Manceau and Hanjalic, and at validating them in several flow configurations. The explicit algebraic methodology consists of projecting the anisotropy tensor equation onto a well chosen basis. The theory of an integrity basis enables the determination of the number of basis tensors necessary to avoid singularities in the model. Under this condition many models can be developed. Due to the accounting of the elliptic blending, an new tensor M appears in the equations in addition to the standard S and W tensors which, on the one hand, makes the theoretical framework more complex, but on the other hand, provides numerous possibilities of approximate models of practical interest. The linear and nonlinear models, validated in channel flow for a wide range of Reynolds numbers, in Couette- Poiseuille flow and in a shearless boundary layer, proved capable of correctly reproducing the near-wall anisotropy. The two-component limit is preserved. The extension of these models to 3D showed the possibility to represent correctly the anisotropy in a 3D case by using a 3-tensor truncated basis, provided that the complete set of invariants appearing in 3D are accounted for.; L'une des difficultés rencontrées dans les écoulements turbulents provient du fait que les statistiques dépendent très fortement de ce qui se passe à proximité des parois. La présence de ces derniers est pour la plupart du temps prise en compte par l'introduction de fonctions d'amortissement dans les équations des modèles. Cette approche particularise d'avantage les modèles et aussi ne favorise pas leur universalité. Le but de cette thèse est de développer des modèles Algébriques Explicites prenant en compte les effets de paroi par l'approche de la pondération elliptique développée par Manceau et Hanjalic, puis de procéder à la validation de ces modèles dans différentes configurations d'écoulements. La méthodologie de la modélisation algébrique explicite consiste à projeter l'équation du tenseur d'anisotropie sur une base convenablement choisie. La théorie des bases d'intégrité permet de déterminer le nombre de tenseurs devant constituer cette base afin d'éviter toute singularité dans les modèles. Dans cette circonstance de nombreux de modèles peuvent alors être développés. La prise en compte de la pondération elliptique à fait apparaître un tenseur supplémentaire M dans les équations en dehors des tenseurs S et W classiquement utilisés. Ce qui d'un côté rend plus complexe le cadre théorique, mais d'un autre côté offre de nombreuses possibilités de modèles approchés utilisables en pratique. Les modèles linéaire et non linéaire développés et testés pour une large gamme de Reynolds en canal, en écoulement de Couette-Poiseuille et en écoulement de couche limite sans cisaillement ont montrés leur capacité à reproduire convenablement l'anisotropie en proche paroi. La limite à deux composantes de la turbulence est préservée. L'extension de ces modèles en 3D a permis de montrer qu'il est possible de représenter correctement l'anisotropie dans un cas 3D en utilisant une base tronquée à trois tenseurs, à condition toutefois de prendre en compte l'intégralité des invariants apparaissant en 3D.

The study described in this thesis is about punctual injection of a saline solution inside a "Hele-Shaw cell" in order to characterize the dispersive behavior of a pollutant in porous media. The chosen experimental approach is based on the setup of an original Laser Induced Fluorescence (LIF) in the Hele-Shaw cell. The setting of the experimental apparatus allows quantitative data reduction of the experimental results. Moreover the "Moments Method" studied precisely the solution mixing dispersive behavior. Using the numerical code FRIPE the same injections have been simulated. The same quantitative data reductions have been applied to the numerical results. This led to an extensive comparison of the numerical and the experimental results, qualitatively but also of the dispersion in the mixing area of the injected solution, L'étude décrite dans cette thèse porte sur l'injection ponctuelle d'une solution saline au sein d'une cellule dite de Hele-Shaw, afin de caractériser le comportement dispersif d'un polluant en milieu poreux. L'approche expérimentale employée est basée sur l'implémentation originale d'un dispositif de Fluorescence Induite par Laser (LIF) dans la cellule. La mise en place d'un protocole de mesure efficace permet de mener une analyse quantitative des résultats expérimentaux. En outre, en appliquant la méthode des moments, il est possible de caractériser avec précision le comportement dispersif de la zone de mélange de la solution injectée. Parallèlement aux expériences, à l'aide du code numérique FRIPE, les injections ont été simulées numériquement. L'analyse quantitative a été appliquée à ces dernières. Une comparaison poussée des résultats expérimentaux et numériques a donc été effectuée, du point de vue qualitatif mais aussi sur l'expression de la dispersion du panache de la zone de mélange de la solution

Knowledge of shale poromechanical behavior is essential for various environmental issues such as CO2 deep geological storage, high level radioactive waste storage, oil field abandonment Furthermore, due to their good confinement properties, there are considered as potential host rocks for high level radioactive waste storage. However, their experimental characterization is very complex because of their very low permeability and their saturating fluid sensitivity. The main objective of this study is to improve knowledge of hydromechanical coupling of such a material. To study the structure of argillites under stress, SEM in-situ uniaxial compression test is performed on parallelepiped sample; compression orientation is perpendicular to bedding planes direction. Correlation techniques are used to estimate strain distributions. Second part underlines the key points of shale experimental characterization using the framework of Biot's mechanics of fluid saturated porous solids. Shales are well known as a more or less transverse isotropy material. The complete methodology to conduct dometric tests on such sensitive and weakly permeable material is described. Measurements realized on Tournemire argillite are proposed and a comparison with the poroelastic parameters of Meuse / Haute-Marne is given. Furthermore, to investigate hydromechanical coupling, a transverse isotropic poroelastic model is used to estimate the influence of anisotropy on diffusion characteristic time. The experimental data are used in numerical simulations to estimate the intrinsic permeability of the samples. In the last part, the caprock behavior of a petroleum field is studied at large scale with a 3D synthetic model (finite element analysis). The purpose of this last chapter is to study changes in stress field during the oil production; Une connaissance plus approfondie du comportement des roches argileuses, et plus particulièrement de leur tenue, s'avère aujourd'hui indispensable dans le cadre de plusieurs thèmes de recherche d'importance croissante : comme couverture des réservoirs pétroliers pour la production d'hydrocarbures et la séquestration de CO2, ou en tant que roche hôte pour le stockage des déchets nucléaires en couches géologiques profondes. Ces roches sont particulièrement difficiles à caractériser expérimentalement du fait de leur très faible perméabilité et de leur sensibilité à l'eau. L'objectif de cette thèse est d'améliorer la compréhension des couplages hydromécaniques intervenant dans ces matériaux. Tout d'abord, une contribution à la caractérisation de l'évolution de la structure des roches argileuses soumises à un chargement mécanique est apportée. A cet effet, un essai de compression uniaxiale sur l'argilite de Tournemire associant des observations au MEB et une méthode de corrélation d'images numériques est proposé. Par la suite, des simulations numériques des essais dométriques ont été réalisées pour appréhender l'influence des propriétés poromécaniques, en particulier de la faible perméabilité sur le comportement de l'échantillon. L'étude des équations de diffusion et les simulations numériques effectuées nous permettent de proposer une méthode de détermination de la perméabilité intrinsèque de ces matériaux. Dans la dernière partie, une étude du comportement d'une roche couverture d'un réservoir pétrolier à grande échelle est réalisée sur un cas synthétique 3D. L'objectif de ce chapitre est d'étudier les variations du champ de contrainte au sein de la couverture au cours de la production du réservoir

Knowledge of shale poromechanical behavior is essential for various environmental issues such as CO2 deep geological storage, high level radioactive waste storage, oil field abandonment Furthermore, due to their good confinement properties, there are considered as potential host rocks for high level radioactive waste storage. However, their experimental characterization is very complex because of their very low permeability and their saturating fluid sensitivity. The main objective of this study is to improve knowledge of hydromechanical coupling of such a material. To study the structure of argillites under stress, SEM in-situ uniaxial compression test is performed on parallelepiped sample; compression orientation is perpendicular to bedding planes direction. Correlation techniques are used to estimate strain distributions. Second part underlines the key points of shale experimental characterization using the framework of Biot's mechanics of fluid saturated porous solids. Shales are well known as a more or less transverse isotropy material. The complete methodology to conduct dometric tests on such sensitive and weakly permeable material is described. Measurements realized on Tournemire argillite are proposed and a comparison with the poroelastic parameters of Meuse / Haute-Marne is given. Furthermore, to investigate hydromechanical coupling, a transverse isotropic poroelastic model is used to estimate the influence of anisotropy on diffusion characteristic time. The experimental data are used in numerical simulations to estimate the intrinsic permeability of the samples. In the last part, the caprock behavior of a petroleum field is studied at large scale with a 3D synthetic model (finite element analysis). The purpose of this last chapter is to study changes in stress field during the oil production; Une connaissance plus approfondie du comportement des roches argileuses, et plus particulièrement de leur tenue, s'avère aujourd'hui indispensable dans le cadre de plusieurs thèmes de recherche d'importance croissante : comme couverture des réservoirs pétroliers pour la production d'hydrocarbures et la séquestration de CO2, ou en tant que roche hôte pour le stockage des déchets nucléaires en couches géologiques profondes. Ces roches sont particulièrement difficiles à caractériser expérimentalement du fait de leur très faible perméabilité et de leur sensibilité à l'eau. L'objectif de cette thèse est d'améliorer la compréhension des couplages hydromécaniques intervenant dans ces matériaux. Tout d'abord, une contribution à la caractérisation de l'évolution de la structure des roches argileuses soumises à un chargement mécanique est apportée. A cet effet, un essai de compression uniaxiale sur l'argilite de Tournemire associant des observations au MEB et une méthode de corrélation d'images numériques est proposé. Par la suite, des simulations numériques des essais dométriques ont été réalisées pour appréhender l'influence des propriétés poromécaniques, en particulier de la faible perméabilité sur le comportement de l'échantillon. L'étude des équations de diffusion et les simulations numériques effectuées nous permettent de proposer une méthode de détermination de la perméabilité intrinsèque de ces matériaux. Dans la dernière partie, une étude du comportement d'une roche couverture d'un réservoir pétrolier à grande échelle est réalisée sur un cas synthétique 3D. L'objectif de ce chapitre est d'étudier les variations du champ de contrainte au sein de la couverture au cours de la production du réservoir

One of the difficulties faced in turbulent flows comes from the fact that the statistics are heavily dependent on phenomena occurring in the vicinity of the walls. The presence of walls is usually accounted for by introducing damping functions in the model equations. This approach make's the model specific to a particular situation and does not favour it's universality. The present thesis aims at developing explicit algebraic models accounting for the wall effects by the elliptic blending approach developed by Manceau and Hanjalic, and at validating them in several flow configurations. The explicit algebraic methodology consists of projecting the anisotropy tensor equation onto a well chosen basis. The theory of an integrity basis enables the determination of the number of basis tensors necessary to avoid singularities in the model. Under this condition many models can be developed. Due to the accounting of the elliptic blending, an new tensor M appears in the equations in addition to the standard S and W tensors which, on the one hand, makes the theoretical framework more complex, but on the other hand, provides numerous possibilities of approximate models of practical interest. The linear and nonlinear models, validated in channel flow for a wide range of Reynolds numbers, in Couette- Poiseuille flow and in a shearless boundary layer, proved capable of correctly reproducing the near-wall anisotropy. The two-component limit is preserved. The extension of these models to 3D showed the possibility to represent correctly the anisotropy in a 3D case by using a 3-tensor truncated basis, provided that the complete set of invariants appearing in 3D are accounted for., L'une des difficultés rencontrées dans les écoulements turbulents provient du fait que les statistiques dépendent très fortement de ce qui se passe à proximité des parois. La présence de ces derniers est pour la plupart du temps prise en compte par l'introduction de fonctions d'amortissement dans les équations des modèles. Cette approche particularise d'avantage les modèles et aussi ne favorise pas leur universalité. Le but de cette thèse est de développer des modèles Algébriques Explicites prenant en compte les effets de paroi par l'approche de la pondération elliptique développée par Manceau et Hanjalic, puis de procéder à la validation de ces modèles dans différentes configurations d'écoulements. La méthodologie de la modélisation algébrique explicite consiste à projeter l'équation du tenseur d'anisotropie sur une base convenablement choisie. La théorie des bases d'intégrité permet de déterminer le nombre de tenseurs devant constituer cette base afin d'éviter toute singularité dans les modèles. Dans cette circonstance de nombreux de modèles peuvent alors être développés. La prise en compte de la pondération elliptique à fait apparaître un tenseur supplémentaire M dans les équations en dehors des tenseurs S et W classiquement utilisés. Ce qui d'un côté rend plus complexe le cadre théorique, mais d'un autre côté offre de nombreuses possibilités de modèles approchés utilisables en pratique. Les modèles linéaire et non linéaire développés et testés pour une large gamme de Reynolds en canal, en écoulement de Couette-Poiseuille et en écoulement de couche limite sans cisaillement ont montrés leur capacité à reproduire convenablement l'anisotropie en proche paroi. La limite à deux composantes de la turbulence est préservée. L'extension de ces modèles en 3D a permis de montrer qu'il est possible de représenter correctement l'anisotropie dans un cas 3D en utilisant une base tronquée à trois tenseurs, à condition toutefois de prendre en compte l'intégralité des invariants apparaissant en 3D.

The study described in this thesis is about punctual injection of a saline solution inside a "Hele-Shaw cell" in order to characterize the dispersive behavior of a pollutant in porous media. The chosen experimental approach is based on the setup of an original Laser Induced Fluorescence (LIF) in the Hele-Shaw cell. The setting of the experimental apparatus allows quantitative data reduction of the experimental results. Moreover the "Moments Method" studied precisely the solution mixing dispersive behavior. Using the numerical code FRIPE the same injections have been simulated. The same quantitative data reductions have been applied to the numerical results. This led to an extensive comparison of the numerical and the experimental results, qualitatively but also of the dispersion in the mixing area of the injected solution, L'étude décrite dans cette thèse porte sur l'injection ponctuelle d'une solution saline au sein d'une cellule dite de Hele-Shaw, afin de caractériser le comportement dispersif d'un polluant en milieu poreux. L'approche expérimentale employée est basée sur l'implémentation originale d'un dispositif de Fluorescence Induite par Laser (LIF) dans la cellule. La mise en place d'un protocole de mesure efficace permet de mener une analyse quantitative des résultats expérimentaux. En outre, en appliquant la méthode des moments, il est possible de caractériser avec précision le comportement dispersif de la zone de mélange de la solution injectée. Parallèlement aux expériences, à l'aide du code numérique FRIPE, les injections ont été simulées numériquement. L'analyse quantitative a été appliquée à ces dernières. Une comparaison poussée des résultats expérimentaux et numériques a donc été effectuée, du point de vue qualitatif mais aussi sur l'expression de la dispersion du panache de la zone de mélange de la solution

The present work deals with the linear and non-linear thermo-hydro-mechanical behaviour of porous rocks such as the argillite by the multiscale modelling approach. Based on microstructure observations, a conceptual model was proposed. In this model, the representative elementary volume of a heterogeneous medium is composed of an argillaceous matrix containing spherical inclusions of minerals quartz and calcite and ellipsoidal inclusions representing the pore space. In a first step, the process of modelling has been exploited by determining the isotropic and transversely isotropic effective properties of the argillite: thermal conductivity and thermo-hydro-mechanical properties. Furthermore, many numerical studies have highlighted the influence of the morphology of the pore space, of the mineralogy and of the estimate schemes to the predictive results. In a second step, we modelled the non linear mechanical behaviour (elasto-plastic, elasto-viscoplastic) of argillaceous rocks. The comparison between numerical simulations and available experimental results (triaxial compression test, creep test) confirmed the validation of the model developed; Le présent travail est consacré à l'étude du comportement thermo-hydro-mécanique linéaire et non linéaire des roches poreuses de type argilites par approche de changement d'échelle. A partir des observations microstructurales de ces matériaux, un modèle conceptuel a été proposé. Dans ce modèle, le volume élémentaire représentatif du milieu hétérogène est composé d'une phase matricielle argileuse contenant des inclusions sphériques de minéraux de quartz et de calcite et des inclusions ellipsoïdales aplaties représentant l'espace poreux. Dans un premier temps, le procédé de la modélisation a été exploité par la détermination des propriétés effectives isotropes et isotrope transverses des argilites : la conductivité thermique et les propriétés thermo-hydro-mécanique croisées. En outre, de nombreuses études numériques ont mis en évidence l'influence de la morphologie de l'espace poreux, de la minéralogie et des schémas d'estimation sur les résultats prédictifs. Dans un deuxième temps, nous avons modélisé le comportement mécanique non linéaire (élasto-plastique, élasto-viscoplastique) des roches argileuses. La comparaison entre les simulations numériques et les résultats expérimentaux disponibles (essai de compression triaxiale, essai de fluage) a confirmé la validation du modèle développé

The present work deals with the linear and non-linear thermo-hydro-mechanical behaviour of porous rocks such as the argillite by the multiscale modelling approach. Based on microstructure observations, a conceptual model was proposed. In this model, the representative elementary volume of a heterogeneous medium is composed of an argillaceous matrix containing spherical inclusions of minerals quartz and calcite and ellipsoidal inclusions representing the pore space. In a first step, the process of modelling has been exploited by determining the isotropic and transversely isotropic effective properties of the argillite: thermal conductivity and thermo-hydro-mechanical properties. Furthermore, many numerical studies have highlighted the influence of the morphology of the pore space, of the mineralogy and of the estimate schemes to the predictive results. In a second step, we modelled the non linear mechanical behaviour (elasto-plastic, elasto-viscoplastic) of argillaceous rocks. The comparison between numerical simulations and available experimental results (triaxial compression test, creep test) confirmed the validation of the model developed, Le présent travail est consacré à l'étude du comportement thermo-hydro-mécanique linéaire et non linéaire des roches poreuses de type argilites par approche de changement d'échelle. A partir des observations microstructurales de ces matériaux, un modèle conceptuel a été proposé. Dans ce modèle, le volume élémentaire représentatif du milieu hétérogène est composé d'une phase matricielle argileuse contenant des inclusions sphériques de minéraux de quartz et de calcite et des inclusions ellipsoïdales aplaties représentant l'espace poreux. Dans un premier temps, le procédé de la modélisation a été exploité par la détermination des propriétés effectives isotropes et isotrope transverses des argilites : la conductivité thermique et les propriétés thermo-hydro-mécanique croisées. En outre, de nombreuses études numériques ont mis en évidence l'influence de la morphologie de l'espace poreux, de la minéralogie et des schémas d'estimation sur les résultats prédictifs. Dans un deuxième temps, nous avons modélisé le comportement mécanique non linéaire (élasto-plastique, élasto-viscoplastique) des roches argileuses. La comparaison entre les simulations numériques et les résultats expérimentaux disponibles (essai de compression triaxiale, essai de fluage) a confirmé la validation du modèle développé

This work has two main objectives. The first one is to develop an integrated methodology allowing to build a 3D geomechanical model and also to image the uncertainties attached to the poro-mechanical properties of the constitutive rocks. This geomechanical model should be based on all related available data and should be consistent with the static and dynamic models, currently built by reservoir geologists and engineers. The second objective is to analyse the impact of geological heterogeneities, which are often neglected, in the mechanical response of the reservoir induced by its exploitation, and furthermore to derive uncertainties on the stress and deformation fields related to the uncertainties on the input properties of the geomechanical model. An integrated methodology based on geostatistical simulations is developed. First, the geometric frame is built; then an approach of embedded stochastic simulations is carried out to infill the different reservoir properties, the lithological description constraining the petrophysical and poro-elastic descriptions. The next step is to generate the mechanical responses of the stochastic realisations, using a finite-element mechanical simulator. The uncertainties on the resulting stress and displacement fields are then deduced from the multiple mechanical responses which are computed. This approach is demonstrated on a real field case, a fluvio-deltaic reservoir in North Sea. It is shown on this example that the reservoir heterogeneities and their uncertainties significantly influence the calculations of stress and strain fields, and also the risks of mechanical failure. Uncertainties on the mechanical quantities under analysis (first invariant of the stress tensor and subsidence) are also derived, L'objectif de ce travail est double : d'une part, il s'agit de développer une méthodologie intégrée pour la construction d'un modèle géomécanique ainsi que la représentation des incertitudes associées aux propriétés poro-élastiques des roches constitutives, en exploitant l'ensemble des données disponibles et en s'appuyant de façon cohérente sur les modèles de réservoir statique et dynamique classiquement utilisés par les géologues et les ingénieurs réservoir ; d'autre part, il s'agit d'analyser quel est l'impact des hétérogénéités géologiques, souvent négligées, dans la réponse mécanique du réservoir sollicité par son exploitation, et d'aboutir à des incertitudes sur les champs de contraintes et de déplacements, issues des incertitudes sur ces hétérogénéités et leurs paramètres mécaniques. Pour ce faire, une méthodologie intégrée s'appuyant sur des simulations géostatistiques a été développée. Après la construction du cadre géométrique 3D, le remplissage des propriétés au sein du réservoir suit une démarche de simulations géostatistiques 3D emboîtées, dans laquelle la représentation des hétérogénéités lithologiques conditionne la génération des propriétés poro-élastiques. La démarche consiste ensuite à représenter les incertitudes sur le modèle géomécanique par des ensembles de réalisations géostatistiques dont la réponse mécanique est alors calculée avec un simulateur mécanique aux éléments finis. Les incertitudes sur les champs de contraintes et de déformations sont déduites ensuite des différentes réponses mécaniques obtenues. La démarche a été mise en oeuvre sur un réservoir réel, dans un environnement fluvio-deltaïque, produisant en Mer du Nord. Dans ce cadre, il a été démontré que les hétérogénéités du réservoir et leurs incertitudes influencent significativement les calculs des champs de contraintes et de déformations, ainsi que les risques mécaniques de rupture. Des incertitudes sur les quantités mécaniques analysées (premier invariant du tenseur des contraintes et subsidence) ont été aussi estimées

Inclusions removal by flotation in mechanically agitated vessels is widely used in liquid aluminium treatments. This process consists of gas injection in a molten aluminium reactor using an impeller. Inclusions are attached to the bubbles while rising into the liquid bulk and released in the froth at the vessel surface.With the aim of both a better understanding of the physical processes acting during flotation and the optimization of the refining process, a 2D mathematical modelling of the flotation tank was set up. Transport phenomena, sedimentation, agglomeration and flotation of inclusions were considered here. The model couples Population Balance with convective transport of the inclusions, in the multiphase flow field. Simulations showed that a good dispersion of the bubbles into the tank, especially in the high turbulence zones, is essential for an efficient process. This situation was achieved with a strong agitation. Gas flow affected the process efficiency weakly.Theoretical and experimental studies on chlorine effects were done. It was noticed that a critical amount of chlorine into the melt is required for the flotation to be efficient. Complete understanding of chlorine role is not achieved yet, but promising ideas for future research have been exhibited, L'élimination des inclusions indésirables par flottation en cuve agitée est un procédé largement utilisé dans l'industrie de l'aluminium. Ce procédé consiste en l'injection de gaz au sein de la cuve par le biais d'un rotor ; les bulles ainsi formées captent, lors de leur ascension dans le métal liquide, les inclusions et les drainent en surface où elles sont éliminées mécaniquement.Dans le but de prédire l'évolution numérale de la population d'inclusions dans la cuve de flottation au cours du temps, un modèle 2D a été développé. Afin de rendre compte du transport convectif des inclusions au sein de l'écoulement diphasique métal liquide bulles, des phénomènes de décantation, d'agglomération et de flottation, ce modèle couple la Mécanique des Fluides Numérique et la technique des Bilans de Population. Une attention particulière a été portée aux modèles d'agglomération et de flottation pour qu'ils soient applicables sur une grande plage de conditions de turbulence.Les simulations numériques réalisées, nous ont notamment permis de mettre en évidence l'importance de veiller à une bonne répartition des bulles dans la cuve, en particulier dans les zones de forte turbulence. En parallèle de ce travail de modélisation, une analyse de la littérature et des travaux expérimentaux sur métal liquide ont été réalisés, afin de mieux comprendre le rôle du chlore dans le procédé de flottation. Un effet de seuil a ainsi pu être mis en évidence : une certaine quantité de chlore injectée dans le métal est nécessaire pour que la flottation devienne efficace. La compréhension du rôle du chlore reste cependant partielle, mais des pistes de recherche prometteuses ont put être identifiées

Inclusions removal by flotation in mechanically agitated vessels is widely used in liquid aluminium treatments. This process consists of gas injection in a molten aluminium reactor using an impeller. Inclusions are attached to the bubbles while rising into the liquid bulk and released in the froth at the vessel surface.With the aim of both a better understanding of the physical processes acting during flotation and the optimization of the refining process, a 2D mathematical modelling of the flotation tank was set up. Transport phenomena, sedimentation, agglomeration and flotation of inclusions were considered here. The model couples Population Balance with convective transport of the inclusions, in the multiphase flow field. Simulations showed that a good dispersion of the bubbles into the tank, especially in the high turbulence zones, is essential for an efficient process. This situation was achieved with a strong agitation. Gas flow affected the process efficiency weakly.Theoretical and experimental studies on chlorine effects were done. It was noticed that a critical amount of chlorine into the melt is required for the flotation to be efficient. Complete understanding of chlorine role is not achieved yet, but promising ideas for future research have been exhibited; L'élimination des inclusions indésirables par flottation en cuve agitée est un procédé largement utilisé dans l'industrie de l'aluminium. Ce procédé consiste en l'injection de gaz au sein de la cuve par le biais d'un rotor ; les bulles ainsi formées captent, lors de leur ascension dans le métal liquide, les inclusions et les drainent en surface où elles sont éliminées mécaniquement.Dans le but de prédire l'évolution numérale de la population d'inclusions dans la cuve de flottation au cours du temps, un modèle 2D a été développé. Afin de rendre compte du transport convectif des inclusions au sein de l'écoulement diphasique métal liquide bulles, des phénomènes de décantation, d'agglomération et de flottation, ce modèle couple la Mécanique des Fluides Numérique et la technique des Bilans de Population. Une attention particulière a été portée aux modèles d'agglomération et de flottation pour qu'ils soient applicables sur une grande plage de conditions de turbulence.Les simulations numériques réalisées, nous ont notamment permis de mettre en évidence l'importance de veiller à une bonne répartition des bulles dans la cuve, en particulier dans les zones de forte turbulence. En parallèle de ce travail de modélisation, une analyse de la littérature et des travaux expérimentaux sur métal liquide ont été réalisés, afin de mieux comprendre le rôle du chlore dans le procédé de flottation. Un effet de seuil a ainsi pu être mis en évidence : une certaine quantité de chlore injectée dans le métal est nécessaire pour que la flottation devienne efficace. La compréhension du rôle du chlore reste cependant partielle, mais des pistes de recherche prometteuses ont put être identifiées

Nature is complex and must be approached in total modesty by engineers seeking to predict the behavior of underground openings. The engineering of industrial projects in underground situations, with high economic and social stakes (Alpine mountain crossings, nuclear waste repository), mean striving to gain better understanding of the behavioral mechanisms of the openings to be designed. This improvement necessarily involves better physical representativeness of macroscopic mechanisms and the provision of prediction tools suited to the expectations and needs of the engineers. The calculation tools developed in this work is in step with this concern for satisfying industrial needs and developing knowledges related to the rheology of geomaterials. These developments led to the proposing of a mechanical constitutive model, suited to lightly fissured rocks, comparable to continuous media, while integrating more particularly the effect of time.Thread of this study, the problematics ensued from the subject of the thesis is precisely about the rock mass delayed behavior in numerical modeling and its consequences on underground openings design.Based on physical concepts of reference, defined in several scales (macro/meso/micro), the developed constitutive model is translated in a mathematical formalism in order to be numerically implemented.Numerical applications presented as illustrations fall mainly within the framework of nuclear waste repository problems. They concern two very different configurations of underground openings: the AECL's underground canadian laboratory, excavated in the Lac du Bonnet granite, and the GMR gallery of Bure's laboratory (Meuse/Haute-Marne), dug in argillaceous rock.In this two cases, this constitutive model use highlights the gains to be obtained from allowing for delayed behavior regarding the accuracy of numerical tunnel behavior predictions in the short, medium and long terms, La nature est complexe et c'est en toute modestie que les ingénieurs doivent chercher à prédire le comportement des ouvrages dans le sous-sol. La réalisation de projets industriels dans le domaine souterrain, à forts enjeux économiques et sociaux (traversées alpines, stockage de déchets nucléaires), nécessite d'évoluer vers une meilleure compréhension des mécanismes comportementaux des ouvrages à concevoir. Cette amélioration passe par une meilleure représentativité physique des mécanismes macroscopiques et par la mise à disposition d'outils de prédiction adaptés aux attentes et aux besoins des ingénieurs. Les outils de calculs développés dans ce travail s'inscrivent dans cette volonté de rapprocher les attentes de l'industrie et les connaissances liées à la rhéologie des géomatériaux. Ces développements ont ainsi débouché sur la proposition d'un modèle de comportement mécanique, adapté aux roches peu fissurées et assimilables à des milieux continus, intégrant, en particulier, l'effet du temps. Fil conducteur de cette étude, la problématique du sujet de thèse concerne précisément la prise en compte du comportement différé des massifs rocheux dans les modélisations et ses conséquences sur les ouvrages souterrains.Fondé sur des concepts physiques de référence, définis à différentes échelles (macro/méso/micro), le modèle rhéologique développé est transcrit dans un formalisme mathématique dans le but d'être mis en oeuvre numériquement.Les applications numériques proposées s'inscrivent principalement dans le contexte du stockage des déchets radioactifs. Elles concernent deux configurations d'ouvrages rigoureusement différentes : l'excavation du laboratoire souterrain canadien de l'AECL, dans le granite du Lac du Bonnet et le creusement de la galerie GMR du laboratoire de Bure (Meuse/Haute-Marne) dans l'argilite de l'Est. Dans les deux cas, l'utilisation du modèle a permis de mettre en évidence l'apport de la prise en compte du comportement différé sur la représentativité des prédictions numériques du comportement à court, moyen et long termes des ouvrages souterrains

This thesis contributes to the analysis of the anisotropic effects in the determination of the poroelastic properties of the porous rocks. After a bibliographical study, the influence of the form and distribution in orientation of the pores on the anisotropy of the rocks were studied. The homogenization scheme in this thesis is that of Mori-Tanaka. The estimate of the effective coefficients poroelastic of heterogeneous materials of type rocks-like composites characterized by a transverse isotropic matrix containing solid mineral inclusions of roughly spherical form was carried out in using the analytical solution of a spherical inclusion isolated in an infinite transverse isotropic matrix. Thereafter, we are interested in the taking into account of spheroidal inclusions of which orientation does not coincid with the orthogonal axis of revolution of the matrix. We propose a numerical approach of integration based on the Green function. Numerical integrate on the unit sphere is carried out using a method of Gauss of which the precision is discussed. The developed numerical tool is applied to a porous rock by considering a homogenization scheme in two stages and various functions of distribution of orientation. Respective influences of pore space anisotropy and matrix anisotropy are clearly distinguished and analyzed.; Cette thèse apporte une contribution à l'analyse des effets d'anisotropie dans la détermination des propriétés poroélastiques des roches poreuses. Après une étude bibliographique, l'influence de la forme et de la distribution en orientation des pores sur l'anisotropie des roches a été étudiée. Le schéma d'estimation pris en compte est celui de Mori-Tanaka. En utilisant la solution analytique d'une inclusion sphérique isolée dans une matrice isotrope transverse infinie, l'estimation des coefficients poroélastiques effectifs des matériaux hétérogènes de type roches-composites a été réalisée. Par la suite, on s'intéresse à la prise en compte des inclusions sphéroidales dont l'orientation ne coincide pas avec l'axe d'orthotropie de révolution de la matrice. On propose une approche d'intégration numérique basée sur la fonction de Green. L'intégration numérique sur la sphère unité est réalisée à l'aide d'une méthode de Gauss dont la précision est discutée. L'outil numérique développé est appliqué à une roche poreuse en considérant un schéma d'homogénéisation en deux étapes et différentes fonctions de distribution d'orientation. Les résultats obtenus mettent en évidence les influences respectives des anisotropies de matrice et de l'espace poreux.

Our aim is here a better understanding of damage process of stressed rocks. We have therefore conducted a micro-macro study of the rock damage in the sense that we have realized microscopic and macroscopic measurements of its evolution during mechanical tests. This study has consisted of : - measurements of elastic modulus before the macroscopic rupture ; - measurements of the velocity of elastic waves and attenuation ; - measurements of acoustic activity (AE events). This study led to the following results : - Results from the three types of measurements are consistent. - Damage of rocks which have different macroscopic behaviours follows from micro-mechanisms which are also different. For dilative rocks, we observe tensile micro-mechanisms (Mode I). For not dilative rocks, damage results from shear microscopic process. - A correlation between micro-mechanisms type (compressive, tensile, shear) and the frequency of AE events is observed. Compressive microscopic process such as microcrack closure ar low frequency. Tensile AE sources are high frequency. Shear micromechanisms are characterized by intermédiate frequency, Afin de mieux comprendre les processus d'endommagement des roches sollicitées mécaniquement, nous avons mené une étude « micro-macro » de l'endommagement au sens où nous avons réalisé des mesures microscopiques et macroscopiques de son évolution durant un essai mécanique. Cette étude a ainsi comporté : - des mesures de modules élastiques avant la rupture macroscopique ; - des mesures de vitesses de propagation des ondes élastiques et d'atténuation ; - des mesures d'activité acoustique (événements d'EA). De cette étude expérimentale, se dégagent les résultats suivants : - Il existe une remarquable cohérence entre les trois types de mesures effectuées. - L'endommagement de roches aux comportements macroscopiques différents découle de micromécanismes eux aussi différents. Pour des roches dilatantes, on observe des micromécanismes en extension (Mode I). L'endommagement de roches non-dilatantes résulte de micromécanismes en cisaillement. - Une corrélation entre le contenu fréquentiel des événements d'EA et le type de micromécanismes est observée. Les micromécanismes en extension (Mode I) ont le contenu fréquentiel le plus élevé. Les micromécanismes en compression (fermeture de microfissures, de vides inter-granulaires) ont un contenu fréquentiel basse fréquence. Les micromécanismes de cisaillement ont un contenu fréquentiel intermédiaire aux deux précédents.

This thesis contributes to the analysis of the anisotropic effects in the determination of the poroelastic properties of the porous rocks. After a bibliographical study, the influence of the form and distribution in orientation of the pores on the anisotropy of the rocks were studied. The homogenization scheme in this thesis is that of Mori-Tanaka. The estimate of the effective coefficients poroelastic of heterogeneous materials of type rocks-like composites characterized by a transverse isotropic matrix containing solid mineral inclusions of roughly spherical form was carried out in using the analytical solution of a spherical inclusion isolated in an infinite transverse isotropic matrix. Thereafter, we are interested in the taking into account of spheroidal inclusions of which orientation does not coincid with the orthogonal axis of revolution of the matrix. We propose a numerical approach of integration based on the Green function. Numerical integrate on the unit sphere is carried out using a method of Gauss of which the precision is discussed. The developed numerical tool is applied to a porous rock by considering a homogenization scheme in two stages and various functions of distribution of orientation. Respective influences of pore space anisotropy and matrix anisotropy are clearly distinguished and analyzed., Cette thèse apporte une contribution à l'analyse des effets d'anisotropie dans la détermination des propriétés poroélastiques des roches poreuses. Après une étude bibliographique, l'influence de la forme et de la distribution en orientation des pores sur l'anisotropie des roches a été étudiée. Le schéma d'estimation pris en compte est celui de Mori-Tanaka.En utilisant la solution analytique d'une inclusion sphérique isolée dans une matrice isotrope transverse infinie, l'estimation des coefficients poroélastiques effectifs des matériaux hétérogènes de type roches-composites a été réalisée.Par la suite, on s'intéresse à la prise en compte des inclusions sphéroidales dont l'orientation ne coïncide pas avec l'axe d'orthotropie de révolution de la matrice. On propose une approche d'intégration numérique basée sur la fonction de Green. L'intégration numérique sur la sphère unité est réalisée à l'aide d'une méthode de Gauss dont la précision est discutée. L'outil numérique développé est appliqué à une roche poreuse en considérant un schéma d'homogénéisation en deux étapes et differentes fonctions de distribution d'orientation. Les résultats obtenus mettent en évidence les influences respectives des anisotropies de matrice et de l'espace poreux.

Natural fractures have a dramatic impact on reservoirs either in fields such as mining and petroleum engineering or hydrogeology. Indeed, acting as conductors or barriers, they often control the flow and also the storage capacity in such reservoirs. The presented work is enclosed in the numerical modelling line of naturally fractured reservoirs, addressing the specific issues of characterizing and modelling fractures in three dimensions. The proposed approaches successively deal with the first logical stages of this modelling: 1. Natural fractures characterization from structural data: a geomechanical and probabilistic approach in three dimensions integrates tensorial strain data (obtained from forward or backward modelling techniques) in order to compute new structural drivers of fracturing intensity and orientation. The proposed mathematical framework also account for petrophysical rock property uncertainties. 2. Numerical fracture network modelling: a discrete fracture network modelling approach in three dimensions was developed because of its flexibility and its realism. A numerical stochastic simulation enables to model any kind of fracture network while honouring field data and predefined geological constraints (such as strata confinement, relative chronology between several fracture systems, etc.). 3. Connectivity analysis of discrete fracture networks: a graph structure highlights the connectivity relationships between the discrete fractures of a three-dimensional network. Several connectivity analysis methods, qualitative and quantitative, are then outlined with an emphasis on their improvement on understanding of the fracturing impact on fluid flows, Les fractures naturelles ont un impact considérable au sein d'un réservoir, que ce soit dans les domaines de l'ingénierie pétrolière et minière, ou de l'hydrogéologie. En effet, qu'elles soient conductrices ou barrières, elles contrôlent généralement les écoulements ainsi que les capacités de stockage des réservoirs. Ces travaux de thèse s'inscrivent au sein de la chaîne de modélisation numérique d'un réservoir naturellement fracturé, leurs objectifs étant de caractériser, puis de modéliser les fractures en trois dimensions. Les méthodologies proposées abordent successivement les premières étapes logiques de cette modélisation: 1. La caractérisation des fractures naturelles à partir de données structurales: les données tensorielles de déformation (résultant d'une modélisation génétique ou inverse) sont intégrées au sein d'une approche tridimensionnelle géomécanique et probabiliste afin de déterminer de nouveaux indicateurs structuraux de fracturation, notamment concernant son intensité et son orientation. La formulation mathématique adoptée prend également en compte les incertitudes sur les propriétés pétrophysiques des roches. 2. La modélisation numérique d'un réseau de fractures: une approche tridimensionnelle par réseau de fractures discrètes a été développée du fait de sa flexibilité et de son réalisme. Une simulation numérique stochastique offre la possibilité de modéliser tout type de réseau de fractures respectant des données de terrains et des contraintes géologiques prédéfinies (confinement aux couches, chronologie relative des fractures, etc.). 3. L'analyse de la connectivité de réseaux de fractures discrètes: une structure de type graphe reflète les relations de connectivité entre les fractures discrètes d'un réseau en trois dimensions. Différentes méthodologies d'analyse qualitative et quantitative de cette connectivité sont ensuite exposées, leur but étant d'améliorer la compréhension de l'impact de la fracturation sur les écoulements fluides

The aim of the paper is to describe the mechanisms which govern the interstitial flow in a semi-solid medium and its influence on the behaviour of the skeleton. A model for the compressible viscoplastic solid coupled with the liquid interstitial flow was studied. A dimensional analysis of this model was made. This analysis shows two dimensionless parameters and which characterize the model ; the first represents compressibility and the second characterizes the liquid influence on the solid and it depends on the characteristic velocity of the problem and on the ratio between the dynamic viscosity of the liquid and the apparent viscosity of the solid. In order to understand the behaviour of these biphasic models, we studied the case of an unidimensional deformation. We have found an analytical solution for the unidimensional problems. These analytical studies show two types of behaviour of the medium according to the order of magnitude of the dimensionless number . At low velocity, the fluid does not have any influence on the response of the solid skeleton. At high velocity, the flow is very localised close to the edges where the fluid can run out. By introducing the elasticity of the skeleton into the preceding model, we found an analytical approached solution. This analysis shows an additional parameter : a relaxation time which characterizes this model. We have showed that this time is proportional to the inverse of the number ., L'objectif de cette thèse est de décrire les mécanismes qui régissent l'écoulement interstitiel dans un milieu semi-solide, ainsi que son influence sur le comportement du squelette. L'étude porte sur un modèle de solide viscoplastique compressible couplé à l'écoulement interstitiel d'un liquide. Une analyse dimensionnelle de ce modèle a été faite. Cette analyse montre deux paramètres sans dimension et qui caractérisent le modèle; le premier représente la compressibilité et le deuxième caractérise l'influence du liquide sur le solide et il dépend de la vitesse caractéristique du problème et du rapport entre la viscosité dynamique du liquide et la viscosité apparente du solide. Pour bien comprendre le comportement de ces modèles biphasiques, nous avons étudié le cas d'une déformation unidimensionnelle. Nous avons montré que les solutions des problèmes unidimensionnels peuvent être menées de manière complètement analytique. Ces études analytiques ont permis de distinguer deux types de comportement du milieu selon l'ordre de grandeur du nombre sans dimension . à faible vitesse, le fluide n'a pas d'influence sur la réponse du squelette solide. à grande vitesse, l'écoulement est très localisé près des bords où le fluide peut s'écouler. En introduisant l'élasticité du squelette dans le modèle précédent, nous avons trouvé une solution analytique approchée. Cette analyse montre un paramètre supplémentaire: un temps de relaxation qui caractérise ce modèle. Ce temps est inversement proportionnel au nombre

This dissertation forms part of a larger research field which deals with using seismic monitoring systems to better understand the behavior of natural rock slope appealing to the laboratory and in-situ experiments and the numerical simulation. For this purpose, a small scale rock mass, former instrumented from the geo-hydro-mechanical point of view, bas been chosen for installing a seismic monitoring system. ln this site, it is possible to easily induce a known hydraulic charging and to measure displacements and pressures in many meaningful points. The results of the in-situ experiments highlight the existence of micro-seismicity induced by a variation of water pressures and the repetitiveness either of the hydro-mechanical and micro-seismic behavior. From results, it bas been hypothesized that the micro-seismic activity onsets when water pressure increases in less permeable zones. The Finite Elements numerical modeling of the site allows to calculate the stress variations due to hydraulic charging and to assess if stress variations could generate seismic sources. Series of cores have been taken from the site in order to tested rock material and rock fractures in laboratory, associated with acoustic emissions. lt results that little stress variations due to water pressure, calculated for the site, do not explain the micro-seismicity observed in laboratory. A micro-mechanical model using Boundary Elements Method (Displacements Discontinuity Method) which simulates the hydro-mechanical behaviour of rock joints sample under normal compression is proposed taking into account the joint topography. Results confirm the laboratory observations and show that the hydrostatic charging can induce instabilities, when contact areas open, comparable to the ill-situ micro-seismicity., Un massif rocheux fracturé, dont le comportement géo-hydromécanique a été étudié depuis plusieurs années, a été instrumenté sismiquement. L'intérêt du site réside dans le fait qu'on peut induire artificiellement une charge hydraulique connue et mesurer des déplacements, des pressions hydrauliques sur plusieurs points significatifs. Les résultats des essais in-situ mettent en évidence l'existence d'une micro-sismicité induite par une variation de pression avec une grande repétitivité. A partir des résultats des essais in-situ, a été formulée l'hypothèse que l'activité micro-sismique est associée à l'augmentation de la pression d'eau dans les zones moins perméables du massif rocheux. Une modélisation numérique du site avec la méthode des Éléments Finis a été réalisée pour déterminer l'amplitude des variations de contraintes induites par la mise en charge hydraulique et vérifier si elles pouvaient engendrer des sources sismiques. Des échantillons prélevés ont permis de déterminer, en laboratoire, le comportement mécanique de la roche intacte et des joints rocheux et la sismicité associée. il s'avère que les variations de contrainte calculées pour le site ne sont pas suffisantes pour générer la sismicité observée au laboratoire. Une modélisation micromécanique à l'échelle du joint rocheux avec la méthode des Éléments de Contour (Discontinuités de Déplacement) a été réalisée pour simuler le comportement hydro-mécanique tenant compte de la morphologie du joint. Les résultats numériques sont en accord avec les obsentations de laboratoire et montrent que la mise en charge hydraulique d'un joint rugueux peut engendrer des instabilités lors de l'ouverture des aires de contact compatibles avec la micro-sismicité observée sur le site.

This dissertation forms part of a larger research field which deals with using seismic monitoring systems to better understand the behavior of natural rock slope appealing to the laboratory and in-situ experiments and the numerical simulation. For this purpose, a small scale rock mass, former instrumented from the geo-hydro-mechanical point of view, bas been chosen for installing a seismic monitoring system. ln this site, it is possible to easily induce a known hydraulic charging and to measure displacements and pressures in many meaningful points. The results of the in-situ experiments highlight the existence of micro-seismicity induced by a variation of water pressures and the repetitiveness either of the hydro-mechanical and micro-seismic behavior. From results, it bas been hypothesized that the micro-seismic activity onsets when water pressure increases in less permeable zones. The Finite Elements numerical modeling of the site allows to calculate the stress variations due to hydraulic charging and to assess if stress variations could generate seismic sources. Series of cores have been taken from the site in order to tested rock material and rock fractures in laboratory, associated with acoustic emissions. lt results that little stress variations due to water pressure, calculated for the site, do not explain the micro-seismicity observed in laboratory. A micro-mechanical model using Boundary Elements Method (Displacements Discontinuity Method) which simulates the hydro-mechanical behaviour of rock joints sample under normal compression is proposed taking into account the joint topography. Results confirm the laboratory observations and show that the hydrostatic charging can induce instabilities, when contact areas open, comparable to the ill-situ micro-seismicity., Un massif rocheux fracturé, dont le comportement géo-hydromécanique a été étudié depuis plusieurs années, a été instrumenté sismiquement. L'intérêt du site réside dans le fait qu'on peut induire artificiellement une charge hydraulique connue et mesurer des déplacements, des pressions hydrauliques sur plusieurs points significatifs. Les résultats des essais in-situ mettent en évidence l'existence d'une micro-sismicité induite par une variation de pression avec une grande repétitivité. A partir des résultats des essais in-situ, a été formulée l'hypothèse que l'activité micro-sismique est associée à l'augmentation de la pression d'eau dans les zones moins perméables du massif rocheux. Une modélisation numérique du site avec la méthode des Éléments Finis a été réalisée pour déterminer l'amplitude des variations de contraintes induites par la mise en charge hydraulique et vérifier si elles pouvaient engendrer des sources sismiques. Des échantillons prélevés ont permis de déterminer, en laboratoire, le comportement mécanique de la roche intacte et des joints rocheux et la sismicité associée. il s'avère que les variations de contrainte calculées pour le site ne sont pas suffisantes pour générer la sismicité observée au laboratoire. Une modélisation micromécanique à l'échelle du joint rocheux avec la méthode des Éléments de Contour (Discontinuités de Déplacement) a été réalisée pour simuler le comportement hydro-mécanique tenant compte de la morphologie du joint. Les résultats numériques sont en accord avec les obsentations de laboratoire et montrent que la mise en charge hydraulique d'un joint rugueux peut engendrer des instabilités lors de l'ouverture des aires de contact compatibles avec la micro-sismicité observée sur le site.

Non disponible / Not available; Actuellement l'état des connaissances sur l'hydrodynamique et la thermique des écoulements de suspensions newtoniennes et non newtoniennes ne permet pas de bien appréhender les problèmes liés à la stérilisation des écoulements chargés en particules dans l'industrie alimentaire. A ce jour, l'influence de la phase des particules solides sur l'écoulement n'est toujours pas très bien connue et pourtant c'est elle qui détermine la rhéologie du mélange et donc toute la Thermomécanique de ces écoulements. Les objectifs de ce travail expérimental sont de vérifier si les modèles rhéologiques établis pour des suspensions de particules de structure microscopique peuvent s'appliquer à des suspensions de "grosses" particules et d'expliquer les problèmes liés aux mécanismes de transfert de chaleur entre la suspension et la paroi. Les études menées, ont montré que le modèle de milieu effectif est efficace dans l'analyse de la perte de charge pour les fluides newtonien et non newtonien, en régimes laminaire et turbulent en situation isotherme. La migration latérale des particules vers l'axe de la conduite contribue à la modification de la rhéologie du mélange. La sédimentation observée et les interactions entre le fluide porteur et les particules sont à l'origine de la non homogénéité des suspensions en écoulement isotherme. Le transfert thermique est influencé par les particules. Le développement de la couche limite thermique est perturbé par la présence des particules. L'objectif est la maîtrise de la stérilisation en continu, des fluides alimentaires chargés de particules, par procédés thermiques. Le fluide porteur s'organise en anneau monophasique qui entoure la zone des particules au centre du tube. Cette organisation de l'écoulement chargé est néfaste puisqu'il fait émigrer les particules vers la zone froide où la stérilisation est difficile.

Non disponible / Not available, Actuellement l'état des connaissances sur l'hydrodynamique et la thermique des écoulements de suspensions newtoniennes et non newtoniennes ne permet pas de bien appréhender les problèmes liés à la stérilisation des écoulements chargés en particules dans l'industrie alimentaire. A ce jour, l'influence de la phase des particules solides sur l'écoulement n'est toujours pas très bien connue et pourtant c'est elle qui détermine la rhéologie du mélange et donc toute la Thermomécanique de ces écoulements. Les objectifs de ce travail expérimental sont de vérifier si les modèles rhéologiques établis pour des suspensions de particules de structure microscopique peuvent s'appliquer à des suspensions de "grosses" particules et d'expliquer les problèmes liés aux mécanismes de transfert de chaleur entre la suspension et la paroi. Les études menées, ont montré que le modèle de milieu effectif est efficace dans l'analyse de la perte de charge pour les fluides newtonien et non newtonien, en régimes laminaire et turbulent en situation isotherme. La migration latérale des particules vers l'axe de la conduite contribue à la modification de la rhéologie du mélange. La sédimentation observée et les interactions entre le fluide porteur et les particules sont à l'origine de la non homogénéité des suspensions en écoulement isotherme. Le transfert thermique est influencé par les particules. Le développement de la couche limite thermique est perturbé par la présence des particules. L'objectif est la maîtrise de la stérilisation en continu, des fluides alimentaires chargés de particules, par procédés thermiques. Le fluide porteur s'organise en anneau monophasique qui entoure la zone des particules au centre du tube. Cette organisation de l'écoulement chargé est néfaste puisqu'il fait émigrer les particules vers la zone froide où la stérilisation est difficile.

Numerical simulation of industrial reactor promises to be a precious tool that could help in prospecting new mixing technologies, in scaling-up and optimising reactors. Fast reactions, such as polymerisation and precipitations, are difficult to handle when mixing processes and reactions strongly interacts. This work aims at evaluating and validating reactive flow simulation for chemical engineering reactors using a special method which is able to take into account mixing reaction interaction for fast reactions. This method is based on solving Navier-Stokes equations in turbulent flow and then on solving the transport equation of the jointtcomposition Probability Density Function (PDF) with a Monte-Carlo solver. The validation is carried-out on a semi-batch Taylor-Couette precipitator. Each calculation step is validated. Finally, comparison between experiment and simulation proved that the method we chose has to be handle with special cares. Calculation has to be validated and 2 model constants have to be fitted and chemical kinetics known quite precisely. Besides, Monte-Carlo algorithm produces numerical noise which can highly affect reliability of the results. The calculation of the specific studied reactor didn'nt give completely good agreement with experiments. Nevertheless, the discrepancy is assigned not to the numerical approach but to the fact that the studied reactive flow is rather tricky to be calculated.; La simulation numérique de réacteurs industriels promet de devenir un outil précieux aidant à prospecter de nouvelles technologies de mélange, à dimensionner et à optimiser les conditions de fonctionnement du réacteur. Ces calculs deviennent difficiles dans le cas de réactions rapides (polymérisations anioniques, précipitations) où existe un fort couplage entre la réaction chimique et le mélange. Ce travail de thèse a pour but d'évaluer et de valider la simulation d'écoulements réactifs pour des applications de génie chimique avec une méthode que nous avons choisie pour sa capacité à tenir compte des interactions entre processus de mélange turbulent et réactions chimiques rapides. La méthode sélectionnée est une méthode basée sur la résolution numérique des équations de Navier-Stokes en milieu turbulent et l'équation de transport de la fonction de densité de probabilité jointe de la composition (PDF) à l'aide d'un algorithme de Monte-Carlo. L'évaluation s'est faite sur un réacteur de précipitation semi-batch de type Taylor-Couette. Chaque étape du calcul a été validée expérimentalement. Au terme de ce travail, il a été montré que l'emploi de ce type d'approche pour le calcul d'un réacteur industriel est hasardeux et délicat. Il demande une validation de l'écoulement, le calage de 2 constantes de modélisation et la connaissance précise des cinétiques de réaction. Enfin, l'utilisation d'un algorithme de type Monte-Carlo engendre un bruit numérique qui selon le type d'écoulement peut nuire à la fiabilité des résultats. Le cas traité n'a pas donné des résultats totalement conformes à l'expérience. La méthode numérique utilisée n'est pas pour autant remise en question car le cas traité était par nature délicat

Non disponible / Not available, La modélisation hydromécanique des milieux fracturés est très complexe. Beaucoup de simplifications sont nécessaires pour modéliser de tels milieux. Dans la plupart des cas, on ne tient compte que des fractures perméables. Le reste du réseau est approché par un milieu continu équivalent. Même si on suppose que cette approche est admissible, les propriétés hydrauliques et mécaniques des fractures et du milieu continu sont rarement connues. On ne connaît pas très bien la nature des mesures qui doivent être effectuées afin de mener à bien une modélisation en milieu discontinu. Afin de mieux comprendre les phénomènes hydromécaniques dans les mieux fracturés, on a travaillé sur deux sites différents. Le premier est le site du Grimsel en Suisse dans lequel un laboratoire souterrain est situé à environ 400 m de profondeur. L?expérience Febex consiste à l?étude in situ des phénomènes consécutifs à la mise en place d?une source de chaleur représentative de déchets radioactifs dans les 17 derniers mètres du tunnel Febex dans le laboratoire de Grimsel. Nous avons seulement modélisé la phase d?excavation de ce tunnel. La modélisation de Febex a permis d?établir une méthodologie de calage des propriétés hydrauliques dans les milieux discontinus. Cependant, ce genre d?étude sur des sites aussi complexes ne permets pas de répondre aux nombreuses questions qui se posent sur le comportement hydromécanique des milieux fracturés. Nous avons donc réalisé des modélisations sur un deuxième site plus accessible et de dimensions plus réduites. Le site expérimental de Coaraze, situé dans les Alpes maritimes, est un petit massif calcaire fracturé vanné au niveau d?une sources principale et muni d?extensomètres et de manomètres pour les mesures simultanées pression-déformation à différents niveaux du massif. La modélisation de Coaraze a permis de reproduire les mesures in situ et de donner un point de vue sur la problématique du couplage hydromécanique dans les milieux fracturés.

During the last decades, the study of coupling phenomena became more and more important in Geomechanics and Civil Engineering. Radioactive wastes repository in deep geological formation is a perfect illustration of this kind of problems. Indeed, due to numerous stakes, all acting phenomena must be considered in order to determine the first fonctionnality of the repository structure: the full safety confining of radioactive wastes. During the building period of the repository structures, the excavation and the ventilation of galleries induce solicitation of the surrounding rock mass. The abject of this work is to characterize the importance .and the nature of this solicitation. It's fundamental to determine the phenomena kinetics (progression of the drying front induced by ventilation, the extent of the perturbated zone) and the strengh of induced couplings (stresses in front of gallery, damage). The aim of all this considerations is to evaluate the eventual modifications of the rock and its confining abilities. The theoretical approach is the Mechanics of Porous Media defined by Coussy (1995), which leads to a formalism at the macroscopic scale and allows to determine the evolutions of the rock, considered as partially saturated porous medium (three-phases medium). The aim of the different models presented in this work is to describe the fluids transfers in the rock mass and to take into account constitutive aspects of the mechanical behavior of the rock (damage, plasticity)., Durant les dernières décennies, l'étude des phénomènes couplés a pris une importance considérable dans de nombreux problèmes de la Géomécanique et du Génie Civil. Le stockage des déchets radioactifs en formations géologiques profondes est une parfaite illustration de ce genre de problème. En effet du fait des nombreux enjeux, il va être nécessaire de prendre en compte tous les phénomènes intervenant, afin de pouvoir tirer des conclusions sur la fonctionnalité première de la structure de stockage à savoir: le confinement en toute sûreté des déchets radioactifs. Durant la phase de construction des structures de stockage, le creusement et la ventilation des galeries vont solliciter le massif rocheux environnant. L'objet de ce travail est de caractériser l'importance et la nature de cette sollicitation. Il est notamment primordial de déterminer la cinétique des phénomènes (progression du front de séchage induit par la ventilation, étendue de la zone perturbée) et la forces des couplages induits (contraintes en paroi des galeries, endommagement, etc.). Toutes ces considérations ont pour but d'évaluer les modifications éventuelles de la roche et ainsi ses capacités de confinement. L'approche théorique retenue est celle de la Mécanique des Milieux Poreux définie par Coussy (1995) permettant, grâce à un formalisme à l'échelle macroscopique, de déterminer les évolutions de la roche en tant que milieu poreux partiellement saturé (formé d'une phase solide, d'une phase liquide et d'une phase gazeuse). Les différents modèles présentés dans ce travail ont pour objectif de traduire au mieux les transferts de fluides au sein du massif rocheux mais également les aspects constitutifs liés au comportement mécanique de la roche (endommagement, plasticité).

Non disponible / Not available, La modélisation hydromécanique des milieux fracturés est très complexe. Beaucoup de simplifications sont nécessaires pour modéliser de tels milieux. Dans la plupart des cas, on ne tient compte que des fractures perméables. Le reste du réseau est approché par un milieu continu équivalent. Même si on suppose que cette approche est admissible, les propriétés hydrauliques et mécaniques des fractures et du milieu continu sont rarement connues. On ne connaît pas très bien la nature des mesures qui doivent être effectuées afin de mener à bien une modélisation en milieu discontinu. Afin de mieux comprendre les phénomènes hydromécaniques dans les mieux fracturés, on a travaillé sur deux sites différents. Le premier est le site du Grimsel en Suisse dans lequel un laboratoire souterrain est situé à environ 400 m de profondeur. L?expérience Febex consiste à l?étude in situ des phénomènes consécutifs à la mise en place d?une source de chaleur représentative de déchets radioactifs dans les 17 derniers mètres du tunnel Febex dans le laboratoire de Grimsel. Nous avons seulement modélisé la phase d?excavation de ce tunnel. La modélisation de Febex a permis d?établir une méthodologie de calage des propriétés hydrauliques dans les milieux discontinus. Cependant, ce genre d?étude sur des sites aussi complexes ne permets pas de répondre aux nombreuses questions qui se posent sur le comportement hydromécanique des milieux fracturés. Nous avons donc réalisé des modélisations sur un deuxième site plus accessible et de dimensions plus réduites. Le site expérimental de Coaraze, situé dans les Alpes maritimes, est un petit massif calcaire fracturé vanné au niveau d?une sources principale et muni d?extensomètres et de manomètres pour les mesures simultanées pression-déformation à différents niveaux du massif. La modélisation de Coaraze a permis de reproduire les mesures in situ et de donner un point de vue sur la problématique du couplage hydromécanique dans les milieux fracturés.

Numerical simulation of industrial reactor promises to be a precious tool that could help in prospecting new mixing technologies, in scaling-up and optimising reactors. Fast reactions, such as polymerisation and precipitations, are difficult to handle when mixing processes and reactions strongly interacts. This work aims at evaluating and validating reactive flow simulation for chemical engineering reactors using a special method which is able to take into account mixing reaction interaction for fast reactions. This method is based on solving Navier-Stokes equations in turbulent flow and then on solving the transport equation of the jointtcomposition Probability Density Function (PDF) with a Monte-Carlo solver. The validation is carried-out on a semi-batch Taylor-Couette precipitator. Each calculation step is validated. Finally, comparison between experiment and simulation proved that the method we chose has to be handle with special cares. Calculation has to be validated and 2 model constants have to be fitted and chemical kinetics known quite precisely. Besides, Monte-Carlo algorithm produces numerical noise which can highly affect reliability of the results. The calculation of the specific studied reactor didn'nt give completely good agreement with experiments. Nevertheless, the discrepancy is assigned not to the numerical approach but to the fact that the studied reactive flow is rather tricky to be calculated., La simulation numérique de réacteurs industriels promet de devenir un outil précieux aidant à prospecter de nouvelles technologies de mélange, à dimensionner et à optimiser les conditions de fonctionnement du réacteur. Ces calculs deviennent difficiles dans le cas de réactions rapides (polymérisations anioniques, précipitations) où existe un fort couplage entre la réaction chimique et le mélange. Ce travail de thèse a pour but d'évaluer et de valider la simulation d'écoulements réactifs pour des applications de génie chimique avec une méthode que nous avons choisie pour sa capacité à tenir compte des interactions entre processus de mélange turbulent et réactions chimiques rapides. La méthode sélectionnée est une méthode basée sur la résolution numérique des équations de Navier-Stokes en milieu turbulent et l'équation de transport de la fonction de densité de probabilité jointe de la composition (PDF) à l'aide d'un algorithme de Monte-Carlo. L'évaluation s'est faite sur un réacteur de précipitation semi-batch de type Taylor-Couette. Chaque étape du calcul a été validée expérimentalement. Au terme de ce travail, il a été montré que l'emploi de ce type d'approche pour le calcul d'un réacteur industriel est hasardeux et délicat. Il demande une validation de l'écoulement, le calage de 2 constantes de modélisation et la connaissance précise des cinétiques de réaction. Enfin, l'utilisation d'un algorithme de type Monte-Carlo engendre un bruit numérique qui selon le type d'écoulement peut nuire à la fiabilité des résultats. Le cas traité n'a pas donné des résultats totalement conformes à l'expérience. La méthode numérique utilisée n'est pas pour autant remise en question car le cas traité était par nature délicat

Non disponible / Not available, Les situations d'envasement (barrages, estuaires), la stabilité de fonds marins, l'évolution de bassins de décantation d'effluents issus de diverses activités humaines constituent des phénomènes ou des procédés fortement conditionnés par le comportement de sédiments fins de teneur en eau élevée. Situés à la frontière entre milieux fluide et solide, ces matériaux sont particulièrement complexes et évolutifs et nécessitent des études interdisciplinaires. La description et la prédiction de leur comportement reposent sur la détermination de lois caractéristiques ; l'expérimentation en laboratoire constitue donc une étape préalable indispensable à la compréhension des phénomènes gouvernant la formation de ces sols. Dans ce contexte, cette étude concerne plus particulièrement le début de la consolidation et repose sur deux approches complémentaires incluant la mise en oeuvre de dispositifs expérimentaux spécifiques : (i) banc gammamétrique et colonnes de sédimentation instrumentées permettant la simulation du processus de dépôt de suspensions de kaolin et la mesure non-destructive de certains paramètres ; (ii) système de type piézocône, adapté à la mesure en place de faibles valeurs de résistance au cisaillement non drainé. Le comportement des suspensions et des dépôts est analysé en terme d'évolution, au cours du temps, de la position de leur surface, de leur densité, de leur cohésion, des surpressions interstitielles, des contraintes effectives. Les résultats font l'objet d'une analyse géotechnique comparative entre différents modes de consolidation, d'une part pour mettre en évidence la continuité du processus sur un domaine étendu de contraintes et d'indice des vides, d'autre part afin d'étudier le comportement mécanique des dépôts obtenus en laboratoire parallèlement à celui observé in situ sur des sédiments marins récents. Nous abordons également la notion de " transition " sédimentation-consolidation, de manière classique et sous l'angle de l'acquisition de résistance.

Non disponible / Not available, Les situations d'envasement (barrages, estuaires), la stabilité de fonds marins, l'évolution de bassins de décantation d'effluents issus de diverses activités humaines constituent des phénomènes ou des procédés fortement conditionnés par le comportement de sédiments fins de teneur en eau élevée. Situés à la frontière entre milieux fluide et solide, ces matériaux sont particulièrement complexes et évolutifs et nécessitent des études interdisciplinaires. La description et la prédiction de leur comportement reposent sur la détermination de lois caractéristiques ; l'expérimentation en laboratoire constitue donc une étape préalable indispensable à la compréhension des phénomènes gouvernant la formation de ces sols. Dans ce contexte, cette étude concerne plus particulièrement le début de la consolidation et repose sur deux approches complémentaires incluant la mise en oeuvre de dispositifs expérimentaux spécifiques : (i) banc gammamétrique et colonnes de sédimentation instrumentées permettant la simulation du processus de dépôt de suspensions de kaolin et la mesure non-destructive de certains paramètres ; (ii) système de type piézocône, adapté à la mesure en place de faibles valeurs de résistance au cisaillement non drainé. Le comportement des suspensions et des dépôts est analysé en terme d'évolution, au cours du temps, de la position de leur surface, de leur densité, de leur cohésion, des surpressions interstitielles, des contraintes effectives. Les résultats font l'objet d'une analyse géotechnique comparative entre différents modes de consolidation, d'une part pour mettre en évidence la continuité du processus sur un domaine étendu de contraintes et d'indice des vides, d'autre part afin d'étudier le comportement mécanique des dépôts obtenus en laboratoire parallèlement à celui observé in situ sur des sédiments marins récents. Nous abordons également la notion de " transition " sédimentation-consolidation, de manière classique et sous l'angle de l'acquisition de résistance.