Un esquema adaptativo semilagrangiano con elementos finitos anisótropos para la resolución de problemas de combustión.

En este trabajo se presenta un método numérico adaptativo para la resolución de problemas temporales de Combustión en la Mecánica de Fluidos. La dis- cretización espacial de las ecuaciones está basada en el método de los elementos finitos y la discretización temporal se realiza con un esquema semilagrangiano para tratar de forma eficiente los términos convectivos. La caracteística principal de la adaptación local, es que el mallado que se construye es anisótropo, lo que le capacita para adaptarse mejor a capas límitee con un reducido coste computacional. Para ello, se define un tensor métrico en cada paso de tiempo, basado en indicadores de error construidos a priori y a posteriori. Ilustraremos el buen comportamiento del código numérico con la modelización de un problema de combustión en 2D y 3D, donde se analilzará la interacción de llamas de difusión de Hidrógeno y vórtices que pueden ser generados en un flujo turbulento.

Modelización de la rotura y fluidificación en geomateriales y geoestructuras

La rotura de las geoestructuras puede ser causada por cambios en las tensiones efectivas debidos bien por cargas externas (terremotos, por ejemplo), bien por variación de las presiones intersticiales (lluvia), o cambios en la geometría (erosión), así como por una disminución de las propiedades resistentes de los materiales (meteorización, ataque químico, etc). El caso particular de los deslizamientos es interesante, existiendo diversas clasificaciones que tienen en cuenta su forma, velocidad de propagación, etc. Dos de estos casos son los deslizamientos propiamente dichos y los flujos. En el primer caso, la deformación se concentra en zonas de pequeño espesor, que se idealiza como una superficie (superficie de rotura). La cinemática de esta rotura se puede considerar como el movimiento relativo de dos masas cuyas deformaciones no son grandes. Este mecanismo está usualmente asociado a materiales sobreconsolidados que presentan reblandecimiento. Los flujos se producen generalmente en materiales de baja densidad y estructura metaestable, con tendencia a compactar, de forma que se generan presiones intersticiales que aumentan el ángulo de rozamiento movilizado, pudiéndose llegar en algunos casos a la licuefacción. Este mecanismo de rotura se conoce como rotura difusa, y no ha sido tan estudiado como el de localización a pesar de se trata frecuentemente de roturas de tipo catastrófico. De hecho, el suelo pasa de un estado sólido a un estado fluidificado, con una gran movilidad. Es importante para el ingeniero predecir tanto el comportamiento de las geoestructuras bajo las cargas de cálculo como las condiciones en las que se producirá la rotura. De esta manera, en algunos casos, se podrán reforzar las zonas más débiles aumentando así su seguridad. En otros casos, no se podrá realizar este refuerzo (grandes deslizamientos como avalanchas, lahares, etc), pero sí se podrán conocer las consecuencias de la rotura: velocidad de propagación, alcance, espesores, etc. La modelización de estos problemas es compleja, ya que aparecen dificultades en los modelos matemáticos, constitutivos o reológicos y numéricos. Dado que en los geomateriales aparece una interacción fuerte entre el esqueleto sólido y los fluidos intersticiales, esto debe ser tenido en cuenta en los modelos matemáticos. En este trabajo se describirán, pues, el desarrollo y aplicación de técnicas avanzadas de modelización; matemática, constitutiva/reológica y numérica. Se dedicará especial atención a los problemas de transición entre suelos y suelos fluidificados, que hoy en día se estudian en una gran mayoría de los casos con modelos diferentes. Así por ejemplo, se emplean modelos constitutivos para el comportamiento previo a la rotura, y reológicos para los materiales fluidificados. En lo que respecta a los modelos matemáticos, existen formulaciones nuevas en velocidades (o desplazamientos), tensiones, y presiones de aire y agua intersticial, de los que se pueden obtener modelos simplificados integrados en profundidad para deslizamientos rápidos. Respecto de los modelos constitutivos, es interesante la teoría de la Plasticidad Generalizada (modelo básico de Pastor-Zienkiewicz y extensiones a suelos no saturados). Se estudiará la extensión de estos modelos elastoplásticos a la viscoplasticidad (Perzyna y consistente), explorando la posibilidad de emplearlos tanto antes como después de la rotura. Finalmente, en lo que a modelos numéricos se refiere, se describirá la implementación de los modelos matemáticos y constitutivos desarrollados en (i) modelos clásicos de elementos finitos, como el GeHoMadrid desarrollado en los grupo investigador M2i al que pertenece el autor de este trabajo, (ii) Métodos de tipo Taylor Galerkin, y (iii) métodos sin malla como el SPH y el Material Point Model. Estos modelos se aplicarán, principalmente a (i) Licuefacción de estructuras cimentadas en el fondo marino (ii) presas de residuos mineros (iii) deslizamientos rápidos de laderas.

Particularidades geológicas del yacimiento aurífero la holguinera, parte noroccidental de la región de Holguín, Cuba.

El yacimiento Holguinera forma parte del campo metalogénico de Aguas Claras, ubicado aproximadamente a 10,5 kilómetros al norte de la ciudad de Holguín, principal núcleo poblacional de la provincia a la que da nombre. Las coordenadas del punto medio del yacimiento son 20º56’18,81’’N y 76º16’59,18’’O; mientras la cota es de 151 metros. Desde el punto de vista geográfico-administrativo, se encuentra al suroeste de la localidad Aguas Claras y colinda al oeste con el caserío de Guajabales. Este yacimiento se ubica en el extremo oriental del conocido yacimiento Reina Victoria, con el cual se conecta mediante pequeñas manifestaciones e indicios auríferos entre los que se mencionan Santiago, Casualidad, Milagro, Relámpago y Nené El yacimiento La Holguinera es uno de los más antiguos de los que se conocen en esta región. Según los informes que aún se conservan en la Oficina Nacional de Recursos Minerales de Cuba. El yacimiento La Holguinera se integra regionalmente en la mega formación conocida como Zona Estructuro-Facial Auras y representa en sí misma una alternancia de formaciones litológicas plegadas y deformadas que proceden del Complejo Ofiolítico y del Arco Volcánico Cretácico. La cobertera que yace sobre ellas consiste en materiales volcánicos y volcano-sedimentarios, que son los productos de la erosión y sedimentación de las principales unidades litológicas. La aflorabilidad del yacimiento La Holguinera y sus flancos es buena, pero se debe más bien a la actividad minera tradicional, gracias a la cual ha quedado al descubierto en superficie gran parte de las menas oxidadas. En la actualidad se observa una gran profusión de trincheras, cortas, sondeos, pozos mineros, desagües, portalones y cortes, que son los testigos de la frenética actividad minera de esos años. Los estudios mineragráficos y petrográficos de las texturas volcánicas relícticas sugieren que la mineralización aurífera del yacimiento La Holguinera se encuentra encajada en una roca de composición andesítica de textura fina microporfirítica, y de forma secundaria en la roca encajante del tipo harzburguita, alterada a listvenitas. Los mismos fluidos que alteraron a los cuerpos intrusivos difundieron hacia el encajante ultramáfico, originando una variedad litológica secundaria conocida como listvenita, compuesta por productos de alteración metasomática del tipo talco-carbonato-cuarzo-sulfuro, con mineralización aurífera acompañante (listvenitas). Las menas de este yacimiento son sulfurosas, sulfuro-polimetálicas, pirítica, cromítica y sulfurosas de hierro.

