Curvature and velocity of methane-air Bunsen flame tips.

PIV and photographic recording are used to measure the velocity of the fresh gas and the shape of the reaction layer in a region around the tip of a methane-air Bunsen flame attached to a cylindrical burner. The results compare well with numerical simulations carried out with an infinite activation energy reaction model. The experimental and numerical results confirm that the well-known linear relation between flame velocity and flame stretch derived from asymptotic theory for weakly curved and strained flames is valid for small and moderate values of the flame stretch if the modified definition of stretch introduced by Echekki and Mungal (Proc Combust Inst 23:455–461, 1990) and Poinsot et al. (Combust Sci Technol 81:45–73, 1992) is used. However, the relation between flame velocity and modified stretch ceases to be linear and approaches a square root law for large values of the stretch, when the curvature of the flame tip becomes large compared to the inverse of the thickness of a planar flame.

Interactive Graphic Simulation: An Advanced Methodology to Improve the Teaching-Learning Process in Nuclear Engineering Education and Training

Nowadays, computer simulators are becoming basic tools for education and training in many engineering fields. In the nuclear industry, the role of simulation for training of operators of nuclear power plants is also recognized of the utmost relevance. As an example, the International Atomic Energy Agency sponsors the development of nuclear reactor simulators for education, and arranges the supply of such simulation programs. Aware of this, in 2008 Gas Natural Fenosa, a Spanish gas and electric utility that owns and operate nuclear power plants and promotes university education in the nuclear technology field, provided the Department of Nuclear Engineering of Universidad Politécnica de Madrid with the Interactive Graphic Simulator (IGS) of “José Cabrera” (Zorita) nuclear power plant, an industrial facility whose commercial operation ceased definitively in April 2006. It is a state-of-the-art full-scope real-time simulator that was used for training and qualification of the operators of the plant control room, as well as to understand and analyses the plant dynamics, and to develop, qualify and validate its emergency operating procedures.

Recuperación secundaria en campos de petróleo y su conversión en almacenamientos subterráneos de gas natural

El objetivo de este proyecto es estudiar la recuperación secundaria de petróleo de la capa sureste Ayoluengo del campo Ayoluengo, Burgos (España), y su conversión en un almacenamiento subterráneo de gas. La capa Ayoluengo se ha considerado como una capa inclinada de 60 km por 10 km de superficie por 30 m de espesor en el que se han perforado 20 pozos, y en donde la recuperación primaria ha sido de un 19%. Se ha realizado el ajuste histórico de la recuperación primaria de gas, petróleo y agua de la capa desde el año 1965 al 2011. La conversión a almacenamiento subterráneo de gas se ha realizado mediante ciclos de inyección de gas, de marzo a octubre, y extracción de gas, de noviembre a febrero, de forma que se incrementa la presión del campo hasta alcanzar la presión inicial. El gas se ha inyectado y extraído por 5 pozos situados en la zona superior de la capa. Al mismo tiempo, se ha realizado una recuperación secundaria debido a la inyección de gas natural de 20 años de duración en donde la producción de petróleo se realiza por 14 pozos situados en la parte inferior de la capa. Para proceder a la simulación del ajuste histórico, conversión en almacenamiento y recuperación secundaria se utilizó el simulador Eclipse100. Los resultados obtenidos fueron una recuperación secundaria de petróleo de un 9% más comparada con la primaria. En cuanto al almacenamiento de gas natural, se alcanzó la presión inicial consiguiendo un gas útil de 300 Mm3 y un gas colchón de 217,3 Mm3. ABSTRACT The aim of this project is to study the secondary recovery of oil from the southeast Ayoluengo layer at the oil field Ayoluengo, Burgos (Spain), and its conversion into an underground gas storage. The Ayoluengo layer is an inclined layer of 60 km by 10km of area by 30 m gross and with 20 wells, which its primary recovery is of 19%. The history matching of the production of oil, gas and water has been carried out from the year 1965 until 2011. The conversion into an underground gas storage has been done in cycles of gas injection from March to October, and gas extraction from November to February, so that the reservoir pressure increases until it gets to the initial pressure. The gas has been injected and extracted through five well situated in the top part of the layer. At the same time, the secondary recovery has occurred due to de injection of natural gas during 20 years where the production of oil has been done through 14 wells situated in the lowest part of the layer. To proceed to the simulation of the history match, the conversion into an underground gas storage and its secondary recovery, the simulator used was Eclipse100. The results were a secondary recovery of oil of 9% more, compared to the primary recovery and concerning the underground gas storage, the initial reservoir pressure was achieved with a working gas of 300 Mm3 and a cushion gas of 217,3 Mm3.

Curvature and velocity of methane-air Bunsen flame tips

PIV and photographic recording are used to measure the velocity of the fresh gas and the shape of the reaction layer in a region around the tip of a methane- air Bunsen flame attached to a cylindrical burner. The results compare well with numerical simulations carried out with an infinite activation energy reaction model. The experimental and numerical results confirm that the well-known linear relation between flame velocity and flame stretch derived from asymptotic theory for weakly curved and strained flames is valid for small and moderate values of the flame stretch if the modified definition of stretch introduced by Echekki and Mungal (Proc Combust Inst 23:455–461,1990) and Poinsot et al. (Combust Sci Technol 81:45–73,1992) is used. However, the relation between flame velocity and modified stretch ceases to be linear and approaches a square root law for large values of the stretch, when the curvature of the flame tip becomes large compared to the inverse of the thickness of a planar flame

Evaluación ambiental de las tecnologías de recuperación de gas en formaciones no convencionales

