Error analysis of a fully discrete structure-preserving finite element scheme for a diffuse-interface model of tumour growth

Resumo Técnico: Análise de Erro de um Esquema de Elementos Finitos Totalmente Discreto e Preservador de Estrutura para um Modelo de Interface Difusa de Crescimento Tumoral

Enunciado do Problema
Este trabalho aborda a aproximação numérica de um modelo de quatro espécies de crescimento tumoral, termodinamicamente consistente, baseado na teoria contínua de misturas. O modelo acopla uma equação do tipo Cahn–Hilliard para a fração volumétrica do tumor ($u$) com uma equação de reação-difusão não linear para a concentração de nutrientes ($n$). O sistema é governado por um funcional de energia livre de Helmholtz e admite uma lei de dissipação de energia no nível contínuo. Embora estudos anteriores tenham proposto vários esquemas numéricos (incluindo divisão convexa não linear e abordagens SAV tempo-semidiscretas), carecia na literatura um esquema totalmente discreto, linear e preservador de estrutura, acompanhado de uma análise rigorosa de erro para este sistema acoplado específico em três dimensões.

Metodologia
Os autores propõem um método de elementos finitos totalmente discreto e linear que preserva as propriedades estruturais do modelo contínuo. A metodologia combina:

1. Formulação de Variável Auxiliar Escalar (SAV): Para lidar com a densidade de energia livre não quadrática $g(u)$, o sistema é reformulado utilizando uma variável auxiliar $r(t) = \sqrt{E_1[u] + B}$, onde $E_1$ é a parte não quadrática da energia livre. Isso permite a construção de um esquema linear mantendo a estabilidade de energia.
2. Discretização Espacial:
   • Um método de elementos finitos mistos é empregado para a parte de Cahn–Hilliard para reduzir a equação de quarta ordem em um sistema acoplado de equações de segunda ordem.
   • Elementos finitos conformes padrão são utilizados para a equação de reação-difusão.
3. Discretização Temporal: Um esquema de passo de tempo semi-implícito de primeira ordem é utilizado.

O sistema algébrico resultante é linear em cada passo de tempo. A implementação do termo SAV é tratada eficientamente via uma correção de posto um à matriz de bloco esparsa de elementos finitos, passível de resolução através da fórmula de Sherman–Morrison, evitando a necessidade de montagem de matrizes densas ou solvers não lineares.

Contribuições Principais e Análise
A principal contribuição deste artigo é a análise de erro rigorosa do esquema proposto sob suposições de regularidade adequadas da solução exata. A análise estabelece:

• Propriedades Estruturais: O esquema é comprovadamente estável em energia de forma incondicional (satisfazendo uma lei de dissipação de energia discreta para um funcional de energia modificado) e preservador de massa.
• Bem-posedade: A existência e unicidade da solução para o sistema algébrico linear em cada passo de tempo são estabelecidas.
• Estimativas de Erro: Os autores derivam estimativas de erro rigorosas nas normas $L^2$, $H^1$ e $L^\infty$. A análise demonstra:
  • Precisão de primeira ordem no tempo ($O(\tau)$).
  • Precisão ótima no espaço ($O(h^{q+1})$ para $L^2$ e $H^1$, onde $q$ é o grau do polinômio).
• Inovações Técnicas na Análise:
  • Ao contrário de análises anteriores tempo-semidiscretas (ex: em 2D), este trabalho aborda o cenário totalmente discreto em um domínio convexo 3D.
  • Os autores superam a falta de estimativas de regularidade $H^2$ para soluções totalmente discretas (que são apenas $H^1$) ao derivar estimativas delicadas para termos não lineares usando imersões de Sobolev discretas e desigualdades de Gagliardo–Nirenberg discretas.
  • Um passo fundamental envolve provar a limitação das soluções numéricas em $L^\infty$ e estabelecer estimativas de superconvergência para a diferença entre a solução de elementos finitos e sua projeção de Ritz.

Resultados
Os resultados teóricos são validados através de simulações numéricas utilizando o pacote FEniCS. Os experimentos incluem:

1. Crescimento de Tumor Único: Demonstrando crescimento constante e conservação de massa.
2. Agregação de Múltiplos Tumores: Simulando a fusão de dois tumores.
3. Taxas de Convergência: Verificando a convergência temporal de primeira ordem teórica e as taxas de convergência espacial ótimas, mesmo para dados iniciais com regularidade inferior à assumida na teoria.
4. Crescimento Quimiotático 3D: Simulando um tumor tridimensional que exibe comportamento quimiotático em direção a fontes de nutrientes. Notavelmente, a simulação permanece estável mesmo quando a condição de coercividade suficiente ($\lambda > \chi_0^2 \delta$) utilizada na análise teórica não é estritamente satisfeita, sugerindo robustez além dos limites provados.

Significância
O artigo afirma preencher uma lacuna na literatura ao fornecer o primeiro esquema numérico totalmente discreto, linear e preservador de estrutura para este modelo específico de crescimento tumoral, completo com uma análise de erro rigorosa em 3D. O método é significativo porque evita o custo computacional de solvers não lineares enquanto garante a conservação de massa e a estabilidade de energia. A análise estende a compreensão de esquemas preservadores de estrutura para cenções totalmente discretos em dimensões superiores, onde o controle de termos não lineares é substancialmente mais delicado do que em casos tempo-semidiscretos ou de menor dimensão. Os experimentos numéricos confirmam que o método é robusto e capaz de capturar fenômenos característicos como agregação e quimiotaxia.