Validación de la malla CFD: los criterios que marcan la diferencia

9 min de lectura Fluidodinámica CFD

Si están aquí, probablemente ya hayan pasado horas lidiando con superficies CAD complejas o intentando que su solucionador procese una geometría complicada.

Hoy hablamos de un tema en el que, durante mi trayectoria como analista, he visto cometer errores con demasiada frecuencia: la validación de la malla. A menudo se tiende a pensar que, una vez pulsado Generate Mesh y tras ver aparecer una densa red de celdas, el trabajo ya está terminado. Como mucho, se llama a un compañero con más experiencia, que revisa la malla y se basa en una inspección visual para evaluar su validez. En realidad, es precisamente en ese punto donde empiezan los problemas (o los buenos resultados). Una malla visualmente atractiva no es necesariamente una malla capaz de producir resultados correctos.

En este artículo exploraremos las métricas fundamentales para determinar si vuestra discretización es una obra maestra de la ingeniería o una bomba de relojería numérica, prestando especial atención a las herramientas que Simcenter STAR-CCM+ pone a nuestra disposición.

¿Por qué validar la malla?

Antes de sumergirnos en los números, demos un paso atrás. La Dinámica de Fluidos Computacional no resuelve las ecuaciones de Navier-Stokes de forma continua en el espacio (simplemente, no es posible), sino que necesita discretizarlas, es decir, transformar un dominio infinito en un número finito de pequeños volúmenes de control, en cada uno de los cuales se resuelven las ecuaciones de Navier-Stokes.

¿El problema? Las ecuaciones son extremadamente exigentes. Si una celda está demasiado deformada, demasiado distorsionada o presenta una diferencia excesiva respecto a sus vecinas, el error numérico aumenta, la convergencia se ralentiza y, en el peor de los casos, el residuo de presión decide despegar hacia el infinito (y más allá).

Validar la malla significa asegurarse de que la cuadrícula que sustenta nuestra física no introduzca errores artificiales capaces de invalidar la solución. No se trata únicamente de una cuestión de precisión, sino de la fiabilidad de los resultados.

Los (principales) criterios geométricos

Las primeras métricas que debemos controlar son puramente geométricas. No nos proporcionan información sobre la física del problema, pero sí nos dicen mucho sobre la calidad matemática de nuestros volúmenes. En STAR-CCM+, estas métricas se calculan automáticamente y están disponibles mediante la herramienta Check Mesh.

Face Validity

La Face Validity mide si el centro de una celda está correctamente posicionado respecto a sus caras. En una celda de buena calidad, las normales de las caras apuntan hacia el exterior, alejándose del centroide de la celda. En una celda con una baja validez de caras, una o más normales apuntan hacia el interior, hacia el centroide de la celda.

Una Face Validity igual a 1.0 significa que todas las normales de las caras apuntan correctamente hacia el exterior, alejándose del centroide de la celda, mientras que valores inferiores a 0.5 indican celdas no convexas o, en el peor de los casos, celdas que se autointersectan (con volumen negativo).

Si tenéis celdas con una validez inferior a 1.0, el solucionador podría tener dificultades; si están por debajo de 0.8, podéis esperar problemas graves de convergencia.

Cell Quality

STAR-CCM+ utiliza un índice de calidad híbrido que combina diferentes factores, como la desviación de los vectores normales entre celdas adyacentes.

Una Cell Quality baja, cercana a cero, suele indicar celdas muy distorsionadas. Un buen analista intenta mantener la gran mayoría de las celdas por encima de 0.1, idealmente con valores mucho más elevados.

Skewness Angle

En términos sencillos, el Skewness Angle (θ) mide cuánto está «inclinada» una celda respecto a sus vecinas. Es el parámetro que nos indica si la información (presión, velocidad) puede transmitirse correctamente de un centro de celda a otro sin «perderse» ni generar inestabilidades matemáticas.

Imaginad dos celdas adyacentes. Para calcular cómo varían las magnitudes entre ellas, el software analiza dos vectores: el vector normal a la cara que las separa y el vector que une los centros de las dos celdas. El ángulo de skewness θ es el ángulo formado entre estos dos vectores. θ = 0° significa que la malla es perfectamente ortogonal. Cuanto mayor es el ángulo, mayor es la distorsión de la malla.

Aspect Ratio

El Aspect Ratio mide cuánto está «alargada» una celda. Imaginad un cubo perfecto: su relación de aspecto es 1:1. Imaginad ahora un espagueti: la relación es muy elevada.

En las zonas centrales del flujo (core flow), buscamos un Aspect Ratio cercano a 1. Sin embargo, en las capas límite, es normal y, a menudo, deseable tener valores elevados de Aspect Ratio (incluso 10:1 o superiores) para capturar los gradientes normales a la pared sin desperdiciar millones de celdas en la dirección paralela al flujo.

Volume Change (Variación de volumen)

Esta es quizás la métrica más subestimada. Representa la relación entre el volumen de una celda y el de sus celdas vecinas.

Los solucionadores odian los cambios bruscos: si se pasa de una celda milimétrica a una centimétrica de forma instantánea, el error de truncamiento en los gradientes aumentará considerablemente.

La física entra en juego: y+ e independencia de la malla

Una vez que la geometría está «sana», debemos preguntarnos: ¿la malla es adecuada para la física que queremos simular? Aquí las métricas ya no son únicamente geométricas, sino que pasan a depender de la solución obtenida.

El parametro y+: el termómetro de la capa límite