Este proyecto pretende ofrecer una visión general de una de las tecnologías más actuales de recuperación de gas en formaciones no convencionales: fracturación hidráulica o “fracking”. El proyecto está motivado por la necesidad de responder a diferentes cuestiones sobre los efectos ambientales, sociales y en la salud humana derivados de la utilización de esa tecnología. Ofrece, además, una descripción del proceso y utilización de la tecnología haciendo especial mención de los riesgos inherentes de su uso, aunque también se intenta establecer una vía de aceptación para su desarrollo cuyo fin último, a parte de los beneficios económicos de quienes la usan, es el de posibilitar la transición hacia el uso de unos recursos (energías fósiles de extracción no convencional) que requieren de dichas técnicas para mantener, a lo largo del tiempo, el suministro de una energía que se supone más respetuosa con el medio ambiente: el gas natural. En primer lugar se expone, a modo introductorio, la necesidad de utilización de nuevas técnicas de estimulación de pozos y su utilización para satisfacer las necesidades energéticas mundiales en los próximos años. A continuación se hace una revisión del marco regulatorio aplicable al gas no convencional. Seguidamente, se hace una descripción de los recursos y fuentes no convencionales de gas y la descripción del proceso de fracturación hidráulica. Se analizan los incidentes relacionados con su desarrollo y las posibilidades y mecanismos que pueden adoptarse para reducirlos. Finalmente, se proponen vías alternativas basadas en las mejores técnicas aplicables al uso de la tecnología, cuya finalidad sea la mayor consideración ambiental posible y el menor riesgo posible en la salud de las personas. ABSTRACT This project aims to provide an overview of the latest technologies in gas recovery unconventional formations: hydraulic fracturing or "fracking". The project is motivated by the need to respond to various questions on the environmental, social and human health arising from the use of this technology. It also offers a description of the process and use of technology with special mention of the inherent risks of their use, but also tries to establish a path of acceptance for development whose ultimate goal, apart from the economic benefits of those who use is of enabling the transition to the use of certain resources (fossil energy extraction unconventional) which require such techniques to maintain, over time, of an energy supply which is more environmentally friendly: natural gas. First discussed the need to use new well stimulation techniques and their use to meet the world's energy needs in the coming years. Below is a review of the regulatory framework applicable to unconventional gas. Next, there is a description of resources and unconventional sources of gas, and the description of the process of hydraulic fracturing. We analyze the incidents related to its development and the possibilities and mechanisms that can be taken to reduce them. Finally, we suggest alternative routes based on the best techniques applicable to the use of technology, aiming at the highest possible environmental consideration and the least possible risk to the health of people.

Estudio experimental sobre la efectividad y la viabilidad de distintas soluciones constructivas para reducir la concentración de gas radón en edificaciones

El trabajo de investigación desarrollado que ha dado lugar a la realización de esta Tesis, aborda la protección de los edificios frente a la entrada de gas radón y su acumulación en los espacios habitados. Dicho gas (isótopo del radón Rn-222) es un elemento radiactivo que se genera, principalmente, en terrenos con altos contenidos de radio (terrenos graníticos por ejemplo). Su alto grado de movilidad permite que penetre en los edificios a través de los materiales de cerramiento del mismo (porosidad de los materiales, fisuras, grietas y juntas) y se acumule en su interior, donde puede ser inhalado en altas concentraciones. La Organización Mundial de la Salud, califica al radón como agente cancerígeno de grado 1. Según este Organismo, el radón es la segunda causa de contracción de cáncer pulmonar detrás del tabaco. Como respuesta a esta alarma, distintos estados ya han elaborado normativas en las que se proponen soluciones para que los niveles de concentración de radón no superen los valores recomendados por los organismos internacionales responsables de la protección radiológica. En España aún no existe normativa de protección frente a este agente cancerígeno causante de numerosas muertes, y es por tal motivo evidente la necesidad de aportar documentación técnica que ayude a las administraciones nacionales y locales a desarrollar dicha normativa para ajustarse a las recomendaciones europeas e internacionales sobre los niveles que no se deben superar y que, por otro lado, ya han contemplado una gran cantidad de países. Como principal aportación de este trabajo se muestran los resultados de reducción de concentración de gas radón de distintas soluciones constructivas enfocadas a frenar la entrada de gas radón al interior de los edificios haciendo uso de técnicas y materiales habituales en el ámbito de la construcción en España. Para ello, se han estudiado las efectividades de dichas soluciones, en lo referente a su capacidad para frenar la inmisión de radón, en un prototipo de vivienda construido al efecto en un terreno con altas concentraciones de radón. Las soluciones propuestas y ensayadas han sido el resultado de una labor de optimización de los sistemas estudiados en la bibliografía con el fin de adaptar las técnicas a los sistemas constructivos habituales en España y en concreto a la situación real del prototipo de vivienda construido en un lugar con contenidos de radón en terreno muy elevados. El trabajo incluye un capítulo inicial con los conceptos básicos necesarios para entender la problemática que supone habitar en espacios con altos contenidos de radón. ABSTRACT The research developed, which has led to the completion of this thesis, deal with the protection of buildings against entry of radon gas and its accumulation in the ocupated spaces. This gas (radon isotope Rn-222) is a radioactive element generated, mainly, in areas with high levels of radio (granitic terrain for example). Its high mobility allows entering in buildings through the enclosure materials of it (porosity of materials, cracks, crevices and joints) and accumulates inside, where it can be inhaled in high concentrations. The World Health Organization describes radon gas as a carcinogen agent in level 1. According to this Agency, radon is the second leading cause of lung cancer behind tobacco. In response to this alarm, some states have developed regulations that propose solutions to reduce radon concentration levels for not exceeding the values recommended by international agencies responsible in radiation protection. In Spain there is still no legislation to protect against this carcinogen element that cause numerous deaths, and for that reason it is evident the need to provide technical documentation to help the national and local governments to develop legislation for reaching the European and international levels recommendations. As the main contribution of this work are the results of reducing radon concentration using different constructive solutions aimed to stop radon entry in buildings, with techniques and materials common in Spain. To do this, effectiveness of such solutions, have been studied in terms of its ability to stop radon entry in a housing prototype built for this purpose in an area with high radon levels. The solutions proposed and tested have been the result of a process of optimization of systems studied in the literature in order to adapt the techniques to Spanish building material and, specifically, to the actual situation of housing prototype built in a place with high contents of radon in soil. The work includes an initial chapter with the basic concepts needed to understand the problem of living in areas with high levels of radon.

