Optimization of multidimensional cross-section tables for few-group core calculations

Multigroup diffusion codes for three dimensional LWR core analysis use as input data pre-generated homogenized few group cross sections and discontinuity factors for certain combinations of state variables, such as temperatures or densities. The simplest way of compiling those data are tabulated libraries, where a grid covering the domain of state variables is defined and the homogenized cross sections are computed at the grid points. Then, during the core calculation, an interpolation algorithm is used to compute the cross sections from the table values. Since interpolation errors depend on the distance between the grid points, a determined refinement of the mesh is required to reach a target accuracy, which could lead to large data storage volume and a large number of lattice transport calculations. In this paper, a simple and effective procedure to optimize the distribution of grid points for tabulated libraries is presented. Optimality is considered in the sense of building a non-uniform point distribution with the minimum number of grid points for each state variable satisfying a given target accuracy in k-effective. The procedure consists of determining the sensitivity coefficients of k-effective to cross sections using perturbation theory; and estimating the interpolation errors committed with different mesh steps for each state variable. These results allow evaluating the influence of interpolation errors of each cross section on k-effective for any combination of state variables, and estimating the optimal distance between grid points.

Assessment of damage domains of the High-Temperature Engineering Test Reactor (HTTR)

This paper presents an assessment analysis of damage domains of the 30 MWth prototype High-Temperature Engineering Test Reactor (HTTR) operated by the Japan Atomic Energy Agency (JAEA). For this purpose, an in-house deterministic risk assessment computational tool was developed based on the Theory of Stimulated Dynamics (TSD). To illustrate the methodology and applicability of the developed modelling approach, assessment results of a control rod (CR) withdrawal accident during subcritical conditions are presented and compared with those obtained by the JAEA.

Revising the Emergency Management Requirements for new generation reactors

The paper presents the application of a new risk-informed methodology for the identification of the Emergency Management Requirements (EMR) to a Generation II, Large size Reactor and a Generation III+ Small Modular Reactor. The results obtained in this test case demonstrate that the actual EMR is conservative, as expected, for the GenII reactor, while the new methodology could be applied for the definition of EMRs for the new generation Nuclear Power Plants, with a possible reduction of the emergency area without loss of safety level. By adopting both probabilistic and deterministic approaches, the study addresses possible accidents and corresponding release scenarios for the two types of reactor, calculates the areas where the accidents have an impact on the population and defines the new EMR considering the health effects on the population.

Code assessment and modelling for Design Basis Accident Analysis of the European sodium fast reactor design. Part I: System description, modelling and benchmarking

The new reactor concepts proposed in the Generation IV International Forum (GIF) are conceived to improve the use of natural resources, reduce the amount of high-level radioactive waste and excel in their reliability and safe operation. Among these novel designs sodium fast reactors (SFRs) stand out due to their technological feasibility as demonstrated in several countries during the last decades. As part of the contribution of EURATOM to GIF the CP-ESFR is a collaborative project with the objective, among others, to perform extensive analysis on safety issues involving renewed SFR demonstrator designs. The verification of computational tools able to simulate the plant behaviour under postulated accidental conditions by code-to-code comparison was identified as a key point to ensure reactor safety. In this line, several organizations employed coupled neutronic and thermal-hydraulic system codes able to simulate complex and specific phenomena involving multi-physics studies adapted to this particular fast reactor technology. In the “Introduction” of this paper the framework of this study is discussed, the second section describes the envisaged plant design and the commonly agreed upon modelling guidelines. The third section presents a comparative analysis of the calculations performed by each organisation applying their models and codes to a common agreed transient with the objective to harmonize the models as well as validating the implementation of all relevant physical phenomena in the different system codes.

Code assessment and modelling for Design Basis Accident analysis of the European Sodium Fast Reactor design. Part II: Optimised core and representative transients analysis

The new reactor concepts proposed in the Generation IV International Forum require the development and validation of computational tools able to assess their safety performance. In the first part of this paper the models of the ESFR design developed by several organisations in the framework of the CP-ESFR project were presented and their reliability validated via a benchmarking exercise. This second part of the paper includes the application of those tools for the analysis of design basis accident (DBC) scenarios of the reference design. Further, this paper also introduces the main features of the core optimisation process carried out within the project with the objective to enhance the core safety performance through the reduction of the positive coolant density reactivity effect. The influence of this optimised core design on the reactor safety performance during the previously analysed transients is also discussed. The conclusion provides an overview of the work performed by the partners involved in the project towards the development and enhancement of computational tools specifically tailored to the evaluation of the safety performance of the Generation IV innovative nuclear reactor designs.

