Herramientas y criterios para mantener bajo control el Wall y+

Si trabajas en el mundo de la dinámica de fluidos computacional, sabes que hay pocas certezas en la vida: una de ellas es que el y⁺ (recuerda que se pronuncia “wai plus”) nunca es exactamente el que esperabas conseguir a la primera.

El Wall y⁺ es el parámetro adimensional que nos indica si la malla cerca de las paredes es adecuada para el modelo de turbulencia seleccionado; es el puente entre la realidad física de la capa límite y la abstracción numérica del software. Pero, ¿cómo se controla realmente en Simcenter STAR-CCM+? No basta con lanzar el cálculo y esperar. Se necesita una estrategia de evaluación rigurosa basada en mapas de contorno, histogramas y umbrales críticos.

Antes de nada, demos un paso atrás para repasar algunos conceptos básicos. La capa límite es donde ocurre toda la magia (y también el sufrimiento): fricción, desprendimiento de la vena fluida e intercambio térmico.

Para decidir cómo tratar la región cercana a la pared, hemos visto que en STAR-CCM+ existen dos enfoques:

• High y+: La primera celda se encuentra en la región turbulenta: en este caso buscamos un y+ > 30 (idealmente entre 30 y 300)
• Low y+: Estamos interesados en una mayor resolución dentro de la subcapa viscosa: en este caso buscamos un y+ ∼ 1

En realidad, existe una tercera opción: un modelo del software capaz de determinar en qué zona nos encontramos entre las dos anteriores y aplicar las Wall Functions adecuadas, denominado All-y+, disponible en STAR-CCM+. En este caso, resulta fundamental encontrarse fuera de la zona buffer, con valores de y+ comprendidos entre 10 y 30.

Equivocarse con el y+ significa equivocarse en el cálculo del drag, de la transferencia de calor o, peor aún, pasar por alto un desprendimiento del flujo que podría cambiar completamente los resultados. Afortunadamente, en STAR-CCM+ disponemos de todas las herramientas necesarias para evitar errores o, al menos, para identificar dónde se ha cometido un fallo antes de entregar los resultados.

Lo primero que hacemos todos es crear una Scalar Scene. En esta fase, la Scalar Scene interpreta el papel de la atalaya. Su función es proporcionar una visión general del problema, como una especie de reconocimiento aéreo que permite identificar a simple vista las macrozonas críticas. Visualizar el mapa permite comprender cuál es el estado global de la malla cerca de las paredes. Puede ocurrir que tengas un y+ excelente en el cuerpo principal del modelo, pero que, cerca de un radio de curvatura o de un borde de salida, la velocidad aumente y tu y+ se dispare.

Si estás buscando un enfoque Low-Re (donde y+ < 1), establece el valor máximo de la escala en 5 o 10. Todo lo que aparezca en rojo intenso estará fuera del objetivo. En cambio, si utilizas Wall Functions logarítmicas, configura un rango entre 30 y 300. Las zonas azules (por debajo de 30) son aquellas donde el modelo matemático puede presentar problemas, ya que la malla es demasiado refinada para las hipótesis en las que se basa la función de pared.

Veamos un caso práctico: imaginemos que queremos calcular la resistencia de un submarino sumergido utilizando un enfoque High y+. Por lo tanto, hemos dimensionado nuestra malla con un objetivo de Wall y+ comprendido entre 30 y 300 en toda la superficie.

Esta primera imagen es un clásico mapa de contornos: muy suave, con muchos niveles de color y transiciones graduales, visualmente agradable… pero no nos dice nada.

Esta segunda imagen, en cambio, presenta una distribución mucho menos suavizada, con zonas más definidas: sin duda es menos atractiva desde el punto de vista visual, pero permite identificar inmediatamente dónde se encuentran los problemas, destacando las celdas con y+ > 30.

Pro Tip: puedes superponer un Data Focus a las Scalar Scenes para resaltar zonas dentro de determinados rangos de una field function mientras visualizas otras magnitudes. Es una herramienta muy útil, por ejemplo, para analizar la velocidad en las zonas con bajo y+ y profundizar aún más en la física que se encuentra detrás de estos fenómenos.

Los mapas son bonitos, pero los revisores técnicos quieren números. ¿Cuánto porcentaje de tu superficie está fuera del objetivo? Para responder a esta pregunta, debemos pasar a herramientas más técnicas, como los Histogram Plots.

Para realizar un análisis objetivo e imparcial que represente el estado real de nuestra malla, necesitamos la inteligencia: el Histogram Plot.

La inteligencia no se deja llevar por las emociones. Toma todos los datos de cada una de las caras de la pared y los transforma en un informe estadístico frío y objetivo. Si la atalaya nos indica dónde se encuentran los problemas, el histograma nos dice hasta qué punto son graves.

El histograma proporciona una visión global de la calidad de la malla. Si observas un pico estrecho alrededor del valor objetivo, significa que has realizado un trabajo impecable. En cambio, si ves una distribución extendida que va desde 0.5 hasta 1000, tu malla no es uniforme y los resultados podrían ser inestables.

Pro Tip: Recuerda ponderar el histograma por área. De esta forma, las caras más grandes tienen un mayor peso en el análisis estadístico respecto a pequeños fragmentos de malla situados en esquinas o zonas localizadas. A la izquierda se muestra un histograma ponderado por área y, a la derecha, uno basado en frecuencia. Es fácil observar cómo la percepción del histograma de la derecha puede dar la impresión de un Wall y+ más extendido de lo que realmente es.

A veces el mapa resulta confuso y el histograma te indica que algo no funciona correctamente, pero no consigues identificar dónde está el problema. Aquí entra en juego nuestro francotirador (el Threshold), una de las herramientas más infravaloradas de STAR-CCM+. Su función es aislar y detectar únicamente las celdas que se comportan como elementos problemáticos, permitiendo identificarlas para su posterior corrección.

Verás iluminarse como por arte de magia las partes de la malla que presentan problemas. Es un método quirúrgico para identificar dónde debes intervenir con un Volumetric Control o dónde necesitas reducir el Base Size de la malla superficial.

Este enfoque es fundamental para los informes de validación: “El 95% de la superficie presenta un y+ > 30; las zonas en las que la primera celda cae dentro de la zona buffer (visibles mediante el Threshold) están localizadas únicamente en correspondencia con puntos de estancamiento y no afectan al drag global.”

Evaluar el y+ no es una tarea que deba realizarse al final del proceso, sino que forma parte de un procedimiento iterativo. A menudo requiere volver al Mesher, ajustar el Prism Layer Thickness o el número de capas prismáticas y volver a calcular.

En STAR-CCM+, gracias a la canalización automatizada, puedes modificar un parámetro del mesher y observar cómo se actualiza el histograma en tiempo real después de unas pocas iteraciones de cálculo. No te conformes con un valor medio global; utiliza mapas, histogramas y thresholds para garantizar que la física de tu simulación sea tan sólida como la geometría de tu modelo.

La buena noticia es que, para geometrías no excesivamente complejas y con la dosis adecuada de experiencia y conocimiento, existen numerosos métodos de cálculo que permiten estimar a priori el espesor de la primera celda. De esta forma, es posible reducir al mínimo (y, con algo de suerte, incluso eliminar) el número de iteraciones necesarias para alcanzar un dimensionamiento adecuado de la capa límite.

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

• Wall y+: cómo controlar realmente la malla cerca de las paredes

• Análisis de Convergencia de Malla: ¿cuándo una simulación CFD se convierte en un resultado fiable?

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