Influencia de las restricciones medioambientales en el sector eléctrico

El objetivo de esta tesis es la caracterización de la generación térmica representativa de la existente en la realidad, para posteriormente proceder a su modelización y simulación integrándolas en una red eléctrica tipo y llevar a cabo estudios de optimización multiobjetivo económico medioambiental. Para ello, en primera instancia se analiza el contexto energético y eléctrico actual, y más concretamente el peninsular, en el que habiendo desaparecido las centrales de fuelóleo, sólo quedan ciclos combinados y centrales de carbón de distinto rango. Seguidamente se lleva a cabo un análisis de los principales impactos medioambientales de las centrales eléctricas basadas en combustión, representados sobre todo por sus emisiones de CO2, SO2 y NOx, de las medidas de control y mitigación de las mismas y de la normativa que les aplica. A continuación, a partir de las características de los combustibles y de la información de los consumos específicos, se caracterizan los grupos térmicos frente a las funciones relevantes que definen su comportamiento energético, económico y medioambiental, en términos de funciones de salida horarias dependiendo de la carga. Se tiene en cuenta la posibilidad de desnitrificación y desulfuración. Dado que las funciones objetivo son múltiples, y que están en conflicto unas con otras, se ha optado por usar métodos multiobjetivo que son capaces de identificar el contorno de puntos óptimos o frente de Pareto, en los que tomando una solución no existe otra que lo mejore en alguna de las funciones objetivo sin empeorarlo en otra. Se analizaron varios métodos de optimización multiobjetivo y se seleccionó el de las ε constraint, capaz de encontrar frentes no convexos y cuya optimalidad estricta se puede comprobar. Se integró una representación equilibrada de centrales de antracita, hulla nacional e importada, lignito y ciclos combinados en la red tipo IEEE-57, en la que se puede trabajar con siete centrales sin distorsionar demasiado las potencias nominales reales de los grupos, y se programó en Matlab la resolución de flujos óptimos de carga en alterna con el método multiobjetivo integrado. Se identifican los frentes de Pareto de las combinaciones de coste y cada uno de los tres tipos de emisión, y también el de los cuatro objetivos juntos, obteniendo los resultados de costes óptimos del sistema para todo el rango de emisiones. Se valora cuánto le cuesta al sistema reducir una tonelada adicional de cualquier tipo de emisión a base de desplazarse a combinaciones de generación más limpias. Los puntos encontrados aseguran que bajo unas determinadas emisiones no pueden ser mejorados económicamente, o que atendiendo a ese coste no se puede reducir más allá el sistema en lo relativo a emisiones. También se indica cómo usar los frentes de Pareto para trazar estrategias óptimas de producción ante cambios horarios de carga. ABSTRACT The aim of this thesis is the characterization of electrical generation based on combustion processes representative of the actual power plants, for the latter modelling and simulation of an electrical grid and the development of economic- environmental multiobjective optimization studies. In this line, the first step taken is the analysis of the current energetic and electrical framework, focused on the peninsular one, where the fuel power plants have been shut down, and the only ones remaining are coal units of different types and combined cycle. Then it is carried out an analysis of the main environmental impacts of the thermal power plants, represented basically by the emissions of CO2, SO2 y NOx, their control and reduction measures and the applicable regulations. Next, based on the combustibles properties and the information about the units heat rates, the different power plants are characterized in relation to the outstanding functions that define their energy, economic and environmental behaviour, in terms of hourly output functions depending on their load. Optional denitrification and desulfurization is considered. Given that there are multiple objectives, and that they go in conflictive directions, it has been decided the use of multiobjective techniques, that have the ability of identifying the optimal points set, which is called the Pareto front, where taken a solution there will be no other point that can beat the former in an objective without worsening it in another objective. Several multiobjective optimization methods were analysed and pondered, selecting the ε constraint technique, which is able to find no convex fronts and it is opened to be tested to prove the strict Pareto optimality of the obtained solutions. A balanced representation of the thermal power plants, formed by anthracite, lignite, bituminous national and imported coals and combined cycle, was integrated in the IEEE-57 network case. This system was selected because it deals with a total power that will admit seven units without distorting significantly the actual size of the power plants. Next, an AC optimal power flow with the multiobjective method implemented in the routines was programmed. The Pareto fronts of the combination of operative costs with each of the three emissions functions were found, and also the front of all of them together. The optimal production costs of the system for all the emissions range were obtained. It is also evaluated the cost of reducing an additional emission ton of any of the emissions when the optimal production mix is displaced towards cleaner points. The obtained solutions assure that under a determined level of emissions they cannot be improved economically or, in the other way, at a determined cost it cannot be found points of lesser emissions. The Pareto fronts are also applied for the search of optimal strategic paths to follow the hourly load changes.

Integración técnica y económica de energías renovables, en un sistema eléctrico liberalizado, a partir de vectores energéticos

