Nueva formulación para la obtención de la energía cinética de la lluvia. Análisis comparativo con ecuaciones empíricas

El presente trabajo describe el estudio de la caracterización cuantitativa de una propiedad de la lluvia muy importante, la energía cinética, y es importante porque es la fuente de entrada de energía en los procesos de erosión hídrica y por tanto, de la pérdida de suelo que sufren nuestros campos.Este trabajo presenta un estudio de la relación existente entre la energía cinética de la lluvia y la intensidad de la lluvia con datos de lluvia registrados de forma automática con el disdrómetro de JOSS y WALDVOGEL (1967). Se ha obtenido una ecuación de regresión entre la energía cinética y la intensi- dad pero con algunas peculiaridades, ya que es una ecuación exponencial que introduce un valor máximo de energía cinética. Esta ecuación exponencial ha proporcionado resultados muy satisfacto- nos, presenta errores relativos bajos y además tiene un significado práctico y físico.Además se ha realizado un análisis comparativo entre los resultados obtenidos con esta ecuación exponencial y ecuaciones empíricas tradicionales en el cálculo de la energía cinética. El criterio de selección de las ecuaciones a comparar fue en primer lugar escoger aquellas ecuaciones que fuesen habituales en los estudios de erosión y también aquellas ecuaciones que por la localización donde fueron obtenidas pudiesen ser interesante para su comparación.Los resultados obtenidos de la comparación conducen a pensar que la relación entre la energía cinética y la intensidad no estaría influenciado por condiciones de localización. Se concluye tambien que el uso de la ecuación exponencial obtenida en este estudio resolvería algunos de los problemas que conlleva el uso de algunas de la ecuaciones empíricas.

Aplicación de la cinética de Arrhenius en la piel del aguacate cv. "Hass" en la maduración durante la atmósfera modificada.

Se presenta el estudio de la pérdida de color verde de la piel de aguacate cv. "Hass en atmósferas modificadas con y sin absorción de etileno, a las temperaturas de +20, +6 y +3°C. Los análisis se basan en el modelo de Arrhenius para la evolución de la coordenada "a" de Hunter evaluando el efecto de la conservación a una temperatura por debajo de la crítica para este cultivar. Los resultados obtenidos en forma global, permiten concluir que la evolución del color de la piel de estos frutos se correlaciona con las concentraciones de etileno, en las atmósferas modificadas y con temperatura según el modelo clásico de cinética. Asimismo los cambios en la linealidad de Arrhenius pueden ser muy útiles en la indicación de las temperaturas que interfieren en el desarrollo normal de los frutos, y en especial para determinar la temperatura más próxima a 1 a crítica, sin causar cambios aparentes en la calidad comercial de los frutos.

Electromagnetic Instability Studies in Fusion Plasmas Edge