Uso del cálculo numérico CFD para el estudio de actuaciones de una motocicleta de competición = CFD analysis of a competition motorbike aerodynamic performance

Esta memoria pretende ilustrar el desarrollo del Proyecto de Fin de Carrera Uso del cálculo numérico CFD para el estudio de actuaciones de una motocicleta de competición de los alumnos Álvaro Ramos Cenzano y Mauro Botella Mompó. Se centra en el estudio de las fuerzas aerodinámicas ejercidas sobre una motocicleta de competición, en concreto de una Moto2 del campeonato de España de velocidad (CEV), por medio de programas de simulación. Los capítulos del 1 al 4, describen por qué surgió el Proyecto, los datos a tener en cuenta del reglamento oficial de Moto2 del CEV (centrándonos sobre todo en los aspectos que afectan a la Aerodinámica, como puede ser el carenado), y los objetivos que se pretenden hallar con este proyecto. Se hará una pequeña introducción a la aerodinámica general para después profundizar en la aerodinámica automovilística y sobretodo en la de las motocicletas. Pasando por el interés que tiene la aerodinámica a diferentes niveles (comercial y de competición), y su evolución a lo largo de la historia. Se analizarán los diferentes métodos de trabajo que se podrían usar para el estudio aerodinámico, centrándose sobre todo en la técnica de estudio informático por Dinámica de Fluidos Computacional o CFD. Dado que el proyecto está basado en unas motocicletas reales, surgirán una serie de problemas en el uso de los métodos numéricos que se irán solventando, intentando mantener la estructura lo más fiel posible a los modelos reales. Todos estos problemas y sus soluciones se detallarán con la intención de que puedan servir de ayuda para futuros estudios. Tras las simulaciones con los programas CFD correspondientes a los distintos modelos de motocicleta, se compararán los resultados para ayudar en la definición del modelo final. El proyecto concluirá con las conclusiones obtenidas por los alumnos.

Evaluación del comportamiento mecánico de estructuras de hormigón armado utilizadas como elementos termoactivos

El objetivo de este trabajo es conocer las posibles modificaciones que puede producir en el comportamiento de las estructuras de hormigón armado (EHA) el hecho de que sean utilizadas como estructuras termoactivas, ya sea como intercambiadores en contacto con el terreno, o como sistema de distribución de calor utilizando la inercia térmica de los elementos de hormigón del edificio, basándose en el uso de energías renovables. Las EHA termoactivas se caracterizan por la incorporación en su interior de tubos de polietileno por los que circulan fluidos a temperaturas medias, que pueden incidir en el comportamiento mecánico de los elementos estructurales debido a dos efectos fundamentales: el incremento de temperatura que se produce en el interior de la EHA y la perturbación provocada por la incorporación de los tubos de polietileno. Con este fin, se ha realizado una campaña experimental de probetas de hormigón, estudiando los dos efectos por separado, por un lado se ha evaluado el comportamiento de probetas de hormigón tipo H-25 y tipo H-30 sometidas a cuatro temperaturas diferentes: 20ºC, 40ºC, 70ºC y 100ºC, ensayando la resistencia a compresión y la resistencia a anclaje/adherencia mediante ensayo “pull-out”; y, por otro lado, se ha evaluado el comportamiento de probetas de hormigón tipo H-25 y tipo H-30, elaboradas con dos tipos de molde (cilíndrico y cúbico), en las que se ha colocado tubos de polietileno en su interior en distintas posiciones, ensayando su resistencia a compresión. Los resultados de los ensayos han puesto de manifiesto que aunque se produce una disminución en la resistencia a compresión, y a arrancamiento, del hormigón, al ser sometido a aumentos de temperatura, esta disminución de la propiedades mecánicas es inferior al 20% al no superar esta tecnología los 70ºC; y respecto a la variación de la resistencia a compresión de probetas cilíndricas y cúbicas, debidas a la incorporación de los tubos de polietileno, se observa que si la posición de los mismos es paralela a la dirección de la carga tampoco se ven comprometidas las propiedad mecánicas del hormigón en valores superiores al 20%. ABSTRACT The aim of this project is to study the effects of using concrete structures as thermo-active constructions, either as energy foundations or other kind of thermo-active ground structures, or as a thermally activated building structure utilizing its own thermal mass conductivity and storage capacity to heat and cool buildings, based on renewable or “free” energy sources. The pipes, filled with a heat carrier fluid, that are embedded into the building´s concrete elements may bring on two different adverse effects on concrete structures. In one hand, the consequence of thermal variations and, on the other hand, because of the fact that the pipes are inside of the concrete mortar and in direct contact with the reinforcing steel bars. For this reason, different types of specimens and testing procedures have been proposed to discuss the effects of temperature (20º, 40ºC, 70ºC y 100ºC) on the performance of two different hardened concrete: H-25 and H-30, and the effects of having the pipes embedded in different positions inside of specimens made of two types of concrete, H-25 and H-30, and with two kind of cast, cylindrical and cubical. The experimental program includes the use of compressive strength test and also pull-out test, in order to investigate the interfacial adhesion quality and interfacial properties between steel bar and concrete. The results of the mechanical test showed that the increase of temperature in hardened concrete specimens lower than 70ºC, and the introduction of embedded pipes placed in parallel to the load, in cylindrical or cubic specimens, does not jeopardize the mechanical properties of concrete with strength decreases higher than 20%.

