Application of Metric-Based Mesh Adaptation to Hypersonic Aerothermal Simulations Using US3D

Este artículo demuestra que la adaptación de malla basada en métricas aplicada a simulaciones aerotérmicas hipersónicas con el código US3D permite predecir con precisión el calentamiento superficial en geometrías complejas, como cápsulas de entrada atmosférica con sistemas de control de reacción, superando las limitaciones de los enfoques de malla estructurada tradicional.

Lo esencial

¡Claro que sí! Imagina que estás diseñando un vehículo espacial que va a viajar a Marte. Cuando este vehículo entra en la atmósfera a velocidades increíbles (hipersónicas), el aire se calienta tanto que podría derretirlo. Para evitar que se funda, los ingenieros necesitan simular exactamente cómo se comporta ese calor en la computadora antes de construirlo.

Este documento de la NASA explica cómo han mejorado una herramienta de simulación para hacer estos cálculos de manera más inteligente y precisa. Aquí te lo explico con analogías sencillas:

1. El Problema: Dibujar un mapa con una cuadrícula rígida

Imagina que quieres dibujar un mapa muy detallado de una montaña con muchos picos y valles.

• El método antiguo (Mallas estructuradas): Era como intentar dibujar esa montaña usando solo papel milimetrado rígido. Si la montaña tiene una forma extraña o un detalle pequeño (como un árbol o una roca), el papel cuadrado no encaja bien. Tienes que deformar el papel o dejar huecos, lo que hace que el dibujo sea borroso y poco preciso. Además, si quieres añadir un detalle nuevo, tienes que redibujar todo el mapa desde cero.
• El problema específico: En el espacio, el "aire" se convierte en un plasma caliente que choca contra la nave. Si la cuadrícula de tu mapa no está perfectamente alineada con ese choque, los cálculos de calor salen mal, como si intentaras medir la lluvia con un cubo que tiene agujeros.

2. La Solución: Un "Mapa Inteligente" que se adapta (MIMIC)

Los autores de este documento (Dirk Ekelschot y su equipo) han desarrollado una herramienta llamada MIMIC.

• La analogía: Imagina que en lugar de usar papel milimetrado fijo, usas una masa de plastilina inteligente.
  • Al principio, la plastilina es gruesa y tosca.
  • La computadora calcula dónde está el "calor" (como si la plastilina se estuviera quemando en ciertos puntos).
  • Entonces, la herramienta estira y aplana la plastilina automáticamente. Donde hay mucho calor o turbulencia, la plastilina se hace muy fina y detallada (como un zoom). Donde no pasa nada importante, la plastilina se queda gruesa para ahorrar espacio.
  • Lo mejor es que esta plastilina puede cambiar de forma infinitamente para encajar en cualquier objeto, por extraño que sea.

3. Dos Experimentos Clave

Experimento A: La Media Esfera (El entrenamiento)
Primero, probaron con una forma simple: una media esfera.

• Descubrieron que si usaban bloques de forma rectangular (hexaedros) en la capa de aire pegada a la superficie, los resultados eran mucho más suaves y precisos que si usaban bloques triangulares (prismas).
• Analogía: Es como pavimentar un camino. Si usas ladrillos rectangulares bien alineados, el camino es liso. Si usas piedras triangulares desordenadas, el camino es irregular y el coche (el calor) rechina.

Experimento B: La Cápsula de Marte (La prueba real)
Luego, probaron con una cápsula real de 70 grados que iría a Marte. Esta cápsula tiene una característica complicada: tiene 8 pequeños motores de control (RCS) en la parte trasera, como pequeñas boquillas.

• El desafío: Con el método antiguo (papel milimetrado), era casi imposible modelar esos 8 motores sin simplificar la geometría o hacer un mapa tan complejo que la computadora se volviera loca.
• El éxito: Con la "plastilina inteligente" (MIMIC), pudieron incluir los 8 motores tal cual son. La herramienta ajustó el mapa automáticamente alrededor de esos motores y en la estela turbulenta detrás de la cápsula.
• Resultado: Los cálculos de calor fueron tan precisos como los métodos antiguos, pero sin tener que redibujar todo el mapa manualmente. Además, pudieron ver detalles del flujo de aire detrás de la cápsula que antes eran invisibles.