Este trabajo esta dedicado al estudio de las estructuras macroscópicas conocidas en la literatura como filamentos o blobs que han sido observadas de manera universal en el borde de todo tipo de dispositivos de fusión por confinamiento magnético. Estos filamentos, celdas convectivas elongadas a lo largo de las líneas de campo que surgen en el plasma fuertemente turbulento que existe en este tipo de dispositivos, parecen dominar el transporte radial de partículas y energía en la región conocida como Scrape-off Layer, en la que las líneas de campo dejan de estar cerradas y el plasma es dirigido hacia la pared sólida que forma la cámara de vacío. Aunque el comportamiento y las leyes de escala de estas estructuras son relativamente bien conocidos, no existe aún una teoría generalmente aceptada acerca del mecanismo físico responsable de su formación, que constituye una de las principales incógnitas de la teoría de transporte del borde en plasmas de fusión y una cuestión de gran importancia práctica en el desarrollo de la siguiente generación de reactores de fusión (incluyendo dispositivos como ITER y DEMO), puesto que la eficiencia del confinamiento y la cantidad de energía depositadas en la pared dependen directamente de las características del transporte en el borde. El trabajo ha sido realizado desde una perspectiva eminentemente experimental, incluyendo la observación y el análisis de este tipo de estructuras en el stellarator tipo heliotrón LHD (un dispositivo de gran tamaño, capaz de generar plasmas de características cercanas a las necesarias en un reactor de fusión) y en el stellarator tipo heliac TJ-II (un dispositivo de medio tamaño, capaz de generar plasmas relativamente más fríos pero con una accesibilidad y disponibilidad de diagnósticos mayor). En particular, en LHD se observó la generación de filamentos durante las descargas realizadas en configuración de alta _ (alta presión cinética frente a magnética) mediante una cámara visible ultrarrápida, se caracterizó su comportamiento y se investigó, mediante el análisis estadístico y la comparación con modelos teóricos, el posible papel de la Criticalidad Autoorganizada en la formación de este tipo de estructuras. En TJ-II se diseñó y construyó una cabeza de sonda capaz de medir simultáneamente las fluctuaciones electrostáticas y electromagnéticas del plasma. Gracias a este nuevo diagnóstico se pudieron realizar experimentos con el fin de determinar la presencia de corriente paralela a través de los filamentos (un parámetro de gran importancia en su modelización) y relacionar los dos tipos de fluctuaciones por primera vez en un stellarator. Así mismo, también por primera vez en este tipo de dispositivo, fue posible realizar mediciones simultáneas de los tensores viscoso y magnético (Reynolds y Maxwell) de transporte de cantidad de movimiento. ABSTRACT This work has been devoted to the study of the macroscopic structures known in the literature as filaments or blobs, which have been observed universally in the edge of all kind of magnetic confinement fusion devices. These filaments, convective cells stretching along the magnetic field lines, arise from the highly turbulent plasma present in this kind of machines and seem to dominate radial transport of particles and energy in the region known as Scrapeoff Layer, in which field lines become open and plasma is directed towards the solid wall of the vacuum vessel. Although the behavior and scale laws of these structures are relatively well known, there is no generally accepted theory about the physical mechanism involved in their formation yet, which remains one of the main unsolved questions in the fusion plasmas edge transport theory and a matter of great practical importance for the development of the next generation of fusion reactors (including ITER and DEMO), since efficiency of confinement and the energy deposition levels on the wall are directly dependent of the characteristics of edge transport. This work has been realized mainly from an experimental perspective, including the observation and analysis of this kind of structures in the heliotron stellarator LHD (a large device capable of generating reactor-relevant plasma conditions) and in the heliac stellarator TJ-II (a medium-sized device, capable of relatively colder plasmas, but with greater ease of access and diagnostics availability). In particular, in LHD, the generation of filaments during high _ discharges (with high kinetic to magnetic pressure ratio) was observed by means of an ultrafast visible camera, and the behavior of this structures was characterized. Finally, the potential role of Self-Organized Criticality in the generation of filaments was investigated. In TJ-II, a probe head capable of measuring simultaneously electrostatic and electromagnetic fluctuations in the plasma was designed and built. Thanks to this new diagnostic, experiments were carried out in order to determine the presence of parallel current through filaments (one of the most important parameters in their modelization) and to related electromagnetic (EM) and electrostatic (ES) fluctuations for the first time in an stellarator. As well, also for the first time in this kind of device, measurements of the viscous and magnetic momentum transfer tensors (Reynolds and Maxwell) were performed.

ALIFE: Caracterización Térmica De La Fermentación Del Café

La fermentación del café es una de las operaciones críticas en el proceso de beneficiado por su impacto en la calidad organoléptica del producto final. De forma general la mayor parte de las explotaciones carecen de tecnificación alguna del fermentador, no teniendo el operario ninguna información objetiva sobre la evolución de los parámetros físico-químicos de control del proceso que le ayuden a la toma de decisión. El objetivo de este trabajo es la integración multi-distribuida de sensores de bajo coste y alta fiabilidad para supervisión de la operación de fermentación del café en origen. Concretamente se ha centrado en la caracterización térmica del proceso de fermentación, mediante la instalación multidistribuida de sensores inalámbricos de temperatura (TurboTag ®) en fermentadores de café en Colombia. El análisis: temporal de la evolución de la temperatura a lo largo del proceso de fermentación y complejo mediante diagrama de fases o de cinética de la evolución de las temperaturas, permiten en este trabajo caracterizar térmicamente fermentadores industriales y prototipos plásticos. Se ha constatado la alta heterogeneidad de las fermentaciones tanto intra-lote (amplitudes térmicas de hasta 8.5 ºC y temperaturas en algunos puntos de los fermentadores de 17 ºC) como inter-lote en fermentaciones sujetas a condiciones exógenas y endógenas muy variables

Modelo de ISING en 3-d. aplicación a sensores magnéticos