In vitro ruminal fermentation of diets containing wheat straw and date pits as forage.

The pH, VFA concentration, total gas and met hane production were determined in the rumen of four Sicilo- Sarde rams fitted with permanent canulas. Rams received a ration that included 1.5 kg DM of oat hay and were supplemented with one of four concentrates: CC (10% barley, 43.3% corn, 25% wheat bran, 17.7% soybean meal, 4% sheep Vitamin and Mineral Mixture (VMM)), SC (66% white sorghum, 30% faba, 4% sheep VMM); TC (71% triticale, 18% faba, 7%, soybean meal, 4% VMM) or BC (71.5% barley, 17.5% faba, 7% soybean meal and 4% VMM). 50 ml samples were taken before, 2, 5 and 8 hours after the morning meal. Total gas was determined on rumen content before the morning meal. The rumen pH was statistically different (P<0.05) before and 2 hours after the morning meal among concentrates feed. It was in favour of TC and BC (P<0.05) concentrates but was comparable at the end of the day. The concentration of VFA was significantly higher (P<0.05) for diets TC and BC following the meal and became comparable among concentrates thereafter. The proportion of acetate and butyrate acids evolved in the same way during the day regardless of the regimen. The total volu me of gas was different (P<0.05) among diets, the BC showed the highest value (87.00±17.29 ml) while the lowest value was found in the TC concentrate (56.58±13.06 ml). The CH4 production for the BC was significantly different (P<0.05) from that of TC. Quantities produced by the CC and SC were similar (22.08±4.18vs . 21.16±3.21).

The influence of diet on the effectiveness of garlic oil and cinnamaldehyde to manipulate in vitro ruminal fermentation and methane production.

The objective of this study was to evaluate the effects of increasing doses [0 (control: CON), 20, 60, 180 and 540 mg/L incubation medium] of garlic oil (GO) and cinnamaldehyde (CIN) on in vitro ruminal fermentation of two diets. Batch cultures of mixed ruminal microorganisms were inoculated with ruminal fluid from four sheep fed a medium-concentrate diet (MC; 50 : 50 alfalfa hay : concentrate) or four sheep fed a high-concentrate diet (HC; 15 : 85 barley straw : concentrate). Diets MC and HC were representative of those fed to dairy and fattening ruminants, respectively. Samples of each diet were used as incubation substrates for the corresponding inoculum, and the incubation was repeated on 4 different days (four replicates per experimental treatment). There were GO × diet-type and CIN × diet-type interactions (P < 0.001–0.05) for many of the parameters determined, indicating different effects of both oils depending on the diet type. In general, effects of GO were more pronounced for MC compared with HC diet. Supplementation of GO did not affect (P > 0.05) total volatile fatty acid (VFA) production at any dose. For MC diet, GO at 60, 180 and 540 mg/L decreased (P < 0.05) molar proportion of acetate (608, 569 and 547 mmol/mol total VFA, respectively), and increased (P < 0.05) propionate proportion (233, 256 and 268 mmol/mol total VFA, respectively), compared with CON values (629 and 215 mmol/mol total VFA for acetate and propionate, respectively). A minimum dose of 180 mg of GO/L was required to produce similar modifications in acetate and propionate proportions with HC diet, but no effects (P > 0.05) on butyrate proportion were detected. Methane/VFA ratio was reduced (P < 0.05) by GO at 60, 180 and 540 mg/L for MC diet (0.23, 0.16 and 0.10 mol/mol, respectively), and by GO at 20, 60, 180 and 540 mg/L for HC diet (0.19, 0.19, 0.16 and 0.08 mol/mol, respectively), compared with CON (0.26 and 0.21 mol/mol for MC and HC diets, respectively). No effects (P = 0.16–0.85) of GO on final pH and concentrations of NH3-N and lactate were detected. For both diet types, the highest CIN dose decreased (P < 0.05) production of total VFA, gas and methane, which would indicate an inhibition of fermentation. Compared with CON, CIN at 180 mg/L increased (P < 0.05) acetate proportion for the MC (629 and 644 mmol/mol total VFA for CON and CIN, respectively) and HC (525 and 540 mmol/mol total VFA, respectively) diets, without affecting the proportions of any other VFA or total VFA production. Whereas for MC diet CIN at 60 and 180 mg/L decreased (P < 0.05) NH3-N concentrations compared with CON, only a trend (P < 0.10) was observed for CIN at 180 mg/L with the HC diet. Supplementation of CIN up to 180 mg/L did not affect (P = 0.18–0.99) lactate concentrations and production of gas and methane for any diet. The results show that effectiveness of GO and CIN to modify ruminal fermentation may depend on diet type, which would have practical implications if they are confirmed in vivo.

Multi-fluid Eulerian model of an electrospray in a host gas.

