Introdução aos caudais Estabilidade da Barragem e Fluxo de Polpa

As barragens de rejeitos de mineração estão entre as maiores estruturas projetadas na Terra, projetadas para armazenar os subprodutos de resíduos de grãos finos (cauda) gerados durante a extração mineral. Esses apreendidos devem permanecer estáveis por décadas, muitas vezes sob condições ambientais severas e cargas sísmicas. A falha de uma barragem de rejeitos pode liberar milhões de metros cúbicos de pastagem tóxica, causando danos ambientais catastróficos e perda de vida. Falhas recentes de alto perfil, como no Brumadinho no Brasil e no Monte Polley no Canadá, têm intensificado a demanda por métodos de análise de estabilidade mais rigorosos.

As avaliações tradicionais de estabilidade geotécnica dependem de métodos de equilíbrio limite e pressupostos simplificados de pressão pore-pressupostos. Embora essas abordagens sejam úteis, muitas vezes não conseguem captar o comportamento complexo e dependente do tempo da lama dentro da barragem. As lamas são fluidos multifásicos, não-newtonianos que evoluem ao longo do tempo à medida que as partículas se estabilizam, consolidam e dessicam. A dinâmica computacional dos fluidos (CFD) oferece uma forma de simular esses processos com alta fidelidade. ANSYS Fluent, uma plataforma líder em CFD, fornece aos engenheiros as ferramentas necessárias para modelar o fluxo de lama, a evolução da pore-pressão e as distribuições de tensão dentro das barragens de caudais. Este artigo fornece um guia abrangente para modelar o fluxo de lama para análise de estabilidade usando o Fluente ANSYS, abrangendo a teoria fundamental, o fluxo de trabalho e considerações práticas.

Compreendendo o fluxo de lama em barragens de caudamentos

As lamas são suspensões de partículas sólidas finas (tipicamente < 200 μm) em água, muitas vezes com uma elevada concentração de sólidos (40%-60% em massa). O comportamento do fluxo destas suspensões é regido pela distribuição do tamanho das partículas, mineralogia e química da água. Em altas concentrações, as lamas apresentam características não-newtonianas: podem ser tixotrópicas (viscosidade dependente do tempo), corta-cisas (viscosidade diminui com a taxa de cisalhamento), ou têm uma tensão de rendimento que deve ser superada antes do início do fluxo. Estas propriedades afectam significativamente a forma como a lama se move através da barragem, como se consolida e como as pressões porosas se desenvolvem.

Modelos reológicos para a caça às vacas

Para modelar o fluxo de lama com precisão no Fluente ANSYS, os engenheiros devem definir um modelo reológico adequado. Os modelos comuns incluem:

  • Plastico de Bingham – uma relação linear entre tensão de cisalhamento e taxa de cisalhamento após uma tensão de escoamento ser ultrapassada.Adequado para pastas de alta concentração que se comportam como um fluido viscoplástico.
  • Herschel-Bulkley – um modelo mais generalizado que explica o corte ou o enfiamento após o rendimento. Amplamente utilizado para rejeitos, porque se encaixa bem em dados experimentais sobre uma ampla gama de taxas de cisalhamento.
  • Lei de poder – útil quando a tensão de rendimento é insignificante, mas muitas vezes inadequada para lamas liquidadas ou altamente concentradas.
  • Casson – historicamente usado para lama e lamas de mineração, embora menos comuns hoje.

A caracterização experimental (por exemplo, utilizando um reômetro rotacional) é essencial para obter os parâmetros para estes modelos. Sem dados reológicos precisos, as simulações CFD podem produzir resultados enganosos.

Ajustamento e consolidação

Com o tempo, as partículas se instalam sob a gravidade, formando um leito denso e compactado no chão da barragem. Esta sedimentação altera a reologia da pasta localmente e afeta a evolução da pressão dos poros. Em muitas barragens de rejeitos, a pasta é depositada em camadas finas (beaching), permitindo que a água se decante e os sólidos se solidifiquem. A ANSYS Fluent pode modelar esta fixação transitória utilizando uma abordagem multifásica euleriana-euleriana com uma fase de sólidos granulares ou uma abordagem euleriana-lagrangiana com partículas discretas. A escolha depende da carga de sólidos e do nível de detalhe necessário.

O papel do CFD na análise de estabilidade da barragem de caudamentos