El modelo de Ising es un problema de física estadística que tiene solución exacta en dos dimensiones, para el caso de tres dimensiones es preciso utilizar procedimientos de simulación. En este trabajo se ha utilizado un método de Monte Carlo para estudiar el comportamiento del sistema en distintas situaciones, siendo de especial interés el estudio del paso por la transición de fase a la temperatura crítica (Temperatura de Curie, Tc). Se ha estudiado la cinética de los dominios magnéticos, considerando la estructura de los dominios desde el punto de vista de la energía, y en consecuencia, hemos tenido en cuenta la energía de canje que tiende a mantener alineados los espines de los electrones en los materiales ferromagnéticos. Este término contribuye a hacer mayor el espesor de la pared, por la tendencia a que los espines de los átomos vecinos se mantengan alineados. Se ha considerado el ferromagnetismo desde el punto de vista cuántico y basado en las propiedades de simetría de las funciones de onda de los electrones, que se manifiestan en variaciones de la energía electrostática de un sistema en función de la orientación de sus espines. Se han estudiado los efectos de histéresis que resultan al aplicar un campo magnético externo a la red y la orientación de los espines de la misma a lo largo de su evolución. Para la determinación de las propiedades de los materiales ferromagnéticos se utiliza el ciclo de histéresis aunque algunas de las propiedades magnéticas, como la dirección de anisotropía, no pueden ser deducidas directamente de esta manera. Se utilizan distintos métodos para la determinación de la anisotropía de las muestras. El acoplamiento entre la magnetización en zonas próximas a la superficie y la magnetización en zonas internas de la muestra puede ser utilizado para obtener un ciclo de histéresis, que permita obtener sensores magnéticos adaptados a las medidas que se quieran realizar. Mediante el control del campo coercitivo y la susceptibilidad se abre una línea de investigación para el desarrollo de sensores magnéticos

Algunos resultados en combustión supersónica

Se discuten algunos aspectos básicos de la combustión supersónica en lo que respecta al uso de motores de aspiración de aire para la propulsión de vehículos espaciales en la atmósfera. Se discuten,asimismo, las ventajas e inconvenientes de la utilización de llamas de difusión en la combustión supersónica. Se incluye, también, un análisis del proceso de mezcla turbulenta supersónica. Se analiza la influencia de la adición de calor debida a la llama de difusión. Y por último, se tratan los efectos de desequilibrio en la cinética química, deteniéndose especialmente en la extinción de la llama,considerando en primer lugar el caso de difusión laminar y a continuación un intento de extensión al caso turbulento.

Modelización de un proceso de biorremediación de suelos contaminados con gasoil

El objeto de esta tesis doctoral es encontrar, mediante herramientas estadísticas, un modelo matemático que prediga la evolución temporal de las concentraciones de los hidrocarburos totales del petróleo (TPH) en un proceso de biorremediación de suelos contaminados con gasoleo. Obtenido el modelo matemático realizar las comparaciones pertinentes para determinar la eficiencia de distintos tratamientos, la influencia tanto del suelo como del nivel de concentración. Para cumplir este objetivo se realizó un diseño de experimentos que tomaba en cuenta dos tipos de suelos, dos niveles de concentración de hidrocarburo y seis tratamientos. Se realizaron en total 1824 ensayos en laboratorio repartidos en 8 campañas de campo durante un año de experiencia tanto en parcelas de 4 x 4 m como en microcosmos. Los resultados obtenidos muestran que el proceso de biorremediación se puede ajustar a un comportamiento exponencial, que es una ecuación de cinética de primer orden, y que las comparaciones realizadas han mostrado resultados satisfactorios de la eficiencia de algunos tratamientos. ABSTRACT The purpose of this thesis is using statistical tools find a mathematical model to predict the time evolution of the total petroleum hydrocarbons (TPH) concentrations in the bioremediation process of diesel contaminated soil. Retrieved mathematical model make relevant comparisons to determine the efficiency of different treatments and the influence of both soil and concentration levels. To achieve this goal a design of experiments was developed, it took into account two types of soil, two concentration levels of hydrocarbon and six treatments. There were a total of 1824 laboratory tests spread over 8 field campaigns during a year of experience in both plots of 4 x 4 m in microcosm. The results obtained show that the bioremediation process can be adjusted to an exponential model, it is an equation of kinetic of first order and that the comparisons have shown satisfactory results of the efficiency of some treatments.

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.

Confined crystallization of nanolayered poly(ethylene terephthalate) using X-ray diffraction methods

