Volumetric Growth in Linear Elasticity Driven by an Optimality Criterion

Este trabajo propone un marco novedoso para modelar el crecimiento volumétrico en elasticidad lineal formulando el tensor de crecimiento como la solución a un problema de optimización con restricciones impulsado por un funcional objetivo, en lugar de depender de leyes fenomenológicas prescritas, el cual se resuelve numéricamente mediante discretización por elementos finitos como un flujo de gradiente proyectado.

Lo esencial

Imagina un organismo vivo, como un árbol o un tumor, creciendo. Por lo general, los científicos intentan predecir este crecimiento escribiendo reglas específicas, como "las células crecen más rápido donde la presión es alta". Este artículo propone una forma diferente de pensar. En lugar de dar al crecimiento un conjunto de instrucciones, los autores sugieren que el crecimiento es como un optimizador inteligente que toma la mejor decisión posible en cada paso individual.

Aquí está la idea central desglosada con analogías simples:

1. El "Constructor Inteligente" frente al "Seguidor de Reglas"

Piensa en un cuerpo en crecimiento como una obra de construcción.

• La Vieja Forma (Seguidor de Reglas): Le das a los trabajadores un manual: "Si sientes un empuje, añade un ladrillo aquí. Si sientes un tirón, añade un ladrillo allá". El crecimiento está dictado por un guion fijo.
• La Forma de Este Artículo (Constructor Inteligente): No les das a los trabajadores un guion. En su lugar, les dices: "Tienes una cantidad limitada de ladrillos nuevos para añadir hoy. Tu objetivo es organizar estos ladrillos para que el edificio sea lo más estable posible (o lo más redondo posible) mientras obedeces las leyes de la física". Los trabajadores deciden dónde poner los ladrillos para lograr ese objetivo. El crecimiento no está "prescrito"; emerge del deseo de ser óptimo.

2. La "Lámina de Goma" y el "Estiramiento Oculto"

Los autores utilizan un modelo simplificado de física (elasticidad linealizada) para describir el cuerpo. Imagina una lámina de goma.

• Elasticidad: Si tiras de la lámina, se estira y quiere volver a su estado original.
• Crecimiento: Ahora, imagina que la lámina puede "crecer" secretamente nuevo material en lugares específicos. Esto es como si la lámina decidiera estirarse permanentemente en un área sin que tú la tires.
• El Conflicto: Si la lámina crece en un lugar pero no en otro, crea tensión interna (esfuerzo), como intentar encajar un clavo cuadrado en un agujero redondo.

3. El "Presupuesto Diario" y la Regla de "Sin Retroceso"

El crecimiento ocurre en pequeños pasos, como días en un calendario.

• El Presupuesto: En cada paso, el cuerpo tiene un "presupuesto" de nueva masa (volumen) para añadir. Puede ser un presupuesto global (todo el cuerpo se hace un poco más grande) o un presupuesto local (lugares específicos reciben más).
• La Regla de Sin Retroceso: El artículo impone una regla de que el cuerpo solo puede añadir material, nunca encogerse ni eliminarlo. Es como una calle de sentido único para el crecimiento; no puedes "descrecer".
• El Objetivo: El cuerpo debe decidir dónde gastar su presupuesto diario.
  • Escenario A (La Viga): Si el cuerpo es una viga que sostiene un peso pesado, el "Constructor Inteligente" añadirá material de una manera que haga la viga más rígida y reduzca el trabajo que el peso ejerce sobre ella. Es como si la viga "pensara": "Necesito engrosarme aquí para dejar de doblarme tanto".
  • Escenario B (La Masa Libre): Si el cuerpo está flotando sin peso, el objetivo podría ser convertirse en un círculo (porque un círculo tiene el perímetro más corto para un área dada). El cuerpo "piensa": "Necesito reorganizar mi crecimiento para volverse más redondo".

4. El "GPS Matemático"

¿Cómo sabe el cuerpo dónde crecer? Los autores muestran que este proceso es matemáticamente equivalente a un flujo de gradiente proyectado.

• Imagina que estás en un paisaje montañoso (que representa la energía o la "mala calidad" de la forma actual). Quieres llegar al valle más bajo (la mejor forma).
• Das un paso cuesta abajo.
• El Giro: Tienes una "restricción de presupuesto" (solo puedes añadir cierta cantidad de masa) y una "regla de sentido único" (no puedes encogerte).
• Las matemáticas actúan como un GPS que te dice: "Da un paso cuesta abajo, pero si ese paso viola tu presupuesto o intenta encogerte, deslízate a lo largo del borde del área permitida hasta encontrar el mejor paso válido".