Passive neutron area monitor with pairs of TLDs as neutron detector

A passive neutron area monitor has been designed using Monte Carlo methods; the monitor is a polyethylene cylinder with pairs of thermoluminescent dosimeters (TLD600 and TLD700) as thermal neutron detector. The monitor was calibrated with a bare and a thermalzed 241AmBe neutron sources and its performance was evaluated measuring the ambient dose equivalent due to photoneutrons produced by a 15 MV linear accelerator for radiotherapy and the neutrons in the output of a TRIGA Mark III radial beam port.

Potential issues related to emergency core cooling system strainers performance at boiling water reactors: application to Cofrentes NPP (Spain)

From the 60s to the 90s, a great number of events related to the Emergency Core Cooling Systems Strainers have been happened in all kind of reactors all over the world. Thus, the Nuclear Regulatory Commission of the USA emitted some Bulletins to address the concerns about the adequacy of Emergency Core Cooling Systems (ECCS) strainer performance at boiling water reactors (BWR). In Spain the regulatory body (Consejo de Seguridad Nuclear, CSN) adopted the USA regulation and Cofrentes NPP installed new strainers with a considerable bigger size than the old strainers. The nuclear industry conducted significant and extensive research, guidance development, testing, reviews, and hardware and procedure changes during the 90s to resolve the issues related to debris blockage of BWR strainers. In 2001 the NRC and CSN closed the Bulletins. Thereafter, the strainers issues were moved to the PWR reactors. In 2004 the NRC issued a Generic Letter (GL). It requested the resolution of several effects which were not noted in the past. The GL regarded to be resolved by the PWR reactors but the NRC in USA and the CSN in Spain have requested that the BWR reactors investigate differences between the methodologies used by the BWRs and PWRs. The developments and improvements done for Cofrentes NPP are detailed. Studies for this plant show that the head loss due to the considered debris is at most half of the limited head loss for the ECCS strainer and the NPSH (Net Positive Suction Head) required for the ECCS pumps is at least three times lower than the NPSH available.

Integrated deterministic-probabilistic safety assessment methodologies

(ENG) IDPSA (Integrated Deterministic-Probabilistic Safety Assessment) is a family of methods which use tightly coupled probabilistic and deterministic approaches to address respective sources of uncertainties, enabling Risk informed decision making in a consistent manner. The starting point of the IDPSA framework is that safety justification must be based on the coupling of deterministic (consequences) and probabilistic (frequency) considerations to address the mutual interactions between stochastic disturbances (e.g. failures of the equipment, human actions, stochastic physical phenomena) and deterministic response of the plant (i.e. transients). This paper gives a general overview of some IDPSA methods as well as some possible applications to PWR safety analyses (SPA)DPSA (Metodologías Integradas de Análisis Determinista-Probabilista de Seguridad) es un conjunto de métodos que utilizan métodos probabilistas y deterministas estrechamente acoplados para abordar las respectivas fuentes de incertidumbre, permitiendo la toma de decisiones Informada por el Riesgo de forma consistente. El punto de inicio del marco IDPSA es que la justificación de seguridad debe estar basada en el acoplamiento entre consideraciones deterministas (consecuencias) y probabilistas (frecuencia) para abordar la interacción mutua entre perturbaciones estocásticas (como por ejemplo fallos de los equipos, acciones humanas, fenómenos físicos estocásticos) y la respuesta determinista de la planta (como por ejemplo los transitorios). Este artículo da una visión general de algunos métodos IDSPA así como posibles aplicaciones al análisis de seguridad de los PWR.

Dynamics of gas bubbles encapsulated by a viscoelastic fluid shell under acoustic fields

