Topology optimization of passively moving rigid bodies in unsteady flows

Este artículo propone un método de optimización topológica para el diseño de cuerpos rígidos de movimiento pasivo en flujos fluidos no estacionarios mediante el acoplamiento de la dinámica de cuerpos rígidos con ecuaciones de fluidos cinéticos de red y la aplicación de un análisis de sensibilidad basado en el adjunto para optimizar las formas para movimientos traslacionales y rotacionales.

Lo esencial

Imagina que estás intentando diseñar la vela perfecta para un barco o la pala perfecta para una turbina eólica. Normalmente, los ingenieros comienzan con una forma, la ponen al viento, ven cómo se mueve, ajustan la forma y vuelven a intentarlo. Este artículo presenta una forma superinteligente y automatizada de hacer eso, pero con un giro: en lugar de que la forma simplemente se quede ahí mientras el viento sopla sobre ella, se le permite moverse porque el viento la está empujando.

Aquí hay un desgero de cómo funciona esta "magia", utilizando analogías sencillas:

1. La gran idea: El "Bailarín Pasivo"

La mayoría de los programas informáticos que diseñan formas asumen que el objeto está pegado al suelo (como un puente o una tubería estacionaria). Si quieres diseñar una pieza móvil, como la pala de un ventilador, normalmente le dices a la computadora: "Gira esta pala a 100 rotaciones por minuto", y la computadora calcula el flujo de aire.

Este artículo cambia el guion. Trata al objeto como un bailarín en una pista de baile.

• La forma antigua (Activa): Le dices al bailarín exactamente cómo moverse y observas cómo se mueve el aire a su alrededor.
• La nueva forma (Pasiva): No le dices al bailarín cómo moverse. Simplemente pones la música (el viento) y le pides a la computadora que diseñe el cuerpo del bailarín para que la música lo empuje naturalmente a girar o deslizarse lo más lejos posible. El movimiento del bailarín es un resultado del viento, no un comando.

2. El truco de las dos mallas: El Mapa vs. El Terreno

Para que esto funcione, la computadora utiliza un truco ingenioso llamado "mallas separadas". Imagina que estás dibujando el mapa de una isla en movimiento sobre una hoja de papel milimetrado.

• La Malla de Diseño (El Mapa): Aquí es donde se dibuja la forma. Es como un cuaderno de bocetos. La computadora decide dónde está el material "sólido" (la isla) y el espacio "vacío" (el agua) aquí.
• La Malla de Análisis (El Terreno): Aquí es donde ocurre la física. Es una malla fija de agua y viento.

Cada diminuta fracción de segundo, el "Mapa" (la forma) se mueve físicamente y rota. Luego, la computadora proyecta ese mapa en movimiento sobre la malla del "Terreno" fija para calcular cómo el viento empuja contra él. Después de que el viento empuja, la computadora calcula cómo debería moverse el objeto a continuación, actualiza el Mapa y repite el ciclo. Es como tomar una foto de un objeto en movimiento, calcular la fuerza del viento, mover el objeto y tomar la siguiente foto instantáneamente.

3. La fuerza "Fantasma" (Fuerza de Brinkman)

¿Cómo sabe la computadora dónde está el objeto sólido? Utiliza un concepto llamado fuerza de Brinkman.
Piensa en el área de diseño como una habitación llena de miel invisible y pegajosa.

• Donde hay material sólido, la miel es súper espesa y pegajosa. El viento no puede pasar a través de ella; simplemente empuja contra la superficie.
• Donde hay espacio vacío, la miel es delgada o inexistente, por lo que el viento fluye libremente.
  La computadora no necesita dibujar una línea dura; simplemente ajusta la "pegajosidad" de la miel en cada punto. Si la pegajosidad es alta, es una pared; si es baja, es aire. Esto permite que la forma se transforme suavemente de una forma a otra.

4. Las matemáticas del "Viaje en el Tiempo" (Método Adjunto)

Para encontrar la forma perfecta, la computadora tiene que saber: "Si cambio este diminuto punto de material aquí, ¿cuánto mejor se moverá el objeto?".
Calcular esto para cada uno de los puntos tardaría una eternidad. Por eso, los autores utilizan un método llamado Método de la Variable Adjunta.

• La Analogía: Imagina que estás intentando encontrar el mejor camino hacia la cima de una montaña en la oscuridad. En lugar de caminar por todos los caminos posibles para ver cuál es el mejor, alumbras con una linterna hacia atrás desde la cima. La luz te muestra exactamente qué pasos te llevan hacia arriba de la manera más eficiente.
• En este artículo, la "linterna" corre hacia atrás en el tiempo, calculando cómo las fuerzas del viento y el movimiento del objeto habrían reaccionado a cada pequeño cambio en la forma. Esto le da a la computadora un "mapa de sensibilidad" que le indica exactamente dónde añadir o quitar material para obtener el mejor resultado.