5. Lo que Mostraron los Experimentos Computacionales

Los autores ejecutaron simulaciones en una computadora para ver qué sucede cuando dejas que este "Constructor Inteligente" tome el control:

• Vigas bajo carga: Las vigas crecieron naturalmente más gruesas en el medio y se afilaron en los extremos, remodelándose efectivamente para volverse más fuertes y rígidas contra la carga. No solo se hicieron más grandes; se volvieron más inteligentes sobre dónde hacerse más grandes.
• Cuerpos flotantes: Evolucionaron naturalmente hacia formas circulares, que es la forma más eficiente para minimizar el perímetro.
• Patrones de Esfuerzo: En el caso de los cuerpos flotantes, el crecimiento creó un patrón de esfuerzo específico: el centro se comprimió (apretado) y el borde exterior se estiró (tensión). Los autores notaron que este patrón específico (núcleo comprimido, cáscara tensa) se observa realmente en tumores biológicos reales, lo que sugiere que esta lógica de "optimización" podría ser un principio fundamental de cómo se forman las formas biológicas.

Resumen

El artículo argumenta que no necesitas reglas biológicas complejas para explicar por qué las cosas crecen en formas específicas. Si simplemente le dices a un objeto en crecimiento: "Aquí tienes tu nuevo material. Úsalo para hacer que el sistema sea lo más eficiente posible (ya sea mecánica o geométricamente) sin encogerte", el objeto evolucionará naturalmente hacia formas complejas, estables y a menudo con apariencia biológica. El crecimiento es el resultado de un problema de optimización, no de un guion preescrito.

Resumen técnico

Resumen Técnico: Crecimiento Volumétrico Impulsado por Optimización en un Marco de Elasticidad Linealizada

Formulación del Problema
El artículo aborda la modelización del crecimiento volumétrico en cuerpos elásticos, un proceso en el que se genera masa a lo largo del volumen, alterando la geometría y la respuesta mecánica. Los enfoques tradicionales de la mecánica de medios continuos suelen prescribir la evolución de un tensor de crecimiento ($E_g$) mediante leyes cinéticas fenomenológicas dependientes del esfuerzo, la deformación o estímulos bioquímicos. En contraste, este trabajo propone un marco en el que el tensor de crecimiento no se prescribe, sino que se determina implícitamente como la solución a un problema de optimización con restricciones en cada paso de tiempo discreto.

El escenario es la elasticidad linealizada en dos dimensiones (con una extensión directa a tres dimensiones). El cuerpo ocupa un dominio $\Omega$ con un borde descompuesto en partes de Dirichlet ($\partial_D \Omega$) y Neumann ($\partial_N \Omega$). El crecimiento se modela mediante una descomposición aditiva del tensor de deformación, donde la deformación total es la suma de la deformación elástica y una deformación de crecimiento inelástica ($E_g$). Las ecuaciones de equilibrio se modifican de modo que $E_g$ actúa como un término fuente interno.

El problema central es una minimización incremental sobre el campo de desplazamientos $u^{(i)}$ y el campo tensorial de crecimiento $E_g^{(i)}$ en cada paso de tiempo $i$. La optimización está sujeta a tres restricciones principales:

1. Equilibrio: Las ecuaciones de equilibrio de la elasticidad linealizada deben satisfacerse.
2. Balance de Masa: El cambio de volumen total (o local) debido al crecimiento debe coincidir con un incremento prescrito $\Gamma^{(i)}$.
3. Irreversibilidad (Acreción): El crecimiento debe ser puramente acretivo, lo que significa que el incremento del tensor de crecimiento debe ser semidefinido positivo ($\lambda_{\min}(E_g^{(i)} - E_g^{(i-1)}) \geq 0$).

El funcional objetivo $\Phi$ codifica el mecanismo impulsor. El artículo investiga dos criterios específicos:

• Minimización del Trabajo Externo: Representa una adaptación impulsada mecánicamente donde el cuerpo busca aumentar la rigidez.
• Minimización del Perímetro: Representa una evolución impulsada geométricamente (principio isoperimétrico).

Se incluye un término de regularización cuadrática, proporcional a $|E_g^{(i)} - E_g^{(i-1)}|^2$, para penalizar grandes desviaciones de la configuración anterior. Este término, inspirado en la teoría de movimientos minimizantes de De Giorgi, asegura la continuidad temporal y conduce a una estructura de flujo de gradiente en el límite continuo.

