[pt] O comportamento sísmico de taludes tem sido um tópico de grande interesse da engenharia geotécnica nos últimos 40 anos. Durante este período, a prática da engenharia nesta área evoluiu do emprego de técnicas elementares para procedimentos numéricos bastante complexos. A abordagem mais simples é a análise pseudo-estática na qual o carregamento do terremoto é simulado por uma aceleração horizontal estática equivalente atuando na massa de solo deslizante, utilizando-se um procedimento de equilíbrio limite (método das fatias), geralmente conservativo. O parâmetro que descreve o comportamento dinâmico do solo é referido como coeficiente sísmico k, e sua seleção depende fortemente da experiência e normas técnicas locais, porque não há maneira simples e segura de se escolher um valor adequado. O segundo procedimento é conhecido como método de Newmark, que envolve o cálculo de uma aceleração de escoamento, definida como a força inercial necessária para o fator de segurança atingir 1 em uma análise pseudo-estática pelo método de equilíbrio limite. O procedimento então usa os registros de aceleração do terremoto de projeto e o integra duplamente no tempo para calcular os deslocamentos permanentes acumulados. O terceiro método é referido como análise de Makdisi- Seed, que procura definir a estabilidade sísmica do talude em termos de deslocamentos aceitáveis em vez de um fator de segurança tradicional através de uma versão modificada do método de Newmark. Esta técnica apresenta uma maneira racional de calcular uma aceleração de escoamento média, necessária para produzir um valor do coeficiente de segurança do talude igual a 1. Gráficos específicos foram também desenvolvidos para estimativa dos deslocamentos permanentes, tendo sido bastante aplicados em aterros rodoviários, barragens e aterros sanitários. Finalmente, o mais sofisticado método para análise de estabilidade sísmica de taludes é conhecido como análise dinâmica, que normalmente incorpora modelos de elementos finitos e relações tensão x deformação complexas numa tentativa de obter melhores representações para o comportamento mecânico de taludes sob cargas cíclicas Os resultados destas análises podem incluir a história no tempo dos deslocamentos e tensões, bem como das freqüências naturais, efeitos de amortecimento, etc. Este trabalho apresenta uma comparação entre os métodos mencionados anteriormente, analisando o comportamento sísmico dos taludes da estrutura de contenção dos resíduos de lixiviação de minério de urânio, na Bahia, e dos taludes do bota-fora sul da mina de cobre Toquepala, situada no Peru. / [en] The seismic stability of slopes has been a topic of considerable interest in geotechnical engineering for the past 40 years. During that period, the state of practice has moved from simples techniques to more complicated numerical procedures. The simplest approach is the pseudo-static analysis in which the earthquake load is simulated by an equivalent static horizontal acceleration acting on the mass of the landslide, according to a generally conservative limit equilibrium analysis. The ground motion parameter used in a pseudo-static analysis is referred to as the seismic coefficient k, and its selection has relied heavily on engineering judgment and local code requirements because there is no simple method for determining an appropriate value. The second main procedure is known as the Newmark displacement analysis which involves the calculation of the yield acceleration, defined as the inertial force required to cause the static factor of safety to reach 1 from the traditional limit equilibrium slope stability analysis. The procedure then uses a design earthquake strong-motion record which is numerically integrated twice for the amplitude of the acceleration above the yield acceleration to calculate the cumulative displacements. These displacements are then evaluated in light of the slope material properties and the requirements of the proposed development. The third method is referred to as the Makdisi-Seed analysis sought to define seismic embankment stability in terms of acceptable deformation instead of conventional factors of safety, using a modified Newmark analysis. Their method presents a rational means to determine yield acceleration, or the average acceleration required to produce a factor of safety of unity. Design curves were developed to estimate the permanent earthquake- induced deformations of embankments, which have since been applied to sanitary landfill and highway embankments. Finally, the most sophisticated method for seismic slope stability calculations is known as the dynamic analysis, which normally incorporates a finite element model and a rather complex stress-strain behavior for geological materials in an attempt to obtain a better representation of the behavior of soils under cyclic loading. The results of the analysis can include a time history of displacements and stresses, as well as natural frequencies, effects of damping, etc. This work presents a comparison of the results obtained by the aforementioned approaches, considering the seismic behavior of the slopes of an uranium lixiviation pad situated in Bahia, Brazil, and the South embankment of the waste landfill of the Toquepala Mine, Peru.

