Aproximación numérica de ecuaciones diferenciales parciales estocásticas no lineales de cuarto orden con ruido blanco aditivo espacio-temporal

Información Básica

• ID del Artículo: 2501.18240
• Título: Numerical approximation of nonlinear fourth-order SPDEs with additive space-time white noise
• Autores: Dirk Blömker, Chengcheng Ling, Johannes Rimmele
• Clasificación: math.NA cs.NA math.PR
• Fecha de Publicación: 14 de octubre de 2025
• Enlace del Artículo: https://arxiv.org/abs/2501.18240

Resumen

Este artículo considera la aproximación numérica fuerte de ecuaciones diferenciales parciales estocásticas no lineales de cuarto orden impulsadas por ruido blanco espacio-temporal en el toro bidimensional. Se utiliza un esquema de Galerkin espectral en el espacio y un esquema de Euler en el tiempo para la discretización completamente discreta. Se demuestra la convergencia que alcanza casi la tasa de convergencia espacial 1 y la tasa de convergencia temporal 1 mediante la técnica de costura estocástica.

Antecedentes de Investigación y Motivación

Contexto del Problema

1. Importancia de las ecuaciones diferenciales parciales estocásticas de cuarto orden: La ecuación de Cahn-Hilliard generalizada estudiada en este artículo proviene del trabajo pionero de Cahn e Hilliard sobre la separación de fases en aleaciones binarias, y se utiliza en ciencia de materiales para modelar cómo el campo de concentración evoluciona en el tiempo formando diferentes regiones de fase.
2. Necesidad del ruido estocástico: Al agregar fluctuaciones estocásticas o ruido térmico a la ecuación de Cahn-Hilliard determinista, se puede describir con mayor precisión la incertidumbre microscópica presente en los materiales reales.
3. Desafíos de los métodos numéricos: La aproximación numérica de ecuaciones diferenciales parciales estocásticas de cuarto orden enfrenta desafíos de limitaciones en las tasas de convergencia, particularmente existe una barrera de orden 1/4 en el tamaño del paso temporal.

Motivación de la Investigación

1. Superar la barrera de la tasa de convergencia: Superar la limitación de la tasa de convergencia de orden 1/4 en el tamaño del paso temporal y lograr tasas de convergencia más altas.
2. Perfeccionamiento teórico: Proporcionar una teoría de análisis numérico rigurosa para ecuaciones diferenciales parciales estocásticas de cuarto orden.
3. Aplicaciones prácticas: Proporcionar métodos de resolución numérica efectivos para modelos físicos relacionados.

Contribuciones Principales

1. Superación de la barrera de la tasa de convergencia: Por primera vez se logra una tasa de convergencia espacial casi de orden 1 y una tasa de convergencia temporal de orden 1 en ecuaciones diferenciales parciales estocásticas de cuarto orden, superando la barrera de convergencia temporal de orden 1/4 que existe en la literatura.
2. Aplicación innovadora de la técnica de costura estocástica: Se aplica exitosamente el lema de costura estocástica, originalmente utilizado para ecuaciones diferenciales estocásticas singulares, al análisis numérico de ecuaciones diferenciales parciales estocásticas de cuarto orden.
3. Marco teórico de análisis completo: Se proporciona un análisis de convergencia completo desde la discretización espacial de Galerkin espectral hasta la discretización temporal de Euler, incluyendo estimaciones de regularidad detalladas.
4. Sistematicidad de los métodos técnicos: Se combinan estimaciones de regularidad del semigrupo generado por $-\Delta^2$ con el lema de costura estocástica, logrando un ajuste preciso de la regularidad espacio-temporal de la solución y el ruido.

Explicación Detallada del Método

Definición del Problema

Se estudia la ecuación diferencial parcial estocástica de cuarto orden tipo Cahn-Hilliard generalizada en el toro bidimensional $T^2$: $\partial_t u = -\delta\Delta^2 u - G(u) + \sigma\xi$

donde:

• $G: \mathbb{R} \to \mathbb{R}$ es un término no lineal que satisface $\|G\|_{\infty}, \|\partial G\|_{\infty} < \infty$
• $\xi$ es ruido blanco espacio-temporal, con intensidad controlada por $\sigma > 0$
• $u$ satisface la condición de conservación de masa: $\int_{T^2} u(t,x)dx = 0$

Diseño del Esquema Numérico

1. Discretización Espacial: Método de Galerkin Espectral

Se utilizan modos de Fourier $e_k(x)$ como funciones base, donde: