Monte Carlo Neutronics and Thermal-hydraulics coupling applied to Fast Reactors

Un escenario habitualmente considerado para el uso sostenible y prolongado de la energía nuclear contempla un parque de reactores rápidos refrigerados por metales líquidos (LMFR) dedicados al reciclado de Pu y la transmutación de actínidos minoritarios (MA). Otra opción es combinar dichos reactores con algunos sistemas subcríticos asistidos por acelerador (ADS), exclusivamente destinados a la eliminación de MA. El diseño y licenciamiento de estos reactores innovadores requiere herramientas computacionales prácticas y precisas, que incorporen el conocimiento obtenido en la investigación experimental de nuevas configuraciones de reactores, materiales y sistemas. A pesar de que se han construido y operado un cierto número de reactores rápidos a nivel mundial, la experiencia operacional es todavía reducida y no todos los transitorios se han podido entender completamente. Por tanto, los análisis de seguridad de nuevos LMFR están basados fundamentalmente en métodos deterministas, al contrario que las aproximaciones modernas para reactores de agua ligera (LWR), que se benefician también de los métodos probabilistas. La aproximación más usada en los estudios de seguridad de LMFR es utilizar una variedad de códigos, desarrollados a base de distintas teorías, en busca de soluciones integrales para los transitorios e incluyendo incertidumbres. En este marco, los nuevos códigos para cálculos de mejor estimación ("best estimate") que no incluyen aproximaciones conservadoras, son de una importancia primordial para analizar estacionarios y transitorios en reactores rápidos. Esta tesis se centra en el desarrollo de un código acoplado para realizar análisis realistas en reactores rápidos críticos aplicando el método de Monte Carlo. Hoy en día, dado el mayor potencial de recursos computacionales, los códigos de transporte neutrónico por Monte Carlo se pueden usar de manera práctica para realizar cálculos detallados de núcleos completos, incluso de elevada heterogeneidad material. Además, los códigos de Monte Carlo se toman normalmente como referencia para los códigos deterministas de difusión en multigrupos en aplicaciones con reactores rápidos, porque usan secciones eficaces punto a punto, un modelo geométrico exacto y tienen en cuenta intrínsecamente la dependencia angular de flujo. En esta tesis se presenta una metodología de acoplamiento entre el conocido código MCNP, que calcula la generación de potencia en el reactor, y el código de termohidráulica de subcanal COBRA-IV, que obtiene las distribuciones de temperatura y densidad en el sistema. COBRA-IV es un código apropiado para aplicaciones en reactores rápidos ya que ha sido validado con resultados experimentales en haces de barras con sodio, incluyendo las correlaciones más apropiadas para metales líquidos. En una primera fase de la tesis, ambos códigos se han acoplado en estado estacionario utilizando un método iterativo con intercambio de archivos externos. El principal problema en el acoplamiento neutrónico y termohidráulico en estacionario con códigos de Monte Carlo es la manipulación de las secciones eficaces para tener en cuenta el ensanchamiento Doppler cuando la temperatura del combustible aumenta. Entre todas las opciones disponibles, en esta tesis se ha escogido la aproximación de pseudo materiales, y se ha comprobado que proporciona resultados aceptables en su aplicación con reactores rápidos. Por otro lado, los cambios geométricos originados por grandes gradientes de temperatura en el núcleo de reactores rápidos resultan importantes para la neutrónica como consecuencia del elevado recorrido libre medio del neutrón en estos sistemas. Por tanto, se ha desarrollado un módulo adicional que simula la geometría del reactor en caliente y permite estimar la reactividad debido a la expansión del núcleo en un transitorio. éste módulo calcula automáticamente la longitud del combustible, el radio de la vaina, la separación de los elementos de combustible y el radio de la placa soporte en función de la temperatura. éste efecto es muy relevante en transitorios sin inserción de bancos de parada. También relacionado con los cambios geométricos, se ha implementado una herramienta que, automatiza el movimiento de las barras de control en busca d la criticidad del reactor, o bien calcula el valor de inserción axial las barras de control. Una segunda fase en la plataforma de cálculo que se ha desarrollado es la simulació dinámica. Puesto que MCNP sólo realiza cálculos estacionarios para sistemas críticos o supercríticos, la solución más directa que se propone sin modificar el código fuente de MCNP es usar la aproximación de factorización de flujo, que resuelve por separado la forma del flujo y la amplitud. En este caso se han estudiado en profundidad dos aproximaciones: adiabática y quasiestática. El método adiabático usa un esquema de acoplamiento que alterna en el tiempo los cálculos neutrónicos y termohidráulicos. MCNP calcula el modo fundamental de la distribución de neutrones y la reactividad al final de cada paso de tiempo, y COBRA-IV calcula las propiedades térmicas en el punto intermedio de los pasos de tiempo. La evolución de la amplitud de flujo se calcula resolviendo las ecuaciones de cinética puntual. Este método calcula la reactividad estática en cada paso de tiempo que, en general, difiere de la reactividad dinámica que se obtendría con la distribución de flujo exacta y dependiente de tiempo. No obstante, para entornos no excesivamente alejados de la criticidad ambas reactividades son similares y el método conduce a resultados prácticos aceptables. Siguiendo esta línea, se ha desarrollado después un método mejorado para intentar tener en cuenta el efecto de la fuente de neutrones retardados en la evolución de la forma del flujo durante el transitorio. El esquema consiste en realizar un cálculo cuasiestacionario por cada paso de tiempo con MCNP. La simulación cuasiestacionaria se basa EN la aproximación de fuente constante de neutrones retardados, y consiste en dar un determinado peso o importancia a cada ciclo computacial del cálculo de criticidad con MCNP para la estimación del flujo final. Ambos métodos se han verificado tomando como referencia los resultados del código de difusión COBAYA3 frente a un ejercicio común y suficientemente significativo. Finalmente, con objeto de demostrar la posibilidad de uso práctico del código, se ha simulado un transitorio en el concepto de reactor crítico en fase de diseño MYRRHA/FASTEF, de 100 MW de potencia térmica y refrigerado por plomo-bismuto. ABSTRACT Long term sustainable nuclear energy scenarios envisage a fleet of Liquid Metal Fast Reactors (LMFR) for the Pu recycling and minor actinides (MAs) transmutation or combined with some accelerator driven systems (ADS) just for MAs elimination. Design and licensing of these innovative reactor concepts require accurate computational tools, implementing the knowledge obtained in experimental research for new reactor configurations, materials and associated systems. Although a number of fast reactor systems have already been built, the operational experience is still reduced, especially for lead reactors, and not all the transients are fully understood. The safety analysis approach for LMFR is therefore based only on deterministic methods, different from modern approach for Light Water Reactors (LWR) which also benefit from probabilistic methods. Usually, the approach adopted in LMFR safety assessments is to employ a variety of codes, somewhat different for the each other, to analyze transients looking for a comprehensive solution and including uncertainties. In this frame, new best estimate simulation codes are of prime importance in order to analyze fast reactors steady state and transients. This thesis is focused on the development of a coupled code system for best estimate analysis in fast critical reactor. Currently due to the increase in the computational resources, Monte Carlo methods for neutrons transport can be used for detailed full core calculations. Furthermore, Monte Carlo codes are usually taken as reference for deterministic diffusion multigroups codes in fast reactors applications because they employ point-wise cross sections in an exact geometry model and intrinsically account for directional dependence of the ux. The coupling methodology presented here uses MCNP to calculate the power deposition within the reactor. The subchannel code COBRA-IV calculates the temperature and density distribution within the reactor. COBRA-IV is suitable for fast reactors applications because it has been validated against experimental results in sodium rod bundles. The proper correlations for liquid metal applications have been added to the thermal-hydraulics program. Both codes are coupled at steady state using an iterative method and external files exchange. The main issue in the Monte Carlo/thermal-hydraulics steady state coupling is the cross section handling to take into account Doppler broadening when temperature rises. Among every available options, the pseudo materials approach has been chosen in this thesis. This approach obtains reasonable results in fast reactor applications. Furthermore, geometrical changes caused by large temperature gradients in the core, are of major importance in fast reactor due to the large neutron mean free path. An additional module has therefore been included in order to simulate the reactor geometry in hot state or to estimate the reactivity due to core expansion in a transient. The module automatically calculates the fuel length, cladding radius, fuel assembly pitch and diagrid radius with the temperature. This effect will be crucial in some unprotected transients. Also related to geometrical changes, an automatic control rod movement feature has been implemented in order to achieve a just critical reactor or to calculate control rod worth. A step forward in the coupling platform is the dynamic simulation. Since MCNP performs only steady state calculations for critical systems, the more straight forward option without modifying MCNP source code, is to use the flux factorization approach solving separately the flux shape and amplitude. In this thesis two options have been studied to tackle time dependent neutronic simulations using a Monte Carlo code: adiabatic and quasistatic methods. The adiabatic methods uses a staggered time coupling scheme for the time advance of neutronics and the thermal-hydraulics calculations. MCNP computes the fundamental mode of the neutron flux distribution and the reactivity at the end of each time step and COBRA-IV the thermal properties at half of the the time steps. To calculate the flux amplitude evolution a solver of the point kinetics equations is used. This method calculates the static reactivity in each time step that in general is different from the dynamic reactivity calculated with the exact flux distribution. Nevertheless, for close to critical situations, both reactivities are similar and the method leads to acceptable practical results. In this line, an improved method as an attempt to take into account the effect of delayed neutron source in the transient flux shape evolutions is developed. The scheme performs a quasistationary calculation per time step with MCNP. This quasistationary simulations is based con the constant delayed source approach, taking into account the importance of each criticality cycle in the final flux estimation. Both adiabatic and quasistatic methods have been verified against the diffusion code COBAYA3, using a theoretical kinetic exercise. Finally, a transient in a critical 100 MWth lead-bismuth-eutectic reactor concept is analyzed using the adiabatic method as an application example in a real system.

Three-dimensional distribution of cortical synapses: a replicated point pattern-based analysis

The biggest problem when analyzing the brain is that its synaptic connections are extremely complex. Generally, the billions of neurons making up the brain exchange information through two types of highly specialized structures: chemical synapses (the vast majority) and so-called gap junctions (a substrate of one class of electrical synapse). Here we are interested in exploring the three-dimensional spatial distribution of chemical synapses in the cerebral cortex. Recent research has showed that the three-dimensional spatial distribution of synapses in layer III of the neocortex can be modeled by a random sequential adsorption (RSA) point process, i.e., synapses are distributed in space almost randomly, with the only constraint that they cannot overlap. In this study we hypothesize that RSA processes can also explain the distribution of synapses in all cortical layers. We also investigate whether there are differences in both the synaptic density and spatial distribution of synapses between layers. Using combined focused ion beam milling and scanning electron microscopy (FIB/SEM), we obtained three-dimensional samples from the six layers of the rat somatosensory cortex and identified and reconstructed the synaptic junctions. A total volume of tissue of approximately 4500μm3 and around 4000 synapses from three different animals were analyzed. Different samples, layers and/or animals were aggregated and compared using RSA replicated spatial point processes. The results showed no significant differences in the synaptic distribution across the different rats used in the study. We found that RSA processes described the spatial distribution of synapses in all samples of each layer. We also found that the synaptic distribution in layers II to VI conforms to a common underlying RSA process with different densities per layer. Interestingly, the results showed that synapses in layer I had a slightly different spatial distribution from the other layers.

