La definición matemática del y+ (por cierto, se pronuncia “wai plus”) tiene su origen en la teoría de la capa límite. Se introduce haciendo adimensional la distancia normal a la pared y respecto a la escala viscosa.
Wall y+: cómo controlar realmente la malla cerca de las paredes
Hoy hablamos de uno de esos parámetros que, cuando se configuran incorrectamente, pueden transformar tu precioso contour de colores en una obra de arte abstracta completamente desprovista de significado físico: el Wall y+, uno de los parámetros más debatidos y, en ocasiones, menos comprendidos.
Comprender en profundidad el y+ permite diseñar mallas que ayuden a los modelos de turbulencia a funcionar según lo previsto, garantizando que la dinámica de la capa límite se capture con la precisión necesaria.
La arquitectura del y+: derivación analítica y escalado adimensional
¿Por qué es tan importante el y+?
Como toda adimensionalización, se trata de una relación entre una magnitud (en nuestro caso y, la distancia a la pared) y una magnitud de referencia coherente: la escala de longitud viscosa v/uτ, donde v es la viscosidad cinemática del fluido y uτ es la velocidad de fricción en la pared, definida como:
con τw el wall shear stress e ρ la densidad del fluido.
Vale, pero ¿no bastaban los milímetros?
El y+ es una distancia adimensional. No nos indica a cuántos milímetros estamos de la pared, sino «qué tan cerca» nos encontramos en relación con las características del fluido (viscosidad) y con el esfuerzo que el propio fluido ejerce sobre la pared.
Echemos otro vistazo al y+: se parece mucho a la definición del número de Reynolds, ¿verdad? Y es precisamente por eso: además de la similitud formal entre ambas definiciones, al igual que el número de Reynolds, el y+ representa la relación entre las fuerzas de inercia y las fuerzas viscosas cerca de la pared.
Cuando hablamos de Wall y+, nos referiremos al valor de y+ correspondiente al centroide de la primera celda (la más cercana a la superficie) de la capa límite.
Dentro la capa límite: la fenomenología del fluido detrás del número
Imagina el aire fluyendo sobre una superficie, como el ala de un avión o el capó de un Ferrari. Justo cerca de la pared, el fluido se ralentiza bruscamente debido a la viscosidad, hasta detenerse completamente sobre la superficie (la famosa condición de no deslizamiento o no-slip condition). La velocidad del fluido varía de forma continua desde un valor nulo en la pared hasta alcanzar la velocidad de la corriente libre U (estrictamente hablando, hasta el 99% del valor de la corriente libre). Esta región es la conocida capa límite, caracterizada por su particular perfil de velocidad.
Dentro de la capa límite, es posible distinguir tres zonas principales en función de la predominancia de las fuerzas viscosas o de las fuerzas inerciales:
- Subcapa viscosa (Viscous Sublayer): en esta región, los efectos de la viscosidad molecular son dominantes. El flujo tiende a estar alineado con la pared y la turbulencia es prácticamente inexistente.
- Capa buffer: una zona de transición compleja en la que ni las leyes viscosas ni las leyes turbulentas son completamente válidas. La turbulencia tiende a disiparse mediante dinámicas complejas.
- Región completamente turbulenta (log-layer): en esta zona, la turbulencia está plenamente desarrollada.
La inteligencia: la relación estadística del histograma
Llegados a este punto, seguramente se habrá encendido una bombilla (o una luz de emergencia con sirena incluida): si tenemos regiones caracterizadas por el predominio de esfuerzos viscosos o fuerzas inerciales, para determinar en qué zona de la capa límite se encuentra un punto concreto necesitamos una herramienta que actúe como puente entre la posición y la relevancia de los esfuerzos viscosos.
El y+ sirve precisamente para indicar en qué zona de la capa límite nos encontramos y, sobre todo, como veremos más adelante, si la malla es lo suficientemente refinada como para comprender qué ocurre en esa región.
- Si y+ < 5, nos encontramos dentro de la subcapa viscosa, donde los perfiles de velocidad siguen una relación lineal respecto a la coordenada normal a la pared (u+ = y+).
Si el objetivo es calcular con precisión el coeficiente de fricción (Cƒ) o el número de Nusselt asociado al intercambio térmico, la malla debe estar resuelta hasta alcanzar valores de y+ ≈ 1.
- Cuando y+ se encuentra aproximadamente entre 10 y 30, la celda se sitúa en la capa buffer, una zona de transición en la que ni los efectos viscosos ni los turbulentos son completamente dominantes.
Es la región más crítica para la estabilidad de los solvers numéricos.
Per y+ > 30, la celda se encuentra en la región logarítmica, donde el perfil medio de velocidad sigue la ley logarítmica de la pared.
En esta zona, es posible utilizar aproximaciones matemáticas para reducir el coste computacional.
Cuando explico el tema del Wall y+ a los estudiantes de nuestros cursos, en este punto siempre hago una pausa durante unos instantes. Tanto para comprobar que todos siguen despiertos, como porque, puntualmente, llega la pregunta habitual y completamente legítima:
Pero, ¿por qué nos interesa saber en qué zona de la capa límite nos encontramos?
Wall Functions: el atajo inteligente para evitar un cálculo infinito
Para calcular correctamente el perfil de velocidad, si quisiéramos resolverlo directamente, necesitaríamos un número extremadamente elevado de celdas muy finas dentro de la capa límite.
El gradiente de velocidad en la pared, del que deriva el esfuerzo cortante y, por tanto, la resistencia por fricción, es muy pronunciado; si las celdas cercanas a la pared son demasiado gruesas, el gradiente de velocidad en la pared no se resuelve correctamente.
Como consecuencia, el esfuerzo cortante en la pared τw se evalúa de forma incorrecta, lo que afecta directamente a la predicción de la resistencia aerodinámica del objeto.
Si, por ejemplo, quisiera calcular de forma directa la capa límite de un vehículo de competición a 200 km/h sin (o casi sin) utilizar funciones de pared, necesitaría una dimensión de celda de 10⁻⁷ m (0,7 μm) para alcanzar un Wall y+ = 0,1. Y no solo eso: el número de celdas aumentaría de forma exponencial también sobre la superficie, donde sería necesario un refinamiento extremo para mantener los Aspect Ratios de las celdas dentro de los valores óptimos.
Si quisiéramos resolver físicamente cada capa límite de una carrocería completa de un vehículo de competición hasta la escala molecular, la simulación probablemente terminaría después de la jubilación del analista.
Y, en realidad, la cuestión sería aún más compleja, mucho más compleja, al introducir el concepto de viscosidad turbulenta, un tema que, sin embargo, va mucho más allá de lo que queremos abordar hoy.
Un razonamiento similar se aplica a la transferencia de calor: cuando la resolución cerca de la pared es insuficiente, el coeficiente local de transferencia de calor también se calcula de forma incorrecta, ya que el gradiente de temperatura en la pared no puede capturarse con precisión.
Entonces, ¿nos rendimos? Evidentemente no.
Aquí es donde entra en juego uno de los trucos más conocidos de la CFD: las Wall Functions.
En lugar de obligar al ordenador a calcular con esfuerzo el perfil de velocidad punto por punto dentro de la subcapa viscosa, decidimos hacer trampa de forma inteligente utilizando un modelo no lineal basado en funciones de pared, que calcula correctamente el gradiente. Esto permite utilizar un menor número de celdas dentro de la capa límite y evita Aspect Ratios excesivamente elevados, garantizando así una mayor estabilidad numérica.
Utilizando las Wall Functions, le decimos al software: «No pierdas tiempo analizando qué ocurre en esa región; sabemos por teoría que, si la primera celda se encuentra a una determinada distancia de la pared, la velocidad seguirá esta curva logarítmica específica». Este enfoque permite «saltar» directamente la región más costosa de la capa límite, aplicando una condición de contorno basada en un perfil de velocidad predefinido.
Es un compromiso brillante: sacrificamos una resolución geométrica infinitesimal para obtener un cálculo miles de veces más rápido, confiando en la validez universal de la física de Prandtl y von Kármán para obtener resultados que siguen siendo increíblemente precisos en la mayoría de las aplicaciones industriales.
Y, por último, aquí es donde entra en juego la última pieza del rompecabezas: un fantástico gráfico que todo analista de fluidodinámica ha visto al menos una vez en su vida y que, espero, ahora tenga mucho más sentido.
En función del valor de y+, es necesario utilizar una ley de tipo lineal (que en el gráfico log-log presenta la forma de una exponencial) o una ley logarítmica (que en el gráfico aparece como una recta).
De forma simplificada, el procedimiento es el siguiente: se calcula la distancia a la pared del centroide de la celda donde se realiza el cálculo y se adimensionaliza para obtener . En función de este valor, se determina si nos encontramos en la subcapa viscosa o en la outer region, y se aplica la ley de pared correspondiente. De esta forma, se calcula y, a partir de él, se obtiene el perfil de velocidad y su derivada, es decir, el esfuerzo cortante en la pared.
De esta forma, en realidad, no estamos resolviendo las ecuaciones de la capa límite, sino que las estamos modelando, utilizando mallas más gruesas y aplicando otros recursos relacionados con la viscosidad, propios de los modelos de turbulencia.
De esta forma, en realidad, no estamos resolviendo las ecuaciones de la capa límite, sino que las estamos modelando, utilizando mallas más gruesas y aplicando otros recursos relacionados con la viscosidad, propios de los modelos de turbulencia.
Si has llegado hasta este punto, una imagen de este tipo debería tener ahora un significado completamente diferente para ti, y habrás comprendido que en un análisis fluidodinámico no se puede prescindir de una correcta modelización de la capa límite basada en el y+.
En este punto, seguramente tendrás muchas preguntas; de hecho, probablemente muchas más que antes si este artículo ha cumplido su objetivo.
¿Cómo dimensiono la capa límite para obtener el valor correcto de y+?
¿Cómo puedo controlarlo en mi simulación?
¿Qué valores puedo aceptar y en qué zonas?
¡Para responder a estas preguntas tendréis que leer el próximo artículo!
Comprueba el producto antes de fabricarlo.
¿Y si pudieras reducir hasta un 80 % los costes de tus prototipos?
