Well-posedness of a first-order formulation for fractionally damped nonlinear acoustics

Resumen Técnico: Bien posed de una formulación de primer orden para la acústica no lineal con amortiguamiento fraccionario

Planteamiento del problema
Este trabajo aborda la bien posed de un sistema hiperbólico cuasilineal de primer orden en el tiempo con memoria, que modela ondas acústicas no lineales con amortiguamiento fraccionario. El sistema general dado es:
$$\begin{cases} A_0 u_t + \sum_{j=1}^d A_j(u, x) u_{x_j} - J * \Delta u = g & \text{en } (0, T) \times \Omega, \\ u(0) = u_0 & \text{en } \Omega, \\ u = 0 & \text{en } (0, T) \times \partial\Omega, \end{cases}$$
donde $u = (p, v)^\top$ representa la presión escalar y el vector de velocidad. El término $J * \Delta u$ denota una convolución temporal con un núcleo de memoria $J$, capturando efectos de amortiguamiento fraccionario. El sistema engloba formulaciones de primer orden de las ecuaciones clásicas de tipo Kuznetsov y Westervelt.

Una dificultad primaria en el análisis de este sistema es que la matriz de flujo $A_j$ es generalmente no simétrica, una propiedad típicamente requerida para la bien posed de sistemas hiperbólicos de primer orden. Además, para núcleos integrables generales $J \in L^1$, el término de amortiguamiento no local no proporciona la mayor regularidad (suavizado parabólico) disponible en los casos donde $J = \delta_0$ (medida de Dirac). Consecuentemente, las estimaciones parabólicas estándar fallan, lo que requiere un enfoque hiperbólico que dependa de estimaciones de energía sin asumir la regularidad inducida por el Laplaciano.

Metodología
Los autores emplean una estrategia de simetrización combinada con nuevas estimaciones de coercividad no lineal para los términos de memoria.

1. Simetrización: Para superar la no simetría de la matriz de flujo, el sistema se multiplica por una matriz multiplicativa adecuadamente elegida $S(u)$. Para el caso de tipo Westervelt ($\mu = \kappa = 0$), se utiliza un simetrizador diagonal. Para el caso de tipo Kuznetsov en el límite de inviscididad ($J=0$), se construye un simetrizador no diagonal para transformar el sistema en un sistema hiperbólico simétrico cuasilineal.
2. Supuestos sobre el núcleo: El análisis se basa en que el núcleo de memoria $J$ sea completamente monótono (Supuesto 1). Esta clase incluye núcleos de Abel, Mittag-Leffler y exponenciales. Se introduce un supuesto más fuerte (Supuesto 2) respecto a la regularidad del resolvente de $J$ para relajar las condiciones de pequeñez de los datos.
3. Coercividad no lineal: Una contribución técnica central es la derivación de estimaciones de coercividad no lineal para el término de memoria. Específicamente, los autores demuestran que, para una matriz $V$ diagonal y definida positiva y un núcleo completamente monótono $J$, la desigualdad
   $$\int_0^t (V(s)(J*y)(s), y(s)) \, ds \geq C_J \int_0^t \|(J*y)(s)\|^2 \, ds$$
   se cumple. Esto permite el control del término de amortiguamiento no local dentro del marco de energía.
4. Argumento de punto fijo: La bien posed del sistema no lineal se establece mediante el teorema del punto fijo de Banach aplicado a un problema linealizado. El problema linealizado se resuelve utilizando aproximaciones de Galerkin y estimaciones de energía de alto orden en espacios de Sobolev ($H^4$ y $H^5$). El mapeo de contracción se realiza en una norma más débil ($L^\infty(0, T; L^2)$) mientras que las cotas a priori se establecen en normas más fuertes, una técnica estándar para problemas hiperbólicos.

Resultados Clave
El artículo establece dos tipos principales de resultados de bien posed:

• Sistema de tipo Westervelt (Dominios acotados, $d \leq 3$):

  • Bajo el Supuesto 1 (núcleo completamente monótono), el sistema admite una solución única local en el tiempo para datos iniciales y términos fuente pequeños. La prueba requiere una condición de pequeñez de los datos relativa a la derivada temporal del simetrizador.
  • Bajo el Supuesto 2 (regularidad adicional en el resolvente de $J$), la condición de pequeñez de los datos se relaja significativamente. Específicamente, la presión inicial solo necesita estar por debajo de un umbral explícito definido por el parámetro de no linealidad, en lugar de requerir que todo el conjunto de datos sea pequeño.
  • Estos resultados se aplican a dominios acotados con fronteras $C^{4,1}$ y condiciones de Dirichlet homogéneas.

• Sistema de tipo Kuznetsov (Caso inviscid, $J=0$):

  • Para el espacio completo $\Omega = \mathbb{R}^d$ y dimensiones arbitrarias, se demuestra que el caso de inviscididad ($J=0$) es localmente bien posed utilizando argumentos hiperbólicos estándar. Este resultado se mantiene sin la restricción $\mu = \kappa = 0$ requerida para el análisis de tipo Westervelt, siempre que el simetrizador produzca autovalores positivos.

Significación y Reivindicaciones
Los autores posicionan este trabajo como una forma de llenar un vacío en el análisis matemático de los modelos acústicos de tiempo fraccionario de primer orden. Mientras que la bien posed para ecuaciones de segundo orden (Westervelt, Kuznetsov) y tercer orden (Jordan–Moore–Gibson–Thompson) con amortiguamiento fraccionario ha sido establecida, existen menos resultados para las formulaciones de primer orden.

El artículo reivindica su importancia en las siguientes áreas:

• Novedad del enfoque: Proporciona los primeros resultados de bien posed para modelos acústicos no lineales de tiempo fraccionario de primer orden utilizando estimaciones de energía que no dependen del suavizado parabólico.
• Cobertura de núcleos generales: El análisis cubre una amplia clase de núcleos físicamente relevantes (Abel, Mittag-Leffler, exponenciales) a través del supuesto de monotonía completa.
• Base para la numeración: Los autores afirman que la formulación de primer orden es particularmente adecuada para el desarrollo de esquemas numéricos eficientes. La bien posed establecida proporciona una base teórica sólida para tales esquemas y aplicaciones potenciales como la identificación de parámetros.

El trabajo es modesto en su alcance, señalando que los resultados son locales en el tiempo y, para el caso de Westervelt, restringidos a datos pequeños (a menos que se cumpla el Supuesto 2) y bajas dimensiones espaciales ($d \leq 3$) debido a los requisitos de incrustación de Sobolev. Se identifica como trabajo futuro la posibilidad de relajar el supuesto de matriz diagonal en los lemas de coercividad para extender los resultados al sistema completo de Kuznetsov con memoria, e investigar comportamientos límite a medida que los núcleos desaparecen.