Neutronics design and radiological aspects of HiPER facilities

En este trabajo se han cubierto diferentes asuntos del diseño neutrónico de los aspectos radiológicos de las dos instalaciones del proyecto HiPER. El proyecto HiPER es un proyecto europeo concebido en el marco del programa ESFRI (European Scientific Facilities Research Infrastructure). Está destinado al desarrollo de la energía de fusión nuclear inercial mediante el uso de láseres y el esquema iluminación directa. Consecuentemente, se trata de una instalación con fines exclusivamente civiles. Se divide en dos fases, correspondientes con dos instalaciones: HiPER Engineering y HiPER Reactor. La instalación HiPER Engineering desarrollará las tecnologías implicadas en la ignición de alta repetición de cápsulas de DT por iluminación directa. El HiPER Reactor será una planta demostradora que produzca electricidad haciendo uso de las tecnologías desarrolladas durante la fase HiPER Engineering. El HiPER Engineering se centrará en las tecnologías relevantes para las igniciones a alta repetición de cápsulas de DT usando la iluminación directa. El principal esfuerzo de desarrollo tecnológico se hará en todos los asuntos directamente relacionados con la ignición: láseres, óptica, inyector, y fabricación masiva de cápsulas entre otros. Se espera una producción de entre 5200 MJ/año y 120000 MJ/año dependiendo del éxito de la instalación. Comparado con la energía esperada en NIF, 1200 MJ/año, se trata de un reto y un paso más allá en la protección radiológica. En este trabajo se ha concebido una instalación preliminar. Se ha evaluado desde el punto de vista de la protección radiológica, siendo las personas y la óptica el objeto de protección de este estudio. Se ha establecido una zonificación durante la operación y durante el mantenimiento de la instalación. Además, se ha llevado a cabo una evaluación de la selección de materiales para la cámara de reacción desde el punto de vista de gestión de residuos radiactivos. El acero T91 se ha seleccionado por, siendo un acero comercial, presentar el mismo comportamiento que el acero de baja activación EUROFER97 al evaluarse como residuo con el nivel de irradiación de HiPER Engineering. Teniendo en cuenta los resultados obtenidos para la instalación preliminar y las modificaciones de la instalación motivadas en otros campos, se ha propuesto una instalación avanzada también en este trabajo. Un análisis más profundo de los aspectos radiológicos, así como una evaluación completa de la gestión de todos los residuos radiactivos generados en la instalación se ha llevado a cabo. La protección radiológica se ha incrementado respecto de la instalación preliminar, y todos los residuos pueden gestionarse en un plazo de 30 sin recurrir al enterramiento de residuos. El HiPER Reactor sera una planta demostradora que produzca electricidad basada en las tecnologías de ignición desarrolladas durante la fase HiPER Engineering. El esfuerzo de desarrollo tecnológico se llevará a cabo en los sistemas relacionados con la generación de electricidad en condiciones económicas: manto reproductor de tritio, ciclos de potencia, vida y mantenimiento de componentes, o sistemas de recuperación de tritio entre otros. En este trabajo la principal contribución a HiPER Reactor está relacionada con el diseño de la cámara de reacción y sus extensiones en la planta. La cámara de reacción es la isla nuclear más importante de la planta, donde la mayoría de las reacciones nucleares tienen lugar. Alberga la primera pared, el manto reproductor de tritio y la vasija de vacío. Todo el trabajo realizado aquí ha pivotado en torno al manto reproductor de tritio y sus interacciones con el resto de componentes de la planta. Tras una revisión profunda de la bibliografía de los diseños recientes de cámaras de reacción con características similares a HiPER Reactor, se ha propuesto y justificado un esquema tecnológico innovador para el manto reproductor de tritio. El material fértil selecconado es el eutéctico 15.7 at.% Litio – 84.3 at.% Plomo, LiPb, evitando el uso de berilio como multiplicador neutrónico mientras se garantiza el ajuste online de la tasa de reproducción de tritio mediante el ajuste en el enriquecimiento en 6Li. Aunque se podría haber elegido Litio purom el LiPb evita problemas relacionados con la reactividad química. El precio a pagar es un reto materializado como inventario radiactivo de Z alto en el lazo de LiPb que debe controlarse. El material estructural seleccionado es el acero de baja activación EUROFER97, que estará en contacto directo con le LiPb fluyendo a alta velocidad. En este esquema tecnológico, el LiPb asegurará la autosuficiente de tritio de la planta mientras el propio LiPb extrae del manto el calor sobre él depositado por los neutrones. Este esquema recibe el nombre de manto de Litio-Plomo auto-refrigerado (SCLL por sus siglas en inglés). Respecto de los conceptos SCLL previos, es destacable que nos e requieren componentes del SiC, puesto que no hay campos magnéticos en la cámara de reacción. Consecuentemente, el manto SCLL propuesto para HiPER presenta riesgo tecnológicos moderados, similares a otros dispositivos de fusión magnética, como el HCLL, e incluso inferiores a los del DCLL, puesto que no se require SiC. Los retos que se deben afrontar son el control del inventario de Z alto así como las tasas de corrosión derivadas de la interacción del LiPb con el EUROFE97. En este trabajo se abordan ambos aspectos, y se presentan los respectivos análisis, junto con otros aspectos neutrónicos y de activación, tales como la protección de la vasija de vacío por parte del material fértil para garantizar la resoldabilidad de por vida en la cara externa de la vasija. También se propone y se estudio un ciclo de potencia de Brayton de Helio para dos configuraciones diferentes de refrigeración del sistema primera pared-manto reproductor. Las principales conclusiones de estos estudios son: i) el inventario de Z alto puede controlarse y es comparable al que se encuentra en dispositivos de fusión similares, ii)la vasija de vacío requiere una mayor protección frente a la radiación neutrónica y iii) las tasas de corrosión son demasiado altas y la temperatura media de salida del LiPb es demasiado baja. Tiendo en cuenta estos resultados juntos con otras consideraciones relacionadas con el mantenimiento de componentes y la viabilidad constructiva, se ha propuesto una evolución de la cámara de reacción. Las evoluciones más destacables son la introducción de un reflector neutrónico de grafito, la modificación de la configuración de la óptica final, la forma y el tamaño de la cámara de vacío y una nueva subdivisión modular del manto. Se ha evaluado desde el punto de vista neutrónico, y su análisis y posterior evolución queda fuera del objeto de este trabajo. Los códigos utilizados en este trabajo son: CATIA para la generación de geometrías 3D complejas MCAM para la traducción de archivos de CATIA a formato de input de MCNP MCNP para el transporte de la radiación (neutrones y gammas) y sus respuestas asociadas ACAB para la evolución del inventario isotópico y sus respuestas asociadas MC2ACAB para acoplar MCNP y ACAB para el cómputo de dosis en parada usando la metodología R2S basada en celda. Moritz para visualizar los reultados de MCNP FLUENT para llevar a cabo cálculos de fluido-dinámica Para llevar a cabo este trabajo, han sido necesarias unas destrezas computacionales. Las más relevantes utilizadas son: generación de geometrás 3D complejas y transmisión a MCNP, diferentes tñecnica de reducción de varianza como importancia por celdas y weight windows basado en malla, metodología Rigorous-two-Steps basada en celdas para el cálculo de dosis en parada y la modificación del código ACAB para el cálculos con múltiples espectros en la misma simulación. Como resumen, la contribución de este trabajo al proyecto HiPER son dos diseños conceptuales de instalación: una para HiPER Engineering y otra para HiPER Reactor. La primera se ha estudio en profundidad desde el punto de vista de protección radiológica y gestión de residuos, mientras que la segunda se ha estudiado desde el punto de vista de operación: seguridad, comportamiento, vida y mantenimiento de componentes y eficiencia del ciclo de potencia.

