Source: http://docplayer.es/2143938-El-cultivo-protegido-bajo-invernadero-tiene-altos-costes-energeticos-derivados-de-la-ventilacion-forzada-en-cambio-la-ventilacion-natural-puede-ser.html
Timestamp: 2017-06-29 00:47:16+00:00

Document:
El cultivo protegido bajo invernadero tiene altos costes energéticos. derivados de la ventilación forzada. En cambio, la ventilación natural puede ser - PDF
El cultivo protegido bajo invernadero tiene altos costes energéticos. derivados de la ventilación forzada. En cambio, la ventilación natural puede ser
Download "El cultivo protegido bajo invernadero tiene altos costes energéticos. derivados de la ventilación forzada. En cambio, la ventilación natural puede ser"
Diego Cordero Quiroga
1 2 3 Resumen (Castellano) El cultivo protegido bajo invernadero tiene altos costes energéticos derivados de la ventilación forzada. En cambio, la ventilación natural puede ser una solución barata que reduzca el consumo energético. No obstante, un diseño que permita la ventilación natural del invernadero supone un reto debido a la complejidad del fenómeno físico. Frente a otros métodos de diseño, como los métodos experimentales o analíticos, en los últimos tiempos la dinámica de fluidos computacional (CFD por sus siglas en inglés) se ha convertido en la herramienta más utilizada para estudiar este tipo de fenómenos, gracias a su relativo bajo coste y a la rapidez en la obtención de resultados. No obstante, los modelos CFD deben refrendarse mediante validaciones realizadas a través de datos experimentales. Un análisis bibliográfico detallado del uso de la CFD aplicada a invernaderos muestra, en general, que los trabajos carecen de suficientes datos experimentales, seguramente debido al alto coste de los sensores para adquirirlos y la dificultad que conlleva el trabajo en condiciones de campo. Además, se observa que en las simulaciones CFD no existe una aproximación común al problema de la ventilación en invernaderos y sería interesante una sistematización. Por ello, la presente tesis, por un lado describe un dispositivo y unos métodos sencillos y baratos para obtener datos atmosféricos, y por otro, propone una visión crítica sobre la investigación realizada hasta el momento, con el fin de sistematizar la manera de generar modelos CFD aplicados a la ventilación natural de invernaderos. Finalmente, la tesis se complementa con un ejemplo sobre un caso práctico. i4 Para ello, en primer lugar, se realizó una revisión bibliográfica de las diferentes guías de buenas prácticas en diferentes campos de la tecnología, principalmente en edificación, para sistematizar y adaptar las recomendaciones para generar modelos CFD en invernaderos. En segundo lugar, se desarrolló un sistema de adquisición de datos sencillo, consistente en una red de sensores, que permite medir simultáneamente la velocidad y dirección de viento en 20 puntos. Este sistema de sensores fue calibrado y probado en campo satisfactoriamente obteniendo una precisión similar a los anemómetros comerciales con un precio 30 veces superior. En tercer lugar, se generaron 24 modelos del flujo de aire alrededor del invernadero, resultantes de la combinación de cuatro modelos de turbulencia (k-ε; RNG k-ε; k-ϖ y RSM); dos esquemas de discretización (primer orden y segundo orden) y tres velocidades de viento exterior (3; 3,5 y 4 m/s ) con el fin de analizar sus diferencias y demostrar sus ventajas e inconvenientes. Por este motivo, en cuarto lugar, se comprobó su capacidad de ajustarse a los datos de campo, validando los modelos con un análisis de regresión lineal sobre los datos experimentales. Con este estudió se reveló que los modelos SST k-ϖ y el RSM (segundo orden) son los que mejor representan el flujo de ventilación y se demostró que el modelo k- ε estándar (primer orden), el más utilizado en la bibliografía, no sólo ofrece resultados diferentes al resto de modelos, sino que su rendimiento es pobre para predecir el flujo de ventilación. ii5 Resúm (Valencià) El cultiu protegit baix hivernacle té alts costos energètics derivats de la ventilació forçada. En canvi, la ventilació natural pot ser una solució barata que reduïsca el consum energètic. No obstant això, un disseny que permeta la ventilació natural de l'hivernacle suposa un repte degut a la complexitat del fenomen físic. Enfront d'altres mètodes de disseny, com els mètodes experimentals o analítics, en els últims temps la dinàmica de fluids computacional (CFD per les seues sigles en anglès) s'ha convertit en la ferramenta més utilitzada per a estudiar este tipus de fenòmens, gràcies al seu relatiu baix cost i a la rapidesa en l'obtenció de resultats. No obstant això, els models CFD han de referendar-se per mitjà de validacions realitzades a través de dades experimentals. Una anàlisi bibliogràfica detallat de l'ús de la CFD aplicada a hivernacles mostra, en general, que els treballs no tenen suficients dades experimentals, segurament a causa de l'alt cost dels sensors per a adquirir-los i la dificultat que comporta el treball en condicions de camp. A més, s'observa que en les simulacions CFD no hi ha una aproximació comuna al problema de la ventilació en hivernacles i seria interessant una sistematització. Per això, la present tesi, per un costat descriu un dispositiu i uns mètodes senzills i barats per a obtindre dades atmosfèriques, i d'un altre, proposa una visió crítica sobre la investigació realitzada fins al moment, a fi de sistematitzar la manera de generar models CFD aplicats a la ventilació natural d'hivernacles. Finalment, la tesi es complementa amb un exemple sobre un cas pràctic. Per a això, en primer lloc, es va realitzar una revisió bibliogràfica de les iii6 diferents guies de bones pràctiques en diferents camps de la tecnologia, principalment en edificació, per a sistematitzar i adaptar les recomanacions per a generar models CFD en hivernacles. En segon lloc, es va desenvolupar un sistema d'adquisició de dades senzilles, consistent en una xarxa de sensors, que permet mesurar simultàniament la velocitat i direcció de vent en 20 punts. Este sistema de sensors va ser calibrat i provat en camp satisfactòriament obtenint una precisió semblant als anemòmetres comercials amb un preu 30 vegades inferior. En tercer lloc, es van generar 24 models del flux d'aire al voltant de l'hivernacle, resultants de la combinació de quatre models de turbulència (k-ε estàndard; RNG k-ε; SST k-ω i RSM); dos esquemes de càlcul (primer i segon ordre) i tres velocitats de vent exterior (3; 3,5 i 4 m/s) a fi d'analitzar les seues diferències i demostrar els seus avantatges i inconvenients. Per este motiu, en quart lloc, es va comprovar la seua capacitat d'ajustar-se a les dades de camp, validant els models amb una anàlisi de regressió lineal sobre les dades experimentals. Amb este estudi es va revelar que els models SST k-ω i el RSM (segon ordre) són els que millor representen el flux de ventilació i es va demostrar que el model k-ω estàndard (primer ordre), el més utilitzat en la bibliografia, no sols ofereix resultats diferents a la resta de models, sinó que el seu rendiment és pobre per a predir el flux de ventilació. iv7 Summary Protected crops, such greenhouses, have high energy costs due to forced ventilation. Hence, natural ventilation can be a reasonable solution to reduce energy consumption. However, a design that allows natural ventilation of the greenhouse is a challenge due to the complexity of the physical phenomena. Contrasting other design methods, such as experimental or analytical methods, in recent times the computational fluid dynamics (CFD) has become the most widely used tool to study the phenomena, thanks to its relative low cost and speed of obtaining results. However, CFD models must be endorsed by validations carried through experimental data. A detailed literature review of the use of CFD applied to greenhouses shows that the work lack sufficient experimental data, probably due to the high cost to acquire experimental data in field conditions. Moreover, it is observed that the CFD simulations are not performed in greenhouses with a systematic procedure. Therefore, this thesis on the one hand describes a device and a simple and inexpensive method to gather atmospheric data, and on the other hand, it proposes a critical view on the research conducted so far, in order to systematize the way to generate CFD models applied to natural ventilation in greenhouses. Finally, the thesis is complemented by an example of a case study. To do this, first, a review of good practices guides in different fields of technology was performed, mainly in environmental building engineering, with the intention to systematize and adapt these recommendations to generate CFD models from build engineering to greenhouses. Second, a simple data acquisition v8 system, consisting of a network of sensors that can simultaneously measure the wind speed and direction at 20 points was developed. This system was calibrated and tested at field conditions successfully with similar accuracy than commercial anemometers (with a price 30 times inferior). Third, 24 model cases were generated to analyse their differences and show their advantages and disadvantages. The cases are the result of the combination of four turbulent models (k-ε; RNG k-ε; SST k-ϖ y RSM); two discretization schemes (first and second order) and three external wind speeds (3, 3.5 and 4 m/s). Fourth, to analyse the performance of the models a validation procedure was perform: computational data and experimental data was compared by using linear regression analysis. This validation revealed that the SST k-ω and RSM models (second order ) are those that best represent the ventilation flow and demonstrated that the k-ε standard (first order), the most widely used in the literature, not only gives different results to other models, but their performance is poor for predicting ventilation flow. vi9 Agradecimentos Por todos los años que he formado parte del Centro de Agroingeniería del Instituto Valenciano de Investigaciones Agrarias y he estado vinculado al Departamento de Ingeniería Agraria de la Universidad Politécnica de Valencia, en los que he recibido el apoyo y ánimo de todos sus miembros para realizar esta tesis doctoral. En especial, en primer lugar, a Dr. Enrique Moltó y Dr. Guillermo Palau, directores de esta tesis doctoral, por haberme apoyado y confiado en mí para el progreso de la investigación en ventilación natural en invernaderos. Me siento afortunado de haber contado con dos personas que le ponen tanto cariño a esto de investigar. Asimismo, a todas las personas ya amigos que han formado parte de esta investigación directa o indirectamente: Rosana, Alejo, Ramón, Patricia, Cruz, Sento, Sergio, Nacho, Santiago, Talo, Cristian Ll., José, e como escribió Shakespeare: the last but not the least Iván. Igualmente, reconocer la buena disponibilidad y colaboración con la empresa Mediterránea de Cactus. De igual manera, a la Generalitat Valenciana por financiar este proyecto y así dotar de medios económicos para llevar a cabo esta investigación. vii10 Finalmente, a la familia no elegida pero tan acertada: a mi Padre, que empezaste esta tesis conmigo pero no tuviste tiempo para verla acabar; a mi Madre, por tu amor incondicional y tu ejemplo diario por el amor al buen trabajo; a mi hermano, por hacerme sentir siempre tan orgulloso. También, a esa familia sí elegida: Rubén y Bea, que me apoyasteis en Inglaterra el tiempo necesario para terminar de escribir este texto; Cristian P, por la paciencia que pusiste en nuestra amistad, entiéndeme, no puedo evitarlo, soy Valenciano; y finalmente a Javitxu, por tu cariño y perseverancia, en nuestro recuerdo quedará ya siempre la Albufereta, lado a lado, y dos tesis por escribir: la tuya y la mía. Y, a tantos otros que sé que me dejo fuera, en especial a los que encontré a partir de mayo de 2011 que consiguieron que mi vida fuera más difícil pero mucho más apasionante. Mejor, cambio el último agradecimiento: A mi Padre, Madre y hermanos. Por hacer que mis días siempre sean más claros, Gracias. viii11 Dedicado a mi Padre La libertad no es nada más que una oportunidad para ser mejor. (Albert Camus) ix12 Índice general 1. INTRODUCCIÓN AL PROBLEMA DE LA VENTILACIÓN NATURAL... 1 Importancia económica... 1 Tipos de ventilación DINÁMICA DE FLUIDOS APLICADA A LA VENTILACIÓN ECUACIONES FÍSICAS QUE DESCRIBEN LOS FLUJOS ATMOSFÉRICOS FLUJO LAMINAR Y FLUJO TURBULENTO EN LA ATMÓSFERA EFECTO DE LA CAPA LÍMITE ATMOSFÉRICA SOBRE EDIFICIOS E INVERNADERO FORMAS CONVENCIONALES DE ABORDAR LA VENTILACIÓN NATURAL: Métodos experimentales Modelos analíticos y empíricos Modelos multizona y de zona DINÁMICA COMPUTACIONAL DE FLUIDOS (CFD) GENERACIÓN DE LA MALLA RESOLUCIÓN DEL FLUJO PRINCIPAL: DNS, LES Y RANS MODELOS DE TURBULENCIA QUE SE EMPLEAN EN LA TÉCNICA RANS Aproximación de Boussinesq (k-ε standard, RNG y SST k-ω) Modelos de transporte de Reynolds (RSM) RESOLUCIÓN DEL FLUJO CERCANO A LA PARED (MODELO DE PARED) ESQUEMAS DE CÁLCULO Y RESIDUOS ANÁLISIS CRÍTICO DE GUÍAS DE BUENAS PRÁCTICAS EN CFD Elección de las variables objetivo Elección de las ecuaciones de aproximación para describir la física del flujo Elección de la representación geométrica de los obstáculos Elección del dominio computacional Elección de las condiciones de contorno x13 Simulación del perfil de viento Simulación del efecto del suelo Condiciones de contorno en la parte superior del dominio Condición de contorno lateral y salida del dominio Condiciones de contorno de la superficie de los edificios Elección de las condiciones iniciales Elección de la malla Elección de las aproximaciones numéricas y de los criterios de convergencia ANÁLISIS CRÍTICO DE LOS ESTUDIOS DE LA VENTILACIÓN NATURAL EN INVERNADEROS MEDIANTE CFD OBJETIVOS MATERIALES Y MÉTODOS EL SISTEMA AUTOMÁTICO DE ADQUISICIÓN SIMULTÁNEA DE DATOS Anemómetros sónicos Anemómetros mecánicos: sistema cazoleta-veleta Registrador de datos de los anemómetros mecánicos Calibración de los anemómetros mecánicos DESCRIPCIÓN DEL INVERNADERO MEDIDAS DE CAMPO SIMULACIONES Variables objetivo Ecuaciones de aproximación Representación geométrica de los obstáculos Dominio computacional Perfil del viento de entrada Modelo de suelo Condiciones de contorno Condiciones iniciales xi14 Malla Aproximaciones numéricas Criterios de convergencia Simulaciones ensayadas y validación RESULTADOS Y DISCUSIÓN CALIBRACIÓN DEL SISTEMA DE MEDIDA ANÁLISIS DE LA MALLA DEL INVERNADERO MEDIDAS EXPERIMENTALES DE VELOCIDAD Y DIRECCIÓN DEL VIENTO RESULTADOS DE LAS SIMULACIONES Y VALIDACIÓN DESCRIPCIÓN DEL FLUJO ALREDEDOR DEL INVERNADERO SUGERENCIA PARA UTILIZAR LOS RESULTADOS DE ESTA TESIS. CASO PRÁCTICO: ESTUDIO DEL FLUJO EN EL INTERIOR DE UN INVERNADERO CONCLUSIONES Y TRABAJO FUTURO CONCLUSIONES PROPUESTA DE TRABAJO FUTURO REFERENCIAS Trabajos y comunicaciones del autor ANEJO A DISCRETIZACIÓN Y LINEALIZACIÓN DE LAS ECUACIONES DE TRANSPORTE ANEJO B ALGORITMO DE RESOLUCIÓN DE LAS ECUACIONES DE NAVIER-STOKES ANEJO C RESULTADOS EXPERIMENTALES Y DE LA SIMULACIÓN ANEJO D RESULTADOS REGRESIÓN ANEJO E FOTOGRAFÍAS DEL INVERNADERO CASO ESTUDIO xii15 Índice de tablas Tabla 1 Principales cultivos y datos de superficie en España (Cotec, 2009) Tabla 2 Estructura de la capa superficial de la atmósfera Tabla 3 Constantes modelo standard k-ε, RNG y SST k-ω Tabla 4 Clasificación revisada de Davenport para obtener la longitud rugosa (zo) en función del tipo de suelo (Wieringa, 1992) Tabla 5 Revisión de las principales guías de buenas prácticas Tabla 6 Comparación de modelos CFD aplicados a invernaderos Tabla 7 Características técnicas de los anemómetros sónicos Tabla 8 Estación meteorológica Red SIAR (Modelo 1) Tabla 9 Valores z0 elegidos para la rugosidad del suelo alrededor del invernadero Tabla 10 Cuadro resumen de los parámetros que se utilizaron en las simulaciones Tabla 11 Malla del ambiente exterior del invernadero Tabla 12 Vectores velocidad alrededor del invernadero (velocidad exterior 4 m/s) Tabla 13 Resultados del análisis de regresión lineal (Caso práctico) Tabla 14 Caso práctico: Resumen de la simulación Tabla 15 Datos observados y predichos por los modelos. Módulo (m/s) Tabla 16 Datos observados y predichos por los modelos. Dirección (º) Tabla 17 Resumen de la calidad de la malla utilizada para el caso práctico Tabla 18 Velocidad estandarizada de las simulaciones y medidas experimentales xiii16 Índice de figuras Figura 1 Esquema del flujo alrededor de un objeto bidimensional con aristas cortantes: (a) líneas de corriente medias y (b) perfiles de la velocidad (Halitsky, 1968) Figura 2 Esquema del flujo alrededor de un objeto tridimensional con aristas cortantes (Halitsky, 1968) Figura 3 Ley de pared para superficies suaves (línea de puntos) y arenosas (línea continua) para diferentes ks+ (Blocken, et al., 2007) Figura 4 Representación esquemática del cálculo de (izquierda) primer orden y de (derecha) segundo orden Figura 5 Diagrama de la primera celda computacional. Perfiles de velocidad del modelo de pared y del perfil logarítmico. (Blocken, et al., 2007) Figura 6 Estación meteorológica PCE-FWS 20 utilizada para la toma de medidas en campo Figura 7 Ejemplo de salida de datos de los anemómetros mecánicos Figura 8 Ejemplo de salida de datos de los anemómetros sónicos 3D Figura 9 Conexión del sistema de medida desarrollado por el IVIA Figura 10 Calibración de los sensores de cazoleta y veleta frente al sensor sónico Figura 11 Emplazamiento invernadero. Catadau (Valencia) Figura 12 Detalle del invernadero y sus alrededores Figura 13 Puesta en campo de los sensores cazoleta, veleta y sónicos Figura 14 Planta del dominio del invernadero dividido mediante: (arriba) la técnica Split y (abajo) la técnica O-grid. Representación de la forma de las celdas en las esquinas: o-grid (a y b) y la técnica split (c y d) Figura 15 Comparación entre los módulos de velocidad los anemómetros de cazoleta y sónico. Ambos anemómetros miden con una precisión similar la velocidad del viento Figura 16 Comparación entre la dirección medida por el anemómetros de cazoleta y sónico. Ambos anemómetros miden con una precisión similar la dirección del viento Figura 17 Velocidad estandarizada en la fachada suroeste del invernadero para los diferentes modelos de turbulencia y datos experimentales (de arriba a abajo): a 3 m y 4 m de altura Figura 18 Velocidad estandarizada en la fachada Sureste (Barlovento) y Noroeste (Sotavento) para los diferentes modelos de turbulencia y experimentales xiv17 Figura 19 Representación de las líneas de corriente en un plano perpendicular al suelo contiene los sensores A. Modelo k-ε (arriba) y modelo RSM (abajo). Los cuadros rojos representan los puntos de medida para sensores A Figura 20 Regresión lineal de los datos experimentales frente a los datos de simulación para (de izquierda a derecha y de arriba abajo) k-ε estándar primer orden vs 3 m/s; k-ε primer orden vs 4 m/s; k-ε segundo orden 4 m/s; k, RSM segundo orden vs 4 m/s Figura 21 Detalle de la energía cinética turbulenta del modelo en la esquina Sureste-Noreste del invernadero k-ε (Izquierda) y RSM (Derecha) para velocidad exterior de 3 m/s Figura 22 Líneas de corriente en un plano longitudinal a 1,5 m de la fachada Noreste del invernadero predichas por el modelo RNG k-ε que contiene los puntos A Figura 23 Contornos de v-velocidad y líneas de corriente en dos planos horizontales (3m y 4m) para una velocidad exterior de 4 m/s y modelo RSM Figura 24 Contorno de la velocidad-v en la sección media del invernadero (a) y líneas de corriente en la misma sección (b) Figura 25 Velocidad estandarizada por la velocidad exterior para un viento (arriba) Sureste-Noroeste y (abajo) Noreste-Suroeste Figura 26 Algoritmo de resolución SIMPLE Figura 27 Representación de las velocidades experimentales estandarizadas por la velocidad exterior en todos los puntos medidos Figura 28 Velocidad estandarizada por la velocidad exterior en la zona de aproximación del viento exterior (sensores B) para los diferentes modelos de turbulencia y media de los datos experimentales (de arriba a abajo): a 3 m y 4 m de altura Figura 29 Contornos de la v-velocidad estandarizada frente a la velocidad exterior para: (de arriba a abajo) k-ε primer orden y RNG, SST, y RSM segundo orden; y (de izquierda a derecha) a 3 m y a 4,5 m de altura Figura 30 Líneas de corriente en un plano paralelo a la fachada Noreste que contiene a los sensores A para (de arriba a abajo): k-ε primer orden; RNG k-ε, SST k-ω y RSM de segundo orden Figura 31 Contornos de v-velocidad en al invernadero (de arriba abajo): centro del vano Nordeste; entre los xv18 dos últimos vanos del Nordeste y a 1,5 metros de la fachada Nordeste; y de (izquierda a derecha): k-ε primer orden; RNG k-ε, SST k-ω y RSM de segundo orden xvi19 Nomenclatura Nomenclatura B Constante para el cáculo del perfil logarítmicos Cs Parámetro de rugosidad de ANSYS Fluent. E Constante empírica de valor aproximado 9,793 f Parámetro de coriolis g Aceleración de la gravedad. Se toma el valor: 9,81(m/s 2 ) H Altura característica de la simulación (edificio considerado) (m) J Operador Jacobiano k Energía cinética turbulenta (m 2 /s 2 ) k s Altura de rugosidad arenosa equivalente (m) L Longitud (m) p punto central de la celda adyacente a la pared Re Número adimensional de Reynolds R 2 Coeficiente de determinación R s Φ Residuo escalado de una variable R Φ Residuo no escalado de una variable T Tiempo (s) xvii20 Nomenclatura U Velocidad media (m/s) u' Velocidad instantánea: Fluctuación con respecto al tiempo de la velocidad de una partícula con respecto a su valor medio (m/s) u;v;w Componentes de la velocidad en la dirección x;y;z (m/s) u p Velocidad en el punto P (m/s) u ref Velocidad de referencia para el cálculo del perfil logarítmico de viento (m/s) u τ 2 Velocidad de fricción para el cálculo de la función de pared (m/s) u Velocidad de fricción de pared (m/s) u ABL Velocidad de fricción para el cálculo del perfil de viento (m/s) y p Distancia de la pared al primer nodo computacional (m) Z Componente dirección vertical (m) z 0 Parámetro experimental llamado longitud rugosa xviii21 Nomenclatura Modelos k-ε G k Representa la generación de energía cinética turbulenta debido a los gradientes de la velocidad media σ k ;σ ε Número turbulentos de Prandtl para k y ε. S k ;S ε Parámetros definidos por el usuario que representa el término fuente para k y ε C μ Constante para la estimación de la viscosidad turbulenta α k ;α ε Números de Prandtl efectivos inversos para k y ε Modelos k-ω G k Generación de energía cinética turbulenta debido a los gradientes de la velocidad media para el modelo SST k-ω G ω Término que representa la generación de ω Γ k ;Γ ω Difusividad efectiva de k y ω en el modelo SST k-ω Y k ;Y ω Disipación de k y ω debido a la turbulencia en el modelo SST k-ω D ω Término que representan el término de difusión cruzada S k ;S ω Parámetros definidos por el usuario que representa el término fuente para k y ω xix22 Nomenclatura Modelos RSM C ij Término que representa el parámetro de la convección D L,ij Término que representa la difusión molecular P ij Término que representa la producción del esfuerzo F ij Término que representa la producción del sistema rotacional D T,ij Término que representa la difusión turbulenta Φ ij Término que representa la deformación de presión ε ij Término que representa la disipación xx23 Nomenclatura Símbolos griegos Φ Valor de una variable cualquiera Φ f Φ interpolada en una cara de una celda Φ nb Valor de una variable cualquiera en el centro de una celda vecina Φ p Valor de una variable cualquiera en el centro de una celda δ ij Función delta de Kronecker ε ij3 Tensor unitario μ Viscosidad dinámica (kg/s m) μ t Viscosidad turbulenta (kg/s m) γ Peso específico (kg/m 3 ) ε Ratio de disipación de turbulencia (m 2 /s 3 ) η Cantidad de una propiedad por unidad de masa κ Constante de Karman ν Viscosidad cinemática (m 2 /s) ρ Densidad kg/m 3 τ Esfuerzo cortante debido a la viscosidad (N/m) ω Disipación específica (s 1 ) xxi24 Nomenclatura Acrónimos ABL Atmospheric boundary layer (Capa Límite atmosférica) AIJ Architectural Institute of Japan (Instituto de arquitectura de Japón) ASCII American Standard Code for Information Interchange (Código internacional estandarizado para intercambio de información) BPG Best practice guidelines (Guía de buenas prácticas) CFD Computational fluid dynamics (Dinámica de fluidos computacional) COST European Cooperation in Science and Technology (Cooperación europea en ciencia y tecnología) DNS Direct Numeric Simulation (Simulación numérica directa) ECORA Evaluation of Computational Fluid Dynamic Methods (Evaluación de métodos para dinámica computacional de fluidos) ERCOFTAC European Research Community on Flow, Turbulence And Combustion (Comunidad de investigación europea en flujo, turbulencia y combustión) FVM Finite volume method (Método de volúmenes finitos) GBP Guías de buenas prácticas HVAC Heating, ventilation and air conditioning (Calefacción, ventilación y aire acondicionado) xxii25 Nomenclatura IP Ingress Protection Rating (Grado de protección) LES Large eddy simulation (Simulación de grandes vórtices) QNET-CFD Quality Net (web que agrupa guías de buenas prácticas) RANS Reynolds average Navier-Stokes (Promediado de Reynolds de las ecuaciones de Navier-Stokes) RNG Renormalization group (Grupo de renormalización) SCM Second Clousre Model (Modelo de cierre de segundo orden) SIAR Sistema de información agroclimática para el regadío SST Shear-Stress Transport (Transporte de esfuerzo de cortadura) VDI Verein Deutscher Ingenieure (Asociación de Ingenieros Alemanes) xxiii26 Introducción al problema de la ventilación 1. Introducción al problema de la ventilación natural El sistema de producción agrícola bajo invernadero asegura el aprovisionamiento continuo de frutas y verduras frescas durante todo el año, en todas las regiones del mundo y bajo todos los climas. Esta independencia geográfica y climática, junto a la posibilidad que ofrecen de gestión de recursos escasos como el agua o suelo, explica el incremento espectacular del uso del invernadero en las últimas décadas, lo que implica una revolución en el volumen de la producción agrícola. Importancia económica La incorporación de este tipo de tecnología en España no se produjo hasta los años 60, es decir, llegó con retraso respecto a países del norte de Europa. Sin embargo, actualmente España es el tercer país en superficie dedicada a invernaderos, después de China y Corea del Sur (International Greenhouse Vegetable Production, 2011). Según el Ministerio de Agricultura, Alimentación y Medio Ambiente (2008) las principales áreas de cultivo se concentran en la cuenca mediterránea: Almería ( hectáreas cultivadas), Murcia (9.000 hectáreas), Canarias (7.500 hectáreas), y Comunidad Valenciana (1.800 hectáreas). En la Tabla 1 se muestra además los principales cultivos y datos de superficie en nuestro país. 127 Introducción al problema de la ventilación Tabla 1 Principales cultivos y datos de superficie en España (FAO, 2011). Cultivo Producción (hectáreas) (toneladas) Berenjenas Calabacín/cal Pimiento Judías V Melones Pepinos Sandías Tomate Junto a esta importancia económica, la horticultura de invernadero tiene gran trascendencia social por el volumen de mano de obra que precisa, empleos entre directos e indirectos, en la fase de producción y comercialización. En nuestro país, aproximadamente la mitad de la superficie bajo invernadero pertenece a pequeñas explotaciones familiares en las que este cultivo es la principal entrada económica. Ministerio de ministerio de Agricultura, Alimentación y Medio Ambiente (2008). Pero este sistema de producción agrícola no está exento de problemas relacionados con el medio ambiente, las condiciones de trabajo, las normas sanitarias, o el aumento del precio de la energía. Actualmente, la crisis energética y económica, unida a la crisis medioambiental, fuerza a las personas y a los gobiernos de todo el mundo a buscar alternativas más sustentables que las formas de producción tradicionales. Así pues, la agricultura bajo invernadero requiere mejorar la ventilación para proporcionar niveles óptimos de temperatura y humedad sobre el cultivo de manera mucho más eficaz, de modo que se reduzcan los actuales consumos de energía. Existen varios estudios en los que se valora el impacto económico del consumo de energía entre los que 228 Introducción al problema de la ventilación destacamos el trabajo de Torrellas et al. (2012). Tipos de ventilación La ventilación es el proceso de renovación del aire en un espacio cerrado. Para ello el aire debe ser continuamente extraído y reemplazado por aire fresco. La ventilación de los invernaderos se usa principalmente para disipar el excesivo calor, mejorar el intercambio de dióxido de carbono y oxígeno de las plantas y mantener niveles aceptables de humedad. En este sentido, la ventilación es un proceso fundamental que influye sobre el clima en el interior del invernadero y afecta fuertemente el rendimiento de los cultivos. En tecnología de invernaderos se utilizan dos sistemas de ventilación: la ventilación forzada y la ventilación natural. La ventilación forzada consiste en inyectar, o extraer, aire del invernadero generalmente a través de sistemas de refrigeración por agua, utilizando ventiladores eléctricos; sistema conocido en la industria como cooling. De esta forma, se puede controlar muchos parámetros climáticos, pero adoptarlo supone una inversión elevada y un elevado consumo energético. Por este motivo se hace imprescindible investigar sobre métodos de ventilación natural en invernaderos (Norton, et al., 2007); (Bournet & Boulard, 2010). La ventilación natural como recurso renovable no consume energía y, por tanto, interesa estudiarla con el fin de implementarla correctamente en el diseño de todo tipo de instalaciones (Tecle, et al., 2012). En los invernaderos normalmente se consigue la ventilación natural mediante aperturas de ventanas permanentes o temporales en el techo y las paredes. Las fuerzas responsables 329 Introducción al problema de la ventilación de la ventilación natural en los invernaderos se deben a dos fenómenos físicos: la flotación, o efecto chimenea y la acción del viento exterior; efectos conocidos en la literatura científica como stack effect y wind driven respectivamente. Las fuerzas de flotación están asociadas a gradientes de densidad producidos por las diferencias de temperatura entre el interior y el exterior del invernadero, que producen un desplazamiento vertical del aire. Tradicionalmente la mayoría de sistemas de ventilación natural, como en el caso de edificación, están diseñados sólo para esta fuerza. En cambio, en tecnología de invernaderos, la flotación se considera despreciable para velocidades medias del viento exterior superiores a 1,5 m/s (Kittas, et al., 1995). Los fenómenos de ventilación producidos por la acción del aire exterior son complejos y se atribuyen a la interacción de dos flujos con diferentes escalas y propiedades: el viento exterior y el interior. El viento exterior incide sobre la fachada a barlovento y crea una diferencia de presiones entre ésta y la fachada a sotavento. Esta diferencia de presiones genera un gradiente que dirige el flujo interior a través del recinto en la dirección de dicho gradiente. El campo de presiones que conduce el flujo el interior depende la estructura y la fluctuación del flujo exterior. 430 Dinámica de fluidos aplicada a la ventilación 2. Dinámica de fluidos aplicada a la ventilación La dinámica de fluidos es la rama de la mecánica que estudia el movimiento de los fluidos así como las fuerzas que provocan este movimiento. Un fluido, desde este punto de vista, es cualquier sustancia cuya estructura molecular no ofrece resistencia a los esfuerzos cortantes. A pesar de que en esta definición entran tanto los gases como los líquidos, en la presente tesis cuando hablemos de fluidos nos referiremos estrictamente a los gases, más concretamente al aire atmosférico. En este epígrafe en primer lugar, se hace una revisión de la física que subyace en la ventilación natural a través de sus ecuaciones de gobierno. En segundo lugar, se describen los fenómenos implicados haciendo hincapié en el flujo turbulento. En tercer lugar, se describe el efecto de la capa límite atmosférica que es la que interviene en los fenómenos de ventilación sobre edificios e invernaderos. Finalmente, se describen los métodos tradicionales para el estudio de los procesos de ventilación en edificios e invernaderos Ecuaciones físicas que describen los flujos atmosféricos El viento se produce como resultado de la transformación de la energía potencial en energía cinética debido a las fuerzas de presión y estas a su vez son generadas por las diferencias de temperatura. Los vientos de interés para la ventilación natural actúan dentro la capa límite superficial de la atmósfera, donde 531 Dinámica de fluidos aplicada a la ventilación los efectos producidos por la rotación de la tierra son despreciables. El tratamiento matemático de estos está basado en las ecuaciones fundamentales del movimiento. La dinámica de fluidos considera los fluidos como medios continuos, es decir no considera las posibles discontinuidades existentes en el nivel microscópico. Estas ecuaciones no son más que expresiones matemáticas deterministas de los principios fundamentales de conservación: de la masa, del momento y de la energía. Aunque esta última, la energía, no será considerada en la presente tesis al no considerarse las diferencias de temperaturas del aire que afectan a las fuerzas de flotación. El principio de conservación de la masa, utilizando el criterio del sumatorio de Einstein, para un volumen finito de fluido, nos lleva a la ecuación de la continuidad para flujos incompresibles en coordenadas cartesianas, ecuación (2.1), o lo que es lo mismo, que la ecuación de la divergencia de la velocidad es cero. Esta ecuación impone la restricción de que la tasa de expansión o compresión de un fluido debe de ser siempre cero para cualquier volumen finito del flujo. Esta hipótesis de flujo incomprensible sólo es correcta para Mach inferior a 0.3. Este número adimensional para flujo atmosférico corresponde aproximadamente a velocidades inferiores a 100 m/s, por lo que los flujos de ventilación en invernaderos se tratan siempre como fluidos no compresibles. u i x i = 0 (2.1) Que representa la variación del vector velocidad, u i, en cada una de las 632 Dinámica de fluidos aplicada a la ventilación direcciones, x i, Además, manteniendo los mismos criterios, la segunda ley de Newton lleva a la llamada ecuación de Navier-Stokes que describe el balance de fuerzas en la atmósfera y que sin considerar las fuerzas centrífugas toma la forma de la ecuación (2.2). u i t = u u i j δ x i3 g + fε ij3 u j 1 p + 1 j ρ x i ρ ( τ ij ) (2.2) x j τ ij = μ ( u i x j + u j x i ) 2 3 μ ( u k x k ) δ ij (2.3) En la ecuación (2.2), u i es la componente horizontal del viento en la dirección i; xi son ejes ortonormales de referencia; δ ij la función delta de Kronecker; g la aceleración de la gravedad que en esta tesis toma el valor de 9,81 m/s; f es el parámetro de Coriolis; ε ij3 es un tensor unitario que para ijk=123, 231 o 312 es +1, para ijk=321, 213 o 132 es -1 y para cualquier otra combinación es 0; ρ es la densidad del aire; p es la presión; τ ij es el tensor cinemático de esfuerzos cortantes (2.3). En la ecuación (2.3), μ es la viscosidad cinemática del aire. De izquierda a derecha, los términos de la ecuación (2.2) son: la tendencia (o término inercial), la advección, la fuerza del gradiente de presión y el esfuerzo cortante molecular (o término viscoso). Las ecuaciones de Navier-Stokes son deterministas, es decir, para unas condiciones iniciales ofrecen una única solución. No obstante, no tienen fácil resolución para flujos complejos como veremos más adelante y su resultado depende fuertemente de las condiciones 733 Dinámica de fluidos aplicada a la ventilación iniciales. En la literatura al conjunto de las tres ecuaciones, (2.1) (2.2) y (2.3), se le denomina ecuaciones de Navier-Stokes (Anderson, 1995) Flujo laminar y flujo turbulento en la atmósfera Tradicionalmente se distinguen dos tipos de flujos: los flujos laminares y los flujo turbulentos. En un flujo laminar el movimiento de las partículas en el fluido es ordenado, estratificado y suave. Se considera que el fluido se mueve en láminas paralelas sin entremezclarse y cada partícula de fluido sigue una trayectoria suave, llamada línea de corriente. Sin embargo, la mayoría de flujos que se encuentran en la naturaleza son turbulentos. La turbulencia no es una característica propia de los fluidos sino del flujo. Resulta difícil dar una definición precisa de este fenómeno y generalmente nos conformamos con realizar una lista con alguna de las características de la turbulencia (Tennekes & Lumely, 1972). Los flujos turbulentos son siempre flujos rotacionales, con rápidas variaciones temporales de sus propiedades en todas las dimensiones del espacio, donde se generan gran cantidad de estructuras generalmente circulares llamadas vórtices. El tamaño de los mayores vórtices coincide con las escalas del flujo principal e interaccionan con éste disminuyendo su energía cinética. Los vórtices son inestables y, por los efectos de la cortadura, o por interacción entre ellos, se dividen en vórtices más pequeños (vortex stretching). El proceso de rotura se produce en cascada: los vórtices más grandes van traspasando su energía cinética a vórtices más pequeños sucesivamente, hasta que son tan pequeños que su energía cinética no es lo suficientemente grande como para que la inestabilidad persista, siendo las 834 Dinámica de fluidos aplicada a la ventilación fuerzas viscosas preponderantes. En estos vórtices más pequeños, la energía cinética contenida se transforma en energía térmica por efecto de la viscosidad. Al proceso completo se le denomina cascada de energía. Este fenómeno es el responsable de una característica importante: los flujos turbulentos son siempre disipativos y, por tanto, la turbulencia requiere un aporte continuo de energía para mantenerse. Finalmente hay que resaltar que los flujos turbulentos son difusivos, produciendo una rápida mezcla de todas sus propiedades y los fenómenos de transporte de masa, cantidad de movimiento y energía se ven notablemente amplificados por efecto de la turbulencia. Esta rápida mezcla de propiedades durante la ventilación es la de mayor interés en el estudio de la ventilación natural en invernaderos. El número de Reynolds (Re) (Reynolds, 1883) permite clasificar el régimen de un fluido, mediante la ecuación (2.4). Re ( ) = vl ν (2.4) Donde v (ms 1 ) es la velocidad media del fluido, L (m) es una longitud característica, que generalmente corresponde con el tamaño del vórtice de mayor tamaño para la escala considerada y ν (m 2 s 1 ) es la viscosidad cinemática. El número de Reynolds establece la relación entre los términos inerciales y los términos viscosos en la ecuación de Navier-Stokes. Los flujos turbulentos ocurren siempre para números de Reynolds altos y, a menudo, se originan como una inestabilidad de los flujos laminares al aumentar el número de Reynolds, es decir, cuando las fuerzas inerciales comienzan a preponderar. 935 Dinámica de fluidos aplicada a la ventilación En la mayoría de los casos, debido a la falta de herramientas matemáticas, resulta complicado realizar un análisis determinista de la turbulencia y es más útil describirla mediante sus propiedades estadísticas; para ello se utiliza la descomposición de Reynolds, que consiste en descomponer una propiedad del flujo en un valor promedio y una fluctuación con respecto al promedio. 1036 Dinámica de fluidos aplicada a la ventilación 2.3. Efecto de la capa límite atmosférica sobre edificios e invernadero La capa superficial de la atmósfera se puede dividir en cuatro capas en función de su régimen: en contacto con el suelo y en los primeros milímetros presenta flujo laminar y a medida que va aumentando la altura el régimen deviene turbulento con una capa donde los dos fenómenos son igual de importantes (Tabla 2). Tabla 2 Estructura de la capa superficial de la atmósfera Altura (m) Nombre Flujo ~ ~ Capa turbulenta Subcapa dinámica Turbulento ~ Subcapa viscosa Laminar/turbulento ~ Capa límite laminar Laminar Al incidir el viento sobre un edificio se produce un cambio en el perfil medio de velocidad del viento y de las líneas de corriente del viento debido a una separación del flujo (Figura 1 y Figura 2). La primera causa de la separación del flujo es su pérdida de energía cinética al encontrarse con el obstáculo, lo que produce un aumento en la energía potencial que se convierte en un aumento de la presión lo largo de la dirección del flujo hasta que éste ya no puede moverse en contra de este incremento de presión. Este gradiente adverso de presión modifica la dirección del flujo y se produce una separación de éste con las capas límite: tanto en la fachada de barlovento como en las fachadas laterales del edificio. Esta separación, junto con la forma del edificio, provoca una estela con una zona de recirculación con dos sistemas de vórtices. Finalmente el flujo se 1137 Dinámica de fluidos aplicada a la ventilación vuelve a unir a sotavento después del edificio. En la estela, el flujo turbulento se caracteriza por tener bajas velocidades pero altas tensiones cortantes y alta intensidad turbulenta. En la parte de la estela más cercana al edificio, el flujo medio y las estructuras turbulentas se ven fuertemente afectadas por dos fenómenos: los esfuerzos cortantes debidos a la separación del flujo y por los dos sistemas de vórtices desprendidos de los ejes del edificio. El primer sistema de vórtices, llamado de herradura horizontal se produce cerca del suelo, a barlovento, y se desvanece a sotavento como pares de remolinos girando en sentido opuesto alrededor del edificio. El segundo sistema consiste en una serie de vórtices periódicos alternados, llamados de Von Karman, que se producen a sotavento. 1238 Dinámica de fluidos aplicada a la ventilación Figura 1 Esquema del flujo alrededor de un objeto bidimensional con aristas cortantes: (a) líneas de corriente medias y (b) perfiles de la velocidad (Halitsky, 1968) Figura 2 Esquema del flujo alrededor de un objeto tridimensional con aristas cortantes (Halitsky, 1968). 1339 Formas tradicionales de abordar la ventilación 2.4. Formas convencionales de abordar la ventilación natural: Tradicionalmente el estudio de la ventilación natural se ha abordado de tres maneras o mediante una combinación de ellas: experimentalmente, mediante modelos analíticos y/o semi-empíricos y mediante los llamados modelos multizona y de zona. A continuación haremos una pequeña revisión de ellas Métodos experimentales La gran mayoría de estudios experimentales del flujo en ventilación se utilizan para obtener información con la que validar modelos numéricos, analíticos, empíricos o semi-empíricos. Estudiar la ventilación natural directamente mediante sensores es difícil, debido al elevado número de sensores que se necesitan para lograr una descripción razonable de las variables de interés tanto en edificios como en invernaderos (Shen, et al., 2012). Aunque a veces las variables del flujo se miden in-situ, a menudo se utilizan modelos a escala en túneles de viento o de agua, lo que permite realizar medidas más reproducibles con un menor coste. La similitud entre el flujo real y el flujo escalado se consigue manteniendo en el modelo a escala los números adimensionales que describen el flujo real, entre otros, los números de Reynolds, Froude, Euler, Mach, Grashof, Arquímedes y Rayleigh. Los túneles de viento permiten simular las condiciones de contorno de la capa atmosférica y medir su interacción con los edificios agrícolas. Los túneles de agua se han utilizado para estudiar las fuerzas de flotación en el interior de los invernaderos producidas por las diferencias de temperatura (Bartzanas, et al., 1440 Formas tradicionales de abordar la ventilación 2007). La mayoría de estudios sobre ventilación natural de edificios se han llevado a cabo en modelos a escala reducida. Sin embargo, los estudios de ventilación en invernaderos se han realizado mayoritariamente in situ (Bournet & Boulard, 2010), pues resulta complicado mantener los números adimensionales (Reichrath & Davies, 2002). Los experimentos de campo incluyen generalmente la medida de la velocidad del flujo, la temperatura, la humedad o la radiación. A veces se usan técnicas indirectas como la liberación de un gas trazador como, por ejemplo, en los trabajos llevados a cabo por Pérez-Parra J. et al. (2004) Modelos analíticos y empíricos Para estudiar la ventilación natural, los modelos analíticos parten de las ecuaciones fundamentales de la dinámica de fluidos y la termodinámica (conservación de la masa, el momento o la cantidad de movimiento, la energía y la concentración de especies químicas), simplificándolas convenientemente para llegar a ecuaciones relativamente sencillas. Los modelos empíricos tienen la misma base teórica, pero en su desarrollo se introducen coeficientes que provienen de datos experimentales o de simulaciones. En general, son modelos poco generalizables, aunque la metodología para su construcción a menudo puede serlo. No obstante, no suelen servir para representar fenómenos de ventilación en los que entran en juego muchas variables físicas (Chen, 2009). Existe un modelo analítico-empírico ampliamente utilizado para invernaderos, basado en la ecuación de Bernouilli. Parte de las suposiciones de que el flujo es permanente, irrotacional e hidrostático, que el fluido no es viscoso, 1541 Formas tradicionales de abordar la ventilación y de que la única acción exterior sobre el mismo es la fuerza de la gravedad. Estas suposiciones violan el principio de la conservación de la energía cinética y por ello en los cálculos se emplea un coeficiente para crearla o destruirla artificialmente. Los valores de este coeficiente dependen de la geometría y del régimen del flujo estudiado y normalmente se ajustan con datos experimentales o simulaciones de ordenador. Esta aproximación limita la precisión de la soluciones de las ecuaciones. Algunos autores como Roy et al. (2002) o Axley y Chung (2005) mejoraron el modelo añadiendo explícitamente ecuaciones que incluyen la conservación de la energía cinética, pero a pesar de estas mejoras el método no permite determinar la distribución espacial de las variables del flujo (Norton, et al., 2007) Modelos multizona y de zona Estos modelos se utilizan principalmente para predecir variables relacionadas con el número de renovaciones horarias (número de veces que el volumen del aire interior es renovado por hora) y la distribución del flujo en edificios. Para crear estos modelos se resuelven las ecuaciones de conservación de la energía, de la masa y de las especies químicas sin tener en cuenta la ecuación de la conservación del momento. Por tanto están restringidos a usos en los que el flujo sea cuasi-estacionario o esté en calma, es decir, en los que el efecto de la conservación del momento pueda ser despreciado. Los modelos multizona asumen la uniformidad de la temperatura y de la concentración de especies químicas. Los modelos de zona dividen el espacio en un número limitado de celdas, resolviendo las ecuaciones en cada una de estas celdas y de 1642 Formas tradicionales de abordar la ventilación este modo permitiendo que dichas variables no se distribuyan uniformemente. No se ha encontrado ningún estudio en el que se utilicen estos modelos para estudiar la ventilación en invernaderos. 1743 Dinámica computacional de fluidos (CFD) 3. Dinámica computacional de fluidos (CFD) Se puede considerar que la CFD surge de los trabajos de Richardson (1910), pero debido a la falta de potencia de cálculo de los ordenadores en la primera mitad del siglo XX y a la necesidad de desarrollo de algunas herramientas matemáticas, no se empezó a utilizar hasta la década de los 70, sobre todo para estudiar flujos simples en aplicaciones aeroespaciales e industriales en las que la predicción del comportamiento del flujo era importante y era complicado realizar medidas experimentales (por ejemplo, la entrada de naves espaciales en la atmósfera). En la década de los 90 se generalizó su uso, y actualmente se emplea en numerosos tipos de aplicaciones en los que intervienen flujos, transferencia de calor, reacciones químicas, dispersión de contaminantes, flujos bifásicos, cambios de fase, transferencia de masa y esfuerzos al interactuar fluidos con sólidos. La dificultad de resolver de forma exacta las ecuaciones diferenciales que describen los flujos lleva a utilizar métodos numéricos para transformarlas en sistemas algebraicos simples. Existe una gran variedad de métodos numéricos, sin embargo, los más utilizados son tres: de diferencias finitas, de elementos finitos y de volúmenes finitos. El método de las diferencias finitas consiste en sustituir las derivadas parciales por diferencias en nudos infinitamente próximos. Así pues, utilizando el sentido físico de la pendiente de una curva en un punto se puede sustituir la derivada por la diferencia entre el valor de la función en puntos (nudos) contiguos, esto dará lugar a diferentes a restas que resultan en las diferencias 1844 Dinámica computacional de fluidos (CFD) finitas. La gran desventaja de este método es que está restringido a mallas simples, ya que tiene dificultades para conservar el momento, la energía y la masa. El método de los elementos finitos consiste en calcular funciones sencillas, generalmente polinomios, cuyo valor en cada punto del dominio se aproxime al valor de la función estudiada en ese punto. Se trata de un método cuyos resultados dependen del tamaño de la celda y es excelente para problemas en el que domine la difusión (flujos viscosos). En cambio, es costoso computacionalmente y no se aconseja para flujos turbulentos. Finalmente, la técnica de los volúmenes finitos consiste en sustituir las ecuaciones diferenciales parciales por ecuaciones algebraicas y aplicarlas en un volumen discreto del espacio (celda). De forma similar a los métodos de elementos finitos y de diferencias finitas, los valores se calcularan en puntos discretos de un domino pero, a diferencia de los anteriores, se discretiza la forma integral de las ecuaciones y no la forma diferencial. Seguidamente, las leyes de conservación se aplican para las variables de flujo en puntos discretos de esas celdas llamados nodos. Se trata de la técnica utilizada en la presente tesis y es la más utilizada en los programas comerciales de CFD. La técnica de los volúmenes finitos se aplica en tres pasos básicos: en primer lugar, se malla el espacio ocupado por el fluido, es decir, se divide el dominio en pequeños volúmenes. En segundo lugar, se integran las ecuaciones de gobierno de la mecánica de fluidos para cada una de las celdas construyendo ecuaciones algebraicas para las variables discretas dependientes (incógnitas) 1945 Dinámica computacional de fluidos (CFD) como la velocidad, presión o temperatura. Finalmente, las ecuaciones algebraicas han de ser linealizadas y resueltas de una forma ordenada para actualizar los valores de las variables independientes. En este capítulo se explicarán las bases de la CFD. En primer lugar se abordará la discretización del dominio computacional a través de la malla. En segundo lugar, se abordarán las diferentes formas de resolver el flujo principal. En tercer lugar, se explicarán diferentes modelos de turbulencia que existen para modelar las ecuaciones de Navier-Stokes. En cuarto lugar, se explicará cómo se resuelve el flujo en las zonas cercanas a la pared. En quinto lugar, se explicarán los esquemas de cálculo disponibles y cómo se obtienen los residuos de las simulaciones para determinar los criterios de convergencia. En sexto lugar, se hará una revisión crítica de las guías de buenas prácticas encontradas en los diferentes campos de la ingeniería. Finalmente, se realizará un análisis crítico de la bibliografía disponible de validación de modelos CFD en invernaderos Generación de la malla La forma de generar la malla influye sobre los resultados de la simulación, por lo que el mallado se ha convertido en una disciplina propia dentro de la CFD (Baker, 2005). Es imprescindible que los resultados del modelo sean independientes de la malla utilizada y sólo dependan de las condiciones de contorno y de los parámetros físicos del problema. También hay que tener en cuenta el tamaño de la malla: a mayor número de celdas en el mismo volumen (mayor resolución), se resuelven con más precisión las ecuaciones de Navier- Stokes. Sin embargo, a mayor resolución, se necesita mayor capacidad de 2046 Dinámica computacional de fluidos (CFD) cálculo; hasta el punto que la simulaciones puedan ser inabordables. Así pues, es necesario un compromiso entre la precisión de los resultados y la resolución de la malla. En la literatura científica se emplean dos métodos para discretizar el espacio en tres dimensiones: las mallas no estructuradas y las estructuradas. En las mallas no estructuradas las celdas se generan automáticamente, conectando arbitrariamente puntos (denominados vértices) mediante líneas con la condición de que éstas no intersecten. Estas líneas definen los vértices de las celdas que, en el caso tridimensional, forman poliedros, generalmente tetraedros. En cambio, las mallas estructuradas se elaboradas manualmente. La idea básica, formulada por Lee, Huang, Yu y Rubbert en (1980), es dividir el dominio físico en varios bloques y posteriormente definir las celdas de forma separada en cada uno de esos bloques generalmente con forma de hexaedro. En general, las mallas no estructuradas se adaptan mejor a dominios físicos complejos y su generación consume mucho menos tiempo. En cambio, las mallas estructuradas resuelven mejor las ecuaciones del flujo, sobre todo cuando dominan los fenómenos viscosos, ya que se puede incrementar la resolución del mallado en las capas límites. La calidad de las mallas se valora por su resolución y a través de factores de forma. A menudo se emplea el factor de simetría, o skewness, que evalúa la forma de las caras o celdas frente a una cara o celda equilátera o equiángula considerada ideal. En este trabajo se ha utilizado el parámetro que se considera más interesante en la literatura, denominado determinante 3x3x3. Éste se 2147 Dinámica computacional de fluidos (CFD) define como el ratio entre el menor determinante dividido por el mayor determinante de la matriz Jacobiana para cada nodo de la malla. El determinante Jacobiano es el determinante del operador Jacobiano (J) que conecta las derivadas de las coordenadas naturales (r,s,t) con las derivadas de las coordenadas locales (x,y,z): J = ((dx/dr dy/dr dz/dr) (dx/ds dy/ds dz/ ds) (dx/dt dy/dt dz/dt)). Así pues, un valor de 1 indica una malla perfectamente regular; un determinante cercano a 0 indica elementos degenerados en cualquiera de las direcciones; y valores negativo indican elementos invertidos. Una malla con calidad insuficiente es inaceptable, porque los métodos numéricos asumen que las celdas son relativamente equiláteras o equiángulas, es decir, con un determinante 3x3x3 cercano a 1. En la actualidad no existe un procedimiento comúnmente aceptado para el diseño de las mallas, hasta el punto de que algunos autores se preguntan si la generación de mallas es un arte o una ciencia (Baker, 2005). No obstante, generalmente se considera que las mallas hexaédricas generan soluciones más precisas que las mallas tetraédricas, pero no existen trabajos que lo demuestren. Sin embargo, existen intentos de sistematizar el proceso de generación de las mallas a través de determinadas guías de buenas prácticas, como veremos más adelante. 2248 Dinámica computacional de fluidos (CFD) 3.2. Resolución del flujo principal: DNS, LES y RANS En CFD las ecuaciones de Navier-Stokes para el tratamiento de la turbulencia se abordan de tres maneras diferentes: la técnica DNS (Direct Numerical Simulation), la técnica RANS (Reynolds-Averaged Navier-Stokes Model) y la técnica LES (Large-Eddy Simulations). Siguiendo la idea de cascada de energía explicada en el epígrafe 2.2, para resolver la enorme cantidad de estructuras que se da en un régimen turbulento, DNS calcula todas las estructuras del flujo. Así pues, obtiene información instantánea del flujo a costa de ser muy exigente en requerimientos computacionales. En sentido opuesto, el método RANS no resuelve directamente las ecuaciones de Navier-Stokes, sino utiliza modelos para calcular el promedio en el tiempo de todo el flujo, con la ventaja de relajar considerablemente las necesidades de cálculo a cambio de perder la información instantánea, ya que sólo calcula el valor medio de las variables. Finalmente, el método LES resuelve los grandes vórtices y simula los pequeños mostrando la tendencia de las variables a lo largo del tiempo con unas necesidades computacionales intermedias entre los dos anteriores. La utilización de cada técnica depende de los requerimientos del caso a resolver. Por ejemplo, DNS debido a su coste computacional, queda fuera de la mayoría de aplicaciones de la ingeniería y se emplea para flujos con números de Reynolds bajos y confinados espacialmente, ya que los requerimientos computacionales aumentan exponencialmente con Re. En este sentido, no se ha encontrado ningún ejemplo en la literatura que estudie los flujos de ventilación a 2349 Dinámica computacional de fluidos (CFD) través de esta técnica. Para reducir los requerimientos computacionales, el método RANS promedia las ecuaciones de Navier-Stokes y utiliza un modelo de turbulencia para resolverlas. Asumiendo la descomposición de Reynolds, separando los términos turbulentos y no turbulentos, y sustituyendo en las ecuaciones de Navier-Stokes se obtienen las ecuaciones de forma cartesiana (3.1) y (3.2)0 que se emplean por RANS para resolver el flujo. ρ t + (ρu i) = 0 (3.1) x i (ρu i ) t + (ρu iu j ) x j = p + [μ ( u i + u j 2 x i x j x j x i 3 δ u l ij )] + (ρu ) i u j (3.2) x l x j Aunque tienen una forma similar a la de las ecuaciones de Navier- Stokes, las variables ahora representan la media espacial o temporal del conjunto. Es importante señalar que aparece un nuevo término, ρu, i u j llamado el tensor de los esfuerzos de Reynolds. El tensor de los esfuerzos de Reynolds introduce nuevas incógnitas, que deben ser resueltas mediante la introducción de nuevas ecuaciones, como detallaremos en el epígrafe 3.3. RANS es la técnica más utilizada en la actualidad para problemas de ingeniería en general, y en tecnología de la ventilación en particular, pues se consiguen resultados aceptables con un bajo consumo computacional. No obstante, se puede cuestionar si estas simplificaciones son adecuadas para el estudio de los flujos 2450 Dinámica computacional de fluidos (CFD) no permanentes que se producen en la capa límite atmosférica, pudiendo enmascarar la dinámica real de las estructuras turbulentas. Por ello, sólo es adecuada para estudios en los que los flujos ambientales sean estables en el tiempo. Por otra parte, aunque de forma minoritaria, existen trabajos que utilizan LES para el estudio ventilación cruzada en edificios, como recogen Rampony y Blocken (2012). No obstante, en invernaderos no se ha encontrado ningún estudio que utilice LES para resolver el flujo de ventilación Modelos de turbulencia que se emplean en la técnica RANS Como se comentó anteriormente, en el método RANS las componentes no permanentes de las variables de flujo se expresan por su valor promedio a través del tensor de esfuerzos de Reynolds, añadiendo seis incógnitas adicionales al sistema de ecuaciones. Para completar el número de ecuaciones y permitir la resolución del sistema generalmente se utilizan dos métodos de modelado, la aproximación de Boussinesq y los modelos de transporte de Reynolds que serán desarrollados a continuación Aproximación de Boussinesq (k-ε standard, RNG y SST k-ω) En la aproximación de Boussinesq se relacionan los esfuerzos turbulentos con el gradiente de velocidad del flujo medio mediante la ecuación (3.3). 2551 Dinámica computacional de fluidos (CFD) ρu i u j = μ t ( u i + u j ) 2 x j x i 3 (ρk + μ u k t ) δ x ij (3.3) k Donde μ t es la llamada viscosidad turbulenta y necesita de ecuaciones adicionales para ser resuelta. Una de las desventajas de la hipótesis de Boussinesq es que asume la viscosidad turbulenta como un escalar isotrópico, lo que no es estrictamente cierto. Para la resolución de la viscosidad turbulenta el modelo más utilizado en ingeniería es el llamado k-ε estándar (Launder & Spalding, 1972). Se trata de un modelo semi-empírico y se describe mediante dos ecuaciones de transporte, una para la energía turbulenta cinética (k), y otra para el ratio de disipación (ε). No obstante, no se trata del único modelo. El modelo k-ɛ sólo es válido para flujos totalmente turbulentos, ya que se desecha el efecto de la viscosidad molecular. La energía cinética turbulenta, k, y su rato de disipación, ɛ, se obtienen las siguientes dos ecuaciones de transporte: (3.4) y (3.5). (ρk) t + (ρku i) x i = x j [(μ + μ t σ k ) (k) x j ] + G k ρε + S k (3.4) (ρε) t + (ρεu i) x i = x j [(μ + μ t σ ε ) (ε) x j ] + C 1ε ε k (G k) C 2ε ρ ε2 k + S ε (3.5) En estas ecuaciones G k representa la generación de energía cinética 2652 Dinámica computacional de fluidos (CFD) turbulenta debido a los gradientes de la velocidad media; σ k y σ ε son los números turbulentos de Prandtl para k y ɛ respectivamente; S k y S ε son términos definidos por el usuario. Finalmente, el resto de parámetros son constantes determinadas experimentales con aire y agua para flujos turbulentos considerados fundamentales, como el flujo de Couette. Sin embargo, debido a su extendido uso se sabe que funcionan razonablemente bien para un gran rango de flujos sin cortadura y/o confinados por paredes. La viscosidad turbulenta, μ t, se computa combinando k y ɛ a través de la ecuación (3.6), donde C μ es una constate: μ t = ρc μ k 2 ε (3.6) El modelo k-ε presenta una serie de desventajas: En primer lugar, la energía cinética turbulenta se sobreestima en las regiones donde el flujo incide sobre obstáculos y se separa; así como en donde se vuelve a unir. Por esta razón estos modelos no son capaces de ajustarse a los flujos reales en zonas angulosas, como el techo de los edificios, o en el exterior de los invernaderos, en los que se produce la separación del flujo. Tampoco reproducen adecuadamente los flujos de recirculación que se dan a sotavento de los obstáculos, debido a que el campo de esfuerzos es muy complejo y está dominado por esfuerzos cortantes, por lo que la asunción de isotropía de la viscosidad turbulenta deja de ser válida 2753 Dinámica computacional de fluidos (CFD) Otras propuesta similar al modelo k-ɛ, que lo intenta mejorar, es el modelo RNG k-ɛ. El modelo RNG k-ɛ deriva de una técnica estadística llamada teoría del grupo de renormalización. Añade varios refinamientos al modelo k-ɛ estándar. Primero, el modelo RNG k-ɛ dispone de términos adicionales para la ecuación ɛ que mejoran significativamente la precisión para flujos rápidamente deformados. Segundo, el efecto de los flujos rotacionales en la turbulencia se incluye y aumenta la precisión de este tipo de flujos. El modelo RNG provee, además, de una formulación analítica para los números de Prandtl, mientras que el modelo k-ɛ estándar utiliza constantes especificadas por el usuario. Finalmente, mientras el modelo k-ɛ estándar es un modelo para números de Reynolds altos, el RNG proporciona una formula diferencial para la viscosidad efectiva, derivada analíticamente, que tiene en cuenta los efectos para flujos con números de Reynolds inferiores. Estas características hacen del modelo RNG k-ɛ un modelo apropiado para una mayor cantidad de flujos que el modelo k-ɛ estándar. Las ecuaciones de transporte son (3.7) y (3.8): (ρk) t + (ρku i) x i = x j (μ eff α ε ε x j ) + G k ρε + S k (3.7) (ρε) t + (ρεu i) x i = x j (μ eff α ε ε x j ) + C 1ε ε k (G k) C 2ε ρ ε2 k R ε + S ε (3.8) En estas ecuaciones, G k, representa la generación energía cinética 2854 Dinámica computacional de fluidos (CFD) turbulenta debido a los gradientes de velocidad media. Los valores α k y α ε son los números de Prandtl efectivos inversos para k y ε, respectivamente. Finalmente, S k y S ε son términos definidos por el usuario. El procedimiento de la eliminación de escalas en la teoría RNG resulta en una ecuación diferencial para la viscosidad turbulenta: Ecuación (3.9). d ( ρ2 k εμ ) = 1.72 v dv (3.9) v C v Donde v = 1,72 μ eff μ y c v 100. La ecuación (3.9) se integra para obtener una descripción precisa de como el transporte turbulento efectivo varía con el número de Reynolds, permitiendo al modelo manejar mejor los flujos con números de Reynolds bajos y los flujos cercanos a la pared. Para elevados números de Reynolds, la ecuación (3.9) se convierte en la ecuación (3.10). μ t = ρc μ k 2 ε (3.10) Con C μ = 0,0845, que se deriva de utilizar la teoría RNG. Es interesante hacer notar que este valor es muy similar al valor empíricamente derivado para el modelo k-ɛ estándar. Por defecto, el programa comercial ANSYS Fluent 2955 Dinámica computacional de fluidos (CFD) (Canonsburg, PA, USA) utiliza la forma para altos números de Reynolds. La principal diferencia de este modelo con el modelo k-ɛ estándar está un término adicional en la disipación de energía, ɛ, expresado por la ecuación (3.11). R ε = C μ ρƞ 3 (1 ƞ ƞ ) 0 ε βƞ 3 k (3.11) Donde ƞ Sk ε; ƞ 0 = 4,38 y β = 0,012. Reordenando la ecuación (3.8) con (3.11) obtenemos la ecuación (3.12) (ρε) t + (ρεu i) = ε ε (μ x i x eff α ε ) + C j x 1ε j k (G k) C 2ε ρ ε2 k (3.12) Donde C 2ε es:: C μ ƞ 3 (1 ƞ ƞ ) C 2ε C 0 (3.13) 2ε 1 + βƞ 3 Es interesante resaltar que en regiones donde η < η 0 el término R ε contribuye positivamente, y el término C 2ε es mayor que el término C 2ε. Como ejemplo, en la capa cercana a la pared, el valor de η 3.0, que es cercano a la magnitud de C 2ε del modelo k-ɛ estándar. Así pues, para flujos moderadamente 3056 Dinámica computacional de fluidos (CFD) deformados el modelo RNG da resultados similares al modelo k-ɛ estándar. En sentido contrario, para η > η 0, es decir regiones con un valor grande de deformación, el término R ε contribuye negativamente y el valor de C 2ε es más pequeño que el valor de C 2ε. Por tanto, comparándolo con el modelo k-ɛ estándar, el modelo RNG genera mayores valores de ɛ, reduce k y, finalmente, la viscosidad efectiva. Es decir, en flujos rápidamente deformados, el modelo RNG produce valores de la viscosidad turbulenta menores que en el modelo k-ɛ estándar. Por esta razón, el modelo RNG responde mejor al efecto de una rápida deformación y curvatura de las líneas de corriente que el modelo k-ɛ estándar, lo que explica que generalmente se ajuste mejor que el modelo k-ɛ estándar a algunas clases de flujos. Tabla 3 Constantes modelo standard k-ε, RNG y SST k-ω Constante Standard k-ε(-) RNG k-ε Constante SST k-ϖ C 1ε 1,44 1,42 σ k,1 1,176 C 2ε 1,92 1,68 σ w,1 2,0 C μ 0,09 0,0845 σ k,2 1,0 σ K 1,0 0,71 σ w,2 1,168 σ ε 1,3 0,71 α 1 0,31 ƞ 0-4,38 β i,1 0,075 β - 0,012 β i,2 0,0828 Otra variación a los modelos de Boussinesque es el modelo SST (shearstress transport) k-ω desarrollado por Menter (1994). Se diferencia del modelo k-ω estándar, desarrolado por Wilcox (1998) en cuatro puntos. Primero, permite un cambio gradual de la formulación que describe el flujo desde la parte interior de la capa límite, donde dominan los efectos viscosos y el número de Reynolds es bajo, donde utiliza la formulación del modelo k-ω estándar. En cambio, en la 3157 Dinámica computacional de fluidos (CFD) zona alejada de la capa límite, para elevados números de Reynolds, utiliza una formulación adaptada del modelo k-ɛ. En segundo lugar, permite modificar la formulación de la viscosidad turbulenta para tener en cuenta los efectos de transporte de los principales esfuerzos cortantes turbulentos. En tercer lugar, utiliza otras constantes (Tabla 3). Finalmente, se incluye la adición de un término de difusión cruzada en la ecuación ω y en la función de mezcla para asegurar que las ecuaciones del modelo se comportan apropiadamente tanto en las zonas cercanas a la capa límite como en las zonas alejadas. Estas dos características hacen del modelo SST k-ω más preciso y fiable para una mayor tipología de flujos, entre ellos los que se dan a través de gradientes adversos de presión. Así pues, las ecuaciones de transporte para k y ω son (3.14) y (3.15): (ρk) t + (ρku i) x i = x j (Γ k k x j ) + G k Y k + S k (3.14) (ρω) t + (ρωu i) x i = x j (Γ ω ω x j ) + G ω Y ω + D ω + S ω (3.15) En estas dos ecuaciones G k representa la generación de energía turbulenta cinética debido a los gradientes de velocidad. G ω representa la generación de ω; Γ k y Γ ω representan la difusividad efectiva de k y ω, respectivamente; Y k y Y ω representan la disipación de k y ω debido a la turbulencia. D ω representa el término de difusión cruzada. Finalmente, S k y S ω 3258 Dinámica computacional de fluidos (CFD) son términos definidos por el usuario Modelos de transporte de Reynolds (RSM) El Modelo de Esfuerzos de Reynolds (RSM, Reynolds Stress Model), o Cierre de segundo orden (SMC, Second Closure Model) fue desarrollado por Launder et al. (1975). RSM utiliza una ecuación para cada uno de los esfuerzos de Reynolds (6 ecuaciones en tres dimensiones), y otra ecuación adicional para la disipación de la energía turbulenta cinética (lo que aumenta la necesidad de recursos computacionales). Por ello, generalmente mejora la descripción de los flujos anisotrópicos. No obstante, la fidelidad de la precisión de los modelos RSM está todavía limitada por las asunciones empleadas para modelar algunos términos en las ecuaciones de transporte de los esfuerzos de Reynolds. En efecto, modelar los campos de presión y los términos del ratio de disipación son particularmente difíciles, y muchas veces se les considera responsables de comprometer la precisión de los modelos RSM. En este sentido, RSM no siempre produce mejores resultados que los modelos de dos ecuaciones, sin embargo, es un modelo que se tiene que tener en cuenta cuando las características estudiadas del flujo de interés son el resultado de fenómenos anisotrópicos, como por ejemplo los flujos muy rotacionales que se dan en las zonas de recirculación de los invernaderos. Finalmente, hay que resaltar que este método no se utiliza de muy a menudo porque que es más complejo y costoso que los modelos basados en la aproximación de Boussinesq. La forma exacta de las ecuaciones de transporte para los esfuerzos de Reynolds pueden escribirse de la siguiente manera: Ecuación (3.16) 3359 Dinámica computacional de fluidos (CFD) (ρu ) i u j t Derivada temporal local + (ρu ku ) i u j x k C ij Convección = [ρu x i u j u k + p(δkj ] u i + δ ik u j ) k D T,ij Difusión turbulenta + [μ (u ) i u j ] x k x k D L,ij Difusión molecular p (u i u u j k + u x j u u i k ) k x k P ij Producción de esfuerzo (3.16) + p ( u i + u j ) x j x i Φ ij Deformación de presión 2pΩ k (u ε j u m ikm + u i u ε m jkm ) 2μ u i u j x k x k F ij Producción del sistema rotacional + S usuario Término definido por el usuario ε ij Disipación De todos los términos definidos en esta ecuación, C ij, D L,ij, P ij y F ij no requieren ningún tipo de modelo y se resuelven directamente. En cambio, D T,ij, Φ ij y ε ij sí necesitan ser modeladas para cerrar las ecuaciones. Finalmente, la ecuación adicional para la disipación de la energía turbulenta cinética es muy similar a la utilizada para el modelo k-ɛ estándar Resolución del flujo cercano a la pared (Modelo de pared) Numerosos experimentos demuestran que la región cercana a una pared puede ser dividida en tres capas: -la capa viscosa, en la que el flujo es laminar y en la que la viscosidad 3460 Dinámica computacional de fluidos (CFD) juega un papel dominante en la transferencia de momento o de masa turbulencia -la capa exterior, o capa turbulenta, dominada por los fenómenos de -capa intermedia, en la que los dos fenómenos son igualmente importantes (Figura 3). Para reducir el número de celdas en la dirección normal de las paredes, donde se dan todos estos fenómenos de transición, se utilizan las funciones de pared como una alternativa a calcular los tensores de esfuerzos cortantes en las celdas, pues en estas zonas se requeriría una malla muy densa que conllevaría grandes necesidades de cálculo. Las funciones de pared (wall functions) son ecuaciones semi-empíricas que describen aproximadamente el comportamiento del flujo dominado por los efectos de la viscosidad y su transición al régimen turbulento dominante. Esta aproximación suele utilizarse porque es económica, robusta y proporciona resultados con una aceptable precisión. La función de pared más utilizada en la CFD es el modelo estándar de Launder y Spalding (1974). Este modelo utiliza tres variables adimensionales: una para describir la velocidad, u + = U/u ; otra para la distancia normal a la pared, y + = u y/θ; y otra que describe la rugosidad en función de los experimentos de Nikuradse (1933), k + s = u k s /θ. Donde U es la velocidad media tangente a la pared, u es una velocidad de fricción de la función de pared, k s es la altura de rugosidad arenosa equivalente, valor obtenido experimentalmente por Nikuradse (1933) 3561 Dinámica computacional de fluidos (CFD) para una gran variedad de superficies y ν es la viscosidad cinemática del fluido. Si representamos la velocidad, u +, obtenida experimentalmente, frente a la distancia, y +, en una escala logarítmica para diferentes rugosidades obtenemos la Figura 3. Figura 3 Ley de pared para superficies suaves (línea de puntos) y arenosas (línea continua) para diferentes ks+ (Blocken, et al., 2007) Observamos que la velocidad media, en la subcapa lineal, aumenta con la distancia a la pared manteniendo la proporción u + = y +, para y + < 5 A medida que nos alejamos todavía más de la pared, en la capa turbulenta, para y + > 30, la velocidad sigue una función logarítmica de pendiente (1/κ) pero con diferentes puntos de corte para diferentes k + s. En esta zona se distinguen tres regímenes en función de la k + s : aerodinámicamente 3662 Dinámica computacional de fluidos (CFD) suave ( k s + < 2,25), transitorio y totalmente rugoso ( k s + > 90 ). En la capa superficial de la atmósfera generalmente se considera un régimen totalmente rugoso ya que los obstáculos son grandes y la capa laminar es ínfima y por tanto el flujo se considera independiente de la viscosidad molecular En la mayoría de códigos comerciales, para resolver el flujo cercano a la pared, la rugosidad de la pared se implementa dando un valor de k s. En ventilación, este parámetro se obtiene a partir de otro parámetro experimental, la longitud rugosa ( z 0 ). La longitud de rugosidad (Roughness length) es un parámetro experimental equivalente a la altura a la que la velocidad del viento se considera nula, que depende de la rugosidad del terreno y se obtiene mediante la clasificación de Davenport (Tabla 4). Generalmente se relaciona con la k s mediante la ecuación (3.17). k s = z 0 exp(κb) (3.17) No obstante, es posible obtener otro tipo de relaciones entre estas dos constantes experimentales como veremos en el epígrafe63 Dinámica computacional de fluidos (CFD) Tabla 4 Clasificación revisada de Davenport para obtener la longitud rugosa (z o ) en función del tipo de suelo (Wieringa, 1992) Clase Nombre Longitud rozamiento 1 Mar 0, Suave 0,005 3 Abierto 0,03 4 Abierto rugosamente 0,1 5 Rugoso 0,25 6 Muy rugoso 0,5 7 Cerrado 1,0 8 Caótico 2 Descripción del paisaje Agua abierta, planicie oceánica, libre de obstáculos >3 km Sin ninguna característica especial apenas cubierto o hielo Terreno plano con césped o vegetación muy baja, con obstáculos ampliamente separados; pista aeropuerto Áreas cultivadas, cultivos bajos, obstáculos ocasionales separados por más de 20 alturas Paisaje abierto, cultivos de alturas variables, setos dispersos etcétera con una separación de 15 alturas Paisajes fuertemente utilizados con espacios abiertos de aproximadamente 10 alturas; setos, frutales bajos, bosques jóvenes densos Cubierta totalmente obstaculizado con espacios abiertos de aproximadamente 1 altura; busques maduros, edificaciones de alturas bajas. Distribución irregular de elementos muy grandes: centro de ciudades, grandes bosques con grandes claros 3864 Dinámica computacional de fluidos (CFD) 3.5. Esquemas de cálculo y residuos Cuando se calculan términos convectivos, los valores de las propiedades en las caras se utilizan en el cálculo junto con los valores centrales de cada celda. Estos valores son interpolados desde los valores centrales ofreciendo un valor por cada una de las caras de la celda (Figura 4). Existen varios esquemas de interpolación. El esquema por defecto utilizado en la mayoría de códigos es el llamado esquema de primer orden (first-order upwind scheme), donde se asume que para variable, Φ, el valor en cada cara es el mismo que en el centro de la celda aguas arriba. Las dos principales ventajas de este método son su facilidad para ser implementado en cualquier código y que evita inestabilidades en los cálculos. No obstante, puede introducir errores en los flujos convectivos debido a un posible incremento en los términos difusivos ya que arrastra los errores a medida que resuelve el flujo aguas abajo. En cambio, en el esquema de segundo orden (second-order upwind scheme), se determina el valor de Φ desde el valor de dos celdas aguas arriba de la cara considerada. Este método suele ser más preciso que el de primer orden, pero en regiones con grandes gradientes puede resultar en valores de cara fuera de los rangos de las celdas, por lo que generalmente se aplican límites a los valores predichos. 3965 Dinámica computacional de fluidos (CFD) Figura 4 Representación esquemática del cálculo de (izquierda) primer orden y de (derecha) segundo orden. En general, la convergencia de las simulaciones de CFD se monitorea mediante los residuos. Para el algoritmo de cálculo basado en la presión que se utiliza en esta tesis, el residuo, R Φ, se define con el desequilibrio en la ecuación (3.18) sumada sobre todas las celdas computacionales del domino (P). Es decir: R Φ = a nb Φ nb + b a p Φ p nb (3.18) Estos residuos de la variable Φ se denominan como no escalados. Los residuos normalmente se escalan por un factor representativo del flujo de la variable Φ a través del dominio. Generalmente se utiliza el valor local de la variable Φ: Ecuación (3.19). R Φ s = p nb a nbφ nb + b a p Φ p a p Φ p (3.19) 4066 Dinámica computacional de fluidos (CFD) Con infinita precisión en el cálculo se alcanzaría una solución que convergería a cero. No obstante, en las computadoras el mínimo residuo que se puede alcanzar es el que ofrece la precisión de cálculo. Así se trabaja generalmente con dos precisiones de cálculo: precisión sencilla que permite caer a los residuos seis órdenes de magnitud, y doble precisión permite reducir hasta doce órdenes de magnitud. En el Anejo A se puede encontrar más información sobre la discretización y linealización de las ecuaciones de transporte. 4167 Dinámica computacional de fluidos (CFD) 3.6. Análisis crítico de guías de buenas prácticas en CFD. Dado que en la revisión bibliográfica no se han encontrado guías de buenas prácticas para simulaciones CFD en invernaderos se ha realizado una búsqueda y análisis crítico de guías de buenas prácticas en otros campos de la tecnología con la intención de adaptarlas para el caso concreto de tecnología de invernaderos. Principalmente, nos hemos centrado en buscar dichas guía en ingeniería de la edificación. Efectivamente, en los últimos años ha habido un interés creciente en utilizar CFD para el estudio del ambiente urbano (ambiente exterior) y los fenómenos de ventilación en los edificios (ambiente interior). El interés en realizar modelos del ambiente exterior en los entornos urbanos está relacionado con la salud de sus habitantes y por ello se estudia la emisión, dispersión y concentración de contaminantes. Los modelos del ambiente interior están relacionados con el confort en el interior de las diferentes dependencias de los edificios. Asimismo, se han añadido guías de ingeniería naval o aeroespacial (Tabla 5). De la misma forma, se incluyen otros autores que no pretenden realizar guías de buenas prácticas pero que sus consideraciones resultan útiles a tal fin. El esquema elegido para sistematizar una simulación CFD es propuesto por la guía COST 732 (Schatzmann, et al., 2010). 4268 Dinámica computacional de fluidos (CFD) Tabla 5 Revisión de las principales guías de buenas prácticas Nombre (Autor) Campo de actuación Característica Guide for the Verification and Validation of Computational Fluid Dynamics Simulation (AIAA, 1998) ERCOFTAC (Casey & Wintergerste, 2000) ECORA (Menter, et al., 2002) QNET-CFD (Bartzis, et al., 2004) VDI guideline: (Technical Division Environmental Meteorology, 2005) AIJ Guidelines for practical applications of CFD (Tominaga, et al., 2008) COST732 (Schatzmann, et al., 2010) Simulaciones CFD aplicada a la industria aeroespacial Simulaciones CFD para aplicaciones generales en la industria Simulaciones CFD para aplicaciones relacionadas con la seguridad de reactores nucleares Simulaciones CFD para la Ingeniería civil (HVAC) Simulaciones para flujo alrededor de edificios Simulaciones CFD para ambientes urbanos a escala de los peatones Simulaciones CFD para ambientes urbanos Se centra en la verificación de códigos CFD y en el proceso de validación Documento más completo hasta la fecha. No se ocupa intencionadamente de problemas relacionados con la capa límite atmosférica Modificación de las guías ERCOFTAC extendidas. Las mejores guías de buenas prácticas en ingeniería eólica hasta la fecha Documento que evalúa los modelos micro-escala que explícitamente resuelven obstáculos en la ABL teniendo en cuenta la orografía Comparación cruzada entre predicciones CFD, túneles de viento y medidas a escala real en 7 casos para diferentes condiciones de contorno A través de la bibliografía genera una guía dirigidas específicamente a validaciones CFD en medios urbanos. 4369 Dinámica computacional de fluidos (CFD) Elección de las variables objetivo Schlünzen (1997) y Menter, et al, (2002) recomiendan que en el primer paso de toda simulaciones CFD se debe definir las variables objetivo. Esto es, se deben incluir aquellas variables que son representativas de los objetivos de la simulación y aquéllas que puedan ser comparadas con sus correspondientes experimentales para la validación del modelo. Las variables más utilizadas, tanto en tecnología de la edificación como en invernaderos, que se miden directamente son la temperatura, la humedad, la velocidad y dirección del viento y la concentración de especies químicas. La tasa de renovación es la variable medida indirectamente más importante en estudios de ventilación de invernaderos Elección de las ecuaciones de aproximación para describir la física del flujo La elección de las ecuaciones básicas influye en enormemente en la precisión de los modelos. Generalmente en los modelos ambientales se utilizan las simplificaciones RANS con el método standard k-ε. Desafortunadamente, no es posible realizar recomendaciones sobre el modelo de turbulencia más apropiado para cada caso, puesto que el rendimiento de la simulación es muy dependiente del modelo y de la malla utilizada. (Franke, et al., 2007). No obstante, repasando lo explicado en el epígrafe 3.3 sí es posible sacar algunas recomendaciones generales para utilizar un modelo de turbulencia. El modelo k-ε estándar está indicado para flujos con altos números de 4470 Dinámica computacional de fluidos (CFD) Reynolds, sin embargo, no es aconsejable en flujos donde el flujo incide sobre obstáculos y se separa o se une. Por esta razón estos modelos deben de evitarse para estudiar zonas angulosas. El modelo k-ε estándar tampoco es aconsejable en los flujos de recirculación que se dan a sotavento de los obstáculos. El modelo RNG k-ɛ está indicado para flujos rápidamente deformados y para los flujos de recirculación. Además, funciona mejor que el modelo k-ε estándar para menores números de Reynolds. Estas características hacen del modelo RNG k-ɛ un modelo apropiado para una mayor cantidad de flujos que el modelo k-ɛ estándar, sobre todo en aquellos en los que existen rápidas deformaciones del flujo o curvaturas de las líneas de corriente. El modelo SST k-ω, también está indicado para flujos rápidamente deformados y para flujos de recirculación. Además, produce buenos resultados para flujos donde se dan gradientes adversos de presión. Finalmente el modelo RSM está indicado para todo tipo de flujos anisotrópicos a costa de ser mucho más costoso computacionalmente. No obstante, la fidelidad de la precisión de los modelos RSM está todavía limitada por las asunciones empleadas para modelar algunos términos en las ecuaciones de transporte de los esfuerzos de Reynolds Elección de la representación geométrica de los obstáculos Las características del área que rodea a un obstáculo (por ejemplo, un invernadero) son determinantes en el comportamiento del flujo. Esto incluye la presencia de otros edificios u obstáculos, así como, la vegetación predominante 4571 Dinámica computacional de fluidos (CFD) en (árboles, césped, etc.) o el tipo de suelo (arenoso, pedregoso, etc.) en dicha área, por lo que la malla en estas zonas debe tener suficiente resolución. COST732 recomienda valorar la influencia de la omisión de detalles en la simulación con la comparación de datos experimentales o la realización de test de sensibilidad si estas características no son conocidas. Tominaga et al. (2008) recomiendan que la región que rodea un edificio, en un radio de 1-2 alturas, esté claramente representada. Además, recomiendan utilizar simplificaciones geométricas de los obstáculos de los alrededores o especificar parámetros apropiados de rozamiento para el suelo. Finalmente, Yoshie et al. (2006) señalan que existe una posibilidad de obtener resultados poco realistas si la extensión computacional se expande sin representar adecuadamente las principales características físicas de los alrededores del edificio estudiado. Es importante indicar que COST732 señala que, en experiencias realizadas en túneles de viento, los edificios con una altura H apenas influyen sobre el viento a distancias superiores a 6 H Elección del dominio computacional COST732 recoge las recomendaciones que, para un único edificio, la extensión vertical del dominio computacional debe de ser al menos 5 veces su altura. Por otro lado, otras fuentes, como VDI (Technical Division Environmental Meteorology, 2005), utilizan el denominado ratio de obstrucción, que se define como la proporción entre el área proyectada del edificio en la dirección del flujo respecto al área de la sección cruzada del domino computacional. VDI sugiere 4672 Dinámica computacional de fluidos (CFD) que la extensión vertical del domino debe de ser entre 4 (para ratios bajos) y 10 (para los mayores ratios) veces la altura del edificio (H), siendo 10 % el ratio máximo aconsejado para una simulación. No obstante, los trabajos realizados por Batke et al. (1990) en túneles de viento sugieren que dicho ratio no exceda de 3%. Las recomendaciones de Tominaga et al. (2008) van en este mismo sentido. Una vez elegida la extensión vertical, y manteniendo el ratio de obstrucción inferior al límite recomendado, se debe calcular la extensión de los laterales. Por ejemplo, para un edificio de igual anchura que altura, una extensión vertical del domino de 6 H, y un ratio de obstrucción 3%, se necesitaría una extensión lateral del dominio de 2,3 H. Sin embargo, Tominaga et al. (2008) consideran las recomendaciones de COST732 insuficientes y recomiendan una extensión lateral mínima de 5 H. En los edificios más extensos que altos, entre los que podrían estar los invernaderos, se recomienda que el ratio de la extensión lateral sobre la extensión vertical del dominio sea similar al de la altura respecto a la longitud del edificio. Si el perfil de viento que se aproxima al edificio es conocido, COST732 recomienda que exista al menos una distancia de 5 H entre el lugar donde se aplican las condiciones de contorno de entrada y el edificio. Bartzis, et al. (2004) recomiendan utilizar una distancia mayor si el perfil no es conocido. VDI recomienda desde una distancia de 2 H para ratios de obstrucción del 3% hasta el 8 H para ratios del 10%. En la zona a sotavento, COST732 recomienda una extensión de 15 H 4773 Dinámica computacional de fluidos (CFD) para permitir la convergencia del flujo Elección de las condiciones de contorno Simulación del perfil de viento A menudo no se dispone de información experimental sobre el perfil de viento que incide sobre los edificios estudiados. En estos casos, Richards y Hoxay (1993) proponen que se un perfil logarítmico de la velocidad (u (z) ; v (z) ; w (z) ); un perfil hiperbólico de la energía de disipación turbulenta (ɛ (z) ) y un perfil constante para la energía cinética turbulenta (k), como se indican en las ecuaciones (3.20), (3.20) y (3.22), siendo (z 0 ) la altura: u(z) = u ABL κ ln (z + z o) z 0 v(z) = 0 w(z)=0 (3.20) k = u 2 ABL (3.21) C μ ε(z) = u 3 ABL κ(z + z 0 ) (3.22) Donde κ es la constante de Karman, y u ABL es la velocidad de fricción, que se define: Ecuación (3.23). κu h u ABL = ln ( (h + z 0) z ) (3.23) 4874 Dinámica computacional de fluidos (CFD) Siendo u h la velocidad a una altura específica h, generalmente 2 metros, obtenida generalmente de estaciones meteorológicas cercanas. En cambio, Tominaga et al. (2008) utilizan una función exponencial. Su propuesta es interesante, ya que permite utilizar información de la distribución vertical de k a través de su relación con la intensidad turbulenta, lo que no ocurre en las guías COST732. Existen otras funciones para la distribución vertical de k propuestas en los trabajos de Yang y Zhang (2009) y Gorlé et al. (2009). No obstante, no se han encontrado trabajos en la literatura donde se apliquen. Simulación del efecto del suelo COST732 indica que para tratar el flujo cercano a las paredes se pueden utilizar los modelos semi-empíricos conocidos como funciones de pared (epígrafe 3.4). No obstante, advierte que estas funciones de pared no son adecuadas para zonas donde se producen fenómenos como la separación o convergencia del flujo. Sin embargo, Castro (2003) muestra que los errores debidos a la mala predicción en estas zonas pueden tener poca influencia en el resultado global de la simulación. Hargreaves y Wright (2007) llaman la atención sobre los errores que se comenten al considerar que el perfil de viento está en equilibrio con el resto de condiciones de contorno y que no sufre variaciones en la zona de aproximación, puesto que el perfil de viento y los parámetros asociados a la turbulencia decaen en la zona de aproximación en la mayoría de las simulaciones. El problema se debe a que la condición de entrada de aire no está balanceada con los modelos 4975 Dinámica computacional de fluidos (CFD) de turbulencia, la función de pared y las condiciones para conducir el flujo. Esto provoca que el perfil de la capa límite atmosférica cambie a medida que se avanza por el dominio y se obtengan perfiles de viento inexactos a barlovento, cerca del edificio. Las diferencias entre los perfiles reales y simulados pueden disminuir el éxito de las simulaciones significativamente (Blocken, et al., 2007). Blocken et al. (2007) establecen cuatro requerimientos para simular correctamente la rugosidad (Figura 5): Primero, se debe de tener una resolución suficientemente fina de malla en la dirección vertical de la pared. Segundo, se debe de mantener un perfil de viento homogéneo en la dirección del flujo en las regiones aguas arriba y aguas abajo del dominio. Tercero, se tiene que mantener una distancia y p desde el punto P, punto central de la celda adyacente a la pared, que sea superior a la altura de rugosidad del terreno k s, esto es y p > k s. Finalmente, se debe de satisfacer la relación entre z 0 y k s explicadas en el epígrafe 3.4. Figura 5 Diagrama de la primera celda computacional. Perfiles de velocidad del modelo de pared y del perfil logarítmico. (Blocken, et al., 2007) No obstante, es prácticamente imposible satisfacer los cuatro 5076 Dinámica computacional de fluidos (CFD) requerimientos. El problema más usual es que mantener la condición y p > k s que, generalmente lleva a grandes celdas en contacto con la pared. Usar valores más pequeños de k s permite una mejor resolución en la pared, pero produce una modificación indeseada del perfil de viento a medida que éste avanza en el dominio. En este sentido, Ramponi y Blocken (2012) realizaron un estudio de sensibilidad de varios parámetros para las simulaciones CFD, y para ventilación natural señalan que la diferencia entre perfiles de entrada en el dominio y el perfil que incide sobre el edificio es uno de los parámetros que más afecta al resultado de las simulaciones, y en particular sobre el parámetro k. En el caso concreto de ANSYS Fluent la condición de pared se define mediante la ecuación (3.24). U P u u τ 2 = 1 κ ln ( Eu y p ν(1 + C s k s + ) ) (3.24) Donde E es una constante empírica con un valor aproximado de 9,793; y u y u τ son dos velocidades de fricción; k p es la energía turbulenta en el punto central P de la celda adyacente a la pared. El término (1 + C s k + s ) permite al usuario introducir el efecto para diferentes rugosidades a través de los parámetro k s y C s, siendo C s una constante de rugosidad de la que no existen guías específicas al respecto y su valor por defecto es 0,5 (-) (ANSYS Fluent, 2010). En comparación con los ambientes urbanos, donde generalmente la k s es muy elevada debido a la presencia de edificios alrededor del edificio de 5177 Dinámica computacional de fluidos (CFD) estudio, en los invernaderos los valores son inferiores. Así, por ejemplo, para una superficie con vegetación baja, muy común alrededor de los invernaderos, y con una separación amplia entre obstáculos, según la clasificación de Davenport, el valor de z o es de 0,015 m, y para mantener el requerimiento mínimo de celda y p > k s, para un valor por defecto de C s = 0,5, obtenemos un valor aproximados k s,abl de 30 cm, es decir que como mínimo, la celda adyacente a la pared tendrá un tamaño mínimo de lado de 60 cm. Que daría como resultado una malla de poca resolución en la zona próxima a la pared. Si aumentamos el parámetro C s a su máximo valor, 1 (-), daría como resultado una celda de 30 cm, que permite una resolución más acorde con las recomendaciones del epígrafe Blocken et al. (2007) indican otras maneras de relajar el requerimiento y p > k s usando otro tipo de funciones de pared, pero esto generalmente no es posible cuando se utilizan programas comerciales. En estos casos sugiere especificar explícitamente la condición τ w = ρ(u ABL ) lo cual fuerza a igualar la fricción de la ABL con la de la pared, lo que reduce considerablemente los errores. No obstante, COST732 recomienda que, debido a la dificultad de satisfacer todas las condiciones en la condición de contorno del suelo, se realice un análisis de sensibilidad en un dominio sin obstáculos, con el fin de estudiar cómo se compatibilizan todas las condiciones de contorno. Por este motivo, no sólo hay que dar información del perfil de entrada al dominio computacional sino también del perfil incidente sobre el edificio. 5278 Dinámica computacional de fluidos (CFD) Finalmente, las diferentes guías advierten de que estos requerimientos se han elaborado para el modelo k-ε estándar, pero que también deben ser tenidos en cuenta en el resto de modelos. Condiciones de contorno en la parte superior del dominio COST732 indica que la elección de esta condición de contorno es importante para mantener el perfil de incidencia de viento igual al de entrada. En este sentido, aconseja utilizar un esfuerzo cortante definido por el usuario en la parte superior del dominio que mantenga iguales el perfil de incidencia y el perfil de entrada. No obstante, COST732 señala que es muy usual utilizar con buenos resultados la condición de deslizamiento libre (sin esfuerzo cortante) en esta zona del dominio. En cambio, indica que no es aconsejable usar la condición de simetría, consistente en generar un flujo paralelo al contorno del dominio, lo que requiere que la componente normal de la velocidad desaparezca, ni la condición de outflow, condición de contorno que permite una componente normal de la velocidad, ya que pueden modificar el perfil de viento atmosférico. No obstante, Tominaga y otros (2008) indican que si el dominio computacional es lo suficientemente alto, es decir, se atiende a las recomendaciones del epígrafe 3.6.4, esta condición de contorno no tiene influencia alrededor del edificio objetivo de estudio. Condición de contorno lateral y salida del dominio En los laterales es común suponer la condición de simetría cuando la dirección del flujo es paralela a los mismos. Se utilizan generalmente las condiciones de contorno llamadas abiertas, que son dos, outflow y pressure 5379 Dinámica computacional de fluidos (CFD) outlet con las que todas las derivadas de las variables del flujo son forzadas a cero. Condiciones de contorno de la superficie de los edificios En edificación generalmente se utiliza la función de pared. Ramponi y Blocken (2012) utilizan k s = 0, considerando que la rugosidad es despreciable Elección de las condiciones iniciales COST732 señala que las condiciones iniciales son más importantes cuanto más largo es el dominio o más pequeña es la velocidad del viento. No obstante, cuando se emplea el método RANS, la influencia de las condiciones iniciales es limitada y el resultado final depende casi únicamente de las condiciones de contorno. En cualquier caso, elegir una condición de entrada lo más parecida posible a la solución final reducirá significativamente los esfuerzos computacionales. En este sentido, no es infrecuente utilizar como condiciones iniciales resultados de simulaciones más sencillas, como simulaciones del flujo laminar o esquemas de cálculo de primer orden Elección de la malla COST732 recomienda realizar un estudio de sensibilidad de la malla a partir de un tamaño mínimo de celda de un décimo de la altura del edificio principal. Sugiere que en el área de interés, la resolución mínima debe de ser al menos de 10 celdas contenidas en 1 m 3. También, puede ser interesante utilizar geometrías simplificadas para valorar las resoluciones iniciales (Blocken, et al., 2004). 5480 Dinámica computacional de fluidos (CFD) Para simulaciones cuyo objetivo sea la validación se utiliza la llamada extrapolación de Richardson, en la que al menos se deben realizar simulaciones con tres mallas con resoluciones que aumentan la resolución 1,5 veces en las tres dimensiones. Tominaga et al. (2008) indican que para predecir el flujo alrededor del edificio con una precisión aceptable, es muy importante reproducir adecuadamente las características de separación del flujo cerca del techo y de las fachadas, pues es donde se dan los fenómenos más complejos, por lo que estas zonas es recomendable que la malla sea más densa. COST732 recomienda que la malla debe de tener un índice de forma adecuado, pero ninguna guía da ninguna recomendación sobre su valor. El ratio de crecimiento entre celdas consecutivas (stretching ratio) debe de ser pequeño en las regiones donde existan mayores gradientes en las propiedades del flujo para mantener los errores de truncamiento pequeños. COST732 en este sentido recomienda que dicho ratio sea menor de 1,3 (-). Otros autores, sin embargo, recomiendan ratios inferiores a 1,2 (-) (Grigoriadis, et al., 2004). Por su parte Tominaga y otros (2008) recomiendan realizar pruebas disminuyendo el stretching ratio, ya que este puede depender de la forma del edificio estudiado Elección de las aproximaciones numéricas y de los criterios de convergencia Tominaga el al. (2008) señalan que los esquemas de solución de primer orden no son apropiados para todas las variables de transporte. Así pues, indican que no se deben de usar. En el mismo sentido se pronuncia COST732, y 5581 Dinámica computacional de fluidos (CFD) recomienda que, si se usa un esquema de simple precisión se debe de demostrar, mediante un análisis de sensibilidad, que los resultados de la simulación no se ven afectados. COST732 indica que las iteraciones deben realizarse al menos hasta que la reducción de residuos sea de cuatro órdenes de magnitud. Además, también recomienda hacer un seguimiento de las variables que forman parte de las ecuaciones de Navier-Stokes utilizadas en la simulación. 5682 Dinámica computacional de fluidos (CFD) 3.7. Análisis crítico de los estudios de la ventilación natural en invernaderos mediante CFD El primer trabajo de CFD aplicado a invernaderos que se tiene constancia fue publicado por Okushima et al. (1988), que realizaron una aproximación de la ventilación natural de un invernadero pequeño de cristal. Desde entonces ha habido un incremento notable en la complejidad así como en el número de los estudios llevados a cabo en invernaderos. La CFD permite de una forma sencilla cuantificar espacialmente la ventilación y por tanto contribuir a la mejora del diseño de las instalaciones y optimizar estrategias encaminadas a incrementar la producción. Además, la CFD soslaya los problemas de semejanza, relacionados con los números adimensionales, y proporciona un control total sobre las condiciones de contorno de una simulación. Esto supone, por ejemplo, realizar estudios paramétricos o de sensibilidad de forma mucho más sencilla que en un túnel de viento o agua. Además, permite combinar los parámetros espacio-temporales y, por tanto, es posible calcular, por ejemplo, los efectos de la variación de la dirección de viento en la ventilación. La ventilación natural es un fenómeno complejo y supone un reto porque depende de las condiciones ambientales exteriores, por lo que es necesario tener un gran conocimiento de física e ingeniería del viento. A pesar su idoneidad, la CFD solo se ha empleado para estudiar de una forma extensiva el clima dentro de los invernaderos con una configuración determinada de ventanas (Bournet & Boulard, 2010). Pese a la importancia de la validación de los modelos, sólo algunos 5783 Dinámica computacional de fluidos (CFD) autores se han propuesto validar el flujo de la ventilación natural, entre los que destacamos los siguientes: Bartzanas et al. (2007) pretenden analizar la sensibilidad del modelo al cambiar el modelo turbulencia; Teitel et al. (2008) estudian el efecto de la dirección del viento en la ventilación; y, finalmente, Majdoubi et al. (2009) estudian los patrones de los campos de velocidad y temperatura (Tabla 6). La validación de los dos primeros estudios se realizó mediante medidas de la dirección y velocidad del viento con anemómetros sónicos. En el trabajo de Bartzanas et al. (2007) se midió la velocidad y dirección del viento mediante un único anemómetro sónico 3D, que se desplazó para obtener información en 15 posiciones; mientras que la información del exterior del invernadero se obtuvo de una estación meteorológica cercana. Por su parte Teitel et al. (2008) utilizaron 3 anemómetros sónicos para medir en el interior del invernadero en 9 posiciones diferentes; el flujo exterior lo midió con 1 sensor cazoleta-veleta cercano. Por otra parte, Majdoubi et al. (2009) utilizaron 20 termo-higrómetros situados simultáneamente en distintas posiciones en el interior del invernadero. La velocidad del aire exterior se midió con un anemómetro de cazoleta emplazado en el mismo invernadero a una altura de 3 metros. Los dos primeros estudios utilizaron el modelo k-ɛ estándar. Ninguno indica el esquema de cálculo utilizado, lo que indica que probablemente se utilizaron los esquemas de cálculo de primer orden. Los tres estudios utilizaron el programa comercial ANSYS Fluent para calcular las aproximaciones numéricas. En relación a la representación geométrica de los obstáculos, ninguno 5884 Dinámica computacional de fluidos (CFD) consideró obstáculos cercanos ni ofrecieron información de por qué no se consideraron. Respecto al dominio computacional, siguiendo sólo las indicaciones COST732, es decir, sin atender las recomendaciones de Tominaga, et al. (2008) y obviando los ratios de obstrucción de VDI, se observa que los dos primeros estudios cumplen con las recomendaciones mínimas de extensión vertical. Respecto a la extensión de entrada, los tres estudios cumplen con las recomendaciones COST732. En cambio, el domino de salida queda insuficientemente definido en el estudio de Bartzanas et al. (2007) y Majdoubi et al. (2009), siendo en éste último caso 3 veces inferior al recomendado. Por otro lado, el trabajo de Teitel et al. (2008) define un dominio de unas 10 veces superior al recomendado. Finalmente, en cuanto al dominio lateral, los tres trabajos cumplen las recomendaciones COST732 y una vez más el trabajo de Teitel et al. (2008) definen un dominio muy superior (18 veces) al recomendado. Respecto a las condiciones de contorno, el trabajo de Bartzanas et al. (2007) no ofrece información sobre ningún parámetro. En los otros dos estudios, se mantienen las recomendaciones COST732 exceptuando el trabajo de Teitel et al. (2008) que definió una condición de simetría en la parte superior del dominio. En cuanto a la condición de entrada del dominio, Bartzanas et al. (2007) definieron un perfil constante de viento, mientras que los otros dos estudios definieron un perfil logarítmico, más acorde con las recomendaciones COST732, que además, utilizaron como condiciones de contorno. Finalmente, ningún estudio menciona el modelo de pared utilizado o sus parámetros; ni el tamaño mínimo de celda o las rugosidades definidas en el inferior del dominio o superficie del invernadero, exceptuando el estudio de Teitel et al. (2008), en el que se consideraron todas 5985 Dinámica computacional de fluidos (CFD) las superficies con la condición de no deslizamiento. Respecto a la definición de la malla, sólo Majdoubi et al. (2009) generan una malla estructurada, pero no ofrecen información sobre el tamaño mínimo de celda, el crecimiento entre dos celdas consecutivas o la calidad de la malla. Respecto a la selección de la malla, ningún estudio sigue las recomendaciones COST732 para la realización de las mallas: Bartzanas et al. (2007) afirman que realizan un estudio de independencia de la malla, pero no indican el procedimiento seguido. Teitel et al. (2008) realizan un estudio de sensibilidad y escogen la de menor resolución que no afecta a los resultados, pero no indican el procedimiento. Majdoubi et al. (2009), indican solamente que escogen la malla que permite mayor rapidez de cálculo. Cabe destacar que los tres modelos declaran validados los modelos CFD sin precisar los errores de cálculo. A pesar de ello, Teitel et al. (2008) avisa de que la CFD no siempre predice los ratios de ventilación adecuadamente. Finalmente cabe resaltar del estudio de Bartzanas et al. (2007) que no muestra diferencias significativas entre los diferentes modelos de turbulencia utilizados para describir el flujo. Así pues, su estudio indica que los modelos de turbulencia no ofrecen resultados que diferencien las principales características del flujo y al validar mediante el ratio de ventilación indica que: el modelo RNG k-ɛ mostró el mejor comportamiento en lo que respecta la precisión del cálculo. Asimismo, en el apartado 8 se añade un apartado con trabajos previos del autor de la presente tesis: Trabajos y comunicaciones del autor 6086 Dinámica computacional de fluidos (CFD) Tabla 6 Comparación de modelos CFD aplicados a invernaderos Invernadero (Autor) Objetivo Variables objetivo Ecuaciones de aproximació n Representa ci Geométrica. Dominio computacion al Condicione s de contorno Condicione s iniciales Malla Aproximacione s Convergenci a Conclusione s principales Invernader o tipo arco (160 m 2 ) (Bartzanas, et al., 2007) Invernader o de 4 módulos (960 m 2 ) (Teitel, et al., 2008) Invernader o canario (1 km 2 ) (Majdoubi, et al., 2009) Validar modelo CFD para analizar la sensibilidad de cuatro modelos de turbulencia Validar modelo CFD para estudiar el efecto de la dirección del viento en la ventilación Validar un modelo CFD para estudiar campos de velocidades y temperatura s Velocidad Interior: 15 puntos con un sensor sónico 3D. Exterior: estación meteorológic a (Bartzanas, et al., 2004) Tasa de ventilación y velocidad Interior: 9 posiciones con tres anemómetro s 3D sónico y dos 2D sónicos Exterior: Veleta y cazoleta Temperatura y humedad Interior 20 puntos. Velocidad Exterior: anemómetro cazoleta a 3 metros de la arista superior del invernadero k-ɛ estándar; RNG k-ε; k-ε realizable; RSM k-ε estándar k-ε estándar Sin datos Sin datos Sin datos Entrada: 7,3 H Salida: 5 H Altura: 5 H Laterales: 5 H Entrada: 50 H Salida: 140 H Altura: 5 H Lateral: 90 H Entrada: 5 H Salida: 5 H Altura: 8 H Lateral: Norte 5 H y Sur 6 H Sin datos Condición de contorno exterior: Sin deslizamiento. Techo y lateral: Simetría Entrada: velocity-ínlet Dominio externo: gradiente de presión Lateral: Conservación del flujo Salida: outlet Inlet: Velocity inlet Perfil de viento entrada constante: v = 3 (ms 1 ) I = 5% Introducción de un perfil de viento logarítmico. (Ecuación ajustada: 2 ms 1 a 5 m de altura. Perfil logarítmico Hexaire (1998) Malla 3D híbrida. No tiene en cuenta el crecimiento. Afirma que realiza un estudio de independenci a de la malla sin especificar cómo No estructurada. Realiza análisis sensibilidad de la malla y escoge la de menor resolución que no afecta a los resultados del modelo. Malla estructurada Después de probar varias mallas se elige la que ofrece mayor rapidez de cálculo ANSYS Fluent ANSYS Fluent ANSYS Fluent Inferior a 1x10 4 ( ) Sin datos Sin datos No existen diferencias significativas entre los resultados al usar diferentes modelos de turbulencia en los patrones del flujo. Cuantitativament e la tasa de ventilación predicha concuerda con datos experimentales. No así las velocidades Valida simulaciones a través de temperaturas 6187 Objetivos 4. Objetivos En la revisión bibliográfica llevada a cabo se observa que la mayoría de los modelos CFD en invernaderos no han estado sujetos a evaluaciones rigurosas. Aunque no es exclusivo del campo de los invernaderos, la aplicación del mismo modelo por diferentes autores tiende a dar resultados diferentes. Además, en los procedimientos de validación que se encuentran en la literatura, a menudo se realiza una mera comparación gráfica entre los resultados CFD y los datos experimentales, lo que no permite estimar el error o cuantificar las incertidumbres debido a condiciones de contorno o a parámetros introducidos incorrectamente. Es importante además, señalar que los resultados individuales de los modelos pueden ser validados, pero no los códigos, o programas. Desgraciadamente, es frecuente encontrar en la literatura afirmaciones consideran un código, o programa, validado únicamente por un experimento. En resumen, los modelos CFD para invernaderos generalmente son creados a partir de muy poca información, debido al alto coste de los datos experimentales, lo que dificulta su validación y generalización por no contar con la suficiente resolución espacial. Se debe de tenerse en cuenta que un anemómetro sónico 3D vale alrededor de euros, lo que sumado a la dificultad de medir en campo, dificulta la toma de medidas en muchos puntos al mismo tiempo. Dado que en los invernaderos generalmente se estudia el flujo atmosférico medio, es decir, sin tener en cuenta las rachas, en rangos de velocidad que van de los 2 a 5 m/s y atendiendo sólo a las 8 direcciones 6288 Objetivos principales de la rosa de los vientos, puede ser interesante utilizar otro tipo de tecnología que ofrezca la posibilidad de aumentar la resolución espacial y temporal de los datos a menor coste. Además, resaltar que no se ha encontrado ninguna guía de buenas prácticas específica para realizar simulaciones CFD en invernadero. Por estos motivos, el objetivo principal de la presente tesis doctoral es proponer un método para generar y validar un modelo de ventilación de un invernadero utilizando la mecánica de fluidos computacional (CFD). Para ello se ha seleccionado un invernadero comercial como ejemplo, sobre el que se han desarrollado diferentes modelos y se expone el procedimiento sistemático que se ha seguido para la selección y validación del modelo final. A lo largo de la tesis se plantearon los siguientes objetivos específicos: Desarrollar y calibrar un sistema automático de medida de la velocidad y dirección del viento de bajo coste y robusto que permita la adquisición de datos en varios puntos simultáneamente. Seguidamente, seleccionar los puntos en los que se deben adquirir dichos datos alrededor del invernadero y realizar las campañas de adquisición. Proponer una sistematización del trabajo de simulación utilizando las guías de buenas prácticas publicadas en otros campos de la tecnología. 6389 Objetivos Generar modelos CFD siguiendo la guía de buenas prácticas propuesta, haciendo especial hincapié en la malla utilizada y el modelo de turbulencia elegido. Validar los resultados experimentales frente a los resultados de las simulaciones realizadas. 6490 Materiales y métodos 5. Materiales y métodos En este capítulo, primero se describen los equipos de medida desarrollados para la adquisición campo y su calibración. En segundo lugar, se presenta el invernadero elegido y se exponen los procedimientos seguidos para la toma de datos experimentales. En tercer lugar, se detalla la generación de los diferentes modelos computacionales para estudiar el flujo ambiental alrededor del invernadero y el flujo de ventilación en el interior. Finalmente, se explican los procedimientos y herramientas estadísticas para la evaluación de los diferentes modelos y su validación El sistema automático de adquisición simultánea de datos. Se ha trabajado con un sistema de adquisición de datos, desarrollado por el Centro de Agroingeniería del Instituto Valenciano de Investigaciones Agrarias. El dispositivo consta de 16 sensores y un sistema de registro de datos. Se complementa con anemómetros convencionales y el conjunto permite medir simultáneamente en 20 puntos la velocidad y dirección de viento. A continuación se describe el sistema y los ensayos realizados para su calibración Anemómetros sónicos El Centro de Agroingeniería dispone de cuatro anemómetros sónicos de la casa comercial Gill Instruments Ltd: dos anemómetros tridimensionales 6591 Materiales y métodos (modelo Windmaster) y dos anemómetros sónicos bidimensionales (modelo Windsonic) cuyas características técnicas principales aparecen en la Tabla 7. Como se verá más adelante, estos anemómetros fueron utilizados como patrones para la calibración del resto del sistema. Estos anemómetros están formados por parejas de transductores opuestos que emiten y reciben ondas ultrasónicas. La medida se basa en que velocidad de propagación de las ondas entre el emisor y el receptor se modifica por la componente de la velocidad de viento en esa dirección. El anemómetro 3D tiene un coste aproximado de frente a de los 2D. Un ordenador se registra directamente las medidas con una frecuencia de 1 Hz a través puerto serie mediante un programa desarrollado por la casa comercial. Tabla 7 Características técnicas de los anemómetros sónicos Windmaster (3D) Windsonic (2D) Anemómetro sónico Velocidad del viento Dirección del viento Rango de medida 0-45 m/s 0-60 m/s Resolución 0,01 m/s 0,01 m/s Precisión Rango de medida <1,5% error medio cuadrático a 12 m/s ±2% a 12 m/s 0-359º 0-359º Resolución 0,1º 1º Precisión 2º a 12 m/s ±3% a 12 m/s Clase de protección IP 65 IP 66 Frecuencia de medida interna 20 o 32 Hz Sin información Frecuencia de salida 1, 2, 4, 8, 10, 16, 20 y 32 Hz 0,25, 0,5, 1, 2 4 Hz Comunicación RS232 RS232 Formato de salida ASCII ASCII 6692 Materiales y métodos Anemómetros mecánicos: sistema cazoleta-veleta En el trabajo se empleó una estación meteorológica comercial de la empresa PCE Ibérica S.L, modelo PCE-FWS 20, que mediante una veleta y un anemómetro de cazoletas ofrece un dato de velocidad y dirección del viento cada 5 minutos (Figura 6). Las cazoletas disponen de dos imanes en su rotor que permiten contar el número de vueltas de las cazoletas por minuto. Este número es proporcional a la velocidad del viento perpendicular al eje del rotor. La veleta tiene instaladas 8 resistencias diferentes que van desde 1kΩ hasta 120kΩ sobre su soporte. Las resistencias son conmutadas por un interruptor que se encuentra en la parte móvil, que gira libremente para alinearse con la dirección del viento. El fabricante no ofrece datos técnicos sobre la precisión resolución o rango de medida. El conjunto dispone de una protección IPX3 y un coste aproximado de 100 euros. Figura 6 Estación meteorológica PCE-FWS 20 utilizada para la toma de medidas en campo. 6793 Materiales y métodos Registrador de datos de los anemómetros mecánicos Para los conjuntos cazoleta-veleta se utilizó un desarrollo propio que adquiere y almacena los datos de velocidad y dirección del viento y se conecta a un ordenador mediante puerto USB. Este circuito está formado por tres bloques. El primer bloque es el encargado recoger los datos de las veletas mediante un microprocesador. Estos datos posteriormente se envían a un micro maestro a través de un bus I2C. El segundo bloque obtiene la señal de los anemómetros de cazoleta y tiene dos microcontroladores. El primero realiza el acondicionamiento de la señal, ya que los anemómetros generan dos pulsos no simétricos por vuelta. El segundo microcontrolador lee los pulsos del acondicionador y los envía al micro maestro por I2C. Éste último solicita los datos de los sensores procesados por los bloques anteriores. La otra tarea que se realiza en este bloque es realizar la comunicación con un ordenador de control a través de USB. Los datos de los sistemas mecánicos se reciben como número de pulsos por segundo y cuadrante de la veleta Y, los datos del registrador se reciben en un ordenador gestionado por un programa de desarrollo propio (Figura 7). Los datos de los anemómetros sónicos se reciben en coordenadas polares, es decir, ángulo de la dirección del viento respecto al Norte y componente horizontal de la velocidad (Figura 8). En total el sistema permite medir simultáneamente la dirección y velocidad del viento en 20 puntos. Por un lado, los datos de 16 conjuntos cazoleta-veleta que provienen de dos registradores de datos de 8 sensores cada uno y, por otro, 68 94 Materiales y métodos los que provienen de los 4 anemómetros sónicos, conectado cada uno un ordenador (Figura 9). Figura 7 Ejemplo de salida de datos de los anemómetros mecánicos Figura 8 Ejemplo de salida de datos de los anemómetros sónicos 3D Figura 9 Conexión del sistema de medida desarrollado por el IVIA 6995 Materiales y métodos Calibración de los anemómetros mecánicos Esta calibración consistió en comparar las medidas de un anemómetro sónico calibrado con las medidas de los anemómetros cazoleta-veleta. Para ello éstos se colocaron junto al anemómetro sónico en una azotea a 15 metros de altura, con el fin de evitar turbulencias del aire debidas a obstáculos cercanos. Los anemómetros se dispusieron a una distancia de 1,5 m entre ellos y se realizaron medidas cada segundo durante un periodo de 18 horas en dos días diferentes (Figura 10). Figura 10 Calibración de los sensores de cazoleta y veleta frente al sensor sónico Se realizó un análisis por separado de la velocidad y de la dirección del viento. En primer lugar, se filtraron los valores del anemómetro de cazoleta inferiores a 2,5 (pulsos/min)x200 ya que se observó que producen distorsiones en la medida debido al comportamiento errático de las cazoletas y de las veletas al iniciarse el movimiento. En segundo lugar, se acondicionó la señal aplicando un filtro de mediana cada 9 datos, es decir, cada 9 segundos. El filtro de mediana 7096 Materiales y métodos permite eliminar el efecto de la inercia que se observa en el anemómetro mecánico y no en el anemómetro sónico, puesto que no es capaz de detectar rápidas variaciones del vector velocidad. En tercer lugar, con la intención de medir el flujo medio, que es el que se empleará en las simulaciones CFD, se promediaron los datos cada 5 minutos, de forma que cada valor de salida corresponde al valor central de la media de cinco minutos, lo que es una frecuencia común para la descripción del viento en simulaciones en invernaderos (Teitel, et al., 2008). Para comparar las medidas de dirección del viento, la señal correspondiente del anemómetro sónico se convirtió a octantes. Al igual que con la velocidad, se promediaron los valores de dirección cada cinco minutos y se realizó un redondeo por truncamiento Descripción del invernadero El invernadero objeto de estudio se encuentra en Catadau (Valencia), pertenece a la empresa Mediterránea de Cactus y está dedicado al cultivo de cactus de pequeño tamaño durante todo el año. Se trata de un invernadero de m 2 de tipo multicapilla 9 capillas de 8 m de luz y 5,5 m de altura y cubierta de plástico, típico en la zona (Figura 11 y Figura 12). El invernadero está en campo abierto, delimitando al Norte y al Este con cultivos de naranjos; al sur con una explanada sin ningún tipo de cultivo ni obstáculos; y al Oeste con una explanada donde se encuentra una nave de pequeño tamaño. Respecto a su orientación, el eje longitudinal del invernadero 7197 Materiales y métodos está en dirección Sureste-Noroeste, que coincide con la dirección de los vientos dominantes en esa zona. Este invernadero tiene dos formas de ventilación natural: a través de las aperturas laterales y a través de las ventanas cenitales, como ocurre en la mayoría de los invernaderos actuales (Anejo E). 7298 Materiales y métodos Figura 11 Emplazamiento invernadero. Catadau (Valencia) Figura 12 Detalle del invernadero y sus alrededores 7399 Materiales y métodos 5.3. Medidas de campo Para estudiar el flujo en el exterior del invernadero se colocaron los sensores mecánicos en 32 puntos, utilizando para ello los 16 sensores colocados alternativamente a 4,5 m y a 3 m sobre mástiles metálicos. Para medir la velocidad de aproximación del viento se situó un anemómetro sónico tridimensional en una zona alejada del invernadero, libre de obstáculos cercanos y a 10 metros de altura. En dos ventanas cenitales del invernadero se tomaron 4 medidas en la cara Sureste y Noroeste, utilizando para ello alternativamente un sensor sónico tridimensional y un sensor sónico bidimensional (Figura 13). Además, se tomaron los datos de una estación meteorológica cercana, situada en Carlet (Valencia), que forma parte de la Red SIAR (Modelo 1). Se trata de las primeras estaciones agroclimáticas de la Comunidad Valenciana que fueron instaladas en 1999 y sus características se recogen en la Tabla 8. Tabla 8 Estación meteorológica Red SIAR (Modelo 1) Instrumento Especificaciones Datalogger Datalogger DL15 de Thies Clima, modelo /.001 Sensor de Temperatura-Humedad Sensor de Radiación Transmisor de humedad y temperatura-compact de Thies Clima, modelo Piranómetro de Kipp & Zonen modelo CMP3 Sensor de velocidad de viento Anemómetro-compact de Thies Clima, modelo Sensor de Dirección de Viento: Anemoscopio-compact de Thies Clima, modelo Pluviómetro Transmisor de precipitación de Thies Clima, modelo100 Materiales y métodos Figura 13 Puesta en campo de los sensores cazoleta, veleta y sónicos. 75101 Materiales y métodos 5.4. Simulaciones Variables objetivo Las variables objetivos fueron la velocidad y la dirección horizontal del viento en 40 puntos: 4 puntos en las ventanas de la fachada a barlovento, 4 puntos en la fachada a sotavento y 32 puntos repartidos en la esquina Sureste- Noroeste del invernadero a 3 metros y 4,5 metros de altura Ecuaciones de aproximación Se utilizaron los siguientes modelos: Standard k-ɛ, k-ɛ RNG, SST K-ω y RSM en esquemas de primer y segundo orden Representación geométrica de los obstáculos Se eligió un invernadero aislado para no que no estuviese afectado por la presencia de otros invernaderos o edificios cercanos. Existe una pequeña nave que se encuentra a una distancia superior de 6H, a 30 m de la zona donde se toman las medidas. Asimismo, no se encuentra a barlovento de la zona de medida y no se espera que su estela pueda afectar a los datos experimentales. No obstante, se realizó una simulación preliminar simple donde se vio que la presencia o su ausencia no modifican los resultados de la simulación en la zona de interés Dominio computacional Para la malla exterior del invernadero se utilizaron extensiones verticales y laterales de 6,5 H. Ambas son superiores a las recomendaciones mínimas que 76102 Materiales y métodos se describen en el apartado y dan como resultado un ratio de obstrucción del 9,7% inferior a la recomendación del 10% de VDI. No obstante, el ratio de obstrucción es superior a la recomendaciones de (Tominaga, et al., 2008) de un máximo del 3% y no se respeta que el ratio del edificio para edificios más largos que anchos. Mantener la primera recomendación hubiera resultado en un dominio de 13 H de alto y 18 H de ancho, lo que sobrepasa la capacidad computacional disponible, y supondría aumentar la superficie transversal simulada de m 2 a m 2. La segunda recomendación hubiese supuesto disminuir el dominio y no mantener las recomendaciones COST732. De todos modos, se realizó una simulación con una malla de un millón de celdas con geometrías simples para comparar los resultados y no se obtuvieron diferencias en las zonas estudiadas. En la zona de aproximación del flujo se mantuvieron todas las recomendaciones de la literatura, con una extensión de 5H a barlovento y 15H a sotavento del invernadero Perfil del viento de entrada Se utilizó el perfil de viento propuesto por Ayra (1988) donde la velocidad de viento (u z ) es una función tipo logarítmico donde se establece una velocidad de referencia (u ref ) a una altura de referencia medida en la estación de enlace, Ecuación (5.1) 77103 Materiales y métodos ln ( z z ) u z = u ref [ 0 ln ( z ] (5.1) ref z ) 0 La intensidad turbulenta se supuso como un 5% de la velocidad de entrada, valor típico para el aire atmosférico Modelo de suelo Se simuló el suelo alrededor del invernadero utilizando la rugosidad ks de las tablas de Davenport. En la cara Sureste y en la cara Noroeste existe una explanada de hormigón. En la cara Suroeste, zona de mayor importancia al tratarse del área de aproximación del viento exterior, existe un solar sin cultivar con presencia de malas hierbas espaciadas. Finalmente, en la cara Noreste, existe una plantación de naranjos con un marco de plantación muy denso a un nivel inferior, de modo que la parte superior de las copas coincide con el nivel del invernadero. Tabla 9 Valores z 0 elegidos para la rugosidad del suelo alrededor del invernadero Cara Características Parámetro z 0 Noreste Sureste Suroeste Noroeste Plantación de naranjos escalonado, las copas de los árboles corresponden al nivel del invernadero Explanada sin cultivar Solar sin cultivar Explanada sin cultivar Rugoso abierto 0.10 Entre terreno suave y abierto Entre terreno suave y abierto Entre terreno suave y abierto Suave104 Materiales y métodos Condiciones de contorno El lateral Sureste se utilizó como condición de entrada la velocidad de entrada, conocida en la literatura como velocity-inlet, mediante un perfil de aire logarítmico. En los laterales Noreste y Suroeste se utilizó la condición de simetría. Finalmente, en el lateral Noroeste se utilizó la condición salida de flujo, conocida en la literatura como outflow Condiciones iniciales Para las simulaciones de con un esquema de cálculo de primer orden se utilizó como condición de contorno el perfil de entrada logarítmico especificado a la entrada del dominio. En cambio, para las simulaciones de segundo orden se utilizaron como condiciones iniciales los resultados de simulaciones anteriores y tan sólo se modificó el esquema de cálculo. 79105 Materiales y métodos Malla Con el objetivo de generar una malla que represente adecuadamente el invernadero y los alrededores se realizó un levantamiento topográfico mediante estación total de la zona, del invernadero y de los obstáculos cercanos ya que no se disponían de planos. En segundo lugar, se definió todo el espacio computacional siguiendo las recomendaciones del apartado y las cotas obtenidas en el levantamiento topográfico. En tercer lugar, se definieron las superficies que representan las superficies del invernadero y las condiciones de contorno. En cuarto lugar, se definió la parte del dominio por donde circula el fluido (cuerpo o body) y por tanto en la que se desea generar las celdas. En quinto lugar, para un mallado más sencillo, el cuerpo se dividió en bloques siguiendo dos técnicas: Split (bloques cartesianos) y 0-grid (malla circular o de tipo 0 ). La primera produce celdas más regulares y es preferible cuando se mallan superficies rectangulares; la segunda se utiliza generalmente para definir con mayor resolución una característica concreta de la geometría. Así pues, se confeccionaron dos mallas, una con cada técnica, con la intención de seleccionar la que produjo celdas más regulares en las esquinas de los invernaderos (Figura 14). 80 106 Materiales y métodos Figura 14 Planta del dominio del invernadero dividido mediante: (arriba) la técnica Split y (abajo) la técnica O-grid. Representación de la forma de las celdas en las esquinas: o-grid (a y b) y la técnica split (c y d). Finalmente, se definieron las celdas en el interior de estos bloques. Siguiendo las recomendaciones del apartado se debe de comenzar al elegir un tamaño mínimo de celda menor de un décimo la altura del invernadero y que permita tener un mínimo de 10 celdas por raíz cúbica. Para ello se generó una primera malla tridimensional estructurada a partir de una celda de 45 cm de lado. Estas celdas se situaron alrededor del invernadero, puesto que es la zona 81107 Materiales y métodos de más interés y necesita tener mayor resolución. Para elaborar la malla se dividió el dominio en un bloque central para el invernadero, con un tamaño uniforme de celdas, y en otros bloques en los que se estableció el criterio de que no hubiese un crecimiento superior al 5% entre dos celdas consecutivas en la dirección exterior al dominio. Este criterio es mucho más restrictivo que los explicados en la revisión bibliográfica, pero se mantuvo porque un estudio de sensibilidad demostró que los ratios superiores al 5% modificaban el resultado de las simulaciones en la zona de interés. Para la selección definitiva de la malla, se generaron otras tres mallas, siguiendo la extrapolación de Richardson, con tamaños de celda mínimos de 45, 30 y 20 cm; manteniendo el resto de criterios constantes. Posteriormente, se realizó una simulación con cada malla para un perfil de velocidad constante de 2 m/s en la dirección longitudinal del flujo, con un modelo k-ɛ estándar y un esquema de discretización de primer orden. Respecto al tamaño de celda cercano al suelo se intentó mantener el criterio de y p > k s para C s = 1. Cuando esto no es posible, ejemplo del tamaño mínimo de celda de 20 cm, es necesario recurrir a las técnicas explicadas en el apartado Se estableció como criterio de generación de la malla que todas las celdas tuvieran una calidad superior a 0.9(-) con el parámetro de calidad determinante 3x3x Aproximaciones numéricas Para todas las simulaciones se utilizaron esquemas de discretización de primer y de segundo orden. 82108 Materiales y métodos El algoritmo de resolución utilizado fue el algoritmo SIMPLE (Anejo B) Criterios de convergencia Se consideró que una simulación convergía cuando los residuos de todas las variables fueron inferiores a 10 4 órdenes de magnitud y la pendiente de los residuos era cercana a Simulaciones ensayadas y validación Se simularon 24 casos en total, resultado de la combinación de cuatro modelos de turbulencia (k-ε; RNG k-ε; k-ϖ y RSM), dos esquemas de discretización (primer orden y segundo orden) y tres velocidades de viento exterior (3; 3,5 y 4 m/s). La Tabla 10 resume los principales parámetros de la simulación. Para realizar la validación se compararon los datos obtenidos con las simulaciones con los datos obtenidos experimentalmente mediante análisis de regresión lineal (Romero & Zúnica, 1993) para perfiles de viento de 3; 3,5 y 4 m/s. 83109 Materiales y métodos Tabla 10 Cuadro resumen de los parámetros que se utilizaron en las simulaciones Apartado Variables objetivo Ecuaciones de aproximación Representación geométrica de los obstáculo Dominio computacional Perfil del viento de entrada Modelo de suelo Condiciones de contorno Condiciones iniciales Malla Aproximaciones numéricas Residuos Elección de parámetros Velocidad y dirección del viento medidos en 40 puntos (simultáneamente 20 puntos) RANS modelos k-ɛ estándar, k-ɛ RNG, SST K-ω y RSM No se ha considerado ningún obstáculo Extensión vertical del dominio: 6,5 H Extensión aproximación del flujo al edificio: 5 H Extensión estela del edificio: 15H Extensión lateral del dominio:5h Perfil logarítmico definiendo la velocidad medida a 10 metros de altura Clasificación de Davenport para la rugosidad del suelo Cara Noreste z 0 : Rugoso abierto 0,10 Cara Noreste z 0 : Entre terreno suave y abierto 0,015 Cara Noreste z 0 : Entre terreno suave y abierto 0,015 Cara Noreste z 0 : Entre terreno suave y abierto 0,015 Laterales Sureste: velocity inlet Lateral Noreste y Suroeste: Simetría Lateral Noroeste: outflow Superficies invernadero: modelo de pared estándar sin rugosidad Simulaciones de primer orden: Perfil de entrada del dominio Simulaciones de segundo orden: Resultado de las simulaciones de primer orden Malla estructurada tridimensionales formadas por hexaedros Crecimiento del 5% entre dos celdas consecutivas Calidad de la malla superior a 0,8(-) Celdas próximas al suelo y p > k s para C s = 1 Extrapolación de Richardson para la elección de la malla SIMPLE Discretización de primer y segundo orden Residuos estabilizados en un valor inferior a para todas las ecuaciones 84110 Resultados y discusión 6. Resultados y discusión Los resultados de esta tesis se presentan en cuatro bloques. En primer lugar, se muestran los resultados de la calibración de los anemómetros de cazoleta y veleta frente al anemómetro sónico 3D utilizado como patrón, demostrando que el sistema de adquisición de datos puede ser utilizado para el estudio. En segundo lugar, se describen los resultados obtenidos alrededor del invernadero y en las ventanas en la campaña de datos experimentales. En tercer lugar, se muestran los resultados de las 24 simulaciones CFD propuestas para estudiar el flujo ambiental alrededor del invernadero junto con la validación de las simulaciones CFD frente a los datos experimentales Calibración del sistema de medida En primer lugar se muestra el gráfico de dispersión de los valores de velocidad de viento de todos los anemómetros de cazoleta frente al anemómetro sónico (Figura 15). La regresión es estadísticamente significativa (F=28,43; p- valor<0,01) y su ecuación es: y ( m ) = 1,077x + 0,351. s El coeficiente de determinación fue R 2 = 0,98 con error absoluto medio de 0,06 m/s. Además, el error estándar de la estimación muestra una desviación típica de los residuos de 0,73 m/s. Este valor podrá ser utilizado para construir límites de predicciones para nuevas observaciones. La recta de calibración es 85111 Resultados y discusión válida para el rango de 2 a 4,5 m/s y cuando las cazoletas están detenidas se supone una velocidad próxima a 0 m/s. En lo que respecta a la dirección del viento, tras analizar más de ocho horas de datos de todas las veletas, las direcciones marcadas por el anemómetro de cazoleta frente al anemómetro sónico coincidieron en el 100% de los casos (Figura 16). Estos datos demuestran que el sistema permite estudiar el vector velocidad media en dos dimensiones en el intervalo de 2 a 4,5 m/s, además de indicar cuando la velocidad del viento es próxima a 0 m/s, con una precisión similar al anemómetro sónico 3D con un conste considerablemente menor, lo que nos permite aumentar el número de medidas simultáneas. El sistema cumple con los requisitos para la adquisición de datos de viento para la generación y validación de modelos CFD en el exterior de los invernaderos. 86112 Cuadrantes Anemómetro sónico (m/s) Resultados y discusión 5 4,5 4 3,5 3 2,5 2 1,5 1 0, Anemómetro Cazoleta (pulsos/min)x200 Figura 15 Comparación entre los módulos de velocidad los anemómetros de cazoleta y sónico. Ambos anemómetros miden con una precisión similar la velocidad del viento Anemómetro de cazoleta Anemómetro sónico Tiempo = 10 h Figura 16 Comparación entre la dirección medida por el anemómetros de cazoleta y sónico. Ambos anemómetros miden con una precisión similar la dirección del viento. 87113 Resultados y discusión 6.2. Análisis de la malla del invernadero Se realizaron tres mallas con diferentes tamaños mínimos de celda: 45 cm de lado, 30 cm y 20 cm. El flujo en las tres mallas fue resuelto utilizando un modelo k-ɛ estándar con un perfil constante de viento de 2 m/s. Para la primera malla, tamaño mínimo de 45 cm, la simulación no convergió. Entre la segunda y tercera malla al comparar los resultados del vector velocidad en la zona de interés no se apreció ninguna diferencia entre ellas. Así pues, se seleccionó la malla con un tamaño mínimo de celda de 30 cm. Esta malla es una malla estructurada, formada por alrededor de 6,1 millones de celdas hexaédricas representando un volumen de más de m 3 con una calidad superior a 0,95(-) en todas las celdas (Tabla 11). Las menores celdas se encuentran en la zona donde se tomaron las medidas experimentales y alrededor del invernadero, que es donde se produce la separación del flujo y los fenómenos son más complejos. Así mismo, la calidad de la malla es superior a 0,99(-) exceptuando 24 celdas, con una calidad superior a 0,95(-), contiguas al vano más pequeño debido al menor espacio disponible para adaptar la forma del invernadero en ese vano. El crecimiento entre celdas consecutivas fue de un máximo del 5%. Esta malla tiene un número de celdas mucho más elevado que lo que se encuentra en la literatura para estudios de invernaderos (alrededor de un millón de celdas). Entre las dos técnicas de generación de la malla, split y 0-grid, se selección la primera técnica ya que se obtienen celdas más regulares en la zona de interés. 88114 Resultados y discusión 1 2 Tabla 11 Malla del ambiente exterior del invernadero Malla Celdas Figuras -66 nodos Dominio a -Menor nodo 30 cm sotavento -Crecimiento 5% (88 m) -Calidad >0.99(-) Dominio barlovento (30 m) -43 nodos -Menor nodo 30 cm -Crecimiento 5% -Calidad >0.99(-) 3 Dominio lateral (30 m) -43 nodos -Menor nodo 30cm -Crecimiento 5% -Calidad >0.99(-) 4 (a) Fachada Sureste y Noroeste (68 y 71 m) -8 vanos x 29 nodos -1 vano x 15 nodos (Suroeste) -Menor nodo 30 cm -Crecimiento 5% -Calidad >0.99(-) 5 (b) Fachada Noreste y Suroeste (95 y 82 m) -319 nodos -Menor nodo 30 cm -Sin crecimiento -Calidad >0.99(-) 6 Altura invernadero (5.5 m) -22 nodos -Menor nodo 30 cm -Crecimiento 5% -Calidad >0.99(-) excepto en vano Suroeste 24 celdas calidad >0.95(-) Esquina superior (c) Canalón (d) Esquina inferior (e) Celdas calidad inferior (f) 89115 Resultados y discusión 6.3. Medidas experimentales de velocidad y dirección del viento Como ejemplo se muestran los resultados experimentales alrededor del invernadero para un viento de dirección Sureste-Noroeste con una velocidad media de 4 m/s. En la zona de aproximación (sensores B) los sensores indican una velocidad promedio aproximadamente un 10% inferior a la velocidad promedio del viento. Sin embargo, para las velocidades de 3 y 3,5 m/s no se observó prácticamente reducción de la velocidad del viento. Además, a pesar de que la velocidad media a 4,5 m de altura es mayor que a 3 m, típico de un perfil logarítmico, estas diferencias no fueron estadísticamente significativas. Esto puede ser debido a la poca rugosidad del terreno a barlovento, donde se desarrolla el perfil de viento que afecta al invernadero. En la esquina Noreste (sensor A1) se observa una aceleración del flujo, ya que se aprecia una velocidad del aire un 30% superior a la velocidad de aproximación, así como un cambio en la dirección del flujo. Además, en las posiciones de los sensores A2 y A3 se midieron velocidades cercanas a 0 m/s y se observó una circulación en sentido contrario al del perfil de aproximación del viento. Más a sotavento (sensores A4 a A8) se recupera el módulo de la velocidad, pero se modifica la dirección del aire respecto al viento de aproximación. En los sensores colocados en las ventanas a barlovento se observa una aceleración del aire de hasta un 40%, sin variar la dirección del flujo de aproximación a la fachada y a sotavento se observa una deceleración. 90116 Resultados y discusión Todo ello se corresponde con los fenómenos de separación de flujo al chocar el aire con un edificio, explicados en el apartado 2.3. En la esquina Sureste-Noreste, se produce la separación del flujo, y por tanto una aceleración que aguas abajo produce recirculaciones del flujo en forma de torbellinos en la fachada Noroeste, uniéndose el flujo de nuevo más a sotavento y recuperándose de nuevo el perfil del viento (Tabla 12). Todos los datos obtenidos se resumen en el Anejo C. Tabla 12 Vectores velocidad alrededor del invernadero (velocidad exterior 4 m/s) Sensor Velocidad Dirección Velocidad Dirección Sensor (m/s) (cuadrante) (m/s) (Cuadrante) A1 4,5 m 5,2 2 4,5 m 3,48 1 B1 3 m 5, m 3,72 1 A2 4,5 m 0 5 4,5 m 2,64 1 B2 3 m m 3,6 1 A3 4,5 m 2,84 8 4,5 m 2,6 1 B3 3 m m 2,52 1 A4 4,5 m 3,44 8 4,5 m 3,68 1 B4 3 m m 3,2 1 A5 4,5 m 3,76 8 4,5 m 3,76 1 B5 3 m 3,4 8 3 m 3,08 1 A6 4,5 m 4,01 8 4,5 m 4,36 1 B6 3 m 4, m 3,52 1 A7 4,5 m 3,32 8 4,5 m 4,32 1 B7 3 m 3, m 4,52 1 A8 4,5 m 4,72 8 4,5 m 3,28 1 B8 3 m 3,6 8 3 m 3,117 Resultados y discusión 6.4. Resultados de las simulaciones y validación Se simularon en total 24 casos: cuatro modelos de turbulencia (k-ε; RNG k-ε; k-ϖ y RSM); para dos esquemas de discretización (primer orden y segundo orden) para una única dirección de viento (Sureste-Noroeste) con tres velocidades de viento exterior (3; 3,5 y 4 m/s). Además se midió la velocidad y dirección del viento en 40 puntos alrededor del invernadero: 16 en la fachada suroeste donde se produce la separación del flujo (a 3 y a 4,5 m de altura); 16 en la explanada de aproximación (a 3 y a 4,5 m de altura); 4 en la fachada de barlovento en dos ventanas cenitales y, finalmente, otros cuatro, en posiciones análogas en la fachada de sotavento. Puesto que se observó que modificar la velocidad exterior del viento no modificó lo resultados de las estructuras del flujo, sólo modificó el módulo de los vectores velocidad proporcionalmente a la velocidad exterior, se decidió representar las velocidades estandarizadas respecto a la velocidad exterior. En la zona de separación del flujo, que se corresponde con el lugar donde se colocaron los sensores A, los diferentes modelos de turbulencia no predicen de la misma forma el flujo de aire en la fachada sudoeste, ni en el módulo, ni en la dirección o sentido. Atendiendo únicamente al módulo de la velocidad media, a 3 m de altura todos los modelos predicen de igual manera una aceleración del flujo de un 30% en A1 respecto a la velocidad exterior. Más a sotavento se observa una deceleración del flujo, que es diferente 92118 Resultados y discusión según el modelo. El modelo que produce resultados más diferentes respecto a los demás y a los datos experimentales es el RNG k-ε, con una deceleración máxima en el punto A2 de un 60%, y con aceleraciones en los puntos A3 y A4. En el resto de simulaciones la aceleración no se produce hasta las posiciones de los sensores A4 y A5. Sin embargo, mientras que los modelos SST k-ω y RSM la velocidad aumenta rápidamente hasta un % de la velocidad exterior, en el resto de modelos la velocidad aumenta más lentamente y en menor medida: entre un 60-70% de la velocidad exterior. En segundo lugar, en las posiciones a 4.5 m, el modelo k-ε estándar de primer orden y de segundo orden captan una deceleración, que es un 30% menor en A2, siendo la deceleración máxima para el modelo k-ε estándar de primer orden en A3. Además, estos dos modelos predicen velocidades a sotavento un 20% menor que los demás modelos (Figura 17). Respecto a la dirección del viento en la fachada Noreste, todos los modelos predicen una recirculación del flujo del invernadero aunque en posiciones diferentes. Si se dibujan las líneas de corriente en un plano perpendicular al suelo que contenga los puntos A, es decir paralelo a la fachada a un metro y medio de distancia de ésta, se observan las diferencias en la predicción de las recirculaciones. El modelo k-ε estándar predice un vórtice aproximadamente la mitad de grande que el vórtice predicho por el modelo RSM, y éste además, está desplazado alrededor de 4,2 metros respecto al modelo RSM (Figura 19). Respecto a las posiciones de los sensores colocados en las fachadas a 93119 Resultados y discusión barlovento y a sotavento los diferentes modelos de turbulencia predijeron con total similitud el flujo de separación en la fachada de barlovento. En cambio, en la fachada a sotavento en tres puntos hubo coincidencia entre las diferentes predicciones mientras que en una de las dos posiciones más elevadas los modelos de turbulencia predijeron los mismos resultados (Figura 18). Así pues, En la fachada a barlovento los diferentes modelos de turbulencia predicen aceleraciones de un 40% en las dos posiciones más elevadas y deceleraciones de un 60% en las otras dos. En la fachada a sotavento, en una de las posiciones elevadas los modelos de turbulencia sí difieren en su predicción: el modelo k-ε estándar resulta en una desaceleración del 30% del flujo mientras que el modelo RNG k-ε predice una aceleración de un 10% de la velocidad exterior. Los modelos SST k-ω y RSM en este punto predicen la misma desaceleración de un 10% respecto la velocidad exterior. Respecto a la dirección del viento en estas posiciones todos los modelos predicen de igual forma. Respecto a la zona de aproximación, equivalente a donde se situaron los sensores B los diferentes modelos de turbulencia produjeron el mismo resultado, tanto en módulo como en la dirección de viento, en los puntos equivalentes a donde se tomaron las medidas experimentales. Por todo lo anterior, se puede decir que los modelos se pueden agrupar en 3 tipos, en función de la similitud de sus predicciones respecto al módulo de la velocidad de viento: por un lado el modelo RNG k-ε, por otro lado los modelos k-ε estándar de primer orden y de segundo orden, y finalmente el modelo SST k- ω y el modelo RNG. 94120 Velocidad estandarizada Velocidad estandarizada Resultados y discusión 1,4 1,2 1 0,8 0,6 0,4 0,2 0 A1 (3 m) A2 (3 m) A3 (3 m) A4 (3 m) A5 (3 m) A6 (3 m) A7 (3 m) A8 (3 m) Velocidad experimental media k-ɛ 1er orden k-ɛ 2do orden RNG k-ɛ 2do orden RSM 2do orden SST K-ω 2do orden 1,4 1,2 1 0,8 0,6 0,4 0,2 0 A1 (4 m) A2 (4 m) A3 (4 m) A4 (4 m) A5 (4 m) A6 (4 m) A7 (4 m) A8 (4 m) Velocidad experimental media k-ɛ 1er orden k-ɛ 2do orden RNG k-ɛ 2do orden RSM 2do orden SST K-ω 2do orden Figura 17 Velocidad estandarizada en la fachada suroeste del invernadero para los diferentes modelos de turbulencia y datos experimentales (de arriba a abajo): a 3 m y 4 m de altura. 95121 Velocidad estandarizada Resultados y discusión 1,6 1,4 1,2 1 0,8 0,6 0,4 0,2 0 Barlovento 1 Barlovento 2 Barlovento 3 Barlovento 4 Sotavento 1 Sotavento 2 Sotavento 3 Sotavento 4 Velocidad experimental media k-ɛ 1er orden k-ɛ 2do orden RNG k-ɛ 2do orden RSM 2do orden SST K-ω 2do orden Figura 18 Velocidad estandarizada en la fachada Sureste (Barlovento) y Noroeste (Sotavento) para los diferentes modelos de turbulencia y experimentales 96122 Resultados y discusión Figura 19 Representación de las líneas de corriente en un plano perpendicular al suelo contiene los sensores A. Modelo k-ε (arriba) y modelo RSM (abajo). Los cuadros rojos representan los puntos de medida para sensores A. 97123 Resultados y discusión Respecto a la validación, el mejor ajuste entre los datos predichos y los experimentales se produjo para los modelos RSM y SST k-ω para todas las velocidades de viento, con coeficientes de determinación R 2 de hasta 0,87 (modelo RSM, velocidad del viento 4 m/s, error estándar de la estimación 0,55 m/s). Por el contrario, el modelo k-ε estándar con el esquema de cálculo de primer orden fue el que peor se ajustó, como ejemplo para 3 m/s (R 2 menor del 0,5, error estándar de estimación 1,01 m/s). El esquema de cálculo de segundo orden, mejoró la regresión (R 2 =0,75, error de estimación 0,95 m/s) pero lejos de los resultados ofrecidos por SST k-ω y RSM. Todo esto muestra diferencias en los patrones del flujo alrededor del invernadero al utilizar diferentes modelos de turbulencia. Estas diferencias necesariamente impactaran al estudiar el flujo en el interior del invernadero. Esto contrasta con las conclusiones del trabajo realizado por Bartzanas et al. (2007), en las que se indicó que no hay diferencias en los patrones del flujo al utilizar los diferentes modelos de turbulencia. Así pues, el modelo k-ε estándar de primer orden, el más utilizado en la literatura científica para ventilación de invernaderos, no sólo ofrece resultados diferentes al resto de modelos, sino que describe peor los flujos que intervienen en la ventilación. El modelo RNG k-ε dio resultados intermedios, pero no se tomaron en consideración debido a problemas de convergencia. Teóricamente el modelo RNG k-ε debería haber mejorado el rendimiento del modelo k-ε estándar en las zonas de recirculación y separación del flujo, pero tampoco predijo 98124 Resultados y discusión correctamente la energía cinética turbulenta. Además, al hacer un análisis más detallado del flujo, se observó que las líneas de corriente de este modelo difieren mucho de los demás modelos estudiados, prediciendo dos vórtices de diferente tamaño (Figura 22). Esto puede ser debido a que el límite de convergencia (10 4 órdenes de magnitud) no es suficiente. Por este motivo se intentó reducir el límite de convergencia un orden de magnitud adicional para el modelo RNG k- ε, pero resultó en un intento infructuoso. La literatura recoge que el pobre rendimiento del modelo k-ε en las zonas de recirculación o separación del flujo está relacionado con la sobrestimación de la energía cinética turbulenta en los lugares donde se producen estos fenómenos. Así pues, se procedió a estudiar este parámetro en la esquina sureste del invernadero. Efectivamente, la energía cinética turbulenta predicha fue un 50% superior en el modelo k-ε estándar. A modo de ejemplo se representa la energía cinética turbulenta del modelo k-ε estándar frente al modelo RSM (Figura 21). Tanto el modelo RSM como el modelo SST k-ω produjeron buenos resultados, probablemente porque incluyen mejoras en la simulación de los flujos de recirculación, donde los esfuerzos cortantes son importantes. Hay que resaltar que el modelo RSM necesita casi 4 veces más de tiempo de computación que el modelo SST k-ω para converger y los resultados son muy similares. Finalmente hay que resaltar que los ajustes de los modelos con los datos experimentales son mejores a mayor velocidad de viento, esto es probable que 99125 Resultados y discusión se deba a que estos vientos fueron menos racheados y, por tanto los fenómenos transitorios, no tenidos en cuenta en los modelos, son menos importantes. Todos los análisis de regresión se recogen en el Anejo D. Plot of Fitted Model Experimental (3ms) = -0, ,1638*k-eps 1er orden Plot of Fitted Model Experimental (4ms) = -0, ,108*k-eps 1er orden 1,6 1,5 Experimental (3ms) 1,2 0,8 0,4 Experimental (4ms) 1,2 0,9 0,6 0, ,3 0,6 0,9 1,2 1,5 k-eps 1er orden 0 0 0,3 0,6 0,9 1,2 1,5 k-eps 1er orden F-ratio:37,5 p-valor: Coeficiente correlación: 0.70 R 2 = Error estándar de estimación:1,01 m/s Error absoluto medio: 0,84 m/s F-ratio:81,08 p-valor: Coeficiente correlación: R 2 = Error estándar de estimación:0,87 m/ Error absoluto medio: 0,65 m 1,5 Plot of Fitted Model Experimental (4ms) = -0, ,11716*k-eps 2do orden 1,5 Plot of Fitted Model Experimental (4ms) = -0, ,13649*RSM 2do Experimental (4ms) 1,2 0,9 0,6 0,3 Experimental (4ms) 1,2 0,9 0,6 0, ,3 0,6 0,9 1,2 1,5 k-eps 2do orden 0 0 0,3 0,6 0,9 1,2 1,5 RSM 2do F-ratio: p-valor: Coeficiente correlación: 0.86 R 2 = Error estándar de estimación:0,76 m Error absoluto medio: 0,59 m F-ratio: p-valor: Coeficiente correlación: 0,93 R 2 = Error estándar de estimación:0,55 m Error absoluto medio: 0,44 m Figura 20 Regresión lineal de los datos experimentales frente a los datos de simulación para (de izquierda a derecha y de arriba abajo) k-ε estándar primer orden vs 3 m/s; k-ε primer orden vs 4 m/s; k-ε segundo orden 4 m/s; k, RSM segundo orden vs 4 m/s. 100 126 Resultados y discusión Figura 21 Detalle de la energía cinética turbulenta del modelo en la esquina Sureste-Noreste del invernadero k-ε (Izquierda) y RSM (Derecha) para velocidad exterior de 3 m/s Figura 22 Líneas de corriente en un plano longitudinal a 1,5 m de la fachada Noreste del invernadero predichas por el modelo RNG k-ε que contiene los puntos A. 101127 Resultados y discusión 6.5. Descripción del flujo alrededor del invernadero Las simulaciones con diferentes velocidades muestran las mismas estructuras y sólo cambia las velocidades medias de los fenómenos. Efectivamente, una vez validado el flujo de ventilación en la esquina Sureste- Noroeste se puede utilizar los resultados de la CFD para describir el flujo. En este sentido, como ejemplo, si consideramos dos planos horizontales: a 3 y a 4 metros de altura, para el modelo RSM de segundo orden y con una velocidad exterior de 4 m/s (Figura 23). Se observa una desaceleración del flujo en la fachada Sureste (barlovento) y una primera recirculación; seguida de una separación prácticamente simétrica del flujo a ambos lados del invernadero, producida en las aristas que genera una segunda recirculación; una aceleración del flujo a ambos lados del invernadero; y una tercera recirculación tras el invernadero. Al incidir el flujo sobre la fachada a barlovento se produce un gradiente negativo de presión que produce primero la desaceleración del flujo y la primera recirculación y un cambio de dirección que se describirá en el siguiente apartado. Luego una separación del flujo a ambos lados del edificio. Esta separación del flujo genera dos recirculaciones alargadas de aproximadamente 3 veces la altura del invernadero Si nos centramos en la esquina Sureste-Noreste del invernadero, donde se ha validado el flujo, las líneas de corriente muestran una recirculación de vórtice único cuyo centro se encuentra aproximadamente 5 m, es decir, a 1H de la pared donde se produce la separación del flujo, desplazado respecto a la fachada aproximadamente 1 m (0.2H), con un tamaño aproximado de 3x15 m. Si continuamos hacia sotavento, el flujo vuelve a 102128 Resultados y discusión pegarse al edificio a una distancia de 20 m del inicio del inicio de la separación, 4H, donde las líneas de corriente vuelven a ser paralelas al flujo principal. En segundo lugar, se observa una aceleración, también simétrica a ambos lado del invernadero. Se trata del sistema de vórtices de herradura horizontal similar al que se produce cerca del suelo, a barlovento, y se desvanece a sotavento como pares de remolinos girando en sentido opuesto alrededor del edificio. Estos remolinos no han sido captados por la simulación debido a que se trata de una simulación estacionaria. Finalmente, de la misma forma, en la estela, se presentarían los típicos sistemas de vórtices periódicos alternados, que tampoco son captados al tratarse de una simulación estacionaria. En vez de ello se presenta un flujo caótico, con las líneas de corriente representando varios vórtices, de diferentes tamaños, siendo mayor el de la esquina Suroeste- Noroeste. Este vórtice puede ser mayor debido a la mayor longitud del edificio en ese lateral, o más probablemente a su mayor altura. En todo caso, la simulación del flujo en esta zona no ha sido validada. En la sección a 4.5 metros de altura, que coincide con la altura de los sensores más elevados, tanto a barlovento como a sotavento se presentan las mismas estructuras comentadas anteriormente, aunque de menor tamaño, ya que nos encontramos a mayor altura respecto del suelo y del techo del edificio, lo que se observa mejor en lateral derecho a sotavento del invernadero. En la sección longitudinal del dominio se observa mejor la recirculación aguas arriba del invernadero. Un segundo vórtice de gran tamaño se genera en el techo, de 5x20 m (Hx4H), al despegarse el flujo debido a la fachada. Este 103129 Resultados y discusión flujo, al contrario de lo que se describe en el apartado 2.3, no vuelve a converger aguas abajo del invernadero, sino sobre el techo, porque el edificio es más largo que ancho. El flujo vuelve a converger 40 m aguas abajo del vértice de la fachada a barlovento, es decir un poco antes de la mitad de la longitud del invernadero. Esto es importante desde el punto de vista de la ventilación, debido a que por las presiones negativas generadas en esta zona, el aire del interior del invernadero tiende a salir. Finalmente, en la estela del invernadero se observa la sección transversal del flujo de recirculación descrito en el apartado anterior, la sección transversal de un vórtice de escala similar al edificio 5x5 m. Como se ha comentado, la separación del flujo producido por el obstáculo en el techo a barlovento no se vuelve a unir aguas abajo del edificio, si no en el techo del invernadero, al ser este más alargado que alto. Por tanto al llegar el final del invernadero, se vuelve a producir una nueva separación del flujo, generando la estela aguas abajo, que produce un nuevo flujo de recirculación dominado por efectos cortantes. 104130 Resultados y discusión Figura 23 Contornos de v-velocidad y líneas de corriente en dos planos horizontales (3m y 4m) para una velocidad exterior de 4 m/s y modelo RSM. 105 131 Resultados y discusión Figura 24 Contorno de la velocidad-v en la sección media del invernadero (a) y líneas de corriente en la misma sección (b) 106 Mostrar más
Introducción a la Ing. Aeroespacial Tema 3 El Campo Fluido Sergio Esteban Roncero Francisco Gavilán Jiménez Departamento de Ingeniería Aeroespacial y Mecánica de Fluidos Escuela Técnica Superior de Ingeniería Más detalles Análisis del diseño de una chimenea solar
CAPÍTULO 4: ENSAYOS DE VALIDACIÓN MECÁNICA (CAE) Diseño, validación y fabricación de un aro protector para envases metálicos mediante el empleo de las tecnologías CAD/CAM/CAE y Rapid Prototyping. 4.1. Más detalles ASADES Avances en Energías Renovables y Medio Ambiente Vol. 8, Nº 2, 2004. Impreso en la Argentina. ISSN 0329-5184
ASADES Avances en Energías Renovables y Medio Ambiente Vol. 8, Nº 2, 2004. Impreso en la Argentina. ISSN 0329-5184 ANALISIS BIDIMENSIONAL DEL EFECTO DEL VIENTO SOBRE UN GRUPO DE VIVIENDAS N. Salvo, V. Más detalles ANSYS-Fluent como herramienta de diseño y evaluación de sistemas auxiliares en Invernaderos
ANSYS-Fluent como herramienta de diseño y evaluación de sistemas auxiliares en Invernaderos J. Flores-Velázquez, F. Villarreal, W. Ojeda y A. Rojano Jorge_flores@tlaloc.imta.mx Introducción Objetivos Contenido Más detalles INTRODUCCIÓN A LA CONVECCIÓN
Introdución a la capa límite laminar bidimensional y estacionaria M. Rodríguez 1 Introducción En los movimientos a altos números de Reynolds (basado en la longitud característica del movimiento), los efectos Más detalles SIMULACIÓN DE FLUIDOS UTILIZANDO COMPUTADORES UNA MODERNA HERRAMIENTA PARA EL ESTUDIO Y ANÁLISIS DE FLUIDOS
Técnicas computacionales aplicadas al tratamiento de aire. Análisis de velocidad y temperatura en la climatización de Edificios a través de CFD Juan Manuel Rodríguez González Esteban Nonay Villalba 1. Más detalles SIMULACIÓN DEL FLUJO TURBULENTO EN UN CODO DE 90 DE TUBERIA DE SECCIÓN CIRCULAR MANUEL ALEJANDRO TABORDA CEBALLOS
SIMULACIÓN DEL FLUJO TURBULENTO EN UN CODO DE 90 DE TUBERIA DE SECCIÓN CIRCULAR MANUEL ALEJANDRO TABORDA CEBALLOS UNIVERSIDAD AUTONOMA DE OCCIDENTE FACULTAD DE INGENIERIA DEPARTAMENTO DE ENERGETICA Y MECANICA Más detalles OPTIMIZACIÓN DE LOS SUBLIMADORES DE AZUFRE
OPTIMIZACIÓN DE LOS SUBLIMADORES DE AZUFRE USADOS EN LA INDUSTRIA FLORICULTORA MEDIANTE SIMULACIÓN NUMÉRICA EN ANSYS/FLOTRAN CFD José Ignacio Huertas*, Juan Manuel Caro* RESUMEN Los sublimadores de azufre Más detalles MODELADO NUMÉRICO DE LA REFRIGERACIÓN REGENERATIVA DE MOTORES COHETE
Tercer Congreso Argentino de Ingeniería Aeronáutica, CAIA 3 - La Plata, 12-14 de Noviembre de 2014 MODELADO NUMÉRICO DE LA REFRIGERACIÓN REGENERATIVA DE MOTORES COHETE E Fernández B. Becario CIC Grupo Más detalles Figura 10.1: Resultante de fuerzas horizontal y vertical sobre un perfil alar.
116 Capítulo 1 Capa Límite n cuerpo que este inmerso en un flujo eperimenta una fuerza resultante debido a la acción entre el flujo y el cuerpo. Esta es la fuerza resultante de los esfuerzos de corte en Más detalles MODELADO Y CONTROL MULTIVARIABLE DE TEMPERATURA Y HUMEDAD EN UN INVERNADERO
MODELADO Y CONTROL MULTIVARIABLE DE TEMPERATURA Y HUMEDAD EN UN INVERNADERO Autor: Cristina Sánchez Pérez Directores: José Luis Guzmán Sánchez y Jorge Antonio Sánchez Molina 19 de septiembre de 2014 2 Más detalles 4. Ecuaciones integrales de la transferencia de calor.
Departamento de Ingeniería Química 76.47 Fenómenos De Transporte -76.03 Operaciones I PROGRAMA ANALÍTICO 1. Nociones fundamentales de la mecánica de fluidos. Concepto de medio continuo. El fluido como Más detalles 2. ACTIVIDAD ACADÉMICA CÁLCULO EXPERIMENTAL DE PÉRDIDAS DE CARGA EN
Nombre de la asignatura: Mecánica de Fluidos I. Carrera : Ingeniería Mecánica Clave de la asignatura: MCB-9330 Clave local: Horas teoría horas practicas créditos: 4-0-8 2. - UBICACIÓN DE LA ASIGNATURA Más detalles La evapotranspiración: concepto y métodos para su determinación. Capítulo I
La evapotranspiración: concepto y métodos para su determinación Capítulo I I. La evapotranspiración: concepto y métodos para su determinación I.1 Evapotranspiración La evaporación es el proceso por el Más detalles VALIDACIÓN DE MODELOS DE VENTILACIÓN EN METROS Y TÚNELES FERROVIARIOS
VALIDACIÓN DE MODELOS DE VENTILACIÓN EN METROS Y TÚNELES FERROVIARIOS Ana Belén Amado García ZITRON Resumen: La precisión en los cálculos teóricos para obras subterráneas se vuelve más y más importante Más detalles TURBULENCIA 2. NATURALEZA Y CARACTERÍSTICAS DE LA TURBULENCIA 3. ORIGEN DE LA TURBULENCIA: INESTABILIDADES
Área de Mecánica de Fluidos. TURBUENCIA 1 Área de Mecánica de Fluidos TURBUENCIA 1. INTRODUCCIÓN. NATURAEZA Y CARACTERÍSTICAS DE A TURBUENCIA 3. ORIGEN DE A TURBUENCIA: INESTABIIDADES 4. A CASCADA DE ENERGÍA Más detalles ANÁLISIS DE LAS CONDICIONES DE FRONTERA EN LA INTERFASE AIRE AGUA (SUPERFICIE) PARA LA SIMULACIÓN HIDRODINÁMICA DE LAGUNAS FACULTATIVAS SECUNDARIAS
ANÁLISIS DE LAS CONDICIONES DE FRONTERA EN LA INTERFASE AIRE AGUA (SUPERFICIE) PARA LA SIMULACIÓN HIDRODINÁMICA DE LAGUNAS FACULTATIVAS SECUNDARIAS A. Aponte*, A. Toro*, L. Dueñas**, S. Laín**, M. R. Peña Más detalles Sistema termodinámico
IngTermica_01:Maquetación 1 16/02/2009 17:53 Página 1 Capítulo 1 Sistema termodinámico 1.1 Introducción En sentido amplio, la Termodinámica es la ciencia que estudia las transformaciones energéticas. Si Más detalles CAPÍTULO V MÉTODO DE ANÁLISIS POR ELEMENTOS FINITOS Y DESCRIPCIÓN DEL SOFTWARE ALGOR.
CAPÍTULO V MÉTODO DE ANÁLISIS POR ELEMENTOS FINITOS Y DESCRIPCIÓN DEL SOFTWARE ALGOR. 5.1 Historia del Método de Elementos Finitos. El Método de Elemento Finito o Finite Element Analysis (FEA), es un método Más detalles Mecánica de Fluidos y Máquinas Hidráulicas
PÉRDIDA DE CARGA Y EFICIENCIA ENERGÉTICA. Con unos costos de la energía en aumento y con unas limitaciones cada vez mayores a la emisión de gases de efecto invernadero, el diseño de equipos e instalaciones Más detalles TUBERIAS. Ricardo García San José Ingeniero Industrial (Noviembre 2.000) TUBERIAS
TUBERIAS Ricardo García San José Ingeniero Industrial (Noviembre 2.000) TUBERIAS INDICE 1.- MATERIALES... 3 2.- PERDIDAS DE CARGA... 4 2.1.- FACTORES QUE INFLUYEN EN LAS PERDIDAS DE CARGA... 4 2.2.- REGIMENES Más detalles Mecánica de Fluidos y Máquinas Hidráulicas
Mecánica de Fluidos y Máquinas Hidráulicas Tema 06. Flujo de Fluidos en Tuberías Severiano F. Pérez Remesal Carlos Renedo Estébanez DPTO. DE INGENIERÍA ELÉCTRICA Y ENERGÉTICA Este tema se publica bajo Más detalles TESIS. Que para obtener el Titulo de INGENIERO MECÁNICO PRESENTA FRANCISCO JAVIER VÁZQUEZ RODRÍGUEZ UNIVERSIDAD NACIONAL AUTÓNOMA DE MÉXICO
UNIVERSIDAD NACIONAL AUTÓNOMA DE MÉXICO FACULTAD DE INGENIERÍA ARRASTRE EN UN CILINDRO CUADRADO PRECEDIDO POR UNA PLACA DE SECCIÓN CUADRADA. TESIS Que para obtener el Titulo de INGENIERO MECÁNICO PRESENTA Más detalles FLUJOS TURBULENTOS BIFASICOS 3D EN LLENADO DE MOLDES DE FUNDICION Y EN IMPACTO DE OLAS EN ESTRUCTURAS Nelson Moraga, Carlos Garrido y Daniel Garrido
Mecánica de Energía Pérdidas de Energía Total Fluidos compresibles e incompresibles Los fluidos incompresibles son aquellos en los que el volumen permanece constante independientemente de las fuerzas aplicadas, Más detalles DIPLOMADO DINÁMICA DE FLUIDOS COMPUTACIONAL CFD
DIPLOMADO DINÁMICA DE FLUIDOS COMPUTACIONAL CFD 2015 INTRODUCCIÓN El Diplomado en Dinámica de Fluidos Computacional (CFD, por sus siglas en inglés Computational Fluid Dynamics), constituye uno de los pilares Más detalles O.C. Zienkiewicz, R.L. Taylor. El Método de los Elementos Finitos. Vols 1 y 2. CIMNE-Mc Graw Hill, 1994.
Índice de la teoría 1. Presentación. Estas lecciones sólo pretenden ser una introducción que sirva para orientar sobre los conceptos, para un estudio más amplio se recomienda leer alguna publicación especializada, Más detalles 1 AYUDA SOLIDWORKS MOTION
ORIENTACIONES PARA LA MATERIA DE FÍSICA Convocatoria 2010 Prueba de Acceso para Mayores de 25 años Para que un adulto mayor de 25 años pueda incorporarse plenamente en los estudios superiores de la Física Más detalles Criterios de cálculo y diseño de flujos laminares
Tecnología Industrial Criterios de cálculo y diseño de flujos laminares Patrones de Flujo del Aire. Condiciones de Flujo laminar José Luis Jiménez Álvarez Gerente / TCI, S.L.L. Los patrones de flujo de Más detalles Energías Renovables. Energía Eólica. Guía Técnica. Fotografía: Ministerio de Educación y Ciencia
Energías Renovables Energía Eólica Guía Técnica Fotografía: Ministerio de Educación y Ciencia INTRODUCCIÓN El objetivo de esta guía técnica es proporcionar una idea del potencial eólico de un emplazamiento Más detalles EL ANÁLISIS POR ELEMENTOS FINITOS: UNA METODOLOGÍA MUY RECIENTE EN ECONOMÍA
EL ANÁLISIS POR ELEMENTOS FINITOS: UNA METODOLOGÍA MUY RECIENTE EN ECONOMÍA Miguel Escribano Ródenas Departamento de Economía Financiera y Contabilidad I Escuela Universitaria de Estudios Empresariales Más detalles 27-02 Simulación CFD de un accidente de caída de un elemento combustible en piscina
27-02 Simulación CFD de un accidente de caída de un elemento combustible en piscina B. Montero 1, R. Corpa 2, G. Jiménez 1, C. Muñoz-Reja 2 1 Departamento de Ingeniería Energética (Área Nuclear), Universidad Más detalles Autor: Beatriz Sánchez Sánchez. Institución: Universidad Complutense de Madrid
Implementación de química reactiva en un modelo de mecánica de fluidos computacional aplicado a la calidad de aire urbano en el marco del proyecto Life Minox Street Autor: Beatriz Sánchez Sánchez Institución: Más detalles ESTUDIO DE MECANISMOS BÁSICOS DE TRANSFERENCIA DE CALOR CON CFD FLUENT
UNIVERSIDAD CARLOS III DE MADRID ESCUELA POLITÉCNICA SUPERIOR ÁREA: MECÁNICA DE FLUIDOS DEPARTAMENTO: INGENIERÍA TÉRMICA Y DE FLUIDOS PROYECTO FIN DE CARRERA INGENIERÍA TÉCNICA INDUSTRIAL, ESPECIALIDAD Más detalles VALIDACIÓN DE METODOLOGÍA CFD PARA SIMULACIÓN DE DISPOSITIVOS HIPERSUSTENTADORES SEGÚN EL 1ST AIAA CFD HIGH LIFT PREDICTION WORKSHOP
VALIDACIÓN DE METODOLOGÍA CFD PARA SIMULACIÓN DE DISPOSITIVOS HIPERSUSTENTADORES SEGÚN EL 1ST AIAA CFD HIGH LIFT PREDICTION WORKSHOP RESUMEN E. Gonzalez Garcia, C. Sacco y G. Weht Departamento de Ingeniería Más detalles CAPÍTULO 4 PRINCIPIOS DE MODELACIÓN DE EVENTOS EN LAGOR.
50 CAPÍTULO 4 PRINCIPIOS DE MODELACIÓN DE EVENTOS EN LAGOR. 4.1 Introducción a Algor simulación de eventos (MES). El futuro de la ingeniería asistida por computadora reposa en la habilidad para representar Más detalles Coordinador de la propuesta: Fernando Varas Mérida
Coordinador de la propuesta: Fernando Varas Mérida Representante de la Empresa: Andrés Gómez Tato TITULO DE LA ACTIVIDAD: JORNADAS DE CONSULTA MATEMÁTICA PARA EMPRESAS E INSTITUCIONES. I-MATH 2008-2011 Más detalles Figura 3.1.1 Vector AP en la trama {A}
3 Desarrollo 3.1 Vector de posición Un punto en el espacio puede ser ubicado utilizando un vector de posición, el cual tiene una dirección y magnitud. Estableciendo un sistema de coordenadas de dos ejes Más detalles Modelación de la combustión y la co-combustión en calderas de carbón pulverizado.
Modelación de la combustión y la co-combustión en calderas de carbón pulverizado. Pérez R. 1, Cornejo P. 2, Flores M. 3, Gordon A. 1, García X. 1 (1). Departamento de Ingeniería Química. Universidad de Más detalles Proyecto Fin de Carrera
Proyecto Fin de Carrera Análisis Computacional de la Dispersión de Contaminantes en Medios Urbanos Autor: Alejandro Sierra Luño Director: Antonio Pascau Benito Escuela de Ingeniería y Arquitectura 2014 Más detalles Sistema de distribución del aire. Cálculo de conductos.
Sistema de distribución del aire. Cálculo de conductos. Objetivos: Que el alumno sea capaz de dimensionar una red de conductos. Se pretende que el alumno pueda identificar los diferentes elementos que Más detalles CAPÍTULO 1 FUNDAMENTOS DEL ESTUDIO DE LA CAPA SUPERFICIAL ATMOSFÉRICA
CAPÍTULO 1 FUNDAMENTOS DEL ESTUDIO DE LA CAPA SUPERFICIAL ATMOSFÉRICA A continuación se exponen los conceptos básicos necesarios para desarrollar el presente trabajo de investigación en micrometeorología. Más detalles Caída de Presión en Tubos de Diferente Diámetro
Caída de Presión en Tubos de Diferente Diámetro Laboratorio de Operaciones Unitarias Equipo 4 Primavera 2008 México D.F., 12 de marzo de 2008 Alumnos: Arlette Mayela Canut Noval arlettecanut@hotmail.com Más detalles Mapas de viento: de la escala sinóptica al aerogenerador. Ignacio Martí imarti@cener.com Centro acional de Energías Renovables www.cener.
Mapas de viento: de la escala sinóptica al aerogenerador Ignacio Martí imarti@cener.com Centro acional de Energías Renovables www.cener.com (CE ER) Wind Power Expo 2009. Zaragoza. Áreas de conocimiento Más detalles Cinemática Inversa del Robot. CI-2657 Robótica M.Sc. Kryscia Ramírez Benavides
Cinemática Inversa del Robot M.Sc. Kryscia Ramírez Benavides Introducción Resuelve la configuración que debe adoptar el robot para una posición y orientación del extremo conocidas. 2 Introducción (cont.) Más detalles Introducción a OpenFoam Presentación general de OpenFoam
Introducción a OpenFoam Presentación general de OpenFoam E. Martín 1, M. Meis 2, F. Varas 1,3 1: Universidad de Vigo, 2: Vicus Desarrollos Tecnológicos, 3: Universidad Politécnica de Madrid Organización Más detalles Capa Límite Superficial
Capa ímite Superficial Física Ambiental. ema 6. ema6. FA (prof. RAMOS 1 ema 6.- Capa ímite Superficial. Capa límite: justificación. Flujos laminar y turbulento, características físicas: números de Reynolds. Más detalles Proyecto Fin de Carrera
Proyecto Fin de Carrera CÁLCULO COMPUTACIONAL DE PERFILES AERODINÁMICOS OS EN AEROTURBINA Víctor M. Martínez Fernández Director: Antonio Pascau Benito Departamento de Ciencia y Tecnología de Materiales Más detalles 9 CAPA LÍMITE. Re = [9.2.1] 9.1 Definición de la capa límite.
9 CAPA LÍMITE. 9.1 Definición de la capa límite. En 194, el Dr. Ludwig Prandtl desarrolló y presentó poco después en las célebres conferencias de la niversidad de Gottingen el concepto de capa límite, Más detalles Air Handling & Climate S. L. Difusión de aire
Air Handling & Climate S. L. Difusión de aire Concepto de difusión de aire El sistema de difusión es la parte terminal y visible de un sistema de climatización, y determina el éxito o el fracaso de la Más detalles Modelización Climática de Alta Resolución en terrenos de orografía compleja
Modelización Climática de Alta Resolución en terrenos de orografía compleja Índice Introducción Modelización del clima Técnicas de downscaling dinámico Participación de ULL en Climatique: Proyecciones Más detalles CONTENIDOS MÍNIMOS BACHILLERATO
SOLUCIONES PRESTACIONALES EN LA PCI DE EDIFICIOS SINGULARES ÍNDICE DE CONTENIDOS 1. Introducción 2. El fuego como fenómeno físico-químico 3. Simulación CFD: La principal herramienta 4. Técnicas de protección Más detalles POR QUÉ SOLIDWORKS FLOW SIMULATION ES LA OPCIÓN PERFECTA PARA LOS INGENIEROS DE PRODUCTO
I N F O R M E T É C N I C O POR QUÉ SOLIDWORKS FLOW SIMULATION ES LA OPCIÓN PERFECTA PARA LOS INGENIEROS DE PRODUCTO La dinámica de fluidos de ingeniería (EFD) es una nueva variedad de software de dinámica Más detalles TEMAS DE ESTUDIO FÍSICA
TEMAS DE ESTUDIO Estimados aspirantes a ingresar a las carreras de: INGENIERÍA AMBIENTAL, ELECTRICIDAD, INGENIERÍA INDUSTRIAL, INGENIERÍA DE SISTEMAS, INGENIERÍA ELECTRÓNICA, INGENIERÍA ELÉCTRICA, INGENIERÍA Más detalles ETAPA DE PREPROCESADO. Generación de geometría. Generación de malla ETAPA DE RESOLUCIÓN. Definición de modelos físicos
La simulación de flujos mediante técnicas computacionales se convertirá en un futuro cercano en una de las herramientas de diseño más valoradas por ingenieros y arquitectos dada su eficacia y versatilidad. Más detalles Capítulo 1 Introducción y planteamiento del proyecto.... 1 1.1 Motivación del proyecto.... 2 1.2 Objetivos del proyecto.... 3 1.3 Metodología de trabajo... 4 1.4 Recursos a utilizar... 4 Capítulo 2 Introducción Más detalles Dinamica de Fluidos: Principio de Bernoulli. Aplicaciones
Dinamica de Fluidos: Principio de Bernoulli. Aplicaciones Cuando un fluido está en movimiento, el flujo se puede clasificar en dos tipos: a) Flujo estacionario o laminar si cada partícula de fluido sigue Más detalles METODOLOGIA DE SUPERFICIES DE RESPUESTA. Esto se logra al determinar las condiciones óptimas de operación del sistema.
37 CAPITULO METODOLOGIA DE SUPERFICIES DE RESPUESTA En este capítulo hablaremos de qué es la Metodología de Superficies de Respuesta, su representación gráfica, el procedimiento a seguir hasta encontrar Más detalles http://saeti.itson.mx/otrosusuarios/plandosmilnueveconsprogamplioimpma.asp?materia...
Page 1 of 5 Departamento: Dpto Cs. Agua y Medio Ambiente Nombre del curso: TRANSFERENCIA DE CALOR CON LABORATORIO Clave: 004269 Academia a la que pertenece: Academia de Ingeniería Química Aplicada en Operaciones Más detalles Carrera: EMM - 0525. Participantes Representante de las academias de ingeniería Electromecánica de los Institutos Tecnológicos.
1.- DATOS DE LA ASIGNATURA Nombre de la asignatura: Carrera: Clave de la asignatura: Horas teoría-horas práctica-créditos Mecánica de Fluidos Ingeniería Electromecánica EMM - 0525 3 2 8 2.- HISTORIA DEL Más detalles Tema 4 Difusión en estado sólido
Tema 4 Difusión en estado sólido Sabemos que los materiales están formados por átomos. Se ha modelado el agrupamiento de los átomos como un conjunto de esferas sólidas ordenadas siguiendo un patrón definido. Más detalles El uso de técnicas de simulación para la predicción de defectos en piezas fundidas.
Ingeniería Mecánica 3 (2006) 47-51 47 El uso de técnicas de simulación para la predicción de defectos en piezas fundidas. T. Rodríguez Moliner, A. Parada Expósito, U. Ordóñez. *Departamento de Tecnología Más detalles Dinámica de Fluidos Computacional para Gráficos
Técnicas Avanzadas de Gráficos 3D Dinámica de Fluidos Computacional para Gráficos Álvaro Pérez Molero Carlos Garre del Olmo Indice 1 Introducción... 3 1.1 Qué es la dinámica de fluidos?... 3 1.2 Qué es Más detalles MANEJO DE HUMO Y CALOR SISTEMAS DE PRESURIZACION PARA CAJA DE ESCALERAS
MANEJO DE HUMO Y CALOR SISTEMAS DE PRESURIZACION PARA CAJA DE ESCALERAS OBJETO Y ALCANCE Es objeto del presente Informe mencionar los fundamentos y conceptos básicos para el Manejo y/o Control de Humo, Más detalles NÚCLEO 4 SISTEMA DE CONDUCCIÓN HIDRÁULICA 4.1 CARÁCTERÍSTICAS HIDRÁULICAS DEL SISTEMA
NÚCLEO 4 SISTEMAS DE CONDUCCIÓN HIDRAÚLICA 4.1 CARÁCTERÍSTICAS HIDRÁULICAS DEL SISTEMA La conducción en un sistema de bombeo es uno de los elementos más importantes, ya que su función es precisamente formar Más detalles 5. PÉRDIDAS DE CARGA EN CONDUCTOS CERRADOS O TUBERIAS
5. PÉRIAS E CARGA EN CONUCTOS CERRAOS O TUBERIAS 5. Perfiles de Velocidad: Laminar y Turbulento 5. Radio Hidráulico para Secciones no Circulares 5.3 Pérdidas Primarias y Secundarias 5.4 Ecuación de arcy Más detalles ANEXO 5 REGLAS BÁSICAS PARA LA IMPLEMENTACIÓN DE MODELOS CFD EN EL MODELO DEL VERTIDO EN CAMPO CERCANO V. 1.1
ANEXO 5 REGLAS BÁSICAS PARA LA IMPLEMENTACIÓN DE MODELOS CFD EN EL MODELO DEL VERTIDO EN CAMPO CERCANO V. 1.1 Este documento se ha obtenido como resultado del proyecto de I+D+i: MEDVSA. Desarrollo e implementación Más detalles CURSO TALLER PROMOTORES DE AHORRO Y EFICIENCIA DE ENERGÍA ELÉCTRICA
PROGRAMA INTEGRAL DE ASISTENCIA TÉCNICA Y CAPACITACIÓN PARA LA FORMACIÓN DE ESPECIALISTAS EN AHORRO Y USO EFICIENTE DE ENERGÍA ELÉCTRICA DE GUATEMALA CURSO TALLER PROMOTORES DE AHORRO Y EFICIENCIA DE ENERGÍA Más detalles Energía eólica. AEROGENERADORES: Funcionamiento, partes y tipos. 1. Funcionamiento. 2. Partes
Energía eólica La energía eólica tiene su origen en el viento, es decir, en el aire en movimiento. El viento se puede definir como una corriente de aire resultante de las diferencias de presión en la atmósfera Más detalles UNIVERSIDAD DE CHILE FACULTAD DE CIENCIAS FÍSICAS Y MATEMÁTICAS DEPARTAMENTO DE INGENIERÍA CIVIL
UNIVERSIDAD DE CHILE FACULTAD DE CIENCIAS FÍSICAS Y MATEMÁTICAS DEPARTAMENTO DE INGENIERÍA CIVIL MODELACIÓN NUMÉRICA 3D APLICADA AL DISEÑO DE LAS OBRAS HIDRÁULICAS DE LA CENTRAL BAKER 1, REGIÓN DE AYSÉN Más detalles INDICE: Concepto y fundamentos de la Ingeniería Química Fenómenos de transporte Balances de materia
INDICE: Concepto y fundamentos de la Ingeniería Química Fenómenos de transporte Balances de materia INTRODUCCIÓN A LA INGENIERÍA DE LOS PROCESOS QUÍMICOS. Tema desarrollado por los profesores Dr. D. Antonio Más detalles Ingeniería de estructuras en situación de Incendio
Ingeniería de estructuras en situación de Incendio Garantizar la resistencia al fuego de la estructura de un edificio en caso de incendio para impedir que este se propague de unos recintos a otros así Más detalles SIMULACIÓN NUMÉRICA CFD APLICADA AL ESTUDIO Y OPTIMIZACIÓN DEL COMPORTAMIENTO TÉRMICO DE PANELES FOTOVOLTAICOS INTEGRADOS EN CUBIERTA
SIMULACIÓN NUMÉRICA CFD APLICADA AL ESTUDIO Y OPTIMIZACIÓN DEL COMPORTAMIENTO TÉRMICO DE PANELES FOTOVOLTAICOS INTEGRADOS EN CUBIERTA Dr. Ing. Marcos Oswaldo Quispe Flores marcos.quispe@teksode.com Dr. Más detalles TÉCNICAS DE CARACTERIZACIÓN DE CAMPO DE VELOCIDADES
TÉCNICAS DE CARACTERIZACIÓN DE CAMPO DE VELOCIDADES La velocimetría de partículas es una de las técnicas más recientes que permite determinar el campo de velocidad del flujo. Los métodos conocidos basados Más detalles EJEMPLO DE CÁLCULO DE UN EQUIPO DE BOMBA DE CALOR POR EL MÉTODO DE LOS SALTOS DE TEMPERATURA. Poveda, M.; Poveda, R; León, M; Poveda, E; Tarifa, M.
X CONGRESO INTERNACIONAL DE INGENIERIA DE PROYECTOS VALENCIA, 13-15 Septiembre, 2006 EJEMPLO DE CÁLCULO DE UN EQUIPO DE BOMBA DE CALOR POR EL MÉTODO DE LOS SALTOS DE TEMPERATURA Poveda, M.; Poveda, R; Más detalles SISTEMAS DE CONTROL I MODELADO DE SISTEMAS FÍSICOS
SISTEMAS DE CONTROL I MODELADO DE SISTEMAS FÍSICOS Ing. Miguel G. Alarcón Agosto de 2011 Temario Sistema Físico. Modelado del Sistema Real. Sistemas Eléctricos. Sistemas Mecánicos. Sistemas Térmicos. Qué Más detalles TIMSS 11.2 DESCRIPCIÓN DE LO EVALUADO EN LOS DOMINIOS DE CONTENIDO MATEMÁTICA Números Incluye la comprensión del proceso de contar, de las maneras de representar los números, de las relaciones entre éstos Más detalles M.C. José Eduardo Frias Chimal Dr. Sergio Alonso Romero M.C. Miguel Ángel Corzo Velázquez M.I. Raúl Lesso Arroyo
Rediseño Geométrico y Validación Numérica de la Optimización de la Transferencia de Calor en Campañas de Cocción de Ladrillo Rojo en el Estado de Guanajuato M.C. José Eduardo Frias Chimal Dr. Sergio Alonso Más detalles TERMODINAMICA 1 Conceptos Basicos
TERMODINAMICA 1 Conceptos Basicos Prof. Carlos G. Villamar Linares Ingeniero Mecánico MSc. Matemáticas Aplicada a la Ingeniería 1 CONTENIDO DEFINICIONES BASICAS Definición de Termodinámica, sistema termodinámico, Más detalles Simulación numérica de la distribución de temperaturas y flujo de aire en el horno de una cocina de uso doméstico
Simulación numérica de la distribución de temperaturas y flujo de aire en el horno de una cocina de uso doméstico UNIVERSIDAD POLITÉCNICA SALESIANA SEDE CUENCA UNIDAD DE POSGRADOS MAESTRÍA EN MÉTODOS Más detalles CAPITULO 6 ANALISIS Y ESTUDIO DE SECADO. El secado de sólidos se puede definir de distintas maneras, según el enfoque que se
52 CAPITULO 6 ANALISIS Y ESTUDIO DE SECADO 6.1 Definición de secado El secado de sólidos se puede definir de distintas maneras, según el enfoque que se desee adoptar. En los estudios más teóricos se pone Más detalles ANÁLISIS DE LOS ESTABILIZADORES DEL TABLERO DE UN PUENTE ATIRANTADO DE GRAN IMPORTANCIA ANTE EL DESPRENDIMIENTO DE VÓRTICES CON CFD RESUMEN ABSTRACT
Sociedad Mexicana de Ingeniería Estructural ANÁLISIS DE LOS ESTABILIZADORES DEL TABLERO DE UN PUENTE ATIRANTADO DE GRAN IMPORTANCIA ANTE EL DESPRENDIMIENTO DE VÓRTICES CON CFD Fabio Alberto Miguel Patiño Más detalles Estudio y simulación de la influencia de la estructura Transformador-Bobina Paralelo en convertidores CC-CC clásicos
ESCUELA POLITÉCNICA SUPERIOR Grupo de Sistemas Electrónicos de Potencia PROYECTO FIN DE CARRERA INGENIERÍA INDUSTRIAL Estudio y simulación de la influencia de la estructura Transformador-Bobina Paralelo Más detalles I.- CONCEPTOS BÁSICOS DEL MÉTODO DE LOS ELEMENTOS FINITOS
I.- CONCEPTOS BÁSICOS DEL MÉTODO DE LOS ELEMENTOS FINITOS 1.1.- Introducción Muchos problemas de importancia práctica que frecuentemente aparecen en ingeniería, resultan de una complejidad matemática tal Más detalles USO DE UN ESCÁNER 3D PARA APLICACIONES EN INGENIERÍA DE RIEGOS
USO DE UN ESCÁNER 3D PARA APLICACIONES EN INGENIERÍA DE RIEGOS Playán, E. 1, Serreta, A. 2, Zapata, N. 1, Burguete, J., Salvador, R., Martínez-Cob, A., Cavero, J. RESUMEN La ingeniería del riego precisa Más detalles Elaboración de un modelo basado en CFD para predecir el comportamiento de un Aerogenerador de Eje Vertical
Máster en Energía: Generación, Gestión y Uso Eficiente. TRABAJO FIN DE MÁSTER Elaboración de un modelo basado en CFD para predecir el comportamiento de un Aerogenerador de Eje Vertical Autora: Carmen Victoria Más detalles MECÁNICA DEL SOLIDO RÍGIDO
MECÁNICA DEL SOLIDO RÍGIDO ESTÁTICA 4.- Estática, Equilibrio La estática es la parte que estudia las interacciones - fuerzas, cargas, momentos- entre cuerpos o entre las diferentes partes de un cuerpo, Más detalles Física P.A.U. ELECTROMAGNETISMO 1 ELECTROMAGNETISMO. F = m a
Física P.A.U. ELECTOMAGNETISMO 1 ELECTOMAGNETISMO INTODUCCIÓN MÉTODO 1. En general: Se dibujan las fuerzas que actúan sobre el sistema. Se calcula la resultante por el principio de superposición. Se aplica Más detalles PROGRAMA DE CURSO. Horas de Trabajo Personal 6 10 3,0 2,0 5,0. Horas de Cátedra
PROGRAMA DE CURSO Código FI1002 Nombre Sistemas Newtonianos Nombre en Inglés Newtonian Systems SCT Unidades Docentes Horas de Cátedra Horas Docenci a Auxiliar Horas de Trabajo Personal 6 10 3,0 2,0 5,0 Más detalles 2017 © DocPlayer.es Política de privacidad | Condiciones del servicio | Feedback

References: RESOLUCIÓN 
 RESOLUCIÓN 
 RESOLUCIÓN 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 Resolución 
 resolución 
 resolución 
 Resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 Resolución 
 Resolución 
 resolución 
 resolución 
 resolución 
 Resolución 
 Resolución