A análise de estabilidade de uma barragem de rejeitos envolve avaliar o equilíbrio entre as forças motrizes (gravidade, pressão de seepage) e forças resistindo (força de corte dos materiais da barragem e os próprios rejeitos). A pressão de pore-água é um fator crítico porque reduz o estresse efetivo, diminuindo a resistência ao cisalhamento. CFD permite aos engenheiros simular a evolução da pressão de pore em toda a barragem, à medida que a lama é depositada, consolida e drena.

As capacidades multifásicas e porosas dos meios ANSYS Fluent permitem tratar os rejeitos como uma mistura fluida-sólida que se move através de um meio poroso deformável (os rejeitos consolidados e a estrutura da barragem). Acoplar CFD com análise geotécnica de elementos finitos (por exemplo, através de uma troca de dados unidirecional ou bidirecional) proporciona ainda mais insight, mas só o Fluent já pode produzir campos de pressão, perfis de velocidade e locais de superfície livre que se alimentam em cálculos de estabilidade convencionais.

Técnicas de modelagem em Fluente ANSYS

A seleção do modelo físico correto em Fluent é a pedra angular de uma simulação bem sucedida. As subseções seguintes detalham as abordagens mais relevantes para o fluxo de lama de rejeito.

Modelos de fluxo multifásico

As lamas de caudas são inerentemente multifase. ANSYS Fluent oferece várias opções:

  • Volume de Fluido (VOF) modelo – melhor quando se rastreia uma interface distinta entre a pasta e o ar (ou água), por exemplo, durante o encalhe ou no lago decantado. VOF captura dinâmicas de superfície livre, mas não modela fases distintas de partículas; a pasta é tratada como um único fluido com uma reologia definida pelo usuário.
  • Modelo euleriano-euleriano – trata cada fase (água, sólidos) como contínua interpenetrante com sua própria velocidade, temperatura e fração volumétrica. Esta é a abordagem mais rigorosa para as lamas densas onde as interações partículas-partícula são significativas. Requer modelos de fechamento para a pressão de arrasto, elevação e granular (por exemplo, a teoria cinética do fluxo granular). É computacionalmente caro, mas fornece informações detalhadas sobre a distribuição de sólidos.
  • Modelo de mistura – um modelo euleriano simplificado que resolve uma equação de momento único para a mistura e usa relações algébricas para velocidades de deslizamento entre fases. Adequado quando velocidades de fixação são bem caracterizadas e a fração de sólidos é moderada.
  • Modelo de Fase Discreta (DPM) – trata partículas como elementos individuais de Lagrangian em um fluido contínuo. Apropriado para baixa carga de sólidos (<10% volume) mas impraticável para rejeitos densos devido ao custo computacional.

Para barragens típicas de rejeitos de alto-sólido, recomenda-se a abordagem Euleriana-Euleriana ou o modelo de mistura (com fechamentos adequados de arrasto e de fixação). O modelo VOF pode ser utilizado em conjunto com um desses modelos para rastreamento de superfície livre.

Modelação de Reologia Não-Newtoniana

A ANSYS Fluent permite aos utilizadores definir viscosidade não-newtoniana através de Funções Definidas pelo Utilizador (UDFs) ou seleccionando um dos modelos incorporados (Bingham, Herschel-Bulkley, etc.). Para uma simulação de caudamentos:

  1. Escolha o modelo reológico adequado com base em dados experimentais. Guarde os parâmetros (stress de rendimento, índice de consistência, índice de fluxo) na base de dados do material.
  2. Se a fração de sólidos varia espacialmente, use uma UDF para fazer os parâmetros reológicos depender da fração de volume local de sólidos, o que é crítico porque a lama perto do ponto de descarga é menos concentrada do que o leito fixo.
  3. Habilite zonas de malha dinâmica ou de mídia porosas para regiões que se consolidaram em um bolo sólido.

Uma armadilha comum é usar uma aproximação de viscosidade Newtoniana, que falha completamente o comportamento de tensão de rendimento. Uma pequena tensão de rendimento (por exemplo, 5 Pa) pode afetar drasticamente a dissipação de pressão poro- e propagação de fluxo dianteiro.

Mídia porosa e consolidação

Uma vez que os rejeitos se instalam, eles se comportam como um meio poroso com permeabilidade decrescente à medida que a consolidação prossegue. Em Fluent, você pode definir zonas porosas com permeabilidade direcional e uma função definida pelo usuário para atualizar a permeabilidade com base na fração sólida local. Isso permite simulação do processo de consolidação: inicialmente, a fração sólida baixa (alta permeabilidade) gradualmente se transiciona para uma torta de baixa permeabilidade. O modelo de mídia porosa Fluente calcula o fluxo através da torta, fornecendo distribuições de poros de pressão necessárias para análise de estabilidade.