[pt] O presente trabalho foca em uma investigação numérica do escoamento através de meios porosos com macroporos, como os carbonatos. A presença de vugs e fraturas tem um forte efeito na característica do escoamento. O fluxo através da matriz porosa é geralmente descrito pela equação de Darcy e o fluxo através dos macroporos pela equação de Stokes. O acoplamento entre essas duas abordagens distintas traz grande complexidade para a modelagem de tais fluxos. Neste trabalho, utilizamos a formulação de Brinkman que é capaz de descrever ambos os fluxos através da matriz porosa e macroporos com uma única equação diferencial. Resolvemos o sistema de equações diferenciais usando o modelo de elementos finitos e implementamos o código com auxílio da plataforma FEniCS. Primeiro, resolvemos o fluxo 1-D através de placas paralelas com uma das paredes sendo um material poroso. O objetivo foi comparar as previsões obtidas pela formulação de Brinkman com aquela obtida usando a condição de contorno de Beavers-Joseph. Então, resolvemos um fluxo 2-D através de um meio poroso com macroporos. A geometria da estrutura porosa foi obtida a partir de cortes 2D de imagens tomográficas de carbonatos. O objetivo desta análise foi avaliar uma permeabilidade equivalente em função da área e estrutura dos macroporos. / [en] The present work focuses on a numerical investigation of the flow through porous media with macropores, such as carbonates. The presence of vugs and fractures have a strong effect on the flow characteristics. The flow through the porous matrix is usually described by the Darcy equation and the flow through the macropores by Stokes equation. The coupling between these two distinct approaches brings great complexity to the modeling of such flows. In this work, we use Brinkman formulation that is able to describe the flow, both in the porous matrix and macropores with a single differential equation. We solved the set of differential equations using the finite element model and implemented the code in the FEniCS platform. We first solved the 1-D flow through parallel plates with one of the walls being a porous material. The goal was to compare the predictions obtained with the Brinkman formulation to that obtained by using the Beavers-Joseph boundary condition. Then, we solved a 2-D flow through a porous medium with macropores. The geometry of the pore structure was obtained from 2D slices of tomographic images of carbonates. The goal of this analysis was to evaluate an equivalent permeability as a function of macropores area and structures.

[pt] Nos últimos anos, os materiais compósitos cimentícios vêm ganhando destaque na indústria da construção civil. Suas excelentes propriedades mecânicas e contribuição para o controle de propagação de fissuras são um atrativo para seu emprego como material de construção. No entanto, normas técnicas para projeto envolvendo estes materiais e estruturas ainda não são consagradas. Uma melhor compreensão do comportamento de materiais cimentícios com adição de fibras requer o estudo de suas fases e da interação entre elas. Análises em diferentes escalas possibilitam esta representação. Tensões e deformações, dano e iniciação de fissuras ocorrem na escala das heterogeneidades e ajudam a explicar e prever o comportamento do concreto em uma escala macroscópica. A modelagem e simulação do comportamento destes compósitos é complexa e desafiadora. Para tal, é necessário definir os principais mecanismos que descrevem o comportamento do material de modo a escolher a descrição matemática adequada. Esta dissertação propõe metodologias para a modelagem numérica multiescala de materiais compósitos cimentícios. A partir de informações obtidas na escala do material, busca-se compreender melhor o comportamento global do compósito. Para isto, serão desenvolvidos métodos numéricos e computacionais baseados no Método dos Elementos Finitos, em técnicas de Inteligência Artificial e nos conceitos da Mecânica do Dano Computacional. Na macroescala, um modelo contínuo equivalente é desenvolvido através de técnicas probabilísticas e de Inteligência Artificial. Na mesoescala, duas abordagens são propostas. A primeira inclui as fibras através de elementos de interface, e a segunda através de um novo elemento compósito fibra-matriz. Os modelos desenvolvidos permitem avaliar a evolução do dano, o processo de propagação de fissuras, e o comportamento global carga-deslocamento do compósito até a ruptura. Resultados experimentais da literatura suportam as conclusões do trabalho. / [en] In recent years, fiber reinforced cement-based materials have gained relevance in the civil engineering industry. Due to its excellent mechanical properties and contribution to crack propagation control, there is a great appeal to its usage as a construction material. However, technical standards for fiber reinforced concrete are still not established. A better understanding of the behavior of cement composite materials requires the representation of the material phases and their interfacial behavior. Stresses and strain distributions, damage evolution and fracture initiation develop at the observation scale of the heterogeneities and help to explain and predict the behavior of concrete at a macroscopic level. The numerical modeling of these composites emerge as challenging and complex problems. For this, it is necessary to define the main mechanisms that describe the material behavior in order to choose the proper mathematical formulation. This dissertation proposes methodologies for the numerical modeling of cement composite materials in a multiscale approach. From the information obtained at the material scale, this work aims at assessing the global behavior of the composite. Numerical and computational procedures will be developed based on the Finite Element Method, Artificial Intelligence techniques and concepts of Computational Damage Mechanics. At the macroscale, an equivalent continuum model is developed through probabilistic and Artificial Intelligence techniques. At the mesoscale, two approaches are proposed. The first includes the fibers through interface elements. The second adopts a new fiber-matrix composite element. With the models developed here, it is possible to evaluate damage evolution, fracture propagation patterns, load-displacement global behavior of the composite upto failure. Experimental results from the literature give support to the conclusions.