The development of crystalline lamellae in ultra-thin layers of poly(ethylene terephthalate) PET confined between polycarbonate (PC) layers in an alternating assembly is investigated as a function of layer thickness by means of X-ray diffraction methods. Isothermal crystallization from the glassy state is in-situ followed by means of small-angle X-ray diffraction. It is found that the reduced size of the PET layers influences the lamellar nanostructure and induces a preferential lamellar orientation. Two lamellar populations, flat-on and edge-on, are found to coexist in a wide range of crystallization temperatures (Tc = 117–150 °C) and within layer thicknesses down to 35 nm. Flat-on lamellae appear at a reduced crystallization rate with respect to bulk PET giving rise to crystals of similar dimensions separated by larger amorphous regions. In addition, a narrower distribution of lamellar orientations develops when the layer thickness is reduced or the crystallization temperature is raised. In case of edge-on lamellae, crystallization conditions also influence the development of lamellar orientation; however, the latter is little affected by the reduced size of the layers. Results suggest that flat-on lamellae arise as a consequence of spatial confinement and edge-on lamellae could be generated due to the interactions with the PC interface. En este trabajo se investiga mediante difracción de rayos X a ángulos bajos (SAXS) y a ángulos altos (WAXS), la cristalización de láminas delgadas de Polietilén tereftalato (PET) confinadas entre láminas de Policarbonato (PC), tomando como referencia PET sin confinar. El espesor de las capas de PET varía entre 35nm y 115 nm. Se realizaron medidas de difracción a tres temperaturas de cristalización (117ºC, 132ºC y 150ºC) encontrándose que el reducido espesor de las capas de PET influye en la estructura lamelar que se desarrolla, induciendo una orientación preferente de las láminas. Se integró la intensidad difractada alrededor del máximo en SAXS para obtener una representación de la intensidad en función del ángulo acimutal. Mediante análisis de mínimos cuadrados se separó la curva experimental obtenida en tres contribuciones diferentes: una función Gausiana que describe la distribución de las orientaciones de las lamelas, una función lorenziana asociada a los máximos meridionales (asociados a las interfases PET-PC) y un background constante. Por otra parte la cantidad de material cristalizado se estimó asumiendo que la intensidad del background en el barrido acimutal, una vez restado el background del primer difractograma (sin máximos en SAXS) se asocia con la contribución del material isotrópico que resta en la muestra cristalizada. Se observa la coexistencia de dos poblaciones de lamelas: flat-on y edge-on. A medida que el espesor de las láminas de PET disminuye la población de las lamelas flat-on experimenta los siguientes cambios: 1) la distribución de orientación se estrecha, 2) la fracción de material cristalizado orientado aumenta, 3) la cinética de cristalización se ralentiza y 4) el largo espaciado aumenta es decir las regiones amorfas entre lamelas aumentan su tamaño. Parece demostrarse que es en las primeras etapas del crecimiento lamelar cuando la restricción espacial fuerza a las lamelas a esta orientación tipo flat-on frente a la orientación edge-on.

Fast measurement of cooling and heating rates in standard cryopreservation protocols based on vitrification

Estudio de la cinética de la crioconservación de tejidos vegetales

Core physics and safety analysis of Generation-IV Sodium Fast Reactors using existing and newly developed computational tools

