A Model-Based Digital Twin for Predicting Deterioration of Mogao Polychrome Sculpture

Resumen Técnico: Un Gemelo Digital Basado en Modelos para la Predicción del Deterioro de la Escultura Policroma de Mogao

Planteamiento del Problema
Las Grutas de Mogao contienen más de 2.400 esculturas de arcilla policromas cuyas superficies pintadas son susceptibles de deterioro impulsado por cambios en la temperatura y la humedad. Si bien existen modelos predictivos para pinturas murales planas y otros materiales patrimoniales (p. ej., celulosa, pintura al óleo, gesso), persiste una brecha fundamental para las esculturas de arcilla exentas de soporte. Estos objetos poseen un sistema de materiales heterogéneo (arcilla con paja, yeso reforzado con fibras, pigmentos minerales unidos con cola animal) que responde de manera distinta al ciclado higrotérmico en comparación con las pinturas de pared. Además, las implementaciones actuales de gemelos digitales en contextos patrimoniales se centran principalmente en la integridad estructural, el rendimiento energético o registros 3D estáticos, careciendo de modelos basados en la física que predigan la degradación de los materiales (específicamente el cambio de color, el moho, el daño por sales, el agrietamiento y el envejecimiento) a nivel de objeto. No existe una función de dosis-respuesta validada experimentalmente que vincule directamente el cambio de color ($\Delta E^*_{ab}$) con la temperatura y la humedad relativa para pigmentos minerales sobre sustratos de arcilla.

Metodología
Los autores presentan un marco de trabajo de gemelo digital basado en modelos que integra cinco mecanismos de deterioro basados en la física sobre una única malla fotogramétrica 3D de una escultura. El sistema se construye utilizando Ingeniería Dirigida por Modelos (MDE) para asegurar que el marco sea genérico y reutilizable.

1. Arquitectura del Marco:

   • Tiempo de Diseño: Un lenguaje de modelado específico de dominio (metamodelo Ecore/EMF) define el grafo de entidades (Cuevas, Estatuas, Sensores). Plantillas Epsilon EGL (aproximadamente 595 líneas) generan automáticamente la capa de datos del backend (esquemas, servicios y controladores Mongoose/Express) para 13 tipos de entidades.
   • Tiempo de Ejecución: El sistema ingiere telemetría de sensores en vivo (Temperatura $T$, Humedad Relativa $RH$, Iluminancia $I$) y ejecuta cinco núcleos de deterioro independientes.
   • Visualización: Un visor 3D basado en navegador renderiza los resultados en dos modos: Simulación Ambiental (escenarios de "qué pasaría si" definidos por el usuario) y Monitoreo Ambiental (reproducción de datos en vivo).

2. Los Cinco Núcleos de Deterioro:

   • Modelo 1 (Segmentación de Pigmentos): Un clasificador determinista de umbral HSV mapea la textura 3D a ocho clases mineralógicas (p. ej., cinabrio, azurita, malaquita, pan de oro) para dirigir la cinética por pigmento.
   • Modelo 2 (Desvanecimiento Químico): Un modelo de dosis-respuesta acoplado a la sorción que combina el escalado de temperatura de Arrhenius con una isoterma de humedad de Paltakari–Karlsson. Calcula la degradación de color por píxel ($\Delta E^*_{ab}$) basada en energías de activación específicas para vías oscuras y claras.
   • Modelo 3 (Multiplicador de Vida Útil): Utiliza el multiplicador de vida útil basado en Arrhenius de Michalski para reescalar la vida de servicio esperada en relación con un entorno de referencia.
   • Modelo 4 (Crecimiento de Moho): Implementa el modelo VTT Hukka–Viitanen para predecir el índice de moho ($M$) basado en umbrales de humedad crítica dependientes de la temperatura.
   • Modelo 5 (Cristalización de Sales): Utiliza el marco termodinámico de Steiger para modelar la presión de cristalización del sulfato de sodio. Crucialmente, tiene en cuenta el entorno sub-deliquescente de las cuevas de Mogao, tratando el daño por sal como una presión estática sostenida en lugar de un ciclo de secado-humedificación, regulado por un factor de disponibilidad capilar ($\alpha_s$).
   • Modelo 6 (Fatiga Higromecánica): Aplica la ley de Basquin y la regla de Miner para modelar el daño por fatiga ($D$) derivado de la amplitud cíclica de la humedad ($\Delta RH$).

3. Evaluación de Riesgo Compuesto:
   Los resultados de los cinco modelos se normalizan a una escala de 0 a 1 y se agregan en un índice compuesto ($R_{composite}$) utilizando un enfoque de máximo riesgo. Este índice dirige una cola de mantenimiento clasificada y superposiciones visuales en la malla 3D.

