Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “COMPUTER GRAPHICS”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 721 records · Page 40

Advances in reactor physics education: Visualization of reactor parameters

Modern computer codes allow detailed neutron transport calculations. In combination with advanced 3D visualization software capable of treating large amounts of data in real time they form a powerful tool that can be used as a convenient modern educational tool for reactor operators, nuclear engineers, students and specialists involved in reactor operation and design. Visualization is applicable not only in education and training, but also as a tool for fuel management, core analysis and irradiation planning. The paper treats the visualization of neutron transport in different moderators, neutron flux and power distributions in two nuclear reactors (TRIGA type research reactor and a typical PWR). The distributions are calculated with MCNP and CORD-2 computer codes and presented using Amira software. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Computation of neutron fluxes in clusters of fuel pins arranged in hexagonal assemblies (2D and 3D)

For computations of fluxes, we have used Carvik's method of collision probabilities. This method requires tracking algorithms. An algorithm to compute tracks (in 2D and 3D) has been developed for seven hexagonal geometries with cluster of fuel pins. This has been implemented in the NXT module of the code DRAGON. The flux distribution in cluster of pins has been computed by using this code. For testing the results, they are compared when possible with the EXCELT module of the code DRAGON. Tracks are plotted in the NXT module by using MATLAB, these plots are also presented here. Results are presented with increasing number of lines to show the convergence of these results. We have numerically computed volumes, surface areas and the percentage errors in these computations. These results show that 2D results converge faster than 3D results. The accuracy on the computation of fluxes up to second decimal is achieved with fewer lines. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Modeling CANDU-6 liquid zone controllers for effects of thorium-based fuels

We use the DRAGON code to model the CANDU-6 liquid zone controllers and evaluate the effects of thorium-based fuels on their incremental cross sections and reactivity worth. We optimize both the numerical quadrature and spatial discretization for 2D cell models in order to provide accurate fuel properties for 3D liquid zone controller supercell models. We propose a low computer cost parameterized pseudo-exact 3D cluster geometries modeling approach that avoids tracking issues on small external surfaces. This methodology provides consistent incremental cross sections and reactivity worths when the thickness of the buffer region is reduced. When compared with an approximate annular geometry representation of the fuel and coolant region, we observe that the cluster description of fuel bundles in the supercell models does not increase considerably the precision of the results while increasing substantially the CPU time. In addition, this comparison shows that it is imperative to finely describe the liquid zone controller geometry since it has a strong impact of the incremental cross sections. This paper also shows that liquid zone controller reactivity worth is greatly decreased in presence of thorium-based fuels compared to the reference natural uranium fuel, since the fission and the fast to thermal scattering incremental cross sections are higher for the new fuels. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A coupled neutronic/thermal-hydraulic scheme between COBAYA3 and SUBCHANFLOW within the NURESIM simulation platform

Multi-scale, multi-physics problems reveal significant challenges while dealing with coupled neutronic/thermal-hydraulic solutions. Current generation of codes applied to Light Water Reactors (LWR) are based on 3D neutronic nodal methods coupled with one or two phase flow thermal-hydraulic system or sub-channel codes. In addition, spatial meshing and temporal schemes are crucial for the proper description of the non-symmetrical core behavior in case of transient and accidents e.g. reactivity insertion accidents. This paper describes the coupling approach between the 3D neutron diffusion code COBAYA3 and the sub-channel code SUBCHANFLOW within SALOME. The coupling is done inside the SALOME open source platform that is characterized by a powerful pre- and post-processing capabilities and a novel functionality for mapping of the neutronic and thermal hydraulic domains. The peculiar functionalities of SALOME and the steps required for the code integration and coupling are presented. The validation of the coupled codes is done based on two benchmarks the PWR MOX/UO{sub 2} RIA and the TMI-1 MSLB benchmark. A discussion of the prediction capability of COBAYA3/SUBCHANFLOW compared to other coupled solutions will be provided too. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Boron injection/dilution capabilities in TRACB/NEM coupled code

