Cette thèse a pour but le développement d'un modèle constitutif hydromécanique prenant en compte le couplage de l'endommagement et de la plasticité dans les géomatériaux non-saturés. Le couplage de ce modèle avec les transferts hydrauliques est rendu possible par son implémentation dans le code aux éléments finis $Theta$-Stock. Un algorithme local a été développé spécifiquement dans ce but. Des problèmes hydromécaniques complètement couplés, tels que la création de la zone endommagée par excavation autour d'un tunnel ou la création d'un endommagement de dessiccation et d'humidification ont été simulés. Une contrainte doublement effective incorporant les effets de la succion et de l'endommagement (considéré isotrope) a été définie en s'appuyant sur des bases thermodynamiques. Cette approche a l'avantage de considérer une unique variable de contrainte étant thermodynamiquement conjuguée aux déformations élastiques. Une formulation hyperélastique dépendante de la pression de confinement est utilisée pour décrire le comportement dans le domaine élastique. L'évolution des rigidités élastiques avec l'endommagement est comparée en considérant deux hypothèses : le principe des déformations équivalentes et le principe de l'énergie élastique équivalente. L'hypothèse d'équivalence des déformations permet d'introduire la contrainte doublement effective dans les équations de plasticité et ainsi de coupler plasticité et endommagement. Les équations de plasticité sont dérivées du modèle de Barcelone (Alonso extit{et al.}, 1990), lui même basé sur le modèle de Cam-Clay. Deux critères distincts d'endommagement et de plasticité sont définis qui peuvent être activés aussi bien indépendamment que simultanément. Les surfaces de charges étant exprimées en fonction de la contrainte effective et de la succion, elles évoluent automatiquement en fonction de la succion et de l'endommagement dans l'espace des contraintes totales. Cela permet de représenter la transition d'un comportement ductile vers un comportement fragile lorsque le matériau est asséché. Un algorithme local explicite a été développé pour gérer le couplage des deux phénomènes dissipatifs. L'implémentation du modèle constitutif dans le code aux éléments finis $Theta$-Stock permet l'étude de problèmes hydromécaniques complètement couplés, les lois de transfert dépendant elles aussi de l'état de saturation du matériau. Ce modèle a d'abord été appliqué à la simulation du développement de micro-fissures lors du séchage d'une éprouvette. L'apparition de l'endommagement est expliquée principalement par le gradient de pression très important créé à la surface de l'échantillon lorsqu'on applique une forte variation de succion. Enfin, un problème à échelle réelle est simulé. L'excavation d'un tunnel, la désaturation du sol environnant dû à la ventilation, ainsi que la création de la zone endommagée par excavation sont étudiés. L'étendue des zones l'endommagement et de déformations plastiques autour du tunnel est étudiée / This thesis work objective is the development of a hydro-mechanical constitutive model which accounts for damage-plasticity couplings in unsaturated geomaterials. The coupling of this model with hydraulic transfers is made possible by its implementation into the Finite Element code Theta-Stock. In order to achieve this implementation, a specific stress-point algorithm has been developed. Fully coupled hydro-mechanical problems have been simulated, such as the creation of the Excavation Damaged Zone around a tunnel and the initiation of damage due to desiccation and humidification. A double effective stress incorporating both the effect of suction and damage (assumed isotropic) is defined based on thermodynamical considerations. The advantage of this approach is that it results in a unique stress variable being thermodynamically conjugated to elastic strains. A pressure-dependent hyperelastic formulation is used to describe the behaviour inside the elastic domain. The evolution of elastic rigidities with damage is then studied. Two hypotheses are compared, the principle of strain equivalence and the principle of equivalent energy. Coupling between damage and plasticity phenomena is achieved by following the principle of strain equivalence and incorporating the effective stress into plasticity equations. The plasticity framework is based on the Barcelona Basic Model (Alonso extit{et al.}, 1990), itself based on the Cam-Clay model. Two distinct criteria are defined for damage and plasticity, which can be activated either independently or simultaneously. Their formulation in terms of effective stress and suction allows them to evolve in the total stress space with suction and damage changes. This leads to a direct coupling between damage and plasticity and allows the model to capture the ductile/brittle behaviour transition occurring when clays are drying. A specific explicit algorithm has been developed to handle the association of the two dissipative phenomena. The implementation of the constitutive model into the Finite Element code Theta-Stock allowed for the simulation of fully coupled hydro-mechanical problems. The hydraulic transfer laws also consider the saturation state. This fully coupled model is first applied to simulate the development of micro-cracks during desiccation of a soil sample. Damage initiation is explained mainly by the important pressure gradient appearing at the boundary when applying a high suction change. Finally, a full-scale problem is simulated. The excavation of a tunnel, the desaturation of the surrounding soil and the creation of the excavation damaged zone are studied