Selection, training and validation process of a sensory panel for bread analysis: Influence of cultivar on the quality of breads made from common wheat and spelta wheat

A protocol of selection, training and validation of the members of the panel for bread sensory analysis is proposed to assess the influence of wheat cultivar on the sensory quality of bread. Three cultivars of bread wheat and two cultivars of spelt wheat organically-grown under the same edaphoclimatic conditions were milled and baked using the same milling and baking procedure. Through the use of triangle tests, differences were identified between the five breads. Significant differences were found between the spelt breads and those made with bread wheat for the attributes ?crumb cell homogeneity? and ?crumb elasticity?. Significant differences were also found for the odor and flavor attributes, with the bread made with ?Espelta Navarra? being the most complex, from a sensory point of view. Based on the results of this study, we propose that sensory properties should be considered as breeding criteria for future work on genetic improvement.

Fluid-solid interaction in arteries incorporating the autoregulation concept in boundary conditions

In pre-surgery decisions in hospital emergency cases, fast and reliable results of the solid and fluid mechanics problems are of great interest to clinicians. In the current investigation, an iterative process based on a pressure-type boundary condition is proposed in order to reduce the computational costs of blood flow simulations in arteries, without losing control of the important clinical parameters. The incorporation of cardiovascular autoregulation, together with the well-known impedance boundary condition, forms the basis of the proposed methodology. With autoregulation, the instabilities associated with conventional pressure-type or impedance boundary conditions are avoided without an excessive increase in computational costs. The general behaviour of pulsatile blood flow in arteries, which is important from the clinical point of view, is well reproduced through this new methodology. In addition, the interaction between the blood and the arterial walls occurs via a modified weak coupling, which makes the simulation more stable and computationally efficient. Based on in vitro experiments, the hyperelastic behaviour of the wall is characterised and modelled. The applications and benefits of the proposed pressure-type boundary condition are shown in a model of an idealised aortic arch with and without an ascending aorta dissection, which is a common cardiovascular disorder.