[pt] O presente trabalho apresenta a Barragem de Terzaghi em um contexto da evolução tecnológica da construção de barragens de terra, colocando em foco sua construção e desempenho ao longo de 65 anos de operação bem como a verificação atual de suas condições físicas, analisando o impacto que a carga de sedimentos causa sobre a face montante da barragem, em termos de percolação, tensões e deformações. Tendo como base os métodos observacionais aplicados ao longo da operação da barragem e modelagem por elementos finitos, foram traçadas expectativas do fator de segurança para a seção principal da barragem ao longo do tempo. / [en] This work presents the Terzaghi Dam in a technological evolutionary context of the construction of earthfill dams, focusing on its own construction a performance after 65 years in operation as well as a scrutinized assessment of its physical conditions, implying the sediment load on its upstream slope, regarding seepage, stresses and strains. Taking into account observational methods applied throughout the operation of the dam and finite element simulations, the evolution of factor of safety over the dam lifetime and on was evaluated for the main section of the dam.

[pt] O design otimizado de estruturas considerando não-linearidades tem sido amplamente pesquisado nas décadas recentes. A análise de elementos finitos aplicada à otimização topológica é prejudicada pela deformação excessiva de elementos de baixa densidade sob alta compressão, o que impede o processo de encontrar uma solução ótima. Dois métodos, o esquema Interpolação de Energia e a técnica de Hiperelasticidade Aditiva, são implementados para superar essa dificuldade no problema de minimização da flexibilidade, e modelos de materiais hiperelásticos são usados para investigar suas influências na topologia otimizada. O Método das Assíntotas Móveis é usado para atualizar as variáves de projeto cujas sensibilidades foram calculadas pelo método adjunto. A equação de estado é resolvida através do método de Newton-Raphson com um incremento de carga ajustável para reduzir o custo computacional. Resultados de dois problemas de referência são comparado com aqueles já estabelecidos na literatura. O uso de diferentes modelos hiperelásticos apresentou pouca influência no design final da estrutura. O método de Interpolação de Energia foi capaz de convergir para cargas muito maiores que o método padrão, enquanto a Hiperelasticidade Aditiva apresentou dificuldades de convergência em estado plano de deformação. / [en] The optimized design of structures considering nonlinearities has been widely researched in the recent decades. The finite element analysis applied to topology optimization is jeopardized by the excessive deformation of low-density elements under high compression, which hinders the process of finding an optimal solution. Two methods, the Energy Interpolation scheme and the Additive Hyperelasticity technique, are implemented to overcome this difficulty in the minimum compliance problem, and hyperelastic material models are used to investigate their influence on the optimized topology. The Method of Moving Asymptotes is used to update the design variables whose sensitivities were calculated from the adjoint method. The state equation is solved through the Newton-Raphson method with an adjusting load step to reduce computational cost. Results for two benchmark problems are compared with those already established in the literature. The use of different hyperelastic models presented little influence on the final design of the structure. The Energy Interpolation method was able to converge for much higher loads than the default method, while the Additive Hyperelasticity presented convergence difficulties in plane strain.