5. Los Resultados: ¿Qué construyeron?

El equipo probó esto en tres escenarios:

1. La Vela 2D: Diseñaron una forma que comienza quieta y es empujada por el viento para deslizarse horizontalmente. El resultado se pareció a un ala de avión curva (perfil alar). El viento empujaba más fuerte en la parte superior que en la inferior, creando una sustentación que tiraba del objeto hacia adelante.
2. La Turbina 2D: Diseñaron una forma que gira. El resultado se pareció a una hélice de cuatro palas. El viento golpeaba las palas curvas, creando un giro que la hacía rotar.
3. La Turbina 3D: Hicieron lo mismo en 3D. El resultado se pareció a una turbina eólica del mundo real.

6. El problema de la "Escala de Grises"

En estos diseños por computadora, los bordes de las formas no siempre son líneas perfectamente nítidas de blanco y negro. A veces son "escala de grises": un poco de ambos, sólido y aire.

• En los ejemplos 2D, los autores descubrieron que incluso si hacían la forma perfectamente nítida (blanco y negro), el rendimiento era casi el mismo. Los bordes "difusos" no afectaron el resultado.
• En el ejemplo 3D, los bordes "difusos" importaron más. Debido a que la malla de la computadora era un poco "tosca" (baja resolución), los bordes difusos cambiaron la forma en que el viento golpeaba las palas. Esto sugiere que para formas 3D complejas, necesitamos un "mapa" más fino para obtener un resultado perfecto.

Resumen

Este artículo presenta una nueva forma de diseñar máquinas móviles (como velas o turbinas) donde la computadora determina la forma y el movimiento simultáneamente. Trata al objeto como un bailarín pasivo empujado por el viento, utiliza un truco de "miel pegajosa" para definir la forma, y ejecuta una simulación matemática hacia atrás en el tiempo para encontrar la forma más eficiente. El resultado son formas que naturalmente parecen alas y hélices, optimizadas para moverse lo más lejos posible o girar lo más rápido posible bajo las fuerzas de los fluidos.

Resumen técnico

Resumen Técnico: Optimización de Topología de Cuerpos Rígidos de Movimiento Pasivo en Flujos No Estacionarios

Planteamiento del Problema
Este estudio aborda la optimización de la topología de cuerpos rígidos que son impulsados por fuerzas inducidas por el flujo de fluido, tales como velas y turbinas. A diferencia de la optimación de topología tradicional en dinámica de fluidos, donde el movimiento del cuerpo rígido es prescrito (activo) para impulsar el flujo, o donde se realiza un seguimiento de partículas en un acoplamiento de una vía donde estas no afectan al fluido, este trabajo se centra en un sistema de acoplamiento bidireccional. Aquí, el movimiento del cuerpo rígido es pasivo, lo que significa que su trayectoria y rotación están determinadas enteramente por las fuerzas del fluido que actúan sobre él. El desafío rector implica resolver simultáneamente las ecuaciones de flujo de fluido no estacionario con las ecuaciones de movimiento del cuerpo rígido, donde las fuerzas del fluido sirven como entradas para las ecuaciones de movimiento, y la posición del cuerpo sirve como condición de contorno para el fluido.

Metodología
El marco propuesto integra varias técnicas numéricas para manejar la compleja interacción entre el sólido en movimiento y el fluido no estacionario:

