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.

Supersonic regime of the Hall-magnetohydrodynamics resistive tearing instability

An earlier analysis of the Hall-magnetohydrodynamics (MHD) tearing instability [E. Ahedo and J. J. Ramos, Plasma Phys. Controlled Fusion 51, 055018 (2009)] is extended to cover the regime where the growth rate becomes comparable or exceeds the sound frequency. Like in the previous subsonic work, a resistive, two-fluid Hall-MHD model with massless electrons and zero-Larmor-radius ions is adopted and a linear stability analysis about a force-free equilibrium in slab geometry is carried out. A salient feature of this supersonic regime is that the mode eigenfunctions become intrinsically complex, but the growth rate remains purely real. Even more interestingly, the dispersion relation remains of the same form as in the subsonic regime for any value of the instability Mach number, provided only that the ion skin depth is sufficiently small for the mode ion inertial layer width to be smaller than the macroscopic lengths, a generous bound that scales like a positive power of the Lundquist number

Neutron features at the UPM neutronics hall

The neutronics hall of the Nuclear Engineering Department at the Polytechnical University of Madrid has been characterized. The neutron spectra and the ambient dose equivalent produced by an 241AmBe source were measured at various source-to-detector distances on the new bench. Using Monte Carlo methods a detailed model of the neutronics hall was designed, and neutron spectra and the ambient dose equivalent were calculated at the same locations where measurements were carried out. A good agreement between measured and calculated values was found.

Teoría y práctica del Generador Hall

El efecto Hall fue descubierto en 1879 por el físico americano del mismo nombre. Durante muchos años ha permanecido en la oscuridad hasta que, históricamente,ha adquirido importancia en el dominio de la física de semiconductores porque condujo a postular la existencia de los huecos y por tanto a construir nuevas teorías no basadas exclusivamente en el flujo electrónico. Actualmente es una prueba corriente en los laboratorios para la determinación de algunas características de los materiales.

Dispositivos basados en el efecto Hall. Parte I

Se presenta aquí un elemento pasivo, de reciente desarrollo, el generador Hall. Este elemento de pequeñísimo tamaño, se basa en un efecto descubierto en 1789 por Hall. En primer lugar, se estudia el efecto, que tiene lugar en los metales, en los semiconductores e incluso en los gases ionizados. El coeficiente de Hall Rh o bien la sensitividad K0 representan en cierto modo la, capacidad de funcionamiento del generador, y por esto resulta muy interesante profundizar un poco en ello y ver qué factores influyen sobre los valores de aquellos parámetros. Se citan algunas de las muchas propiedades de aplicación del elemento. Para terminar esta primera parte, dedicada al efecto en los semiconductores, presentamos un circuito equivalente, cuyo margen de validez vendrá en cada caso determinado por los valores máximos establecidos por el fabricante. Por último, se expone el principio del multiplicador, y sobre un multiplicador comercial se tratan de explicar algunos resultados del mismo que parecen pobres, dada la posibilidad teórica de la cápsula en sí. Aunque comprendemos que gran parte de los inconvenientes provienen de los mismos circuitos, como no conocemos éstos,nos ceñimos más bien a los errores que provienen de la propia cápsula. Parece más honrado y quizás añada algo de luz sobre el conocimiento del elemento considerado

Dispositivos basados en el efecto Hall. Parte II

La medida de un campo magnético con un generador Hall reposa sobre el hecho de que, si se mantiene la corriente de control constante y se consigue una adaptación lineal, la tensión Hall es directamente proporcional a la componente del campo magnético normal a la pastilla semiconductora. La idea de aplicar el efecto Hall a la medida de los campos magnéticos es antigua, pero la resolución práctica sólo ha sido posible a partir de la consecución de las uniones inter metálicas, hacia el año 1952.

Structure of intermediate shocks in collisionsless anisotropic Hall-magnetohydrodynamics plasma models