Pirolisis y gasificación de biomasas: optimización del proceso de transformación energética mediante distintos sistemas de reacción

Las proyecciones muestran un incremento en la demanda de energía primaria a nivel global para los próximos diez años. Para disminuir los efectos adversos sobre el medio ambiente del uso de combustibles fósiles, las políticas energéticas de España y la Unión Europea tienen como uno de sus objetivos principales incrementar el uso fuentes renovables en la producción energética para el año 2020.La biomasa es una de las principales fuentes renovables para la generación de energía primaria en España. Sin embargo su alto contenido en humedad y su baja densidad energética entre otras características de este material, suponen desventajas para su uso directo como combustible. Los procesos de conversión termoquímica son una alternativa para la obtención de productos sólidos, líquidos y gases de mejor potencial energético utilizando biomasa como materia prima. Esta tesis doctoral tiene como objetivo general el estudio y la optimización de procesos de conversión termoquímica de biomasa. Para cumplir con este objetivo el trabajo combina análisis experimentales y de simulación de distintas tecnologías de transformación termoquímica de biomasa. Antes de iniciar los estudios experimentales de los procesos de transformación termoquímica, se realiza una caracterización química, física y combustible de ocho de biomasas de naturaleza lignocelulósica y amilácea. Esta caracterización incluye el análisis del comportamiento termoquímico utilizando las técnicas de termogravimetría (TG-DTG), TG-FTir y pirólisis analítica (Py-GC/MS). Después de la caracterización se seleccionan tres tipos de biomasas (madera de pino, hueso de aceituna y hueso de aguacate) que son utilizadas como materia prima en el estudio experimental de procesos de gasificación, torrefacción y pirólisis en distintos sistemas de reacción de escala de laboratorio. Estos estudios se realizan para optimizar las condiciones de operación de cada uno de los procesos y caracterizar los productos obtenidos. El trabajo realizado incluye el rediseño y puesta en marcha de un sistema de gasificación de lecho fluidizado así como de un sistema de torrefacción y pirólisis en horno rotatorio. Los estudios experimentales incluyen un análisis de distintas condiciones de operación (temperatura, composición atmosférica, tiempo de residencia, condiciones fluidomecánicas, tamaño de partícula del combustible, condiciones de condensación) con el objetivo de determinar las condiciones óptimas de operación de cada uno de los procesos analizados respecto de las características químicas, físicas y combustibles de los productos. Para la optimización de las condiciones de operación del sistema de gasificación de lecho fluidizado se realiza un análisis complementario del comportamiento hidrodinámico utilizando la metodología de análisis computacional fluidodinámico o mecánica de fluidos computacional (CFD). Los resultados obtenidos permiten determinar las condiciones óptimas de operación del sistema de reacción de lecho fluidizado. Finalmente, la tesis incluye el análisis de una planta comercial de gasificación localizada en Júndiz (Álava, España), centrándose en la caracterización del material particulado eliminado de la corriente de gas generado, con el objetivo de determinar su proceso de formación, características y establecer posibles aplicaciones comerciales.

Desarrollo de un sistema de escaneado 3D a partir de un dispositivo de holografía conoscópica

La radioterapia intraoperatoria (RIO) con electrones es una modalidad de tratamiento contra el cáncer que combina cirugía y radiación terapéutica aplicada a un tumor sin resecar o al lecho tumoral después de una resección. Este tratamiento se planifica con una imagen preoperatoria del paciente. Sería adecuado incluir en la planificación información acerca de la superficie del lecho tumoral ya que la acumulación de fluidos, las superficies cóncavas y las irregularidades en la superficie irradiada modifican significativamente la distribución de la dosis. Existen diversos escáneres de superficie pero el que mejor se podría adaptar al entorno de la RIO sería el que utiliza un dispositivo de holografía conoscópica (ConoProbe, Optimet Metrology Ltd.) ya que permite realizar medidas de distancias en cavidades, en superficies reflectantes y en tejidos biológicos. La holografía conoscópica es una técnica de interferometría basada en la propagación de la luz en cristales birrefringentes. Para obtener las coordenadas 3D de los puntos de la superficie barridos por el dispositivo ConoProbe es necesario utilizar un sistema de posicionamiento como el que se utiliza en el quirófano de la RIO para localizar el aplicador de RIO con el que se conduce el haz de electrones (sistema de posicionamiento óptico OptiTrack, NaturalPoint Inc.). El objetivo de este proyecto fin de grado consistió en desarrollar un sistema de escaneado de superficies y realizar diversas pruebas para evaluar la calidad del sistema desarrollado y su viabilidad en el entorno de la RIO. Para integrar la información del dispositivo ConoProbe y del sistema de posicionamiento OptiTrack se realizó una calibración temporal para sincronizar los datos de ambos dispositivos utilizando la función de correlación cruzada y una calibración espacial para transformar la distancia medida por el dispositivo ConoProbe en coordenadas 3D de la superficie del objeto escaneado. Se plantearon dos métodos para realizar esta calibración espacial, por ajuste de pares de puntos y por ajuste a un plano. La calibración espacial elegida fue la primera ya que presentaba menor error. El error del sistema es inferior a 2 mm siendo un error aceptable en el entorno de la RIO. Diferentes pruebas con diversos materiales y formas han permitido comprobar que el sistema de escaneado funciona incluso con líquidos. En un procedimiento de RIO, el escaneado de la superficie se haría después de colocar el aplicador. En este caso, se puede obtener también la superficie del objeto alrededor de centro del aplicador colocando verticalmente el dispositivo ConoProbe. Es complicado obtener la superficie próxima a las paredes del aplicador debido a que estas afectan al cono de luz reflejado. De todas formas, el sistema de escaneado desarrollado proporciona más información en este escenario que un sistema de escáner 3D de luz estructurada (no se podría escanear nada de la superficie con el aplicador colocado). Esta información es útil para la estimación de la distribución de la dosis real que recibe un paciente en un procedimiento de RIO.