Metodología
Los autores emplean un Método de Elementos Finitos (MEF) para discretizar el problema continuo.

• Discretización: El campo de desplazamientos se aproxima utilizando funciones base afines por tramos y continuas, mientras que el tensor de crecimiento se aproxima utilizando funciones base constantes por tramos en elementos triangulares.
• Reducción: Al resolver explícitamente las ecuaciones de equilibrio lineales para el campo de desplazamientos en términos del tensor de crecimiento, los autores eliminan las variables de desplazamiento. Esto reduce el problema a un problema de minimización con restricciones de dimensión finita que involucra únicamente las variables de crecimiento.
• Optimización: El problema discreto resultante se resuelve numéricamente utilizando la rutina fmincon de MATLAB con el algoritmo de Programación Cuadrática Secuencial (SQP). Se proporcionan gradientes analíticos para la función objetivo y las restricciones para garantizar la convergencia.
• Aproximación Analítica: Para revelar la estructura subyacente, los autores derivan una solución analítica aproximada omitiendo temporalmente la restricción no lineal de autovalores (irreversibilidad). Esto permite la derivación de una expresión cerrada para el incremento de crecimiento, demostrando que la evolución sigue un flujo de gradiente proyectado en el espacio de las deformaciones inelásticas.

Resultados Clave
Se realizaron experimentos numéricos en tres casos representativos: una viga doblemente empotrada, una viga en voladizo (ambas bajo carga externa) y un cuerpo libre (minimización del perímetro).

• Adaptación de la Viga: Bajo cargas externas, las vigas evolucionan para aumentar la rigidez estructural. La distribución del crecimiento es no uniforme, concentrándose en regiones donde las fibras de material se alargan. La viga doblemente empotrada desarrolla predominantemente tensiones residuales compresivas, mientras que la viga en voladizo acomoda el crecimiento mediante alargamiento, reduciendo las magnitudes de tensión.
• Minimización del Perímetro: El cuerpo libre evoluciona hacia una forma circular, consistente con el principio isoperimétrico clásico.
• Patrones de Tensión: En el caso impulsado por el perímetro, el campo de tensión residual exhibe un núcleo compresivo y una envoltura exterior tensa. Los autores señalan que este patrón refleja observaciones experimentales y numéricas en tumores murinos y humanos, donde el núcleo del tumor está comprimido y la envoltura exterior está bajo tensión.
• Analítico vs. Numérico: Las soluciones analíticas aproximadas (derivadas sin la restricción de irreversibilidad) coinciden estrechamente con los resultados numéricos completos para los casos de viga, con diferencias neglijibles en los valores del funcional objetivo y los minimizadores. Para el caso del perímetro, la solución analítica converge al resultado numérico a pesar de las violaciones iniciales de la restricción de autovalores, siempre que el parámetro de regularización sea suficientemente grande.

Significado y Afirmaciones
El artículo afirma demostrar que el crecimiento volumétrico puede formularse como un proceso impulsado por optimización donde el tensor de crecimiento es seleccionado por un principio de selección mecánica en lugar de una ley fenomenológica.

• Cambio Teórico: El trabajo se distingue al tratar el crecimiento no como una entrada cinemática descriptiva, sino como el resultado de la minimización de un funcional específico sujeto a restricciones mecánicas.
• Insight Estructural: El modelo revela que la evolución puede interpretarse como un flujo de gradiente proyectado en el espacio de las deformaciones inelásticas. Los autores destacan que la convexidad del funcional de complianza es crucial; cuando se pierde la convexidad, pueden surgir no unicidad y localización.
• Prueba de Concepto: La teoría linealizada sirve como prueba de concepto, mostrando que la optimización con restricciones mecánicas por sí sola es suficiente para generar morfologías complejas inducidas por el crecimiento (como la adaptación de la forma y los patrones de tensión residual) sin asumir leyes de crecimiento biológico específicas.
• Perspectiva Futura: Los autores declaran que este marco allana el camino para extender el enfoque a la elasticidad totalmente no lineal y a sistemas más realistas que involucran objetivos energéticos competitivos, aunque no afirman haber resuelto estos sistemas biológicos complejos en el trabajo actual.

El artículo concluye que el marco propuesto unifica con éxito el equilibrio mecánico, el balance de masa y la irreversibilidad dentro de un marco variacional, proporcionando una herramienta robusta para predecir morfologías impulsadas por el crecimiento.