Three-dimensional numerical simulation of fluid flow and heat transfer in fin-and-tube heat exchangers at different flow regimes

This thesis aims at unifying two distinct branches of work within the Heat Transfer Technological Center (CTTC). On one side, extensive experimental work has been done during the past years by the researchers of the laboratory. This experimental work has been complemented with numerical models for the calculation of fin and tube heat exchangers thermal and fluid dynamic behavior. Such numerical models can be referred to as fast numerical tool which can be used for industrial rating and design purposes. On the other hand, the scientists working at the research center have successfully developed a general purpose multi-physics Computational Fluid Dynamics (CFD) code (TermoFluids). This high performance CFD solver has been extensively used by the co-workers of the group mainly to predict complex flows of great academic interest. The idea of bringing together this two branches, comes from the necessity of a reliable numerical platform with detailed local data of the flow and heat transfer on diverse heat exchanger applications. Being able to use local heat transfer coefficients as an input on the rating and design tool will lead to affordable and accurate prediction of industrial devices performance, by which the center can propose enhanced alternatives to its industrial partners. To accomplish these goals, several contributions have been made to the existing TermoFluids software which is in continuous evolution in order to meet the competitive requirements. The most significant problematics to adequately attack this problem are analyzed and quite interesting recommendations are given. Some of the challenging arising issues involve the generation of suitable and affordable meshes, the implementation and validation of three dimensional periodic boundary condition and coupling of different domains with important adjustments for the study of cases with different flow physics like time steps and thermal development. Turbulence is present in most of engineering flows, and refrigeration evaporator heat exchangers are not an exception. The presence of many tubes (acting like bluff bodies for the flow) arranged in different configurations and the fact that the flow is also confined by fins, create complex three dimensional flow features that have usually turbulent or transition to turbulent regime. Therefore, three dimensional turbulent forced convection in a matrix of wall-bounded pins is analyzed. Large Eddy Simulations (LES) are performed in order to assess the performance of three different subgrid-scale models, namely WALE, QR and VMS. The Reynolds numbers of the study were set to 3000, 10000 and 30000. Some of the main results included are the pressure coefficient around the cylinders, the averaged Nusselt number at the endwalls and vorticity of the flow. The final part of the thesis is devoted to study the three dimensional fluid flow and conjugated heat transfer parameters encountered in a plate fin and tube heat exchanger used for no-frost refrigeration. The numerical code and post processing tools are validated with a very similar but smaller case of a heat exchanger with two rows of tubes at low Reynolds for which experimental data is available. The next analysis presented is a typical configuration for no-frost evaporators with double fin spacing (for which very few numerical data is reported in the scientific literature). Conjugated convective heat transfer in the flow field and heat conduction in the fins are coupled and considered. The influence of some geometrical and flow regime parameters is analyzed for design purposes. In conclusion, the implementations and general contributions of the present thesis together with the previous existent multi-physics computational code, has proved to be capable to perform successful top edge three dimensional simulations of the flow features and heat transfer mechanisms observed on heat exchanger devices.Esta tesis tiene como objetivo unificar dos ramas de trabajo dentro del Centro Tecnológico de Transferencia de Calor (CTTC). Por un lado, se ha realizado un amplio trabajo experimental durante los últimos años. Este trabajo experimental se ha complementado con modelos numéricos para el estudio de intercambiadores de calor de tipo aleta-tubo. Tales modelos numéricos pueden considerarse una herramienta numérica de bajo coste empleada con propósitos de diseño principalmente. Por otro lado, los científicos que trabajan en el centro han desarrollado con éxito un código de Dinámica de Fluidos Computacionales (TermoFluids). Este código de alto rendimiento ha sido ampliamente utilizado principalmente para predecir flujos complejos de gran interés académico. La idea de unir a estas dos ramas, proviene de la necesidad de una plataforma numérica fiable con datos locales propios del flujo y de la transferencia de calor en diversas aplicaciones de intercambiadores de calor. Ser capaz de generar coeficientes locales de transferencia de calor para abastecer con datos propios los modelos existentes de bajo coste, permitirá la correcta predicción del rendimiento de dichos dispositivos. Para lograr estos objetivos, se han hecho varias contribuciones al código TermoFluids que está en continua evolución. Algunas de las mayores cuestiones que se plantean implican la generación de mallas adecuadas y asequibles, la implementación y validación de la condición de contorno periódica tridimensional y el acoplamiento de los diferentes dominios para el estudio de casos con diferentes comportamientos físicos, como desarrollo transitorio e inercia térmica. La turbulencia está presente en la mayoría de los flujos de ingeniería, y los intercambiadores de calor de evaporadores para refrigeración no son una excepción. La presencia de muchos tubos (que actúan como obstáculos para el fluido) colocados en diferentes configuraciones y el hecho de que el flujo también está confinado por aletas, crean características de flujo tridimensionales complejas que tienen generalmente régimen turbulento o en transición. Por lo tanto, se analiza la convección forzada turbulenta en una matriz de pines delimitados por paredes. simulando las grandes escalas de turbulencia y modelando las pequeñas (LES) con el fin de evaluar el desempeño de los tres modelos seleccionados, a saber WALE, QR y VMS. Los números de Reynolds establecidos para el estudio son 3000, 10000 y 30000. Algunos de los principales resultados que se incluyen son el coeficiente de presión alrededor los cilindros, el número de Nusselt promedio en las paredes y la vorticidad del flujo. La parte final de la tesis se dedica a estudiar el flujo tridimensional y los parámetros de transferencia de calor encontrados en un intercambiador de calor de tipo aleta-tubo utilizado para la refrigeración doméstica en equipos de 'no-escarcha'. Las implementaciones del código y el postproceso numéricos se validan en un caso muy similar para un intercambiador de calor con dos filas de tubos a bajos Reynolds para el cual se dispone de datos experimentales. El siguiente análisis que se presenta es una configuración típica para evaporadores 'no-escarcha' con paso de aleta doble (para el que se tiene muy poca información numérica en la literatura). Se considera el acoplamiento conjugado de la transferencia de calor convectiva entre fluido y sólido y conductiva dentro de la aleta. La influencia de algunos parámetros geométricos y de régimen de flujo se analizan con propósitos de diseño. En conclusión, las contribuciones generales de esta tesis junto con el código computacional ya existente, ha demostrado ser capaz de realizar con éxito simulaciones tridimensionales para predecir las características del flujo y los mecanismos responsables de la transferencia de calor en intercambiadores de calor de tipo aleta-tuboPostprint (published version