Nuclear powerplant standardization : light water reactors.

"April 1981" (v. 1); "June 1981" (v. 2)

CFD-modeling and experiments of insulation debris transport phenomena in water flow

The investigation of insulation debris generation, transport, and sedimentation becomes more important with regard to reactor safety research for pressurized water reactors and boiling water reactors when considering the long-term behavior of emergency core coolant systems during all types of loss-of-coolant accidents (LOCAs). The insulation debris released near the break during a LOCA incident consists of a mixture of disparate particle populations that varies with size, shape, consistency, and other properties. Some fractions of the released insulation debris can be transported into the reactor sump, where it may perturb/impinge on the emergency core cooling systems. Open questions of generic interest are, for example, the particle load on strainers and corresponding pressure drop, the sedimentation of the insulation debris in a water pool, and its possible resuspension and transport in the sump water flow. A joint research project on such questions is being performed in cooperation with the University of Applied Sciences Zittau/Görlitz. The project deals with the experimental investigation and the development of computational fluid dynamics (CFD) models for the description of particle transport phenomena in coolant flow. While the experiments are performed at the University of Applied Sciences Zittau/Görlitz, the theoretical work is concentrated at Forschungszentrum Dresden-Rossendorf. In the current paper the basic concepts for CFD modeling are described and feasibility studies including the conceptual design of the experiments are presented.

Numerical and experimental investigations for insulation particle transport phenomena in water flow

The investigation of insulation debris generation, transport and sedimentation becomes more important with regard to reactor safety research for pressurized and boiling water reactors, when considering the long-term behaviour of emergency core coolant systems during all types of loss of coolant accidents (LOCA). The insulation debris released near the break during a LOCA incident consists of a mixture of a disparate particle population that varies with size, shape, consistency and other properties. Some fractions of the released insulation debris can be transported into the reactor sump, where it may perturb or impinge on the emergency core cooling systems. Open questions of generic interest are for example the particle load on strainers and corresponding pressure-drop, the sedimentation of the insulation debris in a water pool, its possible re-suspension and transport in the sump water flow. A joint research project on such questions is being performed in cooperation with the University of Applied Science Zittau/Görlitz and the Forschungszentrum Dresden-Rossendorf. The project deals with the experimental investigation and the development of computational fluid dynamic (CFD) models for the description of particle transport phenomena in coolant flow. While the experiments are performed at the University Zittau/Görlitz, the theoretical work is concentrated at Forschungszentrum Dresden-Rossendorf. In the present paper, the basic concepts for computational fluid dynamic (CFD) modelling are described and experimental results are presented. Further experiments are designed and feasibility studies were performed.

Degradation of compounds present in cork boiling water by gamma radiation

Cork boiling water is an aqueous and complex dark liquor with high concentration of phenolic compounds such as phenolic acids and tannins [1, 2], which are considered biorecalcitrants [2]. Ionizing radiation has been widely studied as an alternative technology for the degradation of organic contaminants without the addition of any other (e.g.: Fenton technologies). The aim of this work was to identify the compounds present in cork boiling water and further evaluate the resulting stable degradation products after gamma irradiation. The irradiation experiments of standard solutions were carried out at room temperature using a Co-60 experimental equipment. The applied absorbed doses were 20 and 50 kGy at a dose rate of 1.5 kGy/h, determined by routine dosimeters [3]. The identification of radiolytic products was carried out by HPLC-DAD-ESI/MS. The phenolic compounds were identified by comparing their retention times and UV–vis and mass spectra with those obtained from standard compounds, when available, as well as by comparing the obtained information with available data reported in the literature. Concerning the obtained results and the literature review, the main cork wastewater components are: quinic, gallic, protocatechuic, vanillic, syringic and ellagic acids. Based on this, we used protocatechuic, vanillic and syringic acids as model compounds to study their degradation by gamma radiation in order to identify the corresponding radiolytic products. Standard aqueous solutions were irradiated and the derivatives of each model compound are represented in figure 1. The obtained results seem to demonstrate that the derivatives of the parent compounds could also be phenolic acids, since it was observed the loss of 44 u (CO2) from the [M-H]- ions. Gallic and protocatechuic acids are identified as derivatives of vanillic and syringic acids, and gallic acid as a protocatechuic acid derivative. Compound 5 ([M-H]- at m/z 169) was tentatively identified as 2,4,6-trihydroxybenzoic acid, since its fragmentation pattern (m/z 151, 125 and 107) is similar to that previously reported in literature [4]. The structure of compound 7 was proposed based on the molecular ion and its fragmentation and compound 6 remains unknown.

Calculation of the eigenfunctions of the two-group neutron diffusion equation and application to modal decomposition of BWR instabilities