El objetivo de la tesis es estudiar la bondad del almacenamiento de energía en hidrógeno para minorar los desvíos de energía respecto a su previsión de parques eólicos y huertas solares. Para ello se ha partido de datos de energías horarias previstas con 24 h de antelación y la energía real generada. Se ha procedido a dimensionar la planta de hidrógeno, a partir de una modelización de la operación de la misma, teniendo siempre como objetivo la limitación de los desvíos. Posteriormente, se ha procedido a simular la operación de la planta con dos objetivos en mente, uno limitar los desvíos y por otro lado operar la planta como una central de bombeo, generando hidrógeno en horas valle y generando electricidad en horas punta. Las dos simulaciones se han aplicado a tres parques eólicos de diferentes potencias, y a una huerta solar fotovoltaica. Se ha realizado un estudio económico para determinar la viabilidad de las plantas dimensionadas, obteniendo como resultado que no son viables a día de hoy y con la estimación de precios considerada, necesitando disminuir considerablemente los costes, dependiendo fuertemente de la bondad de los métodos de previsión de viento. Por último se ha estudiado la influencia de la disminución de los desvíos generados sobre una red tipo de 30 nudos, obteniendo como resultado, que si bien no disminuyen sensiblemente los extra costes generados en regulación, sí que mejora la penetración de las energías renovables no despachables en la red. Se observa disminuyen los vertidos eólicos cuando se usa la planta de hidrógeno. ABSTRACT The aim of this thesis is to study the benefit of hydrogen energy storage to minimize energy deviations of Wind Power and Solar Photovoltaic (PV) Power Plants compared to its forecast. To achieve this goal, first of all we have started with hourly energy data provided 24 h in advance (scheduled energy), and real generation (measured energy). Secondly, It has been sized the hydrogen plant, from a modeling of its working mode, always keeping the goal in mind of limiting energy imbalances. Subsequently, It have been simulated the plant working mode following two goals, one, to limit energy imbalances and secondly to operate the plant as a pumping power plant, generating hydrogen-in valley hours and generating electricity at peak hours. The two simulations have been applied to three wind power plants with different installed power capacities, and a photovoltaic solar power plant. It has been done an economic analysis in order to determine the viability of this sized plants, turning out not viable plants today with the estimated prices considered, requiring significantly lower costs, depending heavily on the reliability of the Wind Power forecast methods. Finally, It has been studied the influence of decreasing measured imbalances (of energy) in a 30 grid node, resulting that, while it not reduces significantly the extra costs generated by reserve power, it does improve the penetration of non-manageable renewable energy on the grid, by reducing the curtailments of power of these plants.

Software de diseño eléctrico para instalaciones fotovoltaicas

En este documento se describen los requisitos y/o especificaciones que servirán como base para una correcta definición de la aplicación a desarrollar. Esta aplicación se encargará de gestionar el diseño eléctrico de instalaciones fotovoltaicas. Este software será diseñado especialmente como herramienta de ayuda para los Ingenieros del Departamento de Construcción de Solaria Energía y Medio Ambiente, S.A.

Análisis de topología y restricciones de seguridad en mercados eléctricos competitivos: aplicación al sistema eléctrico de República Dominicana

Las técnicas de resolución de restricciones son mecanismos que garantizan que todas las transacciones en un mercado competitivo de electricidad se pueden lograr sin la violación de los límites de operación. Es un tema importante en un mercado competitivo de la electricidad, ya que un manejo inadecuado puede segmentar el sistema y hacer que algunos de los participantes puedan ejercer poder de mercado. El propósito principal de este artículo es presentar una metodología que permite garantizar un despacho seguro, tanto en estado normal como después de contingencias, utilizando la técnica del re-despacho tras el análisis topológico. El artículo describe paso a paso el proceso de aplicación de la metodología propuesta e ilustra su aplicación al caso del sistema eléctrico de la República Dominicana y cómo puede ayudar a solventar el problema mencionado cuando el ahorro durante los periodos de contingencia supera al extra-coste durante los periodos de funcionamiento normal

Contaminación cero

Devolver la calle al peatón, limitar el uso del coche e, incluso, eliminarlo de las urbes; aprovechar los avances tecnológicos para que sean más “inteligentes” y optimizar al máximo los recursos, cada vez más escasos, para que sean más sostenibles. Así sería la ciudad perfecta para los arquitectos

Estado del Arte del Uso del Gas de Gasificación Termoquímica de Biomasa(GG), en Motores de Combustión Interna Alternativos