The coupled code TRAC-BF1/NEM is a thermal-hydraulic-neutronic code which allows transient simulations considering neutronic 3D and thermal-hydraulic process in multiple channels with one-dimensional geometry. TRAC-BF1 and NEM can be executed either in stand-alone mode, i.e. without coupling, as well as coupled. In stand-alone calculations NEM code is used without coupling and the thermal-hydraulic conditions (fuel temperature, moderator density and boron concentration) and xenon concentration for each node are taken from the SIMULATE3 output files. The NEM's source code has been modified to be able to read these conditions from external files when it is executed without being coupled. The coupling between TRAC-BF1 and NEM follows an integration scheme in which the thermal-hydraulic solution of TRAC-BF1 is sent to NEM to incorporate the feedback effects through the cross sections. TRAC-BF1 solves heat conduction equations inside of the heat structures using the 3D power distribution from NEM. The coupling is carried out through the communication protocol functions of PVM (Parallel Virtual Machine). The present article presents a study which constitutes an advance in the simulation of injection, transport and mix of boron in the reactor, increasing the capabilities of TRAC-BF1/NEM coupled code. This article shows the modifications introduced in the TRAC-BF1/NEM's source code to allow a more realistic simulation of boron injection transients. The qualification of these improvements in both codes is performed simulating a steady state of a generic BWR at nominal power. The results have been compared with SIMULATE3 which is used as a reference to obtain the cross sections through the SIMTAB methodology. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Development and preliminary verification of the 3D core neutronic code: COCO

As the recent blooming economic growth and following environmental concerns (China)) is proactively pushing forward nuclear power development and encouraging the tapping of clean energy. Under this situation, CGNPC, as one of the largest energy enterprises in China, is planning to develop its own nuclear related technology in order to support more and more nuclear plants either under construction or being operation. This paper introduces the recent progress in software development for CGNPC. The focus is placed on the physical models and preliminary verification results during the recent development of the 3D Core Neutronic Code: COCO. In the COCO code, the non-linear Green's function method is employed to calculate the neutron flux. In order to use the discontinuity factor, the Neumann (second kind) boundary condition is utilized in the Green's function nodal method. Additionally, the COCO code also includes the necessary physical models, e.g. single-channel thermal-hydraulic module, burnup module, pin power reconstruction module and cross-section interpolation module. The preliminary verification result shows that the COCO code is sufficient for reactor core design and analysis for pressurized water reactor (PWR). (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

CASMO5/TSUNAMI-3D spent nuclear fuel reactivity uncertainty analysis

The CASMO5 lattice physics code is used in conjunction with the TSUNAMI-3D sequence in ORNL's SCALE 6 code system to estimate the uncertainties in hot-to-cold reactivity changes due to cross-section uncertainty for PWR assemblies at various burnup points. The goal of the analysis is to establish the multiplication factor uncertainty similarity between various fuel assemblies at different conditions in a quantifiable manner and to obtain a bound on the hot-to-cold reactivity uncertainty over the various assembly types and burnup attributed to fundamental cross-section data uncertainty. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

3D analysis of the reactivity insertion accident in VVER-1000

Fuel parameters such as peak enthalpy and temperature during rod ejection accident are calculated. The calculations are performed by 3D neutron kinetics code NESTLE and 3D thermal-hydraulic code VIPRE-W. Both hot zero power and hot full power cases were studied for an equilibrium cycle with Westinghouse hex fuel in VVER-1000. It is shown that the use of 3D methodology can significantly increase safety margins for current criteria and met future criteria. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Identification of limiting case between DBA and SBDBA (CL break area sensitivity): A new model for the boron injection system

Atucha-2 is a Siemens-designed PHWR reactor under construction in the Republic of Argentina. Its geometrical complexity and (e.g., oblique Control Rods, Positive Void coefficient) required a developed and validated complex three dimensional (3D) neutron kinetics (NK) coupled thermal hydraulic (TH) model. Reactor shut-down is obtained by oblique CRs and, during accidental conditions, by an emergency shut-down system (JDJ) injecting a highly concentrated boron solution (boron clouds) in the moderator tank, the boron clouds reconstruction is obtained using a CFD (CFX) code calculation. A complete LBLOCA calculation implies the application of the RELAP5-3D{sup C} system code. Within the framework of the third Agreement 'NA-SA - Univ. of Pisa' a new RELAP5-3D control system for the boron injection system was developed and implemented in the validated coupled RELAP5-3D/NESTLE model of the Atucha 2 NPP. The aim of this activity is to find out the limiting case (maximum break area size) for the Peak Cladding Temperature for LOCAs under fixed boundary conditions. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Full core reactor analysis: Running Denovo on Jaguar

Fully-consistent, full-core, 3D, deterministic neutron transport simulations using the orthogonal mesh code Denovo were run on the massively parallel computing architecture Jaguar XT5. Using energy and spatial parallelization schemes, Denovo was able to efficiently scale to more than 160 k processors. Cell-homogenized cross sections were used with step-characteristics, linear-discontinuous finite element, and trilinear-discontinuous finite element spatial methods. It was determined that using the finite element methods gave considerably more accurate eigenvalue solutions for large-aspect ratio meshes than using step-characteristics. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Hybrid parallel code acceleration methods in full-core reactor physics calculations

When dealing with nuclear reactor calculation schemes, the need for three dimensional (3D) transport-based reference solutions is essential for both validation and optimization purposes. Considering a benchmark problem, this work investigates the potential of discrete ordinates (Sn) transport methods applied to 3D pressurized water reactor (PWR) full-core calculations. First, the benchmark problem is described. It involves a pin-by-pin description of a 3D PWR first core, and uses a 8-group cross-section library prepared with the DRAGON cell code. Then, a convergence analysis is performed using the PENTRAN parallel Sn Cartesian code. It discusses the spatial refinement and the associated angular quadrature required to properly describe the problem physics. It also shows that initializing the Sn solution with the EDF SPN solver COCAGNE reduces the number of iterations required to converge by nearly a factor of 6. Using a best estimate model, PENTRAN results are then compared to multigroup Monte Carlo results obtained with the MCNP5 code. Good consistency is observed between the two methods (Sn and Monte Carlo), with discrepancies that are less than 25 pcm for the k{sub eff}, and less than 2.1% and 1.6% for the flux at the pin-cell level and for the pin-power distribution, respectively. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

ABWR start-up test analysis using BWR core simulator with three-dimensional direct response matrix method

The ABWR start-up test analysis has been done with the BWR core simulator using the three--dimensional direct response matrix (3D-DRM) method. The Monte Carlo code VMONT made the sub-response matrices for the 3D-DRM method. Each boundary surface was subdivided by 4 x 4 for transverse segments, by 4 for angular segments and by 4 for axial zones in a node. For the calculation speedup, the 3D-DRM code used the divided sub-response matrices data set. The code used the MPI and OpenMP for the parallelized method. The median value is set as the average critical eigenvalues. The changes from the maximum value to the minimum value are 0.34 %{Delta}k with the spectral history method and 0.40 %{Delta}k without it, and the respective standard deviations were 0.12 % and 0.14 %. Using the spectral history method decreased the variation by 0.06 %{Delta}k. The root mean square differences of the axial power distribution were about 6 % between the analysis results and the plant data. Using the currents which converged in the previous exposure step reduced the number of iterations when the CR pattern changed only slightly. The averaged calculation time for each exposure step was about 5 hours on 12 PC Linux cluster servers with Core 2 Quad 3 GHz. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Benchmark of Atucha-2 PHWR RELAP5-3D control rod model by Monte Carlo MCNP5 core calculation

Atucha-2 is a Siemens-designed PHWR reactor under construction in the Republic of Argentina. Its geometrical complexity and peculiarities require the adoption of advanced Monte Carlo codes for performing realistic neutronic simulations. Therefore core models of Atucha-2 PHWR were developed using MCNP5. In this work a methodology was set up to collect the flux in the hexagonal mesh by which the Atucha-2 core is represented. The scope of this activity is to evaluate the effect of obliquely inserted control rod on neutron flux in order to validate the RELAP5-3D{sup C}/NESTLE three dimensional neutron kinetic coupled thermal-hydraulic model, applied by GRNSPG/UNIPI for performing selected transients of Chapter 15 FSAR of Atucha-2. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

An example of neutronic penalizations in reactivity transient analysis using 3D coupled chain HEMERA

HEMERA (Highly Evolutionary Methods for Extensive Reactor Analyses), is a fully coupled 3D computational chain developed jointly by IRSN and CEA. It is composed of CRONOS2 (core neutronics, cross sections library from APOLLO2), FLICA4 (core thermal-hydraulics) and the system code CATHARE. Multi-level and multi-dimensional models are developed to account for neutronics, core thermal-hydraulics, fuel thermal analysis and system thermal-hydraulics, dedicated to best-estimate, conservative simulations and sensitivity analysis. In IRSN, the HEMERA chain is widely used to study several types of reactivity accidents and for sensitivity studies. Just as an example of the HEMERA possibilities, we present here two types of neutronic penalizations and their impact on a power transient due to a REA (Rod Ejection Accident): in the first one, we studied a bum-up distribution modification and in the second one, a delayed-neutron fraction modification. Both modifications are applied to the whole core or localized in a few assemblies. Results show that it is possible to use global or local changes but 1) in case of bum-up modification, the total core power can increase when assembly peak power decrease so, care has to be taken if the goal is to maximize a local power peak and 2) for delayed-neutron fraction, a local modification can have the same effect as the one on the whole core, provided that it is large enough. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Interpretation of 3D void measurements with Tripoli4.6/JEFF3.1.1 Monte Carlo code