4. Piloto Experimental (Envejecimiento Acelerado):
   Para fundamentar el modelo de desvanecimiento, los autores fabricaron réplicas a escala 1:5 de un pedestal de loto de Mogao utilizando materiales tradicionales (arcilla chengban, paja, cáñamo, cola animal y pigmentos minerales). Estas fueron envejecidas durante 61 días bajo dos condiciones: un brazo de "sala" (23°C, 40% RH, luz ambiental) y un brazo de "cámara" (40°C, 10% RH, oscuridad). Las mediciones colorimétricas ($\Delta E^*_{ab}$) se tomaron en 21 puntos temporales.

Resultados Clave

• Cinética del Piloto: La campaña de envejecimiento acelerado midió los contrastes de la tasa de desvanecimiento entre los dos brazos. El cinabrio, la azurita y la malaquita envejecieron de 1.5 a 2.0 veces más rápido en el brazo de la sala (más fresca, húmeda y con luz) en comparación con el brazo de la cámara (caliente, seca y oscura). Debido a las variables de confusión (la luz, la temperatura y la humedad variando simultáneamente), el estudio reporta estas relaciones de tasas en lugar de las energías de activación deconvolucionadas.
• Aplicación de Caso de Estudio: El marco se aplicó al Bodhisattva de Rodillas (Mogao Cueva 328).
  • Clima Actual: Bajo el clima monitoreado de las cuevas de Mogao (media de 12°C, RH < 80%), la escultura se encuentra en un régimen benigno. El crecimiento de moho está latente (índice = 0) y la fatiga higromecánica es insignificante. La cristalización de sales ejerce una presión estática sostenida pero está espacialmente limitada a la base alimentada por capilaridad.
  • Escenarios Contrafactuales: El gemelo demuestra que si se pierde el control climático (p. ej., la RH aumenta al 85–96% y la temperatura aumenta), el riesgo de cristalización de sales y de moho se activa rápidamente. Un escenario de "Almacenamiento Deficiente" proyecta que el índice de riesgo compuesto subirá a niveles de falla ($R=1.0$) dentro de cuatro décadas, mientras que un amortiguamiento activo mantiene la línea base de bajo riesgo.
• Visualización: El sistema logra renderizar firmas de deterioro distintas (p. ej., desvanecimiento de pigmento no uniforme, eflorescencia de sal en la base, redes de craquelado) directamente sobre la malla 3D, permitiendo a los conservadores inspeccionar los mecanismos de forma individual o compuesta.

Significancia y Reivindicaciones
El artículo reivindica cuatro contribuciones primarias:

1. Marco Genérico: Introduce un enfoque de Ingeniería Dirigida por Modelos donde el backend se genera a partir de un modelo de dominio específico, lo que hace que el gemelo digital sea adaptable a diferentes objetos patrimoniales sin necesidad de reescribir la capa de datos central.
2. Renderizado de "Qué pasaría si": Permite la visualización interactiva y basada en la física de las firmas de deterioro en una malla 3D, permitiendo a los conservadores simular trayectorias ambientales y observar cambios de textura en tiempo real.
3. Integración de Múltiples Mecanismos: Operacionaliza cinco núcleos de degradación física distintos en un solo objeto, haciendo que los procesos heterogéneos sean comparables y reduciéndolos a una cola de mantenimiento derivada de la física y clasificada.
4. Cinética del Piloto: Proporciona los primeros contrastes de tasa de desvanecimiento controlados para los sistemas de arcilla policroma de Mogao, estableciendo que las condiciones más frescas, húmedas y con luz pueden ser más dañinas para ciertos pigmentos que las condiciones calientes, secas y oscuras, incluso si el cambio de color absoluto es pequeño durante periodos cortos.

Limitaciones y Modestia
Los autores declaran explícitamente que el piloto actual no proporciona energías de activación independientes y deconvolucionadas debido a la confusión de la luz y la temperatura en el diseño de dos brazos. Las proyecciones de largo horizonte (p. ej., simulaciones de 200 años) están calibradas para coincidir con los patrones de degradación arqueológica (aprox. 48% de pérdida de color en 1200 años) en lugar de ser validadas contra registros independientes de varias décadas. Los autores reconocen que la extrapolación de tasa lineal del piloto de 61 días muestra un bajo rendimiento fuera de la muestra ($Q^2 = -16.7$), lo que requiere el uso de un modelo de saturación para las proyecciones a largo plazo. Identifican la necesidad de diseños factoriales de luz emparejada y la validación externa contra datos de archivo como prioridades para trabajos futuros.