Preconditioning for a Cahn-Hilliard-Navier-Stokes model for morphology formation in organic solar cells

Resumo Técnico: Pré-condicionamento para um Modelo de Cahn-Hilliard-Navier-Stokes para a Formação de Morfologia em Células Solares Orgânicas

1. Definição do Problema

A fabricação de fotovoltaicos orgânicos (OPVs) depende da secagem de uma mistura à base de solvente contendo um polímero (doador de elétrons), um aceitador não-fulereno (NFA) e um solvente. O desempenho do dispositivo resultante é criticamente dependente da morfologia do filme fino, que se forma através de uma interação complexa entre a separação de fases espontânea e a evaporação do solvente. Otimizar experimentalmente as condições de produção para alcançar uma morfologia desejada é demorado e muitas vezes inviável devido à sensibilidade do processo às propriedades dos materiais e aos parâmetros de processamento.

O desafio matemático reside em modelar essa evolução morfológica. O sistema envolve:

• Dinâmica de fluidos macroscópica: O movimento da mistura de fluido.
• Separação de fases microscópica: A separação dos componentes do doador, aceitador e solvente.
• Evaporação: A transição do solvente de líquido para gás e sua difusão no ar circundante.

Esses processos físicos são acoplados, levando a um sistema de equações diferenciais parciais (PDEs) que, quando discretizadas, resultam em sistemas lineares e não lineares de grande escala e mal condicionados. Resolver esses sistemas de forma eficiente, particularmente em relação às mudanças nos parâmetros de discretização (tamanho da malha, passo de tempo), é o principal obstáculo numérico abordado neste trabalho.

2. Metodologia

Modelagem Matemática
Os autores propõem um modelo de contínuo acoplando três componentes principais:

1. Equações de Navier-Stokes: Descrevem a velocidade macroscópica ($v$) e a pressão ($p$) da mistura de fluido.
2. Equações de Cahn-Hilliard: Descrevem a evolução das frações volumétricas ($\phi_i$) para o polímero, NFA, solvente e ar, impulsionadas por potenciais químicos ($\mu_i$) derivados de um funcional de energia livre.
3. Equação de Allen-Cahn: Descreve a evolução de um parâmetro de ordem de vapor ($\Phi_v$) para modelar a mudança de fase do solvente de líquido para gás. Esta equação não conserva a massa, refletindo a perda de solvente para o ambiente.

O modelo incorpora um funcional de energia livre baseado na teoria de Flory-Huggins, interpolando entre as fases líquida e gasosa usando uma função suave do parâmetro de ordem de vapor. As condições de contorno incluem no-slip (sem deslizamento) para a velocidade e condições de fluxo específicas na fronteira superior para modelar a saída do solvente e a entrada de ar.

Discretização

• Espacial: As equações são discretizadas utilizando o Método de Elementos Finitos (FEM). Os autores empregam uma formulação dividida (split formulation) para lidar com os requisitos de continuidade das funções de base. Elementos Taylor-Hood são usados para as equações de Navier-Stokes (velocidade de grau $\ell+1$, pressão de grau $\ell$), enquanto elementos Lagrange contínuos são usados para os campos de fase e potenciais químicos.
• Temporal: Um esquema de passo de tempo semi-implícito é utilizado. Termos lineares são tratados implicitamente, enquanto termos não lineares provenientes do potencial de energia livre são tratados explicitamente.

Estratégia de Pré-condicionamento
A contribuição central do artigo é o desenvolvimento de um solver iterativo pré-condicionado robusto para os sistemas de ponto de sela resultantes.

• Pré-condicionamento em Bloco: Os autores constroem pré-condicionadores de bloco diagonal que exploram a estrutura de ponto de sela das equações discretizadas.
• Aproximação do Complemento de Schur: Para os blocos de Cahn-Hilliard e Navier-Stokes, o pré-condicionador aproxima o complemento de Schur.
  • Para Cahn-Hilliard, o complemento de Schur é aproximado por um produto de matrizes (abordagem de correspondência ou matching approach) para evitar a inversão direta de uma soma de grandes matrizes.
  • Para Navier-Stokes, um pré-condicionador de convecção-difusão-pressão (PCD) é usado, baseado no comutador dos operadores.