The dynamics of a gas-filled microbubble encapsulated by a viscoelastic fluid shell immersed in a Newtonian liquid and subject to an external pressure field is theoretically studied. The problem is formulated by considering a nonlinear Oldroyd type constitutive equation to model the rheological behavior of the fluid shell. Heat and mass transfer across the surface bubble have been neglected but radiation losses due to the compressibility of the surrounding liquid have been taken into account. Bubble collapse under sudden increase of the external pressure as well as nonlinear radial oscillations under ultrasound fields are investigated. The numerical results obtained show that the elasticity of the fluid coating intensifies oscillatory collapse and produces a strong increase of the amplitudes of radial oscillations which may become chaotic even for moderate driving pressure amplitudes. The role played by the elongational viscosity has also been analyzed and its influence on both, bubble collapse and radial oscillations, has been recognized. According to the theoretical predictions provided in the present work, a microbubble coated by a viscoelastic fluid shell is an oscillating system that, under acoustic driving, may experience volume oscillations of large amplitude, being, however, more stable than a free bubble. Thus, it could be expected that such a system may have a suitable behavior as an echogenic agent.

Desarrollo de una herramienta de verificación para cálculos de difusión mediante COBAYA

El código COBAYA4 es un simulador de núcleo multi-escala que resuelve la ecuación de difusión 3D en multigrupos en geometría cartesiana y hexagonal[3]. Este código ha sido desarrollado en el Departamento de Ingeniería Nuclear desde los años 80[2] ampliando su alcance y funcionalidades de forma continua. Como parte de estos desarrollos es necesaria la verificación continua de que el código sigue teniendo al menos las mismas capacidades que tenía anteriormente. Además es necesario establecer casos de referencia que nos permitan confirmar que los resultados son comparables a los obtenidos con otros códigos con modelos de mayor precisión. El desarrollo de una herramienta informática que automatice la comparación de resultados con versiones anteriores del código y con resultados obtenidos mediante modelos de mayor precisión es crucial para implementar en el código nuevas funcionalidades. El trabajo aquí presentado ha consistido en la generación de la mencionada herramienta y del conjunto de casos de referencia que han constituido la matriz mencionada.

Metodología de simulación numérica de la transferencia de calor por ebullición en película en torno a superficies esféricas refrigeradas por líquido saturado en convección libre: aplicación a reactores nucleares de lecho de bolas (PBR)

En la presente Tesis se realizó un análisis numérico, usando el código comercial Ansys-Fluent, de la refrigeración de una bola de combustible de un reactor de lecho de bolas (PBR, por sus siglas en inglés), ante un escenario de emergencia en el cual el núcleo sea desensamblado y las bolas se dejen caer en una piscina de agua, donde el mecanismo de transferencia de calor inicialmente sería la ebullición en película, implicando la convección y la radiación al fluido. Previamente se realizaron pruebas de validación, comparando los resultados numéricos con datos experimentales disponibles en la literatura para tres geometrías diferentes, lo cual permitió seleccionar los esquemas y modelos numéricos con mejor precisión y menor costo computacional. Una vez identificada la metodología numérica, todas las pruebas de validación fueron ampliamente satisfactorias, encontrándose muy buena concordancia en el flujo de calor promedio con los datos experimentales. Durante estas pruebas de validación se lograron caracterizar numéricamente algunos parámetros importantes en la ebullición en película con los cuales existen ciertos niveles de incertidumbre, como son el factor de acoplamiento entre convección y radiación, y el factor de corrección del calor latente de vaporización. El análisis térmico de la refrigeración de la bola del reactor por ebullición en película mostró que la misma se enfría, a pesar del calor de decaimiento, con una temperatura superficial de la bola que desciende de forma oscilatoria, debido al comportamiento inestable de la película de vapor. Sin embargo, la temperatura de esta superficie tiene una buena uniformidad, notándose que las áreas mejor y peor refrigeradas están localizadas en la parte superior de la bola. Se observó la formación de múltiples domos de vapor en diferentes posiciones circunferenciales, lo cual causa que el área más caliente de la superficie se localice donde se forman los domos más grandes. La separación entre los domos de vapor fue consistente con la teoría hidrodinámica, con la adición de que la separación entre domos se reduce a medida que evolucionan y crecen, debido a la curvatura de la superficie. ABSTRACT A numerical cooling analysis of a PBR fuel pebble, after an emergency scenario in which the nucleus disassembly is made and the pebbles are dropped into a water pool, transmitting heat by film boiling, involving convection and radiation to the fluid, is carried out in this Thesis. First, were performed validation tests comparing the numerical results with experimental works available for three different geometries, which allowed the selection of numerical models and schemes with better precision and lower computational cost. Once identified the numerical methodology, all validation tests were widely satisfactory, finding very good agreement with experimental works in average heat flux. During these validation tests were achieved numerically characterize some important parameters in film boiling with which there are certain levels of uncertainty, such as the coupling factor between convection and radiation, and the correction factor of the latent heat of vaporization. The thermal analysis of pebble cooling by film boiling shows that despite its decay heat, cooling occurs, with pebble surface temperature descending from an oscillatory manner, due to the instability of the vapor film. However, the temperature of this surface has a good uniformity, noting that the best and worst refrigerated area is located at the top of the pebble. The formation of multiple vapor domes at different circumferential positions is observed, which cause that the hottest area of the surface was located where biggest vapor domes were formed. The separation between vapor domes was consistent with the hydrodynamic theory, with the addition that the separation is reduced as the vapor dome evolves and grows, due to the surface curvature.

