Safe Trajectory Tracking of the Stefan Problem with Second-Order Moving Boundary Dynamics

Resumen Técnico: Seguimiento de Trayectoria Seguro del Problema de Stefan con Dinámica de Frontera Móvil de Segundo Orden

Planteamiento del Problema
Este artículo aborda el control de seguimiento de trayectorias seguro del problema de Stefan de una fase, un modelo clásico para fenómenos de cambio de fase térmica como la fusión y la solidificación. El desafío específico que se aborda es el control de un sistema donde la frontera móvil (la interfaz entre fases) está gobernada por dinámicas de segundo orden, que representan la inercia térmica, en lugar de las dinámicas de primer orden más comunes. El sistema se modela mediante una Ecuación Diferencial Parcial (EDP) parabólica para el perfil de temperatura definido en un dominio espacial variable en el tiempo, acoplada con una Ecuación Diferencial Ordinaria (EDO) de segundo orden para la posición de la frontera.

El objetivo de control es diseñar una entrada de flujo de calor en la frontera, $q_c(t)$, que impulse la frontera móvil y el perfil de temperatura del sistema para seguir una trayectoria de referencia deseada, $s_r(t)$ y $T^r(x,t)$. Crucialmente, el controlador debe garantizar la seguridad, lo que significa que las restricciones físicas del modelo deben satisfacerse en todo momento: la temperatura debe permanecer por encima del punto de fusión ($T(x,t) \ge T_m$) para evitar que la fase líquida se congele en la entrada de calor, y la posición de la frontera debe permanecer dentro del dominio físico ($0 < s(t) < L$).

Metodología
La estrategia de control propuesta combina un componente de alimentación directa con un controlador de retroalimentación diseñado mediante conformación de energía y retroceso (backstepping) de EDP.

1. Control de Alimentación Directa (Problema de Stefan Inverso):
   Para generar la trayectoria de referencia, los autores resuelven el problema de Stefan inverso para un movimiento de frontera de referencia dado $s_r(t)$. Emplean un enfoque de expansión en serie, representando el perfil de temperatura de referencia como una serie de potencias en la coordenada espacial relativa a la frontera móvil. Los coeficientes de esta serie se determinan recursivamente. El artículo demuestra la convergencia de esta serie infinita bajo el supuesto de que la trayectoria de referencia es de orden Gevrey, asegurando la existencia de un perfil de temperatura de referencia bien definido y el flujo de calor de alimentación directa correspondiente $q_c^r(t)$.

2. Control de Retroalimentación (Conformación de Energía y Seguridad):
   Se deriva un controlador de seguimiento de trayectorias utilizando un enfoque de conformación de energía. La energía total del sistema se define para incluir tanto la energía térmica como la energía cinética/potencial de la frontera móvil. La ley de control se diseña como:
   $$q_c(t) = q_c^r(t) - c(E(t) - E_r(t))$$
   donde $E(t)$ y $E_r(t)$ son las energías real y de referencia, respectivamente, y $c$ es una ganancia positiva.
   Para garantizar la seguridad (específicamente, que el flujo de calor permanezca no negativo y la temperatura se mantenga por encima del punto de fusión), los autores utilizan el Principio del Máximo. Imponen suposiciones específicas sobre la energía inicial y la trayectoria de referencia para asegurar que el sistema en lazo cerrado satisfaga las restricciones físicas. Este enfoque actúa efectivamente como un filtro de seguridad, asegurando que incluso si la trayectoria de referencia exige enfriamiento (flujo de calor negativo), la entrada de control real permanezca no negativa, evitando la violación de la condición de cambio de fase.

3. Análisis de Estabilidad:
   La estabilidad exponencial global del error de seguimiento se establece utilizando análisis de Lyapunov combinado con una transformación de retroceso (backstepping) de EDP. El sistema de error de seguimiento se transforma en un sistema objetivo con propiedades de estabilidad deseables. Mediante la construcción de una función de Lyapunov que tenga en cuenta la energía de los estados de error y la dinámica de la frontera móvil, los autores demuestran que el error de seguimiento converge exponencialmente a cero. El análisis maneja cuidadosamente el dominio variable en el tiempo y el acoplamiento entre la EDP y la EDO de segundo orden.

Contribuciones Clave

• Dinámicas de Segundo Orden: El artículo extiende el control del problema de Stefan para incluir dinámicas de frontera móvil de segundo orden, lo que captura mejor la inercia térmica en los procesos de cambio de fase.
• Prueba Rigurosa de Seguridad: A diferencia de trabajos anteriores que a menudo asumían la seguridad o referencias estáticas, este artículo proporciona una prueba rigurosa de que la ley de control propuesta garantiza el cumplimiento de las restricciones físicas (no congelación y límites del dominio) para una trayectoria de referencia dinámica.
• Convergencia de la Solución en Serie: Los autores demuestran la convergencia de la solución de expansión en serie para el problema de Stefan inverso bajo dinámicas de frontera de segundo orden, un paso necesario para generar el control de alimentación directa.
• Estabilidad Exponencial Global: El artículo establece la estabilidad exponencial global para el error de seguimiento en lazo cerrado utilizando una combinación novedosa de conformación de energía, retroceso (backstepping) y análisis de Lyapunov, sin depender de linealizaciones o suposiciones simplificadas sobre la dinámica del error.

Resultados
Se realizaron simulaciones numéricas utilizando parámetros para la fusión de Zinc. La trayectoria de referencia se diseñó como una función sinusoidal no trivial que decae con el tiempo, asegurando que la frontera se mueva monótonamente pero con velocidad variable.

• Rendimiento de Seguimiento: El sistema en lazo cerrado siguió con éxito la posición de la frontera de referencia y el perfil de temperatura, con el error de seguimiento convergiendo a cero dentro del tiempo de simulación (aproximadamente 60 minutos).
• Verificación de Seguridad: La simulación confirmó que, aunque el flujo de calor de referencia ocasionalmente se volvió negativo (indicando una necesidad teórica de enfriamiento), la entrada de control real permaneció no negativa, satisfaciendo la restricción de seguridad. El perfil de temperatura se mantuvo estrictamente por encima del punto de fusión durante toda la simulación.
• Comportamiento Transitorio: Los resultados mostraron que el controlador aplicó inicialmente un flujo de calor grande para compensar la menor energía inicial del sistema en comparación con la referencia, lo que llevó a un breve sobrepaso antes de asentarse en un seguimiento preciso.

Significado
Los autores presentan este trabajo como el primero en proporcionar una prueba rigurosa tanto de seguridad como de estabilidad de seguimiento para el problema de Stefan con control de seguimiento de trayectorias. Al abordar la no linealidad geométrica y el acoplamiento entre la EDP y la EDO de segundo orden, el artículo ofrece un diseño de control sistemático que garantiza la validez física de la solución (seguridad) mientras logra un seguimiento preciso. La metodología cierra la brecha entre el diseño de control teórico para EDPs y la necesidad práctica de satisfacer restricciones físicas rígidas en aplicaciones de cambio de fase. El código para el método propuesto se pone a disposición para facilitar investigaciones y validaciones futuras.