El futuro de la energía nuclear de fisión dependerá, entre otros factores, de la capacidad que las nuevas tecnologías demuestren para solventar los principales retos a largo plazo que se plantean. Los principales retos se pueden resumir en los siguientes aspectos: la capacidad de proporcionar una solución final, segura y fiable a los residuos radiactivos; así como dar solución a la limitación de recursos naturales necesarios para alimentar los reactores nucleares; y por último, una mejora robusta en la seguridad de las centrales que en definitiva evite cualquier daño potencial tanto en la población como en el medio ambiente como consecuencia de cualquier escenario imaginable o más allá de lo imaginable. Siguiendo estas motivaciones, la Generación IV de reactores nucleares surge con el compromiso de proporcionar electricidad de forma sostenible, segura, económica y evitando la proliferación de material fisible. Entre los sistemas conceptuales que se consideran para la Gen IV, los reactores rápidos destacan por su capacidad potencial de transmutar actínidos a la vez que permiten una utilización óptima de los recursos naturales. Entre los refrigerantes que se plantean, el sodio parece una de las soluciones más prometedoras. Como consecuencia, esta tesis surgió dentro del marco del proyecto europeo CP-ESFR con el principal objetivo de evaluar la física de núcleo y seguridad de los reactores rápidos refrigerados por sodio, al tiempo que se desarrollaron herramientas apropiadas para dichos análisis. Efectivamente, en una primera parte de la tesis, se abarca el estudio de la física del núcleo de un reactor rápido representativo, incluyendo el análisis detallado de la capacidad de transmutar actínidos minoritarios. Como resultado de dichos análisis, se publicó un artículo en la revista Annals of Nuclear Energy [96]. Por otra parte, a través de un análisis de un hipotético escenario nuclear español, se evalúo la disponibilidad de recursos naturales necesarios en el caso particular de España para alimentar una flota específica de reactores rápidos, siguiendo varios escenarios de demanda, y teniendo en cuenta la capacidad de reproducción de plutonio que tienen estos sistemas. Como resultado de este trabajo también surgió una publicación en otra revista científica de prestigio internacional como es Energy Conversion and Management [97]. Con objeto de realizar esos y otros análisis, se desarrollaron diversos modelos del núcleo del ESFR siguiendo varias configuraciones, y para diferentes códigos. Por otro lado, con objeto de poder realizar análisis de seguridad de reactores rápidos, son necesarias herramientas multidimensionales de alta fidelidad específicas para reactores rápidos. Dichas herramientas deben integrar fenómenos relacionados con la neutrónica y con la termo-hidráulica, entre otros, mediante una aproximación multi-física. Siguiendo este objetivo, se evalúo el código de difusión neutrónica ANDES para su aplicación a reactores rápidos. ANDES es un código de resolución nodal que se encuentra implementado dentro del sistema COBAYA3 y está basado en el método ACMFD. Por lo tanto, el método ACMFD fue sometido a una revisión en profundidad para evaluar su aptitud para la aplicación a reactores rápidos. Durante ese proceso, se identificaron determinadas limitaciones que se discutirán a lo largo de este trabajo, junto con los desarrollos que se han elaborado e implementado para la resolución de dichas dificultades. Por otra parte, se desarrolló satisfactoriamente el acomplamiento del código neutrónico ANDES con un código termo-hidráulico de subcanales llamado SUBCHANFLOW, desarrollado recientemente en el KIT. Como conclusión de esta parte, todos los desarrollos implementados son evaluados y verificados. En paralelo con esos desarrollos, se calcularon para el núcleo del ESFR las secciones eficaces en multigrupos homogeneizadas a nivel nodal, así como otros parámetros neutrónicos, mediante los códigos ERANOS, primero, y SERPENT, después. Dichos parámetros se utilizaron más adelante para realizar cálculos estacionarios con ANDES. Además, como consecuencia de la contribución de la UPM al paquete de seguridad del proyecto CP-ESFR, se calcularon mediante el código SERPENT los parámetros de cinética puntual que se necesitan introducir en los típicos códigos termo-hidráulicos de planta, para estudios de seguridad. En concreto, dichos parámetros sirvieron para el análisis del impacto que tienen los actínidos minoritarios en el comportamiento de transitorios. Concluyendo, la tesis presenta una aproximación sistemática y multidisciplinar aplicada al análisis de seguridad y comportamiento neutrónico de los reactores rápidos de sodio de la Gen-IV, usando herramientas de cálculo existentes y recién desarrolladas ad' hoc para tal aplicación. Se ha empleado una cantidad importante de tiempo en identificar limitaciones de los métodos nodales analíticos en su aplicación en multigrupos a reactores rápidos, y se proponen interesantes soluciones para abordarlas. ABSTRACT The future of nuclear reactors will depend, among other aspects, on the capability to solve the long-term challenges linked to this technology. These are the capability to provide a definite, safe and reliable solution to the nuclear wastes; the limitation of natural resources, needed to fuel the reactors; and last but not least, the improved safety, which would avoid any potential damage on the public and or environment as a consequence of any imaginable and beyond imaginable circumstance. Following these motivations, the IV Generation of nuclear reactors arises, with the aim to provide sustainable, safe, economic and proliferationresistant electricity. Among the systems considered for the Gen IV, fast reactors have a representative role thanks to their potential capacity to transmute actinides together with the optimal usage of natural resources, being the sodium fast reactors the most promising concept. As a consequence, this thesis was born in the framework of the CP-ESFR project with the generic aim of evaluating the core physics and safety of sodium fast reactors, as well as the development of the approppriated tools to perform such analyses. Indeed, in a first part of this thesis work, the main core physics of the representative sodium fast reactor are assessed, including a detailed analysis of the capability to transmute minor actinides. A part of the results obtained have been published in Annals of Nuclear Energy [96]. Moreover, by means of the analysis of a hypothetical Spanish nuclear scenario, the availability of natural resources required to deploy an specific fleet of fast reactor is assessed, taking into account the breeding properties of such systems. This work also led to a publication in Energy Conversion and Management [97]. In order to perform those and other analyses, several models of the ESFR core were created for different codes. On the other hand, in order to perform safety studies of sodium fast reactors, high fidelity multidimensional analysis tools for sodium fast reactors are required. Such tools should integrate neutronic and thermal-hydraulic phenomena in a multi-physics approach. Following this motivation, the neutron diffusion code ANDES is assessed for sodium fast reactor applications. ANDES is the nodal solver implemented inside the multigroup pin-by-pin diffusion COBAYA3 code, and is based on the analytical method ACMFD. Thus, the ACMFD was verified for SFR applications and while doing so, some limitations were encountered, which are discussed through this work. In order to solve those, some new developments are proposed and implemented in ANDES. Moreover, the code was satisfactorily coupled with the thermal-hydraulic code SUBCHANFLOW, recently developed at KIT. Finally, the different implementations are verified. In addition to those developments, the node homogenized multigroup cross sections and other neutron parameters were obtained for the ESFR core using ERANOS and SERPENT codes, and employed afterwards by ANDES to perform steady state calculations. Moreover, as a result of the UPM contribution to the safety package of the CP-ESFR project, the point kinetic parameters required by the typical plant thermal-hydraulic codes were computed for the ESFR core using SERPENT, which final aim was the assessment of the impact of minor actinides in transient behaviour. All in all, the thesis provides a systematic and multi-purpose approach applied to the assessment of safety and performance parameters of Generation-IV SFR, using existing and newly developed analytical tools. An important amount of time was employed in identifying the limitations that the analytical nodal diffusion methods present when applied to fast reactors following a multigroup approach, and interesting solutions are proposed in order to overcome them.