In this thesis, numerical methods aiming at determining the eigenfunctions, their adjoint and the corresponding eigenvalues of the two-group neutron diffusion equations representing any heterogeneous system are investigated. First, the classical power iteration method is modified so that the calculation of modes higher than the fundamental mode is possible. Thereafter, the Explicitly-Restarted Arnoldi method, belonging to the class of Krylov subspace methods, is touched upon. Although the modified power iteration method is a computationally-expensive algorithm, its main advantage is its robustness, i.e. the method always converges to the desired eigenfunctions without any need from the user to set up any parameter in the algorithm. On the other hand, the Arnoldi method, which requires some parameters to be defined by the user, is a very efficient method for calculating eigenfunctions of large sparse system of equations with a minimum computational effort. These methods are thereafter used for off-line analysis of the stability of Boiling Water Reactors. Since several oscillation modes are usually excited (global and regional oscillations) when unstable conditions are encountered, the characterization of the stability of the reactor using for instance the Decay Ratio as a stability indicator might be difficult if the contribution from each of the modes are not separated from each other. Such a modal decomposition is applied to a stability test performed at the Swedish Ringhals-1 unit in September 2002, after the use of the Arnoldi method for pre-calculating the different eigenmodes of the neutron flux throughout the reactor. The modal decomposition clearly demonstrates the excitation of both the global and regional oscillations. Furthermore, such oscillations are found to be intermittent with a time-varying phase shift between the first and second azimuthal modes.

Leak Detection In Pressure Tubes Of A Pressurized Heavy-Water Reactor By Acoustic-Emission Technique

Leak detection in the fuel channels is one of the challenging problems during the in-service inspection (ISI) of Pressurised Heavy Water Reactors (PHWRs). In this paper, the use of an acoustic emission (AE) technique together with AE signal analysis is described, to detect a leak that was ncountered in one (or more) of the 306 fuel channels of the Madras Atomic Power Station (PHWR), Unit I. The paper describes the problems encountered during the ISI, the experimental methods adopted and the results obtained. Results obtained using acoustic emission signal analysis are compared with those obtained from other leak detection methods used in such cases.

Near-complete transuranic waste incineration in a thorium fuelled pressurised water reactor

The production of long-lived transuranic (TRU) waste is a major disadvantage of fission-based nuclear power. Incineration, and virtual elimination, of waste stockpiles is possible in a thorium (Th) fuelled critical or subcritical fast reactor. Fuel cycles producing a net decrease in TRUs are possible in conventional pressurised water reactors (PWRs). However, minor actinides (MAs) have a detrimental effect on reactivity and stability, ultimately limiting the quality and quantity of waste that can be incinerated. In this paper, we propose using a thorium-retained-actinides fuel cycle in PWRs, where the reactor is fuelled with a mixture of thorium and TRU waste, and after discharge all actinides are reprocessed and returned to the reactor. To investigate the feasibility and performance of this fuel cycle an assembly-level analysis for a one-batch reloading strategy was completed over 125 years of operation using WIMS 9. This one-batch analysis was performed for simplicity, but allowed an indicative assessment of the performance of a four-batch fuel management strategy. The build-up of 233U in the reactor allowed continued reactive and stable operation, until all significant actinide populations had reached pseudo-equilibrium in the reactor. It was therefore possible to achieve near-complete transuranic waste incineration, even for fuels with significant MA content. The average incineration rate was initially around 330 kg per GW th year and tended towards 250 kg per GW th year over several decades: a performance comparable to that achieved in a fast reactor. Using multiple batch fuel management, competitive or improved end-of-cycle burn-up appears achievable. The void coefficient (VC), moderator temperature coefficient (MTC) and Doppler coefficient remained negative. The quantity of soluble boron required for a fixed fuel cycle length was comparable to that for enriched uranium fuel, and acceptable amounts can be added without causing a positive VC or MTC. This analysis is limited by the consideration of a single fuel assembly, and it will be necessary to perform a full core coupled neutronic-thermal-hydraulic analysis to determine if the design in its current form is feasible. In particular, the potential for positive VCs if the core is highly or locally voided is a cause for concern. However, these results provide a compelling case for further work on concept feasibility and fuel management, which is in progress. © 2011 Elsevier Ltd. All rights reserved.

Decay heat in fast reactors with transuranic fuels