[pt] Na análise da predição de vida de estruturas a presença de defeitos por trincas é de fundamental importância. Estudos iniciais utilizando a Mecânica da Fratura Linear Elástica mostram que o Fator de Intensidade de Tensão controla o crescimento de trinca por fadiga. Contudo, sobrecargas podem induzir efeitos de memória no material ao reduzir, parar ou acelerar taxas de propagação, comportamento não previsto considerando-se um único parâmetro elástico. Devido às hipóteses consideradas, as soluções fornecidas pelos métodos analíticos não verificam condições importantes de equilíbrio e compatibilidade, especialmente quando são usados campos singulares idealizados de tensão e de deformação na região da ponta da trinca. Neste trabalho uma extensa revisão dos resultados analíticos recentemente publicados para materiais nos regimes elástico e elastoplástico são apresentadas: Williams e HRR com campo de distribuição singular das tensões e Creager-Paris com campo não-singular das tensões com o material elástico. Neste estudo, estas soluções são comparadas com resultados da análise numérica utilizando a discretização por elementos finitos tridimensionais elastoplásticos com o critério de escoamento do material de von Mises na representação do cegamento da ponta da trinca e importantes conclusões são apresentadas com referência aos limites de aplicação dos modelos teóricos bem como a extensão dos requerimentos para o modelo numérico. Além de carregamentos monotonicamente crescentes, o descarregamento foi também considerado na análise numérica, apesar desta condição não estar presente nas referências consideradas para as soluções analíticas. Para as simulações numéricas consideradas algoritmos de elastoplasticidade foram implementados em um framework com arquitetura baseada em plugins. / [en] In structure life prediction analysis the occurrence of crack defects are of paramount importance to be considered. Basic studies using Linear Elastic Fracture Mechanics approach shows that the Stress Intensity Factor (SIF) parameter controls Fatigue crack growth (FCG). However, overloads may induce material memory effects that delay, arrest or accelerate the FCG rate, a behavior that is not described by considering a single elastic parameter. To account for service variable amplitude loadings, the use of a prescribed stress-strain distribution has been proposed in recent research studies, as the driving force of FCG, using the critical damage approach. Due to the assumptions considered in the analytical derivations, important equilibrium and compatibility conditions are violated in the obtained solutions, since an idealized singular stress-strain field at the crack front region is used. In this work a comprehensive review of the recently published analytical results for solutions under the elastic and elastoplastic material behavior regimens are presented: Williams and HRR with singular stress distribution field and Creager-Paris with a non singular stress field, but with the material elastic approach. These solutions are compared, throughout the study, to results obtained from the numerical analysis using a 3D finite element discretization with elastoplastic von Mises yielding criteria employed for the modeling of the crack tip blunt and some comprehensive conclusions are derived regarding application limits of the theoretical models as well as the required extent of the numerical model representation. In addition to monotonic increasing applied loads, unloading conditions were also considered in the numerical analysis, although no reference to this condition is available in the literature, considering the analytical approach. For the numerical simulations presented elastoplastic algorithms were implemented in a plugin based framework.

[pt] Entre materiais renováveis de baixo impacto ambiental o bambu se destaca devido ao seu crescimento rápido, alta resistência à tração, baixo peso específico e flexibilidade. É um material compósito natural com fibras longas alinhadas unidirecionalmente, anisotrópico, que não está sendo usado amplamente nas construções civis por razão de falta de informações técnico-científicas. O bambu é um material inteligente com a distribuição das fibras longitudinais como functionally graded composite material (FGM). O objetivo deste trabalho é continuação dos estudos sobre o comportamento de bambu para aplicação nas obras de engenharia, que inicio na PUC-Rio em 1979. Neste estudo específico foi observado comportamento de bambu da espécie Phyllostachys pubescens (Mosô) (através dos ensaios experimentais) sujeito a carga de torção e flexão pura (realizados com utilização do sistema de reconhecimento de imagem) e simulações computacionais em programa de elementos finitos Abaqus. Os ensaios foram feitos para os corpos de prova inteiros (de 1 metro de comprimento) e os mesmos corpos de prova foram posteriormente fissurados no meio, alargando a fissura para cada ensaio até a ruptura. / [en] Among renewable materials with low environmental impact bamboo stands out because of rapid growth, high tensile strength, low specific weight and flexibility. It is a natural composite anisotropic material with unidirectionally aligned long fibers, which is not widely used in civil construction yet due to lack of technical and scientific information. Bamboo is an intelligent material with the distribution of longitudinal fibers as functionally graded composite material (FGM). The objective of this work is to continue the study the behavior of bamboo for application to engineering works, which started at PUC-Rio in 1979. In this specific study was observed the behavior of Moso bamboo (Phyllostachys pubescens) through the experimental observation of pure torsion and bending (realized with utilization of image recognition system) and computer simulations in finite element program Abaqus. The tests were done for the whole bamboo specimens (1 meter long) and the same specimens were subsequently cracked in the middle extending the crack for each test until failure.