An Eulerian multifluid model is used to describe the evolution of an electrospray plume and the flow induced in the surrounding gas by the drag of the electrically charged spray droplets in the space between an injection electrode containing the electrospray source and a collector electrode. The spray is driven by the voltage applied between the two electrodes. Numerical computations and order-of-magnitude estimates for a quiescent gas show that the droplets begin to fly back toward the injection electrode at a certain critical value of the flux of droplets in the spray, which depends very much on the electrical conditions at the injection electrode. As the flux is increased toward its critical value, the electric field induced by the charge of the droplets partially balances the field due to the applied voltage in the vicinity of the injection electrode, leading to a spray that rapidly broadens at a distance from its origin of the order of the stopping distance at which the droplets lose their initial momentum and the effect of their inertia becomes negligible. The axial component of the electric field first changes sign in this region, causing the fly back. The flow induced in the gas significantly changes this picture in the conditions of typical experiments. A gas plume is induced by the drag of the droplets whose entrainment makes the radius of the spray away from the injection electrode smaller than in a quiescent gas, and convects the droplets across the region of negative axial electric field that appears around the origin of the spray when the flux of droplets is increased. This suppresses fly back and allows much higher fluxes to be reached than are possible in a quiescent gas. The limit of large droplet-to-gas mass ratio is discussed. Migration of satellite droplets to the shroud of the spray is reproduced by the Eulerian model, but this process is also affected by the motion of the gas. The gas flow preferentially pushes satellite droplets from the shroud to the core of the spray when the effect of the inertia of the droplets becomes negligible, and thus opposes the well-established electrostatic/inertial mechanism of segregation and may end up concentrating satellite droplets in an intermediate radial region of the spray.

Modelling and forecasting fossil fuels, CO2 and electricity prices and their volatilities

In the current uncertain context that affects both the world economy and the energy sector, with the rapid increase in the prices of oil and gas and the very unstable political situation that affects some of the largest raw materials’ producers, there is a need for developing efficient and powerful quantitative tools that allow to model and forecast fossil fuel prices, CO2 emission allowances prices as well as electricity prices. This will improve decision making for all the agents involved in energy issues. Although there are papers focused on modelling fossil fuel prices, CO2 prices and electricity prices, the literature is scarce on attempts to consider all of them together. This paper focuses on both building a multivariate model for the aforementioned prices and comparing its results with those of univariate ones, in terms of prediction accuracy (univariate and multivariate models are compared for a large span of days, all in the first 4 months in 2011) as well as extracting common features in the volatilities of the prices of all these relevant magnitudes. The common features in volatility are extracted by means of a conditionally heteroskedastic dynamic factor model which allows to solve the curse of dimensionality problem that commonly arises when estimating multivariate GARCH models. Additionally, the common volatility factors obtained are useful for improving the forecasting intervals and have a nice economical interpretation. Besides, the results obtained and methodology proposed can be useful as a starting point for risk management or portfolio optimization under uncertainty in the current context of energy markets.

Continuum models of the mechanical behavior of rolled and die-cast magnesium alloys

En los últimos años ha habido una fuerte tendencia a disminuir las emisiones de CO2 y su negativo impacto medioambiental. En la industria del transporte, reducir el peso de los vehículos aparece como la mejor opción para alcanzar este objetivo. Las aleaciones de Mg constituyen un material con gran potencial para el ahorro de peso. Durante la última década se han realizado muchos esfuerzos encaminados a entender los mecanismos de deformación que gobiernan la plasticidad de estos materiales y así, las aleaciones de Mg de colada inyectadas a alta presión y forjadas son todavía objeto de intensas campañas de investigación. Es ahora necesario desarrollar modelos que contemplen la complejidad inherente de los procesos de deformación de éstos. Esta tesis doctoral constituye un intento de entender mejor la relación entre la microestructura y el comportamiento mecánico de aleaciones de Mg, y dará como resultado modelos de policristales capaces de predecir propiedades macro- y microscópicas. La deformación plástica de las aleaciones de Mg está gobernada por una combinación de mecanismos de deformación característicos de la estructura cristalina hexagonal, que incluye el deslizamiento cristalográfico en planos basales, prismáticos y piramidales, así como el maclado. Las aleaciones de Mg de forja presentan texturas fuertes y por tanto los mecanismos de deformación activos dependen de la orientación de la carga aplicada. En este trabajo se ha desarrollado un modelo de plasticidad cristalina por elementos finitos con el objetivo de entender el comportamiento macro- y micromecánico de la aleación de Mg laminada AZ31 (Mg-3wt.%Al-1wt.%Zn). Este modelo, que incorpora el maclado y tiene en cuenta el endurecimiento por deformación debido a las interacciones dislocación-dislocación, dislocación-macla y macla-macla, predice exitosamente las actividades de los distintos mecanismos de deformación y la evolución de la textura con la deformación. Además, se ha llevado a cabo un estudio que combina difracción de electrones retrodispersados en tres dimensiones y modelización para investigar el efecto de los límites de grano en la propagación del maclado en el mismo material. Ambos, experimentos y simulaciones, confirman que el ángulo de desorientación tiene una influencia decisiva en la propagación del maclado. Se ha observado que los efectos no-Schmid, esto es, eventos de deformación plástica que no cumplen la ley de Schmid con respecto a la carga aplicada, no tienen lugar en la vecindad de los límites de baja desorientación y se hacen más frecuentes a medida que la desorientación aumenta. Esta investigación también prueba que la morfología de las maclas está altamente influenciada por su factor de Schmid. Es conocido que los procesos de colada suelen dar lugar a la formación de microestructuras con una microporosidad elevada, lo cuál afecta negativamente a sus propiedades mecánicas. La aplicación de presión hidrostática después de la colada puede reducir la porosidad y mejorar las propiedades aunque es poco conocido su efecto en el tamaño y morfología de los poros. En este trabajo se ha utilizado un enfoque mixto experimentalcomputacional, basado en tomografía de rayos X, análisis de imagen y análisis por elementos finitos, para la determinación de la distribución tridimensional (3D) de la porosidad y de la evolución de ésta con la presión hidrostática en la aleación de Mg AZ91 (Mg- 9wt.%Al-1wt.%Zn) colada por inyección a alta presión. La distribución real de los poros en 3D obtenida por tomografía se utilizó como input para las simulaciones por elementos finitos. Los resultados revelan que la aplicación de presión tiene una influencia significativa tanto en el cambio de volumen como en el cambio de forma de los poros que han sido cuantificados con precisión. Se ha observado que la reducción del tamaño de éstos está íntimamente ligada con su volumen inicial. En conclusión, el modelo de plasticidad cristalina propuesto en este trabajo describe con éxito los mecanismos intrínsecos de la deformación de las aleaciones de Mg a escalas meso- y microscópica. Más especificamente, es capaz de capturar las activadades del deslizamiento cristalográfico y maclado, sus interacciones, así como los efectos en la porosidad derivados de los procesos de colada. ---ABSTRACT--- The last few years have seen a growing effort to reduce CO2 emissions and their negative environmental impact. In the transport industry more specifically, vehicle weight reduction appears as the most straightforward option to achieve this objective. To this end, Mg alloys constitute a significant weight saving material alternative. Many efforts have been devoted over the last decade to understand the main mechanisms governing the plasticity of these materials and, despite being already widely used, high pressure die-casting and wrought Mg alloys are still the subject of intense research campaigns. Developing models that can contemplate the complexity inherent to the deformation of Mg alloys is now timely. This PhD thesis constitutes an attempt to better understand the relationship between the microstructure and the mechanical behavior of Mg alloys, as it will result in the design of polycrystalline models that successfully predict macro- and microscopic properties. Plastic deformation of Mg alloys is driven by a combination of deformation mechanisms specific to their hexagonal crystal structure, namely, basal, prismatic and pyramidal dislocation slip as well as twinning. Wrought Mg alloys present strong textures and thus specific deformation mechanisms are preferentially activated depending on the orientation of the applied load. In this work a crystal plasticity finite element model has been developed in order to understand the macro- and micromechanical behavior of a rolled Mg AZ31 alloy (Mg-3wt.%Al-1wt.%Zn). The model includes twinning and accounts for slip-slip, slip-twin and twin-twin hardening interactions. Upon calibration and validation against experiments, the model successfully predicts the activity of the various deformation mechanisms and the evolution of the texture at different deformation stages. Furthermore, a combined three-dimensional electron backscatter diffraction and modeling approach has been adopted to investigate the effect of grain boundaries on twin propagation in the same material. Both experiments and simulations confirm that the misorientation angle has a critical influence on twin propagation. Non-Schmid effects, i.e. plastic deformation events that do not comply with the Schmid law with respect to the applied stress, are absent in the vicinity of low misorientation boundaries and become more abundant as misorientation angle increases. This research also proves that twin morphology is highly influenced by the Schmid factor. Finally, casting processes usually lead to the formation of significant amounts of gas and shrinkage microporosity, which adversely affect the mechanical properties. The application of hydrostatic pressure after casting can reduce the porosity and improve the properties but little is known about the effects on the casting’s pores size and morphology. In this work, an experimental-computational approach based on X-ray computed tomography, image analysis and finite element analysis is utilized for the determination of the 3D porosity distribution and its evolution with hydrostatic pressure in a high pressure diecast Mg AZ91 alloy (Mg-9wt.%Al-1wt.%Zn). The real 3D pore distribution obtained by tomography is used as input for the finite element simulations using an isotropic hardening law. The model is calibrated and validated against experimental stress-strain curves. The results reveal that the pressure treatment has a significant influence both on the volume and shape changes of individuals pores, which have been precisely quantified, and which are found to be related to the initial pore volume. In conclusion, the crystal plasticity model proposed in this work successfully describes the intrinsic deformation mechanisms of Mg alloys both at the mesoscale and the microscale. More specifically, it can capture slip and twin activities, their interactions, as well as the potential porosity effects arising from casting processes.