El presente trabajo, desarrollado en el marco del Convenio de Cooperación educativa entre la ETSII - UPM y el Ciemat, se realiza con el fin de determinar líneas futuras de investigación y/o aplicación de la tecnología de gasificación termoquímica de biomasa integrada a motores de combustión interna alternativos (MCIA) para generación de potencia, motivados por la necesidad de reducir las emisiones contaminantes, aumentar el uso de las fuentes renovables de energía, reducir la dependencia económica de los combustibles fósiles, aprovechar energéticamente infinidad de residuos del sector agroindustrial y por la necesidad de generar energía a base de combustibles autóctonos que permitan resolver los problemas de suministro eléctrico en zonas no interconectadas eléctricamente en países en vía de desarrollo. En el capítulo 1 se comienza por presentar los objetivos y la justificación del presente trabajo, se enmarca esta tecnología desde el punto de vista histórico, de los combustibles y de los gasógenos, que fue la primera aplicación extendida masivamente durante las dos guerras mundiales en Europa. Además se hace un breve recuento de los principales grupos de investigación y fabricantes a nivel comercial que actualmente trabajan en el desarrollo de esta tecnología. En el capítulo 2 se hacer una breve descripción del proceso de gasificación termoquímica, mostrando los diferentes tipos de gasificadores existentes, las propiedades del combustible primario usado y los factores que afectan la eficiencia de este proceso. En el capítulo 3 se estudia el gas de gasificación desde el punto de vista de la composición, propiedades como combustible motor, requisitos, tratamiento necesario para su uso como combustible en motores de combustión interna, otros usos del GG y riesgos que conlleva su utilización. En el capítulo 4 se presentan un estudio general de los motores de gas, donde se presenta una clasificación, se estudia la manera de regular la operación de estos motores, se hace una descripción cualitativa de la combustión, se muestran algunas aplicaciones y se estudia la combustión en los motores a gas desde el punto de vista de los factores que la afectan. En el capítulo 5 se entra en profundidad sobre el uso del gas de gasificación -GG en MCIA, inicialmente estudiado desde el punto de vista teórico, luego presentando los resultados de varias investigaciones realizadas y por último mostrando algunos de las aplicaciones comerciales actualmente en el mercado. Finalmente se presentan las conclusiones, la bibliografía consultada y un glosario de términos.

Sistema de conversión electromecánica de alta potencia específica para generación eléctrica de origen renovable

La presente tesis estudia, analiza y explora la problemática de la conversión electromecánica a baja velocidad y propone soluciones para mejorarla en términos de prestaciones y coste. La baja velocidad es una característica común de la mayor parte de sistemas de conversión de fuentes renovables. Dicha baja velocidad se debe a la manera en que se establece la interacción entre la fuente y el sistema de generación. Se aborda la problemática del accionamiento directo entre el absorbedor y el generador eléctrico en los casos de energía eólica y energía con olas. A continuación se profundiza en las máquinas de imanes permanentes, las cuales son de gran interés en este caso por las prestaciones que permiten en términos de eficiencia y densidad de potencia. Se presentan los modelos físicos empleados y en base a ellos se proponen metodologías de diseño, una de las cuales maximiza el par en relación a las pérdidas, y otra que minimiza el efecto de la incertidumbre existente sobre el precio de los imanes permanentes, problemática enorme en el caso de las tierras raras. Adicionalmente, y como contribución a la técnica más importante de esta tesis, se propone y justifica el uso en aplicaciones de accionamiento directo de una configuración de máquina que mejora las características de las ya existentes.

Hidrógeno electrolítico para la gestión de la red eléctrica

La variabilidad de la demanda a lo largo de las horas del día y las estaciones del año, junto a la generación aleatoria de la mayoría de las fuentes de energía renovable, sitúan al sistema eléctrico ante un nuevo y fascinante paradigma; ello resulta especialmente relevante en redes relativamente aisladas y con una notable penetración de los parques eólicos, como sucede actualmente en España. La figura 1 muestra una curva típica, con dos puntas que reflejan sobretodo las necesidades de climatización durante las horas cálidas del estío e iluminación en los días cortos del invierno.

Sistemas de Gestión de la Calidad. Actividades del proceso de mejora continua: Estudio y Análisis de los defectos de construcción en edificación de viviendas

Uno de los pilares fundamentales que defienden el grupo de Normas ISO 9000 para garantizar la eficacia y eficiencia de los Sistemas de Gestión de la Calidad es la Mejora Continua. Para llevar a cabo el proceso de Mejora Continua es preciso empezar detectando y analizando donde estamos fallando y qué podemos mejorar. Organizaciones como PMI® asignan a esta actividad un proceso propio, que llaman también Lecciones Aprendidas y debe estar presente siempre, en todos los proyectos, ya que sólo conociendo nuestros fallos y analizando sus causas podemos tomar acciones preventivas que eviten que vuelvan a repetirse. Desafortunadamente, la experiencia práctica nos dice que, en la mayoría de los casos, esta actividad o no se ejecuta, o no se realiza de forma eficaz. La primera actividad que debería llevarse a cabo para detectar y analizar donde estamos fallando es la inspección del producto, en nuestro caso el edificio de viviendas, y realizar un listado de los defectos encontrados. Esta tarea, que debería realizarla la empresa constructora y ser supervisada por la Dirección de Ejecución de Obra, generalmente, no se realiza con el rigor suficiente que permita entregar al usuario su vivienda sin defectos aparentes. Se consideran repasos de postventa, dejando en el usuario final la responsabilidad de verificar el correcto acabado del producto que acaba de adquirir y que en la mayoría de los casos será la compra de mayor importe que realice en su vida. Cuando adquirimos cualquier objeto, un coche, un televisor, un teléfono móvil, etc., el fabricante es el responsable de verificar que el producto que está entregando al cliente o al suministrador no tiene daños o defectos aparentes. Sin embargo, cuando compramos una vivienda es el comprador o usuario final el que realiza, en la mayoría de los casos, esta verificación.

