Interface problems of mixed spatial order

Resumen Técnico: Problemas de Interfaz de Orden Mixto mediante el Método de la Transformada Unificada

Planteamiento del Problema
Este trabajo aborda el problema de valor inicial en la frontera (IBVP) para ecuaciones diferenciales parciales (PDE) de evolución lineal definidas en un dominio espacial dividido en dos semirrectas, $(-\infty, 0)$ y $(0, \infty)$, separadas por una interfaz en $x=0$. A diferencia de estudios previos donde las ecuaciones gobernantes a ambos lados de la interfaz eran idénticas o diferían solo por coeficientes escalares, este artículo investiga problemas de interfaz de "orden mixto". Específicamente, considera la interacción entre ecuaciones de evolución lineal de órdenes espaciales y comportamientos físicos fundamentalmente diferentes (dispersivos frente a disipativos).

Los ejemplos prototípicos analizados involucran:

1. Semirrecta Izquierda ($x < 0$): La ecuación de Schrödinger lineal (lS), la ecuación del calor (h), o la ecuación de Schrödinger biarmónica (biS).
2. Semirrecta Derecha ($x > 0$): La ecuación de Airy o la ecuación de Korteweg de Vries linealizada (lKdV).

Estos sistemas están acoplados mediante condiciones de interfaz en $x=0$. El artículo considera dos conjuntos de condiciones:

• Condición (IfC.1): Continuidad de la solución y de su primera derivada espacial ($C^1$ continuidad).
• Condición (IfC.2): Continuidad de la solución y de sus primeras dos derivadas espaciales ($C^2$ continuidad), requerida para la ecuación de Schrödinger biarmónica de cuarto orden.

Los valores de la interfaz (datos de frontera en $x=0$) no están prescritos a priori; son funciones desconocidas determinadas por el acoplamiento de las dos PDE distintas.

Metodología
Los autores emplean el Método de la Transformada Unificada (UTM), también conocido como el método de Fokas, que sirve como una generalización de la transformada de Fourier para IBVPs en dominios con fronteras. El procedimiento de solución implica los siguientes pasos:

1. Relaciones Globales: Para cada problema de semirrecta, los autores derivan "relaciones globales" (identidades espectrales) que vinculan los datos iniciales, la solución y los valores de frontera desconocidos en el plano espectral complejo.
2. Cambio de Variables: Para manejar los órdenes mixtos, los autores introducen funciones específicas de tipo raíz ($\sigma, \rho, \tau, \gamma$) y un biholomorfismo $\nu$ (para el caso de lKdV), que mapean el parámetro espectral $\lambda$ al plano complejo $k$ de manera que alineen las relaciones de dispersión de las diferentes ecuaciones (por ejemplo, relacionando $k^2, k^3$ y $k^4$).
3. Construcción del Sistema Lineal: Al aplicar estos cambios de variables a las relaciones globales, los autores construyen un sistema lineal (típicamente de $3\times3$ o $4\times4$) para las transformadas de los valores de interfaz desconocidos.
4. Representación de la Solución: Resolver este sistema lineal permite la eliminación de los valores de interfaz desconocidos. Los autores luego sustituyen estas expresiones en las representaciones integrales de las soluciones para cada semirrecta.
5. Deformación de Contornos y Asintótica: Una contribución técnica crítica es la justificación rigurosa para deformar los contornos de integración (por ejemplo, desde el eje real hacia contornos como $\tilde{R}$ en el semiplano inferior) para evitar singularidades (específicamente en $\lambda=0$) y para eliminar las integrales que involucran los términos de la solución dependientes del tiempo ($u(x,t), v(x,t)$). Esto se logra mediante argumentos personalizados que utilizan el comportamiento asintótico específico de las funciones raíz y los núcleos exponenciales, en lugar del lema de Jordan estándar, debido al desajuste en los términos de orden principal.

Contribuciones Clave y Resultados
El artículo proporciona fórmulas de solución explícitas para cuatro problemas de interfaz distintos, presentadas en las Proposiciones 1 a 4:

• Proposición 1: Resuelve la interfaz entre la ecuación de Schrödinger lineal y la ecuación de Airy bajo condiciones de interfaz $C^1$.
• Proposición 2: Resuelve la interfaz entre la ecuación de Schrödinger lineal y la ecuación de KdV linealizada (incluyendo un término de orden inferior $av_x$) bajo condiciones $C^1$. Esto requiere la introducción del biholomorfismo $\nu$ para manejar la relación de dispersión no homogénea.
• Proposición 3: Resuelve la interfaz entre la ecuación del calor (disipativa) y la ecuación de Airy (dispersiva) bajo condiciones $C^1$.
• Proposición 4: Resuelve la interfaz entre la ecuación de Schrödinger biarmónica (dispersiva de cuarto orden) y la ecuación de Airy bajo condiciones $C^2$. Este caso involucra un sistema lineal de $4\times4$ debido al orden superior de la PDE.

Las fórmulas explícitas se expresan mediante integrales de contorno que involucran las transformadas de los datos iniciales y funciones núcleo específicas derivadas de los mapeos de tipo raíz.

Evaluación Numérica
Para demostrar la utilidad de las fórmulas derivadas, los autores realizan evaluaciones numéricas utilizando la cuadratura de Filon. Este método es elegido para manejar la naturaleza altamente oscilatoria de los integrandos. El artículo presenta gráficos de las soluciones para diversas configuraciones de datos iniciales:

• Schrödinger-Airy: Muestra la transición de ondas dispersivas desde el régimen unidireccional de Airy hacia el régimen bidireccional de Schrödinger, y viceversa, resaltando los patrones de tipo "ruleta" y los efectos de arrastre causados por la interfaz.
• Biharmonic Schrödinger-Airy: Ilustra interacciones de ondas similares pero con una dinámica gobernada por una relación de dispersión cuártica, lo que conduce a la aparición más rápida de patrones complejos.
• Heat-Airy: Demuestra la interacción entre ondas difusivas y dispersivas, mostrando un decaimiento exponencial en el lado del calor y dispersión en el lado de Airy.

Significancia
El artículo afirma su importancia al extender el Método de la Transformada Unificada a una nueva clase de problemas donde las ecuaciones gobernantes a ambos lados de una interfaz son de tipos fundamentalmente diferentes (diferentes órdenes espaciales y naturalezas físicas). Las aplicaciones previas del UTM a problemas de interfaz se limitaron en gran medida a casos donde las ecuaciones eran idénticas o diferían solo por coeficientes escalares.

Al resolver con éxito estos problemas de orden mixto, los autores demuestran la robustez del UTM para manejar estructuras espectrales complejas que surgen de relaciones de dispersión dispares. Este trabajo proporciona el primer marco de solución analítica para tales interacciones, las cuales son relevantes para modelar fenómenos físicos que involucran medios heterogéneos, como el crecimiento de tumores en el tejido cerebral (como se cita en la introducción respecto a los gliomas) o materiales compuestos con propiedades físicas variables. Las fórmulas explícitas derivadas sirven como base para análisis futuros, incluyendo posibles extensiones a problemas no lineales mediante técnicas de mapeo de contracción.