• Ecuaciones Rectores: El flujo de fluido se modela utilizando las ecuaciones de Navier-Stokes incomprimibles, resueltas mediante el Esquema Cinético de Red (LKS), una versión extendida del Método de Lattice Boltzmann (LBM) adecuada para flujos no estacionarios. La dinámica del cuerpo rígido está gobernada por las ecuaciones de Newton-Euler, donde las fuerzas y torques se computan como la reacción de la fuerza de Brinkman dentro del dominio del fluido.
• Grillas de Diseño y Análisis Separadas: Un componente central de la metodología es la separación de la grilla de diseño (donde se define la forma) de la grilla de análisis (donde se computan las variables de estado del fluido). La grilla de diseño experimenta un movimiento de cuerpo rígido en cada paso de tiempo y se mapea sobre la grilla de análisis mediante una función de peso. Esto permite la representación directa del movimiento del cuerpo rígido sin necesidad de remallar el dominio del fluido.
• Análisis de Sensibilidad Adjunto: Para computar las sensibilidades de diseño de manera eficiente, el estudio emplea el Esquema Cinético de Red Adjunto (ALKS). Los autores derivan ecuaciones adjuntas continuas que acoplan las ecuaciones de flujo de fluido con las ecuaciones de movimiento del cuerpo rígido. Esta formulación trata las variables de movimiento del cuerpo rígido (posición, velocidad, ángulo de rotación y velocidad angular) como campos de estado adicionales dentro del sistema adjunto.
• Algoritmo de Optimización: La forma se actualiza utilizando el Método de las Asíntotas Móviles (MMA). Para asegurar diseños de tipo binario y fabricables, se utiliza un filtro de densidad con proyección de Heaviside. Se aplica un esquema de continuación al parámetro de convexidad y al parámetro de estiramiento de la proyección para guiar la optimización desde un estado inicial difuso hacia una distribución sólido-fluido clara.
• Aproximación Periódica: Para problemas que involucran movimiento cíclico (por ejemplo, turbinas), se emplea una aproximación de problema periódico. Los valores iniciales para los campos de estado y adjunto al inicio de una iteración de optimización se establecen con los valores terminales de la iteración anterior, reduciendo el costo computacional al evitar la simulación de fases de arranque transitorias en cada iteración.

Contribuciones Clave

1. Formulación de Movimiento Pasivo: El artículo presenta un novedoso marco de optimización de topología específicamente para cuerpos rígidos impulsados por fuerzas de fluido, distinguiéndolo del movimiento activo o del seguimiento de partículas de una vía.
2. Derivación de Adjunto Acoplado: Proporciona una derivación rigurosa de las ecuaciones adjuntas y las sensibilidades de diseño que acoplan completamente la dinámica de fluidos con las ecuaciones de movimiento del cuerpo rígido, tratando las variables cinemáticas del cuerpo como campos de estado.
3. Implementación de Separación de Grillas: El estudio demuestra la efectividad del enfoque de separación de grillas para problemas de cuerpos rígidos en movimiento y no estacionarios, facilitando la transferencia de cantidades como el coeficiente de Brinkman y la velocidad del sólido entre grillas móviles y estacionarias.
4. Validación Numérica: El método se verifica contra aproximaciones de diferencias finitas tanto para casos de traslación como de rotación, confirmando la precisión de las sensibilidades derivadas.

Resultos
El método se aplicó a tres ejemplos numéricos:

• Vela 2D: La optimización maximizó el desplazamiento horizontal de un cuerpo rígido. La forma resultante se asemeja a un perfil alar, generando fuerzas de sustentación que propulsan el cuerpo hacia adelante. La forma optimizada superó significativamente a una forma de referencia con una superficie superior plana.
• Turbina 2D: El objetivo fue maximizar el ángulo de rotación. Al emplear una condición periódica cíclica, la optimización produjo una forma con cuatro álabes similares a perfiles alares que rotan efectivamente bajo las fuerzas del fluido. Los resultados indicaron que la presencia de regiones de escala de grises en la forma optimizada tuvo un efecto insignificante en el rendimiento al ser binarizada.
• Turbina 3D: Se optimizó una forma rotatoria tridimensional, resultando en una estructura tipo turbina. Sin embargo, debido a limitaciones computacionales resultantes de una resolución de malla más gruesa, la versión binarizada mostró diferencias no despreciables en fase y periodo en comparación con la forma optimizada, resaltando la sensibilidad de las estructuras delgadas 3D a la resolución de la malla.

Significancia y Reclamaciones
Los autores afirman que el método propuesto optimiza eficazmente las formas de dispositivos como velas y turbinas donde el rendimiento está dictado por la interacción entre las fuerzas del fluido y el movimiento del cuerpo rígido. El estudio demuestra que un análisis no estacionario es necesario para tales problemas, ya que el ángulo de ataque efectivo y las fuerzas cambian dinámicamente con el movimiento del cuerpo.

El artículo señala modestamente que, si bien el método es efectivo, existen limitaciones, particularmente en casos tridimensionales donde las restricciones de resolución de malla pueden conducir a discrepancias entre las formas optimizadas y las binarizadas. Los autores sugieren que el trabajo futuro debería centrarse en mejorar la resolución de la malla (por ejemplo, mediante métodos adaptativos) y extender el marco a flujos turbulentos. El estudio concluye que este enfoque de optimización de topología acoplado proporciona una herramienta robusta para diseñar componentes móviles impulsados por fluidos, ofreciendo perspectivas sobre los mecanismos físicos de generación de sustentación y arrastre en regímenes no estacionarios.