Improved modeling of photoluminescent and electroluminescent coupling in multijunction solar cells

The performance of tandem stacks of Group III?V multijunction solar cells continues to improve rapidly, both through improved performance of the individual cells in the stack and throughi ncrease in the number of stacked cells. As the radiative efficiency of these individual cells increases, radiative coupling between the stacked cells becomes an increasingly important factor not only in cell design, but also in accurate efficiency measurement and in determining performance of cells and systems under varying spectral conditions in the field. Past modeling has concentrated on electroluminescent coupling between the cells, although photoluminescent coupling is shown to be important for cells operating near their maximum power point voltage or below or when junction defect recombination is significant. Extension of earlier models i sproposed to allow this non-negligible component of luminescent coupling to be included. Therefined model is validated by measurement of the closely related external emission from both single and double junction cells.

A point mutation in the γ2 subunit of γ-aminobutyric acid type A receptors results in altered benzodiazepine binding site specificity

Benzodiazepines allosterically modulate γ-aminobutyric acid (GABA) evoked chloride currents of γ-aminobutyric acid type A (GABAA) receptors. Coexpression of either rat γ2 or γ3, in combination with α1 and β2 subunits, results both in receptors displaying high [3H]Ro 15-1788 affinity. However, receptors containing a γ3 subunit display a 178-fold reduced affinity to zolpidem as compared with γ2-containing receptors. Eight chimeras between γ2 and γ3 were constructed followed by nine different point mutations in γ2, each to the homologous amino acid residue found in γ3. Chimeric or mutant γ subunits were coexpressed with α1 and β2 in human embryonic kidney 293 cells to localize amino acid residues responsible for the reduced zolpidem affinity. Substitution of a methionine-to-leucine at position 130 of γ2 (γ2M130L) resulted in a 51-fold reduction in zolpidem affinity whereas the affinity to [3H]Ro 15-1788 remained unchanged. The affinity for diazepam was only decreased by about 2-fold. The same mutation resulted in a 9-fold increase in Cl 218872 affinity. A second mutation (γ2M57I) was found to reduce zolpidem affinity by about 4-fold. Wild-type and γ2M130L-containing receptors were functionally expressed in Xenopus oocytes. Upon mutation allosteric coupling between agonist and modulatory sites is preserved. Dose–response curves for zolpidem and for diazepam showed that the zolpidem but not the diazepam apparent affinity is drastically reduced. The apparent GABA affinity is not significantly affected by the γ2M130L mutation. The identified amino acid residues may define part of the benzodiazepine binding pocket of GABAA receptors. As the modulatory site in the GABAA receptor is homologous to the GABA site, and to all agonist sites of related receptors, γ2M130 may either point to a homologous region important for agonist binding in all receptors or define a new region not underlying this principle.

Defect in IgV gene somatic hypermutation in Common Variable Immuno-Deficiency syndrome

Common Variable Immuno-Deficiency (CVID) is the most common symptomatic primary antibody-deficiency syndrome, but the basic immunologic defects underlying this syndrome are not well defined. We report here that among eight patients studied (six CVID and two hypogammaglobulinemic patients with recurrent infections), there is in two CVID patients a dramatic reduction in Ig V gene somatic hypermutation with 40–75% of IgG transcripts totally devoid of mutations in the circulating memory B cell compartment. Functional assays of the T cell compartment point to an intrinsic B cell defect in the process of antibody affinity maturation in these two cases.

Symmetries in bacterial motility