Libertad, responsabilidad y técnica

¿Qué han hecho los ciudadanos que se desplazan en coche para que tantos grupúsculos, grupos, organismos y medios estén empeñados en restringir este estupendo modo de desplazamiento? ¿Están justificados tantos ataques? ¿Hay recursos técnicos para ordenar la situación, ajustando la libertad a sus límites y dando y pidiendo responsabilidades a los conductores, a otros actores y a la Administración? Para empezar, necesitamos términos nuevos. Acompañante no es un término correcto para un viajero que no conduce. Parecería que, no teniendo nada mejor que hacer, se mete en el coche con el conductor (un estudio reseñado en La Vanguardia afirma que el 33,3% % de los desplazamientos de las mujeres catalanas se hacen en coche, pero ellas prefieren viajar sin conducir). Por su parte, conductores son unos cuantos, pero esto queda bastante lejos de todos los que manejan un coche. La misma denominación de conductor trae a la cabeza a alguien cuyo primer objetivo es conducir, y esto no suele ser cierto: es desplazarse con otro motivo, no es pasear en coche sin más. Cochero no parece un término adecuado. Viajero individual tampoco. Pero necesitamos que el lenguaje refleje la realidad: los que llamamos conductores son viajeros, o pasajeros (aunque no compren un pasaje).

La absorción acústica a ruido de impacto en placas de aislamiento térmico de suelo radiante

Actualmente en el mercado existen una reducida variedad de materiales específicos de ruido de impactos. Los materiales que han de ser utilizados como material resiliente, han de cumplir unos valores muy bajos del módulo de elasticidad E. Mientras el valor del caucho blando está en el orden de los 100 kg/cm2, el valor del acero es de 2.100.000 kg/cm2 (21.000 veces mayor). Estos valores tan bajos hacen que diferencias de apenas 50 kg/cm2 produzcan variaciones muy importantes en las ganancias de aislamiento acústico. La realización de los ensayos de caracterización acústica a ruido de impacto en laboratorio es un proceso largo y caro. En primer lugar se necesita elaborar un forjado de referencia, que ha de ser ensayado de forma independiente. Con posterioridad, se ha de colocar la losa flotante, y volver a realizar la prueba. Con la diferencia de las medidas realizadas, se determinan los valores de módulo de elasticidad. Como alternativa a este método, en la ponencia del congreso “TECNIACÚSTICA GANDÍA 2006” sobre “estudio de cinco métodos para determinar las propiedades dinámicas de capas elásticas para la mejora del aislamiento a ruido de impactos”, presentada por Francisco Simón, David K. Anthony, Mª José Fernández, proponen un total de cinco alternativas al ensayo indicado en la norma. Al margen de las condiciones acústicas, las placas de suelo radiante necesitan garantizar unos valores mínimos de aislamiento térmico. Con las actuales necesidades normativas, la solución habitual consiste en la superposición de dos materiales: en primer lugar, inmediatamente sobre el forjado, se coloca el material encargado de proporcionar el aislamiento acústico, generalmente una lámina flexible de polietileno. Encima de ella se coloca la placa de aislamiento térmico, generalmente una plancha de poliestireno expandido con resaltes en las que se encajan las tuberías que conducen el agua calefactada. El poliestireno expandido convencional, no tiene unos valores de módulo elástico suficientemente bajo para poder actuar como material resiliente. Para conseguir estos valores, es necesario elastificar el material mediante su introducción en una prensa (Poliestireno expandido elastificado EEPS). La utilización de esta técnica es incompatible con el proceso de fabricación de planchas moldeadas con tetones para la integración de los soportes de las tuberías. La lámina flexible de polietileno, si tienes valores aceptables, tanto de aislamiento térmico como acústico, en cambio no es posible la integración de los tetones durante su proceso de fabricación para soportar las canalizaciones de agua. El objetivo de la investigación consiste en integrar en una única placa las características tanto térmicas como acústicas necesarias para satisfacer todas las necesidades. Como se ha indicado más arriba, los materiales resilientes han de cumplir unos valores muy bajos del módulo de elasticidad E. Mientras el valor del caucho blando está en el orden de los 100 kg/cm2, el valor del acero es de 2.100.000 kg/cm2 (21.000 veces mayor). Pero no es menos cierto que con acero, dando la forma adecuada (un helicoide) se consigue un elemento muy elástico, un muelle. Tal es así, que variando la forma se puede conseguir el módulo de elásticidad E deseado, resultando utilidades tan dispares como el resorte de un bolígrafo o la suspensión de un coche. Por otro lado, y como se ha indicado anteriormente, la determinación de los valores mecánicos de los materiales aislantes térmicos es cara y complicada. Como objetivo secundario de la investigación, se establece el proponer un ensayo alternativo, más rápido y económico que el propuesto por la norma para facilitar el estudio de nuevas alternativas, sin tener que hacer grandes desembolsos económicos en