Proyecto de gestión de residuos sólidos urbanos en la Villa de Ibo

El objetivo de la Gestión de Residuos Sólidos Urbanos (GRSU) es paliar los problemas que representan los residuos sólidos urbanos. Entre ellos destacan2: el volumen y espacio que ocupan, molestias que generan, el riesgo para la salud pública, las cuestiones económicas sobre su gestión y efectos en otros campos, la contaminación del medio ambiente y el consumo de recursos. Se trata por tanto de un problema de gestión, con vistas en el medio y largo plazo, que abarca diversas perspectivas como la tecnológica, la política, la medioambiental o la económica.

Impacto y protección radiológica asociados al transporte de residuos de alta actividad (RAA) al almacén temporal centralizado (ATC) en España

Las cuestiones relacionadas con el transporte de residuos radiactivos de alta actividad (RAA) al previsto almacén temporal centralizado (ATC) en Villar de Cañas (Cuenca) están de actualidad, debido a la movilidad que se espera en un futuro próximo, el compromiso con el medio ambiente, la protección de las personas, así, como la normativa legal reguladora. En esta tesis se ha evaluado el impacto radiológico asociado a este tipo de transportes mediante una nueva herramienta de procesamiento de datos, que puede ser de utilidad y servir como documentación complementaria a la recogida en el marco legal del transporte. Además puede facilitar el análisis desde una perspectiva más científica, para investigadores, responsables públicos y técnicos en general, que pueden utilizar dicha herramienta para simular distintos escenarios de transportes radiactivos basados únicamente en datos de los materiales de entrada y las rutas elegidas. Así, conociendo el nivel de radiación a un metro del transporte y eligiendo una ruta, obtendremos los impactos asociados, tales como las poblaciones afectadas, la dosis recibida por la persona más expuesta, el impacto radiológico global, las dosis a la población en el trayecto y el posible detrimento de su salud. En España se prevé una larga “ruta radiactiva” de más de 2.000 kilómetros, por la que el combustible nuclear gastado se transportará presumiblemente por carretera desde las centrales nucleares hasta el ATC, así como los residuos vitrificados procedentes del reprocesado del combustible de la central nuclear Vandellos I, que en la actualidad están en Francia. Como conclusión más importante, se observa que la emisión de radiaciones ionizantes procedentes del transporte de residuos radiactivos de alta actividad en España, en operación normal, no es significativa a la hora de generar efectos adversos en la salud humana y su impacto radiológico puede considerarse despreciable. En caso de accidente, aunque la posibilidad del suceso es remota, las emisiones, no serán determinantes a la hora de generar efectos adversos en la salud humana. Issues related to the transport of high level radioactive wastes (HLW) to the new centralised temporary storage facility to be built in Villar de Cañas (Cuenca) are attracting renewed attention due to the mobility expected in the near future for these materials, the commitment to the environment, the protection of persons and the legal regulatory standards. This study assesses the radiological impacts associated with this type of transport by means of a new dataprocessing tool, which may be of use and serve as documentation complementary to that included in the legal framework covering transport. Furthermore, it may facilitate analysis from a more scientific perspective for researchers, public servants and technicians in general, who may use the tool to simulate different radioactive transport scenarios based only on input materials data and the routes selected. Thus, by knowing the radiation level at a distance of one metre from the transport and selecting a route, it is possible to obtain the associated impacts, such as the affected populations, the dose received by the most exposed individual, the overall radiological impact and the doses to the public en route and the possible detriment to their health. In Spain a long “radioactive route” of more than 2,000 kilometres is expected, along which spent nuclear fuels will be transported – foreseeably by road – from the nuclear power plants to the CTS facility. The route will also be used for the vitrified wastes from fuel reprocessing of the fuel from Vandellós I nuclear power plant, which are currently in France. In conclusion, it may be observed that the emission of ionising radiations from transport of high level radioactive wastes in Spain is insignificant, in normal operations, as regards the generation of adverse effects for human health, and that the radiological impact may be considered negligible. In the event of an accident, the possibility of which is remote, the emissions will not be also a very determining factor as regards adverse effects for human health.

