Optimización del peaje en autopistas interurbanas que compiten con carreteras convencionales libres de pago : una perspectiva de eficiencia y equidad

En algunos países, como por ejemplo España, es común que entre un origen y un destino existan dos carreteras paralelas en las que existen ciertas diferencias. La más importante es que una de las vías es una autopista que ofrece a los usuarios una mayor comodidad y un menor tiempo de viaje a cambio del pago de un peaje, el cual no es necesario abonar en la carretera convencional paralela. Así, el problema de la tarificación vial ha sido estudiado en diversas ocasiones. Existe un amplio consenso en que para lograr el máximo bienestar social, los usuarios deben internalizar las externalidades que producen y no perciben a través de un peaje. Sin embargo, dicho peaje puede perjudicar a los usuarios con bajos ingresos. Dependiendo del objetivo (por ejemplo, maximizar el bienestar, maximizar la equidad social, la amortización de la construcción de la carretera, etc) el peaje óptimo podría variar sustancialmente. La literatura académica acerca de los peajes, la eficiencia y la equidad es vasta y diversa. Sin embargo, hemos encontrado una deficiencia en dicha literatura sobre el peaje óptimo, en corredores donde una carretera y una autopista con diferentes características de calidad compiten para capturar el tráfico. Particularmente no se ha encontrado ninguna investigación acerca del establecimiento del peaje que maximice el bienestar social o la equidad para distintas distribuciones del valor del tiempo de viaje (VTT), caracterizadas por su media y varianza. Por ello, el principal objetivo de la investigación es estimar la influencia que tiene la distribución de la renta de una sociedad sobre el peaje óptimo. La presente tesis doctoral trata de obtener, por medio de una metodología robusta, las diferentes políticas de peajes que los planificadores de transporte deberían llevar a cabo según la riqueza y cohesión social de los potenciales usuarios del corredor, la demanda y el objetivo que se busque con dicha tarificación, esto es, maximizar el bienestar social o la equidad. Adicionalmente también se obtienen los peajes óptimos dependiendo de si el corredor se encuentra totalmente tarificado o únicamente se debe pagar un peaje por circular en la autopista. In some countries, such as Spain, it is very common that in the same corridor there are two roads with the same origin and destination but with some differences. The most important contrast is that one is a toll highway which offers a better quality than the parallel road in exchange of a price. The users decide if the price of the toll worth to pay for the advantages offered. The problem of road pricing has been largely studied. It is well acknowledged that in order to achieve the maximum social welfare, users must internalize the externalities they produce and do not perceive through a toll. However, that toll can harm users with low income. Depending on the objective (e.g. maximize welfare, maximize social equity, amortize the construction of the road, etc) the optimal toll might vary substantially. The academic literature about pricing, efficiency and equity is vast and diverse. However, as far as we have found, there is a gap in the literature regarding the optimal price where a road and a highway with different quality characteristics compete for capturing the traffic in a corridor. Particularly we did not find any research estimating the optimal welfare price or the optimal equity price for different Value of Travel Time (VTT) distributions characterized by different VTT average and variance. The objective of the research is to fill this gap. In this research a theoretical model in order to obtain the optimal price in a toll highway that competes for capturing the traffic with a conventional road is developed. This model is done from the welfare and equity perspective and for non‐usual users who decide over the expectation of free flow conditions. The model is finally applied to the variables we want to focus on: average value of travel time (VTT) which is strongly related with income, dispersion of this VTT, different kind of distributions of VTT and traffic levels, from free flow to congestion. Furthermore, we also obtain the optimal tolls with the corridor completely charged or with untolled alternative.

Estudio paramétrico del conjunto tractor-trailer mediante aplicación de mecánica de fluidos computacional