Desarrollo y caracterización de un material celular de alta porosidad con base cementicia activada mediante agentes inorgánicos

Esta tesis ha estudiado los morteros celulares, centrándose en la experimentación con pastas de cemento aireadas (PCA) con polvo de aluminio como agente expansor. El objetivo es el desarrollo de un material cementicio con una baja conductividad térmica que sirva como aislamiento térmico. La naturaleza inorgánica del material lo hace incombustible, en contraste con las espumas poliméricas existentes en el mercado, cuya aplicación en cámaras ventiladas ha sido prohibida por normativas de construcción tanto a nivel nacional como internacional. Las posibles aplicaciones son con proyección neumática o en paneles prefabricados. Se han ensayado dos series de pastas de cemento con polvo de aluminio: - Serie WPC/CAC/CH. Mezcla de referencia con cemento blanco (WPC), cemento de aluminato cálcico (CAC) y cal aérea (CH) en proporción 5:1:4. - Serie OPC/CH. Mezcla de referencia con cemento portland con cenizas volantes (OPC) y cal aérea (CH) en proporción OPC/CH de 4:1 A las mezclas de referencia se le han añadido adiciones de metacaolín (MK) (10 y 20%) o sepiolita (SP) (1 y 2%) para observar el efecto que producen tanto en el mortero fresco como en el mortero endurecido. Se ha estudiado la reología de las pastas en estado fresco, analizando el proceso de expansión de las pastas, registrando los valores de tensión de fluencia, aire ocluido y temperatura durante la expansión. Con los valores obtenidos se ha discutido la influencia de las adiciones utilizadas en la cinética de corrosión del polvo de aluminio que genera la expansión, concluyendo que las adiciones puzolánicas (CV y MK) y la SP reducen mucho el periodo de inducción, lo que provoca poros más grandes y mayor cantidad de aire ocluido. Asimismo se ha analizado la relación entre la tensión de fluencia y el contenido de aire ocluido, deduciendo que a mayor tensión de fluencia en el momento de iniciarse la expansión, menor tamaño de poros y contenido de aire ocluido. Finalmente, se han obtenido las densidades y capacidades de retención de agua de los morteros frescos. Para caracterizar la red porosa de las pastas aireadas endurecidas, se obtuvieron tanto las densidades reales, netas, aparentes y relativas como las porosidades abiertas, cerradas y totales con ensayos hídricos. Finalmente se obtuvieron imágenes de los poros con tomografía axial computerizada para obtener las porosimetrías de las muestras. La caracterización de la red porosa ha servido para terminar de analizar lo observado en la evolución de la expansión del mortero fresco. Se ha analizado la influencia de la red porosa en la conductividad térmica, obtenida con caja caliente, comparándola con la existente en la literatura existente tanto de morteros celulares como de espumas poliméricas. Se concluye que los valores de conductividad térmica conseguida están en el mínimo posible para un material celular de base cementicia. La microestructura se ha estudiado con microscopía electrónica de barrido, difracción de rayos X y ensayos térmicos TG/ATD, observando que los productos de hidratación encontrados coinciden con los que se producen en morteros sin airear. Las imágenes SEM y los resultados de ultrasonidos han servido para analizar la presencia de microfisuras de retracción en las pastas aireadas, observando que en las muestras con adiciones de MK y SP, se reduce la presencia de microfisuras. ABSTRACT This thesis has studied cellular mortars, focusing in testing aerated cement pastes with aluminum powder as expansive agent. The purpose is the development of a cementitious material with low thermal conductivity that can be used as thermal isolation. Inorganic nature of this material makes it non-combustible, in contrast with polymeric foams in market, whose application in ventilated double skin façade systems has been banned by building standards, both domestically and internationally. Possible uses for this material are pneumatically sprayed applications and precast panels. Two series of batches with aluminum powder have been tested: - WPC/CAC/CH series. Reference paste with white portland cement (WPC), calcium aluminate cement (CAC) and lime (CH) with 5:1:4 ratio. - OPC/CH series. Reference paste with portland cement with fly ash (OPC) and lime (CH) with 4:1 ratio. Metakaolin (MK) (10 and 20%) or sepiolite (SP) (1 and 2%) additions were used in reference pastes to characterize the effect in fresh and hardened mortar. Rheology in fresh pastes was studied, expansion process of pastes was analyzed, recording yield stress, entrained air and temperature values during expansion. Recorded values were used to discuss influence of additions on reaction kinetics of aluminum powder corrosion, that produces expansion.. Conclusion is that pozzolanic additions (FA and MK) and SP greatly reduce induction period, producing bigger pores and more entrained air. Relation between yield stress and entrained air has been also analyzed, observing that the bigger yield stress at beginning of expansion, the smaller pores size and the lower entrained air values. Finally density and water retention of fresh mortars were obtained. Pore network in hardened aerated cement pastes was characterized by imbibition methods providing true, bulk and relative density, and providing also open, closed and total porosity. Finally, pore system imaging were obtained with computerized axial tomography to study porosimetry of specimens. Pore network characterization was useful to complete facts analysis observed in expansion of fresh mortars. Influence of pore network in thermal conductivity, checked in hot box, was analyzed comparing with those existing values in cellular mortar and polymeric foams researches. It was concluded that thermal conductivity values achieved are close to minimum possible in a cementitious cellular material. Microstructure was studied with Scanning Electron Microscopy, X-Ray Diffractometry and TG-DTA analysis, observing that hydration phases found, are those produced in non aerated mortar. SEM imaging and ultrasound results were useful to analyze shrinkage microcracks in aerated cement pastes, concluding that microcrack presence in specimens with MK and SP additions were reduced.

