Shear Banding in Simulations of Polymer Melts

Este estudo compara simulações numéricas de polímeros sob fluxo de cisalhamento com um modelo teórico acoplado à dinâmica de tubos e à liberação de restrições convectadas (CCR), revelando que o cisalhamento bandado ocorre quando o número de emaranhamentos excede um valor crítico dependente da rigidez à flexão, embora existam desvios devido às limitações dos modelos de tubo em líquidos parcialmente desenredados sob fluxo intenso.

O essencial

Imagine que você está tentando misturar um pote enorme de mel muito grosso e cheio de fios de lã emaranhados. Se você mexer devagar, tudo se move junto. Mas, se você começar a mexer muito rápido, algo estranho acontece: em vez de todo o mel girar uniformemente, ele se divide em camadas. Algumas camadas giram como loucas, enquanto outras ficam quase paradas. Na física, chamamos isso de "bandamento de cisalhamento" (shear banding).

Este artigo é uma investigação científica sobre por que isso acontece em plásticos derretidos (polímeros) e se conseguimos prever quando vai ocorrer. Os autores usaram supercomputadores para simular esses "fios de lã" e compararam os resultados com uma teoria matemática.

Aqui está a explicação simplificada, usando analogias do dia a dia:

1. O Cenário: A "Festa" dos Fios Embaralhados

Pense em um polímero (como plástico derretido) como uma panela cheia de espaguete cozido.

• Em repouso: Os fios estão todos emaranhados uns nos outros. É difícil se mover porque você está preso aos vizinhos.
• Ao mexer (cisalhamento): Quando você aplica força para misturar, os fios tentam se alinhar e esticar.

2. O Problema: O "Efeito Rebote" (CCR)

Aqui entra o conceito chave do artigo: Liberação de Restrição Convectiva (CCR).
Imagine que, enquanto você puxa os fios de espaguete para esticá-los, eles também começam a se soltar dos nós que os prendiam.

• A Teoria: Se os fios se soltarem muito rápido, o material fica "mais mole" e flui de forma estranha, criando essas faixas (bandas) de velocidade diferente.
• O Parâmetro Beta ($\beta$): Os cientistas criaram um "botão de controle" chamado $\beta$ na matemática. Esse botão mede quão rápido os fios se soltam quando são puxados.
  • Se o plástico é muito rígido (como um espaguete duro), ele se solta de um jeito.
  • Se é flexível (como um espaguete mole), se solta de outro.

3. A Grande Descoberta: Quantos Fios são Necessários?

A teoria matemática (o modelo DO) previa uma regra simples:

"Para que o efeito de divisão em faixas (bandamento) aconteça, você precisa de muitos fios emaranhados (chamados de $Z$)."

É como se dissessem: "Se você tiver menos de 10 fios, o espaguete vai girar uniformemente. Se tiver mais de 10, ele vai começar a se dividir em faixas rápidas e lentas."

O que os computadores mostraram?
Os autores rodaram simulações com milhões de "fios" virtuais.

• A boa notícia: A teoria acertou em cheio! Eles viram que, de fato, só acontece o bandamento quando há muitos fios emaranhados.
• A nuance: A teoria às vezes exagera um pouco na quantidade de força necessária para ver isso acontecer, mas a direção geral está correta.

4. O Mistério dos Plásticos Reais

Aqui está a parte mais interessante e um pouco frustrante:

• Nas simulações: O bandamento acontece claramente.
• Na vida real (experimentos com plásticos): É muito difícil ver isso.

Por que?
Os autores sugerem que os plásticos reais têm uma estrutura molecular mais complexa (como pequenos "joelhos" ou dobras nas cadeias) que faz com que eles se soltem (o parâmetro $\beta$) de uma forma diferente das simulações simples.

• Analogia: Imagine que nas simulações, os fios são feitos de um material que se solta facilmente quando puxado. Nos plásticos reais, eles podem ter "travas" ou "ganchos" internos que impedem que eles se soltem tão rápido, evitando que o material se divida em faixas.

5. Conclusão: O Que Aprendemos?

1. A Teoria Funciona: O modelo matemático é uma ferramenta poderosa para prever quando um plástico vai começar a se comportar de forma estranha (dividir em faixas) sob pressão.
2. O Limite: Existe um número mágico de emaranhamentos. Abaixo dele, o plástico flui bem. Acima dele, o caos (bandamento) pode começar.
3. O Futuro: Isso ajuda os engenheiros a escolherem quais plásticos usar. Se você quer um plástico que não quebre ou flua de forma irregular em máquinas de alta velocidade, precisa evitar aqueles que têm muitos emaranhamentos e que se soltam muito rápido.