La mejora del comportamiento aerodinámico de vehículos pesados en la carretera ha adquirido una mayor importancia estos últimos años debido a la crisis y su consecuente aumento del precio de los combustibles. Dado que una reducción en la resistencia aerodinámica del vehículo conlleva un menor empleo de combustible, el objetivo de este proyecto es realizar un estudio paramétrico del conjunto tractor-tráiler mediante la mecánica de fluidos computacional (CFD) para así obtener la geometría que proporciona un menor gasto de combustible cuando ésta es expuesta a viento frontal. La influencia de 3 parámetros, que son la separación existente entre cabina y remolque, la altura del remolque y el radio de curvatura de las aristas frontales de la cabina, es analizada en este estudio dividido en varias etapas que implican el uso de programas específicos como son: la parametrización y creación de las geometrías en 3D que es llevada a cabo mediante CATIA, el mallado del dominio (realizado con Gambit) y resolución de las ecuaciones mediante FLUENT. Finalmente se obtendrá una relación entre la resistencia aerodinámica (representada mediante el coeficiente de arrastre) y la combinación de los 3 parámetros, que nos permitirá decidir que geometría es la óptima. Abstract The improvement of the heavy vehicle’s aerodynamic behavior on the road has gained a great importance for these last years because of the economic crisis and the consequent increase of the price of the fuels. Due to the fact that a reduction in the aerodynamic resistance of the vehicle involves using a smaller amount of fuel, the objective of this project is to carry out a parametric study about the ensemble tractor-trailer by computational fluid dynamics methods (CFD) in order to obtain the geometry which expense of fuel is the lowest when it’s exposed to frontal wind. The influence of the three parameters, which are the space between cab and trailer, the height of the trailer and the curvature of the frontal cab edges, is analysed in this study which is divided into different parts involving specific programs: choosing the parameters and building the geometries, which is done by using CATIA, the mesh is built by Gambit, and the program equations-solver is FLUENT. Finally a ratio between aerodynamic resistance and a combination of the three parameters will be obtained and it will allow us to choose the best geometry.

Desarrollo de familias jerárquicas de elementos finitos de clase CK

Existen numerosas situaciones en la Técnica para las que es preciso resolver problemas de condiciones iniciales y contorno con una elevada exigencia de continuidad en las soluciones. Algunos ejemplos no exhaustivos se citan a continuación: flexión de vigas y placas en el análisis de las estructuras, problemas de láminas, topografía, trazado de vías de comunicación, definición geométrica de estructuras o reconocimiento caligráfico. Para ello se utilia el método de los elementos finitos que constituye en la actualidad una técnica matemática de discretización bien establecida . Su eficiencia se manifiesta en particular en la resolución de problemas de contorno planteados en su formulación débil y en los problemas de determinación de extrema les de funcionales.

Un modelo no paramétrico de evaluación de la eficiencia y la gestión de las redes sociales virtuales : una aplicación a las empresas del sector de las telecomunicaciones en España

Este artículo analiza la relación entre la eficiencia productiva y las Redes Sociales Virtuales (RSV) en las empresas de telecomunicaciones en España. En una primera etapa, se aplica el análisis envolvente de datos (DEA) incorporando varios indicadores de actividad ?Social Media?. En una segunda etapa, se utiliza una regresión logística para caracterizar las empresas eficientes. Los resultados muestran que la capacidad de absorción y utilización de las RSV es un factor determinante en la mejora de la eficiencia productiva. La utilización combinada y las distintas capacidades de gestión de las RSV permiten identificarlas como un recurso heterogéneo. Este trabajo presenta un modelo para la evaluación del desempeño estratégico al abordar su presencia y actividad en RSV. ABSTRACT. This paper analyzes the relationship between the productive efficiency and the Online Social Networks - OSN in the Spanish telecommunications firms. First, a data envelopment analysis (DEA) is used and several indicators of business "Social Media" activities are incorporated. In a second stage, a logistic regression model regression is applied to characteri ze the efficient enterprises. Results show that the company's ability to absorb and utilize this OSN is a key factor in improving the productive efficiency. These results on the combined use and different management capabilities of OSN point to a definitio n of OSN as a heterogeneous resource. This paper presents a model for assessing the strategic performance to address their presence and activity in OSN.

Dos propuestas legislativas para la eficiencia de los parques tecnológicos.

Se plantea la conveniencia de una norma con rango de Directiva comunitaria que limite la competencia entre administraciones públicas de áreas geográficas diversas para atraer la implantación de empresas en sus parques tecnológicos. También se expone la necesidad de desarrollar instrumentos de inversión colectiva que favorezcan a la PYME. Con estos fines, se conceptúa el significado de parque tecnológico, a la vez que se especula sobre sus principales efectos en el tejido industrial y sobre la actividad investigadora de la empresa.