Development and implementation of an air quality integrated assessment modelling system for the Iberian Peninsula

La mejora de la calidad del aire es una tarea eminentemente interdisciplinaria. Dada la gran variedad de ciencias y partes involucradas, dicha mejora requiere de herramientas de evaluación simples y completamente integradas. La modelización para la evaluación integrada (integrated assessment modeling) ha demostrado ser una solución adecuada para la descripción de los sistemas de contaminación atmosférica puesto que considera cada una de las etapas involucradas: emisiones, química y dispersión atmosférica, impactos ambientales asociados y potencial de disminución. Varios modelos de evaluación integrada ya están disponibles a escala continental, cubriendo cada una de las etapas antesmencionadas, siendo el modelo GAINS (Greenhouse Gas and Air Pollution Interactions and Synergies) el más reconocido y usado en el contexto europeo de toma de decisiones medioambientales. Sin embargo, el manejo de la calidad del aire a escala nacional/regional dentro del marco de la evaluación integrada es deseable. Esto sin embargo, no se lleva a cabo de manera satisfactoria con modelos a escala europea debido a la falta de resolución espacial o de detalle en los datos auxiliares, principalmente los inventarios de emisión y los patrones meteorológicos, entre otros. El objetivo de esta tesis es presentar los desarrollos en el diseño y aplicación de un modelo de evaluación integrada especialmente concebido para España y Portugal. El modelo AERIS (Atmospheric Evaluation and Research Integrated system for Spain) es capaz de cuantificar perfiles de concentración para varios contaminantes (NO2, SO2, PM10, PM2,5, NH3 y O3), el depósito atmosférico de especies de azufre y nitrógeno así como sus impactos en cultivos, vegetación, ecosistemas y salud como respuesta a cambios porcentuales en las emisiones de sectores relevantes. La versión actual de AERIS considera 20 sectores de emisión, ya sea equivalentes a sectores individuales SNAP o macrosectores, cuya contribución a los niveles de calidad del aire, depósito e impactos han sido modelados a través de matrices fuentereceptor (SRMs). Estas matrices son constantes de proporcionalidad que relacionan cambios en emisiones con diferentes indicadores de calidad del aire y han sido obtenidas a través de parametrizaciones estadísticas de un modelo de calidad del aire (AQM). Para el caso concreto de AERIS, su modelo de calidad del aire “de origen” consistió en el modelo WRF para la meteorología y en el modelo CMAQ para los procesos químico-atmosféricos. La cuantificación del depósito atmosférico, de los impactos en ecosistemas, cultivos, vegetación y salud humana se ha realizado siguiendo las metodologías estándar establecidas bajo los marcos internacionales de negociación, tales como CLRTAP. La estructura de programación está basada en MATLAB®, permitiendo gran compatibilidad con software típico de escritorio comoMicrosoft Excel® o ArcGIS®. En relación con los niveles de calidad del aire, AERIS es capaz de proveer datos de media anual y media mensual, así como el 19o valor horario más alto paraNO2, el 25o valor horario y el 4o valor diario más altos para SO2, el 36o valor diario más alto para PM10, el 26o valor octohorario más alto, SOMO35 y AOT40 para O3. En relación al depósito atmosférico, el depósito acumulado anual por unidad de area de especies de nitrógeno oxidado y reducido al igual que de azufre pueden ser determinados. Cuando los valores anteriormente mencionados se relacionan con características del dominio modelado tales como uso de suelo, cubiertas vegetales y forestales, censos poblacionales o estudios epidemiológicos, un gran número de impactos puede ser calculado. Centrándose en los impactos a ecosistemas y suelos, AERIS es capaz de estimar las superaciones de cargas críticas y las superaciones medias acumuladas para especies de nitrógeno y azufre. Los daños a bosques se calculan como una superación de los niveles críticos de NO2 y SO2 establecidos. Además, AERIS es capaz de cuantificar daños causados por O3 y SO2 en vid, maíz, patata, arroz, girasol, tabaco, tomate, sandía y trigo. Los impactos en salud humana han sido modelados como consecuencia de la exposición a PM2,5 y O3 y cuantificados como pérdidas en la esperanza de vida estadística e indicadores de mortalidad prematura. La exactitud del modelo de evaluación integrada ha sido contrastada estadísticamente con los resultados obtenidos por el modelo de calidad del aire convencional, exhibiendo en la mayoría de los casos un buen nivel de correspondencia. Debido a que la cuantificación de los impactos no es llevada a cabo directamente por el modelo de calidad del aire, un análisis de credibilidad ha sido realizado mediante la comparación de los resultados de AERIS con los de GAINS para un escenario de emisiones determinado. El análisis reveló un buen nivel de correspondencia en las medias y en las distribuciones probabilísticas de los conjuntos de datos. Las pruebas de verificación que fueron aplicadas a AERIS sugieren que los resultados son suficientemente consistentes para ser considerados como razonables y realistas. En conclusión, la principal motivación para la creación del modelo fue el producir una herramienta confiable y a la vez simple para el soporte de las partes involucradas en la toma de decisiones, de cara a analizar diferentes escenarios “y si” con un bajo coste computacional. La interacción con políticos y otros actores dictó encontrar un compromiso entre la complejidad del modeladomedioambiental con el carácter conciso de las políticas, siendo esto algo que AERIS refleja en sus estructuras conceptual y computacional. Finalmente, cabe decir que AERIS ha sido creado para su uso exclusivo dentro de un marco de evaluación y de ninguna manera debe ser considerado como un sustituto de los modelos de calidad del aire ordinarios. ABSTRACT Improving air quality is an eminently inter-disciplinary task. The wide variety of sciences and stakeholders that are involved call for having simple yet fully-integrated and reliable evaluation tools available. Integrated AssessmentModeling has proved to be a suitable solution for the description of air pollution systems due to the fact that it considers each of the involved stages: emissions, atmospheric chemistry, dispersion, environmental impacts and abatement potentials. Some integrated assessment models are available at European scale that cover each of the before mentioned stages, being the Greenhouse Gas and Air Pollution Interactions and Synergies (GAINS) model the most recognized and widely-used within a European policy-making context. However, addressing air quality at the national/regional scale under an integrated assessment framework is desirable. To do so, European-scale models do not provide enough spatial resolution or detail in their ancillary data sources, mainly emission inventories and local meteorology patterns as well as associated results. The objective of this dissertation is to present the developments in the design and application of an Integrated Assessment Model especially conceived for Spain and Portugal. The Atmospheric Evaluation and Research Integrated system for Spain (AERIS) is able to quantify concentration profiles for several pollutants (NO2, SO2, PM10, PM2.5, NH3 and O3), the atmospheric deposition of sulfur and nitrogen species and their related impacts on crops, vegetation, ecosystems and health as a response to percentual changes in the emissions of relevant sectors. The current version of AERIS considers 20 emission sectors, either corresponding to individual SNAP sectors or macrosectors, whose contribution to air quality levels, deposition and impacts have been modeled through the use of source-receptor matrices (SRMs). Thesematrices are proportionality constants that relate emission changes with different air quality indicators and have been derived through statistical parameterizations of an air qualitymodeling system (AQM). For the concrete case of AERIS, its parent AQM relied on the WRF model for meteorology and on the CMAQ model for atmospheric chemical processes. The quantification of atmospheric deposition, impacts on ecosystems, crops, vegetation and human health has been carried out following the standard methodologies established under international negotiation frameworks such as CLRTAP. The programming structure isMATLAB ® -based, allowing great compatibility with typical software such as Microsoft Excel ® or ArcGIS ® Regarding air quality levels, AERIS is able to provide mean annual andmean monthly concentration values, as well as the indicators established in Directive 2008/50/EC, namely the 19th highest hourly value for NO2, the 25th highest daily value and the 4th highest hourly value for SO2, the 36th highest daily value of PM10, the 26th highest maximum 8-hour daily value, SOMO35 and AOT40 for O3. Regarding atmospheric deposition, the annual accumulated deposition per unit of area of species of oxidized and reduced nitrogen as well as sulfur can be estimated. When relating the before mentioned values with specific characteristics of the modeling domain such as land use, forest and crops covers, population counts and epidemiological studies, a wide array of impacts can be calculated. When focusing on impacts on ecosystems and soils, AERIS is able to estimate critical load exceedances and accumulated average exceedances for nitrogen and sulfur species. Damage on forests is estimated as an exceedance of established critical levels of NO2 and SO2. Additionally, AERIS is able to quantify damage caused by O3 and SO2 on grapes, maize, potato, rice, sunflower, tobacco, tomato, watermelon and wheat. Impacts on human health aremodeled as a consequence of exposure to PM2.5 and O3 and quantified as losses in statistical life expectancy and premature mortality indicators. The accuracy of the IAM has been tested by statistically contrasting the obtained results with those yielded by the conventional AQM, exhibiting in most cases a good agreement level. Due to the fact that impacts cannot be directly produced by the AQM, a credibility analysis was carried out for the outputs of AERIS for a given emission scenario by comparing them through probability tests against the performance of GAINS for the same scenario. This analysis revealed a good correspondence in the mean behavior and the probabilistic distributions of the datasets. The verification tests that were applied to AERIS suggest that results are consistent enough to be credited as reasonable and realistic. In conclusion, the main reason thatmotivated the creation of this model was to produce a reliable yet simple screening tool that would provide decision and policy making support for different “what-if” scenarios at a low computing cost. The interaction with politicians and other stakeholders dictated that reconciling the complexity of modeling with the conciseness of policies should be reflected by AERIS in both, its conceptual and computational structures. It should be noted however, that AERIS has been created under a policy-driven framework and by no means should be considered as a substitute of the ordinary AQM.