Resumo em uma frase:
O artigo mostra que, ao mexer plásticos derretidos muito rápido, eles podem se dividir em camadas de velocidade diferente, e isso depende de quantos "nós" existem entre as cadeias de plástico e de quão rápido esses nós se desfazem quando puxados. A matemática prevê isso bem, mas a realidade dos plásticos é um pouco mais teimosa do que os computadores imaginam.

Resumo técnico

1. Problema e Contexto

O fenômeno de bandamento de cisalhamento (shear banding) em polímeros entrelaçados ocorre quando o fluxo de um fluido complexo se torna heterogêneo, dividindo-se em bandas com taxas de cisalhamento diferentes sob uma tensão de cisalhamento constante.

• O Dilema Experimental vs. Teórico: Embora teorias de modelos de tubo (como o modelo Doi-Edwards) prevejam curvas constitutivas não monótonas (com um platô de tensão) que levam ao bandamento, a maioria dos experimentos em fusos poliméricos (polymer melts) mostra apenas um leve aumento na tensão, sem platôs claros. Isso é frequentemente atribuído a efeitos de deslizamento na parede (wall slip) e fratura de borda (edge fracture) que mascaram o bandamento real.
• O Papel das Simulações: Simulações de dinâmica molecular (DM) oferecem um ambiente controlado sem essas instabilidades experimentais, permitindo observar o bandamento em estado estacionário. No entanto, é necessário validar se os modelos teóricos atuais conseguem prever corretamente quando e sob quais condições esse bandamento ocorre em fusos.
• Objetivo: O artigo compara as previsões do modelo DO (Dolata-Olmsted), que acopla a dinâmica do tubo (tipo Rolie-Poly) com a liberação de restrições convectada (CCR), com resultados de simulações de dinâmica molecular de polímeros do tipo Kremer-Grest (KG). O foco é determinar como parâmetros físicos (como rigidez da cadeia e número de entrelaçamentos) influenciam a transição para o bandamento.

2. Metodologia

Os autores utilizaram uma abordagem combinada de modelagem teórica e simulação computacional:

• Modelo Teórico (Modelo DO):

  • Baseia-se em equações de movimento que acoplam o tensor de conformação do tubo ($A$) e o número reduzido de "dobras" (kinks) da primitiva ($\nu$), representando a perda de entrelaçamentos sob fluxo.
  • O modelo inclui o mecanismo de Liberação de Restrições Convectada (CCR), controlado pelo parâmetro $\beta$. O CCR acelera o relaxamento das cadeias esticadas pelo fluxo, reduzindo o número de entrelaçamentos.
  • O modelo possui parâmetros adimensionais: número de entrelaçamentos ($Z$), parâmetro de anisotropia de relaxamento ($\alpha$), parâmetro CCR ($\beta$), estiramento máximo ($\lambda_{max}$) e razão de viscosidade ($\epsilon$).
  • O objetivo é mapear o diagrama de fases onde a curva constitutiva (tensão vs. taxa de cisalhamento) se torna não monótona, indicando a possibilidade de bandamento.

• Simulações de Dinâmica Molecular (DM):

  • Utilizaram o software LAMMPS para simular polímeros semi-flexíveis do tipo Kremer-Grest (esferas conectadas por molas FENE com potencial Lennard-Jones).
  • Variação de Parâmetros: A rigidez da cadeia foi controlada variando o potencial de ângulo de torção ($k_\theta$), criando cadeias flexíveis ($k_\theta=0$) e semi-flexíveis ($k_\theta=1.5, 2.0$).
  • Tamanho da Cadeia: Foram estudadas cadeias com comprimentos ($N_{mon}$) variando de 200 a 2000 monômeros, resultando em diferentes números de entrelaçamentos ($Z$).
  • Algoritmos de Cisalhamento:
    1. SLLOD: Impõe um perfil de velocidade linear uniforme, suprimindo o bandamento (usado para medir propriedades de equilíbrio e determinar $\beta$).
    2. Deformação da Caixa (fix deform): Permite que o perfil de velocidade se desenvolva livremente, permitindo a formação de bandas (usado para observar o bandamento).
  • Análise de Entrelaçamentos: Utilizaram o código Z1+ para analisar as trajetórias das cadeias e calcular o número de entrelaçamentos ($Z_k$) em função da taxa de cisalhamento, permitindo ajustar o parâmetro $\beta$ do modelo teórico.

3. Contribuições Principais