Passos no processo de modelagem

Um fluxo de trabalho robusto garante que a simulação seja precisa e computacionalmente viável.

Geometria e geração de malha

A geometria de uma barragem de rejeitos inclui o aterro da barragem, o piso da bacia, o lago de decantação, e o gasoduto de descarga de chorume ou espigota. Em muitos casos, a barragem cresce com o tempo; um modelo 3D pode incorporar a construção em fase de estágio usando múltiplas regiões ou movimento dinâmico de malha. A malha deve capturar características críticas:

  • Superfície livre (tipicamente refinada perto da interface ar-ar).
  • Zona de descarga de lama (gradientes de alta velocidade).
  • Represa interior e fundação (onde as pressões de poros são calculadas).

As malhas híbridas (tetraédricas em regiões complexas, hexaédricas em secções rectas) são comuns. Um estudo de independência de grelhas deve ser realizado para pelo menos duas malhas (coarse/fina) para garantir que as pressões de poros e as frentes de fluxo calculadas não são sensíveis à malha.

Condições de Limite e Propriedades do Material

As condições de fronteira chave incluem:

  • Inlet – taxa de entrada de lamas, fração de volume sólido, parâmetros reológicos (ou um UDF para variar).
  • Outlet – decantar superfície do lago ou dreno do dedo do pé. Normalmente definido como saída de pressão com uma cabeça hidrostática fixa.
  • Muros – face da barragem e piso da bacia: o deslize é adequado para as fases de fluido, com tensão de cisalhamento calculada a partir da reologia da suspensão.
  • Simetria – utilizada se a geometria da barragem for simétrica.

As propriedades do material devem ser definidas para cada fase. Para a fase dos sólidos num modelo euleriano, especifique densidade, diâmetro de partículas e modelos de viscosidade granular. A fase fluida (água) utiliza propriedades padrão. Lembre-se de especificar um tamanho de passo de tempo apropriado: esquemas explícitos podem exigir passos de tempo muito pequenos (por exemplo, 0,001 s) para manter a estabilidade em regiões de alto rendimento-stress.

Resolver configurações e simulação

Para simulações transitórias de deposição e consolidação de lamas:

  • Use o solucionador baseado em pressão com o algoritmo PISO ou SIMPLE para acoplamento de velocidade de pressão.
  • Active a equação da fração de volume para o modelo Euleriano, ou defina o esquema explícito do VOF se rastrear a superfície livre.
  • Selecione uma discretização espacial de segunda ordem para a fração de volume e momento para reduzir a difusão numérica.
  • Defina cuidadosamente os factores de sub-relaxamento — os fluxos de tensão podem ser sensíveis a factores de relaxamento inferiores a 0,3.
  • Inicialize o domínio com um pequeno volume de ar acima do leito de caudas iniciais. Use regiões de patch para definir a altura inicial da pasta.

Execute a simulação por tempo físico suficiente para capturar o ciclo de deposição completo (horas a dias de tempo real). Verifique a convergência em cada passo temporal: os resíduos devem cair pelo menos três ordens de magnitude.

Pós-Processo e Validação

Após a simulação, as ferramentas de pós-processamento da ANSYS Fluent podem visualizar:

  • Contornos de pressão de poros e estresse eficaz.
  • Fração de volume de sólidos para mostrar o perfil de praia e frente de consolidação.
  • Lotes vetoriais de velocidade de suspensão, indicando zonas de fluxo preferencial que podem levar à tubulação.
  • Localização de superfície livre para comparação com fotos aéreas ou imagens de drones.

A validação é essencial. Compare pressões de poros computadas em várias profundidades com leituras de piezômetros da barragem real. Se a simulação corresponder aos dados de campo dentro de tolerância aceitável (por exemplo, ±15%), o modelo pode ser usado para prever o comportamento futuro ou testar modificações de projeto.

Considerações Práticas para Análise de Estabilidade