Sistemas de salvaguardia en reactores EPR

En este documento se describe brevemente el funcionamiento de los diversos sistemas de una planta nuclear operada con un reactor de tipo PWR. Más concretamente, el proyecto se centra en una descripción exhaustiva de los sistemas de salvaguardia y seguridad que regulan el funcionamiento de un reactor de tipo EPR, así como la central nuclear que contiene a dicho reactor. El proceso ha consistido en clasificar y resumir los distintos sistemas que operan en dicha planta, estudiando sus características y parámetros de funcionamiento. También se han estudiado los accidentes más comunes que pueden tener lugar en este tipo de centrales nucleares. Tras el análisis y estudio realizado acerca del reactor EPR, se puede concluir que las centrales nucleares que operan con este tipo de reactor experimentan una serie de mejoras en cuanto a la prevención de accidentes, así como una serie de mejoras de diseño en una gran variedad de sistemas y elementos del reactor, como pueden ser la vasija, los SG, etc. ABSTRACT This document gives a brief description of the operation of several systems of a nuclear power plant operating with a PWR reactor. More specifically, the project focuses on a thorough description of the safety and security systems that govern the operation of an EPR reactor and its plant. The process consisted on classify and summarize the different operating systems of this nuclear plant, studying its characteristics and operating parameters. We have also studied the most common accidents that can occur in this type of nuclear power plants. After the analysis and study on the EPR reactor, it can be concluded that nuclear power plants operating with this type of reactor undergo a series of improvements in the prevention of accidents, as well as a number of design improvements in several reactor systems and components, such as the vessel, the SG, etc.

Análisis de la distribución de hidrógeno en contención y pozo seco de C.N. Cofrentes.

El desarrollo del modelo de la Contención de C.N. Cofrentes mediante GOTHIC se ha llevado a cabo introduciendo todos los datos geométricos y de estructuras de la Contención, pudiendo así modelar todos los recintos interiores y habitaciones que la componen. De esta forma se ha obtenido un modelo 3D detallado y con la precisión suficiente para el estudio global de la gestión del hidrógeno, permitiendo tener en cuenta, a la hora de la distribución del hidrógeno, la asimetría tanto de la contención como de las descargas de masa y en energía que en ella se realizan, permitiendo simular la distribución del vapor y el hidrógeno presentes en el accidente severo para poder determinar las zonas de mayor riesgo de deflagración o detonación durante la evolución del accidente

Estudio de la distribución de hidrógeno en una contención PWR con códigos CFD

Durante el desarrollo de un accidente severo en un reactor PWR, se pueden generar grandes cantidades de hidrógeno por la oxidación de los metales presentes en el núcleo, principalmente el zirconio de las vainas del combustible. Este hidrógeno, junto con vapor y otros gases, puede ser liberado a la atmósfera de la contención por una fuga o rotura en el circuito primario y alcanzar condiciones en las que pueda darse combustión. La combustión provoca cargas térmicas y de presión que pueden dañar los sistemas de seguridad y la integridad del edificio de contención, última barrera de confinamiento de los materiales radiactivos. La principal condición que define las características de la combustión es la concentración de especies, por lo que el conocimiento detallado de la distribución de hidrógeno resulta muy importante para predecir correctamente los posibles daños en la contención en el caso de que se produjera combustión.

Simulación CFD de un accidente de caída de un elemento combustible en piscina

It is intended to provide a methodology of analysis more realistic this accident referred to in calculations of the license that requires fuel catastrophic break regardless of the height of the fall, with the consequent release of inventory analysers. Accidents that occurred in the past indicate that this hypothesis could be too conservative.