La dispersion dans des milieux poreux homogènes (empilements de grains) a été étudiée par des mesures par résonance magnétique nucléaire (RMN) et des simulations de marches aléatoires dans un réseau de pores. La RMN permet de mesurer l’ensemble des déplacements des molécules d’eau durant un temps tΔ, et d’obtenir propagateurs et moments caractéristiques. L’évolution temporelle du second moment σ (σ2 ∝ taΔ) permet de caractériser de manière précise le régime de dispersion des molécules (Gaussien ou anormal). Des mesures pour des écoulements de 15 < Pe < 45 dans un empilement de grains de 30μm ont permis d’observer une dispersion anormale faiblement super-dispersive (a = 1.17) en écoulement saturé et une augmentation progressive du caractère super-dispersif avec la diminution de la saturation en eau (jusqu’à a = 1.5 pour 42 %)lors d’une co-injection stationnaire eau-huile. En écoulement saturé, les propagateurs et courbes de percée sont quasi-gaussiennes, tandis qu’en écoulement insaturé, les propagateurs sont asymétriques et les courbes de percée présentent des trainées aux grands temps. Dans ces conditions, on montre que la dispersion anormale observée est mieux décrite par des lois stables de Lévy que par des lois gaussiennes. Des simulations de marche aléatoire ont été réalisées dans un réseau de pores extrait d’un milieu poreux réel par imagerie microscanner.Elles permettent d’obtenir les mêmes informations que la RMN, les marcheurs se déplaçant par advection et diffusion. Ces simulations montrent l’existence d’une stagnation non observée dans les expériences, montrant que la simplification du réseau poreux est trop importante et empêche de reproduire certains aspects du champ de vitesses détecté par la RMN. Toutefois, l’évolution temporelle du second moment a également un caractère super-dispersif à temps long à 100 % de saturation / We investigated dispersion in homogeneous porous media (grain packs) by nuclear magnetic resonance (NMR) measurements and random walk simulations in pore networks. We measured water molecules displacements during a time interval tΔ by NMR measurements, which allows us to obtain propagators and charateristic cumulants of displacements such as the mean square displacement σ. The evolution of the cumulant σ as a function of time tΔ (σ2 ∝ taΔ) is a very sensitive test of Gaussian behaviour compared to the analysis of the shape of propagators. In a homogeneous 30μm grain pack and low Peclet numbers (15 < Pe < 45), we observed weak super dispersion in saturated conditions (a = 1.17) and gradually stronger super-dispersionas the water saturation decreases (up to a = 1.5 for 42 %) during steady-state oil-water two phase flow. Insaturated conditions, propagators and breakthrough curves are Gaussian or nearly Gaussian, whereas in two phase conditions, propagators are non symmetric and breakthrough curves show thick tails at long time. Weshow that the anomalous dispersion observed is better explained by Lévy stable laws (asymetric for longitudina ldispersion, and symetric for transverse dispersion) than by Gaussian laws. Random walk simulations were performed in a pore network constructed using high resolution images of the grain pack. They allow us to obtain the same informations than the NMR, with walkers submitted to diffusive and advective effects. The simulations show the existence of an anomalous stagnation not observed in experiments, highlighting the oversimplification of the pore network that prevent reproducing some aspects of the velocity field detected by NMR. However, the simulations indicate similarly a super-dispersion at long time in saturated conditions