[pt] Este trabalho apresenta o desenvolvimento/montagem de um sistema computacional para análise de fluxo em meios porosos, meios fraturados, porosos fraturados e em combinações destes meios, considerando regime permanente ou transiente, sob condições saturadas e não saturadas. O sistema consiste de quatro programas, três programas de funções específicas interligadas por rotinas de programação feitas na linguagem Cmaismais e o quarto é um visualizador de resultados. O FracGen 3D (Telles, 2006) gera fraturas ou famílias de fraturas de forma determinística ou probabilística. O programa ICEM CFD v.14 divide o domínio de interesse em sub-dominios, através da geração de malha de elementos finitos. O programa FTPF-3D (Telles, 2006) utiliza o método de elementos finitos para discretizar as equações governantes no espaço e em diferenças finitas no tempo, e para resolver a não linearidade, utiliza o método iterativo de Picard ou o método iterativo BFGS e finalmente O Pos3D é o responsável pela visualização dos resultados. Neste trabalho foram desenvolvidos cinco exemplos, dois deles para a validação deste procedimento, e três aplicados a um talude típico do Rio de Janeiro, os quais incluem fraturas verticais e juntas de alívio. Estes casos estudados verificam a influência das fraturas nos meios porosos em termos de carga de pressão, totais e campo de velocidades, para a verificação do comportamento hidráulico dos maciços e de eventuais instabilidades. / [en] This work presents the development/assembly of a computational system for flow analysis in porous media, fractured and fractured porous media and in combination of both media, considering steady or transient states under saturated and unsaturated conditions. The system comprehends four computational programs, three of them of specific functions interconnected by Cplusplus programing routines and the last program is an output viewer. FracGen 3D program (Telles, 2006) generates fractures or fracture families in a determinist or probabilistic way. ICEM CFD v.14 program divides the interest domain in sub-domains by means of the element finite mesh generation. FTPF-3D program (Telles, 2006) uses the element finite method to discretize the governing equations in the space domain and the difference finite method for the time domain and for solving the nonlinearity is used the iterative Picard or BFGS method, so that, finally, Pos3D viewer program is answerable by visualization of the results. In the present dissertation five examples were developed, two of them for the validation of this procedure and the three others applied to a typical slope in Rio de Janeiro, which include vertical fractures and relief joints on their slopes. All those studied cases evaluate the influence of the fractures on porous media in terms of pressure and total heads and velocity fields for verifying of the hydraulic behavior of solid masses and eventual instabilities.

[pt] Emulsões podem ser usadas como agentes de controle de mobilidade em métodos de recuperação avançada de óleo de forma a conseguir um varrido mais eficiente do reservatório. A aplicação de tal técnica requer um entendimento completo de como as emulsões escoam em um material poroso. O comportamento macroscópico do escoamento deve ser determinado baseado no escoamento na escala de poros, em que as emulsões não podem ser tratadas como líquidos não Newtonianos monofásicos, uma vez que os diâmetros das gotas são da mesma ordem de magnitude das gargantas de poros. Experimentos realizados por Cobos et al.(1) sobre o escoamento de emulsões através de micro-capilares com garganta, que serviu como modelo para a geometria de uma garganta de poro conectando dois poros adjacentes, mostrou o efeito da fase dispersa na queda de pressão para diferentes condições de escoamento e características das emulsõs. De forma a ampliar a faixa de número de capilaridade e tamanhos de gota explorados nos experimentos, determinar o efeito da tensão interfacial, razão de viscosidades, geometria do capilar e determinar características do escoamento com mais de uma gota, é estudado o escoamento de gotas suspensas em uma fase contínua escoando através de um micro-capilar. A evolução da interface das gotas é determinada pelo método de level-set, que é acoplado a solução das equações de Navier Stokes, baseada no método de elementos finitos. Foi investigado o efeito da distância inicial entre duas gotas, assim como o efeito da tensão interfacial na velocidade de cada uma das gotas. Os resultados obtidos fornecem uma descrição mais detalhada do escoamento de emulsões em gargantas de poros. / [en] Emulsions can be used as mobility control agent in enhanced oil recovery methods in order to achieve a more efficient sweep of the reservoir. The application of such technique requires full understanding of how emulsions flow in a porous material. The macroscopic flow behavior can be determined based on the pore scale flow, at which emulsions cannot be treated as a single phase non Newtonian liquid, since the drop diameters are in the same order of magnitude of the pore throats. Experiments by Cobos et al.(1) of flow of emulsions through constricted micro capillaries, wich served as a model to the geometry of a pore throat connecting two adjacent pore bodies, have shown the effect of the dispersed phase on the pressure drop for different flow conditions and emulsions characteristics. In order to widen the range of capillary number and drop size explored in the experiments and to determine the effect of interfacial tension, viscosity ratio, capillary geometry and determine the characteristics of the flow with more the one drop, we study the flow of a drops suspended in a continuous phase flowing through a constricted micro capillary. The evolution of the drop interface is determined by the level set method which is incorporated to a fully coupled implicit Navier Stokes solver based finite element method. We investigate the effect of an initial distance between the two drops, as well as the effects of the interfacial tension on the drops speed. The results provide a more detailed description of the flow of emulsions through pore throats.