Al Wall y+ dedicaremos no uno, sino dos artículos específicos dentro de esta serie. Aquí resumimos sus principales efectos, pero os invitamos a leerlos en cuanto sean publicados.

Si estáis simulando turbulencia (y con toda probabilidad es así), el y+ es vuestro mejor aliado.

O vuestra peor pesadilla.

Es una coordenada adimensional que mide el nivel de refinamiento de la malla cerca de las paredes.

  • High-y+ (Wall Functions): 30 < y+ < 100 (ideal para ahorrar recursos computacionales).
  • Low-y+ (All y+ Wall Treatment): y+ ≈ 1 (necesario para capturar la subcapa viscosa).

Equivocarse con el y+ significa equivocarse en el cálculo de la resistencia y del intercambio de calor. STAR-CCM+ permite visualizar el y+ como una escena escalar sobre la superficie después de unas pocas iteraciones: ¡hacedlo siempre!

Mesh Sensitivity Study (Estudio de convergencia de la malla)

Ningún artículo sobre la validación estaría completo sin mencionar el estudio de independencia de la malla. El concepto es sencillo: si duplico el número de celdas y mi resultado (por ejemplo, la sustentación de un ala) cambia un 10%, significa que mi malla original no era lo suficientemente refinada.

El estándar de referencia es el uso del Grid Convergence Index (GCI), basado en la extrapolación de Richardson.

Donde Fs es un factor de seguridad, ε es el error relativo entre dos mallas, ρ es la relación de refinamiento y  es el orden de precisión del esquema numérico.

Antes de pulsar Run: la lista de comprobación "anti-crash" para no quemar el cálculo

Ya casi estamos. Las métricas geométricas están en verde, el y+ está dentro del rango adecuado, pero antes de dejar que el solucionador procese datos durante toda la noche (si todo va bien, porque existe el riesgo de que diverja después de 5 iteraciones y no os deis cuenta hasta la mañana siguiente), hay un último paso fundamental. Debemos asegurarnos de que la malla sea «topológicamente sana». Omitir estos tres controles finales significa exponerse a crashes inesperados o a una pérdida de precisión sutil pero devastadora.

  • La búsqueda de volúmenes negativos (elementos invertidos): Es el control más sencillo, pero también el más crítico. Si durante la generación de la malla (por ejemplo, para seguir un radio de redondeo muy estrecho o una característica geométrica de escala milimétrica) el mesher genera elementos con volumen cero o negativo, tendréis una bomba de relojería dentro del dominio. En cuanto el solucionador intente dividir por ese volumen en las ecuaciones discretizadas, se encontrará ante una división por cero y detendrá el cálculo con un error crítico. Elementos negativos = crash instantáneo.
  • La trampa de las interfaces no conformes: Cuando se conectan diferentes dominios — por ejemplo, la zona estacionaria de un compartimento motor con el dominio rotativo de un ventilador, o simplemente zonas con diferentes topologías de malla (por ejemplo, de hexaédricas a poliédricas) — STAR-CCM+ crea interfaces. Verificad siempre que no haya «huecos» o acoplamientos parciales. Si las superficies no coinciden geométricamente o si la tolerancia de acoplamiento es incorrecta, el flujo quedará literalmente bloqueado o, peor aún, se generarán recirculaciones artificiales que invalidarán la física del problema.
  • La alineación con el flujo (no ortogonalidad numérica): Esta, permitidme decirlo, es una consideración propia de analistas experimentados. En las zonas donde el fluido se desplaza a alta velocidad siguiendo una dirección preferente (como en el canal de un radiador o a lo largo de la carrocería), la orientación de las celdas es tan importante como su tamaño. Si las caras de los elementos están fuertemente desalineadas respecto a las líneas de flujo reales, los esquemas de discretización (incluso los de segundo orden o los métodos upwind avanzados) introducirán una difusión numérica Orientar la malla (utilizando, por ejemplo, celdas direccionales o extrusiones) de forma que las caras sean lo más paralelas o perpendiculares posible a los vectores de velocidad reduce el error artificial, proporcionando una solución mucho más nítida y cercana a la realidad.

Ingenieros y analistas pasan horas debatiendo sobre los detalles más avanzados de los modelos de turbulencia, las configuraciones de los solucionadores segregados o acoplados y las condiciones de contorno más complejas. Pero la realidad es que podéis configurar la física más sofisticada del mundo, pero si la malla que la sustenta es deficiente, simplemente estaréis acelerando hacia un resultado incorrecto.

Recordad: una malla deficiente con un solucionador de segundo orden siempre producirá resultados peores que una malla excelente con un solucionador de primer orden. Invertid tiempo en la malla, y la física os lo agradecerá.

Validar la malla no es un ejercicio académico para hacer más pesados los informes técnicos; es el primer paracaídas que tenéis para aseguraros de que los milisegundos consumidos por la CPU (y el coste de las licencias) estén calculando la realidad y no el error numérico del software.

Antes de lanzar el próximo cálculo pesado y abandonar la estación de trabajo durante la noche, abrid el panel de diagnóstico de STAR-CCM+, analizad los gráficos de distribución de calidad y haced una pregunta honesta: ¿esta malla realmente merece el tiempo de mi solucionador? Invertid hoy ese 10% de tiempo adicional en la malla. La convergencia de vuestros residuos, pero sobre todo la credibilidad de vuestros resultados mañana, os lo agradecerán.

Buen meshing.

Lee los demás artículos de la serie:

Comprueba el producto antes de fabricarlo.

¿Y si pudieras reducir hasta un 80 % los costes de tus prototipos?