Preparados enológicos comerciales a base de levaduras secas inactivas : caracterización, modo de acción e influencia en la composición y características sensoriales de los vino

El empleo de preparados comerciales de levaduras secas inactivas (LSI) se ha generalizado en la industria enológica con el fin de mejorar los procesos tecnológicos, así como las características sensoriales y la eliminación de compuestos indeseables en los vinos. Estos preparados se obtienen por crecimiento de la levadura vínica Saccharomyces cerevisiae en un medio rico en azúcares, y posterior inactivación y secado para obtener un producto en forma de polvo. Su composición puede incluir una fracción soluble procedente del citoplasma de la levadura (péptidos, aminoácidos, proteínas) y una fracción insoluble compuesta principalmente por las paredes celulares. Dependiendo de su composición, los preparados de LSI están indicados para distintas aplicaciones (activadores de la fermentación alcohólica y maloláctica, mejorar el color y el aroma de los vinos, etc). Sin embargo, a pesar de que el empleo de preparados de LSI está cada vez más extendido en la industria enológica, existen pocos estudios científicos encaminados a conocer sus efectos en el vino así como a discernir su modo de acción. Por ello, se ha planteado esta Tesis Doctoral con el principal objetivo de caracterizar la composición química de algunos de los preparados comerciales de LSI que actualmente más se están empleando en la elaboración de los vinos, así como determinar su modo de acción y su efecto tanto en la composición como en las características organolépticas de los vinos. Para la caracterización de preparados de LSI, en primer lugar se ha estudiado la cesión de compuestos nitrogenados y polisacáridos a medios vínicos, así como las diferencias entre preparados procedentes de diferentes casas comerciales e indicados para diferentes aplicaciones durante la vinificación. Los resultados de este estudio mostraron diferencias en la cesión de estos compuestos dependiendo del tipo de vino al que van dirigidos y de la casa comercial de la procedencia. Así, los preparados indicados para la elaboración de vinos tintos liberan una mayor concentración de polisacáridos, que están asociados a una modulación de la astringencia y de la protección del color de los vinos tintos, comparado con los indicados para vinos blancos. Sin embargo en estos últimos se determinó una mayor cesión de aminoácidos libres. Posteriormente, se estudió el efecto de componentes específicos presentes en preparados de LSI en el crecimiento de distintas especies de bacterias lácticas del vino. Se obtuvieron extractos de diferente composición empleando la extracción con fluidos presurizados, y se ensayó su actividad en el crecimiento de cepas bacterianas pertenecientes a las especies Lactobacillus hilgardii, Pediococcus pentosaceus y Oenococcus oeni. Los resultados de este trabajo mostraron que los distintos componentes presentes en los preparados de LSI pueden ejercer un efecto estimulante o inhibidor en el crecimiento de las bacterias lácticas, que depende del tipo de compuesto, de su concentración y de la especie bacteriana a ensayar. Entre los compuestos responsables del efecto estimulante se encuentran tanto aminoácidos esenciales como monosacáridos libres, mientras que los ácidos grasos de cadena corta y larga así como compuestos heterocíclicos volátiles originados por la reacción de Maillard parecen estar directamente implicados en el efecto inhibidor observado. En la tercera parte de este trabajo, se ha evaluado el papel de los preparados comerciales de LSI en el aroma y en las características organolépticas de los vinos. Para ello, empleando vinos sintéticos, se estudiaron por un lado, el efecto que los componentes cedidos por los preparados de LSI podrían tener en la volatilidad de compuestos del aroma ya presentes en el vino, y por otro lado, en la posible cesión al vino de compuestos odorantes presentes en los preparados y originados durante su fabricación. Los resultados de estos estudios pusieron de manifiesto que los preparados de LSI pueden modificar la volatilidad de compuestos del aroma del vino, dependiendo principalmente de las características físico-químicas del compuesto del aroma, del tiempo de permanencia del preparado en el vino y del tipo de preparado de LSI. Por otro lado, gracias al empleo de diferentes técnicas de extracción combinadas, la extracción con fluidos presurizados (PLE) y microextracción en fase sólida en el espacio de cabeza (HS-SPME) y posterior análisis por GO-MS, se caracterizó el perfil volátil de un preparado de LSI representativo del resto de los estudiados y se comprobó que estaba formado tanto por componentes procedentes de la levadura, como ácidos grasos de cadena corta y larga, así como por compuestos volátiles originados durante las etapas de fabricación de los preparados y resultantes de la reacción de Maillard. Se determinó a su vez, que algunos de estos compuestos volátiles, pueden ser cedidos tanto a medios vínicos, como a vinos reales. Concretamente, se comprobó que los preparados indicados para vinificación en blanco liberan una menor concentración de compuestos volátiles, y especialmente, de productos de la reacción de Maillard, que los indicados para vinificación en tinto. Estos compuestos, y principalmente los pertenecientes al grupo de las pirazinas, podrían modificar el perfil sensorial de los vinos, ya que son compuestos con, en general, bajos umbrales de percepción. Para intentar eliminar estos compuestos de los preparados, es decir, para su desodorización, en esta parte del trabajo, también se desarrolló una metodología basada en la extracción con CO2 supercrítico, que permitió la eliminación selectiva de este tipo de compuestos sin provocar modificaciones en su composición no volátil. Debido a la cantidad de preparados de LSI ricos en glutatión (GSH) disponibles en el mercado para mejorar las características sensoriales de los vinos, se estudió la cesión de este compuesto a medios vínicos. Para ello, se puso a punto una metodología que permitiera la cuantificación del GSH en sus diferentes formas (reducido y total), en vinos sintéticos y reales. El método se aplicó para cuantificar la cantidad de GSH aportada por los LSI a los vinos. Los resultados mostraron que los preparados de LSI indicados para vinificación en blanco y rosado liberaron, inmediatamente tras su adición a medios vínicos sintéticos, ♯y-glutamil-cisteína y GSH principalmente en su forma reducida a mayores concentraciones que los preparados indicados para otras aplicaciones. Por otro lado, se aislaron las fracciones de peso molecular <3000 Da de preparados de LSI y se emplearon para suplementar vinos sintéticos aromatizados con terpenos y sometidos a condiciones de envejecimiento acelerado. Los resultados revelaron una reducción en la oxidación por parte de las fracciones aisladas tanto de preparados de LSI ricos en glutatión como de preparados empleados como nutrientes de fermentación. Esta reducción era incluso mayor que cuando se empleaba el glutatión comercial, lo que indicaba la presencia de otros compuestos con propiedades antioxidantes. Finalmente, para conocer el efecto de los preparados de LSI a escala de bodega, se elaboraron a escala industrial dos tipos de vinos, con y sin la adición de un preparado de LSI rico en GSH, a escala industrial. En estos vinos, se determino el contenido en glutatión, y la evolución del color y de la composición fenólica y volátil durante la vinificación y a lo largo del envejecimiento en botella. Los resultados de este trabajo mostraron la capacidad del preparado de LSI de ceder glutatión en su forma reducida, aunque esté es rápidamente oxidado durante la fermentación alcohólica. Sin embargo, se comprobó que los vinos elaborados con el preparado presentaron mayor estabilidad de color y aroma a lo largo del envejecimiento en botella. De hecho, a los 9 meses de envejecimiento, se encontraron diferencias sensoriales entre los vinos con y sin el preparado. Sin embargo, estas diferencias sensoriales no fueron suficientes para que en el estudio de preferencia que se llevó a cabo con los mismos vinos, los consumidores mostraran preferencia hacia los vinos elaborados con el preparado de LSI.