Optimization of urban mobility measures to achieve win-win strategies

European cities are essential in the development of Europe as they constitute the living environment of more than 60% of the population in the European Union and are drivers of the European economy – just under 85% of the EU’s gross domestic product is produced in urban areas (EC, 2007a). The car has been one of the main factors of development during the 20th century, but it is at the same time the origin of the key problems cities have to face: traffic increase. This has resulted in chronic congestion with many adverse consequences such as air pollution and noise. This loss of environmental quality is one of the reasons for urban sprawl in European cities during recent decades. But this urban sprawl at the same time worsens the environmental conditions. We must return to the dense city, but clean and competitive, and this implies reducing car use yet provides quality transport alternatives sufficient to recover and maintain the competitiveness of cities (EC, 2007a). Consequently, European cities need to establish an urban transport strategy which helps reduce their environmental problems –mainly emissions and noise – but without decreasing their trip attraction. This aspect is very important because a loss of trip attraction would result in an increase of people moving to more disperse areas, contributing towards worsening the current situation. This thesis is an attempt to contribute solutions to this problem in two ways: 1) The first is to analyze the complementarity and possible synergies of several urban transport measures aimed at improving a modal split to a more sustainable means of transport. This analysis will focus on the three aspects already mentioned: emissions, noise and attractiveness or competitiveness. 2) Once possible synergies and complementarities have been analyzed, the second objective is to propose the best combination of these measures, in terms of level of implementation, to achieve the maximum benefit with respect to the three aspects previously established: emissions, noise and attractiveness or competitiveness. Therefore, within the wide range of measures enhancing sustainable urban transport, three of them have been be selected in this thesis to establish a methodology for achieving these objectives. The analysis will be based on the region of Madrid, which is also the case study selected for this research. Las ciudades europeas son piezas fundamentales para el desarrollo europeo, ya que son el lugar de residencia de más del 60% de la población de la unión europea así como los motores de su economía – casi el 85% del PIB europeo se produce en áreas urbanas (EC, 2007a). El coche ha sido uno de los principales motores de desarrollo de las ciudades durante el siglo XX, pero se ha terminado por convertir a su vez en uno de los principales problemas con los que tiene que lidiar las ciudades: el aumento del tráfico. Esto ha derivado en unos niveles crónicos de congestión, con multitud de efectos adversos, entre los que cabe destacar la contaminación del aire y el ruido. Esta pérdida de calidad ambiental es una de las razones que ha propiciado la dispersión urbana que han experimentado las ciudades europeas en las últimas décadas. Pero esta dispersión urbana a su vez contribuye a empeorar las condiciones ambientales de las ciudades. Debemos retornar a la ciudad densa, pero limpia y competitiva, y esto implica reducir el uso del coche, pero proporcionando alternativas de transporte que permitan recuperar y mantener la competitividad de las ciudades (EC, 2007a). Por lo tanto, las ciudades europeas necesitan encontrar una estrategia de transporte urbano que ayude a reducir sus problemas medio ambientales – principalmente ruido y emisiones – pero sin hacerlas perder atractividad o competitividad. Este aspecto tiene gran importancia porque una pérdida de la misma se traduciría en un aumento de dispersión de la población hacia áreas periféricas, contribuyendo a empeorar la situación actual. Esta tesis contribuye a solucionar este problema de dos maneras: 1) La primera, analizando la complementariedad y posibles sinergias de diferentes medidas de transporte urbano orientadas a promover un reparto modal hacia modos más sostenibles. Este análisis se centrará en los tres aspectos anteriormente citados: emisiones, ruido y atractividad o competitividad. 2) Una vez las posibles sinergias y complementariedades se han analizado, el segundo objetivo es proponer la mejor combinación de estas medidas – en términos de grado de aplicación - para lograr el máximo beneficio en lo que respecta a los tres objetivos previamente establecidos. Para ello, en esta tesis se han seleccionado una serie de medidas que permitan establecer una metodología para alcanzar estos objetivos previamente definidos. El análisis se centra en la ciudad de Madrid y su área metropolitana, la cual se ha escogido como caso de estudio para realizar esta investigación.

Asistencia técnica e informe final sobre la demolición del gran contrafuerte añadido a la iglesia de Melón en el antiguo claustro del Monasterio (Melón, Ourense)

