Bulldozing an immersed granular material in a confined channel

Resumo Técnico: Escavação de Material Granular Imerso em um Canal Confinado

Definição do Problema
Este estudo investiga a mecânica de um material granular imerso sendo "escavado" (bulldozing) por um pistão dentro de um canal horizontal confinado. O sistema consiste em uma camada granular densa, saturada de fluido, inicialmente repousando sobre a placa inferior de um canal com uma altura de fenda $b$. Um pistão rígido move-se com uma velocidade constante, deslocando tanto os grãos quanto o fluido intersticial. O problema é caracterizado por interações complexas entre o atrito grão-grão, o arrasto grão-fluido e o atrito grão-parede. O foco específico é a formação de uma região de "ponte" (bridged region) onde os grãos se acumulam para preencher a fenda do canal, uma zona de transição onde a altura da camada granular aumenta do seu valor inicial para a altura da fenda, e o campo distante não perturbado. Este trabalho generaliza estudos anteriores sobre escavação de granulares secos (que ignoravam o fluido) e deslocamento granular impulsionado por capilaridade (que frequentemente dependiam de simulações baseadas em regras) ao desenvolver um modelo de contínuo para o sistema acoplado fluido-grão.

Metodologia
Os autores desenvolvem um modelo de contínuo de ordem reduzida baseado em uma abordagem de duas fases e duas camadas com média de profundidade. O sistema é tratado como um conjunto acoplado de filmes finos: uma camada superior de fluido Newtoniano incompressível e uma camada inferior de material granular saturado de fluido.

• Equações Governantes: Conservação de massa e momento são aplicadas sob uma aproximação de lubrificação (assumindo que o canal é fino e o movimento do pistão é lento, de modo que a inércia macroscópica é desprezível). A fase granular é modelada como um material visco-plástico, denso e poroso.
• Reologia: A mecânica interna da fase granular é descrita usando reologias do tipo $\mu(I)$. Os autores consideram dois regimes distinguidos pelo número de Stokes ($St$):
  • Limite Seco ($St \gg 1$): Efeitos inerciais dominam; a viscosidade granular é descrita pelo número de inércia $I$.
  • Limite Úmido ($St \ll 1$): Efeitos viscosos dominam; a viscosidade granular é descrita pelo número viscoso $I_v$.
• Condições de Contorno: O modelo incorpora atrito do tipo Coulomb nas paredes superior e inferior. Os grãos deslizam apenas quando a tensão de cisalhamento local excede um limiar definido pela tensão normal local e um coeficiente de atrito de parede ($\mu_w$). Um desafio fundamental abordado é a transição entre estados não pontuados (dominados por cisalhamento de fluido) e pontuados (dominados por atrito de parede).
• Regularização: Para lidar com a natureza descontínua de fluidos de tensão de escoamento (yield-stress) e deslizamento de parede (onde regiões de deslizamento e aderência coexistem), os autores introduzem um novo esquema de regularização usando funções tangente hiperbólica. Isso permite uma solução numérica única e contínua sem a necessidade de rastrear explicitamente fronteiras livres entre regiões escoadas e não escoadas.
• Validação Numérica e Experimental: O modelo é resolvido numericamente em um referencial que se move com o pistão. Os resultados são comparados contra uma série de experimentos de laboratório conduzidos em um tubo de vidro quadrado usando esferas de vidro e fluidos intersticiais (água e soluções de glicerol) em várias velocidades de pistão.

Principais Contribuições e Resultados

1. Regime Seco ($St \gg 1$):

   • Na região de ponte, o perfil de velocidade granular é não-monotônico, apresentando uma região central de "plugue" (onde os grãos movem-se como um corpo rígido) entre sanduíche de regiões escoadas próximas às paredes.
   • A tensão normal no pistão cresce exponencialmente com o comprimento da região de ponte, consistente com o efeito Janssen clássico e simulações anteriores de método de elementos discretos (DEM).
   • A forma da zona de transição (a cunha de grãos ascendentes) depende do número de Froude ($Fr$). À medida que $Fr$ aumenta, a zona de transição encurta e a inclinação da pilha granular torna-se mais íngreme, aproximando-se do ângulo de atrito da parede.

2. Regime Úmido ($St \ll 1$):

   • A presença do fluido introduz um comportamento não-monotônico no comprimento da zona de transição ($\lambda$) como uma função do parâmetro de Shields modificado ($\theta$), que mede o cisalhamento do fluido em relação ao peso dos grãos.
   • Em $\theta$ baixo, o fluido tem impacto desprezível, e o sistema comporta-se de forma semelhante ao caso seco.
   • À medida que $\theta$ aumenta, a resistência viscosa da camada de fluido sobrejacente alonga a zona de transição.
   • Além de um $\theta$ crítico, o modelo prevê que uma solução de onda viajante estacionária para a zona de transição deixa de existir; a zona se alonga continuamente e, para $\theta$ suficientemente alto, os grãos falham em formar a ponte na fenda completamente.

3. Comparação Experimental:

   • O modelo apresenta concordância qualitativa ampla com os experimentos em relação à forma da pilha granular e à variação não-monotônica do comprimento da zona de transição com $\theta$.
   • Discrepâncias quantitativas no comprimento da zona de transição foram observadas. Os autores atribuem isso à omissão do atrito das paredes laterais no modelo base. Quando o atrito das paredes laterais é incluído no modelo, a concordância com os dados experimentais melhora significativamente.
   • Perfis de velocidade medidos via Velocimetria de Imagem de Partículas (PIV) geralmente coincidem com a previsão do modelo de aumento monotônico da velocidade com a altura na zona de transição, embora o modelo preveja uma curvatura convexa ($\partial^2_z u_s < 0$) não totalmente capturada nos experimentos.

Significância e Alegações
O artigo afirma fornecer o primeiro modelo de contínuo quasi-1D/2D derivado de mecânica fundamental para a escavação de materiais granulares imersos em geometrias confinadas. Ele consegue unir com sucesso as teorias de fluxo granular seco e problemas de interação fluido-grão.

Os autores afirmam modestamente que seu modelo captura a física essencial das interações entre atrito de parede, tensões internas do tipo viscoso e fluxo de fluido. Eles destacam que a técnica de regularização permite a simulação eficiente de padrões complexos de fluxo friccional sem a necessidade de rastreamento de fronteira livre. Embora o modelo reproduza as principais tendações experimentais, os autores reconhecem limitações, incluindo a suposição de fração sólida constante (ignorando compactação/dilatação), a negligência da inércia macroscópica e o tratamento simplificado dos efeitos das paredes laterais. Eles sugerem que trabalhos futuros podem estender este arcabouço para incorporar fluidos não molhantes para melhor modelar o deslocamento granular impulsionado por capilaridade em células de Hele-Shaw e tubos capilares.