Caracterización mecánica de tablas de madera de pino silvestre, usadas como barandillas en sistemas provisionales de protección de borde

El empleo de piezas de madera de pino silvestre de pequeña escuadría es una práctica muy habitual en el sector de la Construcción en la zona central de la Península Ibérica. Existen numerosos estudios sobre este tipo de material pero todos ellos se han realizado sobre piezas de gran escuadría. No obstante, se constata un vacío absoluto de informes técnicos en lo que se refiere al comportamiento mecánico de elementos de pequeña escuadría, que son los que normalmente se utilizan en protecciones colectivas (PC) y medios auxiliares de edificación (MAE). La norma UNE-EN 13374 “Sistemas provisionales de protección de borde. Especificaciones del producto, métodos de ensayo” especifica que la madera empleada en sistemas provisionales de protección de borde debe poseer como mínimo la clase resistente C14. La obtención de la clase resistente de la madera empleada en construcción para PC y MAE mediante el empleo de la norma UNE 56544 “Clasificación visual de la madera aserrada para uso estructural. Maderas de coníferas ”, dejaría fuera de uso prácticamente a todos los elementos empleados en la práctica, que suelen ser de madera de cuarta, a la que no se le puede asignar una clase resistente y además la norma supone que su resistencia es inferior a la de la clase C18 en elementos de pino silvestre. En este trabajo se ha realizado una caracterización mecánica de tablas de madera de pino silvestre de las diferentes calidades que existen en el mercado (especial, primera, segunda, tercera y cuarta) en los grosores comerciales 22, 27, 30 y 40 mm y ancho 150 mm., que son los utilizados habitualmente en (PC) y (MAE). En todos los elementos se ha realizado una clasificación visual y resistente previa según la norma UNE 56544. Posteriormente se han ensayado las piezas para determinar sus características mecánicas reales. Se han comparado los resultados con los obtenidos por otros investigadores sobre piezas de gran escuadría, buscando las posibles correlaciones entre algunas de sus propiedades.

Aceros al molibdeno sinterizados con gradiente funcional

Las propiedades que se solicitan a las piezas y productos exigen, la mayoría de las veces, una combinación de materiales que respondan a requerimientos específicos en función de su posición. La fabricación de materiales compuestos y materiales con gradiente funcional proporciona grandes posibilidades de respuesta a las necesidades demandas, tanto de satisfacción de requerimientos como de reducción de materias primas necesarias. En este sentido, satisfacer las propiedades demandadas, de manera localizada geométricamente, puede suponer un ahorro de materias primas y de los residuos relacionados en su fabricación. En el presente estudio de la sinterización a 1120 ºC de un compacto de dos aceros de diferentes composiciones (un acero con un 0,3 % C y un 2 % de cobre y un acero al molibdeno con un 0,7 % carbono y un 2 % de cobre) ofrece unas posibilidades muy atractivas de eficiencia, al reducir en el interior la composición que se requiere en la superficie de las piezas. Se han caracterizado las propiedades mecánicas de este acero sinterizado con gradiente funcional: densidad, variación dimensional, durezas y resistencia mecánica; ilustradas con el análisis de imagen mediante microscopía óptica y electrónica de barrido.

Multiscale modeling of the plastic behaviour in single crystal tungsten: from atomistic to crystal plasticity simulations