El monasterio de Santa María de Melón, fundado en el siglo XII, fue abandonado tras la desamortización de Mendizábal en 1835. Las inclemencias del tiempo y el expolio lo han convertido en una ruina. La ruina ha sido consolidada y es visitable gracias a dos intervenciones realizadas en 2003 y 2010, respectivamente. La iglesia, sin embargo, ha continuado en uso hasta la actualidad. La cabecera es románica y siguió una ampliación gótica, que añadió el crucero y las tres naves. Distintas intervenciones se sucedieron hasta el siglo XVI (para una historia detallada del templo, véase Fernández Rodríguez 2010). A finales del siglo XIX por motivos que se desconocen se hundió toda la parte de los pies de la iglesia, esto es, las tres naves desde el crucero hasta los pies. Se desconocen las causas del hundimiento (no podemos estar de acuerdo con la interpretación de Fernández Rodríguez (2010, 45) que atribuye el hundimiento a una tormenta con aparato eléctrico y vendaval). Sea como fuere, con fecha desconocida se decidió dejar la iglesia ápoda, sin naves, y rematarla con un cuerpo de fachada que está formado por dos cubos macizos entre los que se sitúa la fachada y, en el interior, el coro, Figura 1, A y B.. (La bóveda sobre el coro situada entre ellos debió ser reconstruida.) La particularidad más interesante de esta intervención es que buena parte de los restos de la ruina fueron almacenados en dichos cubos a los pies de la iglesia. Esto se puede verificar en el cubo A de la izquierda mirando a través de las ventanas del husillo de la escalera, cegadas al construir el cubo. Se hizo una cata a unos 4 m de altura, desde el coro, al cubo B de la derecha con el mismo resultado. El cubo A tiene un volumen aproximado de 6×5×12 m3 de volumen. El cubo B tiene una planta algo menor y un volumen de 6×4×12 m3. En la hipótesis de que ambos cubos están rellenos hasta su coronación esto daría un volumen de fábrica de unos 600 m3 (esto parece probable a la luz de la información disponible sobre el cubo izquierdo).

Código híbrido avanzado de motores de plasma de efecto Hall