Urbanismo aerodinámico: entre lo sólido y lo fluido

Dentro de las variables meteorológicas menos estudiadas en relación a las ciudades se encuentra el viento. Las constantes variaciones de cambio a nivel de gradiente vertical, de la meso-escala a la micro-escala junto a los cambios a nivel del cañón urbano en velocidad, dirección y fuerza lo vuelven complejo de investigar. Para poder estudiar el viento se planteó una estrategia metodológica que abordará la multidimensionalidad del fenómeno del viento, aplicándolo en la ciudad más austral del mundo, la ciudad de Punta Arenas donde los vientos no solo alcanzar los 183 Km por hora sino que además vienen desde la misma dirección en más de un 90% de las veces. Este comportamiento constante en la direccionalidad del viento permite reconocer los fenómenos aerodinámicos en primera instancia producidos por la rugosidad urbana en sus líneas de viento y luego a nivel morfológico en el comportamiento aerodinámico del fluido laminar sobre los cuerpos edificados, generándose un patrón regular entre lo sólidos y lo fluidos. Los fuertes vientos en ciudades de climas fríos inciden finalmente en el uso estancial de sus espacios públicos, los cuales al no estar diseñados bajo estas condiciones no desarrollan apropiados niveles de confort térmico de sus espacios exteriores.

Estudio matemático del comportamiento de un fluido incompresible en flujo laminar

El objetivo de este trabajo es realizar un estudio del comportamiento de un fluido en flujo laminar. En primer lugar se realizará una breve introducción Histórica, con el fin de situar al lector en la época donde suceden los acontecimientos más importantes relativos a la mecánica de fluidos. Se presentarán a las figuras más importantes de la mecánica de fluidos moderna, como Sir Gabriel Stokes o Osborne Reynolds, así como los acontecimientos más importantes. A continuación se expondran las ecuaciones que rigen el comportamiento de un fluido, las ecuaciones de Navier-Stokes y la ecuación de Reynolds, con el objetivo de ayudar al lector a entender los análisis posteriores que se realizarán sobre dichas ecuaciones. En tercer lugar se analizará la ecuación de Reynolds, la existencia y unicidad de soluciones, para seguidamente, realizar una simulación del problema. Dicha simulación se ha realizado en un script, bajo la herramienta Matlab, se explicará como se ha realizado la simulación y se expondrán ejemplos de los casos de un fluido compresible y uno incompresible para diversas ecuaciones. Seguidamente se describe el fenómeno de la capa límite y se explica con un ejemplo en un fluido compresible. Para finalizar se exponen las conclusiones obtenidas en la realización de este trabajo.

NURBS-Based Geometry Parameterization for Aerodynamic Shape Optimization