Descriptions are given of three kinds of symmetries encountered in studies of bacterial locomotion, and of the ways in which they are circumvented or broken. A bacterium swims at very low Reynolds number: it cannot propel itself using reciprocal motion (by moving through a sequence of shapes, first forward and then in reverse); cyclic motion is required. A common solution is rotation of a helical filament, either right- or left-handed. The flagellar rotary motor that drives each filament generates the same torque whether spinning clockwise or counterclockwise. This symmetry is broken by coupling to the filament. Finally, bacterial populations, grown in a nutrient medium from an inoculum placed at a single point, usually move outward in symmetric circular rings. Under certain conditions, the cells excrete a chemoattractant, and the rings break up into discrete aggregates that can display remarkable geometric order.

Coupling the phosphotransferase system and the methyl-accepting chemotaxis protein-dependent chemotaxis signaling pathways of Escherichia coli.

Chemotactic responses in Escherichia coli are typically mediated by transmembrane receptors that monitor chemoeffector levels with periplasmic binding domains and communicate with the flagellar motors through two cytoplasmic proteins, CheA and CheY. CheA autophosphorylates and then donates its phosphate to CheY, which in turn controls flagellar rotation. E. coli also exhibits chemotactic responses to substrates that are transported by the phosphoenolpyruvate (PEP)-dependent carbohydrate phosphotransferase system (PTS). Unlike conventional chemoreception, PTS substrates are sensed during their uptake and concomitant phosphorylation by the cell. The phosphoryl groups are transferred from PEP to the carbohydrates through two common intermediates, enzyme I (EI) and phosphohistidine carrier protein (HPr), and then to sugar-specific enzymes II. We found that in mutant strains HPr-like proteins could substitute for HPr in transport but did not mediate chemotactic signaling. In in vitro assays, these proteins exhibited reduced phosphotransfer rates from EI, indicating that the phosphorylation state of EI might link the PTS phospho-relay to the flagellar signaling pathway. Tests with purified proteins revealed that unphosphorylated EI inhibited CheA autophosphorylation, whereas phosphorylated EI did not. These findings suggest the following model for signal transduction in PTS-dependent chemotaxis. During uptake of a PTS carbohydrate, EI is dephosphorylated more rapidly by HPr than it is phosphorylated at the expense of PEP. Consequently, unphosphorylated EI builds up and inhibits CheA autophosphorylation. This slows the flow of phosphates to CheY, eliciting an up-gradient swimming response by the cell.

Numerical modeling of electromagnetic coupling phenomena in the vicinities of overhead power transmission lines.

Electromagnetic coupling phenomena between overhead power transmission lines and other nearby structures are inevitable, especially in densely populated areas. The undesired effects resulting from this proximity are manifold and range from the establishment of hazardous potentials to the outbreak of alternate current corrosion phenomena. The study of this class of problems is necessary for ensuring security in the vicinities of the interaction zone and also to preserve the integrity of the equipment and of the devices there present. However, the complete modeling of this type of application requires the three- -dimensional representation of the region of interest and needs specific numerical methods for field computation. In this work, the modeling of problems arising from the flow of electrical currents in the ground (the so-called conductive coupling) will be addressed with the finite element method. Those resulting from the time variation of the electromagnetic fields (the so-called inductive coupling) will be considered as well, and they will be treated with the generalized PEEC (Partial Element Equivalent Circuit) method. More specifically, a special boundary condition on the electric potential is proposed for truncating the computational domain in the finite element analysis of conductive coupling problems, and a complete PEEC formulation for modeling inductive coupling problems is presented. Test configurations of increasing complexities are considered for validating the foregoing approaches. These works aim to provide a contribution to the modeling of this class of problems, which tend to become common with the expansion of power grids.

Amination of enzymes to improve biocatalyst performance: coupling genetic modification and physicochemical tools