Investigación de la relación entre el microfracturamiento producido por el explosivo y la disolución del metal en una mena oxidada de cobre

El efecto de la cantidad de explosivos es investigado con el propósito de conocer cuál es la influencia del explosivo con respecto a la cinética de la lixiviación de un mineral oxidado y conocer cuál es la recuperación del metal que contiene. Si la cantidad de explosivo es más grande nos permite obtener un mayor grado de microfracturamiento que permite una lixiviación más eficiente. Este microfracturamiento se conoce a través de medir la superficie específica del mineral, permitiendo conocer cuál es la acción de la onda explosiva. Se extraen tres muestras de mineral oxidado de cobre y se les lixivia en pruebas de laboratorio de lixiviación columnar. Los resultados de la cinética de la lixiviación y la recuperación del metal fueron positivos, es decir, fueron superiores cuando mayor es la superficie específica o la cantidad de explosivo. ABSTRACT The effect of the quantity of explosives is investigated with the purpose to know which is the influence of the explosive with regard to the kinetic of the leaching of a mineral one oxidized, and which the recovery of the metal that contains is. The greater quantity of explosives permits to obtain a greater one microcracking that permits a lixiviation more efficient. This microcracking knows him through measuring the specific surface of the mineral, permitting to know which the action of the explosive wave is. Three samples of mineral oxidized of Cu are extracted, and tests were performed of Iixiviation columnar. The results of the kinetic of the lixiviation and of the recovery of the metal they were positive, that is to say, they enlarged when enlarged the specific surface or the quantity of explosive key.

Simulación numérica de un problema de contaminación atmosférica

