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

Resumen Técnico: Solución Exacta a una Clase de Problemas para la Ecuación de Burgers en Intervalos Acotados

Enunciado del Problema
El artículo aborda la ecuación de Burgers viscosa, una ecuación diferencial parcial no lineal que modela el equilibrio entre convección y difusión, definida en intervalos acotados genéricos $[l_1, l_2]$. El problema específico consiste en determinar la función $w(x, t)$ sujeta a:

1. Ecuación Gobernante: $w_t - a^2 w_{xx} + w w_x = 0$, donde $a^2$ representa el coeficiente de viscosidad cinemática ($1/Re$).
2. Condición Inicial: Un perfil inicial no homogéneo $w(x, 0) = w_0(x)$.
3. Condiciones de Contorno: Condiciones de contorno de Dirichlet fijas $w(l_1, t) = \alpha_1$ y $w(l_2, t) = \alpha_2$.

Los autores señalan que, aunque los métodos numéricos (diferencias finitas, elementos finitos, métodos espectrales) y las técnicas semianalíticas (descomposición de Adomian, perturbación de homotopía) son prevalentes, a menudo carecen de soluciones en forma cerrada o requieren recursos computacionales significativos. El objetivo es derivar una solución operacional exacta que ofrezca una expresión en forma cerrada en el dominio de Laplace y una solución precisa en el dominio del tiempo mediante transformadas inversas.

Metodología
La estrategia de solución propuesta se basa en un proceso de transformación de dos pasos:

1. Transformación de Hopf-Cole: La ecuación de Burgers no lineal se linealiza utilizando la transformación clásica de Hopf-Cole:
   $$w(x, t) = -2a^2 \frac{u_x(x, t)}{u(x, t)}$$
   Esto convierte el problema no lineal en un problema de valor en la frontera lineal para la ecuación de difusión de calor unidimensional (difusión-reacción con término de reacción cero):
   $$u_t - a^2 u_{xx} = 0$$
   La condición inicial para $u$ se deriva de la integral de la condición inicial $w_0$, y las condiciones de Dirichlet para $w$ se transforman en condiciones de contorno de Robin homogéneas para $u$.

2. Transformada de Laplace y Solución Operacional: La ecuación de calor lineal se resuelve en el dominio de Laplace (dominio $p$). Los autores utilizan una solución operacional exacta previamente establecida para ecuaciones lineales de difusión-reacción con coeficientes constantes en intervalos finitos.

   • El problema se formula como un problema de valor en la frontera de dos puntos en el dominio de Laplace.
   • Se deriva una expresión exacta para la transformada de Laplace $U(x, p)$, que involucra términos exponenciales dependientes de $\sqrt{p}$ y un sistema lineal de cuatro ecuaciones para determinar los valores en la frontera ($U(l_1, p), U_x(l_1, p), U(l_2, p), U_x(l_2, p)$).
   • El término fuente en la ecuación linealizada es cero, lo que simplifica la solución operacional para que dependa de la transformada de la condición inicial $R(x, p)$.

3. Recuperación en el Dominio del Tiempo: La solución exacta para la variable original $w(x, t)$ se obtiene implícitamente mediante la transformada inversa de Laplace:
   $$w(x, t) = -2a^2 \frac{\mathcal{L}^{-1}\{U_x(x, p)\}}{\mathcal{L}^{-1}\{U(x, p)\}}$$
   El artículo señala que, aunque las inversas analíticas pueden existir en forma cerrada (potencialmente mediante transformadas de Mellin), se emplean algoritmos numéricos altamente eficientes (específicamente el algoritmo de de Hoog et al. implementado como invlap.m en Matlab) para calcular la solución en el dominio del tiempo, independientemente de la complejidad de las expresiones de Laplace.

Contribuciones y Resultados Clave
El artículo presenta dos ejemplos ilustrativos para validar la metodología:

• Ejemplo 1: Un problema con una condición inicial $w_0(x) = \frac{2a^2\pi \sin(\pi x)}{\sigma + \cos(\pi x)}$ y condiciones de contorno nulas. Los autores derivaron la solución en el dominio de Laplace y aplicaron la transformada inversa. La solución resultante en el dominio del tiempo coincidió perfectamente con la solución clásica exacta en forma cerrada conocida (Wood, 2014).
• Ejemplo 2: Un problema con una condición inicial $w_0(x) = \sin(\pi x)$. Se sabe que la solución exacta clásica para este caso es una serie infinita. El método de los autores produjo una solución exacta implícita en el dominio de Laplace. Cuando se invirtió numéricamente, los resultados se alinearon estrechamente con la solución de la serie infinita truncada ("Solución Teórica").

Rendimiento Comparativo
El estudio destaca dos ventajas principales del método propuesto sobre los enfoques existentes:

1. Eficiencia Computacional: En el Ejemplo 2, se encontró que el procedimiento numérico de inversa de Laplace fue más de 100 veces más rápido en tiempo de ejecución en comparación con la evaluación de la solución de la serie infinita truncada.
2. Precisión: Las soluciones de inversa de Laplace numérica mostraron errores despreciables en comparación con las soluciones de series teóricas. Por el contrario, las soluciones de diferencias finitas (calculadas mediante pdepe() de Matlab) mostraron una discrepancia creciente a medida que disminuía el coeficiente de viscosidad, particularmente en los primeros pasos de tiempo.

Significado y Afirmaciones
Los autores afirman que su trabajo proporciona una solución operacional en forma cerrada en el dominio de Laplace y una solución exacta implícita en el dominio del tiempo. Sostienen que este enfoque ofrece una "alternativa novedosa" a las soluciones basadas en series infinitas o en discretización puramente numérica.

El significado del trabajo se enmarca de la siguiente manera:

• Punto de Referencia: La solución exacta derivada sirve como un referente robusto para validar y probar otros métodos numéricos y semianalíticos.
• Eficiencia: El método es computacionalmente superior a las expansiones en series, que pueden ser intensivas en recursos.
• Generalidad: La metodología maneja tanto soluciones suaves como soluciones generalizadas con discontinuidades en las derivadas.

El artículo concluye modestamente, indicando que la metodología de solución será una referencia para trabajos futuros, incluidas posibles extensiones a ecuaciones de Burgers forzadas en dimensiones superiores y comparaciones con otras técnicas de ecuaciones de evolución no lineales como el análisis de Painlevé y la técnica bilineal de Hirota.