A few notes about viscoplastic rheologies

El artículo utiliza herramientas rigurosas de análisis convexo para examinar combinaciones en serie y paralelo de viscosidad lineal y plasticidad perfecta, con el objetivo de sintetizar potenciales de disipación viscoplásticos y comparar modelos viscosos en serie rigurosos con modelos empíricos basados en medias armónicas.

Lo esencial

¡Hola! Vamos a desglosar este artículo científico de una manera sencilla, como si estuviéramos contando una historia sobre cómo se comportan las cosas cuando las empujamos o estiramos.

Imagina que el mundo está lleno de materiales: desde el hielo de un glaciar hasta la roca fundida bajo la Tierra, pasando por la masa de pan o la miel. Todos estos materiales tienen una personalidad especial: a veces se comportan como sólidos rígidos, a veces como líquidos viscosos (como la miel) y a veces como plásticos que se deforman y no vuelven a su forma original.

El autor, Tomás Roubíček, quiere crear un "manual de instrucciones" matemático para predecir exactamente cómo se comportarán estos materiales cuando los sometemos a fuerzas.

Aquí tienes la explicación paso a paso:

1. Los dos personajes principales: El "Gelatina" y el "Chicle"

Para entender los materiales complejos, el autor los descompone en dos tipos básicos de comportamiento:

• El "Gelatina" (Viscosidad): Imagina un fluido espeso como la miel o el aceite. Si lo mueves lento, ofrece resistencia. Si lo mueves rápido, ofrece más resistencia. Es como un amortiguador de coche: cuanto más rápido intentas moverlo, más se opone. En física, esto es la viscosidad.
• El "Chicle" o "Goma Elástica Rota" (Plasticidad Perfecta): Imagina un trozo de chicle o arcilla. Si lo empujas un poquito, no pasa nada. Pero si superas cierta fuerza (un umbral), se deforma y nunca vuelve a su forma original. No importa cuánto empujes después, ya se rompió su "resistencia inicial". Esto es la plasticidad.

2. El gran problema: ¿Cómo se juntan?

El artículo se centra en una pregunta clave: ¿Qué pasa si un material tiene ambas propiedades a la vez? (Como el hielo de un glaciar, que es duro pero fluye lentamente, o la roca del manto terrestre).

El autor explora dos formas principales de combinar estos dos "personajes":

Opción A: En Paralelo (El equipo de trabajo)

Imagina que tienes un amortiguador de gelatina y un bloque de chicle pegados uno al lado del otro, y ambos tienen que moverse juntos.

• Lo que pasa: Para mover el sistema, tienes que vencer la resistencia del chicle (si la fuerza es suficiente) Y al mismo tiempo vencer la resistencia de la gelatina.
• Resultado: Es un modelo muy "duro". Se usa mucho para describir fluidos como el fluido de Bingham (como la pasta de dientes: no sale del tubo hasta que le das suficiente presión, y luego fluye).

Opción B: En Serie (La cadena de montaje)

Ahora imagina que el amortiguador y el chicle están uno detrás del otro, como eslabones de una cadena.

• Lo que pasa: La fuerza que aplicas es la misma para ambos, pero la deformación se reparte.
  • Si la fuerza es baja, solo el chicle se resiste y nada se mueve (o se mueve muy poco).
  • Si la fuerza supera el umbral del chicle, este se "rompe" (se activa) y permite que el amortiguador de gelatina haga todo el trabajo de fluir.
• Resultado: Este modelo es más suave. Permite que el material se mueva lentamente incluso con fuerzas pequeñas (creep o fluencia lenta), y luego acelere si la fuerza aumenta. Es ideal para modelar terremotos o el movimiento lento de placas tectónicas.

3. La Magia Matemática: La "Cocina" de las Fórmulas

El autor usa herramientas matemáticas avanzadas (llamadas análisis convexo) para mezclar estas dos propiedades.

• La Metáfora de la Receta: Imagina que tienes dos recetas (una para la gelatina y otra para el chicle). El autor no quiere simplemente ponerlas una al lado de la otra; quiere crear una nueva receta única que describa el comportamiento del material combinado.
• El Truco del "Promedio": A veces, los ingenieros usan una fórmula rápida y aproximada (como un promedio armónico) para mezclar estas propiedades. El autor dice: "Cuidado, esa fórmula rápida a veces miente". Él demuestra que si usas las herramientas matemáticas correctas (llamadas convolución infimal), obtienes una receta exacta y rigurosa.
• El Resultado: A veces, la fórmula exacta es muy compleja (como resolver una ecuación cúbica), pero a veces se simplifica en algo elegante y suave que los ingenieros pueden usar sin miedo a errores.