The present work details the first analysis of the 3D void phase conducted during the EPICURE/UM17x17/7% mixed UOX/MOX configuration. This configuration is composed of a homogeneous central 17x17 MOX-7% assembly, surrounded by portions of 17x17 1102 assemblies with guide-tubes. The void bubble is modelled by a small waterproof 5x5 fuel pin parallelepiped box of 11 cm height, placed in the centre of the MOX assembly. This bubble, initially placed at the core mid-plane, is then moved in different axial positions to study the evolution in the core of the axial perturbation. Then, to simulate the growing of this bubble in order to understand the effects of increased void fraction along the fuel pin, 3 and 5 bubbles have been stacked axially, from the core mid-plane. The C/E comparison obtained with the Monte Carlo code Tripoli4 for both radial and axial fission rate distributions, and in particular the reproduction of the very important flux gradients at the void/water interfaces, changing as the bubble is displaced along the z-axis are very satisfactory. It demonstrates both the capability of the code and its library to reproduce this kind of situation, as the very good quality of the experimental results, confirming the UM-17x17 as an excellent experimental benchmark for 3D code validation. This work has been performed within the frame of the V and V program for the future APOLL03 deterministic code of CEA starting in 2012, and its V and V benchmarking database. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Development of 3D full-core ERANOS-2.2/MCNPX-2.7.0 models and neutronic analysis of the BFS-2 zero-power facility

The present paper is addressing the development and validation against experimental data of 3D full-core models of the BFS-2 zero-power fast-reactor using both the deterministic system code ERANOS-2.2 and the stochastic code MCNPX-2.7.0. The model configuration of BFS considered for analysis is the BFS-62-3A benchmark. To extend the - deterministic/stochastic - code-to-code comparison, neutronic parameters, i.e. reactivity, neutron spectrum and reaction rates, were also simulated at the cell level with the Monte Carlo code SERPENT-1.1.7 with two modern data libraries, ENDF-B/VII and JEFF-3.1.1. The BFS-2 critical zero-power facility at the Inst. of Physics and Power Engineering (IPPE) was designed for simulations of the core and shielding of sodium-cooled, fast reactors, for neutron data validation and comparison with experimental results. At the BFS-2 facility, the BFS-62-3A critical benchmark experiment was set-up as a mock-up of the BN-600 reactor, with hybrid MOX fuel and stainless steel reflectors. A UO{sub 2} blanket and a large non-homogeneous stainless-steel reflector surround the core. The lattice is hexagonal of pitch 5.1 cm and metallic dowels are used to keep in central position cylindrical rods made of different types of material (fissile, fertile, blanket, plenum, shielding and absorber). A typical subassembly is formed in piling up various pellets of about 1 cm in height and 4.6 cm in diameter, conferring large heterogeneity in the axial direction. The full-core model development was a complex task due to the large number of subassemblies and the axial subassembly heterogeneity. In ERANOS-2.2, it was necessary to homogenize axially per region the pellets used to form the subassembly. The self-shielded macroscopic cross-sections were calculated using the cell code ECCO in association with JEFF-3.1 and ENDF/B-VI.8 data libraries. The core calculations were performed with broad cross-sections data in 33 neutron energy groups with the solver AVNM in the diffusion approximation, mostly. In MCNPX-2.7.0, a step-by-step approach was used, starting with a model in which the fissile rods were simulated on a homogeneous level, to finally integrate the actual heterogeneous description of the subassemblies. The code-to-code cell analysis performed between ECCO, SERPENT and MCNPX with different modern nuclear data library revealed that the results for the infinite multiplication factor between Monte Carlo and deterministic analysis are in good agreement ({Delta}p < 100 pcm). The differences between the results were observed to be larger for the neutron data libraries, with reactivity differences up to 350 pcm. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Application of the JENDL-4.0 nuclear data set for uncertainty analysis of the prototype FBR Monju