En una planta de fusión, los materiales en contacto con el plasma así como los materiales de primera pared experimentan condiciones particularmente hostiles al estar expuestos a altos flujos de partículas, neutrones y grandes cargas térmicas. Como consecuencia de estas diferentes y complejas condiciones de trabajo, el estudio, desarrollo y diseño de estos materiales es uno de los más importantes retos que ha surgido en los últimos años para la comunidad científica en el campo de los materiales y la energía. Debido a su baja tasa de erosión, alta resistencia al sputtering, alta conductividad térmica, muy alto punto de fusión y baja retención de tritio, el tungsteno (wolframio) es un importante candidato como material de primera pared y como posible material estructural avanzado en fusión por confinamiento magnético e inercial. Sin embargo, el tiempo de vida del tungsteno viene controlado por diversos factores como son su respuesta termo-mecánica en la superficie, la posibilidad de fusión y el fallo por acumulación de helio. Es por ello que el tiempo de vida limitado por la respuesta mecánica del tungsteno (W), y en particular su fragilidad, sean dos importantes aspectos que tienes que ser investigados. El comportamiento plástico en materiales refractarios con estructura cristalina cúbica centrada en las caras (bcc) como el tungsteno está gobernado por las dislocaciones de tipo tornillo a escala atómica y por conjuntos e interacciones de dislocaciones a escalas más grandes. El modelado de este complejo comportamiento requiere la aplicación de métodos capaces de resolver de forma rigurosa cada una de las escalas. El trabajo que se presenta en esta tesis propone un modelado multiescala que es capaz de dar respuestas ingenieriles a las solicitudes técnicas del tungsteno, y que a su vez está apoyado por la rigurosa física subyacente a extensas simulaciones atomísticas. En primer lugar, las propiedades estáticas y dinámicas de las dislocaciones de tipo tornillo en cinco potenciales interatómicos de tungsteno son comparadas, determinando cuáles de ellos garantizan una mayor fidelidad física y eficiencia computacional. Las grandes tasas de deformación asociadas a las técnicas de dinámica molecular hacen que las funciones de movilidad de las dislocaciones obtenidas no puedan ser utilizadas en los siguientes pasos del modelado multiescala. En este trabajo, proponemos dos métodos alternativos para obtener las funciones de movilidad de las dislocaciones: un modelo Monte Cario cinético y expresiones analíticas. El conjunto de parámetros necesarios para formular el modelo de Monte Cario cinético y la ley de movilidad analítica son calculados atomísticamente. Estos parámetros incluyen, pero no se limitan a: la determinación de las entalpias y energías de formación de las parejas de escalones que forman las dislocaciones, la parametrización de los efectos de no Schmid característicos en materiales bcc,etc. Conociendo la ley de movilidad de las dislocaciones en función del esfuerzo aplicado y la temperatura, se introduce esta relación como ecuación de flujo dentro de un modelo de plasticidad cristalina. La predicción del modelo sobre la dependencia del límite de fluencia con la temperatura es validada experimentalmente con ensayos uniaxiales en tungsteno monocristalino. A continuación, se calcula el límite de fluencia al aplicar ensayos uniaxiales de tensión para un conjunto de orientaciones cristalográticas dentro del triángulo estándar variando la tasa de deformación y la temperatura de los ensayos. Finalmente, y con el objetivo de ser capaces de predecir una respuesta más dúctil del tungsteno para una variedad de estados de carga, se realizan ensayos biaxiales de tensión sobre algunas de las orientaciones cristalográficas ya estudiadas en función de la temperatura.-------------------------------------------------------------------------ABSTRACT ----------------------------------------------------------Tungsten and tungsten alloys are being considered as leading candidates for structural and functional materials in future fusion energy devices. The most attractive properties of tungsten for the design of magnetic and inertial fusion energy reactors are its high melting point, high thermal conductivity, low sputtering yield and low longterm disposal radioactive footprint. However, tungsten also presents a very low fracture toughness, mostly associated with inter-granular failure and bulk plasticity, that limits its applications. As a result of these various and complex conditions of work, the study, development and design of these materials is one of the most important challenges that have emerged in recent years to the scientific community in the field of materials for energy applications. The plastic behavior of body-centered cubic (bcc) refractory metals like tungsten is governed by the kink-pair mediated thermally activated motion of h¿ (\1 11)i screw dislocations on the atomistic scale and by ensembles and interactions of dislocations at larger scales. Modeling this complex behavior requires the application of methods capable of resolving rigorously each relevant scale. The work presented in this thesis proposes a multiscale model approach that gives engineering-level responses to the technical specifications required for the use of tungsten in fusion energy reactors, and it is also supported by the rigorous underlying physics of extensive atomistic simulations. First, the static and dynamic properties of screw dislocations in five interatomic potentials for tungsten are compared, determining which of these ensure greater physical fidelity and computational efficiency. The large strain rates associated with molecular dynamics techniques make the dislocation mobility functions obtained not suitable to be used in the next steps of the multiscale model. Therefore, it is necessary to employ mobility laws obtained from a different method. In this work, we suggest two alternative methods to get the dislocation mobility functions: a kinetic Monte Carlo model and analytical expressions. The set of parameters needed to formulate the kinetic Monte Carlo model and the analytical mobility law are calculated atomistically. These parameters include, but are not limited to: enthalpy and energy barriers of kink-pairs as a function of the stress, width of the kink-pairs, non-Schmid effects ( both twinning-antitwinning asymmetry and non-glide stresses), etc. The function relating dislocation velocity with applied stress and temperature is used as the main source of constitutive information into a dislocation-based crystal plasticity framework. We validate the dependence of the yield strength with the temperature predicted by the model against existing experimental data of tensile tests in singlecrystal tungsten, with excellent agreement between the simulations and the measured data. We then extend the model to a number of crystallographic orientations uniformly distributed in the standard triangle and study the effects of temperature and strain rate. Finally, we perform biaxial tensile tests and provide the yield surface as a function of the temperature for some of the crystallographic orientations explored in the uniaxial tensile tests.