La motivación de esta tesis es el desarrollo de una herramienta de optimización automática para la mejora del rendimiento de formas aerodinámicas enfocado en la industria aeronáutica. Este trabajo cubre varios aspectos esenciales, desde el empleo de Non-Uniform Rational B-Splines (NURBS), al cálculo de gradientes utilizando la metodología del adjunto continuo, el uso de b-splines volumétricas como parámetros de diseño, el tratamiento de la malla en las intersecciones, y no menos importante, la adaptación de los algoritmos de la dinámica de fluidos computacional (CFD) en arquitecturas hardware de alto paralelismo, como las tarjetas gráficas, para acelerar el proceso de optimización. La metodología adjunta ha posibilitado que los métodos de optimización basados en gradientes sean una alternativa prometedora para la mejora de la eficiencia aerodinámica de los aviones. La formulación del adjunto permite calcular los gradientes de una función de coste, como la resistencia aerodinámica o la sustentación, independientemente del número de variables de diseño, a un coste computacional equivalente a una simulación CFD. Sin embargo, existen problemas prácticos que han imposibilitado su aplicación en la industria, que se pueden resumir en: integrabilidad, rendimiento computacional y robustez de la solución adjunta. Este trabajo aborda estas contrariedades y las analiza en casos prácticos. Como resumen, las contribuciones de esta tesis son: • El uso de NURBS como variables de diseño en un bucle de automático de optimización, aplicado a la mejora del rendimiento aerodinámico de alas en régimen transónico. • El desarrollo de algoritmos de inversión de punto, para calcular las coordenadas paramétricas de las coordenadas espaciales, para ligar los vértices de malla a las NURBS. • El uso y validación de la formulación adjunta para el calculo de los gradientes, a partir de las sensibilidades de la solución adjunta, comparado con diferencias finitas. • Se ofrece una estrategia para utilizar la geometría CAD, en forma de parches NURBS, para tratar las intersecciones, como el ala-fuselaje. • No existen muchas alternativas de librerías NURBS viables. En este trabajo se ha desarrollado una librería, DOMINO NURBS, y se ofrece a la comunidad como código libre y abierto. • También se ha implementado un código CFD en tarjeta gráfica, para realizar una valoración de cómo se puede adaptar un código sobre malla no estructurada a arquitecturas paralelas. • Finalmente, se propone una metodología, basada en la función de Green, como una forma eficiente de paralelizar simulaciones numéricas. Esta tesis ha sido apoyada por las actividades realizadas por el Área de Dinámica da Fluidos del Instituto Nacional de Técnica Aeroespacial (INTA), a través de numerosos proyectos de financiación nacional: DOMINO, SIMUMAT, y CORESFMULAERO. También ha estado en consonancia con las actividades realizadas por el departamento de Métodos y Herramientas de Airbus España y con el grupo Investigación y Tecnología Aeronáutica Europeo (GARTEUR), AG/52. ABSTRACT The motivation of this work is the development of an automatic optimization strategy for large scale shape optimization problems that arise in the aeronautics industry to improve the aerodynamic performance; covering several aspects from the use of Non-Uniform Rational B-Splines (NURBS), the calculation of the gradients with the continuous adjoint formulation, the development of volumetric b-splines parameterization, mesh adaptation and intersection handling, to the adaptation of Computational Fluid Dynamics (CFD) algorithms to take advantage of highly parallel architectures in order to speed up the optimization process. With the development of the adjoint formulation, gradient-based methods for aerodynamic optimization become a promising approach to improve the aerodynamic performance of aircraft designs. The adjoint methodology allows the evaluation the gradients to all design variables of a cost function, such as drag or lift, at the equivalent cost of more or less one CFD simulation. However, some practical problems have been delaying its full implementation to the industry, which can be summarized as: integrability, computer performance, and adjoint robustness. This work tackles some of these issues and analyse them in well-known test cases. As summary, the contributions comprises: • The employment of NURBS as design variables in an automatic optimization loop for the improvement of the aerodynamic performance of aircraft wings in transonic regimen. • The development of point inversion algorithms to calculate the NURBS parametric coordinates from the space coordinates, to link with the computational grid vertex. • The use and validation of the adjoint formulation to calculate the gradients from the surface sensitivities in an automatic optimization loop and evaluate its reliability, compared with finite differences. • This work proposes some algorithms that take advantage of the underlying CAD geometry description, in the form of NURBS patches, to handle intersections and mesh adaptations. • There are not many usable libraries for NURBS available. In this work an open source library DOMINO NURBS has been developed and is offered to the community as free, open source code. • The implementation of a transonic CFD solver from scratch in a graphic card, for an assessment of the implementability of conventional CFD solvers for unstructured grids to highly parallel architectures. • Finally, this research proposes the use of the Green's function as an efficient paralellization scheme of numerical solvers. The presented work has been supported by the activities carried out at the Fluid Dynamics branch of the National Institute for Aerospace Technology (INTA) through national founding research projects: DOMINO, SIMUMAT, and CORESIMULAERO; in line with the activities carried out by the Methods and Tools and Flight Physics department at Airbus and the Group for Aeronautical Research and Technology in Europe (GARTEUR) action group AG/52.

Urban wind energy: empirical optimization of high-rise building roof shape for the wind energy exploitation