Co-simulation Methodologies for Hybrid and Electric Vehicle Dynamics

In recent decades, full electric and hybrid electric vehicles have emerged as an alternative to conventional cars due to a range of factors, including environmental and economic aspects. These vehicles are the result of considerable efforts to seek ways of reducing the use of fossil fuel for vehicle propulsion. Sophisticated technologies such as hybrid and electric powertrains require careful study and optimization. Mathematical models play a key role at this point. Currently, many advanced mathematical analysis tools, as well as computer applications have been built for vehicle simulation purposes. Given the great interest of hybrid and electric powertrains, along with the increasing importance of reliable computer-based models, the author decided to integrate both aspects in the research purpose of this work. Furthermore, this is one of the first final degree projects held at the ETSII (Higher Technical School of Industrial Engineers) that covers the study of hybrid and electric propulsion systems. The present project is based on MBS3D 2.0, a specialized software for the dynamic simulation of multibody systems developed at the UPM Institute of Automobile Research (INSIA). Automobiles are a clear example of complex multibody systems, which are present in nearly every field of engineering. The work presented here benefits from the availability of MBS3D software. This program has proven to be a very efficient tool, with a highly developed underlying mathematical formulation. On this basis, the focus of this project is the extension of MBS3D features in order to be able to perform dynamic simulations of hybrid and electric vehicle models. This requires the joint simulation of the mechanical model of the vehicle, together with the model of the hybrid or electric powertrain. These sub-models belong to completely different physical domains. In fact the powertrain consists of energy storage systems, electrical machines and power electronics, connected to purely mechanical components (wheels, suspension, transmission, clutch…). The challenge today is to create a global vehicle model that is valid for computer simulation. Therefore, the main goal of this project is to apply co-simulation methodologies to a comprehensive model of an electric vehicle, where sub-models from different areas of engineering are coupled. The created electric vehicle (EV) model consists of a separately excited DC electric motor, a Li-ion battery pack, a DC/DC chopper converter and a multibody vehicle model. Co-simulation techniques allow car designers to simulate complex vehicle architectures and behaviors, which are usually difficult to implement in a real environment due to safety and/or economic reasons. In addition, multi-domain computational models help to detect the effects of different driving patterns and parameters and improve the models in a fast and effective way. Automotive designers can greatly benefit from a multidisciplinary approach of new hybrid and electric vehicles. In this case, the global electric vehicle model includes an electrical subsystem and a mechanical subsystem. The electrical subsystem consists of three basic components: electric motor, battery pack and power converter. A modular representation is used for building the dynamic model of the vehicle drivetrain. This means that every component of the drivetrain (submodule) is modeled separately and has its own general dynamic model, with clearly defined inputs and outputs. Then, all the particular submodules are assembled according to the drivetrain configuration and, in this way, the power flow across the components is completely determined. Dynamic models of electrical components are often based on equivalent circuits, where Kirchhoff’s voltage and current laws are applied to draw the algebraic and differential equations. Here, Randles circuit is used for dynamic modeling of the battery and the electric motor is modeled through the analysis of the equivalent circuit of a separately excited DC motor, where the power converter is included. The mechanical subsystem is defined by MBS3D equations. These equations consider the position, velocity and acceleration of all the bodies comprising the vehicle multibody system. MBS3D 2.0 is entirely written in MATLAB and the structure of the program has been thoroughly studied and understood by the author. MBS3D software is adapted according to the requirements of the applied co-simulation method. Some of the core functions are modified, such as integrator and graphics, and several auxiliary functions are added in order to compute the mathematical model of the electrical components. By coupling and co-simulating both subsystems, it is possible to evaluate the dynamic interaction among all the components of the drivetrain. ‘Tight-coupling’ method is used to cosimulate the sub-models. This approach integrates all subsystems simultaneously and the results of the integration are exchanged by function-call. This means that the integration is done jointly for the mechanical and the electrical subsystem, under a single integrator and then, the speed of integration is determined by the slower subsystem. Simulations are then used to show the performance of the developed EV model. However, this project focuses more on the validation of the computational and mathematical tool for electric and hybrid vehicle simulation. For this purpose, a detailed study and comparison of different integrators within the MATLAB environment is done. Consequently, the main efforts are directed towards the implementation of co-simulation techniques in MBS3D software. In this regard, it is not intended to create an extremely precise EV model in terms of real vehicle performance, although an acceptable level of accuracy is achieved. The gap between the EV model and the real system is filled, in a way, by introducing the gas and brake pedals input, which reflects the actual driver behavior. This input is included directly in the differential equations of the model, and determines the amount of current provided to the electric motor. For a separately excited DC motor, the rotor current is proportional to the traction torque delivered to the car wheels. Therefore, as it occurs in the case of real vehicle models, the propulsion torque in the mathematical model is controlled through acceleration and brake pedal commands. The designed transmission system also includes a reduction gear that adapts the torque coming for the motor drive and transfers it. The main contribution of this project is, therefore, the implementation of a new calculation path for the wheel torques, based on performance characteristics and outputs of the electric powertrain model. Originally, the wheel traction and braking torques were input to MBS3D through a vector directly computed by the user in a MATLAB script. Now, they are calculated as a function of the motor current which, in turn, depends on the current provided by the battery pack across the DC/DC chopper converter. The motor and battery currents and voltages are the solutions of the electrical ODE (Ordinary Differential Equation) system coupled to the multibody system. Simultaneously, the outputs of MBS3D model are the position, velocity and acceleration of the vehicle at all times. The motor shaft speed is computed from the output vehicle speed considering the wheel radius, the gear reduction ratio and the transmission efficiency. This motor shaft speed, somehow available from MBS3D model, is then introduced in the differential equations corresponding to the electrical subsystem. In this way, MBS3D and the electrical powertrain model are interconnected and both subsystems exchange values resulting as expected with tight-coupling approach.When programming mathematical models of complex systems, code optimization is a key step in the process. A way to improve the overall performance of the integration, making use of C/C++ as an alternative programming language, is described and implemented. Although this entails a higher computational burden, it leads to important advantages regarding cosimulation speed and stability. In order to do this, it is necessary to integrate MATLAB with another integrated development environment (IDE), where C/C++ code can be generated and executed. In this project, C/C++ files are programmed in Microsoft Visual Studio and the interface between both IDEs is created by building C/C++ MEX file functions. These programs contain functions or subroutines that can be dynamically linked and executed from MATLAB. This process achieves reductions in simulation time up to two orders of magnitude. The tests performed with different integrators, also reveal the stiff character of the differential equations corresponding to the electrical subsystem, and allow the improvement of the cosimulation process. When varying the parameters of the integration and/or the initial conditions of the problem, the solutions of the system of equations show better dynamic response and stability, depending on the integrator used. Several integrators, with variable and non-variable step-size, and for stiff and non-stiff problems are applied to the coupled ODE system. Then, the results are analyzed, compared and discussed. From all the above, the project can be divided into four main parts: 1. Creation of the equation-based electric vehicle model; 2. Programming, simulation and adjustment of the electric vehicle model; 3. Application of co-simulation methodologies to MBS3D and the electric powertrain subsystem; and 4. Code optimization and study of different integrators. Additionally, in order to deeply understand the context of the project, the first chapters include an introduction to basic vehicle dynamics, current classification of hybrid and electric vehicles and an explanation of the involved technologies such as brake energy regeneration, electric and non-electric propulsion systems for EVs and HEVs (hybrid electric vehicles) and their control strategies. Later, the problem of dynamic modeling of hybrid and electric vehicles is discussed. The integrated development environment and the simulation tool are also briefly described. The core chapters include an explanation of the major co-simulation methodologies and how they have been programmed and applied to the electric powertrain model together with the multibody system dynamic model. Finally, the last chapters summarize the main results and conclusions of the project and propose further research topics. In conclusion, co-simulation methodologies are applicable within the integrated development environments MATLAB and Visual Studio, and the simulation tool MBS3D 2.0, where equation-based models of multidisciplinary subsystems, consisting of mechanical and electrical components, are coupled and integrated in a very efficient way.