Las ordenanzas ambientales como instrumento para la eficiencia energertica de los barrios residenciales

(re) inventar a cidade em tempos de midança. 6º Congreso Luso-Brasileiro para o Planeamiento Urbano, Tegional. Integrado e Sustentável

Uva mecânica : A colheita da uva no Brasil ainda é feita de forma manual na maioria dos parreirais. Mas neste cenário já é possível ver máquinas realizando a tarefa antes feita pelas mãos dos trabalhadores, operação que em outros países já possui uma taxa de mecanização maior.

Uva mecânica A colheita da uva no Brasil ainda é feita de forma manual na maioria dos parreirais. Mas neste cenário já é possível ver máquinas realizando a tarefa antes feita pelas mãos dos trabalhadores, operação que em outros países já possui uma taxa de mecanização maior.

La rehabilitación de Ciudad de los Ángeles en Madrid. Mejorar la habitabilidad de un barrio de 1950 con criterios de eficiencia energética

Regeneración urbana, habitabilidad, rehabilitación, eficiencia energética, calidad ambiental

Análisis de la eficiencia de los sistemas de colonias de hormigas para problemas de accesibilidad

A pesar de los avances en materia de predicción, los desastres naturales siguen teniendo consecuencias devastadoras. Entre los principales problemas a los que se enfrentan los equipos de ayuda y rescate después de un desastre natural o provocado por el hombre se encuentra la planificación de las tareas de reparación de carreteras para conseguir la máxima ventaja de los limitados recursos económicos y humanos. En la presente Tesis Fin de Máster se intenta dar solución al problema de la accesibilidad, es decir, maximizar el número de supervivientes que consiguen alcanzar el centro regional más cercano en un tiempo mínimo mediante la planificación de qué carreteras rurales deberían ser reparadas dados unos recursos económicos y humanos limitados. Como se puede observar, es un problema combinatorio ya que el número de planes de reparación y conexiones entre las ciudades y los centros regionales crece de forma exponencial con el tamaño del problema. Para la resolución del problema se comienza analizando una adaptación básica de los sistemas de colonias de hormigas propuesta por otro autor y se proponen múltiples mejoras sobre la misma. Posteriormente, se propone una nueva adaptación más avanzada de los sistemas de colonias de hormiga al problema, el ACS con doble hormiga. Este sistema hace uso de dos tipos distintos de hormigas, la exploradora y la trabajadora, para resolver simultáneamente el problema de encontrar los caminos más rápidos desde cada ciudad a su centro regional más cercano (exploradora), y el de obtener el plan óptimo de reparación que maximice la accesibilidad de la red (trabajadora). El algoritmo propuesto se ilustra por medio de un ejemplo de gran tamaño que simula el desastre natural ocurrido en Haití, y su rendimiento es comparado con la combinación de dos metaheurísticas, GRASP y VNS.---ABSTRACT---In spite of the advances in forecasting, natural disaster continue to ocasionate devastating consequences. One of the main problems relief teams face after a natural or man-made disaster is how to plan rural road repair work to take maximum advantage of the limited available financial and human resources. In this Master´s Final Project we account for the accesability issue, that is, to maximize the number of survivors that reach the nearest regional center in a minimum time by planning whic rural roads should be repaired given the limited financial and human resources. This is a combinatorial problem since the number of possible repairing solutions and connections between cities and regional centers grows exponentially with the size of the problem. In order to solve the problem, we analyze the basic ant colony system adaptation proposed by another author and point out multiple improvements on it. Then, we propose a novel and more advance adaptation of the ant colony systems to the problem, the double- ant ACS. This system makes use of two diferent type of ants, the explorer and the worker, to simultaneously solve the problem of finding the shorthest paths from each city to their nearest regional center (explorer), and the problem of identifying the optimal repairing plan that maximize the network accesability (worker). The proposed algorithm is illustrated by means of a big size example that simulates the natural disaster occurred in Haiti, and its performance is compared with a combination of two metaheuristics, GRASP and VNS.