• Solvers Lineares: Os blocos diagonais e as aproximações do complemento de Schur são invertidos usando métodos de Multigrid Algébrico (AMG) (especificamente Ruge-Stüben e Smoothed Aggregation). O sistema global é resolvido usando o método do Gradiente Conjugado Generalizado (GMRES).

3. Principais Contribuições

• Nova Formulação de Modelo: O artigo apresenta um modelo acoplado ligando dinâmica de fluidos, separação de fases e evaporação especificamente para o contexto da formação de morfologia de células solares orgânicas, estendendo modelos anteriores de Cahn-Hilliard-Navier-Stokes para incluir efeitos de evaporação via uma equação de Allen-Cahn.
• Pré-condicionamento Robusto: Os autores projetam e implementam uma técnica de pré-condicionamento que garante que o número de iterações do GMRES permaneça pequeno e independente do tamanho da malha ($h$) e do passo de tempo ($\tau$). Esse comportamento "robusto ao parâmetro" é crítico para a escalabilidade do solver em simulações 3D de grande escala.
• Melhorias na Estabilidade Numérica: Para abordar instabilidades numéricas decorrentes do funcional de energia livre logarítmico, os autores propõem uma aproximação polinomial (inspirada em potenciais de poço duplo) que mantém as características físicas desejadas (mínimos de separação de fase) enquanto melhora a estabilidade numérica.

4. Resultados Numéricos

Os autores validam o modelo e o solver através de extensos experimentos numéricos em domínios 1D, 2D e 3D:

• Desempenho do Solver: Em todos os cenários testados (variando tamanhos de malha de $50^2$ a $400^2$ em 2D, e de $10^3$ a $50^3$ em 3D, e variando passos de tempo), as contagens de iteração do GMRES permaneceram baixas e estáveis. O tempo de execução escalou linearmente com os graus de liberdade, confirmando a eficiência do pré-condicionador.
• Comportamento Físico:
  • Evaporação: Simulações de evaporação de solvente puro mostraram a depleção esperada do solvente e a formação de uma camada de vapor.
  • Separação de Fases: No modelo completo (polímero, NFA, solvente), as simulações capturaram com sucesso a separação de fases das camadas ativas. O uso da aproximação de energia livre polinomial permitiu simulações estáveis de morfologias complexas em 3D, que anteriormente eram propensas a instabilidades numéricas com o potencial logarítmico.
• Morfologia 3D: O artigo demonstra com sucesso a formação de morfologias 3D, mostrando a evolução dos domínios de polímero e NFA durante o processo de secagem.

5. Significância e Alegações

O artigo alega que a modelagem precisa da formação de morfologia é uma ferramenta indispensável para otimizar a produção de OPV, já que o método de tentativa e erro experimental é impraticável. A significância deste trabalho reside em:

1. Viabilização de Simulações 3D: Ao fornecer um esquema de pré-condicionamento robusto, os autores permitem a simulação do processo de formação de morfologia em três dimensões, o que é um passo necessário para a modelagem realista de dispositivos.
2. Robustez de Parâmetros: O solver proposto é robusto em relação aos parâmetros de discretização, tornando-o adequado para refinamento de malha adaptativo e estratégias de passo de tempo variáveis, exigidos para simulações físicas complexas.
3. Captura Qualitativa: O modelo captura qualitativamente o comportamento da morfologia do filme durante a secagem, incluindo a interação entre a separação de fases e a evaporação.

Os autores concluem que, embora o modelo capture a física essencial, a capacidade de realizar essas simulações em 3D fornece uma ferramenta poderosa para investigar condições de produção otimizadas, potencialmente levando a um melhor desempenho de OPV. Eles não alegam ter resolvido o problema inverso de encontrar as condições ideais, mas sim fornecem o arcabouço numérico necessário para estudar o problema direto da evolução da morfologia.