The existence of discontinuities within the double-adiabatic Hall-magnetohydrodynamics (MHD) model is discussed. These solutions are transitional layers where some of the plasma properties change from one equilibrium state to another. Under the assumption of traveling wave solutions with velocity C and propagation angle θ with respect to the ambient magnetic field, the Hall-MHD model reduces to a dynamical system and the waves are heteroclinic orbits joining two different fixed points. The analysis of the fixed points rules out the existence of rotational discontinuities. Simple considerations about the Hamiltonian nature of the system show that, unlike dissipative models, the intermediate shock waves are organized in branches in parameter space, i.e., they occur if a given relationship between θ and C is satisfied. Electron-polarized (ion-polarized) shock waves exhibit, in addition to a reversal of the magnetic field component tangential to the shock front, a maximum (minimum) of the magnetic field amplitude. The jumps of the magnetic field and the relative specific volume between the downstream and the upstream states as a function of the plasma properties are presented. The organization in parameter space of localized structures including in the model the influence of finite Larmor radius is discussed

Extension of Euskalduna Conference Centre and Concert Hall: A Contemporary Application of Irregular Reciprocal Frames

Reciprocal frame structures, formed by a set of self-supported elements in a closed circuit, have long been used since antiquity to cover large spans with small elements. The roof structure of the Euskalduna conference centre and concert hall extension in Bilbao, covering an irregu- lar geometry of 3000 m2 with a maximum span of 45 m, presented an interesting opportunity to revisit the concept and to apply these classical systems. Furthermore, its analysis and develop- ment led to an interesting discussion on reciprocal frames. They showed great sensitivity of these systems to the local modification of a particular element, establishment of irregular load paths, mobilisation of almost the entire sys- tem when locally applying a punctual load and, finally, its large deformability. Besides, reciprocal frames present particular construction complexities and possibilities due to the moderate length of the structural elements, the predominance of shear-only connec- tions and the necessity of the entire system to be completely erected to guarantee its stability. Euskalduna extension, completed in 2012, is one of the largest and a very par- ticular case of irregular reciprocal frame structures built in the world. It shows the formal possibilities and potentiality of reciprocal frames to respond to free and irregular geometries.

City Hall, Fene, La Coruña (1981)

Ayuntamiento de Fene, La Coruña

SOPHIA CPV module round robin: power rating at CSOC

Within the European funded project SOPHIA, a Round Robin measurement on CPV module has been initiated. Seven different test laboratories located in Europe between 48°N and 37°N perform measurements of four SOITEC CPV modules. The modules are electrically characterized with different measurement equipment under various climatic conditions. One pyrheliometer and one spectral sensor based on component cells are shipped together with the modules. This ensures that the irradiance and spectrum, two factors with high impact on CPV module performance, are measured with the identical equipment at each site. The round robin activity is performed in closeco-operation with the IEC TC82 WG7 power rating team in order to support the work on the CPV module power rating draft standard 62670-3. The resultingrated module power outputs at CSOC (Concentrator Standard Operating Conditions) are compared amongst the power rating methods and amongst the test labs. In this manner, a deviation in rated power output between different test labs and power rating methods is determined.

Electron transport and azimuthal oscillations in hall thrusters

The Hall Effect Thruster (HET) is a type of satellite electric propulsion device initially developed in the 1960’s independently by USA and the former USSR. The development continued in the shadow during the 1970’s in the Soviet Union to reach a mature status from the technological point of view in the 1980’s. In the 1990’s the advanced state of this Russian technology became known in western countries, which rapidly restarted the analysis and development of modern Hall thrusters. Currently, there are several companies in USA, Russia and Europe manufacturing Hall thrusters for operational use. The main applications of these thrusters are low-thrust propulsion of interplanetary probes, orbital raising of satellites and stationkeeping of geostationary satellites. However, despite the well proven in-flight experience, the physics of the Hall Thruster are not completely understood yet. Over the last two decades large efforts have been dedicated to the understanding of the physics of Hall Effect thrusters. However, the so-called anomalous diffusion, short name for an excessive electron conductivity along the thruster, is not yet fully understood as it cannot be explained with classical collisional theories. One commonly accepted explanation is the existence of azimuthal oscillations with correlated plasma density and electric field fluctuations. In fact, there is experimental evidence of the presence of an azimuthal oscillation in the low frequency range (a few kHz). This oscillation, usually called spoke, was first detected empirically by Janes and Lowder in the 1960s. More recently several experiments have shown the existence of this type of oscillation in various modern Hall thrusters. Given the frequency range, it is likely that the ionization is the cause of the spoke oscillation, like for the breathing mode oscillation. In the high frequency range (a few MHz), electron-drift azimuthal oscillations have been detected in recent experiments, in line with the oscillations measured by Esipchuk and Tilinin in the 1970’s. Even though these low and high frequency azimuthal oscillations have been known for quite some time already, the physics behind them are not yet clear and their possible relation with the anomalous diffusion process remains an unknown. This work aims at analysing from a theoretical point of view and via computer simulations the possible relation between the azimuthal oscillations and the anomalous electron transport in HET. In order to achieve this main objective, two approaches are considered: local linear stability analyses and global linear stability analyses. The use of local linear stability analyses shall allow identifying the dominant terms in the promotion of the oscillations. However, these analyses do not take into account properly the axial variation of the plasma properties along the thruster. On the other hand, global linear stability analyses do account for these axial variations and shall allow determining how the azimuthal oscillations are promoted and their possible relation with the electron transport.