O objetivo final da simulação CFD é fornecer dados para uma análise quantitativa de estabilidade. As seguintes saídas são particularmente úteis:

  • Campo de pressão poro – directamente ingerido no software de equilíbrio limite (por exemplo, SLIDE, GeoStudio) através de um mapa de pressão. Altas pressões de poros na massa de rejeitos reduzem a resistência efectiva e aumentam o risco de liquefação estática.
  • Distribuição de tensão de corte – Fluent pode calcular a tensão de cisalhamento exercida pela pasta no chão e paredes da barragem. Comparando-os à resistência de cisalhamento dos materiais de fundação indica onde é provável a falha.
  • Localização de superfície livre – a altura da superfície fréática dentro da barragem é uma entrada clássica para gráficos de estabilidade.Uma superfície fréática excessivamente alta sugere drenagem inadequada.
  • Caminhos de fluxo desleixados – vetores de velocidade podem revelar erosão interna (piping) se zonas de alta velocidade forem encontradas perto do núcleo da barragem ou fundação.

Os engenheiros frequentemente combinam os resultados do CFD com uma análise de tensão de elementos finitos (por exemplo, usando Ansys Mechanical) para calcular o fator de segurança contra falha de deslizamento e rotação. Esta abordagem acoplada fornece uma imagem mais completa do que qualquer dos métodos isoladamente.

Estudos de Caso e Aplicações

Várias operações de mineração têm usado o ANSYS Fluent para analisar o gerenciamento de rejeitos. Por exemplo:

  • Uma mina de cobre no Chile modelou a deposição sequencial de rejeitos em uma barragem de vale-preenchimento para otimizar a localização da espiga e reduzir o risco de sobre-recipiente. A simulação mostrou que mover o ponto de descarga 20 m acima reduziu a altura da superfície freática em 3 m, melhorando o fator de segurança de 1,2 para 1,5.
  • Uma mina de ouro na África Ocidental utilizou um modelo Eulerian-Eulerian com a reologia de Herschel-Bulkley para avaliar a taxa de dissipação de poros em uma barragem de rejeitos de paredes grossas. Os resultados levaram à instalação de camadas de drenagem horizontais adicionais, que encurtaram o tempo de consolidação de 18 meses para 10 meses.

Embora os detalhes sejam proprietários, tais exemplos demonstram o valor prático do CFD na engenharia de barragens de rejeitos.

Limitações e desafios

Nenhuma simulação é perfeita. As principais limitações da abordagem descrita incluem:

  • Custo computacional – simulações Eulerian-Eulerian de alta resolução com modelos acoplados a DEM podem levar dias em uma estação de trabalho multi-core. Escalar para modelos transientes 3D em escala completa pode exigir computação de cluster.
  • Incerteza reológica – as propriedades de rejeitos mudam com o teor mineral, a forma das partículas e a química da água.A simulação é tão boa quanto os dados reológicos de entrada, que podem não ser representativos de todas as fases de deposição.
  • Leis de constituição geotécnica simplificadas – Fluent não modela nativamente o comportamento de tensão da estrutura de barragem sólida. Acoplamento de CFD com FEA geotécnica ainda é uma tarefa especializada que requer mapeamento cuidadoso de dados.
  • Validação de dados – muitas minas carecem de instrumentação abrangente (piezômetros, inclinômetros) para comparar. Sem validação, a simulação continua sendo uma ferramenta preditiva com precisão desconhecida.

Apesar desses desafios, o uso de CFD para a estabilidade de rejeitos tem crescido rapidamente à medida que o poder computacional aumenta e como a pressão regulatória demanda avaliações de risco mais quantitativas.

Conclusão

Modelando o fluxo de lamas em barragens de rejeitos de mineração usando a ANSYS Fluent, é uma base poderosa para a análise de estabilidade. Ao capturar a natureza multifásica da pasta, sua reologia não-newtoniana e a evolução da pressão de poros ao longo do tempo, os engenheiros podem identificar mecanismos de falha que são mais simples. O fluxo de trabalho – criação de geometria, malhagem, condições de contorno, configuração de solucionadores, pós-processamento e validação – exige atenção cuidadosa aos detalhes, mas o pagamento é uma compreensão mais profunda do comportamento da barragem e a capacidade de projetar impundments mais seguros e confiáveis. À medida que a indústria se move para uma gestão de riscos mais rigorosa, ferramentas de CFD como a ANSYS Fluent desempenharão um papel cada vez mais central para garantir que as barragens de rejeito permaneçam estáveis ao longo de toda a sua vida útil e além.

Para mais informações sobre a física subjacente e as melhores práticas, consultar ANSYS Fluent documentation, o Conselho Internacional de Mineração e Metais (ICMM) orientações de gestão de rejeitos, e os Geoenginer.org recursos de rejeitos de barragens[.Insight adicional sobre reologia de chorume pode ser encontrado através de ]Resunção de aplicações de lamas de mineração.