El programa Europeo HORIZON2020 en Futuras Ciudades Inteligentes establece como objetivo que el 20% de la energía eléctrica sea generada a partir de fuentes renovables. Este objetivo implica la necesidad de potenciar la generación de energía eólica en todos los ámbitos. La energía eólica reduce drásticamente las emisiones de gases de efecto invernadero y evita los riesgos geo-políticos asociados al suministro e infraestructuras energéticas, así como la dependencia energética de otras regiones. Además, la generación de energía distribuida (generación en el punto de consumo) presenta significativas ventajas en términos de elevada eficiencia energética y estimulación de la economía. El sector de la edificación representa el 40% del consumo energético total de la Unión Europea. La reducción del consumo energético en este área es, por tanto, una prioridad de acuerdo con los objetivos "20-20-20" en eficiencia energética. La Directiva 2010/31/EU del Parlamento Europeo y del Consejo de 19 de mayo de 2010 sobre el comportamiento energético de edificaciones contempla la instalación de sistemas de suministro energético a partir de fuentes renovables en las edificaciones de nuevo diseño. Actualmente existe una escasez de conocimiento científico y tecnológico acerca de la geometría óptima de las edificaciones para la explotación de la energía eólica en entornos urbanos. El campo tecnológico de estudio de la presente Tesis Doctoral es la generación de energía eólica en entornos urbanos. Específicamente, la optimization de la geometría de las cubiertas de edificaciones desde el punto de vista de la explotación del recurso energético eólico. Debido a que el flujo del viento alrededor de las edificaciones es exhaustivamente investigado en esta Tesis empleando herramientas de simulación numérica, la mecánica de fluidos computacional (CFD en inglés) y la aerodinámica de edificaciones son los campos científicos de estudio. El objetivo central de esta Tesis Doctoral es obtener una geometría de altas prestaciones (u óptima) para la explotación de la energía eólica en cubiertas de edificaciones de gran altura. Este objetivo es alcanzado mediante un análisis exhaustivo de la influencia de la forma de la cubierta del edificio en el flujo del viento desde el punto de vista de la explotación energética del recurso eólico empleando herramientas de simulación numérica (CFD). Adicionalmente, la geometría de la edificación convencional (edificio prismático) es estudiada, y el posicionamiento adecuado para los diferentes tipos de aerogeneradores es propuesto. La compatibilidad entre el aprovechamiento de las energías solar fotovoltaica y eólica también es analizado en este tipo de edificaciones. La investigación prosigue con la optimización de la geometría de la cubierta. La metodología con la que se obtiene la geometría óptima consta de las siguientes etapas: - Verificación de los resultados de las geometrías previamente estudiadas en la literatura. Las geometrías básicas que se someten a examen son: cubierta plana, a dos aguas, inclinada, abovedada y esférica. - Análisis de la influencia de la forma de las aristas de la cubierta sobre el flujo del viento. Esta tarea se lleva a cabo mediante la comparación de los resultados obtenidos para la arista convencional (esquina sencilla) con un parapeto, un voladizo y una esquina curva. - Análisis del acoplamiento entre la cubierta y los cerramientos verticales (paredes) mediante la comparación entre diferentes variaciones de una cubierta esférica en una edificación de gran altura: cubierta esférica estudiada en la literatura, cubierta esférica integrada geométricamente con las paredes (planta cuadrada en el suelo) y una cubierta esférica acoplada a una pared cilindrica. El comportamiento del flujo sobre la cubierta es estudiado también considerando la posibilidad de la variación en la dirección del viento incidente. - Análisis del efecto de las proporciones geométricas del edificio sobre el flujo en la cubierta. - Análisis del efecto de la presencia de edificaciones circundantes sobre el flujo del viento en la cubierta del edificio objetivo. Las contribuciones de la presente Tesis Doctoral pueden resumirse en: - Se demuestra que los modelos de turbulencia RANS obtienen mejores resultados para la simulación del viento alrededor de edificaciones empleando los coeficientes propuestos por Crespo y los propuestos por Bechmann y Sórensen que empleando los coeficientes estándar. - Se demuestra que la estimación de la energía cinética turbulenta del flujo empleando modelos de turbulencia RANS puede ser validada manteniendo el enfoque en la cubierta de la edificación. - Se presenta una nueva modificación del modelo de turbulencia Durbin k — e que reproduce mejor la distancia de recirculación del flujo de acuerdo con los resultados experimentales. - Se demuestra una relación lineal entre la distancia de recirculación en una cubierta plana y el factor constante involucrado en el cálculo de la escala de tiempo de la velocidad turbulenta. Este resultado puede ser empleado por la comunidad científica para la mejora del modelado de la turbulencia en diversas herramientas computacionales (OpenFOAM, Fluent, CFX, etc.). - La compatibilidad entre las energías solar fotovoltaica y eólica en cubiertas de edificaciones es analizada. Se demuestra que la presencia de los módulos solares provoca un descenso en la intensidad de turbulencia. - Se demuestran conflictos en el cambio de escala entre simulaciones de edificaciones a escala real y simulaciones de modelos a escala reducida (túnel de viento). Se demuestra que para respetar las limitaciones de similitud (número de Reynolds) son necesarias mediciones en edificaciones a escala real o experimentos en túneles de viento empleando agua como fluido, especialmente cuando se trata con geometrías complejas, como es el caso de los módulos solares. - Se determina el posicionamiento más adecuado para los diferentes tipos de aerogeneradores tomando en consideración la velocidad e intensidad de turbulencia del flujo. El posicionamiento de aerogeneradores es investigado en las geometrías de cubierta más habituales (plana, a dos aguas, inclinada, abovedada y esférica). - Las formas de aristas más habituales (esquina, parapeto, voladizo y curva) son analizadas, así como su efecto sobre el flujo del viento en la cubierta de un edificio de gran altura desde el punto de vista del aprovechamiento eólico. - Se propone una geometría óptima (o de altas prestaciones) para el aprovechamiento de la energía eólica urbana. Esta optimización incluye: verificación de las geometrías estudiadas en el estado del arte, análisis de la influencia de las aristas de la cubierta en el flujo del viento, estudio del acoplamiento entre la cubierta y las paredes, análisis de sensibilidad del grosor de la cubierta, exploración de la influencia de las proporciones geométricas de la cubierta y el edificio, e investigación del efecto de las edificaciones circundantes (considerando diferentes alturas de los alrededores) sobre el flujo del viento en la cubierta del edificio objetivo. Las investigaciones comprenden el análisis de la velocidad, la energía cinética turbulenta y la intensidad de turbulencia en todos los casos. ABSTRACT The HORIZON2020 European program in Future Smart Cities aims to have 20% of electricity produced by renewable sources. This goal implies the necessity to enhance the wind energy generation, both with large and small wind turbines. Wind energy drastically reduces carbon emissions and avoids geo-political risks associated with supply and infrastructure constraints, as well as energy dependence from other regions. Additionally, distributed energy generation (generation at the consumption site) offers significant benefits in terms of high energy efficiency and stimulation of the economy. The buildings sector represents 40% of the European Union total energy consumption. Reducing energy consumption in this area is therefore a priority under the "20-20-20" objectives on energy efficiency. The Directive 2010/31/EU of the European Parliament and of the Council of 19 May 2010 on the energy performance of buildings aims to consider the installation of renewable energy supply systems in new designed buildings. Nowadays, there is a lack of knowledge about the optimum building shape for urban wind energy exploitation. The technological field of study of the present Thesis is the wind energy generation in urban environments. Specifically, the improvement of the building-roof shape with a focus on the wind energy resource exploitation. Since the wind flow around buildings is exhaustively investigated in this Thesis using numerical simulation tools, both computational fluid dynamics (CFD) and building aerodynamics are the scientific fields of study. The main objective of this Thesis is to obtain an improved (or optimum) shape of a high-rise building for the wind energy exploitation on the roof. To achieve this objective, an analysis of the influence of the building shape on the behaviour of the wind flow on the roof from the point of view of the wind energy exploitation is carried out using numerical simulation tools (CFD). Additionally, the conventional building shape (prismatic) is analysed, and the adequate positions for different kinds of wind turbines are proposed. The compatibility of both photovoltaic-solar and wind energies is also analysed for this kind of buildings. The investigation continues with the buildingroof optimization. The methodology for obtaining the optimum high-rise building roof shape involves the following stages: - Verification of the results of previous building-roof shapes studied in the literature. The basic shapes that are compared are: flat, pitched, shed, vaulted and spheric. - Analysis of the influence of the roof-edge shape on the wind flow. This task is carried out by comparing the results obtained for the conventional edge shape (simple corner) with a railing, a cantilever and a curved edge. - Analysis of the roof-wall coupling by testing different variations of a spherical roof on a high-rise building: spherical roof studied in the litera ture, spherical roof geometrically integrated with the walls (squared-plant) and spherical roof with a cylindrical wall. The flow behaviour on the roof according to the variation of the incident wind direction is commented. - Analysis of the effect of the building aspect ratio on the flow. - Analysis of the surrounding buildings effect on the wind flow on the target building roof. The contributions of the present Thesis can be summarized as follows: - It is demonstrated that RANS turbulence models obtain better results for the wind flow around buildings using the coefficients proposed by Crespo and those proposed by Bechmann and S0rensen than by using the standard ones. - It is demonstrated that RANS turbulence models can be validated for turbulent kinetic energy focusing on building roofs. - A new modification of the Durbin k — e turbulence model is proposed in order to obtain a better agreement of the recirculation distance between CFD simulations and experimental results. - A linear relationship between the recirculation distance on a flat roof and the constant factor involved in the calculation of the turbulence velocity time scale is demonstrated. This discovery can be used by the research community in order to improve the turbulence modeling in different solvers (OpenFOAM, Fluent, CFX, etc.). - The compatibility of both photovoltaic-solar and wind energies on building roofs is demonstrated. A decrease of turbulence intensity due to the presence of the solar panels is demonstrated. - Scaling issues are demonstrated between full-scale buildings and windtunnel reduced-scale models. The necessity of respecting the similitude constraints is demonstrated. Either full-scale measurements or wind-tunnel experiments using water as a medium are needed in order to accurately reproduce the wind flow around buildings, specially when dealing with complex shapes (as solar panels, etc.). - The most adequate position (most adequate roof region) for the different kinds of wind turbines is highlighted attending to both velocity and turbulence intensity. The wind turbine positioning was investigated for the most habitual kind of building-roof shapes (flat, pitched, shed, vaulted and spherical). - The most habitual roof-edge shapes (simple edge, railing, cantilever and curved) were investigated, and their effect on the wind flow on a highrise building roof were analysed from the point of view of the wind energy exploitation. - An optimum building-roof shape is proposed for the urban wind energy exploitation. Such optimization includes: state-of-the-art roof shapes test, analysis of the influence of the roof-edge shape on the wind flow, study of the roof-wall coupling, sensitivity analysis of the roof width, exploration of the aspect ratio of the building-roof shape and investigation of the effect of the neighbouring buildings (considering different surrounding heights) on the wind now on the target building roof. The investigations comprise analysis of velocity, turbulent kinetic energy and turbulence intensity for all the cases.