Ingeniería sísmica : discurso leído en la solemne entrega de los diplomas académicos a los ingenieros industriales de la promoción 134

En estas páginas se ha tomado el tema de la Ingeniería Sísmica como ejemplo y pretexto para mostrar la actitud del ingeniero: observación, comprensión, reflexión y generalización son virtudes comunes en el establecimiento de los paradigmas tanto científicos como técnicos. El ingeniero arranca de ellos, pero debe prolongar su actividad en la elaboración de artefactos que aprovechen las leyes establecidas para extraer un beneficio o conseguir un fin. Por eso el premio Nobel Simon habla de la ingeniería como la ciencia de lo artificial.El éxito o el fracaso del artefacto permite detectar nuevos problemas o establecer un Código de buenas prácticas al que atenerse. Este equilibrio entre práctica profesional e investigación de nuevas posibilidades es el que mantiene activa la llama de esta Escuela y mide la vitalidad de una especialidad. Si falla la creatividad se cae en la rutina y la obsolescencia, si falta la aplicación profesional se cae en el aislamiento. En cualquiera de los casos se pierde la competitividad, el mercado, y en definitiva, las señas de identidad de la profesión.

Fukushima: el plan de estabilización y control de la central

Se describen de forma panorámica los pasos e hitos principales que comprende el plan de estabilización y recuperación de la central de Fukushima-Daiichi, con los logros alcanzados hasta conseguir una refrigeración estable y llegar a la parada fría de los reactores, mantener la refrigeración de las piscinas, detener prácticamente por completo de las descargas radiactivas, gestionar y reutilizar cantidades enormes de agua contaminada, limpiar el emplazamiento y aislar la unidad 1. También se presentan sucintamente las etapas previstas posteriormente hasta llegar al desmantelamiento de la central, pasando por la extracción del combustible de las piscinas, y en una fase posterior mucho más delicada, del material fundido de los núcleos de los reactores y el resto de materiales y residuos altamente radiactivos.