[pt] O rebaixamento e controle das águas subterrâneas são atividades implementadas e monitoradas de forma contínua ao longo da vida dos projetos de mineração subterrânea ou a céu aberto. A implementação apropriada e eficiente destas atividades depende de estudos hidrogeológicos de grande porte, que permitem avaliar os sistemas de controle mais adequados. A procura da eficiência técnico-econômica destes processos demanda análises numéricas de fluxo tridimensionais de toda a região em estudo, caracterizada por profundas e complexas estratificações de materiais permeáveis abaixo do lençol freático, como normalmente abrangem projetos de mineração, onde as soluções analíticas não podem mais ser aplicadas. O presente trabalho de pesquisa contribui na melhor compreensão das formulações numéricas que representam o comportamento do fluxo subterrâneo, através de dois estudos de caso, o primeiro em uma mina subterrânea e segundo em uma mina superficial. No caso da mina subterrânea foram incorporadas feições cársticas através de elementos discretos 1D dentro de um modelo tridimensional de elementos finitos com o intuito de representar caminhos preferenciais de fluxo. Foram discutidas as vantagens de incorporar tais feições de forma explicita, quantificando o fluxo que passam por estas, que alimentam à mina através de conexões diretas com um rio adjacente. Estes elementos discretos permitem uma representação mais realista do meio hidrogeológico e ao mesmo tempo, uma avaliação mais aprimorada dos efeitos no comportamento do fluxo subterrâneo devido à impermeabilização superficial do rio, como a solução mais coerente para este problema de infiltração. Também foi elaborado um modelo hidrogeológico conceitual para representar o comportamento hidrogeológico de uma mina a céu aberto, desenvolvendo uma sistemática de uso adequado das condições de contorno e de restrição, a calibração deste modelo e a verificação de diferentes cenários de fluxo, como resultado da incorporação das diferentes técnicas de controle das águas avaliadas em regime transiente, que abrange poços de bombeamento, paredes cut-off assim como ponteiras filtrantes e drenos horizontais. As diferentes técnicas modeladas mostraram resultados satisfatórios, sendo que arranjos de várias técnicas, configuradas de forma localizada, resultam ser mais recomendáveis e eficientes para tratar problemas particulares. / [en] Dewatering and groundwater control are activities continuously implemented and monitored throughout the duration of underground mining or open pit projects. The proper and efficient implementation of these activities depends on large hydrogeological studies, for assessing the most appropriate control systems. The demand for technical and economic efficiency of these processes requires three-dimensional flow numerical analysis of the entire study area, characterized by deep and complex stratifications of permeable materials below the water table, as usually cover mining projects, where analytical solutions cannot be applied. The present research contributes to a better understanding of the numerical formulations that represent the behavior of groundwater flow via two case studies—the first in an underground mine and the second in a surface mine. In the case of the underground mine, 1D discrete elements were incorporated within a finite-three dimensional model in order to represent preferential flow paths. The advantages of incorporating such features explicitly to quantify the flow passing through them, that feed the mine through direct connections with an adjacent river, were discussed. The use of discrete elements allows for a more realistic representation of the hydrogeological environment and, at the same time, a more refined assessment of the effects on the behavior of groundwater flow due to surface sealing of the river, as the most coherent solution to this infiltration problem. Furthermore, a conceptual hidrogeologic model representing the hydrogeological behavior of an open pit mine was created, developing a system of the appropriate use of boundary and constraint conditions, calibrating the model and verifying different flow scenarios, as a result of the incorporation of different water control techniques assessed in transient regime, such as covering pumping wells, cut-off walls as well as wellpoints and horizontal drains. The different techniques modeled showed satisfactory results, with various arrangements of techniques configured in a localized form resulting in the most desirable and efficient treatments for particular problems.