La propulsión eléctrica constituye hoy una tecnología muy competitiva y de gran proyección de futuro. Dentro de los diversos motores de plasma existentes, el motor de efecto Hall ha adquirido una gran madurez y constituye un medio de propulsión idóneo para un rango amplio de misiones. En la presente Tesis se estudian los motores Hall con geometría convencional y paredes dieléctricas. La compleja interacción entre los múltiples fenómenos físicos presentes hace que sea difícil la simulación del plasma en estos motores. Los modelos híbridos son los que representan un mejor compromiso entre precisión y tiempo de cálculo. Se basan en utilizar un modelo fluido para los electrones y algoritmos de dinámica de partículas PIC (Particle-In- Cell) para los iones y los neutros. Permiten hacer uso de la hipótesis de cuasineutralidad del plasma, a cambio de resolver separadamente las capas límite (o vainas) que se forman en torno a las paredes de la cámara. Partiendo de un código híbrido existente, llamado HPHall-2, el objetivo de la Tesis doctoral ha sido el desarrollo de un código híbrido avanzado que mejorara la simulación de la descarga de plasma en un motor de efecto Hall. Las actualizaciones y mejoras realizadas en las diferentes partes que componen el código comprenden tanto aspectos teóricos como numéricos. Fruto de la extensa revisión de la algoritmia del código HPHall-2 se han conseguido reducir los errores de precisión un orden de magnitud, y se ha incrementado notablemente su consistencia y robustez, permitiendo la simulación del motor en un amplio rango de condiciones. Algunos aspectos relevantes a destacar en el subcódigo de partículas son: la implementación de un nuevo algoritmo de pesado que permite determinar de forma más precisa el flujo de las magnitudes del plasma; la implementación de un nuevo algoritmo de control de población, que permite tener suficiente número de partículas cerca de las paredes de la cámara, donde los gradientes son mayores y las condiciones de cálculo son más críticas; las mejoras en los balances de masa y energía; y un mejor cálculo del campo eléctrico en una malla no uniforme. Merece especial atención el cumplimiento de la condición de Bohm en el borde de vaina, que en los códigos híbridos representa una condición de contorno necesaria para obtener una solución consistente con el modelo de interacción plasma-pared, y que en HPHall-2 aún no se había resuelto satisfactoriamente. En esta Tesis se ha implementado el criterio cinético de Bohm para una población de iones con diferentes cargas eléctricas y una gran dispersión de velocidades. En el código, el cumplimiento de la condición cinética de Bohm se consigue por medio de un algoritmo que introduce una fina capa de aceleración nocolisional adyacente a la vaina y mide adecuadamente el flujo de partículas en el espacio y en el tiempo. Las mejoras realizadas en el subcódigo de electrones incrementan la capacidad de simulación del código, especialmente en la región aguas abajo del motor, donde se simula la neutralización del chorro del plasma por medio de un modelo de cátodo volumétrico. Sin abordar el estudio detallado de la turbulencia del plasma, se implementan modelos sencillos de ajuste de la difusión anómala de Bohm, que permiten reproducir los valores experimentales del potencial y la temperatura del plasma, así como la corriente de descarga del motor. En cuanto a los aspectos teóricos, se hace especial énfasis en la interacción plasma-pared y en la dinámica de los electrones secundarios libres en el interior del plasma, cuestiones que representan hoy en día problemas abiertos en la simulación de los motores Hall. Los nuevos modelos desarrollados buscan una imagen más fiel a la realidad. Así, se implementa el modelo de vaina de termalización parcial, que considera una función de distribución no-Maxwelliana para los electrones primarios y contabiliza unas pérdidas energéticas más cercanas a la realidad. Respecto a los electrones secundarios, se realiza un estudio cinético simplificado para evaluar su grado de confinamiento en el plasma, y mediante un modelo fluido en el límite no-colisional, se determinan las densidades y energías de los electrones secundarios libres, así como su posible efecto en la ionización. El resultado obtenido muestra que los electrones secundarios se pierden en las paredes rápidamente, por lo que su efecto en el plasma es despreciable, no así en las vainas, donde determinan el salto de potencial. Por último, el trabajo teórico y de simulación numérica se complementa con el trabajo experimental realizado en el Pnnceton Plasma Physics Laboratory, en el que se analiza el interesante transitorio inicial que experimenta el motor en el proceso de arranque. Del estudio se extrae que la presencia de gases residuales adheridos a las paredes juegan un papel relevante, y se recomienda, en general, la purga completa del motor antes del modo normal de operación. El resultado final de la investigación muestra que el código híbrido desarrollado representa una buena herramienta de simulación de un motor Hall. Reproduce adecuadamente la física del motor, proporcionando resultados similares a los experimentales, y demuestra ser un buen laboratorio numérico para estudiar el plasma en el interior del motor. Abstract Electric propulsion is today a very competitive technology and has a great projection into the future. Among the various existing plasma thrusters, the Hall effect thruster has acquired a considerable maturity and constitutes an ideal means of propulsion for a wide range of missions. In the present Thesis only Hall thrusters with conventional geometry and dielectric walls are studied. The complex interaction between multiple physical phenomena makes difficult the plasma simulation in these engines. Hybrid models are those representing a better compromise between precision and computational cost. They use a fluid model for electrons and Particle-In-Cell (PIC) algorithms for ions and neutrals. The hypothesis of plasma quasineutrality is invoked, which requires to solve separately the sheaths formed around the chamber walls. On the basis of an existing hybrid code, called HPHall-2, the aim of this doctoral Thesis is to develop an advanced hybrid code that better simulates the plasma discharge in a Hall effect thruster. Updates and improvements of the code include both theoretical and numerical issues. The extensive revision of the algorithms has succeeded in reducing the accuracy errors in one order of magnitude, and the consistency and robustness of the code have been notably increased, allowing the simulation of the thruster in a wide range of conditions. The most relevant achievements related to the particle subcode are: the implementation of a new weighing algorithm that determines more accurately the plasma flux magnitudes; the implementation of a new algorithm to control the particle population, assuring enough number of particles near the chamber walls, where there are strong gradients and the conditions to perform good computations are more critical; improvements in the mass and energy balances; and a new algorithm to compute the electric field in a non-uniform mesh. It deserves special attention the fulfilment of the Bohm condition at the edge of the sheath, which represents a boundary condition necessary to match consistently the hybrid code solution with the plasma-wall interaction, and remained as a question unsatisfactory solved in the HPHall-2 code. In this Thesis, the kinetic Bohm criterion has been implemented for an ion particle population with different electric charges and a large dispersion in their velocities. In the code, the fulfilment of the kinetic Bohm condition is accomplished by an algorithm that introduces a thin non-collisional layer next to the sheaths, producing the ion acceleration, and measures properly the flux of particles in time and space. The improvements made in the electron subcode increase the code simulation capabilities, specially in the region downstream of the thruster, where the neutralization of the plasma jet is simulated using a volumetric cathode model. Without addressing the detailed study of the plasma turbulence, simple models for a parametric adjustment of the anomalous Bohm difussion are implemented in the code. They allow to reproduce the experimental values of the plasma potential and the electron temperature, as well as the discharge current of the thruster. Regarding the theoretical issues, special emphasis has been made in the plasma-wall interaction of the thruster and in the dynamics of free secondary electrons within the plasma, questions that still remain unsolved in the simulation of Hall thrusters. The new developed models look for results closer to reality, such as the partial thermalization sheath model, that assumes a non-Maxwellian distribution functions for primary electrons, and better computes the energy losses at the walls. The evaluation of secondary electrons confinement within the chamber is addressed by a simplified kinetic study; and using a collisionless fluid model, the densities and energies of free secondary electrons are computed, as well as their effect on the plasma ionization. Simulations show that secondary electrons are quickly lost at walls, with a negligible effect in the bulk of the plasma, but they determine the potential fall at sheaths. Finally, numerical simulation and theoretical work is complemented by the experimental work carried out at the Princeton Plasma Physics Laboratory, devoted to analyze the interesting transitional regime experienced by the thruster in the startup process. It is concluded that the gas impurities adhered to the thruster walls play a relevant role in the transitional regime and, as a general recomendation, a complete purge of the thruster before starting its normal mode of operation it is suggested. The final result of the research conducted in this Thesis shows that the developed code represents a good tool for the simulation of Hall thrusters. The code reproduces properly the physics of the thruster, with results similar to the experimental ones, and represents a good numerical laboratory to study the plasma inside the thruster.