Modelo de gestión logística para pymes industriales

This paper presents the design and results of applying a model for logistics management in industrial SMEs. To identify the variables in the model, we conducted a thorough review of the state of the art logistics management; to characterize SMEs, developed a Likert questionnaire with the variables collected in the previous step. Once validated the questionnaire, was applied the same to a group of seventy-five (75) SMEs in the industrial sector, located in Bolivar State, Venezuela. To determine statistically the most relevant variables of management was used exploratory factor analysis technique applied to the data collected. The qualification obtained for all companies evaluated (47% compliance), highlights the weakness of logistics management in industrial SME.

Gestión docente de la creatividad en la Escuela Técnica Superior de Ingenieros Industriales (ETSII) de Madrid: análisis de los procesos de generación, transmisión y absorción de conocimiento

El análisis de la metodología docente empleada en la ETSII se ha llevado a cabo con encuestas realizadas a profesores y alumnos del Centro. Las conclusiones de las encuestas apuntan a que hay un amplio campo de mejoras encaminadas a mejorar la transmisión de habilidades creativas a los futuros ingenieros. La aplicación de las TICs, potenciación del trabajo en equipo, mayor control de la calidad de la enseñanza impartida, desarrollo de grupos de investigación de las metodologías docentes, reducción del número de alumnos por clase y contratación de profesores con experiencia profesional, son aspectos a potenciar.

Gestión de seguridad y salud en el trabajo: aplicación a las pymes industriales

En este artículo se presentan el diseño y los resultados de la aplicación de un modelo para la evaluación y mejora de la gestión de la seguridad y salud en el trabajo; esta investigación forma parte de un proyecto más amplio dirigido al diseño de un modelo de sistema inteligente para la evaluación y mejora de la gestión empresarial, que soporte la toma de decisión en las Pymes industriales.

Gestión estratégica. Aplicación a las PYMES industriales

En este artículo se presentan el diseño y los resultados de la aplicación de un modelo para la evaluación y mejora de la gestión estratégica; esta investigación forma parte de un proyecto más amplio dirigido al diseño de un modelo de sistema inteligente para la evaluación y mejora de la gestión empresarial, que soporte la toma de decisión en las Pymes industriales. Para la determinación de las variables, elementos y evidencias, se hizo una revisión minuciosa de la bibliografía existente sobre gestión estratégica; para hallar las relaciones entre las variables se conformó un equipo de expertos, el cual a través de la matriz de análisis estructural, determinó los valores de motricidad y dependencia de las variables en estudio. Una vez diseñado el sistema de evaluación, se aplicó el instrumento a un grupo de diecinueve (19) Pymes industriales, ubicadas en Ciudad Guayana, Estado Bolívar, Venezuela, las cuales obtuvieron una calificación de 39 puntos en promedio, sobre una base de medición de 100, lo que evidencia que estas empresas no diseñan, implementan y evalúan estrategias, lo cual les impide llevar a cabo una gestión según un plan estratégico y no a merced de los vaivenes del mercado en el cual se desenvuelven.

