A numerical study of turbulent heat transfer in pipes

Este trabajo está orientado al estudio detallado del problema de transferencia de calor turbulenta en un tubo sujeto a un flujo de calor no uniforme en la superficie y recorrido por un fluido con propiedades dependientes de la temperatura. El estudio de este problema busca ser representativo de las condiciones de trabajo de los tubos en un receptor solar de torre. La operación de estas plantas de energía solar térmica por concentración depende en gran manera de obtener una solución precisa del problema térmico, muchas veces por encima de lo que proporcionan las correlaciones actualmente en uso. La resolución de este problema es, además, de utilidad para otras muchas aplicaciones en ingeniería El principal objetivo de esta tesis es contribuir al avance en el entendimiento de la influencia de las condiciones de contorno térmicas no uniformes y de las propiedades del fluido variables en las estadísticas del flujo turbulento dentro del tubo. Con este fin, se ha generado una base de datos numérica modificando los parámetros del problema: el número de Reynolds de fricción (Re_ = 180 − 360), el número de Prandtl (Pr = 0.7 − 4), la distribución de radiación en la superficie y la sensibilidad de la viscosidad y la difusividad del fluido con la temperatura. En el estudio hemos observado que, mientras que el número de Reynolds tiene poco efecto en la distribución de temperatura, el impacto del número de Prandtl es significativo. Este resultado es de importancia para el diseño de los receptores solares ya que su operación tiene limitado el rango de temperaturas de funcionamiento. Otro resultado notable del estudio es la observación del efecto que tiene el calentamiento o enfriamiento de la pared en un caso con propiedades variable del fluido. Debido al cambio de viscosidad con la temperatura, se obtiene que la fluctuaciones de temperatura incrementan cerca de la superficie fría y se amortiguan en la superficie caliente. La variación de propiedades con la temperatura produce además la generación de pequeñas pero distinguibles velocidades medias en el plano del tubo. Estas velocidades secundarias son responsables de una contribución no desdeñable del flujo de calor total. Con objeto de profundizar en el fenómeno del flujo de calor turbulento en el tubo, se ha realizado un descomposición modal de los campos de velocidad y temperatura instantáneos, extrayendo las estructuras más energéticas en el flujo. Usando la técnica de EPOD (Extended Proper Orthogonal Decomposition), se han identificado las estructuras turbulentas responsables de la transferencia de calor. Se ha observado que, en el rango considerado, el tamaño de estas estructuras es poco dependiente del número de Reynolds pero está muy afectado por el número de Prandtl.The present work is devoted to the detailed study of the turbulent heat transfer problem in a pipe flow with variable fluid properties subjected to highly non-uniform heat flux distribution on the pipe surface. This problem aims to be representative of the conditions found in the tubes of heat receivers in Solar Power Tower (SPT) plants where sun radiation is converted into electric energy. The operational and economical performance of SPT plants rely on an accurate prediction of the thermal field, sometimes beyond the current semi-empirical approach to the heat transfer problem. Moreover, this study might be of interest for many other heat transfer applications in engineering. The primary intention of this thesis is to shed some light on the influence of the circumferentially varying heat flux conditions and of the temperature-dependent fluid properties in the statistics of a fully-developed turbulent flow. To that end, we create and analyze a numerical database modifying the main flow parameters of interest: the friction Reynolds number (Re_ = 180 − 360), the Prandtl number (Pr = 0.7 − 4), the heat flux distribution and the sensitivity to temperature of the fluid variables. From the analysis, we observe that, while the friction Reynolds number has little effect on the temperature distribution on the pipe wall, the impact of the Prandtl number is significant. This is of importance for the heat receiver application, as the maximum flow temperature at the inner wall, the so-called film temperature, might be always below certain limit. Another important outcome from the analysis of a case with a heated half of the pipe but cooled in the other half is that, because of the influence of the variable fluid properties, the temperature turbulent fluctuations are enhanced near the cold wall but damped near the hot wall. This behavior is seen to reduce the maximum temperature peaks in the flow while maintaining a similar wall temperature distribution. The variation of the fluid properties are also shown to induce small but discernible secondary velocities in the pipe cross-plane. These velocities are responsible of a non negligible contribution to the heat flux from the heated wall to the cooled wall. To gain a deeper insight in the turbulent heat transfer phenomena, we perform a modal decomposition of the instantaneous velocity and temperature fields to extract the information on the most energetic coherent structures in the flow. The use of a relatively new technique, Extended Proper Orthogonal Decomposition, allows us to discern how the velocity fluctuations are correlated to the temperature fluctuations, hence increasing our knowledge in the heat transfer process. From this analysis, we characterize the modes that are responsible of the heat transport in the pipe crossplane, obtaining that the scales of the structures bringing hot fluid from the heated wall to the pipe core are strongly affected by the Prandtl number.Programa Oficial de Doctorado en Ingeniería MatemáticaPresidente: Manuel Rodríguez Fernández.- Secretario: Carolina Marugán Cruz.- Vocal: Sergio Hoyas Calv