El lenguaje de la estructura: Mies van der Rohe y la construcción con acero y hormigón

John Summerson, en El Lenguaje Clásico de la Arquitectura, defiende que los órdenes clásicos que empleaban los arquitectos romanos para decorar sus edificios, no tienen una función estructural pero hacen expresivos a los edificios. Les hacen hablar. Arthur Schopenhauer afirma que “el destino de la Bella Arquitectura es poner de manifiesto la lucha entre el peso y la rigidez de los elementos estructurales”. Y Auguste Perret define la Arquitectura como “el arte de hacer cantar al punto de apoyo”. El objetivo de esta Tesis Doctoral es profundizar en la capacidad de expresión de la estructura. A través del estudio de las estructuras históricas, que se realiza en la primera parte de la Tesis, podemos concluir que existen tres categorías, tres maneras de expresión de la Estructura. Estructuras Vistas, que hablan, Estructuras Ocultas, que se esconden y Estructuras Ilusorias, que fingen. El Partenón y la Sainte Chapelle de París se estudian en el apartado de Estructuras Vistas. El Panteón, el Palacio de Carlos V, la Catedral de San Pablo en Londres, y otras arquitecturas renacentistas y romanas, en el apartado de Estructuras Ocultas. Y como Estructuras Ilusorias, la Alhambra (Dos Hermanas, Comares, el Patio de los Leones), Santa Sofía, y otras arquitecturas del barroco italiano. En la segunda parte de la Tesis se analiza la obra completa de Mies van der Rohe desde el punto de vista de esas tres categorías. Lo visto, lo oculto y lo ilusorio en las estructuras de Mies. Se estudia la evolución en la estructura de la casa, desde las primeras casas con Estructura Oculta de muro de ladrillo, hasta las últimas casas con Estructura Vista y columnas adelantadas, pasando por una etapa intermedia de casas con estructura mixta de muro de ladrillo en la que el acero comienza a hacer su aparición. Se analizan también seis soluciones estructurales en los Bloques y en las Torres: Estructura Vista reverberante, expresiva o inexpresiva vs Estructura Oculta con vestido horizontal, vestido reticular o vestido vertical. Y por último, se estudian las tres soluciones de Estructura Ilusoria que emplea Mies en sus Pabellones. La metodología de trabajo que se ha empleado se divide en cuatro apartados: El análisis bibliográfico; el análisis in situ de los edificios, que nos permite comprobar, por ejemplo, los efectos lumínicos de la columna acanalada del Partenón, o el efecto reverberante de las columnas de la Weissenhofsiedlung; el análisis crítico de planos y detalles constructivos, que nos lleva a concluir que la disposición de pantallas del Pabellón de Barcelona anula la lectura de la crujía estructural, y que la columna del restaurante Cantor se dispone con su alma perpendicular a la cercha, y no paralela, como cabría suponer si se quisiera aprovechar toda la capacidad portante del perfil en H; y por último, el análisis numérico y estructural, que nos lleva a confirmar el sobredimensionado de la Estructura del Patio de los Leones de la Alhambra o el sobredimensionado de la Estructura de la Casa Farnsworth. Lo que se confirma con esta Tesis Doctoral es que la Estructura es algo más, mucho más, que sólo transmisión de las cargas. ABSTRACT John Summerson, in The Classic Language of Architecture, argues that the classic orders used by Roman architects in the decoration of their buildings did not have a structural function, but made buildings expressive. They make them speak. Arthur Schopenhauer affirms that “the goal of Great Architecture is to highlight the struggle between the gravity and rigidity of structural elements”. And Auguste Perret defines Architecture as “the art of making the points of support sing”. The objective of this Doctoral Thesis is to examine the expressive capacity of structure. Following a study of historic structures in the first part of the thesis, we conclude that three categories exist, three ways of expressing Structure. Visible Structures that speak, Concealed Structures that are hidden and Illusory Structures that pretend. The Parthenon and the Sainte Chapelle in Paris are studied in the section on Visible Structures. The Pantheon, the Palace of Charles V, Saint Paul’s Cathedral in London, and other Renaissance and Roman architectures are dealt with in the Concealed Structures section. And, as examples of Illusory Structures, we focus on the Alhambra (The Hall of the Two Sisters, the Comares and the Court of the Lions), Saint Sophia, and other Italian Baroque architectures. In the second part of the Thesis the complete work of Mies van der Rohe is analysed from the perspective of these three categories. The visible, the concealed and the illusory, in the structures of Mies. We study how the structure of the house evolves, from the first houses with the Hidden Structure of the brick wall, to the later houses with Visible Structures and columns, via an intermediate phase of mixed-structure houses with brick walls, where steel first began to make its appearance. We also analyse six structural solutions in the Blocks and Towers: reverberant, expressive or inexpressive Visible Structure vs Concealed Structure with horizontal cladding, reticular cladding or vertical cladding. And finally, we look at the three Illusory Structure solutions that Mies employs in his Pavilions. The methodology employed is divided into four sections: a bibliographic analysis; an analysis in situ of the buildings, which allows us to test, for example, the lighting effects of the fluted column in the Parthenon, or the reverberant effect of the Weissenhofsiedlung columns; a critical analysis of plans and constructive details, which leads us to conclude that the arrangement of panels in the Barcelona Pavilion cancels out the structural centreline, and that the column in the Cantor restaurant is placed with its web perpendicular to the truss, and not parallel to it, as one might expect if one wanted to avail of all the load-bearing capacity of the H beam; and lastly, a numeric and structural analysis, which confirms the oversizing of the Court of the Lions structure in the Alhambra or the oversized structure of Farnsworth House. All of which confirms in this Doctoral Thesis that structure is something more, much more, than a mere conveyor of loads.