Simulación tridimensional de una cámara de combustión experimental para obtener los campos de velocidad, presión y temperatura, mediante el software CFX 5.6

El objetivo de la presente investigación es predecir los campos de velocidad, presión y temperatura en una cámara de combustión experimental, mediante la técnica de la simulación numérica de flujo de fluidos. Para ello se revisa el procedimiento de solución numérica de las ecuaciones de transporte, aplicadas a la cámara de combustión experimental. La simulación está basada en el software CFX 5.6, el cual fue adquirido por la universidad nacional experimental del Táchira por medio del Decanato de Post-grado y de Investigación. Se hace un estudio de la sensibilidad de malla para adecuar el criterio de convergencia que el software requiere. La cámara de combustión experimental empleada para éste estudio es una cámara de combustión diseñada por estudiantes de Pre-grado para determinar la temperatura de flama adiabática, aunque el diseño de esta cámara no es estándar, es útil para medir la temperatura en tiempo real. El combustible empleado para éste análisis es propano (C3H8) el cual es inyectado a la cámara de combustión por una tubería concéntrica al flujo de aire. En la solución de la simulación computacional se aprecia, a través del perfil de temperatura, la envolvente de la llama, formada por el contorno de temperatura máxima, la cual es similar a la observada en cualquier cámara de combustión.

Simulación transitoria de una mezcla no reactiva de aire, combustible y gases de recirculación dentro de una precámara de combustión de un motor

El objetivo del presente trabajo fue predecir, en función del tiempo, los campos de velocidad, presión, temperatura y concentración de especies químicas de una mezcla no reactiva de aire, gases de recirculación y combustible (iso-octano), en una pre-cámara de combustión de un motor CFR, mediante la técnica de la simulación numérica de flujo de fluidos. Para ello se revisó el procedimiento de la solución numérica de las ecuaciones de transporte, aplicadas a la pre-cámara de combustión del motor CFR. La simulación se basó en el software CFX 10. Se hizo un estudio de la sensibilidad de malla para adecuar el criterio de·convergencia que el software requiere. La pre-cámara de combustión empleada para este estudio fue una pre-cámara de combustión modificada de un motor CFR para hacer que el mismo funcione como un motor HCCI. El motor CFR está ubicado en la Universidad Politécnica de Madrid.

Pattern formation of parametrically driven surface waves induced by vibration

Esta tesis se centra en la generación de ondas superficiales subarmónicas en fluidos sometidos a vibración forzada en el régimen gravitatorio capilar con líquidos de baja viscosidad. Tres problemas diferentes han sido estudiados: un contenedor rectangular con vibración horizontal, la misma geometría pero con una combinación de vibración vertical y horizontal y un obstáculo completamente sumergido vibrado verticalmente en un contenedor grande. Se deriva una ecuación de amplitud desde primeros principios para describir las ondas subarmónicas con forzamiento parámetrico inducido por la vibración. La ecuación es bidimensional mientras que el problema original es tridimensional y admite un forzamiento espacial no uniforme. Usando esta ecuación los tres sistemas han sido analizados, centrándose en calcular la amplitud crítica, la orientación de los patrones y el carácter temporal de los patrones espaciotemporales, que pueden ser estrictamente subarmónicos o cuasiperiodicos con una frecuencia de modulación temporal. La dependencia con los parámetros adimensionales también se considera. La teoría será comparada con los experimentos disponibles en la literatura. Abstract This thesis focus on the generation of subharmonic surface waves on fluids subject to forced vibration in the gravity-capillary regime with liquids of small viscosity. Three different problems have been considered: a rectangular container under horizontal vibration; the same geometry but under a combination of horizontal and vertical vibration; and a fully submerged vertically vibrated obstacle in a large container. An amplitude equation is derived from first principles that fairly precisely describes the subharmonic surfaces waves parametrically driven by vibration. That equation is two dimensional while the underlying problem is three-dimensional and permits spatially nonuniform forcing. Using this equation, the three systems have been analyzed, focusing on the calculation of the threshold amplitude, the pattern orientation, and the temporal character of the spatio-temporal patterns, which can be either strictly subharmonic or quasi-periodic, showing an additional modulation frequency. Dependence on the non-dimensional parameters is also considered. The theory is compared with the experiments available in the literature.