Cultivos hortícolas industriales en las alternativas de regadío: Productividades y mano de obra empleada.

El cultivo de especies hortícolas con destino industrial se considera como el elemento que puede mejorar a corto plazo la productividad económica y de la mano de obra de las explotaciones de regadío de las Vegas del Guadiana. Se estudian cinco alternativas que incluyen en total siete cultivos diferentes y para dos niveles de explotación (5 ha y 20 ha de superficie). Se calculan las necesidades y distribución anuales de la mano de obra, las productividades y la incidencia de la mecanización de la recolección en las especies hortícolas en que ello es hoy posible. Las diferencias en mano de obra empleada entre las distintas alternativas estudiadas son muy importantes, así como las productividades. La incidencia de la mecanización de la recolección tiende a favorecer la introducción de cultivos hortícolas industriales en explotaciones de mayor dimensión, hoy dedicadas a la producción de grano.

Selección de alternativas para el tratamiento de depósitos de residuos radiactivos bajo una respectiva multiatributo borrosa.

La planta química Pridneprovsky fue una de las empresas mas grandes de procesamiento de uranio de la antigua URSS. Los depósitos de Zapadnoe contienen la mayoría de estos residuos. Proponemos un marco teórico basado en el Análisis de Decisión Multicriterio y la Lógica Borrosa para analizar diferentes alternativas de remediación de los residuos Zapadnoe, en las que objetivos económicos, sociales y ambientales potencialmente conflictivos entre sí son simultáneamente tenidos en cuenta. Se ha construido una jerarquía de objetivos, que incluye todos los aspectos relevantes del problema, en la que se observa claramente los objetivos y atributos que son necesarios conseguir si se quieren alcanzar objetivos de un nivel superior. Se han usado conjuntos borrosos en lugar de valores exactos, para evaluar las alternativas de remediación en los diferentes criterios y cuanticar las preferencias del decisor. Para ello, se propone que las alternativas de remediación deben ser evaluadas por medio de una función de utilidad de multiatributo aditiva. Este método producirá una utilidad esperada para cada alternativa representada por medio de un numero borroso trapezoidal. Finalmente se usará un algoritmo de ordenamiento de números borrosos para elegir la mejor alternativa.

Evolución del flujo de residuos generados en obras de edificación residencial

En la actualidad, conocer los residuos de construcción y demolición (RCD) generados en una obra es indispensable para optimizar su gestión. Por tanto, cualquier herramienta que establezca una estimación de los residuos generados debe ser considerada como un instrumento más para plantear soluciones reales en el campo de la sostenibilidad. La presente comunicación da a conocer la evolución del flujo de residuos generados en la construcción de edificios residenciales, identificando la fase de obra que más RCD produce. Para ello, se han analizado varias obras reales construidas en la Comunidad de Madrid. Los resultados obtenidos permiten conocer, con la antelación suficiente, la cantidad de RCD que se va a generar, y por tanto planificar la cantidad y dimensión de los contenedores necesarios en un determinado momento de la obra, así como el espacio preciso en la obra para su correcta gestión.

Simulación en tiempo real de vehículos industriales con modelos multicuerpo de gran complejidad

