Exact solution to a class of problems for the Burgers' equation on bounded intervals

Resumo Técnico: Solução Exata para uma Classe de Problemas da Equação de Burgers em Intervalos Limitados

Enunciado do Problema
O artigo aborda a equação de Burgers viscosa, uma equação diferencial parcial não linear que modela o equilíbrio entre convecção e difusão, definida em intervalos limitados genéricos $[l_1, l_2]$. O problema específico envolve determinar a função $w(x, t)$ sujeita a:

1. Equação Governante: $w_t - a^2 w_{xx} + w w_x = 0$, onde $a^2$ representa o coeficiente de viscosidade cinemática ($1/Re$).
2. Condição Inicial: Um perfil inicial não homogêneo $w(x, 0) = w_0(x)$.
3. Condições de Contorno: Condições de contorno de Dirichlet fixas $w(l_1, t) = \alpha_1$ e $w(l_2, t) = \alpha_2$.

Os autores observam que, embora métodos numéricos (diferenças finitas, elementos finitos, métodos espectrais) e técnicas semianalíticas (decomposição de Adomian, perturbação de homotopia) sejam prevalentes, eles frequentemente carecem de soluções de forma fechada ou exigem recursos computacionais significativos. O objetivo é derivar uma solução operacional exata que ofereça uma expressão de forma fechada no domínio de Laplace e uma solução precisa no domínio do tempo por meio de transformadas inversas.

Metodologia
A estratégia de solução proposta baseia-se em um processo de transformação em duas etapas:

1. Transformação de Hopf-Cole: A equação não linear de Burgers é linearizada usando a clássica transformação de Hopf-Cole:
   $$w(x, t) = -2a^2 \frac{u_x(x, t)}{u(x, t)}$$
   Isso converte o problema não linear em um problema de valor de contorno linear para a equação de difusão de calor unidimensional (reação-difusão com termo de reação nulo):
   $$u_t - a^2 u_{xx} = 0$$
   A condição inicial para $u$ é derivada da integral da condição inicial $w_0$, e as condições de Dirichlet para $w$ transformam-se em condições de contorno de Robin homogêneas para $u$.

2. Transformada de Laplace e Solução Operacional: A equação de calor linear é resolvida no domínio de Laplace (domínio $p$). Os autores utilizam uma solução operacional exata previamente estabelecida para equações de reação-difusão lineares com coeficientes constantes em intervalos finitos.

   • O problema é formulado como um problema de valor de contorno de dois pontos no domínio de Laplace.
   • Uma expressão exata para a transformada de Laplace $U(x, p)$ é derivada, envolvendo termos exponenciais dependentes de $\sqrt{p}$ e um sistema linear de quatro equações para determinar os valores de contorno ($U(l_1, p), U_x(l_1, p), U(l_2, p), U_x(l_2, p)$).
   • O termo fonte na equação linearizada é zero, simplificando a solução operacional para depender da transformada da condição inicial $R(x, p)$.

3. Recuperação no Domínio do Tempo: A solução exata para a variável original $w(x, t)$ é obtida implicitamente por meio da transformada inversa de Laplace:
   $$w(x, t) = -2a^2 \frac{\mathcal{L}^{-1}\{U_x(x, p)\}}{\mathcal{L}^{-1}\{U(x, p)\}}$$
   O artigo observa que, embora inversos analíticos possam existir de forma fechada (potencialmente por meio de transformadas de Mellin), algoritmos numéricos altamente eficientes (especificamente o algoritmo de de Hoog et al., implementado como invlap.m no Matlab) são empregados para calcular a solução no domínio do tempo, independentemente da complexidade das expressões de Laplace.

Contribuições e Resultados Principais
O artigo apresenta dois exemplos ilustrativos para validar a metodologia:

• Exemplo 1: Um problema com condição inicial $w_0(x) = \frac{2a^2\pi \sin(\pi x)}{\sigma + \cos(\pi x)}$ e condições de contorno nulas. Os autores derivaram a solução no domínio de Laplace e aplicaram a transformada inversa. A solução resultante no domínio do tempo coincidiu perfeitamente com a conhecida solução exata clássica de forma fechada (Wood, 2014).
• Exemplo 2: Um problema com condição inicial $w_0(x) = \sin(\pi x)$. Sabe-se que a solução exata clássica para este caso é uma série infinita. O método dos autores produziu uma solução exata implícita no domínio de Laplace. Quando invertida numericamente, os resultados alinharam-se estreitamente com a solução da série infinita truncada ("Solução Teórica").

Desempenho Comparativo
O estudo destaca duas vantagens primárias do método proposto sobre abordagens existentes:

1. Eficiência Computacional: No Exemplo 2, o procedimento numérico de inversão de Laplace mostrou-se mais de 100 vezes mais rápido em tempo de execução em comparação com a avaliação da solução da série infinita truncada.
2. Precisão: As soluções de inversão numérica de Laplace apresentaram erro negligenciável em comparação com as soluções de série teóricas. Em contraste, as soluções por diferenças finitas (computadas via pdepe() do Matlab) mostraram discrepância crescente à medida que o coeficiente de viscosidade diminuía, particularmente nos primeiros passos de tempo.

Significado e Alegações
Os autores afirmam que seu trabalho fornece uma solução operacional de forma fechada no domínio de Laplace e uma solução exata implícita no domínio do tempo. Eles afirmam que essa abordagem oferece uma "alternativa nova" a soluções baseadas em séries infinitas ou discretização puramente numérica.

O significado do trabalho é enquadrado da seguinte forma:

• Ponto de Referência: A solução exata derivada serve como um benchmark robusto para validar e testar outros métodos numéricos e semianalíticos.
• Eficiência: O método é computacionalmente superior a expansões em série, que podem ser intensivas em recursos.
• Generalidade: A metodologia lida com soluções suaves, bem como soluções generalizadas com descontinuidades nas derivadas.

O artigo conclui modestamente, afirmando que a metodologia de solução será uma referência para trabalhos futuros, incluindo potenciais extensões para equações de Burgers forçadas em dimensões superiores e comparações com outras técnicas de equações de evolução não lineares, como a análise de Painlevé e a técnica bilinear de Hirota.