Air-steam gasification of sewage sludge in a bubbling bed reactor: Effect of alumina as a primary catalyst

Numerous references can be found in scientific literature regarding biomass gasification. However, there are few works related to sludge gasification. A study of sewage sludge gasification process in a bubbling fluidised bed gasifier on a laboratory scale is here reported. The aim was to find the optimum conditions for reducing the production of tars and gain more information on the influx of different operating variables in the products resulting from the gasification of this waste. The variables studied were the equivalence ratio (ER), the steam-biomass ratio (SB) and temperature. Specifically, the ER was varied from 0.2 to 0.4, the SB from 0 to 1 and the temperature from 750 °C (1023 K) to 850 °C (1123 K). Although it was observed that tar production could be considerably reduced (up to 72%) by optimising the gasification conditions, the effect of using alumina (aluminium oxide, of proven efficacy in destroying the tar produced in biomass gasification) as primary catalyst in air and air-steam mixture tests was also verified. The results show that by adding small quantities of alumina to the bed (10% by weight of fed sludge) considerable reductions in tar production can be obtained (up to 42%) improving, at the same time, the lower heating value (LHV) of the gas and carbon conversion.

Use of alternative fuels obtained from renewable sources in brayton cycles

Analysis and simulation of the behaviour of gas turbines for power generation using different nonconventional fuels obtained from different renewable sources are presented. Three biomass-tobiofuel processes are considered: anaerobic digestion of biomass (biogas), biomass gasification (synthesis gas) and alcoholic fermentation of biomass and dehydration (bioethanol), each of them with two different biomass substrates (energy crops and municipal solid waste) as input. The gas turbine behaviour in a Brayton cycle is simulated both in an isolated operation and in combined cycle. The differences in gas turbine performance when fired with the considered biofuels compared to natural gas are studied from different points of view related with the current complex energetic context: energetic and exergetic efficiency of the simple/combined cycle and CO2 emissions. Two different tools have been used for the simulations, each one with a different approach: while PATITUG (own software) analyses the behaviour of a generic gas turbine allowing a total variability of parameters, GT-PRO (commercial software) is more rigid, albeit more precise in the prediction of real gas turbine behaviour. Different potentially interesting configurations and its thermodynamic parameters have been simulated in order to obtain the optimal range for all of them and its variation for each fuel.