In this work, we performed an evaluation of decay heat power of advanced, fast spectrum, lead and molten salt-cooled reactors, with flexible conversion ratio. The decay heat power was calculated using the BGCore computer code, which explicitly tracks over 1700 isotopes in the fuel throughout its burnup and subsequent decay. In the first stage, the capability of the BGCore code to accurately predict the decay heat power was verified by performing a benchmark calculation for a typical UO2 fuel in a Pressurized Water Reactor environment against the (ANSI/ANS-5.1-2005, "Decay Heat Power in Light Water Reactors," American National Standard) standard. Very good agreement (within 5%) between the two methods was obtained. Once BGCore calculation capabilities were verified, we calculated decay power for fast reactors with different coolants and conversion ratios, for which no standard procedure is currently available. Notable differences were observed for the decay power of the advanced reactor as compared with the conventional UO2 LWR. The importance of the observed differences was demonstrated by performing a simulation of a Station Blackout transient with the RELAP5 computer code for a lead-cooled fast reactor. The simulation was performed twice: using the code-default ANS-79 decay heat curve and using the curve calculated specifically for the studied core by BGCore code. The differences in the decay heat power resulted in failure to meet maximum cladding temperature limit criteria by ∼100 °C in the latter case, while in the transient simulation with the ANS-79 decay heat curve, all safety limits were satisfied. The results of this study show that the design of new reactor safety systems must be based on decay power curves specific to each individual case in order to assure the desired performance of these systems. © 2009 Elsevier B.V. All rights reserved.

Light water breeding with nitride fuel

This scoping study proposes using mixed nitride fuel in Pu-based high conversion LWR designs in order to increase the breeding ratio. The higher density fuel reduces the hydrogen-to-heavy metal ratio in the reactor which results in a harder spectrum in which breeding is more effective. A Resource-renewable Boiling Water Reactor (RBWR) assembly was modeled in MCNP to demonstrate this effect in a typical high conversion LWR design. It was determined that changing the fuel from (U,TRU)O2 to (U,TRU)N in the assembly can increase its fissile inventory ratio (fissile Pu mass divided by initial fissile Pu mass) from 1.04 to up to 1.17. © 2011 Elsevier Ltd. All rights reserved.

Rethinking Rice Preparation for Highly Efficient Removal of Inorganic Arsenic Using Percolating Cooking Water

A novel way of cooking rice to maximize the removal of the carcinogen inorganic arsenic (Asi) is presented here. In conventional rice cooking water and grain are in continuous contact, and it is known that the larger the water:rice cooking ratio, the more Asi removed by cooking, suggesting that the Asi in the grain is mobile in water. Experiments were designed where rice is cooked in a continual stream of percolating near boiling water, either low in Asi, or Asi free. This has the advantage of not only exposing grain to large volumes of cooking water, but also physically removes any Asi leached from the grain into the water receiving vessel. The relationship between cooking water volume and Asi removal in conventional rice cooking was demonstrated for the rice types under study. At a water-to-rice cooking ratio of 12:1, 57±5% of Asi could be removed, average of 6 wholegrain and 6 polished rice samples. Two types of percolating technology were tested, one where the cooking water was recycled through condensing boiling water steam and passing the freshly distilled hot water through the grain in a laboratory setting, and one where tap water was used to cook the rice held in an off-the-shelf coffee percolator in a domestic setting. Both approaches proved highly effective in removing Asi from the cooking rice, with up to 85% of Asi removed from individual rice types. For the recycled water experiment 59±8% and 69±10% of Asi was removed, on average, compared to uncooked rice for polished (n=27) and wholegrain (n=13) rice, respectively. For coffee percolation there was no difference between wholegrain and polished rice, and the effectiveness of Asi removal was 49±7% across 6 wholegrain and 6 polished rice samples. The manuscript explores the potential applications and further optimization of this percolating cooking water, high Asi removal, discovery.

Epoxide ring opening in a zinc(II) complex in water without any Lewis acid catalyst: formation of only one diastereomer out of 2