This paper deals with uncertainty analysis of the Monju reactor using JENDL-4.0 and the ERANOS code 1. In 2010 the Japan Atomic Energy Agency - JAEA - released the JENDL-4.0 nuclear data set. This new evaluation contains improved values of cross-sections and emphasizes accurate covariance matrices. Also in 2010, JAEA restarted the sodium-cooled fast reactor prototype Monju after about 15 years of shutdown. The long shutdown time resulted in a build-up of {sup 241}Am by natural decay from the initially loaded Pu. As well as improved covariance matrices, JENDL-4.0 is announced to contain improved data for minor actinides 2. The choice of Monju reactor as an application of the new evaluation seems then even more relevant. The uncertainty analysis requires the determination of sensitivity coefficients. The well-established ERANOS code was chosen because of its integrated modules that allow users to perform sensitivity and uncertainty analysis. A JENDL-4.0 cross-sections library is not available for ERANOS. Therefor a cross-sections library had to be made from the original ENDF files for the ECCO cell code (part of ERANOS). For confirmation of the newly made library, calculations of a benchmark core were performed. These calculations used the MZA and MZB benchmarks and showed consistent results with other libraries. Calculations for the Monju reactor were performed using hexagonal 3D geometry and PN transport theory. However, the ERANOS sensitivity modules cannot use the resulting fluxes, as these modules require finite differences based fluxes, obtained from RZ SN-transport or 3D diffusion calculations. The corresponding geometrical models have been made and the results verified with Monju restart experimental data 4. Uncertainty analysis was performed using the RZ model. JENDL-4.0 uncertainty analysis showed a significant reduction of the uncertainty related to the fission cross-section of Pu along with an increase of the uncertainty related to the capture cross-section of {sup 238}U compared with the previous JENDL-3.3 version. Covariance data recently added in JENDL-4.0 for {sup 241}Am appears to have a non-negligible contribution. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

The developments and verifications of trace model for IIST LOCA experiments

The test facility IIST (INER Integral System Test) is a Reduced-Height and Reduced-Pressure (RHRP) integral test loop, which was constructed for the purposes of conducting thermal hydraulic and safety analysis of the Westinghouse three-loop PWR Nuclear Power Plants. The main purpose of this study is to develop and verify TRACE models of IIST through the IIST small break loss of coolant accident (SBLOCA) experiments. First, two different IIST TRACE models which include a pipe-vessel model and a 3-D vessel component model have been built. The steady state and transient calculation results show that both TRACE models have the ability to simulate the related IIST experiments. Comparing with IIST SBLOCA experiment data, the 3-D vessel component model has shown better simulation capabilities so that it has been chosen for all further thermal hydraulic studies. The second step is the sensitivity studies of two phase multiplier and subcooled liquid multiplier in choked flow model; and two correlation constants in CCFL model respectively. As a result, an appropriate set of multipliers and constants can be determined. In summary, a verified IIST TRACE model with 3D vessel component, and fine-tuned choked flow model and CCFL model is established for further studies on IIST experiments in the future. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