4. ¿Por qué nos importa esto? (El mundo real)

Este no es solo un ejercicio de matemáticas aburridas. Sirve para entender cosas gigantes:

• Los Glaciares: El hielo no es solo hielo duro; fluye lentamente como un líquido muy espeso. Este modelo ayuda a predecir cuánto tardará un glaciar en llegar al mar.
• El Interior de la Tierra: Bajo nuestros pies, la roca es sólida, pero a lo largo de millones de años, se comporta como un líquido muy viscoso, permitiendo que las placas tectónicas se muevan.
• Volcanes y Magma: El magma es una mezcla de roca fundida y cristales. Entender si se comporta más como "gelatina" o como "chicle" ayuda a predecir erupciones.
• Medicina e Industria: Desde cómo fluye la sangre (que a veces se espesa) hasta cómo se procesan los plásticos.

5. La Conclusión en una frase

El autor nos dice: "Podemos crear un modelo matemático perfecto y único para materiales que son a la vez líquidos y sólidos, pero debemos tener cuidado de no usar atajos matemáticos que parezcan fáciles pero que nos den resultados incorrectos. Usar las herramientas correctas nos permite predecir el futuro de los glaciares, los volcanes y los materiales del mañana."

En resumen, es como si el autor hubiera creado un traductor universal que nos permite hablar el idioma de la Tierra y entender cómo se mueven sus entrañas, mezclando la lógica de los sólidos y los líquidos en una sola ecuación maestra.

Resumen técnico

Aquí presento un resumen técnico detallado del artículo "A few notes about viscoplastic rheologies" (Algunas notas sobre reologías viscoplásticas) de Tomáš Roubíček, estructurado según los puntos solicitados.

1. Problema y Contexto

El artículo aborda la complejidad de modelar matemáticamente el comportamiento reológico de materiales que exhiben tanto comportamiento viscoso (dependiente de la tasa de deformación) como plástico (deformación permanente tras superar un umbral).

• El desafío: Existen numerosos modelos empíricos para fluidos no newtonianos y materiales geofísicos (rocas, hielo, magma). Sin embargo, a menudo carecen de una base rigurosa cuando se combinan elementos en serie o en paralelo.
• La necesidad: Se requiere un marco unificado para sintetizar un único potencial de disipación convexo ("viscoplástico") a partir de los potenciales de elementos individuales (viscosidad lineal, plasticidad perfecta, viscosidades no lineales).
• El objetivo: Utilizar herramientas rigurosas del análisis convexo para derivar modelos matemáticamente consistentes y compararlos con los modelos empíricos comunes (basados en medias armónicas) utilizados en geofísica e ingeniería.

2. Metodología

El autor emplea exclusivamente el análisis convexo para construir y analizar los modelos reológicos. Las herramientas matemáticas clave incluyen:

• Potenciales de Disipación ($\zeta_{vp}$): Se asume que la tensión disipativa $\sigma$ es la derivada (o subdiferencial) de un potencial convexo $\zeta_{vp}$ respecto a la tasa de deformación $\varepsilon$.
• Conjugación Convexa ($^*$): Se utiliza la transformada de Legendre-Fenchel para pasar del potencial de disipación al potencial dual, facilitando el manejo de elementos en serie.
• Convolución Infimal ($\square$): Esta es la operación central para combinar elementos en serie. Si dos elementos están en serie, el potencial total es la convolución infimal de sus potenciales individuales ($\zeta_{vp} = \zeta_1 \square \zeta_2$).
• Suma de Potenciales: Para elementos en paralelo, el potencial total es simplemente la suma de los potenciales individuales ($\zeta_{vp} = \zeta_1 + \zeta_2$).
• Comparación Analítica: Se contrastan las derivaciones rigurosas (basadas en la convolución infimal) con las aproximaciones empíricas que asumen una "media armónica" de las viscosidades, reemplazando las tasas de deformación parciales por la tasa total.

3. Contribuciones Clave y Resultados

A. Escenarios Básicos (Sección 2)

El autor analiza dos configuraciones fundamentales combinando plasticidad perfecta (potencial homogéneo de grado 1) y viscosidad lineal (potencial cuadrático):

1. Modelo en Paralelo (Fluido de Bingham):
   • La tensión es la suma de las tensiones de cada elemento.
   • Resultado: $\zeta_{vp} = \zeta_1 + \zeta_2$.
   • Comportamiento: La viscosidad efectiva es la suma de la viscosidad lineal y un término que decae con la tasa de deformación ($D + \sigma_a/|\varepsilon|$).
2. Modelo en Serie (Viscoplasticidad estricta):
   • La tasa de deformación es la suma de las tasas de cada elemento bajo una tensión común.
   • Resultado: $\zeta_{vp} = \zeta_1 \square \zeta_2$.
   • Hallazgo importante: Para la plasticidad perfecta y viscosidad cuadrática, este potencial resultante es la aproximación de Yosida de la función de plasticidad, conocida como función de Huber.
   • Ventaja: A diferencia del modelo de Bingham, este modelo en serie produce un potencial $\zeta_{vp}$ suave (diferenciable), eliminando la discontinuidad en la derivada del potencial conjugado, lo cual es crucial para la estabilidad numérica y la modelización de deslizamientos rápidos y fluencia lenta.

B. Escenarios de Tres Elementos (Sección 3)

Se exploran combinaciones de dos viscosidades lineales y un elemento plástico para suavizar aún más las transiciones:

• Se demuestra que dos variantes topológicas diferentes (una con viscosidad en serie y otra con viscosidad en paralelo al bloque plástico) son matemáticamente equivalentes si se ajustan adecuadamente los parámetros (ecuación 7).
• Estos modelos generan una viscosidad efectiva que es una suma de una media armónica y una viscosidad constante, permitiendo modelar fenómenos de "adelgazamiento por cizalla" (shear-thinning) de manera rigurosa.

C. Viscosidades No Lineales y Potencias (Sección 4)

El artículo extiende el análisis a fluidos con leyes de potencia (modelos Norton-Hoff):

• Combinación Serie (Plasticidad + Ley de Potencia): Esto conduce al modelo de fluido Herschel-Bulkley.
• Dificultad de Cálculo: El autor demuestra que obtener una forma explícita del potencial de disipación resultante ($\zeta_{vp} = \zeta_{plástico} \square \zeta_{no-lineal}$) es extremadamente difícil o imposible analíticamente para exponentes generales $n$.
  • Se resuelven casos específicos ($n=2, n=3$) utilizando fórmulas algebraicas (cuadrática y de Cardano para cúbicas).
• Crítica a los Modelos Empíricos: Se compara la derivación rigurosa (convolución infimal) con la aproximación empírica común (media armónica de viscosidades, ecuación 15).
  • Resultado: Los modelos empíricos no son equivalentes a los modelos rigurosos. La aproximación empírica falla en capturar la física correcta, especialmente en regímenes de baja y alta tasa de deformación, y produce curvas de viscosidad efectiva con formas distintas (ver Figura 6).

D. Generalización a Grandes Deformaciones (Sección 4, Rem. 6)

Se propone una formulación para grandes deformaciones utilizando la descomposición multiplicativa de Kröner-Lee-Liu.

• Contribución: Se muestra que trabajar con el potencial conjugado ($\zeta^*$) en lugar del potencial directo es ventajoso en sistemas viscoelastoplásticos en serie. En este enfoque, la combinación en serie se traduce simplemente en una suma de potenciales conjugados, evitando la compleja construcción de la convolución infimal en el espacio de deformaciones.

4. Significado e Impacto

• Rigor Matemático: El artículo establece una base teórica sólida para la reología viscoplástica, demostrando que muchos modelos empíricos usados en geofísica (manto terrestre, glaciología) son aproximaciones que pueden desviarse significativamente de la realidad física descrita por la termodinámica de procesos irreversibles.
• Suavidad Numérica: La identificación de la función de Huber como el potencial correcto para la combinación serie (plasticidad + viscosidad) ofrece una alternativa suave y diferenciable a los modelos de Bingham, lo cual es vital para la convergencia en simulaciones numéricas de elementos finitos.
• Aplicaciones Geofísicas: Proporciona herramientas para modelar con mayor precisión fenómenos como la fluencia de dislocaciones en rocas, el flujo de magma y el movimiento de placas tectónicas, donde la transición entre comportamiento plástico y viscoso es crítica.
• Advertencia sobre Modelos Empíricos: Advierte explícitamente contra el uso indiscriminado de la "media armónica" de viscosidades en regímenes no lineales, mostrando que esto altera la física del modelo (diferencia entre ecuaciones 14 y 15).

En resumen, el paper actúa como un puente entre la teoría abstracta del análisis convexo y la modelización práctica de materiales complejos, ofreciendo derivaciones exactas y advirtiendo sobre las limitaciones de las simplificaciones empíricas comunes.