1. Determinação de $\beta$ a partir de Simulações: Os autores desenvolveram um método robusto para extrair o parâmetro de CCR ($\beta$) diretamente das simulações de DM, correlacionando-o com a rigidez da cadeia ($k_\theta$) e a razão entre o comprimento de Kuhn e o comprimento de empacotamento ($b_K/p$).
2. Validação do Diagrama de Fases: Compararam as regiões de bandamento observadas nas simulações com as previsões teóricas do modelo DO, gerando um diagrama de fases ($Z$ vs. $\beta$) que valida semi-quantitativamente a teoria.
3. Análise da Dependência Histórica: Investigaram como a história de cisalhamento (início a partir do repouso vs. variação gradual da taxa) afeta a formação de bandas, confirmando a natureza metastável do fenômeno.
4. Discussão sobre a Falta de Bandamento em Experimentos: Oferecem uma explicação plausível para a ausência de bandamento em fusos poliméricos experimentais, sugerindo que polímeros reais com rotações de diedro (como polietileno) podem ter valores de $\beta$ muito altos, elevando o número crítico de entrelaçamentos ($Z_c$) necessário para o bandamento a valores inatingíveis em experimentos comuns.

4. Resultados Chave

• Correlação entre Rigidez e $\beta$: O parâmetro $\beta$ aumenta com a rigidez da cadeia ($k_\theta$). Para cadeias flexíveis (KG padrão), $\beta \approx 0.03$, enquanto para cadeias semi-flexíveis, $\beta \approx 0.17 - 0.19$.
• Condição para Bandamento: O modelo prevê que o bandamento ocorre quando o número de entrelaçamentos $Z$ excede um valor crítico $Z_c(\beta, \alpha)$.
  • Simulações confirmam que sistemas com $Z < Z_c$ exibem fluxo uniforme, enquanto $Z > Z_c$ exibem bandamento.
  • Há um acordo semi-quantitativo notável entre a teoria e as simulações, apesar do modelo DO ser uma versão simplificada de modelos microscópicos mais complexos (como o GLaMM).
• Desvios e Limitações:
  • O modelo superestima a faixa de taxas de cisalhamento onde o bandamento ocorre.
  • Existe uma tendência inversa nas simulações: o valor crítico $Z_c$ diminui ligeiramente com o aumento de $\beta$, enquanto o modelo prevê que $Z_c$ deve aumentar. Os autores atribuem isso a aproximações do modelo em fluxos fortes e à possível inadequação da descrição de "tubo" para líquidos parcialmente desenredados.
• Influência da História: O estado final (bandado ou não) depende da história de deformação. Iniciar a partir de uma taxa de cisalhamento alta (não bandada) e reduzir pode impedir a formação de bandas, indicando que o estado não bandado pode ser metaestável.
• Implicação para Polímeros Reais: A análise sugere que polímeros com detalhes atômicos (como polietileno, UA-PE) podem ter valores de $\beta$ muito maiores ($\sim 1.25$) devido a mecanismos de relaxamento específicos. Se $Z_c$ escala com $\beta$, isso explicaria por que fusos poliméricos reais (com $Z$ tipicamente entre 10 e 50) não mostram bandamento, enquanto simulações de modelos simplificados (com $\beta$ baixo) o mostram facilmente.

5. Significado e Conclusão

O trabalho fornece uma ponte crucial entre a teoria de polímeros entrelaçados e a realidade das simulações de dinâmica molecular.

• Validação Teórica: Confirma que o mecanismo de CCR é fundamental para descrever a desentrelaçamento induzido por fluxo e a estabilidade do fluxo homogêneo.
• Explicação para a Discrepância Experimental: Sugere que a falta de evidências experimentais de bandamento em fusos não é necessariamente uma falha da teoria de bandamento, mas sim uma consequência de que polímeros reais possuem parâmetros de relaxamento ($\beta$) que elevam o limiar de entrelaçamento necessário para a instabilidade acima do que é acessível experimentalmente.
• Direção Futura: O estudo destaca a necessidade de modelos mais precisos para descrever a dinâmica em altas taxas de cisalhamento e a importância de considerar a química específica (rigidez, rotações de diedro) ao prever o comportamento reológico de novos materiais poliméricos.

Em resumo, o artigo demonstra que o bandamento de cisalhamento é um fenômeno real em simulações de fusos, controlado pela competição entre o número de entrelaçamentos e a eficiência da liberação de restrições convectada, e oferece uma hipótese sólida para a discrepância entre simulações e experimentos em polímeros reais.