El objetivo de esta Tesis ha sido la consecución de simulaciones en tiempo real de vehículos industriales modelizados como sistemas multicuerpo complejos formados por sólidos rígidos. Para el desarrollo de un programa de simulación deben considerarse cuatro aspectos fundamentales: la modelización del sistema multicuerpo (tipos de coordenadas, pares ideales o impuestos mediante fuerzas), la formulación a utilizar para plantear las ecuaciones diferenciales del movimiento (coordenadas dependientes o independientes, métodos globales o topológicos, forma de imponer las ecuaciones de restricción), el método de integración numérica para resolver estas ecuaciones en el tiempo (integradores explícitos o implícitos) y finalmente los detalles de la implementación realizada (lenguaje de programación, librerías matemáticas, técnicas de paralelización). Estas cuatro etapas están interrelacionadas entre sí y todas han formado parte de este trabajo. Desde la generación de modelos de una furgoneta y de camión con semirremolque, el uso de tres formulaciones dinámicas diferentes, la integración de las ecuaciones diferenciales del movimiento mediante métodos explícitos e implícitos, hasta el uso de funciones BLAS, de técnicas de matrices sparse y la introducción de paralelización para utilizar los distintos núcleos del procesador. El trabajo presentado en esta Tesis ha sido organizado en 8 capítulos, dedicándose el primero de ellos a la Introducción. En el Capítulo 2 se presentan dos formulaciones semirrecursivas diferentes, de las cuales la primera está basada en una doble transformación de velocidades, obteniéndose las ecuaciones diferenciales del movimiento en función de las aceleraciones relativas independientes. La integración numérica de estas ecuaciones se ha realizado con el método de Runge-Kutta explícito de cuarto orden. La segunda formulación está basada en coordenadas relativas dependientes, imponiendo las restricciones por medio de penalizadores en posición y corrigiendo las velocidades y aceleraciones mediante métodos de proyección. En este segundo caso la integración de las ecuaciones del movimiento se ha llevado a cabo mediante el integrador implícito HHT (Hilber, Hughes and Taylor), perteneciente a la familia de integradores estructurales de Newmark. En el Capítulo 3 se introduce la tercera formulación utilizada en esta Tesis. En este caso las uniones entre los sólidos del sistema se ha realizado mediante uniones flexibles, lo que obliga a imponer los pares por medio de fuerzas. Este tipo de uniones impide trabajar con coordenadas relativas, por lo que la posición del sistema y el planteamiento de las ecuaciones del movimiento se ha realizado utilizando coordenadas Cartesianas y parámetros de Euler. En esta formulación global se introducen las restricciones mediante fuerzas (con un planteamiento similar al de los penalizadores) y la estabilización del proceso de integración numérica se realiza también mediante proyecciones de velocidades y aceleraciones. En el Capítulo 4 se presenta una revisión de las principales herramientas y estrategias utilizadas para aumentar la eficiencia de las implementaciones de los distintos algoritmos. En primer lugar se incluye una serie de consideraciones básicas para aumentar la eficiencia numérica de las implementaciones. A continuación se mencionan las principales características de los analizadores de códigos utilizados y también las librerías matemáticas utilizadas para resolver los problemas de álgebra lineal tanto con matrices densas como sparse. Por último se desarrolla con un cierto detalle el tema de la paralelización en los actuales procesadores de varios núcleos, describiendo para ello el patrón empleado y las características más importantes de las dos herramientas propuestas, OpenMP y las TBB de Intel. Hay que señalar que las características de los sistemas multicuerpo problemas de pequeño tamaño, frecuente uso de la recursividad, y repetición intensiva en el tiempo de los cálculos con fuerte dependencia de los resultados anteriores dificultan extraordinariamente el uso de técnicas de paralelización frente a otras áreas de la mecánica computacional, tales como por ejemplo el cálculo por elementos finitos. Basándose en los conceptos mencionados en el Capítulo 4, el Capítulo 5 está dividido en tres secciones, una para cada formulación propuesta en esta Tesis. En cada una de estas secciones se describen los detalles de cómo se han realizado las distintas implementaciones propuestas para cada algoritmo y qué herramientas se han utilizado para ello. En la primera sección se muestra el uso de librerías numéricas para matrices densas y sparse en la formulación topológica semirrecursiva basada en la doble transformación de velocidades. En la segunda se describe la utilización de paralelización mediante OpenMP y TBB en la formulación semirrecursiva con penalizadores y proyecciones. Por último, se describe el uso de técnicas de matrices sparse y paralelización en la formulación global con uniones flexibles y parámetros de Euler. El Capítulo 6 describe los resultados alcanzados mediante las formulaciones e implementaciones descritas previamente. Este capítulo comienza con una descripción de la modelización y topología de los dos vehículos estudiados. El primer modelo es un vehículo de dos ejes del tipo chasis-cabina o furgoneta, perteneciente a la gama de vehículos de carga medianos. El segundo es un vehículo de cinco ejes que responde al modelo de un camión o cabina con semirremolque, perteneciente a la categoría de vehículos industriales pesados. En este capítulo además se realiza un estudio comparativo entre las simulaciones de estos vehículos con cada una de las formulaciones utilizadas y se presentan de modo cuantitativo los efectos de las mejoras alcanzadas con las distintas estrategias propuestas en esta Tesis. Con objeto de extraer conclusiones más fácilmente y para evaluar de un modo más objetivo las mejoras introducidas en la Tesis, todos los resultados de este capítulo se han obtenido con el mismo computador, que era el top de la gama Intel Xeon en 2007, pero que hoy día está ya algo obsoleto. Por último los Capítulos 7 y 8 están dedicados a las conclusiones finales y las futuras líneas de investigación que pueden derivar del trabajo realizado en esta Tesis. Los objetivos de realizar simulaciones en tiempo real de vehículos industriales de gran complejidad han sido alcanzados con varias de las formulaciones e implementaciones desarrolladas. ABSTRACT The objective of this Dissertation has been the achievement of real time simulations of industrial vehicles modeled as complex multibody systems made up by rigid bodies. For the development of a simulation program, four main aspects must be considered: the modeling of the multibody system (types of coordinates, ideal joints or imposed by means of forces), the formulation to be used to set the differential equations of motion (dependent or independent coordinates, global or topological methods, ways to impose constraints equations), the method of numerical integration to solve these equations in time (explicit or implicit integrators) and the details of the implementation carried out (programming language, mathematical libraries, parallelization techniques). These four stages are interrelated and all of them are part of this work. They involve the generation of models for a van and a semitrailer truck, the use of three different dynamic formulations, the integration of differential equations of motion through explicit and implicit methods, the use of BLAS functions and sparse matrix techniques, and the introduction of parallelization to use the different processor cores. The work presented in this Dissertation has been structured in eight chapters, the first of them being the Introduction. In Chapter 2, two different semi-recursive formulations are shown, of which the first one is based on a double velocity transformation, thus getting the differential equations of motion as a function of the independent relative accelerations. The numerical integration of these equations has been made with the Runge-Kutta explicit method of fourth order. The second formulation is based on dependent relative coordinates, imposing the constraints by means of position penalty coefficients and correcting the velocities and accelerations by projection methods. In this second case, the integration of the motion equations has been carried out by means of the HHT implicit integrator (Hilber, Hughes and Taylor), which belongs to the Newmark structural integrators family. In Chapter 3, the third formulation used in this Dissertation is presented. In this case, the joints between the bodies of the system have been considered as flexible joints, with forces used to impose the joint conditions. This kind of union hinders to work with relative coordinates, so the position of the system bodies and the setting of the equations of motion have been carried out using Cartesian coordinates and Euler parameters. In this global formulation, constraints are introduced through forces (with a similar approach to the penalty coefficients) are presented. The stabilization of the numerical integration is carried out also by velocity and accelerations projections. In Chapter 4, a revision of the main computer tools and strategies used to increase the efficiency of the implementations of the algorithms is presented. First of all, some basic considerations to increase the numerical efficiency of the implementations are included. Then the main characteristics of the code’ analyzers used and also the mathematical libraries used to solve linear algebra problems (both with dense and sparse matrices) are mentioned. Finally, the topic of parallelization in current multicore processors is developed thoroughly. For that, the pattern used and the most important characteristics of the tools proposed, OpenMP and Intel TBB, are described. It needs to be highlighted that the characteristics of multibody systems small size problems, frequent recursion use and intensive repetition along the time of the calculation with high dependencies of the previous results complicate extraordinarily the use of parallelization techniques against other computational mechanics areas, as the finite elements computation. Based on the concepts mentioned in Chapter 4, Chapter 5 is divided into three sections, one for each formulation proposed in this Dissertation. In each one of these sections, the details of how these different proposed implementations have been made for each algorithm and which tools have been used are described. In the first section, it is shown the use of numerical libraries for dense and sparse matrices in the semirecursive topological formulation based in the double velocity transformation. In the second one, the use of parallelization by means OpenMP and TBB is depicted in the semi-recursive formulation with penalization and projections. Lastly, the use of sparse matrices and parallelization techniques is described in the global formulation with flexible joints and Euler parameters. Chapter 6 depicts the achieved results through the formulations and implementations previously described. This chapter starts with a description of the modeling and topology of the two vehicles studied. The first model is a two-axle chassis-cabin or van like vehicle, which belongs to the range of medium charge vehicles. The second one is a five-axle vehicle belonging to the truck or cabin semi-trailer model, belonging to the heavy industrial vehicles category. In this chapter, a comparative study is done between the simulations of these vehicles with each one of the formulations used and the improvements achieved are presented in a quantitative way with the different strategies proposed in this Dissertation. With the aim of deducing the conclusions more easily and to evaluate in a more objective way the improvements introduced in the Dissertation, all the results of this chapter have been obtained with the same computer, which was the top one among the Intel Xeon range in 2007, but which is rather obsolete today. Finally, Chapters 7 and 8 are dedicated to the final conclusions and the future research projects that can be derived from the work presented in this Dissertation. The objectives of doing real time simulations in high complex industrial vehicles have been achieved with the formulations and implementations developed.

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.

Criterios de sostenibilidad en la localización de instalaciones industriales

La incorporación de nuevos valores del comportamiento empresarial condiciona los criterios para adoptar decisiones estratégicas ya que coexisten los clásicos principios de inversión y beneficio con las obligaciones de respeto a la calidad medioambiental y los criterios de comportamiento sostenible. Todo ello justifica una investigación orientada a conseguir el objetivo de proporcionar una herramienta de ayuda en la toma de decisiones en materia de localización, ampliación, traslado o cierre para jerarquizar objetiva y eficazmente las alternativas evaluadas. Se ha decidido evaluar la incorporación de factores de sostenibilidad cuyo empleo esté apoyado en métodos de análisis multicriterio. Así de la intersección entre la primera etapa de elección de criterios de carácter eminentemente técnicos con la segunda donde se establecen las preferencias reales del centro decisor mediante la denominada función criterio o de utilidad (Romero, 1996) se puede obtener una solución óptima y sostenible al problema. 1-27 La investigación se ha centrado en la etapa de definición de criterios con el objetivo adicional de que buena parte de ellos puedan ser criterios utilizados frecuentemente en los Sistemas de Gestión Medioambiental y en la evaluación de comportamientos sobre sostenibilidad y Responsabilidad Social Corporativa. A su vez, se ha previsto alcanzar cinco objetivos específicos: 1. Evaluar el concepto de desarrollo sostenible para su adaptación a los múltiples ámbitos en los que actualmente se aplica. 2. Revisar la aplicación del concepto de sostenibilidad al ámbito empresarial, incidiendo en los instrumentos y metodologías desarrolladas para ello. 3. Investigar y seleccionar los aspectos frecuentemente considerados en la evaluación de la sostenibilidad aplicada al ámbito empresarial y explorar su utilidad como factor de sostenibilidad en la toma de decisiones estratégicas. 4. Establecer una metodología para la selección de indicadores de sostenibilidad y un modelo que permita su aplicación en la toma de decisiones estratégicas. 5. Evaluar y validar el comportamiento del modelo mediante su aplicación a situaciones reales. Los tres primeros objetivos tienen un carácter exclusivamente teórico orientado a lograr una adecuada identificación de los planteamientos realizados en torno al concepto general de sostenibilidad y las particularidades de su aplicación en el ámbito empresarial. El cumplimiento del cuarto y quinto objetivos supone evaluar la utilidad del modelo y analizar las dificultades en su utilización. El logro de los objetivos de la investigación puede contribuir al desarrollo de nuevos instrumentos que permitan incluir la sostenibilidad como criterio en la toma de decisiones relativa a la implantación, ampliación, cierre o traslado de instalaciones industriales. Adicionalmente, el resultado de esta investigación puede ofrecer la utilidad añadida de mejorar o potenciar los instrumentos actuales de evaluación de la responsabilidad social empresarial y de gestión medioambiental.