SOPHIA CPV module round robin: power rating at CSOC

In the frame of the European project SOPHIA a concentrator photovoltaic (CPV) module measurement round robin has been initiated. The round robin includes measurements of four CPV modules at seven different test laboratories located in Europe. IV curves of the modules are measured with different measurement equipment under various climatic conditions. The aim of this activity is to perform at each site a rating of the modules at concentrator standard operating conditions CSOC according to IEC 62670-1. The outcome of the round robin is intended for direct feedback to the current draft standard IEC 62670-3 “Concentrator Photovoltaic (CPV) Performance Testing - Performance Measurements and Power Rating”. The paper discusses initial results from the first three partners that have already finished the measurements up to now.

Induced radioisotopes in a linac treatment hall

When linacs operate above 8MV an undesirable neutron field is produced whose spectrum has three main components: the direct spectrum due to those neutrons leaking out from the linac head, the scattered spectrum due to neutrons produced in the head that collides with the nuclei in the head losing energy and the third spectrum due to room-return effect. The third category of spectrum has mainly epithermal and thermal neutrons being constant at any location in the treatment hall. These neutrons induce activation in the linac components, the concrete walls and in the patient body. Here the induced radioisotopes have been identified in concrete samples located in the hall and in one of the wedges. The identification has been carried out using a gamma-ray spectrometer.