4. ¿Por qué es importante esto?

1. Ahorro de tiempo y esfuerzo: Antes, los ingenieros pasaban semanas creando mallas perfectas para formas simples. Ahora, la computadora hace el trabajo sucio de ajustar la malla automáticamente.
2. Geometrías complejas: Permite simular naves con detalles reales (como motores, antenas o formas irregulares) que antes eran "demasiado difíciles" de simular con precisión.
3. Seguridad: Al tener mapas de calor más precisos, podemos diseñar escudos térmicos más seguros para los astronautas y las naves, asegurando que no se quemen al entrar en la atmósfera.

En resumen

Este documento nos dice que la NASA ha aprendido a usar una herramienta de "auto-ajuste" para sus simulaciones de vuelo espacial. En lugar de forzar la geometría a encajar en una cuadrícula rígida, ahora dejan que la cuadrícula se adapte inteligentemente a la forma del objeto y al comportamiento del calor. Es como pasar de dibujar con regla y compás a usar un pincel mágico que sabe exactamente dónde poner cada pincelada para capturar la realidad con perfección.

Resumen técnico

Resumen Técnico: Aplicación de Adaptación de Malla Basada en Métricas a Simulaciones Aerotérmicas Hipersónicas Usando US3D

1. Problema y Contexto
Las simulaciones de flujo hipersónico para vehículos de entrada atmosférica presentan desafíos significativos en la predicción precisa del calentamiento superficial (flujo de calor). Los códigos de Dinámica de Fluidos Computacional (CFD) de última generación, como DPLR y US3D, suelen depender de mallas estructuradas de bloques compuestas por elementos hexaédricos alineados con la onda de choque para garantizar la precisión. Sin embargo, la generación de estas mallas es laboriosa y poco flexible, especialmente para geometrías complejas que incluyen características como sistemas de control de reacción (RCS), pasos hacia atrás o estructuras de múltiples escalas.

Por otro lado, las mallas no estructuradas (tetraédricas) son más fáciles de generar pero tienden a producir predicciones de calentamiento inexactas debido a la generación de entropía numérica en la onda de choque, que se convecta hacia la superficie y contamina los resultados. Además, las mallas estructuradas tradicionales tienen dificultades para capturar dinámicas de estela complejas y fenómenos no estacionarios en la parte trasera (backshell) del vehículo.

2. Metodología
El trabajo presenta un flujo de trabajo que integra el solucionador de flujo compresible no estructurado US3D con la capacidad de mejora de mallas informada por métricas MIMIC (Metric-Informed Mesh Improvement Capability).

• Estrategia de Adaptación: Se utiliza una adaptación de malla anisotrópica basada en métricas. El tensor métrico se calcula utilizando el Hessiano (segundas derivadas) del campo de temperatura de la solución del flujo. Este tensor actúa como un indicador de error local que dicta dónde debe refinarse o alisarse la malla.
• Tipos de Elementos: La metodología soporta mallas híbridas. La capa límite cerca de la pared se mantiene fija y puede consistir en prismas o hexaedros (extruidos desde triángulos o cuadriláteros en la superficie), mientras que el volumen restante se tesseliza con tetraedros no estructurados que se adaptan dinámicamente.
• Proceso Iterativo:
  1. Se genera una malla base híbrida gruesa.
  2. Se resuelve el flujo (inicialmente con esquemas de primer orden para evitar efectos de "carbuncle" en la onda de choque).
  3. Se calcula el tensor métrico basado en la solución de temperatura.
  4. Se utiliza la biblioteca paralela ParMMG para adaptar los tetraedros del volumen, alineándolos con la métrica.
  5. El proceso se repite iterativamente (5 iteraciones en los casos de estudio) hasta converger en una solución precisa.
• Casos de Estudio:
  • Caso 1: Flujo supersónico sobre un hemisferio (gas perfecto) para comparar la influencia de usar capas límite de prismas vs. hexaedros.
  • Caso 2: Flujo hipersónico (mezcla reactiva de gases simulando la atmósfera marciana: $CO_2$-$N_2$) sobre una cápsula de entrada de 70° con 8 inyectores de RCS en la parte trasera.