L'objectif de ce travail est d'obtenir, par homogénéisation périodique, des modèles macroscopiques de transfert hydrique dans les milieux poreux partiellement saturés à partir des équations de transfert de l'eau liquide et de vapeur d'eau écrites à une échelle microscopique. La dimensionnalisation des équations fait apparaître naturellement des nombres sans dimension caractérisant les problèmes de transfert hydriques dans les milieux partiellement saturés. Nous nous sommes intéressés à trois différents régimes de transfert (diffusion de vapeur prédominante, couplage diffusion/convection, convection de l'eau liquide prédominante). Pour chaque modèle homogénéisé, nous avons obtenu une expression différente du tenseur de diffusion hydrique homogénéisé. Nous avons ensuite calculé les tenseurs de diffusion hydrique homogénéisés obtenus dans les deux régions hygroscopique et super-hygroscopique, sur des géométries plus ou moins complexes décrivant la microstructure en 2D et 3D. Des comparaisons avec des valeurs expérimentales ont été ensuite effectuées. Pour finir, une résolution numérique de l'équation de transfert hydrique macroscopique homogénéisée a été effectuée en se basant sur les données expérimentales d'un béton BHP. / We propose in this work to construct, by periodic homogenization, macroscopic models of moisture transfer in unsaturated porous media. To do this, the liquid water and water vapor transport equations are averaged from the microscopic scale. The dimensional analysis of transport equations naturally lets appear dimensionless numbers characterizing the moisture transfer in unsaturated porous media. Three different transfer regimes are addressed (predominant water vapor diffusion, coupling diffusion / convection, predominant liquid water convection). For each transfer regime, the associated homogenized moisture diffusion tensor has a different expression. Then, the homogenized moisture diffusion tensors are calculated in both hygroscopic and super-hygroscopic regions on several geometries with varying complexity, describing 2D and 3D microstructures. Comparisons with experimental values are also addressed. Finally, based on experimental data of a BHP concrete, a numerical resolution of the homogenized macroscopic moisture transfer equation is performed.

On s'intéresse plus particulièrement à la modélisation du comportement d'un matériau hétérogène méso-fissuré (béton, roche,...), soumis à une sollicitation thermo-hydro-mécanique avec prise en compte du couplage géométrique. Pour conduire cette étude, on s'appuie notamment sur les approches micro-mécaniques du comportement des milieux méso-fissurés non saturés développées depuis quelques années au Laboratoire des Matériaux et des Structures du Génie Civil - Ur Navier - Université Paris Est. Le milieu fissuré non saturé traité ici est constitué d'une matrice solide homogène élastique linéaire et de fissures connectées saturées par deux fluides immiscibles : un liquide et un gaz séparés par une surface capillaire. La fissure est traditionnellement considérée comme une cavité ellipsoïdale (cas 3D) ou elliptique (cas 2D) dont le rapport d'aspect tend vers zéro. Deux morphologies typiques de matériau sont considérés dans ce travail : la situation où les fissures sont toutes orientées dans la même direction et la situation où les fissures possèdent des orientations aléatoires. Dans une première étape, on rappelle brièvement les résultats disponibles concernant la modélisation des fissures non saturées par des cavités ellipsoïdales aplaties. A la fin de cette première partie, on complète les résultats déjà disponibles en étudiant l'influence de l'histoire de chargement sur la réponse de matériau. Dans une deuxième étape, on s'attache à valider une partie des résultats obtenus en utilisant une description des efforts capillaires dans les fissures par une précontrainte homogène en seréférant aux solutions analytiques exactes disponibles dans la littérature permettant de décrire le comportement d'une fissure isolée au sein d'une matrice élastique. Dans une troisième étape, on s'intéresse aux phénomènes de propagation des fissures en condition non saturée. Les lois de propagation sous critique et le phénomène de branchement des fissures sont également prises en compte dans cette approche. La dernière partie de la thèse concerne l'influence de la température sur le comportement des milieux poreux non saturés / The main topic of my work is the development of a micromechanical model for the behaviour of unsaturated mesocracks in media (concrete, rock...) in which the thermo-hydro-mechanical loadingsand thermo-hydro-mechanical couplings are taken into account. For this, we used the micromechanical approach model of behaviour of cracked porous media recently developed at LMSGC. My thesis is focused on the equilibrium configurations of a porous material whose pore space is saturated by a vapour and a liquid phase. The behaviour of an elastic medium containing unsaturated mesocracks is studied in the framework of a micromechanical approach. The cracks are filled by two immiscible fluids, namely a liquid and a gas, separated by a capillary interface. Furthermore, it is assumed that the set of cracks constitutes a connected network ; the capillary pressure is uniform over a representative elementary volume. The cracks are modelled as flat oblate spheroid cavities. Several geometrical configurations of cracks in porous media are considered in the framework of Eshelby-based homogenization methods (parallel cracks, randomly oriented cracks). First, a previously developed model showed that when coupling between the deformation of the cracks and the capillary forces is taken into account, there is no more a one-to-one relationship between the loading parameters and the state-variables. Thus, we describe the loading history prescribed to the material in order to compute its response. Second, we validate these results referring to the exact solutions available in the literature to describe the behaviour of a unsaturated crack within an elastic matrix. Third, the description of crack propagation in unsaturated media is considered in the framework of linear elastic fracture mechanics. The phenomenon of subcritical crack growth due to stress corrosion cracking is taken into account in this approach. Mixed mode fracture in the plane is also examined. Finally, we are interested in the influence of the temperature on the behavior of unsaturated porous media in the framework of the micromechanical approach

Au cours de la dernière décennie, la méthode d’analyse isogéométrique (AIG) a attiré l’attention des chercheurs grâce à ses capacités supérieures à la méthode standard des éléments finis (MEF). Le concept AIG utilise les mêmes fonctions de base que celles utilisées dans la conception assistée par ordinateur (CAO) pour l’approximation des champs inconnus tels que les déplacements, pression interstitielle ou la température dans la solution des éléments finis d’un problème thermo–hydro–mécanique (éventuellement couplé). Parmi les caractéristiques les plus importantes d’AIG, la régularité, le taux de convergence et surtout sa continuité intrinsèque d’ordre supérieur représentent une nette amélioration par rapport à la méthode standard des éléments finis, permettant d’obtenir des avantages computationnels significatifs en termes de précision de la solution et de efficacité.Ce travail tente d’exploiter les caractéristiques d’AIG pour la résolution numérique des problèmes hydromécaniques (HM) couplés dans les géomatériaux de second gradient de type poro–élastoplastiques partiellement saturés. D’une part, le modèle second gradient appartenant à la théorie des milieux continus avec microstructure assure l’objectivité des résultats en présence de phénomènes de localisation de la déformation en termes d’indépendance de maillage de la solution numérique, ce qui ne peut être réalisé avec des modèles constitutifs classiques qui n’implique pas l’intervention d’une longueur interne. D’autre part, la continuité C1 réalisable au moyen de fonctions de base AIG permet une implémentation directe de tels modèles constitutifs d’ordre supérieur, dans une formulation HM dérivée de l’approche de mélange classique. De plus, la régularité des fonctions de base AIG s’est révélée très efficace dans la modélisation de processus couplés caractérisés par de forts gradients hydrauliques – comme la simulation de la propagation d’un front de saturation dans une pente partiellement saturée. Dernier point, mais non des moindres, il convient de noter que, par rapport aux approches existantes basées sur les multiplicateurs de Lagrange, la méthode AIG pour résoudre les problèmes hydromécaniques (HM) couplés dans les matériaux du second gradient saturé et partiellement saturé permet une réduction considérable du nombre de degrés de libertés requis pour atteindre le même niveau de précision. Cela entraîne non seulement une augmentation significative de l’efficacité de calcul, mais permet également d’étendre la formulation du second gradient à l’analyse de problèmes réalistes en 3D, dont la solution a été présentée pour la première fois dans ce travail.La formulation poro–élastoplastique du second gradient développée dans ce travail est mise en œuvre dans le code orienté vers la recherche GeoPDEs, un code IAG–MEF open source écrit en Matlab et développé à l’Université de Pavia. Sur la base des résultats obtenus dans une large série de problèmes aux limites en 2D et 3D analysées dans ce travail, on peut conclure que la combinaison de AIG et d’élastoplasticité du second gradient représente un outil puissant pour la simulation numérique de problèmes géotechniques caractérisés par de forts couplages multiphysiques, un comportement fortement non linéaire du sol, et des gradients de déplacement et de pression interstitielle fortement localisés. / During the last decade, Isogeometric Analysis (IGA) has drawn the attention of the Finite Element community to its superior capabilities over the standard Finite Element Method (FEM). The IGA concept uses the same basis functions used in Computed Aided Design (CAD) for the approximation of the unknown fields such as displacements, pore pressure or temperature in the Finite Element solution of a (possibly coupled) thermo– hydro–mechanical problem. Among the most relevant features of IGA, its smoothness, its convergence rate and particularly its intrinsic higher–order continuity between elements represent a definite improvement over the standard FEM, which allow to obtain significant computational advantages in terms of accuracy of the solution and computa- tional efficiency.This work attempts to exploit the characteristics of IGA for the numerical solution of coupled hydro–mechanical (HM) problems in saturated and partially saturated second gradient poro–elastoplastic geomaterials. On one hand, the second gradient model belonging to the theory of continua with microstructure ensures the objectivity of the results in presence of strain localization phenomena in terms of mesh independence of the numerical solution, which cannot be achieved with classical constitutive models without an internal length scale. On the other hand, the C1–continuity achievable by means of IGA basis functions allows a straightforward implementation of such higher order constitutive models, within a HM formulation derived from the classical mixture approach. In addition, the smoothness of the IGA basis functions proved to be very efficient in the modeling of coupled processes characterized by strong hydraulic gradients – such as the simulation of the downward propagation of a saturation front in a partially saturated slope subject to rainfall infiltration. Last but not least, it is worth noting that, as compared to the existing approaches based on Lagrange multipliers, the IGA approach to the solution of coupled hydro-mechanical (HM) problems in saturated and partially saturated second gradient materials allows a dramatic reduction in the number of degrees of freedoms required to achieve the same level of accuracy. This not only results in a significant increase of the computational efficiency, but also allows to extend the complete second gradient formulation to the analysis of realistic 3D problems, the solution of which has been presented in this work for the first time.The local second gradient poro–elastoplastic formulation developed in this work is implemented in the research-oriented code GeoPDEs, a Matlab open source IGA–FEM code developed at the University of Pavia. Based on the results obtained in a large series of representative 2D and 3D initial–boundary value problems analyzed in this work, it can be concluded that the combination of IGA and the second gradient elastoplasticity represents a powerful tool for the numerical simulation of geotechnical problems characterized by strong multiphysics couplings, highly nonlinear behavior of the soil, and strongly localized displacement and pore pressure gradients.

A l’heure actuelle, des sédiments fluviaux et lacustres se trouvent pollués par les polychlorobiphényles (PCB), à des degrés de gravité divers. Les PCB se présentent comme des micropolluants persistants dans l’environnement et nuisibles à la santé. Un programme de recherche « PCB-Axelera » est conduit afin de pencher sur cette problématique et explorer les solutions de diminution de cette pollution. En France, l'ensemble des arrêtés préfectoraux (30/10/2002 et 14/06/2000) et le décret du 18/04/2002 considèrent les sédiments pollués comme déchets et autorisent leur stockage dans les installations de stockage de déchets si leurs teneurs en PCB sont comprises entre 0,68 ppm et 50 ppm. L’objectif principal de cette étude est de concevoir une méthodologie adaptée pour une installation de déchets non-dangereux ISDND. Les géomatériaux dans ce stockage seront utilisés comme barrière d’étanchéité pour les sédiments pollués. L’étude est composée de trois parties : Dans la première partie la caractérisation hydromécanique chimique des sédiments pollués par les PCB est présentée. La classification des sédiments est effectuée. Le comportement mécanique du sédiment est déterminé à l’aide de l’œdomètre. Quant au comportement hydromécanique il est étudié à l’aide du filtre presse et du perméamètre à paroi rigide (PPR). Cette caractérisation est indispensable pour proposer la meilleure mise en œuvre des sédiments dans le centre de stockage. Dans la deuxième partie, les performances hydraulique (perméabilité) et mécanique (gonflement) des différents géomatériaux en contact avec un fluide extrait du sédiment et d’autres fluides synthétiques contenant des PCB sont explorées. L’étude du gonflement des géomatériaux est réalisée à l’échelle du laboratoire tandis que, l’étude de sa performance hydrique est réalisée sur deux échelles centimétrique et métrique. L’étude menée sur un pilote de 1,2 m3 de volume comprenant un système de pluviation automatique a permis de suivre les tassements à l’aide d’un capteur ultrason et de mesurer la perméabilité des sédiments et des géomatériaux sous conditions réelles de stockage. Dans la troisième partie, un modèle hydraulique 1D est développé afin de suivre le front d’eau au sein d’un échantillon de sédiment soumis à une charge hydraulique représentant l’expérimentation dans un PPR. La méthode VER est utilisée. Les résultats expérimentaux montrent que le degré de contamination en PCB de ces sédiments n’influence pas le comportement hydromécanique et donc la qualité d’étanchéité des géomatériaux étudiés. Les résultats numériques de simulation du comportement hydraulique du sédiment sont en concordance avec les résultats expérimentaux. La solution de stockage hors site des sédiments pollués est alors envisageable puisqu’elle confine la pollution en éliminant le risque de contamination du milieu environnant. / Currently, river and lake sediments are contaminated with polychlorinated biphenyls (PCBs), to varying degrees of severity. PCBs are persistent micropollutants in the environment and harmful to health. A research program "PCB-Axelera" is conducted to address this issue and explore solutions to eliminate or decrease the pollution. In France, all prefectural orders (30/10/2002 and 14/06/2000) and the Decree of 18/04/2002 considered contaminated sediments as waste and allow their landfilling in waste storage facility if their PCB concentrations are between 0.68 ppm and 50 ppm. The main objective of this study is to develop a methodology suitable for installation of non-hazardous waste landfill sites. Geomaterials in this storage will be used as a barrier for confinement of contaminated sediment. The study consists of three parts: 1. In the first part, hydro-chemical characterization of sediments contaminated with PCBs is presented. The sediment classification is performed. The mechanical behavior of the sediment is determined using oedometer tests. As for the hydro-mechanical behavior, it is studied using the filter press and the rigid wall permeameter (PPR). This characterization is essential to provide the best implementation of sediment in the storage center. 2. In the second part, the hydraulic performance (permeability) and mechanical (swelling) of various geomaterials in contact with a fluid extracted from the sediment and other synthetic fluids containing PCBs are being explored. The study of swelling of geomaterials is performed on a laboratory scale while studying its water performance was achieved on two centimetric and metric scales. The pilot scale study was conducted on a 1.2 m3 volume tank. The pilot test was equipped with an automatic rain system and an ultrasonic sensor to measure to measure the follow the settlements of the sediment. The permeability of sediments and geomaterials were measured using the pilot test under actual storage conditions . 3. In the third part, a 1D hydraulic model is developed to follow the waterfront in a sediment sample subjected to a hydraulic load representing a PPR experiment. The VER method is used. A water model of the sediment is proposed to reproduce the experimental results. The water retention curve of the sediment is calculated using Durner’s model (1994). The experimental results show that the degree of PCB contamination of the sediments does not affect the hydro-mechanical behavior and therefore the hydro-mechanical performance of the studied geomaterials studied. Numerical simulation results of the hydraulic behavior of sediment are consistent with the experimental results. The solution of landfilling the PCB contaminated sediment is possible since it assure the confinement of the pollution and eliminates the risk of environmental contamination.

[pt] Este trabalho apresenta o estágio atual do desenvolvimento de uma ferramenta numérica para simulação e análise do fluxo bifásico e do transporte multicomponente e de calor em meios porosos com interação química em duas dimensões. O programa desenvolvido foi baseado no UNSATCHEM_2D escrito em linguagem Fortran. A discretização das equações diferenciais parciais é feita em elementos finitos no espaço e em diferenças finitas no tempo. Uma interface grá- fica foi criada em linguagem Lua para a entrada dos dados e visualização dos resultados. Esta interface foi baseada no programa Mtool do TecGraf. O objetivo do trabalho foi implementar uma ferramenta que acoplasse diferentes processos em meios porosos: fluxo bifásico, transporte de solutos, processos químicos e temperatura. A ferramenta desenvolvida deverá ser útil em análises de problemas ambientais. Em particular, o enfoque deste trabalho é no acoplamento de processos químicos. Para isto, uma rotina de especiação química foi acoplada à cinética química com relevância na oxidação mineral. Para resolver o problema da especiação química, o programa Phreeqc foi utilizado. O trabalho apresenta detalhes dos processos de análise, de suas equações relevantes e da implementação numérica. No estágio atual de desenvolvimento, o programa é capaz de avaliar o transporte de metais e, em particular, a simulação da drenagem ácida em barragens de rejeito e pilhas de estéril sob condições complexas. Exemplos de validação e ilustrativos são mostrados e discutidos. São sugeridas algumas possibilidades de futuros acoplamentos à ferramenta numérica desenvolvida. / [en] This work presents the present stage of development of a numerical platform for simulation and analysis of two-phase flow and of, non isotermic, multicomponent transport in porous media with chemical interaction in two dimensions. The developed computer code is based on the UNSATCHEM_2D code written in Fortran language. The discretization of the partial differential equations is done by finite elements at the space and in finite difference at the time. A graphical interface was created in language Lua for the input of the data and result visualization. This interface was based on the program TecGraf’s Mtool. The purpose of the work was to construct a numerical tool that couples different processes in porous media: flow (two-phase), transport of solutes, chemical processes and temperature. The developed tool should be useful in the analysis of environmental problems. In particular, the work focused on the coupling of chemical processes. For this, a chemical speciation routine was coupled to chemical kinetics with relevance to mineral oxidation. In order to solve the chemical speciation problem, the program Phreeqc was used. The work presents details of the analyzed processes, their relevant equations and of the numerical implementation. At the present stage of development, the computer program is able to evaluate the transport of metals and, at particular, the simulation of the acid mine drainage at tailings discharge and waste rock piles under complex conditions. Validation and illustrative examples are shown and discussed. Possibilities of future additions to the developed tool are suggested.