Improvement of the features of an enzyme is in many instances a pre-requisite for the industrial implementation of these exceedingly interesting biocatalysts. To reach this goal, the researcher may utilize different tools. For example, amination of the enzyme surface produces an alteration of the isoelectric point of the protein along with its chemical reactivity (primary amino groups are the most widely used to obtain the reaction of the enzyme with surfaces, chemical modifiers, etc.) and even its “in vivo” behavior. This review will show some examples of chemical (mainly modifying the carboxylic groups using the carbodiimide route), physical (using polycationic polymers like polyethyleneimine) and genetic amination of the enzyme surface. Special emphasis will be put on cases where the amination is performed to improve subsequent protein modifications. Thus, amination has been used to increase the intensity of the enzyme/support multipoint covalent attachment, to improve the interaction with cation exchanger supports or polymers, or to promote the formation of crosslinkings (both intra-molecular and in the production of crosslinked enzyme aggregates). In other cases, amination has been used to directly modulate the enzyme properties (both in immobilized or free form). Amination of the enzyme surface may also pursue other goals not related to biocatalysis. For example, it has been used to improve the raising of antibodies against different compounds (both increasing the number of haptamers per enzyme and the immunogenicity of the composite) or the ability to penetrate cell membranes. Thus, amination may be a very powerful tool to improve the use of enzymes and proteins in many different areas and a great expansion of its usage may be expected in the near future.

Improvements in DFT Calculations of Spin–Spin Coupling Constants

Different types of spin–spin coupling constants (SSCCs) for several representative small molecules are evaluated and analyzed using a combination of 10 exchange functionals with 12 correlation functionals. For comparison, calculations performed using MCSCF, SOPPA, other common DFT methods, and also experimental data are considered. A detailed study of the percentage of Hartree–Fock exchange energy in SSCCs and in its four contributions is carried out. From the above analysis, a combined functional formed with local Slater (34%), Hartree–Fock exchange (66%), and P86 correlation functional (S66P86) is proposed in this paper. The accuracy of the values obtained with this hybrid functional (mean absolute deviation of 4.5 Hz) is similar to that of the SOPPA method (mean absolute deviation of 4.6 Hz).

Assessing the Single Supervisory Mechanism: passing the point of no return for Europe's banking union. Egmont Paper No. 58, June 2013

From the Executive Summary. Europe’s financial and sovereign debt crises have become increasingly interconnected. In order to break the negative feedback loop between the two, the EU has decided to create a common supervisory framework for the banking sector: the Single Supervisory Mechanism (SSM). The SSM will involve a supervisory system including both the national supervisors and the European Central Bank (ECB). By endowing the ECB with supervisory authority over a major part of the European banking sector, the SSM’s creation will result in a shake-up of the way in which the European financial sector is being supervised. Under the right circumstances, this could be a major step forward in addressing Europe’s interconnected crises.

Historic turning point or just another chapter? The results of a foreign-policy summit. EPC Post-Summit Analysis, 24 March 2014

The March 2014 European Council could enter the history books as a turning point, not only in the EU’s relations with Russia but also in its role as a foreign policy actor. Events in Ukraine inevitably dominated the Summit, with EU leaders adopting a balanced approach aimed at achieving three key objectives – de-escalation, containment/deterrence and cooperation – based on political and economic support for Ukraine, increased but limited pressure on Russia, and moves to strengthen ties with other EU neighbours. The Summit also discussed a range of economic and environmental policy issues, with the situation in Ukraine casting a long shadow over the discussion on energy policy, but failed to reach agreement on the EU’s climate goals to 2030, or to put more flesh on the bones of calls for a European “industrial renaissance”. However, two other developments were particularly significant: the creation of the second pillar of the future banking union, establishing a single regime for winding down failing banks; and changes to the savings tax directive, bringing years of negotiation to an end.

The Berlin Republic. Evolution of Germany's politics of memory and German patriotism. OSW Point of View, August 2009

Important changes have occurred in recent years in the attitude of a majority of the German elite towards the history of the 20th century and the political identity built on collective memory. Until recently, the sense of guilt for the crimes of the Third Reich and the obligation to remember were prevalent. While these two elements of Germany's memory of World War II are still important, currently the focus increasingly shifts to the German resistance against Nazism and the fate of the Germans who suffered in the war. Positive references to Germany's post-war history also occupy more and more space in the German memory. In 2009, i.e. the year of the 60th anniversary of the Federal Republic of Germany and the 20th anniversary of the fall of Communism, the efforts of German public institutions concentrate on promoting a new canon of history built around the successful democratisation and Germany's post-war economic success. The purpose behind these measures is to build a common historical memory that could be shared by the eastern and western parts of Germany and appeal to Germany's immigrants, who account for a growing proportion of the society.