3. Contribuciones Clave

• Flexibilidad Geométrica: Demostración de que las mallas no estructuradas adaptadas pueden manejar geometrías complejas (como los inyectores de RCS en la parte trasera) que son difíciles o imposibles de modelar con mallas estructuradas de bloques sin simplificaciones excesivas.
• Mejora de la Capa Límite: Extensión de la capacidad de MIMIC para soportar superficies cuadriláteras extruidas en capas de hexaedros, no solo triángulos extruidos en prismas. Se demuestra que los hexaedros en la capa límite producen distribuciones de flujo de calor más suaves y simétricas.
• Precisión con Mallas No Estructuradas: Validación de que, mediante la adaptación basada en métricas, las mallas no estructuradas pueden alcanzar niveles de precisión en el flujo de calor superficial comparables a los de las mallas estructuradas tradicionales (DPLR), siempre que la onda de choque esté adecuadamente resuelta.
• Automatización: Simplificación del proceso de generación de mallas, permitiendo la automatización de la refinación en regiones de interés (ondas de choque, capas de corte, estelas) sin intervención manual intensiva.

4. Resultados

• Caso Hemisferio:
  • Las mallas iniciales gruesas mostraron distribuciones de flujo de calor ruidosas y distorsionadas.
  • Tras 5 iteraciones de adaptación, la definición de la onda de choque mejoró drásticamente.
  • Se observó una alta sensibilidad al espesor de la capa límite: las mallas con 110 capas de hexaedros produjeron resultados más suaves y simétricos que las de 80 capas o las de prismas.
  • Los perfiles de presión y flujo de calor de las mallas adaptadas con hexaedros coincidieron dentro de los márgenes de error (±7.5%) con los datos de referencia de DPLR.
• Caso Cápsula de Entrada (70°):
  • La adaptación de la malla permitió capturar características de flujo complejas en la estela, como la interacción entre la capa de corte y la burbuja de recirculación, que serían imprácticas de resolver con mallas estructuradas de un solo bloque.
  • La inclusión de la geometría real de los inyectores de RCS en la malla no estructurada fue exitosa, algo que el enfoque estructurado no pudo hacer sin simplificar la geometría a una superficie lisa.
  • El flujo de calor superficial en el morro (forebody) convergió hacia una distribución simétrica y libre de "rayas" (streaks) numéricas tras las iteraciones de adaptación.
  • La comparación cualitativa con DPLR mostró un acuerdo razonable en el calentamiento del morro, aunque se notaron diferencias en la distancia de separación de la onda de choque en el hombro superior, atribuidas a diferencias en la resolución de la malla y el dominio de simulación (medio dominio vs. dominio completo).

5. Significado y Conclusiones
Este trabajo establece un flujo de trabajo viable para el análisis aerotérmico de diseños de entrada atmosférica complejos utilizando mallas no estructuradas adaptadas. Las conclusiones principales son:

1. Viabilidad: La adaptación de malla basada en métricas es esencial para simular flujos hipersónicos sobre geometrías complejas, superando las limitaciones de las mallas estructuradas tradicionales.
2. Precisión: Es posible obtener predicciones de flujo de calor de alta fidelidad con mallas no estructuradas, siempre que se utilicen capas límite de hexaedros y se realice una adaptación suficiente para alinear la malla con la onda de choque.
3. Eficiencia: El proceso reduce significativamente la carga de trabajo de generación de mallas, permitiendo a los ingenieros simular características geométricas detalladas (como propulsores RCS) sin sacrificar la precisión aerotérmica.
4. Futuro: El enfoque abre la puerta a estudios futuros que involucren interacciones de estela de chorro de RCS (plume-wake interaction) y el uso de esquemas numéricos avanzados (como HyperSolve) que funcionan mejor en mallas anisotrópicas.

En resumen, el documento demuestra que la combinación de US3D y MIMIC ofrece una alternativa robusta, flexible y precisa a los métodos estructurados tradicionales para la ingeniería aerotérmica hipersónica.