La estructura asintótica de la capa límite determina una malla fuertemente no uniforme: utilice celdas muy pequeñas normales a las paredes sólidas, donde los gradientes de velocidad y escalares son pronunciados, y celdas progresivamente mayores en el flujo exterior, donde los gradientes son más débiles. La resolución necesaria cerca de la pared depende del solucionador y del tratamiento de pared: la DNS resuelve directamente las escalas viscosas y turbulentas, la LES requiere una resolución adecuada de las estructuras cercanas a la pared o un modelo de pared, y la RANS suele establecer como objetivo una coordenada de pared específica, como (y^+). En una planta piloto, los caudales, las caídas de presión, las temperaturas y las condiciones de pared medidas proporcionan la evidencia necesaria para validar esa estrategia de mallado.
A números de Reynolds elevados, el flujo no presenta la misma dificultad de resolución en todas partes. La capa límite crea una región viscosa delgada que controla el esfuerzo cortante en la pared, la pérdida de presión, la separación y la transferencia de calor o masa; por tanto, el diseño de la malla debe asignar la resolución de acuerdo con las escalas locales de la capa límite y verificarla mediante mediciones de la planta piloto.
Por qué la capa límite controla el diseño de la malla
El flujo exterior no necesita una resolución fina uniforme
A números de Reynolds elevados, el flujo puede dividirse conceptualmente en una región exterior, aproximadamente no viscosa, y una delgada capa límite viscosa adyacente a la pared.
Los gradientes de velocidad normales a la pared son pequeños en gran parte de la región exterior. Una malla uniformemente fina en esa zona desperdicia recursos computacionales sin mejorar necesariamente las magnitudes de interés.
La región de la pared contiene los gradientes importantes
La condición de no deslizamiento exige que la velocidad alcance el valor cero en la pared. Después, esta aumenta rápidamente hasta aproximarse a la velocidad del flujo exterior a través de la capa límite.
Esto genera valores elevados de (\partial u/\partial y), que determinan directamente el esfuerzo cortante en la pared:
[ \tau_w = \mu \left.\frac{\partial u}{\partial y}\right|_w ]
Una malla demasiado gruesa en la dirección normal a la pared suaviza numéricamente este gradiente. Los errores resultantes pueden manifestarse como factores de fricción, caídas de presión, arrastre, distribuciones del tiempo de residencia o coeficientes de transferencia de calor incorrectos.
El espesor de la capa límite disminuye con el número de Reynolds
Una relación de escala característica es
[ \frac{\delta}{L} \sim \frac{1}{\sqrt{Re}} ]
para una capa límite laminar bajo condiciones apropiadas. Por tanto, al aumentar el número de Reynolds, la capa límite generalmente se vuelve más delgada en relación con la longitud del equipo.
La consecuencia práctica es directa: las simulaciones con números de Reynolds más elevados generalmente requieren un refinamiento más agresivo en la dirección normal a la pared, a menos que se utilice una estrategia adecuada de funciones de pared o de modelo de pared.
Traducir la estructura a una malla
Refinar primero en dirección normal a la pared
La primera prioridad es el espaciado normal a la pared. Varias celdas deben quedar distribuidas a través de la capa límite, con una colocación especialmente cuidadosa cerca de la pared, donde el gradiente de velocidad es máximo.
Una malla típica cerca de la pared utiliza capas de inflación o prismáticas cuyo espesor aumenta gradualmente al alejarse de la pared. Los saltos bruscos entre celdas diminutas en la pared y celdas grandes en el exterior pueden reducir la precisión y crear problemas de convergencia.
Estirar la malla al alejarse de la pared
La malla debe expandirse suavemente desde la pared hacia el flujo exterior. Esto refleja la menor necesidad de resolución a medida que el perfil de velocidad se aproxima al valor del flujo externo.
La razón de expansión debe mantenerse suficientemente moderada para preservar la calidad de la discretización. Un estiramiento excesivo puede producir errores de interpolación, estimaciones deficientes de los gradientes y sensibilidad del solucionador.
Resolver el desarrollo en la dirección del flujo
Las capas límite crecen aguas abajo, por lo que la extensión normal a la pared y la distribución de celdas necesarias cambian a lo largo del equipo.
La malla en la dirección del flujo también debe resolver:
- El desarrollo de la región de entrada.
- Los cambios en el gradiente de presión.
- La separación y la readherencia.
- La mezcla, los codos, las contracciones y las expansiones.
- Los cambios locales en la transferencia de calor o masa.
Una malla diseñada únicamente a partir de la capa límite de entrada puede resultar inadecuada cerca de un codo, una válvula, un impulsor, un deflector o un cambio brusco de área.
Añadir resolución transversal cuando importen los efectos tridimensionales
En flujos rectos e idealizados, la resolución transversal puede ser menos exigente que la resolución normal a la pared. Sin embargo, en las plantas piloto son habituales los efectos tridimensionales.
Los codos, colectores, placas distribuidoras, haces tubulares, equipos rotativos y accesorios asimétricos pueden crear flujos secundarios y gradientes transversales. Estos requieren una resolución suficiente en las direcciones transversales, en lugar de depender únicamente de una suposición bidimensional de capa límite.
Adaptar la malla al método CFD
Simulaciones RANS
Los modelos RANS suelen utilizar mallas con paredes resueltas o funciones de pared.
Una malla RANS con paredes resueltas coloca la primera celda suficientemente cerca de la pared para resolver el perfil de velocidad cercano a ella. Una malla con funciones de pared coloca la primera celda en un intervalo apropiado de la capa logarítmica y modela la región más próxima a la pared.
El objetivo correcto no es simplemente «lo más fino posible». Es el intervalo de (y^+) requerido por el modelo de turbulencia y el tratamiento de pared seleccionados. Una malla que no sea coherente con dicho tratamiento puede producir predicciones poco fiables del esfuerzo cortante en la pared y de la transferencia de calor.
Simulaciones LES
La LES resuelve las estructuras turbulentas grandes y modela las más pequeñas. Cerca de las paredes sólidas, las escalas turbulentas más pequeñas resultan especialmente costosas de resolver.
Por consiguiente, la LES puede requerir:
- Una resolución muy fina cerca de la pared para una LES con paredes resueltas.
- Un modelo de pared validado para simulaciones de equipos a escala de planta y con números de Reynolds elevados.
La malla también debe ser suficientemente fina en las direcciones del flujo y transversal para representar las estructuras turbulentas resueltas. Refinar únicamente la dirección normal a la pared no es suficiente.
Simulaciones DNS
La DNS resuelve todas las escalas espaciales y temporales relevantes sin utilizar un modelo de turbulencia. Sus requisitos de malla son mucho más exigentes, especialmente cerca de las paredes y a números de Reynolds elevados.
La DNS suele ser adecuada para estudios fundamentales o geometrías representativas pequeñas, no para equipos de plantas piloto a escala completa en condiciones industriales de operación. No obstante, puede proporcionar datos de referencia para validar modelos de orden reducido o tratamientos de pared.
Utilizar las mediciones de la planta piloto para definir y validar la malla
Comenzar con condiciones de operación medibles
El diseño de la malla debe basarse en el intervalo real de operación, no únicamente en la geometría nominal. Entre las entradas importantes se incluyen:
- Caudal y velocidad.
- Densidad, viscosidad y propiedades térmicas del fluido.
- Niveles de turbulencia o perturbación en la entrada.
- Temperatura de la pared o flujo de calor.
- Rugosidad superficial.
- Presión y temperatura de operación.
- Condiciones transitorias de puesta en marcha o ciclado.
Estos factores determinan los números de Reynolds locales, los espesores de las capas límite y las regiones donde el refinamiento es más importante.
Validar magnitudes que dependan de la resolución cerca de la pared
Los datos de la planta piloto deben compararse con resultados de simulación sensibles a la capa límite, incluidos:
- Caída de presión a lo largo de tuberías y equipos.
- Factor de fricción.
- Esfuerzo cortante en la pared.
- Perfiles de velocidad.
- Ubicaciones de separación y readherencia.
- Coeficientes de transferencia de calor.
- Perfiles de concentración o transferencia de masa.
- Comportamiento de la mezcla y del tiempo de residencia.
La coincidencia únicamente en el caudal global no demuestra que la malla cercana a la pared sea adecuada.
Utilizar estudios de convergencia de malla
Un estudio de malla debe comparar al menos tres mallas relacionadas sistemáticamente, por ejemplo, gruesa, media y fina. La comparación debe centrarse en los resultados de ingeniería, no solo en la reducción de los residuos.
Una simulación puede converger numéricamente y seguir dependiendo de la malla. Si las predicciones de caída de presión, esfuerzo cortante en la pared o transferencia de calor cambian sustancialmente con el refinamiento, la malla aún no demuestra independencia de la solución.
Combinar la incertidumbre numérica y experimental
Las mediciones de la planta piloto contienen incertidumbre procedente de los sensores, la calibración, las perturbaciones del flujo y la variabilidad operativa. Los resultados numéricos contienen incertidumbre de discretización, de modelado y de las entradas.
La pregunta adecuada no es si la simulación y el experimento son idénticos punto por punto. Es si la diferencia es coherente con la incertidumbre combinada y si las mismas tendencias se mantienen en distintas condiciones de operación.
Considerar las fronteras, la separación y los transitorios
Distinguir entre regiones de flujo adherido y separado
Las suposiciones de capa límite son más fiables para flujos predominantemente adheridos con cambios geométricos graduales. Los fuertes gradientes de presión adversos pueden provocar separación, en la que la capa límite pierde su adherencia a la pared y el flujo se vuelve altamente inestable.
Las zonas de separación requieren refinamiento local en más de una dirección. La estela, la capa de cizallamiento, la región de recirculación y el punto de readherencia pueden controlar tanto la pérdida de presión como la mezcla.
Refinar alrededor de las perturbaciones geométricas
Las plantas piloto contienen características que invalidan la representación simple de una «capa delgada a lo largo de una pared plana». Se debe refinar alrededor de:
- Entradas y salidas.
- Codos y derivaciones en T.
- Válvulas y accesorios.
- Expansiones y contracciones bruscas.
- Deflectores y elementos internos.
- Entradas de tubos.
- Impulsores y componentes rotativos.
- Distribuidores porosos.
Estas regiones generan gradientes intensos, flujo secundario o recirculación que pueden dominar el resultado a escala del equipo.
Resolver las escalas temporales relevantes
Durante la puesta en marcha, la parada, la pulsación o la operación periódica, la capa límite evoluciona con el tiempo. El paso temporal debe resolver el transitorio relevante y ser compatible con la malla espacial.
Magnitudes adimensionales como el número de Strouhal y el número de Reynolds ayudan a determinar si la variación del proceso es lenta en relación con el desarrollo del flujo o suficientemente rápida para producir una respuesta transitoria diferenciada de la capa límite.
Comprender las compensaciones
Las mallas más finas aumentan rápidamente el coste
El refinamiento cerca de la pared incrementa el número de celdas, la demanda de memoria y el tiempo de cálculo. Las capas de inflación tridimensionales pueden hacer que este coste sea especialmente elevado en tuberías largas, intercambiadores de calor y equipos de proceso complejos.
Por tanto, el objetivo es una resolución específica, no la máxima resolución en todas partes.
Las funciones de pared reducen el coste, pero añaden supuestos
Las funciones de pared pueden hacer viables las simulaciones a números de Reynolds elevados al modelar la región cercana a la pared que no se resuelve. Son menos fiables cuando el flujo presenta una separación intensa, cambios rápidos del gradiente de presión, efectos de rugosidad, curvatura pronunciada o una transferencia significativa de calor y masa.
Su uso debe ser coherente con la ubicación de la malla y el régimen físico.
Las simplificaciones de capa límite no son universalmente válidas
Las ecuaciones reducidas de la capa límite pueden tratarse como parabólicas en flujos adecuados con números de Reynolds elevados y predominantemente aguas abajo. Esto permite procedimientos eficientes de avance, ya que la información se propaga principalmente aguas abajo.
Sin embargo, se requiere un tratamiento completo de Navier–Stokes cuando son importantes la influencia aguas arriba, la recirculación, la fuerte tridimensionalidad, la separación o el acoplamiento elíptico de la presión. Una planta piloto no debe simplificarse únicamente porque su número de Reynolds global sea elevado.
El orden formal no garantiza la precisión de ingeniería
Una discretización de orden superior puede reducir el error de truncamiento, pero no puede compensar una colocación inadecuada cerca de la pared, una representación geométrica deficiente, un modelado inapropiado de la turbulencia o una resolución insuficiente de la separación.
Del mismo modo, añadir celdas sin realizar un estudio de convergencia no demuestra que el resultado sea fiable.
Cómo aplicar esto a su proyecto
Un flujo de trabajo práctico consiste en diseñar la malla en torno al intervalo de operación y a los resultados sensibles a la pared que sean relevantes para la operación unitaria.
- Si su objetivo principal es la caída de presión o el esfuerzo cortante en la pared: Priorice el refinamiento normal a la pared, seleccione una estrategia coherente de (y^+) y valide los resultados frente a las pérdidas de presión y los factores de fricción medidos.
- Si su objetivo principal es la transferencia de calor o masa: Resuelva tanto las capas límite hidrodinámica y térmica como la de concentración, que pueden tener espesores diferentes, y compare los coeficientes de transferencia previstos con los datos de la planta piloto.
- Si su objetivo principal es la mezcla o el tiempo de residencia: Refine las estelas, las zonas de recirculación, las entradas, los elementos internos y las regiones tridimensionales de flujo transversal, en lugar de concentrar la resolución únicamente en las paredes rectas.
- Si su objetivo principal es la separación o la operación no estacionaria: Utilice un refinamiento local multidireccional, pasos temporales suficientemente pequeños y mediciones que capturen la presión, la velocidad o el comportamiento transitorio de readherencia.
- Si su objetivo principal es la eficiencia computacional: Utilice una malla estratificada y refinada localmente, y seleccione RANS, LES o DNS según la física requerida, en lugar de aplicar el método más costoso en todas partes.
- Si su objetivo principal es la validación del modelo: Realice estudios de convergencia de malla y compare múltiples condiciones de operación, no solo un punto de datos de la planta piloto.
Una malla basada en la estructura asintótica de la capa límite concentra el esfuerzo computacional donde la física —y las decisiones de ingeniería— son más sensibles.
Tabla resumen:
| Prioridad del diseño de la malla | Enfoque recomendado | Resultados clave que deben validarse |
|---|---|---|
| Resolución cerca de la pared | Utilice capas prismáticas o de inflación y seleccione un tratamiento de pared coherente con el objetivo de y+ |
Esfuerzo cortante en la pared, factor de fricción, caída de presión |
| Malla del flujo exterior | Aumente gradualmente el tamaño de las celdas donde los gradientes de velocidad y escalares sean más débiles | Desarrollo global del flujo y eficiencia |
| Perturbaciones geométricas | Refine codos, válvulas, entradas, salidas, elementos internos, estelas y zonas de recirculación | Separación, mezcla, tiempo de residencia |
| Método CFD | Adapte la resolución de la malla a los requisitos de RANS, LES o DNS | Fiabilidad de la turbulencia y del modelo de pared |
| Transferencia de calor y masa | Resuelva las capas límite hidrodinámica, térmica y de concentración | Coeficientes y perfiles de transferencia |
| Validación | Compare mallas gruesas, medias y finas con las mediciones de la planta piloto | Independencia de la malla y concordancia de las incertidumbres |
¿Está planificando un estudio CFD fiable para sus equipos de proceso? Póngase en contacto con LABPARK para analizar sus necesidades. LABPARK proporciona plantas piloto educativas y de formación para operaciones unitarias en ingeniería química, bioprocesos y biotecnología, así como en aplicaciones ambientales y de tratamiento de aguas. Nuestras plantas piloto ayudan a universidades, institutos de investigación y empresas a evaluar la caída de presión, la transferencia de calor y masa, la mezcla, el tiempo de residencia, la separación y el comportamiento operativo transitorio en condiciones medibles. Podemos ayudarle a seleccionar los equipos pertinentes, definir las condiciones de operación, apoyar estrategias de mallado basadas en la capa límite y conectar los resultados de simulación con una validación experimental práctica.