Electromagnetic analysis of periodic structures under a waveguide viewpoint

The fundamental objective of this Ph. D. dissertation is to demonstrate that, under particular circumstances which cover most of the structures with practical interest, periodic structures can be understood and analyzed by means of closed waveguide theories and techniques. To that aim, in the first place a transversely periodic cylindrical structure is considered and the wave equation, under a combination of perfectly conducting and periodic boundary conditions, is studied. This theoretical study runs parallel to the classic analysis of perfectly conducting closed waveguides. Under the light shed by the aforementioned study it is clear that, under certain very common periodicity conditions, transversely periodic cylindrical structures share a lot of properties with closed waveguides. Particularly, they can be characterized by a complete set of TEM, TE and TM modes. As a result, this Ph. D. dissertation introduces the transversely periodic waveguide concept. Once the analogies between the modes of a transversely periodic waveguide and the ones of a closed waveguide have been established, a generalization of a well-known closed waveguide characterization method, the generalized Transverse Resonance Technique, is developed for the obtention of transversely periodic modes. At this point, all the necessary elements for the consideration of discontinuities between two different transversely periodic waveguides are at our disposal. The analysis of this type of discontinuities will be carried out by means of another well known closed waveguide method, the Mode Matching technique. This Ph. D. dissertation contains a sufficient number of examples, including the analysis of a wire-medium slab, a cross-shaped patches periodic surface and a parallel plate waveguide with a textured surface, that demonstrate that the Transverse Resonance Technique - Mode Matching hybrid is highly precise, efficient and versatile. Thus, the initial statement: ”periodic structures can be understood and analyzed by means of closed waveguide theories and techniques”, will be corroborated. Finally, this Ph. D. dissertation contains an adaptation of the aforementioned generalized Transverse Resonance Technique by means of which the analysis of laterally open periodic waveguides, such as the well known Substrate Integrated Waveguides, can be carried out without any approximation. The analysis of this type of structures has suscitated a lot of interest in the recent past and the previous analysis techniques proposed always resorted to some kind of fictitious wall to close the structure. vii Resumen El principal objetivo de esta tesis doctoral es demostrar que, bajo ciertas circunstancias que se cumplen para la gran mayoría de estructuras con interés práctico, las estructuras periódicas se pueden analizar y entender con conceptos y técnicas propias de las guías de onda cerradas. Para ello, en un primer lugar se considera una estructura cilíndrical transversalmente periódica y se estudia la ecuación de onda bajo una combinación de condiciones de contorno periódicas y de conductor perfecto. Este estudio teórico y de caracter general, sigue el análisis clásico de las guías de onda cerradas por conductor eléctrico perfecto. A la luz de los resultados queda claro que, bajo ciertas condiciones de periodicidad (muy comunes en la práctica) las estructuras cilíndricas transversalmente periódicas guardan multitud de analogías con las guías de onda cerradas. En particular, pueden ser descritas mediante un conjunto completo de modos TEM, TE y TM. Por ello, ésta tesis introduce el concepto de guía de onda transversalmente periódica. Una vez establecidas las similitudes entre las soluciones de la ecuación de onda, bajo una combinación de condiciones de contorno periódicas y de conductor perfecto, y los modos de guías de onda cerradas, se lleva a cabo, con éxito, la adaptación de un conocido método de caracterización de guías de onda cerradas, la técnica de la Resonancia Transversal Generalizada, para la obtención de los modos de guías transversalmente periódicas. En este punto, se tienen todos los elementos necesarios para considerar discontinuidades entre guías de onda transversalmente periódicas. El analisis de este tipo de discontinuidades se llevará a cabo mediante otro conocido método de análisis de estructuras cerradas, el Ajuste Modal. Esta tesis muestra multitud de ejemplos, como por ejemplo el análisis de un wire-medium slab, una superficie de parches con forma de cruz o una guía de placas paralelas donde una de dichas placas tiene cierta textura, en los que se demuestra que el método híbrido formado por la Resonancia Transversal Generalizada y el Ajuste Modal, es tremendamente preciso, eficiente y versátil y confirmará la validez de el enunciado inicial: ”las estructuras periódicas se pueden analizar y entender con conceptos y técnicas propias de las guías de onda cerradas” Para terminar, esta tésis doctoral incluye también una modificación de la técnica de la Resonancia Transversal Generalizada mediante la cual es posible abordar el análisis de estructuras periódica abiertas en los laterales, como por ejemplo las famosas guías de onda integradas en sustrato, sin ninguna aproximación. El análisis de este tipo de estructuras ha despertado mucho interés en los últimos años y las técnicas de análisis propuestas hasta ix el momento acostumbran a recurrir a algún tipo de pared ficticia para simular el carácter abierto de la estructura.