Caracterización físico-mecánica del yeso de construcción con adición de residuos de poliestireno extruído y residuos cerámicos

Este trabajo estudia el comportamiento de una matriz de yeso de construcción a la que se le han añadido residuos de construcción y demolición (RCD), residuos de poliestireno extruido (XPS) y residuos cerámicos respectivamente, combinados y en diferentes porcentajes en función del peso del yeso. Los residuos de XPS son producto de una obra en Madrid donde el material fue utilizado como aislamiento térmico y los residuos cerámicos corresponden a trozos de ladrillos toscos encontrados en una obra paralizada en la ciudad de Ávila. Se confeccionaron probetas con porcentajes hasta 3% de XPS y hasta 50% de cerámicos en función del peso del yeso utilizado, como referencia se confeccionaron probetas sin adición de RCD. Fueron ensayadas en laboratorio y se determinaron las características físicas y mecánicas de las mismas. Tras un análisis comparativo se evidencian que la adición de residuos de XPS y residuos cerámicos en conjunto disminuye la densidad seca del material y la absorción de agua por capilaridad, en algunos casos disminuye la conductividad térmica y aumenta la dureza superficial. ABSTRACT: This paper studies the behavior of a building gypsum matrix to which have been added Construction and Demolition waste (CDW), residues of extruded polystyrene (XPS) and ceramic waste respectively, and combined in different percentages depending on the weight of gypsum. XPS waste are the product of a work in Madrid where the material was used as thermal insulation and ceramic waste correspond to pieces of rough bricks found in a paralyzed work in the city of Avila. Specimens were prepared with percentages up to 3% of XPS and up to 50% depending on the weight ceramic gypsum used as reference samples were prepared without addition of CDW. They were tested in laboratory and determined the physical and mechanical characteristics thereof. After a comparative analysis they show that the addition of ceramic waste and waste XPS decreases together dry material density and water absorption by capillary action, in some cases the thermal conductivity decreases and increases surface hardness.

Investigaciones realizadas en la Universidad de Aveiro sobre caracterización mecánica de las construcciones existentes en adobe en Portugal y propuestas de rehabilitación y refuerzo. Resultados alcanzados

En Portugal, sobre todo en la región de Aveiro, las estructuras de fábrica de adobe fueron el sistema constructivo predominante durante la primera mitad del siglo XX. En la actualidad es aún significativo el legado del patrimonio – cerca de 30% de las construcciones existentes – el cual está formado por un conjunto significativo de construcciones de elevado valor histórico y arquitectónico. En las últimas décadas se ha descuidado su rehabilitación, que ha redundado en el estado actual de daño en muchas de estas construcciones, con desfavorables condiciones de seguridad, y en algunos casos comprometiendo su uso. No obstante, en los últimos años, se observa un creciente interés, asociado a factores ambientales, económicos y de salvaguarda del patrimonio, tanto por parte de las administraciones locales como en el ámbito de los inversores privados, por la conservación y rehabilitación de estas construcciones. Ante esta nueva realidad, que pone de manifiesto la necesidad de desarrollar soluciones y técnicas de intervención que permitan, a través de la rehabilitación y/o refuerzo, prolongar la vida útil de estas estructuras, en la Universidad de Aveiro un grupo de investigadores viene desarrollando un trabajo multidisciplinar sobre las construcciones de adobe en Portugal. El objetivo principal del mismo se centra en la creación de una base de resultados experimentales que apoye las intervenciones de rehabilitación y/o refuerzo del patrimonio edificado, y el desarrollo de soluciones de mejora de su comportamiento estructural. Esta investigación se constituye como un guía que procura contribuir al conocimiento de estas fábricas, enfocado sobre todo a la rehabilitación de las mismas, aportando información que permita a los diferentes profesionales que trabajan en este ámbito actuar en la conservación de este patrimonio

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.