La ecuación en derivadas parciales de advección difusión con reacción química es la base de los modelos de dispersión de contaminantes en la atmósfera, y los diferentes métodos numéricos empleados para su resolución han sido objeto de amplios estudios a lo largo de su desarrollo. En esta Tesis se presenta la implementación de un nuevo método conservativo para la resolución de la parte advectiva de la ecuación en derivadas parciales que modela la dispersión de contaminantes dentro del modelo mesoescalar de transporte químico CHIMERE. Este método está basado en una técnica de volúmenes finitos junto con una interpolación racional. La ventaja de este método es la conservación exacta de la masa transportada debido al empleo de la ley de conservación de masas. Para ello emplea una formulación de flujo basado en el cálculo de la integral ponderada dentro de cada celda definida para la discretización del espacio en el método de volúmenes finitos. Los resultados numéricos obtenidos en las simulaciones realizadas (implementando el modelo conservativo para la advección en el modelo CHIMERE) se han comparado con los datos observados de concentración de contaminantes registrados en la red de estaciones de seguimiento y medición distribuidas por la Península Ibérica. Los datos estadísticos de medición del error, la media normalizada y la media absoluta normalizada del error, presentan valores que están dentro de los rangos propuestos por la EPA para considerar el modelo preciso. Además, se introduce un nuevo método para resolver la parte advectivadifusiva de la ecuación en derivadas parciales que modeliza la dispersión de contaminantes en la atmósfera. Se ha empleado un método de diferencias finitas de alto orden para resolver la parte difusiva de la ecuación de transporte de contaminantes junto con el método racional conservativo para la parte advectiva en una y dos dimensiones. Los resultados obtenidos de la aplicación del método a diferentes situaciones incluyendo casos académicos y reales han sido comparados con la solución analítica de la ecuación de advección-difusión, demostrando que el nuevo método proporciona un resultado preciso para aproximar la solución. Por último, se ha desarrollado un modelo completo que contempla los fenómenos advectivo y difusivo con reacción química, usando los métodos anteriores junto con una técnica de diferenciación regresiva (BDF2). Esta técnica consiste en un método implícito multipaso de diferenciación regresiva de segundo orden, que nos permite resolver los problemas rígidos típicos de la química atmosférica, modelizados a través de sistemas de ecuaciones diferenciales ordinarias. Este método hace uso de la técnica iterativa Gauss- Seidel para obtener la solución de la parte implícita de la fórmula BDF2. El empleo de la técnica de Gauss-Seidel en lugar de otras técnicas comúnmente empleadas, como la iteración por el método de Newton, nos proporciona rapidez de cálculo y bajo consumo de memoria, ideal para obtener modelos operativos para la resolución de la cinética química atmosférica. ABSTRACT Extensive research has been performed to solve the atmospheric chemicaladvection- diffusion equation and different numerical methods have been proposed. This Thesis presents the implementation of an exactly conservative method for the advection equation in the European scale Eulerian chemistry transport model CHIMERE based on a rational interpolation and a finite volume algorithm. The advantage of the method is that the cell-integrated average is predicted via a flux formulation, thus the mass is exactly conserved. Numerical results are compared with a set of observation registered at some monitoring sites in Spain. The mean normalized bias and the mean normalized absolute error present values that are inside the range to consider an accurate model performance. In addition, it has been introduced a new method to solve the advectiondiffusion equation. It is based on a high-order accurate finite difference method to solve de diffusion equation together with a rational interpolation and a finite volume to solve the advection equation in one dimension and two dimensions. Numerical results obtained from solving several problems include academic and real atmospheric problems have been compared with the analytical solution of the advection-diffusion equation, showing that the new method give an efficient algorithm for solving such problems. Finally, a complete model has been developed to solve the atmospheric chemical-advection-diffusion equation, adding the conservative method for the advection equation, the high-order finite difference method for the diffusion equation and a second-order backward differentiation formula (BDF2) to solve the atmospheric chemical kinetics. The BDF2 is an implicit, second order multistep backward differentiation formula used to solve the stiff systems of ordinary differential equations (ODEs) from atmospheric chemistry. The Gauss-Seidel iteration is used for approximately solving the implicitly defined BDF solution, giving a faster tool than the more commonly used iterative modified Newton technique. This method implies low start-up costs and a low memory demand due to the use of Gauss-Seidel iteration.

Análisis cinético y cinemático de la carrera de velocidad 100 metros lisos en sus diferentes fases

El presente estudio, se realizó en las instalaciones del CSD con los velocistas del Centro de Alto Rendimiento, se analizó en él una carrera de velocidad de 100m utilizando instrumentos novedosos como son el láser para la obtención de la velocidad instantánea y plantillas podológicas instrumentadas para la medición de la fuerza vertical de los apoyos. Estos instrumentos nos aportan una información imposible de conseguir hasta la fecha fuera de un contexto de laboratorio permitiendo que podamos analizar con detalle en cada fase de la carrera de velocidad tanto la cinemática como la cinética de la misma. Además, a partir de los resultados que nos aportan estos instrumentos y de su posterior análisis, nos encontramos en disposición de entregar a los entrenadores con los que colaboramos un informe detallado acerca de las características específicas de sus corredores en cada tramo de la carrera que les permitirá hacer los ajustes necesarios en su entrenamiento con los que optimizar el rendimiento de sus atletas.