The epoxide ring in 5,6-dihydro-5,6-epoxy-1,10-phenanthroline (L) opens up in its reaction with 4-methylaniline and 4-methoxyaniline in water in equimolar proportion at room temperature without any Lewis acid catalyst to give a monohydrate of 6-(4-methyl-phenylamino)-5,6-dihydro-1,10-phenanthrolin-5-ol (L′·H2O) and 6-(4-methoxyphenyl-amino)-5,6-dihydro-1,10-phenanthrolin-5-ol (L″) respectively. Reaction time decreases from 72 to 14 h in boiling water. But the yields become less. Reaction of L with Zn(ClO4)2·6H2O in methanol in 3:1 molar ratio at room temperature affords white [ZnL3](ClO4)2·H2O. The X-ray crystal structure of the acetonitrile solvate [ZnL3](ClO4)2·MeCN has been determined which shows that the metal has a distorted octahedral N6 coordination sphere. [ZnL3](ClO4)2·2H2O reacts with 4-methylaniline and 4-methoxyaniline in boiling water in 1:3 molar proportion in the absence of any Lewis acid catalyst to produce [ZnL′3](ClO4)2·4H2O and [ZnL″3](ClO4)2·H2O, respectively in 1–4 h time in somewhat low yield. In the 1H NMR spectra of [ZnL′3](ClO4)2·4H2O and [ZnL″3](ClO4)2·H2O, only one sharp methyl signal is observed implicating that only one diastereomer out of the 23 possibilities is formed. The same diastereomers are obtained when L′·H2O and L″ are reacted directly with Zn(ClO4)2·6H2O in tetrahydrofuran at room temperature in very good yields. Reactions of L′·H2O and L″ with Ru(phen)2Cl2·2H2O (phen = 1,10-phenanthroline) in equimolar proportion in methanol–water mixture under refluxing condition lead to the isolation of two diastereomers of [Ru(phen)2L′](ClO4)2·2H2O and [Ru(phen)2L″](ClO4)2·2H2O.

Assessment of CTF boiling transition and critical heat flux modeling capabilities using the OECD/NRC BFBT and PSBT benchmark databases

The need to refine models for best-estimate calculations, based on good-quality experimental data, has been expressed in many recent meetings in the field of nuclear applications. The modeling needs arising in this respect should not be limited to the currently available macroscopic methods but should be extended to next-generation analysis techniques that focus on more microscopic processes. One of the most valuable databases identified for the thermalhydraulics modeling was developed by the Nuclear Power Engineering Corporation (NUPEC), Japan. From 1987 to 1995, NUPEC performed steady-state and transient critical power and departure from nucleate boiling (DNB) test series based on the equivalent full-size mock-ups. Considering the reliability not only of the measured data, but also other relevant parameters such as the system pressure, inlet sub-cooling and rod surface temperature, these test series supplied the first substantial database for the development of truly mechanistic and consistent models for boiling transition and critical heat flux. Over the last few years the Pennsylvania State University (PSU) under the sponsorship of the U.S. Nuclear Regulatory Commission (NRC) has prepared, organized, conducted and summarized the OECD/NRC Full-size Fine-mesh Bundle Tests (BFBT) Benchmark. The international benchmark activities have been conducted in cooperation with the Nuclear Energy Agency/Organization for Economic Co-operation and Development (NEA/OECD) and Japan Nuclear Energy Safety (JNES) organization, Japan. Consequently, the JNES has made available the Boiling Water Reactor (BWR) NUPEC database for the purposes of the benchmark. Based on the success of the OECD/NRC BFBT benchmark the JNES has decided to release also the data based on the NUPEC Pressurized Water Reactor (PWR) subchannel and bundle tests for another follow-up international benchmark entitled OECD/NRC PWR Subchannel and Bundle Tests (PSBT) benchmark. This paper presents an application of the joint Penn State University/Technical University of Madrid (UPM) version of the well-known subchannel code COBRA-TF, namely CTF, to the critical power and departure from nucleate boiling (DNB) exercises of the OECD/NRC BFBT and PSBT benchmarks

Analysis of a Spanish energy scenario with Generation IV nuclear reactors

The advantages of fast-spectrum reactors consist not only of an efficient use of fuel through the breeding of fissile material and the use of natural or depleted uranium, but also of the potential reduction of the amount of actinides such as americium and neptunium contained in the irradiated fuel. The first aspect means a guaranteed future nuclear fuel supply. The second fact is key for high-level radioactive waste management, because these elements are the main responsible for the radioactivity of the irradiated fuel in the long term. The present study aims to analyze the hypothetical deployment of a Gen-IV Sodium Fast Reactor (SFR) fleet in Spain. A nuclear fleet of fast reactors would enable a fuel cycle strategy different than the open cycle, currently adopted by most of the countries with nuclear power. A transition from the current Gen-II to Gen-IV fleet is envisaged through an intermediate deployment of Gen-III reactors. Fuel reprocessing from the Gen-II and Gen-III Light Water Reactors (LWR) has been considered. In the so-called advanced fuel cycle, the reprocessed fuel used to produce energy will breed new fissile fuel and transmute minor actinides at the same time. A reference case scenario has been postulated and further sensitivity studies have been performed to analyze the impact of the different parameters on the required reactor fleet. The potential capability of Spain to supply the required fleet for the reference scenario using national resources has been verified. Finally, some consequences on irradiated final fuel inventory are assessed. Calculations are performed with the Monte Carlo transport-coupled depletion code SERPENT together with post-processing tools.

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.

The slowing down distribution from a point source of fission neutrons in light water /

Work performed at the Physics Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee.