Interface requirements to couple thermal-hydraulic codes to 3D neutronic codes
Langenbuch, S.; Austregesilo, H.; Velkov, K.
1997-07-01
The present situation of thermalhydraulics codes and 3D neutronics codes is briefly described and general considerations for coupling of these codes are discussed. Two different basic approaches of coupling are identified and their relative advantages and disadvantages are discussed. The implementation of the coupling for 3D neutronics codes in the system ATHLET is presented. Meanwhile, this interface is used for coupling three different 3D neutronics codes.
3D neutronic codes coupled with thermal-hydraulic system codes for PWR, and BWR and VVER reactors
Langenbuch, S.; Velkov, K.; Lizorkin, M.
1997-07-01
This paper describes the objectives of code development for coupling 3D neutronics codes with thermal-hydraulic system codes. The present status of coupling ATHLET with three 3D neutronics codes for VVER- and LWR-reactors is presented. After describing the basic features of the 3D neutronic codes BIPR-8 from Kurchatov-Institute, DYN3D from Research Center Rossendorf and QUABOX/CUBBOX from GRS, first applications of coupled codes for different transient and accident scenarios are presented. The need of further investigations is discussed.
TART97 a coupled neutron-photon 3-D, combinatorial geometry Monte Carlo transport code
Cullen, D.E.
1997-11-22
TART97 is a coupled neutron-photon, 3 Dimensional, combinatorial geometry, time dependent Monte Carlo transport code. This code can on any modern computer. It is a complete system to assist you with input preparation, running Monte Carlo calculations, and analysis of output results. TART97 is also incredibly FAST; if you have used similar codes, you will be amazed at how fast this code is compared to other similar codes. Use of the entire system can save you a great deal of time and energy. TART97 is distributed on CD. This CD contains on- line documentation for all codes included in the system, the codes configured to run on a variety of computers, and many example problems that you can use to familiarize yourself with the system. TART97 completely supersedes all older versions of TART, and it is strongly recommended that users only use the most recent version of TART97 and its data riles.
Podlazov, L. N.
1998-07-29
Two realistic benchmark problems are defined and used to assess the performance of coupled thermal-hydraulic and neutronic codes used in simulating dynamic processes in VVER-1000 and RBMK reactor systems. One of the problems simulates a design basis accident involving the ejection of three control and protection system rods from a VVER-1000 reactor. The other is based on a postulated rod withdrawal from an operating RBMK reactor. Preliminary results calculated by various codes are compared. While these results show significant differences, the intercomparisons performed so far provide a basis for further evaluation of code limitations and modeling assumptions.
Parameterization of nuclear cross-sections for coupled neutronic- thermalhydraulic codes
Miro, R.; Verdu, G.; Barrachina, T.; Rosello, O.
2006-07-01
The present work consists of developing an in-house methodology, called SIMTAB, to characterize, in a simplified way, the reactor core of LWR Nuclear Power Plants. Specifically, a cross-sections and kinetic parameters set are obtained as a function of the prompt and control variables. So that, the core can be modeled using a limited number of neutronic regions, in such a way that the reactor kinetic behavior is properly characterized. This simplification of the reactor core permits, from an operative point of view, the use of few cross sections data sets in coupled 3D neutronic-thermalhydraulic codes. (authors)
Energy Science and Technology Software Center (ESTSC)
2013-06-24
Version 07 TART2012 is a coupled neutron-photon Monte Carlo transport code designed to use three-dimensional (3-D) combinatorial geometry. Neutron and/or photon sources as well as neutron induced photon production can be tracked. It is a complete system to assist you with input preparation, running Monte Carlo calculations, and analysis of output results. TART2012 is also incredibly FAST; if you have used similar codes, you will be amazed at how fast this code is compared tomore » other similar codes. Use of the entire system can save you a great deal of time and energy. TART2012 extends the general utility of the code to even more areas of application than available in previous releases by concentrating on improving the physics, particularly with regard to improved treatment of neutron fission, resonance self-shielding, molecular binding, and extending input options used by the code. Several utilities are included for creating input files and displaying TART results and data. TART2012 uses the latest ENDF/B-VI, Release 8, data. New for TART2012 is the use of continuous energy neutron cross sections, in addition to its traditional multigroup cross sections. For neutron interaction, the data are derived using ENDF-ENDL2005 and include both continuous energy cross sections and 700 group neutron data derived using a combination of ENDF/B-VI, Release 8, and ENDL data. The 700 group structure extends from 10-5 eV up to 1 GeV. Presently nuclear data are only available up to 20 MeV, so that only 616 of the groups are currently used. For photon interaction, 701 point photon data were derived using the Livermore EPDL97 file. The new 701 point structure extends from 100 eV up to 1 GeV, and is currently used over this entire energy range. TART2012 completely supersedes all older versions of TART, and it is strongly recommended that one use only the most recent version of TART2012 and its data files. Check authors homepage for related information: http
A Coupled Neutron-Photon 3-D Combinatorial Geometry Monte Carlo Transport Code
Energy Science and Technology Software Center (ESTSC)
1998-06-12
TART97 is a coupled neutron-photon, 3 dimensional, combinatorial geometry, time dependent Monte Carlo transport code. This code can run on any modern computer. It is a complete system to assist you with input preparation, running Monte Carlo calculations, and analysis of output results. TART97 is also incredibly fast: if you have used similar codes, you will be amazed at how fast this code is compared to other similar codes. Use of the entire system canmore » save you a great deal of time and energy. TART 97 is distributed on CD. This CD contains on-line documentation for all codes included in the system, the codes configured to run on a variety of computers, and many example problems that you can use to familiarize yourself with the system. TART97 completely supersedes all older versions of TART, and it is strongly recommended that users only use the most recent version of TART97 and ist data files.« less
Cynod: A Neutronics Code for Pebble Bed Modular Reactor Coupled Transient Analysis
Hikaru Hiruta; Abderrafi M. Ougouag; Hans D. Gougar; Javier Ortensi
2008-09-01
The Pebble Bed Reactor (PBR) is one of the two concepts currently considered for development into the Next Generation Nuclear Plant (NGNP). This interest is due, in particular, to the concept’s inherent safety characteristics. In order to verify and confirm the design safety characteristics of the PBR computational tools must be developed that treat the range of phenomena that are expected to be important for this type of reactors. This paper presents a recently developed 2D R-Z cylindrical nodal kinetics code and shows some of its capabilities by applying it to a set of known and relevant benchmarks. The new code has been coupled to the thermal hydraulics code THERMIX/KONVEK[1] for application to the simulation of very fast transients in PBRs. The new code, CYNOD, has been written starting with a fixed source solver extracted from the nodal cylindrical geometry solver contained within the PEBBED code. The fixed source solver was then incorporated into a kinetic solver.. The new code inherits the spatial solver characteristics of the nodal solver within PEBBED. Thus, the time-dependent neutron diffusion equation expressed analytically in each node of the R-Z cylindrical geometry sub-domain (or node) is transformed into one-dimensional equations by means of the usual transverse integration procedure. The one-dimensional diffusion equations in each of the directions are then solved using the analytic Green’s function method. The resulting equations for the entire domain are then re-cast in the form of the Direct Coarse Mesh Finite Difference (D-CMFD) for convenience of solution. The implicit Euler method is used for the time variable discretization. In order to correctly treat the cusping effect for nodes that contain a partially inserted control rod a method is used that takes advantage of the Green’s function solution available in the intrinsic method. In this corrected treatment, the nodes are re-homogenized using axial flux shapes reconstructed based on the
Petruzzi, A.; D'Auria, F.; Giannotti, W.; Ivanov, K.
2005-02-15
The best-estimate calculation results from complex system codes are affected by approximations that are unpredictable without the use of computational tools that account for the various sources of uncertainty.The code with (the capability of) internal assessment of uncertainty (CIAU) has been previously proposed by the University of Pisa to realize the integration between a qualified system code and an uncertainty methodology and to supply proper uncertainty bands each time a nuclear power plant (NPP) transient scenario is calculated. The derivation of the methodology and the results achieved by the use of CIAU are discussed to demonstrate the main features and capabilities of the method.In a joint effort between the University of Pisa and The Pennsylvania State University, the CIAU method has been recently extended to evaluate the uncertainty of coupled three-dimensional neutronics/thermal-hydraulics calculations. The result is CIAU-TN. The feasibility of the approach has been demonstrated, and sample results related to the turbine trip transient in the Peach Bottom NPP are shown. Notwithstanding that the full implementation and use of the procedure requires a database of errors not available at the moment, the results give an idea of the errors expected from the present computational tools.
Cullen, D E
1998-11-22
TART98 is a coupled neutron-photon, 3 Dimensional, combinatorial geometry, time dependent Monte Carlo radiation transport code. This code can run on any modern computer. It is a complete system to assist you with input preparation, running Monte Carlo calculations, and analysis of output results. TART98 is also incredibly FAST; if you have used similar codes, you will be amazed at how fast this code is compared to other similar codes. Use of the entire system can save you a great deal of time and energy. TART98 is distributed on CD. This CD contains on-line documentation for all codes included in the system, the codes configured to run on a variety of computers, and many example problems that you can use to familiarize yourself with the system. TART98 completely supersedes all older versions of TART, and it is strongly recommended that users only use the most recent version of TART98 and its data files.
Code System for 2-Group, 3D Neutronic Kinetics Calculations Coupled to Core Thermal Hydraulics.
Energy Science and Technology Software Center (ESTSC)
2000-05-12
Version 00 QUARK is a combined computer program comprising a revised version of the QUANDRY three-dimensional, two-group neutron kinetics code and an upgraded version of the COBRA transient core analysis code (COBRA-EN). Starting from either a critical steady-state (k-effective or critical dilute Boron problem) or a subcritical steady-state (fixed source problem) in a PWR plant, the code allows one to simulate the neutronic and thermal-hydraulic core transient response to reactivity accidents initiated both inside themore » vessel (such as a control rod ejection) and outside the vessel (such as the sudden change of the Boron concentration in the coolant). QUARK output can be used as input to PSR-470/NORMA-FP to perform a subchannel analysis from converged coarse-mesh nodal solutions.« less
Leap Frog and Time Step Sub-Cycle Scheme for Coupled Neutronics and Thermal-Hydraulic Codes
Lu, S.
2002-07-01
As the result of the advancing TCP/IP based inter-process communication technology, more and more legacy thermal-hydraulic codes have been coupled with neutronics codes to provide best-estimate capabilities for reactivity related reactor transient analysis. Most of the coupling schemes are based on closely coupled serial or parallel approaches. Therefore, the execution of the coupled codes usually requires significant CPU time, when a complicated system is analyzed. Leap Frog scheme has been used to reduce the run time. The extent of the decoupling is usually determined based on a trial and error process for a specific analysis. It is the intent of this paper to develop a set of general criteria, which can be used to invoke the automatic Leap Frog algorithm. The algorithm will not only provide the run time reduction but also preserve the accuracy. The criteria will also serve as the base of an automatic time step sub-cycle scheme when a sudden reactivity change is introduced and the thermal-hydraulic code is marching with a relatively large time step. (authors)
Bingham, Philip R; Santos-Villalobos, Hector J
2011-01-01
Coded aperture techniques have been applied to neutron radiography to address limitations in neutron flux and resolution of neutron detectors in a system labeled coded source imaging (CSI). By coding the neutron source, a magnified imaging system is designed with small spot size aperture holes (10 and 100 m) for improved resolution beyond the detector limits and with many holes in the aperture (50% open) to account for flux losses due to the small pinhole size. An introduction to neutron radiography and coded aperture imaging is presented. A system design is developed for a CSI system with a development of equations for limitations on the system based on the coded image requirements and the neutron source characteristics of size and divergence. Simulation has been applied to the design using McStas to provide qualitative measures of performance with simulations of pinhole array objects followed by a quantitative measure through simulation of a tilted edge and calculation of the modulation transfer function (MTF) from the line spread function. MTF results for both 100um and 10um aperture hole diameters show resolutions matching the hole diameters.
Curca-Tivig, Florin; Merk, Stephan; Pautz, Andreas; Thareau, Sebastien
2007-07-01
Anticipating future needs of our customers and willing to concentrate synergies and competences existing in the company for the benefit of our customers, AREVA NP decided in 2002 to develop the next generation of coupled neutronics/ core thermal-hydraulic (TH) code systems for fuel assembly and core design calculations for both, PWR and BWR applications. The global CONVERGENCE project was born: after a feasibility study of one year (2002) and a conceptual phase of another year (2003), development was started at the beginning of 2004. The present paper introduces the CONVERGENCE project, presents the main feature of the new code system ARCADIA{sup R} and concludes on customer benefits. ARCADIA{sup R} is designed to meet AREVA NP market and customers' requirements worldwide. Besides state-of-the-art physical modeling, numerical performance and industrial functionality, the ARCADIA{sup R} system is featuring state-of-the-art software engineering. The new code system will bring a series of benefits for our customers: e.g. improved accuracy for heterogeneous cores (MOX/ UOX, Gd...), better description of nuclide chains, and access to local neutronics/ thermal-hydraulics and possibly thermal-mechanical information (3D pin by pin full core modeling). ARCADIA is a registered trademark of AREVA NP. (authors)
Lee, Deokjung; Downar, Thomas J.; Ulses, Anthony; Akdeniz, Bedirhan; Ivanov, Kostadin N.
2004-10-15
An analysis of the Peach Bottom Unit 2 Turbine Trip 2 (TT2) experiment has been performed using the U.S. Nuclear Regulatory Commission coupled thermal-hydraulics and neutronics code TRAC-M/PARCS. The objective of the analysis was to assess the performance of TRAC-M/PARCS on a BWR transient with significance in two-phase flow and spatial variations of the neutron flux. TRAC-M/PARCS results are found to be in good agreement with measured plant data for both steady-state and transient phases of the benchmark. Additional analyses of four fictitious extreme scenarios are performed to provide a basis for code-to-code comparisons and comprehensive testing of the thermal-hydraulics/neutronics coupling. The obtained results of sensitivity studies on the effect of direct moderator heating on transient simulation indicate the importance of this modeling aspect.
Parallelization of the MAAP-A code neutronics/thermal hydraulics coupling
Froehle, P.H.; Wei, T.Y.C.; Weber, D.P.; Henry, R.E.
1998-12-31
A major new feature, one-dimensional space-time kinetics, has been added to a developmental version of the MAAP code through the introduction of the DIF3D-K module. This code is referred to as MAAP-A. To reduce the overall job time required, a capability has been provided to run the MAAP-A code in parallel. The parallel version of MAAP-A utilizes two machines running in parallel, with the DIF3D-K module executing on one machine and the rest of the MAAP-A code executing on the other machine. Timing results obtained during the development of the capability indicate that reductions in time of 30--40% are possible. The parallel version can be run on two SPARC 20 (SUN OS 5.5) workstations connected through the ethernet. MPI (Message Passing Interface standard) needs to be implemented on the machines. If necessary the parallel version can also be run on only one machine. The results obtained running in this one-machine mode identically match the results obtained from the serial version of the code.
Russell, G.J.; Pitcher, E.J.; Ferguson, P.D.
1995-12-01
Optimizing the neutronic performance of a coupled-moderator system for a Long-Pulse Spallation Source is a new and challenging area for the spallation target-system designer. For optimal performance of a neutron source, it is essential to have good communication with instrument scientists to obtain proper design criteria and continued interaction with mechanical, thermal-hydraulic, and materials engineers to attain a practical design. A good comprehension of the basics of coupled-moderator neutronics will aid in the proper design of a target system for a Long-Pulse Spallation Source.
Horiguchi, Hironori; Sato, Tatsuhiko; Kumada, Hiroaki; Yamamoto, Tetsuya; Sakae, Takeji
2015-01-01
The absorbed doses deposited by boron neutron capture therapy (BNCT) can be categorized into four components: α and 7Li particles from the 10B(n, α)7Li reaction, 0.54-MeV protons from the 14N(n, p)14C reaction, the recoiled protons from the 1H(n, n) 1H reaction, and photons from the neutron beam and 1H(n, γ)2H reaction. For evaluating the irradiation effect in tumors and the surrounding normal tissues in BNCT, it is of great importance to estimate the relative biological effectiveness (RBE) for each dose component in the same framework. We have, therefore, established a new method for estimating the RBE of all BNCT dose components on the basis of the microdosimetric kinetic model. This method employs the probability density of lineal energy, y, in a subcellular structure as the index for expressing RBE, which can be calculated using the microdosimetric function implemented in the particle transport simulation code (PHITS). The accuracy of this method was tested by comparing the calculated RBE values with corresponding measured data in a water phantom irradiated with an epithermal neutron beam. The calculation technique developed in this study will be useful for biological dose estimation in treatment planning for BNCT. PMID:25428243
NASA Astrophysics Data System (ADS)
Ichikawa, Ryoko; Masuhara, Yasuhiro; Kasahara, Fumio
The Best Estimate Plus Uncertainty (BEPU) method has been prepared for the regulatory cross-check analysis at Japan Nuclear Energy Safety Organization (JNES) on base of the three-dimensional neutron-kinetics/thermal- hydraulics coupled code SKETCH-INS/TRACE5.0. In the preparation, TRACE5.0 is verified against the large-scale thermal-hydraulic tests carried out with NUPEC facility. These tests were focused on the pressure drop of steam-liquid two phase flow and void fraction distribution. From the comparison of the experimental data with other codes (RELAP5/MOD3.3 and TRAC-BF1), TRACE5.0 was judged better than other codes. It was confirmed that TRACE5.0 has high reliability for thermal hydraulics behavior and are used as a best-estimate code for the statistical safety evaluation. Next, the coupled code SKETCH-INS/TRACE5.0 was applied to turbine trip tests performed at the Peach Bottom-2 BWR4 Plant. The turbine trip event shows the rapid power peak due to the voids collapse with the pressure increase. The analyzed peak value of core power is better simulated than the previous version SKETCH-INS/TRAC-BF1. And the statistical safety evaluation using SKETCH-INS/TRACE5.0 was applied to the loss of load transient for examining the influence of the choice of sampling method.
Coupled Neutron Transport for HZETRN
NASA Technical Reports Server (NTRS)
Slaba, Tony C.; Blattnig, Steve R.
2009-01-01
Exposure estimates inside space vehicles, surface habitats, and high altitude aircrafts exposed to space radiation are highly influenced by secondary neutron production. The deterministic transport code HZETRN has been identified as a reliable and efficient tool for such studies, but improvements to the underlying transport models and numerical methods are still necessary. In this paper, the forward-backward (FB) and directionally coupled forward-backward (DC) neutron transport models are derived, numerical methods for the FB model are reviewed, and a computationally efficient numerical solution is presented for the DC model. Both models are compared to the Monte Carlo codes HETC-HEDS, FLUKA, and MCNPX, and the DC model is shown to agree closely with the Monte Carlo results. Finally, it is found in the development of either model that the decoupling of low energy neutrons from the light particle transport procedure adversely affects low energy light ion fluence spectra and exposure quantities. A first order correction is presented to resolve the problem, and it is shown to be both accurate and efficient.
Morgan C. White
2000-07-01
The fundamental motivation for the research presented in this dissertation was the need to development a more accurate prediction method for characterization of mixed radiation fields around medical electron accelerators (MEAs). Specifically, a model is developed for simulation of neutron and other particle production from photonuclear reactions and incorporated in the Monte Carlo N-Particle (MCNP) radiation transport code. This extension of the capability within the MCNP code provides for the more accurate assessment of the mixed radiation fields. The Nuclear Theory and Applications group of the Los Alamos National Laboratory has recently provided first-of-a-kind evaluated photonuclear data for a select group of isotopes. These data provide the reaction probabilities as functions of incident photon energy with angular and energy distribution information for all reaction products. The availability of these data is the cornerstone of the new methodology for state-of-the-art mutually coupled photon-neutron transport simulations. The dissertation includes details of the model development and implementation necessary to use the new photonuclear data within MCNP simulations. A new data format has been developed to include tabular photonuclear data. Data are processed from the Evaluated Nuclear Data Format (ENDF) to the new class ''u'' A Compact ENDF (ACE) format using a standalone processing code. MCNP modifications have been completed to enable Monte Carlo sampling of photonuclear reactions. Note that both neutron and gamma production are included in the present model. The new capability has been subjected to extensive verification and validation (V&V) testing. Verification testing has established the expected basic functionality. Two validation projects were undertaken. First, comparisons were made to benchmark data from literature. These calculations demonstrate the accuracy of the new data and transport routines to better than 25 percent. Second, the ability to
Iterative Reconstruction of Coded Source Neutron Radiographs
Santos-Villalobos, Hector J; Bingham, Philip R; Gregor, Jens
2012-01-01
Use of a coded source facilitates high-resolution neutron imaging but requires that the radiographic data be deconvolved. In this paper, we compare direct deconvolution with two different iterative algorithms, namely, one based on direct deconvolution embedded in an MLE-like framework and one based on a geometric model of the neutron beam and a least squares formulation of the inverse imaging problem.
Miro, R.; Maggini, F.; Barrachina, T.; Verdu, G.; Gomez, A.; Ortego, A.; Murillo, J. C.
2006-07-01
The Reactor Ejection Accident (REA) belongs to the Reactor Initiated Accidents (RIA) category of accidents and it is part of the licensing basis accident analyses required for pressure water reactors (PWR). The REA at hot zero power (HZP) is characterized by a single rod ejection from a core position with a very low power level. The evolution consists basically of a continuous reactivity insertion. The main feature limiting the consequences of the accident in a PWR is the Doppler Effect. To check the performance of the coupled code RELAP5/PARCS2.5 and RELAP5/VALKIN a REA in Trillo NPP is simulated. These analyses will allow knowing more accurately the PWR real plant phenomenology in the RIA most limiting conditions. (authors)
Secondary neutron source modelling using MCNPX and ALEPH codes
NASA Astrophysics Data System (ADS)
Trakas, Christos; Kerkar, Nordine
2014-06-01
Monitoring the subcritical state and divergence of reactors requires the presence of neutron sources. But mainly secondary neutrons from these sources feed the ex-core detectors (SRD, Source Range Detector) whose counting rate is correlated with the level of the subcriticality of reactor. In cycle 1, primary neutrons are provided by sources activated outside of the reactor (e.g. Cf252); part of this source can be used for the divergence of cycle 2 (not systematic). A second family of neutron sources is used for the second cycle: the spontaneous neutrons of actinides produced after irradiation of fuel in the first cycle. Both families of sources are not sufficient to efficiently monitor the divergence of the second cycles and following ones, in most reactors. Secondary sources cluster (SSC) fulfil this role. In the present case, the SSC [Sb, Be], after activation in the first cycle (production of Sb124, unstable), produces in subsequent cycles a photo-neutron source by gamma (from Sb124)-neutron (on Be9) reaction. This paper presents the model of the process between irradiation in cycle 1 and cycle 2 results for SRD counting rate at the beginning of cycle 2, using the MCNPX code and the depletion chain ALEPH-V1 (coupling of MCNPX and ORIGEN codes). The results of this simulation are compared with two experimental results of the PWR 1450 MWe-N4 reactors. A good agreement is observed between these results and the simulations. The subcriticality of the reactors is about at -15,000 pcm. Discrepancies on the SRD counting rate between calculations and measurements are in the order of 10%, lower than the combined uncertainty of measurements and code simulation. This comparison validates the AREVA methodology, which allows having an SRD counting rate best-estimate for cycles 2 and next ones and optimizing the position of the SSC, depending on the geographic location of sources, main parameter for optimal monitoring of subcritical states.
Iterative Reconstruction of Coded Source Neutron Radiographs
Santos-Villalobos, Hector J; Bingham, Philip R; Gregor, Jens
2013-01-01
Use of a coded source facilitates high-resolution neutron imaging through magnifications but requires that the radiographic data be deconvolved. A comparison of direct deconvolution with two different iterative algorithms has been performed. One iterative algorithm is based on a maximum likelihood estimation (MLE)-like framework and the second is based on a geometric model of the neutron beam within a least squares formulation of the inverse imaging problem. Simulated data for both uniform and Gaussian shaped source distributions was used for testing to understand the impact of non-uniformities present in neutron beam distributions on the reconstructed images. Results indicate that the model based reconstruction method will match resolution and improve on contrast over convolution methods in the presence of non-uniform sources. Additionally, the model based iterative algorithm provides direct calculation of quantitative transmission values while the convolution based methods must be normalized base on known values.
Fast-neutron, coded-aperture imager
NASA Astrophysics Data System (ADS)
Woolf, Richard S.; Phlips, Bernard F.; Hutcheson, Anthony L.; Wulf, Eric A.
2015-06-01
This work discusses a large-scale, coded-aperture imager for fast neutrons, building off a proof-of concept instrument developed at the U.S. Naval Research Laboratory (NRL). The Space Science Division at the NRL has a heritage of developing large-scale, mobile systems, using coded-aperture imaging, for long-range γ-ray detection and localization. The fast-neutron, coded-aperture imaging instrument, designed for a mobile unit (20 ft. ISO container), consists of a 32-element array of 15 cm×15 cm×15 cm liquid scintillation detectors (EJ-309) mounted behind a 12×12 pseudorandom coded aperture. The elements of the aperture are composed of 15 cm×15 cm×10 cm blocks of high-density polyethylene (HDPE). The arrangement of the aperture elements produces a shadow pattern on the detector array behind the mask. By measuring of the number of neutron counts per masked and unmasked detector, and with knowledge of the mask pattern, a source image can be deconvolved to obtain a 2-d location. The number of neutrons per detector was obtained by processing the fast signal from each PMT in flash digitizing electronics. Digital pulse shape discrimination (PSD) was performed to filter out the fast-neutron signal from the γ background. The prototype instrument was tested at an indoor facility at the NRL with a 1.8-μCi and 13-μCi 252Cf neutron/γ source at three standoff distances of 9, 15 and 26 m (maximum allowed in the facility) over a 15-min integration time. The imaging and detection capabilities of the instrument were tested by moving the source in half- and one-pixel increments across the image plane. We show a representative sample of the results obtained at one-pixel increments for a standoff distance of 9 m. The 1.8-μCi source was not detected at the 26-m standoff. In order to increase the sensitivity of the instrument, we reduced the fastneutron background by shielding the top, sides and back of the detector array with 10-cm-thick HDPE. This shielding configuration led
The CoCoNuT code: from neutron star oscillations to supernova explosions
NASA Astrophysics Data System (ADS)
Cerdá-Durán, P.; Gabler, M.; Müller, E.; Font, J. A.; Stergioulas, N.; Obergaulinger, M.; Aloy, M. A.; DeBrye, N.; Cordero-Carrión, I.; Ibáñez, J. M.
2013-05-01
CoCoNuT is a numerical code, that evolves the General relativistic magneto-hydrodynamics equations coupled to the Einstein equations in the CFC approximation. Its main purpose is to simulate astrophysical scenarios in which strong gravity is important such as the collapse of massive stars and the evolution of neutron stars. I review recent results of the numerical code regarding neutron star oscillations and core collapse supernova and its observational consequences.
3D neutronic/thermal-hydraulic coupled analysis of MYRRHA
Vazquez, M.; Martin-Fuertes, F.
2012-07-01
The current tendency in multiphysics calculations applied to reactor physics is the use of already validated computer codes, coupled by means of an iterative approach. In this paper such an approach is explained concerning neutronics and thermal-hydraulics coupled analysis with MCNPX and COBRA-IV codes using a driver program and file exchange between codes. MCNPX provides the neutronic analysis of heterogeneous nuclear systems, both in critical and subcritical states, while COBRA-IV is a subchannel code that can be used for rod bundles or core thermal-hydraulics analysis. In our model, the MCNP temperature dependence of nuclear data is handled via pseudo-material approach, mixing pre-generated cross section data set to obtain the material with the desired cross section temperature. On the other hand, COBRA-IV has been updated to allow for the simulation of liquid metal cooled reactors. The coupled computational tool can be applied to any geometry and coolant, as it is the case of single fuel assembly, at pin-by-pin level, or full core simulation with the average pin of each fuel-assembly. The coupling tool has been applied to the critical core layout of the SCK-CEN MYRRHA concept, an experimental LBE cooled fast reactor presently in engineering design stage. (authors)
Energy Science and Technology Software Center (ESTSC)
1982-04-30
Version 00 As a part of the measurement and analysis plan for the Dosimetry Experiment at the "JOYO" experimental fast reactor, neutron flux spectral analysis is performed using the NEUPAC (Neutron Unfolding Code Package) code system. NEUPAC calculates the neutron flux spectra and other integral quantities from the activation data of the dosimeter foils.
Multigroup Complex Geometry Neutron Diffusion Code System.
Energy Science and Technology Software Center (ESTSC)
2002-12-18
Version 01 SNAP-3D is based on SNAP2 and is a one- two- or three-dimensional multigroup diffusion code system. It is primarily intended for neutron diffusion calculations, but it can also carry out gamma-ray calculations if the diffusion approximation is accurate enough. It is suitable for fast and thermal reactor core calculations and for shield calculations. SNAP-3D can solve the multi-group neutron diffusion equations using finite difference methods in (x,y,z), (r,theta,z), (TRI,z), (HEX,z) or (spherical) coordinates.more » The one-dimensional slab and cylindrical geometries and the two-dimensional (x,y), (r,z), (r,theta), (HEX) and (TRI) are all treated as simple special cases of three-dimensional geometries. Numerous reflective and periodic symmetry options are available and may be used to reduce the number of mesh points necessary to represent the system. Extrapolation lengths can be specified at internal and external boundaries. The problem classes are: 1) eigenvalue search for critical k-effective, 2) eigenvalue search for critical buckling, 3) eigenvalue search for critical time-constant, 4) fixed source problems in which the sources are functions of regions, 5) fixed source problems in which the sources are provided, on disc, for every mesh point and group.« less
TORT/MCNP coupling method for the calculation of neutron flux around a core of BWR.
Kurosawa, Masahiko
2005-01-01
For the analysis of BWR neutronics performance, accurate data are required for neutron flux distribution over the In-Reactor Pressure Vessel equipments taking into account the detailed geometrical arrangement. The TORT code can calculate neutron flux around a core of BWR in a three-dimensional geometry model, but has difficulties in fine geometrical modelling and lacks huge computer resource. On the other hand, the MCNP code enables the calculation of the neutron flux with a detailed geometry model, but requires very long sampling time to give enough number of particles. Therefore, a TORT/MCNP coupling method has been developed to eliminate the two problems mentioned above in each code. In this method, the TORT code calculates angular flux distribution on the core surface and the MCNP code calculates neutron spectrum at the points of interest using the flux distribution. The coupling method will be used as the DOT-DOMINO-MORSE code system. This TORT/MCNP coupling method was applied to calculate the neutron flux at points where induced radioactivity data were measured for 54Mn and 60Co and the radioactivity calculations based on the neutron flux obtained from the above method were compared with the measured data. PMID:16604689
Calculations of the giant-dipole-resonance photoneutrons using a coupled EGS4-morse code
Liu, J.C.; Nelson, W.R.; Kase, K.R.; Mao, X.S.
1995-10-01
The production and transport of the photoneutrons from the giant-dipoleresonance reaction have been implemented in a coupled EGS4-MORSE code. The total neutron yield (including both the direct neutron and evaporation neutron components) is calculated by folding the photoneutron yield cross sections with the photon track length distribution in the target. Empirical algorithms based on the measurements have been developed to estimate the fraction and energy of the direct neutron component for each photon. The statistical theory in the EVAP4 code, incorporated as a MORSE subroutine, is used to determine the energies of the evaporation neutrons. These represent major improvements over other calculations that assumed no direct neutrons, a constant fraction of direct neutrons, monoenergetic direct neutron, or a constant nuclear temperature for the evaporation neutrons. It was also assumed that the slow neutrons (< 2.5 MeV) are emitted isotropically and the fast neutrons are emitted anisotropically in the form of 1+Csin{sup 2}{theta}, which have a peak emission at 900. Comparisons between the calculated and the measured photoneutron results (spectra of the direct, evaporation and total neutrons; nuclear temperatures; direct neutron fractions) for materials of lead, tungsten, tantalum and copper have been made. The results show that the empirical algorithms, albeit simple, can produce reasonable results over the interested photon energy range.
Coupling procedure for TRANSURANUS and KTF codes
Jimenez, J.; Alglave, S.; Avramova, M.
2012-07-01
The nuclear industry aims to ensure safe and economic operation of each single fuel rod introduced in the reactor core. This goal is even more challenging nowadays due to the current strategy of going for higher burn-up (fuel cycles of 18 or 24 months) and longer residence time. In order to achieve that goal, fuel modeling is the key to predict the fuel rod behavior and lifetime under thermal and pressure loads, corrosion and irradiation. In this context, fuel performance codes, such as TRANSURANUS, are used to improve the fuel rod design. The modeling capabilities of the above mentioned tools can be significantly improved if they are coupled with a thermal-hydraulic code in order to have a better description of the flow conditions within the rod bundle. For LWR applications, a good representation of the two phase flow within the fuel assembly is necessary in order to have a best estimate calculation of the heat transfer inside the bundle. In this paper we present the coupling methodology of TRANSURANUS with KTF (Karlsruhe Two phase Flow subchannel code) as well as selected results of the coupling proof of principle. (authors)
Coded source neutron imaging with a MURA mask
NASA Astrophysics Data System (ADS)
Zou, Y. B.; Schillinger, B.; Wang, S.; Zhang, X. S.; Guo, Z. Y.; Lu, Y. R.
2011-09-01
In coded source neutron imaging the single aperture commonly used in neutron radiography is replaced with a coded mask. Using a coded source can improve the neutron flux at the sample plane when a very high L/ D ratio is needed. The coded source imaging is a possible way to reduce the exposure time to get a neutron image with very high L/ D ratio. A 17×17 modified uniformly redundant array coded source was tested in this work. There are 144 holes of 0.8 mm diameter on the coded source. The neutron flux from the coded source is as high as from a single 9.6 mm aperture, while its effective L/ D is the same as in the case of a 0.8 mm aperture. The Richardson-Lucy maximum likelihood algorithm was used for image reconstruction. Compared to an in-line phase contrast neutron image taken with a 1 mm aperture, it takes much less time for the coded source to get an image of similar quality.
User Guide for the R5EXEC Coupling Interface in the RELAP5-3D Code
Forsmann, J. Hope; Weaver, Walter L.
2015-04-01
This report describes the R5EXEC coupling interface in the RELAP5-3D computer code from the users perspective. The information in the report is intended for users who want to couple RELAP5-3D to other thermal-hydraulic, neutron kinetics, or control system simulation codes.
Monte Carlo code for neutron scattering instrumentation design and analysis
Daemen, L.; Fitzsimmons, M.; Hjelm, R.; Olah, G.; Roberts, J.; Seeger, P.; Smith, G.; Thelliez, T.
1996-09-01
This is the final report of a one-year, Laboratory-Directed Research and Development (LDRD) at the Los Alamos National Laboratory (LANL). The development of next generation, accelerator based neutron sources calls for the design of new instruments for neutron scattering studies of materials. It will be necessary, in the near future, to evaluate accurately and rapidly the performance of new and traditional neutron instruments at short- and long-pulse spallation neutron sources, as well as continuous sources. We have developed a code that is a design tool to assist the instrument designer model new or existing instruments, test their performance, and optimize their most important features.
Wake coupling to full potential rotor analysis code
NASA Technical Reports Server (NTRS)
Torres, Francisco J.; Chang, I-Chung; Oh, Byung K.
1990-01-01
The wake information from a helicopter forward flight code is coupled with two transonic potential rotor codes. The induced velocities for the near-, mid-, and far-wake geometries are extracted from a nonlinear rigid wake of a standard performance and analysis code. These, together with the corresponding inflow angles, computation points, and azimuth angles, are then incorporated into the transonic potential codes. The coupled codes can then provide an improved prediction of rotor blade loading at transonic speeds.
An approach for coupled-code multiphysics core simulations from a common input
Schmidt, Rodney; Belcourt, Kenneth; Hooper, Russell; Pawlowski, Roger P.; Clarno, Kevin T.; Simunovic, Srdjan; Slattery, Stuart R.; Turner, John A.; Palmtag, Scott
2014-12-10
This study describes an approach for coupled-code multiphysics reactor core simulations that is being developed by the Virtual Environment for Reactor Applications (VERA) project in the Consortium for Advanced Simulation of Light-Water Reactors (CASL). In this approach a user creates a single problem description, called the “VERAIn” common input file, to define and setup the desired coupled-code reactor core simulation. A preprocessing step accepts the VERAIn file and generates a set of fully consistent input files for the different physics codes being coupled. The problem is then solved using a single-executable coupled-code simulation tool applicable to the problem, which ismore » built using VERA infrastructure software tools and the set of physics codes required for the problem of interest. The approach is demonstrated by performing an eigenvalue and power distribution calculation of a typical three-dimensional 17 × 17 assembly with thermal–hydraulic and fuel temperature feedback. All neutronics aspects of the problem (cross-section calculation, neutron transport, power release) are solved using the Insilico code suite and are fully coupled to a thermal–hydraulic analysis calculated by the Cobra-TF (CTF) code. The single-executable coupled-code (Insilico-CTF) simulation tool is created using several VERA tools, including LIME (Lightweight Integrating Multiphysics Environment for coupling codes), DTK (Data Transfer Kit), Trilinos, and TriBITS. Parallel calculations are performed on the Titan supercomputer at Oak Ridge National Laboratory using 1156 cores, and a synopsis of the solution results and code performance is presented. Finally, ongoing development of this approach is also briefly described.« less
An approach for coupled-code multiphysics core simulations from a common input
Schmidt, Rodney; Belcourt, Kenneth; Hooper, Russell; Pawlowski, Roger P.; Clarno, Kevin T.; Simunovic, Srdjan; Slattery, Stuart R.; Turner, John A.; Palmtag, Scott
2014-12-10
This study describes an approach for coupled-code multiphysics reactor core simulations that is being developed by the Virtual Environment for Reactor Applications (VERA) project in the Consortium for Advanced Simulation of Light-Water Reactors (CASL). In this approach a user creates a single problem description, called the “VERAIn” common input file, to define and setup the desired coupled-code reactor core simulation. A preprocessing step accepts the VERAIn file and generates a set of fully consistent input files for the different physics codes being coupled. The problem is then solved using a single-executable coupled-code simulation tool applicable to the problem, which is built using VERA infrastructure software tools and the set of physics codes required for the problem of interest. The approach is demonstrated by performing an eigenvalue and power distribution calculation of a typical three-dimensional 17 × 17 assembly with thermal–hydraulic and fuel temperature feedback. All neutronics aspects of the problem (cross-section calculation, neutron transport, power release) are solved using the Insilico code suite and are fully coupled to a thermal–hydraulic analysis calculated by the Cobra-TF (CTF) code. The single-executable coupled-code (Insilico-CTF) simulation tool is created using several VERA tools, including LIME (Lightweight Integrating Multiphysics Environment for coupling codes), DTK (Data Transfer Kit), Trilinos, and TriBITS. Parallel calculations are performed on the Titan supercomputer at Oak Ridge National Laboratory using 1156 cores, and a synopsis of the solution results and code performance is presented. Finally, ongoing development of this approach is also briefly described.
Modular Code and Data System for Fast Reactor Neutronics Analyses
RIMPAULT, G.
2008-06-30
Version 00. The European Reactor ANalysis Optimized calculation System, ERANOS, has been developed and validated with the aim of providing a suitable basis for reliable neutronic calculations of current as well as advanced fast reactor cores. It consists of data libraries, deterministic codes and calculation procedures which have been developed within the European Collaboration on Fast Reactors over the past 20 years or so, in order to answer the needs of both industrial and R&D organisations. The whole system counts roughly 250 functions and 3000 subroutines totalling 450000 lines of FORTRAN-77 and ESOPE instructions. ERANOS is written using the ALOS software which requires only standard FORTRAN compilers and includes advanced programming features. A modular structure was adopted for easier evolution and incorporation of new functionalities. Blocks of data (SETs) can be created or used by the modules themselves or by the user via the LU control language. Programming, and dynamic memory allocation, are performed by means of the ESOPE language. External temporary storage and permanent storage capabilities are provided by the GEMAT and ARCHIVE functions, respectively. ESOPE, LU, GEMAT and ARCHIVE are all part of the ALOS software. This modular structure allows different modules to be linked together in procedures corresponding to recommended calculation routes ranging from fast-running and moderately-accurate 'routine' procedures to slow-running but highly-accurate 'reference' procedures. The main contents of the ERANOS-2.0 package are: nuclear data libraries (multigroup cross-sections from the JEF-2.2 evaluated nuclear data file, and other specific data files), a cell and lattice code (ECCO), reactor flux solvers (diffusion, Sn transport, nodal variational transport), a burn-up module, various processing modules (material and neutron balance, breeding gains,...), tools related to perturbation theory and sensitivity analysis, core follow-up modules (connected in the
Modular Code and Data System for Fast Reactor Neutronics Analyses
Energy Science and Technology Software Center (ESTSC)
2008-06-30
Version 00. The European Reactor ANalysis Optimized calculation System, ERANOS, has been developed and validated with the aim of providing a suitable basis for reliable neutronic calculations of current as well as advanced fast reactor cores. It consists of data libraries, deterministic codes and calculation procedures which have been developed within the European Collaboration on Fast Reactors over the past 20 years or so, in order to answer the needs of both industrial and R&Dmore » organisations. The whole system counts roughly 250 functions and 3000 subroutines totalling 450000 lines of FORTRAN-77 and ESOPE instructions. ERANOS is written using the ALOS software which requires only standard FORTRAN compilers and includes advanced programming features. A modular structure was adopted for easier evolution and incorporation of new functionalities. Blocks of data (SETs) can be created or used by the modules themselves or by the user via the LU control language. Programming, and dynamic memory allocation, are performed by means of the ESOPE language. External temporary storage and permanent storage capabilities are provided by the GEMAT and ARCHIVE functions, respectively. ESOPE, LU, GEMAT and ARCHIVE are all part of the ALOS software. This modular structure allows different modules to be linked together in procedures corresponding to recommended calculation routes ranging from fast-running and moderately-accurate 'routine' procedures to slow-running but highly-accurate 'reference' procedures. The main contents of the ERANOS-2.0 package are: nuclear data libraries (multigroup cross-sections from the JEF-2.2 evaluated nuclear data file, and other specific data files), a cell and lattice code (ECCO), reactor flux solvers (diffusion, Sn transport, nodal variational transport), a burn-up module, various processing modules (material and neutron balance, breeding gains,...), tools related to perturbation theory and sensitivity analysis, core follow-up modules (connected
Greene, N.M.; Ford, W.E. III; Petrie, L.M.; Arwood, J.W.
1992-10-01
AMPX-77 is a modular system of computer programs that pertain to nuclear analyses, with a primary emphasis on tasks associated with the production and use of multigroup cross sections. AH basic cross-section data are to be input in the formats used by the Evaluated Nuclear Data Files (ENDF/B), and output can be obtained in a variety of formats, including its own internal and very general formats, along with a variety of other useful formats used by major transport, diffusion theory, and Monte Carlo codes. Processing is provided for both neutron and gamma-my data. The present release contains codes all written in the FORTRAN-77 dialect of FORTRAN and wig process ENDF/B-V and earlier evaluations, though major modules are being upgraded in order to process ENDF/B-VI and will be released when a complete collection of usable routines is available.
Current status of the PSG Monte Carlo neutron transport code
Leppaenen, J.
2006-07-01
PSG is a new Monte Carlo neutron transport code, developed at the Technical Research Centre of Finland (VTT). The code is mainly intended for fuel assembly-level reactor physics calculations, such as group constant generation for deterministic reactor simulator codes. This paper presents the current status of the project and the essential capabilities of the code. Although the main application of PSG is in lattice calculations, the geometry is not restricted in two dimensions. This paper presents the validation of PSG against the experimental results of the three-dimensional MOX fuelled VENUS-2 reactor dosimetry benchmark. (authors)
A COUPLED TH/NEUTRONICS/CRUD FRAMEWORK IN PREDICTION OF CIPS PHENOMENON
Ling Zou; Hongbin Zhang; Jess Gehin; Brendan Kochunas
2012-04-01
A coupled TH/Neutronics/CRUD framework, which is able to simulate the CRUD deposits impact on CIPS phenomenon, was described in this paper. This framework includes the coupling among three essential physics, thermal-hydraulics, CRUD and neutronics. The overall framework was implemented by using the CFD software STAR-CCM+, developing CRUD codes, and using the neutronics code DeCART. The coupling was implemented by exchanging data between softwares using intermediate exchange files. A typical 3 by 3 PWR fuel pin problem was solved under this framework. The problem was solved in a 12 months length period of time. Time-dependent solutions were provided, including CRUD deposits inventory and their distributions on fuels, boron hideout amount inside CRUD deposits, as well as power shape changing over time. The results clearly showed the power shape suppression in regions where CRUD deposits exist, which is a strong indication of CIPS phenomenon.
CALTRANS: A parallel, deterministic, 3D neutronics code
Carson, L.; Ferguson, J.; Rogers, J.
1994-04-01
Our efforts to parallelize the deterministic solution of the neutron transport equation has culminated in a new neutronics code CALTRANS, which has full 3D capability. In this article, we describe the layout and algorithms of CALTRANS and present performance measurements of the code on a variety of platforms. Explicit implementation of the parallel algorithms of CALTRANS using both the function calls of the Parallel Virtual Machine software package (PVM 3.2) and the Meiko CS-2 tagged message passing library (based on the Intel NX/2 interface) are provided in appendices.
Coded source neutron imaging at the PULSTAR reactor
Xiao, Ziyu; Mishra, Kaushal; Hawari, Ayman; Bingham, Philip R; Bilheux, Hassina Z; Tobin Jr, Kenneth William
2011-01-01
A neutron imaging facility is located on beam-tube No.5 of the 1-MW PULSTAR reactor at North Carolina State University. An investigation of high resolution imaging using the coded source imaging technique has been initiated at the facility. Coded imaging uses a mosaic of pinholes to encode an aperture, thus generating an encoded image of the object at the detector. To reconstruct the image data received by the detector, the corresponding decoding patterns are used. The optimized design of coded mask is critical for the performance of this technique and will depend on the characteristics of the imaging beam. In this work, a 34 x 38 uniformly redundant array (URA) coded aperture system is studied for application at the PULSTAR reactor neutron imaging facility. The URA pattern was fabricated on a 500 ?m gadolinium sheet. Simulations and experiments with a pinhole object have been conducted using the Gd URA and the optimized beam line.
Magnified Neutron Radiography with Coded Sources
NASA Astrophysics Data System (ADS)
Bingham, P.; Santos-Villalobos, H.; Lavrik, N.; Gregor, J.; Bilheux, H.
A coded source imaging (CSI) system has been developed and tested at the High Flux Isotope Reactor (HFIR) CG-1D beamline at Oak Ridge National Laboratory (ORNL). The goal of this system is to use magnification to improve resolution of the imaging system beyond the detector resolution. For this system, coded masks have been manufactured at 10 μm resolution with 9 μm thick Gd patterned on Si wafers, a system model base iterative reconstruction code developed, and experiments have been performed at resolutions of 200 μm, 100 μm, 50 μm, 20 μm, and 10 μm with the object place greater than 5.5m from the detector giving magnifications up to 25 times.
Boron injection/dilution capabilities in TRACB/NEM coupled code
Jambrina, A.; Barrachina, T.; Miro, R.; Verdu, G.
2012-07-01
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)
A New Code for Proto-neutron Star Evolution
NASA Astrophysics Data System (ADS)
Roberts, L. F.
2012-08-01
A new code for following the evolution and emissions of proto-neutron stars during the first minute of their lives is developed and tested. The code is one dimensional, fully implicit, and general relativistic. Multi-group, multi-flavor neutrino transport is incorporated that makes use of variable Eddington factors obtained from a formal solution of the static general relativistic Boltzmann equation with linearized scattering terms. The timescales of neutrino emission and spectral evolution obtained using the new code are broadly consistent with previous results. Unlike other recent calculations, however, the new code predicts that the neutrino-driven wind will be characterized, at least for part of its existence, by a neutron excess. This change, potentially consequential for nucleosynthesis in the wind, is due to an improved treatment of the charged current interactions of electron-flavored neutrinos and anti-neutrinos with nucleons. A comparison is also made between the results obtained using either variable Eddington factors or simple equilibrium flux-limited diffusion. The latter approximation, which has been frequently used in previous studies of proto-neutron star cooling, accurately describes the total neutrino luminosities (to within 10%) for most of the evolution, until the proto-neutron star becomes optically thin.
A NEW CODE FOR PROTO-NEUTRON STAR EVOLUTION
Roberts, L. F.
2012-08-20
A new code for following the evolution and emissions of proto-neutron stars during the first minute of their lives is developed and tested. The code is one dimensional, fully implicit, and general relativistic. Multi-group, multi-flavor neutrino transport is incorporated that makes use of variable Eddington factors obtained from a formal solution of the static general relativistic Boltzmann equation with linearized scattering terms. The timescales of neutrino emission and spectral evolution obtained using the new code are broadly consistent with previous results. Unlike other recent calculations, however, the new code predicts that the neutrino-driven wind will be characterized, at least for part of its existence, by a neutron excess. This change, potentially consequential for nucleosynthesis in the wind, is due to an improved treatment of the charged current interactions of electron-flavored neutrinos and anti-neutrinos with nucleons. A comparison is also made between the results obtained using either variable Eddington factors or simple equilibrium flux-limited diffusion. The latter approximation, which has been frequently used in previous studies of proto-neutron star cooling, accurately describes the total neutrino luminosities (to within 10%) for most of the evolution, until the proto-neutron star becomes optically thin.
Study of an External Neutron Source for an Accelerator-Driven System using the PHITS Code
Sugawara, Takanori; Iwasaki, Tomohiko; Chiba, Takashi
2005-05-24
A code system for the Accelerator Driven System (ADS) has been under development for analyzing dynamic behaviors of a subcritical core coupled with an accelerator. This code system named DSE (Dynamics calculation code system for a Subcritical system with an External neutron source) consists of an accelerator part and a reactor part. The accelerator part employs a database, which is calculated by using PHITS, for investigating the effect related to the accelerator such as the changes of beam energy, beam diameter, void generation, and target level. This analysis method using the database may introduce some errors into dynamics calculations since the neutron source data derived from the database has some errors in fitting or interpolating procedures. In this study, the effects of various events are investigated to confirm that the method based on the database is appropriate.
Benchmarking of Neutron Production of Heavy-Ion Transport Codes
Remec, Igor; Ronningen, Reginald M.; Heilbronn, Lawrence
2012-01-01
Accurate prediction of radiation fields generated by heavy ion interactions is important in medical applications, space missions, and in design and operation of rare isotope research facilities. In recent years, several well-established computer codes in widespread use for particle and radiation transport calculations have been equipped with the capability to simulate heavy ion transport and interactions. To assess and validate these capabilities, we performed simulations of a series of benchmark-quality heavy ion experiments with the computer codes FLUKA, MARS15, MCNPX, and PHITS. We focus on the comparisons of secondary neutron production. Results are encouraging; however, further improvements in models and codes and additional benchmarking are required.
Coupled Monte Carlo neutronics and thermal hydraulics for power reactors
Bernnat, W.; Buck, M.; Mattes, M.; Zwermann, W.; Pasichnyk, I.; Velkov, K.
2012-07-01
The availability of high performance computing resources enables more and more the use of detailed Monte Carlo models even for full core power reactors. The detailed structure of the core can be described by lattices, modeled by so-called repeated structures e.g. in Monte Carlo codes such as MCNP5 or MCNPX. For cores with mainly uniform material compositions, fuel and moderator temperatures, there is no problem in constructing core models. However, when the material composition and the temperatures vary strongly a huge number of different material cells must be described which complicate the input and in many cases exceed code or memory limits. The second problem arises with the preparation of corresponding temperature dependent cross sections and thermal scattering laws. Only if these problems can be solved, a realistic coupling of Monte Carlo neutronics with an appropriate thermal-hydraulics model is possible. In this paper a method for the treatment of detailed material and temperature distributions in MCNP5 is described based on user-specified internal functions which assign distinct elements of the core cells to material specifications (e.g. water density) and temperatures from a thermal-hydraulics code. The core grid itself can be described with a uniform material specification. The temperature dependency of cross sections and thermal neutron scattering laws is taken into account by interpolation, requiring only a limited number of data sets generated for different temperatures. Applications will be shown for the stationary part of the Purdue PWR benchmark using ATHLET for thermal- hydraulics and for a generic Modular High Temperature reactor using THERMIX for thermal- hydraulics. (authors)
Interface requirements to couple thermal hydraulics codes to severe accident codes: ICARE/CATHARE
Camous, F.; Jacq, F.; Chatelard, P.
1997-07-01
In order to describe with the same code the whole sequence of severe LWR accidents, up to the vessel failure, the Institute of Protection and Nuclear Safety has performed a coupling of the severe accident code ICARE2 to the thermalhydraulics code CATHARE2. The resulting code, ICARE/CATHARE, is designed to be as pertinent as possible in all the phases of the accident. This paper is mainly devoted to the description of the ICARE2-CATHARE2 coupling.
Calleja, M.; Stieglitz, R.; Sanchez, V.; Jimenez, J.; Imke, U.
2012-07-01
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)
Development of a coupling code for PWR reactor cavity radiation streaming calculation
Zheng, Z.; Wu, H.; Cao, L.; Zheng, Y.; Zhang, H.; Wang, M.
2012-07-01
PWR reactor cavity radiation streaming is important for the safe of the personnel and equipment, thus calculation has to be performed to evaluate the neutron flux distribution around the reactor. For this calculation, the deterministic codes have difficulties in fine geometrical modeling and need huge computer resource; and the Monte Carlo codes require very long sampling time to obtain results with acceptable precision. Therefore, a coupling method has been developed to eliminate the two problems mentioned above in each code. In this study, we develop a coupling code named DORT2MCNP to link the Sn code DORT and Monte Carlo code MCNP. DORT2MCNP is used to produce a combined surface source containing top, bottom and side surface simultaneously. Because SDEF card is unsuitable for the combined surface source, we modify the SOURCE subroutine of MCNP and compile MCNP for this application. Numerical results demonstrate the correctness of the coupling code DORT2MCNP and show reasonable agreement between the coupling method and the other two codes (DORT and MCNP). (authors)
3D Multigroup Sn Neutron Transport Code
Energy Science and Technology Software Center (ESTSC)
2001-02-14
ATTILA is a 3D multigroup transport code with arbitrary order ansotropic scatter. The transport equation is solved in first order form using a tri-linear discontinuous spatial differencing on an arbitrary tetrahedral mesh. The overall solution technique is source iteration with DSA acceleration of the scattering source. Anisotropic boundary and internal sources may be entered in the form of spherical harmonics moments. Alpha and k eigenvalue problems are allowed, as well as fixed source problems. Forwardmore » and adjoint solutions are available. Reflective, vacumn, and source boundary conditions are available. ATTILA can perform charged particle transport calculations using slowing down (CSD) terms. ATTILA can also be used to peform infra-red steady-state calculations for radiative transfer purposes.« less
3D Multigroup Sn Neutron Transport Code
McGee, John; Wareing, Todd; Pautz, Shawn
2001-02-14
ATTILA is a 3D multigroup transport code with arbitrary order ansotropic scatter. The transport equation is solved in first order form using a tri-linear discontinuous spatial differencing on an arbitrary tetrahedral mesh. The overall solution technique is source iteration with DSA acceleration of the scattering source. Anisotropic boundary and internal sources may be entered in the form of spherical harmonics moments. Alpha and k eigenvalue problems are allowed, as well as fixed source problems. Forward and adjoint solutions are available. Reflective, vacumn, and source boundary conditions are available. ATTILA can perform charged particle transport calculations using slowing down (CSD) terms. ATTILA can also be used to peform infra-red steady-state calculations for radiative transfer purposes.
Interim report on fuel cycle neutronics code development.
Rabiti, C; Smith, M. A.; Kaushik, D.; Yang, W. S.
2008-05-13
As part of the Global Nuclear Energy Partnership (GNEP), a fast reactor simulation program was launched in April 2007 to develop a suite of modern simulation tools specifically for the analysis and design of sodium cooled fast reactors. The general goal of the new suite of codes is to reduce the uncertainties and biases in the various areas of reactor design activities by enhanced prediction capabilities. Under this fast reactor simulation program, a high-fidelity deterministic neutron transport code named UNIC is being developed. The final objective is to produce an integrated, advanced neutronics code that allows the high fidelity description of a nuclear reactor and simplifies the multi-step design process by direct coupling with thermal-hydraulics and structural mechanics calculations. Currently there are three solvers for the neutron transport code incorporated in UNIC: PN2ND, SN2ND, and MOCFE. PN2ND is based on a second-order even-parity spherical harmonics discretization of the transport equation and its primary target area of use is the existing homogenization approaches that are prevalent in reactor physics. MOCFE is based upon the method of characteristics applied to an unstructured finite element mesh and its primary target area of use is the fine grained nature of the explicit geometrical problems which is the long term goal of this project. SN2ND is based on a second-order, even-parity discrete ordinates discretization of the transport equation and its primary target area is the modeling transition region between the PN2ND and MOCFE solvers. The major development goal in fiscal year 2008 for the MOCFE solver was to include a two-dimensional capability that is scalable to hundreds of processors. The short term goal of this solver is to solve two-dimensional representations of reactor systems such that the energy and spatial self-shielding are accounted for and reliable cross sections can be generated for the homogeneous calculations. In this report we
A new pad-based neutron detector for stereo coded aperture thermal neutron imaging
NASA Astrophysics Data System (ADS)
Dioszegi, I.; Yu, B.; Smith, G.; Schaknowski, N.; Fried, J.; Vanier, P. E.; Salwen, C.; Forman, L.
2014-09-01
A new coded aperture thermal neutron imager system has been developed at Brookhaven National Laboratory. The cameras use a new type of position-sensitive 3He-filled ionization chamber, in which an anode plane is composed of an array of pads with independent acquisition channels. The charge is collected on each of the individual 5x5 mm2 anode pads, (48x48 in total, corresponding to 24x24 cm2 sensitive area) and read out by application specific integrated circuits (ASICs). The new design has several advantages for coded-aperture imaging applications in the field, compared to the previous generation of wire-grid based neutron detectors. Among these are its rugged design, lighter weight and use of non-flammable stopping gas. The pad-based readout occurs in parallel circuits, making it capable of high count rates, and also suitable to perform data analysis and imaging on an event-by-event basis. The spatial resolution of the detector can be better than the pixel size by using a charge sharing algorithm. In this paper we will report on the development and performance of the new pad-based neutron camera, describe a charge sharing algorithm to achieve sub-pixel spatial resolution and present the first stereoscopic coded aperture images of thermalized neutron sources using the new coded aperture thermal neutron imager system.
Bragg optics computer codes for neutron scattering instrument design
Popovici, M.; Yelon, W.B.; Berliner, R.R.; Stoica, A.D.
1997-09-01
Computer codes for neutron crystal spectrometer design, optimization and experiment planning are described. Phase space distributions, linewidths and absolute intensities are calculated by matrix methods in an extension of the Cooper-Nathans resolution function formalism. For modeling the Bragg reflection on bent crystals the lamellar approximation is used. Optimization is done by satisfying conditions of focusing in scattering and in real space, and by numerically maximizing figures of merit. Examples for three-axis and two-axis spectrometers are given.
Code Development in Coupled PARCS/RELAP5 for Supercritical Water Reactor
Hu, Po; Wilson, Paul
2014-01-01
The new capability is added to the existing coupled code package PARCS/RELAP5, in order to analyze SCWR design under supercritical pressure with the separated water coolant and moderator channels. This expansion is carried out on both codes. In PARCS, modification is focused on extending the water property tables to supercritical pressure, modifying the variable mapping input file and related code module for processing thermal-hydraulic information from separated coolant/moderator channels, and modifying neutronics feedback module to deal with the separated coolant/moderator channels. In RELAP5, modification is focused on incorporating more accurate water properties near SCWR operation/transient pressure and temperature in themore » code. Confirming tests of the modifications is presented and the major analyzing results from the extended codes package are summarized.« less
A new code for the Hall-driven magnetic evolution of neutron stars
NASA Astrophysics Data System (ADS)
Viganò, D.; Pons, J. A.; Miralles, J. A.
2012-10-01
Over the past decade, the numerical modeling of the magnetic field evolution in astrophysical scenarios has become an increasingly important field. In the crystallized crust of neutron stars the evolution of the magnetic field is governed by the Hall induction equation. In this equation the relative contribution of the two terms (Hall term and Ohmic dissipation) varies depending on the local conditions of temperature and magnetic field strength. This results in the transition from the purely parabolic character of the equations to the hyperbolic regime as the magnetic Reynolds number increases, which presents severe numerical problems. Up to now, most attempts to study this problem were based on spectral methods, but they failed in representing the transition to large magnetic Reynolds numbers. We present a new code based on upwind finite differences techniques that can handle situations with arbitrary low magnetic diffusivity and it is suitable for studying the formation of sharp current sheets during the evolution. The code is thoroughly tested in different limits and used to illustrate the evolution of the crustal magnetic field in a neutron star in some representative cases. Our code, coupled to cooling codes, can be used to perform long-term simulations of the magneto-thermal evolution of neutron stars.
Coupled receiver-decoders for low rate turbo codes
NASA Technical Reports Server (NTRS)
Hamkins, J.; Divsalar, D.
2003-01-01
We present a rate 1/31 turbo code that achieves a bit error rate of 10(sup -6) at Eb/No = -0.9 db. At this coding gain, insufficient energy per symbol is present for a conventional receiver to recover the carrier phase properly. We present a method to overcome this problem by coupling the receiver and decoder functions.
Neutron Limit on the Strongly-Coupled Chameleon Field
NASA Astrophysics Data System (ADS)
Pushin, Dmitry
2016-03-01
One of the major open questions of cosmology is the physical origin of the dark energy. There are a few sets of theories which might explain this origin that could be tested experimentally. The chameleon dark energy theory postulates self-interacting scalar field that couples to matter. This coupling induces a screening mechanism chosen so that the field amplitude is nonzero in empty space but is greatly suppressed in regions of terrestrial matter density. On behalf of the INDEX collaboration, I will report the most stringent upper bound on the free neutron-chameleon coupling in the strongly-coupled limit of the chameleon theory using neutron interferometric techniques. In our experiment we measure neutron phase induced by chameleon field. We report a 95 % confidence level upper bound on the neutron-chameleon coupling ranging from β < 4 . 7 ×106 for a Ratra-Peebles index of n = 1 in the nonlinear scalar field potential to β < 2 . 4 ×107 for n = 6 , one order of magnitude more sensitive than the most recent free neutron limit for intermediate n. This work was supported by NIST; NSF Grants: PHY-1205342, PHY-1068712, PHY-1307426; DOE award DE-FG02-97ER41042; NSERC CREATE and DISCOVERY programs; CERC; IUCSS and IU FRS program.
Adaptation of a neutron diffraction detector to coded aperture imaging
Vanier, P.E.; Forman, L.
1997-02-01
A coded aperture neutron imaging system developed at Brookhaven National Laboratory (BNL) has demonstrated that it is possible to record not only a flux of thermal neutrons at some position, but also the directions from whence they came. This realization of an idea which defied the conventional wisdom has provided a device which has never before been available to the nuclear physics community. A number of potential applications have been explored, including (1) counting warheads on a bus or in a storage area, (2) investigating inhomogeneities in drums of Pu-containing waste to facilitate non-destructive assays, (3) monitoring of vaults containing accountable materials, (4) detection of buried land mines, and (5) locating solid deposits of nuclear material held up in gaseous diffusion plants.
Valtavirta, V.
2013-07-01
This paper describes the unique way of simultaneously solving the power and temperature distributions of a nuclear system with the Monte Carlo neutron transport code Serpent 2. The coupled solution is achieved through the implementation of an internal temperature solver and material property correlations in the code. The program structure is reviewed concerning the temperature solver and the internal correlations as well as the internal coupling between these two and the neutron transport part. To estimate the reactivity differences resulting from correlation choices a simple pin-cell case has been calculated. It is established, that some correlation choices may result in difference in reactivity of approximately 100 pcm. (authors)
Bousbia-Salah, Anis; Vedovi, Juswald; D'Auria, Francesco; Ivanov, Kostadin; Galassi, Giorgio
2004-10-15
Thanks to continuous progress in computer technology, it is now possible to perform best-estimate simulations of complex scenarios in nuclear power plants. This method is carried out through the coupling of three-dimensional (3-D) neutron modeling of a reactor core into system codes. It is particularly appropriate for transients that involve strong interactions between core neutronics and reactor loop thermal hydraulics. For this purpose, the Peach Bottom boiling water reactor turbine trip test was selected to challenge the capability of such coupled codes. The test is characterized by a power excursion induced by rapid core pressurization and a self-limiting course behavior. In order to perform the closest simulation, the coupled thermal-hydraulic system code RELAP5 and 3-D neutron kinetic code PARCS were used. The obtained results are compared to those available from experimental data. Overall, the coupled code calculations globally predict the most significant observed aspects of the transient, such as the pressure wave amplitude across the core and the power course, with an acceptable agreement. However, sensitivity studies revealed that more-accurate code models should be considered in order to better match the void dynamic and the cross-section variations during transient conditions.
Coupled-Channel Computation of Direct Neutron Capture on Non-Spherical Nuclei
NASA Astrophysics Data System (ADS)
Arbanas, Goran; Thompson, Ian; Escher, Jutta; Nunes, Filomena; Elster, Charlotte; Zhang, Shi-Sheng
2014-09-01
Models of direct neutron capture of neutrons have so far accounted for the effects of non-spherical nuclei either in the incoming wave functions (via non-spherical optical model potentials), or in the final bound states (via non-spherical real potential wells), but not in both. Since it is known that spherical optical potentials do not give a good reproduction of low energy neutron-scattering observables of deformed nuclei, we have performed calculations in which the initial and final states are both treated in a self-consistent, non-spherical-nucleus picture. We have done this in the coupled-channels model of nuclear reactions implemented in the FRESCO code by using the same deformation-length for the couplings to the rotational-band states in the incoming and the final state configurations. We compute direct capture using this method for even-mass calcium isotopes 40 , 42 , 44 , 46 , 48Ca to study the effect across the two closed neutron shells, for neutron-rich even-mass tin isotopes relevant to models of astrophysical nucleosynthesis, and for 56Fe that is an important structural material used in nuclear applications. Models of direct neutron capture of neutrons have so far accounted for the effects of non-spherical nuclei either in the incoming wave functions (via non-spherical optical model potentials), or in the final bound states (via non-spherical real potential wells), but not in both. Since it is known that spherical optical potentials do not give a good reproduction of low energy neutron-scattering observables of deformed nuclei, we have performed calculations in which the initial and final states are both treated in a self-consistent, non-spherical-nucleus picture. We have done this in the coupled-channels model of nuclear reactions implemented in the FRESCO code by using the same deformation-length for the couplings to the rotational-band states in the incoming and the final state configurations. We compute direct capture using this method for even
Liu Yueqiang
2008-07-15
In order to model a magnetohydrodynamic (MHD) instability that strongly couples to external conducting structures (walls and/or coils) in a fusion device, it is often necessary to combine a MHD code solving for the plasma response, with an eddy current code computing the fields and currents of conductors. We present a rigorous proof of the coupling schemes between these two types of codes. One of the coupling schemes has been introduced and implemented in the CARMA code [R. Albanese, Y. Q. Liu, A. Portone, G. Rubinacci, and F. Villone, IEEE Trans. Magn. 44, 1654 (2008); A. Portone, F. Villone, Y. Q. Liu, R. Albanese, and G. Rubinacci, Plasma Phys. Controlled Fusion 50, 085004 (2008)] that couples the MHD code MARS-F[Y. Q. Liu, A. Bondeson, C. M. Fransson, B. Lennartson, and C. Breitholtz, Phys. Plasmas 7, 3681 (2000)] and the eddy current code CARIDDI[R. Albanese and G. Rubinacci, Adv. Imaging Electron Phys. 102, 1 (1998)]. While the coupling schemes are described for a general toroidal geometry, we give the analytical proof for a cylindrical plasma.
Neutron limit on the strongly-coupled chameleon field
NASA Astrophysics Data System (ADS)
Li, K.; Arif, M.; Cory, D. G.; Haun, R.; Heacock, B.; Huber, M. G.; Nsofini, J.; Pushin, D. A.; Saggu, P.; Sarenac, D.; Shahi, C. B.; Skavysh, V.; Snow, W. M.; Young, A. R.; Index Collaboration
2016-03-01
The physical origin of the dark energy that causes the accelerated expansion rate of the Universe is one of the major open questions of cosmology. One set of theories postulates the existence of a self-interacting scalar field for dark energy coupling to matter. In the chameleon dark energy theory, this coupling induces a screening mechanism such that the field amplitude is nonzero in empty space but is greatly suppressed in regions of terrestrial matter density. However measurements performed under appropriate vacuum conditions can enable the chameleon field to appear in the apparatus, where it can be subjected to laboratory experiments. Here we report the most stringent upper bound on the free neutron-chameleon coupling in the strongly coupled limit of the chameleon theory using neutron interferometric techniques. Our experiment sought the chameleon field through the relative phase shift it would induce along one of the neutron paths inside a perfect crystal neutron interferometer. The amplitude of the chameleon field was actively modulated by varying the millibar pressures inside a dual-chamber aluminum cell. We report a 95% confidence level upper bound on the neutron-chameleon coupling β ranging from β <4.7 ×106 for a Ratra-Peebles index of n =1 in the nonlinear scalar field potential to β <2.4 ×107 for n =6 , one order of magnitude more sensitive than the most recent free neutron limit for intermediate n . Similar experiments can explore the full parameter range for chameleon dark energy in the foreseeable future.
Neutron-induced background in charge-coupled device detectors
Jaanimagi, P. A.; Boni, R.; Keck, R. L.
2001-01-01
The inertial confinement fusion (ICF) community must become more cognizant of the neutron-induced background levels in charge-coupled device (CCD) detectors that are replacing film as the recording medium in many ICF diagnostics. This background degrades the signal-to-noise ratio (SNR) of the recorded signals and for the highest-yield shots comprises a substantial fraction of the pixel's full well capacity. CCD detectors located anywhere in the OMEGA Target Bay are precluded from recording high precision signals (SNR>30) for deuterium--tritium neutron yields greater than 10{sup 13}. CCDs make excellent calibrated neutron detectors. The average CCD background level is proportional to the neutron yield, and we have measured a linear response over four decades. The spectrum of deposited energy per pixel is heavily weighted to low energies, <50 keV, with a few isolated saturated pixels. Most of the background recorded by the CCDs is due to secondary radiation produced by interactions of the primary neutrons with all the materials in the Target Bay as well as the shield walls and the floor. Since the noise source comes from all directions it is very difficult to shield. The fallback position of using film instead of CCD cameras for high-neutron-yield target shots is flawed, as we have observed substantially increased fog levels on our x-ray recording film as a function of the neutron yield.
Neutron Activation Analysis PRognosis and Optimization Code System.
Energy Science and Technology Software Center (ESTSC)
2004-08-20
Version 00 NAAPRO predicts the results and main characteristics (detection limits, determination limits, measurement limits and relative precision of the analysis) of neutron activation analysis (instrumental and radiochemical). Gamma-ray dose rates for different points of time after sample irradiation and input count rate of the spectrometry system are also predicted. The code uses standard Windows user interface and extensive graphical tools for the visualization of the spectrometer characteristics (efficiency, response and background) and simulated spectrum.more » Optimization part is not included in the current version of the code. This release is designated NAAPRO, Version 01.beta. The MCNP code was used for generating detector responses. PREPRO-2000 and FCONV programs were used at the preparation of the program nuclear databases. A special program was developed for viewing, editing and updating of the program databases (not included into the present program package). The MCNP, PREPRO-2000 and FCONV software packages are not included in the NAAPRO package.« less
TOPLHA and ALOHA: comparison between Lower Hybrid wave coupling codes
NASA Astrophysics Data System (ADS)
Meneghini, Orso; Hillairet, J.; Goniche, M.; Bilato, R.; Voyer, D.; Parker, R.
2008-11-01
TOPLHA and ALOHA are wave coupling simulation tools for LH antennas. Both codes are able to account for realistic 3D antenna geometries and use a 1D plasma model. In the framework of a collaboration between MIT and CEA laboratories, the two codes have been extensively compared. In TOPLHA the EM problem is self consistently formulated by means of a set of multiple coupled integral equations having as domain the triangles of the meshed antenna surface. TOPLHA currently uses the FELHS code for modeling the plasma response. ALOHA instead uses a mode matching approach and its own plasma model. Comparisons have been done for several plasma scenarios on different antenna designs: an array of independent waveguides, a multi-junction antenna and a passive/active multi-junction antenna. When simulating the same geometry and plasma conditions the two codes compare remarkably well both for the reflection coefficients and for the launched spectra. The different approach of the two codes to solve the same problem strengthens the confidence in the final results.
ALOHA: an Advanced LOwer Hybrid Antenna coupling code
NASA Astrophysics Data System (ADS)
Hillairet, J.; Voyer, D.; Ekedahl, A.; Goniche, M.; Kazda, M.; Meneghini, O.; Milanesio, D.; Preynas, M.
2010-12-01
The Advanced LOwer Hybrid Antenna (ALOHA) code, has been developed to improve the modelling of the coupling of lower hybrid (LH) waves from the antenna to a cold inhomogeneous plasma while keeping a fast tool. In contrast to the previous code Slow Wave ANtenna (SWAN) (that only described the interaction of the slow wave between the waveguides and the plasma in a 1D model), the equations are now solved in 2D including the contribution of both the slow and fast waves, with a low computational cost. This approach is completed either by a full-wave computation of the antenna that takes into account its detailed geometry or by a mode-matching code dedicated to multijunctions modelling, which is convenient in preliminary design phases. Moreover, ALOHA can treat more realistic scrape-off layers in front of the antenna, by using a two-layer electron density profile. The ALOHA code has been compared with experimental results from Tore Supra LH antennas of different geometries, as well as benchmarked against other LH coupling codes, with very good results. Once validated, ALOHA has been used as a support for the design of COMPASS and ITER LH antennas and has shown to be a fast and reliable tool for LH antenna design.
Coded aperture Fast Neutron Analysis: Latest design advances
NASA Astrophysics Data System (ADS)
Accorsi, Roberto; Lanza, Richard C.
2001-07-01
Past studies have showed that materials of concern like explosives or narcotics can be identified in bulk from their atomic composition. Fast Neutron Analysis (FNA) is a nuclear method capable of providing this information even when considerable penetration is needed. Unfortunately, the cross sections of the nuclear phenomena and the solid angles involved are typically small, so that it is difficult to obtain high signal-to-noise ratios in short inspection times. CAFNAaims at combining the compound specificity of FNA with the potentially high SNR of Coded Apertures, an imaging method successfully used in far-field 2D applications. The transition to a near-field, 3D and high-energy problem prevents a straightforward application of Coded Apertures and demands a thorough optimization of the system. In this paper, the considerations involved in the design of a practical CAFNA system for contraband inspection, its conclusions, and an estimate of the performance of such a system are presented as the evolution of the ideas presented in previous expositions of the CAFNA concept.
Coupled full core neutron transport/CFD simulations of pressurized water reactors
Kochunas, B.; Stimpson, S.; Collins, B.; Downar, T.; Brewster, R.; Baglietto, E.; Yan, J.
2012-07-01
Recently as part of the CASL project, a capability to perform 3D whole-core coupled neutron transport and computational fluid dynamics (CFD) calculations was demonstrated. This work uses the 2D/1D transport code DeCART and the commercial CFD code STAR-CCM+. It builds on previous CASL work demonstrating coupling for smaller spatial domains. The coupling methodology is described along with the problem simulated and results are presented for fresh hot full power conditions. An additional comparison is made to an equivalent model that uses lower order T/H feedback to assess the importance and cost of high fidelity feedback to the neutronics problem. A simulation of a quarter core Combustion Engineering (CE) PWR core was performed with the coupled codes using a Fixed Point Gauss-Seidel iteration technique. The total approximate calculation requirements are nearly 10,000 CPU hours and 1 TB of memory. The problem took 6 coupled iterations to converge. The CFD coupled model and low order T/H feedback model compared well for global solution parameters, with a difference in the critical boron concentration and average outlet temperature of 14 ppm B and 0.94 deg. C, respectively. Differences in the power distribution were more significant with maximum relative differences in the core-wide pin peaking factor (Fq) of 5.37% and average relative differences in flat flux region power of 11.54%. Future work will focus on analyzing problems more relevant to CASL using models with less approximations. (authors)
A new neutron energy spectrum unfolding code using a two steps genetic algorithm
NASA Astrophysics Data System (ADS)
Shahabinejad, H.; Hosseini, S. A.; Sohrabpour, M.
2016-03-01
A new neutron spectrum unfolding code TGASU (Two-steps Genetic Algorithm Spectrum Unfolding) has been developed to unfold the neutron spectrum from a pulse height distribution which was calculated using the MCNPX-ESUT computational Monte Carlo code. To perform the unfolding process, the response matrices were generated using the MCNPX-ESUT computational code. Both one step (common GA) and two steps GAs have been implemented to unfold the neutron spectra. According to the obtained results, the new two steps GA code results has shown closer match in all energy regions and particularly in the high energy regions. The results of the TGASU code have been compared with those of the standard spectra, LSQR method and GAMCD code. The results of the TGASU code have been demonstrated to be more accurate than that of the existing computational codes for both under-determined and over-determined problems.
Impact of the phonon coupling on the radiative neutron capture
Avdeenkov, A. V.; Goriely, S.; Kamerdzhiev, S. P.
2010-07-15
Inclusion of the coupling of quasiparticle degrees of freedom with phonon degrees is a natural extention of the standard QRPA approach. The paper presents the quantitative impact of this phonon coupling on the dipole strength and radiative neutron capture for the stable {sup 124}Sn and very exotic {sup 150}Sn isotopes, as an illustration, using the self-consistent version of the Extended Theory of Finite Fermi Systems. It was found that the phonon contribution to the pygmy-dipole resonance and radiative neutron capture cross section is increased with the (N - Z) difference growth. The results show that the self-consistent nuclear structure calculations are important for unstable nuclei, where phenomenological approaches do not work.
Measurement of the spin-rotation coupling in neutron polarimetry
NASA Astrophysics Data System (ADS)
Demirel, Bülent; Sponar, Stephan; Hasegawa, Yuji
2015-02-01
The effect of spin-rotation coupling is measured for the first time with neutrons. The coupling of spin with the angular velocity of a rotating spin turner can be observed as a phase shift in a neutron polarimeter set-up. After the neutron’s spin is rotated by 2π through a rotating magnetic field, different phase shifts are induced for ‘up’ and ‘down’ spin eigenstates. This phase difference results in the rotation of the neutron’s spin-vector, which turns out to depend solely on the frequency of the rotation of the magnetic field. The experimental results agree well with the solutions acquired by the Pauli-Schrödinger equation.
Validation of a coupled reactive transport code in porous media
NASA Astrophysics Data System (ADS)
Mugler, C.; Montarnal, P.; Dimier, A.
2003-04-01
The safety assessment of nuclear waste disposals needs to predict the migration of radionuclides and chemical species through a geological medium. It is therefore necessary to develop and assess qualified and validated tools which integrate both the transport mechanisms through the geological media and the chemical mechanisms governing the mobility of radionuclides. In this problem, both geochemical and hydrodynamic phenomena are tightly linked together. That is the reason why the French Nuclear Energy Agency (CEA) and the French Agency for the Management of Radioactive Wastes (ANDRA) are conjointly developping a coupled reactive transport code that solves simultaneously a geochemical model and a transport model. This code, which is part of the software project ALLIANCES, leans on the libraries of two geochemical codes solving the complex ensemble of reacting chemical species: CHESS and PHREEQC. Geochemical processes considered here include ion exchange, redox reactions, acid-base reactions, surface complexation and mineral dissolution and/or precipitation. Transport is simulated using the mixed-hybrid finite element scheme CAST3M or the finite volume scheme MT3D. All together solve Darcy's law and simulate several hydrological processes such as advection, diffusion and dispersion. The coupling algorithm is an iterative sequential algorithm. Several analytical test cases have been defined and used to validate the reactive transport code. Numerical results can be compared to analytical solutions.
Limits on tensor coupling from neutron β decay
NASA Astrophysics Data System (ADS)
Pattie, R. W., Jr.; Hickerson, K. P.; Young, A. R.
2013-10-01
Limits on the tensor couplings generating a Fierz interference term b in mixed Gamow-Teller Fermi decays can be derived by combining data from measurements of angular correlation parameters in neutron decay, the neutron lifetime, and GV=GFVud as extracted from measurements of the Ft values from the 0+→0+ superallowed decay data set. These limits are derived by comparing the neutron β-decay rate as predicted in the standard model with the measured decay rate while allowing for the existence of beyond the standard model (BSM) couplings. We analyze limits derived from the electron-neutrino asymmetry a, or the beta asymmetry A, finding that the most stringent limits for CT/CA under the assumption of no right-handed neutrinos is -0.0026
An example of neutronic penalizations in reactivity transient analysis using 3D coupled chain HEMERA
Dubois, F.; Normand, B.; Sargeni, A.
2012-07-01
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)
Energy Science and Technology Software Center (ESTSC)
1991-08-01
Version: 00 The original MORSE code was a multipurpose neutron and gamma-ray transport Monte Carlo code. It was designed as a tool for solving most shielding problems. Through the use of multigroup cross sections, the solution of neutron, gamma-ray, or coupled neutron-gamma-ray problems could be obtained in either the forward or adjoint mode. Time dependence for both shielding and criticality problems is provided. General three-dimensional geometry could be used with an albedo option available atmore » any material surface. Isotropic or anisotropic scattering up to a P16 expansion of the angular distribution was allowed. MORSE-CG incorporated the Mathematical Applications, Inc. (MAGI) combinatorial geometry routines. MORSE-B modifies the Monte Carlo neutron and photon transport computer code MORSE-CG by adding routines which allow various flexible options.« less
Interface requirements to couple thermal-hydraulic codes to severe accident codes: ATHLET-CD
Trambauer, K.
1997-07-01
The system code ATHLET-CD is being developed by GRS in cooperation with IKE and IPSN. Its field of application comprises the whole spectrum of leaks and large breaks, as well as operational and abnormal transients for LWRs and VVERs. At present the analyses cover the in-vessel thermal-hydraulics, the early phases of core degradation, as well as fission products and aerosol release from the core and their transport in the Reactor Coolant System. The aim of the code development is to extend the simulation of core degradation up to failure of the reactor pressure vessel and to cover all physically reasonable accident sequences for western and eastern LWRs including RMBKs. The ATHLET-CD structure is highly modular in order to include a manifold spectrum of models and to offer an optimum basis for further development. The code consists of four general modules to describe the reactor coolant system thermal-hydraulics, the core degradation, the fission product core release, and fission product and aerosol transport. Each general module consists of some basic modules which correspond to the process to be simulated or to its specific purpose. Besides the code structure based on the physical modelling, the code follows four strictly separated steps during the course of a calculation: (1) input of structure, geometrical data, initial and boundary condition, (2) initialization of derived quantities, (3) steady state calculation or input of restart data, and (4) transient calculation. In this paper, the transient solution method is briefly presented and the coupling methods are discussed. Three aspects have to be considered for the coupling of different modules in one code system. First is the conservation of masses and energy in the different subsystems as there are fluid, structures, and fission products and aerosols. Second is the convergence of the numerical solution and stability of the calculation. The third aspect is related to the code performance, and running time.
Cooperative solutions coupling a geometry engine and adaptive solver codes
NASA Technical Reports Server (NTRS)
Dickens, Thomas P.
1995-01-01
Follow-on work has progressed in using Aero Grid and Paneling System (AGPS), a geometry and visualization system, as a dynamic real time geometry monitor, manipulator, and interrogator for other codes. In particular, AGPS has been successfully coupled with adaptive flow solvers which iterate, refining the grid in areas of interest, and continuing on to a solution. With the coupling to the geometry engine, the new grids represent the actual geometry much more accurately since they are derived directly from the geometry and do not use refits to the first-cut grids. Additional work has been done with design runs where the geometric shape is modified to achieve a desired result. Various constraints are used to point the solution in a reasonable direction which also more closely satisfies the desired results. Concepts and techniques are presented, as well as examples of sample case studies. Issues such as distributed operation of the cooperative codes versus running all codes locally and pre-calculation for performance are discussed. Future directions are considered which will build on these techniques in light of changing computer environments.
The Application Programming Interface for the PVMEXEC Program and Associated Code Coupling System
Walter L. Weaver III
2005-03-01
This report describes the Application Programming Interface for the PVMEXEC program and the code coupling systems that it implements. The information in the report is intended for programmers wanting to add a new code into the coupling system.
Valentine, T.E.; Rugama, Y. Munoz-Cobos, J.; Perez, R.
2000-10-23
The design of reactivity monitoring systems for accelerator-driven systems must be investigated to ensure that such systems remain subcritical during operation. The Monte Carlo codes LAHET and MCNP-DSP were combined together to facilitate the design of reactivity monitoring systems. The coupling of LAHET and MCNP-DSP provides a tool that can be used to simulate a variety of subcritical measurements such as the pulsed neutron, Rossi-{alpha}, or noise analysis measurements.
A novel approach to correct the coded aperture misalignment for fast neutron imaging
NASA Astrophysics Data System (ADS)
Zhang, F. N.; Hu, H. S.; Zhang, T. K.; Jia, Q. G.; Wang, D. M.; Jia, J.
2015-12-01
Aperture alignment is crucial for the diagnosis of neutron imaging because it has significant impact on the coding imaging and the understanding of the neutron source. In our previous studies on the neutron imaging system with coded aperture for large field of view, "residual watermark," certain extra information that overlies reconstructed image and has nothing to do with the source is discovered if the peak normalization is employed in genetic algorithms (GA) to reconstruct the source image. Some studies on basic properties of residual watermark indicate that the residual watermark can characterize coded aperture and can thus be used to determine the location of coded aperture relative to the system axis. In this paper, we have further analyzed the essential conditions for the existence of residual watermark and the requirements of the reconstruction algorithm for the emergence of residual watermark. A gamma coded imaging experiment has been performed to verify the existence of residual watermark. Based on the residual watermark, a correction method for the aperture misalignment has been studied. A multiple linear regression model of the position of coded aperture axis, the position of residual watermark center, and the gray barycenter of neutron source with twenty training samples has been set up. Using the regression model and verification samples, we have found the position of the coded aperture axis relative to the system axis with an accuracy of approximately 20 μm. Conclusively, a novel approach has been established to correct the coded aperture misalignment for fast neutron coded imaging.
A novel approach to correct the coded aperture misalignment for fast neutron imaging
Zhang, F. N.; Hu, H. S. Wang, D. M.; Jia, J.; Zhang, T. K.; Jia, Q. G.
2015-12-15
Aperture alignment is crucial for the diagnosis of neutron imaging because it has significant impact on the coding imaging and the understanding of the neutron source. In our previous studies on the neutron imaging system with coded aperture for large field of view, “residual watermark,” certain extra information that overlies reconstructed image and has nothing to do with the source is discovered if the peak normalization is employed in genetic algorithms (GA) to reconstruct the source image. Some studies on basic properties of residual watermark indicate that the residual watermark can characterize coded aperture and can thus be used to determine the location of coded aperture relative to the system axis. In this paper, we have further analyzed the essential conditions for the existence of residual watermark and the requirements of the reconstruction algorithm for the emergence of residual watermark. A gamma coded imaging experiment has been performed to verify the existence of residual watermark. Based on the residual watermark, a correction method for the aperture misalignment has been studied. A multiple linear regression model of the position of coded aperture axis, the position of residual watermark center, and the gray barycenter of neutron source with twenty training samples has been set up. Using the regression model and verification samples, we have found the position of the coded aperture axis relative to the system axis with an accuracy of approximately 20 μm. Conclusively, a novel approach has been established to correct the coded aperture misalignment for fast neutron coded imaging.
A novel approach to correct the coded aperture misalignment for fast neutron imaging.
Zhang, F N; Hu, H S; Zhang, T K; Jia, Q G; Wang, D M; Jia, J
2015-12-01
Aperture alignment is crucial for the diagnosis of neutron imaging because it has significant impact on the coding imaging and the understanding of the neutron source. In our previous studies on the neutron imaging system with coded aperture for large field of view, "residual watermark," certain extra information that overlies reconstructed image and has nothing to do with the source is discovered if the peak normalization is employed in genetic algorithms (GA) to reconstruct the source image. Some studies on basic properties of residual watermark indicate that the residual watermark can characterize coded aperture and can thus be used to determine the location of coded aperture relative to the system axis. In this paper, we have further analyzed the essential conditions for the existence of residual watermark and the requirements of the reconstruction algorithm for the emergence of residual watermark. A gamma coded imaging experiment has been performed to verify the existence of residual watermark. Based on the residual watermark, a correction method for the aperture misalignment has been studied. A multiple linear regression model of the position of coded aperture axis, the position of residual watermark center, and the gray barycenter of neutron source with twenty training samples has been set up. Using the regression model and verification samples, we have found the position of the coded aperture axis relative to the system axis with an accuracy of approximately 20 μm. Conclusively, a novel approach has been established to correct the coded aperture misalignment for fast neutron coded imaging. PMID:26724035
FRAPTRAN Fuel Rod Code and its Coupled Transient Analysis with the GENFLO Thermal-Hydraulic Code
Valtonen, Keijo; Hamalainen, Anitta; Cunningham, Mitchel E. )
2002-05-01
The FRAPTRAN computer code has been developed for the U.S. Nuclear Regulatory Commission (NRC) to calculate fuel behavior during power and/or cooling transients at burnup levels up to 65 MWd/kgM. FRAPTRAN has now been assessed and peer reviewed. STUK/VTT have coupled GENFLO to FRAPTRAN for calculations with improved coolant boundary conditions and prepared example calculations to show the effect of improving the coolant boundary conditions.
A neutron spectrum unfolding computer code based on artificial neural networks
NASA Astrophysics Data System (ADS)
Ortiz-Rodríguez, J. M.; Reyes Alfaro, A.; Reyes Haro, A.; Cervantes Viramontes, J. M.; Vega-Carrillo, H. R.
2014-02-01
The Bonner Spheres Spectrometer consists of a thermal neutron sensor placed at the center of a number of moderating polyethylene spheres of different diameters. From the measured readings, information can be derived about the spectrum of the neutron field where measurements were made. Disadvantages of the Bonner system are the weight associated with each sphere and the need to sequentially irradiate the spheres, requiring long exposure periods. Provided a well-established response matrix and adequate irradiation conditions, the most delicate part of neutron spectrometry, is the unfolding process. The derivation of the spectral information is not simple because the unknown is not given directly as a result of the measurements. The drawbacks associated with traditional unfolding procedures have motivated the need of complementary approaches. Novel methods based on Artificial Intelligence, mainly Artificial Neural Networks, have been widely investigated. In this work, a neutron spectrum unfolding code based on neural nets technology is presented. This code is called Neutron Spectrometry and Dosimetry with Artificial Neural networks unfolding code that was designed in a graphical interface. The core of the code is an embedded neural network architecture previously optimized using the robust design of artificial neural networks methodology. The main features of the code are: easy to use, friendly and intuitive to the user. This code was designed for a Bonner Sphere System based on a 6LiI(Eu) neutron detector and a response matrix expressed in 60 energy bins taken from an International Atomic Energy Agency compilation. The main feature of the code is that as entrance data, for unfolding the neutron spectrum, only seven rate counts measured with seven Bonner spheres are required; simultaneously the code calculates 15 dosimetric quantities as well as the total flux for radiation protection purposes. This code generates a full report with all information of the unfolding in
Neutron Transport Models and Methods for HZETRN and Coupling to Low Energy Light Ion Transport
NASA Technical Reports Server (NTRS)
Blattnig, S.R.; Slaba, T.C.; Heinbockel, J.H.
2008-01-01
Exposure estimates inside space vehicles, surface habitats, and high altitude aircraft exposed to space radiation are highly influenced by secondary neutron production. The deterministic transport code HZETRN has been identified as a reliable and efficient tool for such studies, but improvements to the underlying transport models and numerical methods are still necessary. In this paper, the forward-backward (FB) and directionally coupled forward-backward (DC) neutron transport models are derived, numerical methods for the FB model are reviewed, and a computationally efficient numerical solution is presented for the DC model. Both models are compared to the Monte Carlo codes HETCHEDS and FLUKA, and the DC model is shown to agree closely with the Monte Carlo results. Finally, it is found in the development of either model that the decoupling of low energy neutrons from the light ion (A<4) transport procedure adversely affects low energy light ion fluence spectra and exposure quantities. A first order correction is presented to resolve the problem, and it is shown to be both accurate and efficient.
Interface requirements for coupling a containment code to a reactor system thermal hydraulic codes
Baratta, A.J.
1997-07-01
To perform a complete analysis of a reactor transient, not only the primary system response but the containment response must also be accounted for. Such transients and accidents as a loss of coolant accident in both pressurized water and boiling water reactors and inadvertent operation of safety relief valves all challenge the containment and may influence flows because of containment feedback. More recently, the advanced reactor designs put forth by General Electric and Westinghouse in the US and by Framatome and Seimens in Europe rely on the containment to act as the ultimate heat sink. Techniques used by analysts and engineers to analyze the interaction of the containment and the primary system were usually iterative in nature. Codes such as RELAP or RETRAN were used to analyze the primary system response and CONTAIN or CONTEMPT the containment response. The analysis was performed by first running the system code and representing the containment as a fixed pressure boundary condition. The flows were usually from the primary system to the containment initially and generally under choked conditions. Once the mass flows and timing are determined from the system codes, these conditions were input into the containment code. The resulting pressures and temperatures were then calculated and the containment performance analyzed. The disadvantage of this approach becomes evident when one performs an analysis of a rapid depressurization or a long term accident sequence in which feedback from the containment can occur. For example, in a BWR main steam line break transient, the containment heats up and becomes a source of energy for the primary system. Recent advances in programming and computer technology are available to provide an alternative approach. The author and other researchers have developed linkage codes capable of transferring data between codes at each time step allowing discrete codes to be coupled together.
Rod Ejection Accident by the Coupled System Code ATHLET-QUABOX/CUBBOX
NASA Astrophysics Data System (ADS)
Perin, Yann; Velkov, Kiril; Pasichnyk, Igor; Langenbuch, Siegfried
The paper considers a Rod Ejection Accident (REA) which has been calculated by the coupled-code system ATHLET-QUABOX/CUBBOX. For the present study, a MOX/UOX mixed core loading was developed on the basis of a generic PWR. The results are particularly focused on the fuel enthalpy rise which is the main safety criterion for such transient. A parametric REA study has been performed, showing the influence of some important thermal-hydraulic and neutron-physical parameters. Simulations have been performed using realistic or artificially decreased delayed neutron fractions for two different core states (HZP and 30% of the nominal power). Effective fuel rod temperature influence (i.e. Doppler coefficient) has been studied by using different correlations (0.5/0.5 weighting factors or the typical TDoppler = 0.7 TSurface + 0.3 TCenter) or by changing the fuel gap conductance. It is shown that the maximum enthalpy (and enthalpy increase) does not always appear in the affected fuel assembly but can also appear in the neighboring ones. This result is a direct consequence of the burn up dependence of the enthalpy. The paper also considers the case of local delayed neutron parameters and briefly describes the future REA studies foreseen at GRS such as an investigation of quantitative uncertainty propagation from the nuclear data to the transient behavior.
Cahalan, J. E.; Ama, T.; Palmiotti, G.; Taiwo, T. A.; Yang, W. S.
2000-03-09
The VARIANT-K and DIF3D-K nodal spatial kinetics computer codes have been coupled to the SAS4A and SASSYS-1 liquid metal reactor accident and systems analysis codes. SAS4A and SASSYS-1 have been extended with the addition of heavy liquid metal (Pb and Pb-Bi) thermophysical properties, heat transfer correlations, and fluid dynamics correlations. The coupling methodology and heavy liquid metal modeling additions are described. The new computer code suite has been applied to analysis of neutron source and thermal-hydraulics transients in a model of an accelerator-driven minor actinide burner design proposed in an OECD/NEA/NSC benchmark specification. Modeling assumptions and input data generation procedures are described. Results of transient analyses are reported, with emphasis on comparison of P1 and P3 variational nodal transport theory results with nodal diffusion theory results, and on significance of spatial kinetics effects.
Myers, W.L.
1996-12-31
Analysis of a criticality accident scenario occuring in a simple fast metal system using the coupled neutronic-hydrodynamic method is demonstrated by examining the last Godiva-I criticality accident. The basis tools and information for creating a coupled neutronic-hydrodynamic code are presented. Simplifying assumptions and approximations for creating an idealized model for the Godiva-I system are discussed. Estimates of the total energy generation and the maximum attainable kinetic energy yield are the most important results that are obtained from the code. With modifications, the methodology presented in this paper can be extended to analyze criticality accident excursions in other kinds of nuclear systems.
Preliminary study of coded-source-based neutron imaging at the CPHS
NASA Astrophysics Data System (ADS)
Li, Yuanji; Huang, Zhifeng; Chen, Zhiqiang; Kang, Kejun; Xiao, Yongshun; Wang, Xuewu; Wei, Jie; Loong, C.-K.
2011-09-01
A cold neutron radiography/tomography instrument is under construction at the Compact Pulsed Hadron Source (CPHS) at Tsinghua University, China. The neutron flux is so low that an acceptable neutron radiographic image requires a long exposure time in the single-hole imaging mode. The coded-source-based imaging technique is helpful to increase the utilization of neutron flux to reduce the exposure time without loss in spatial resolution and provides high signal-to-noise ratio (SNR) images. Here we report a preliminary study on the feasibility of coded-source-based technique applied to the cold neutron imaging with a low-brilliance neutron source at the CPHS. A proper coded aperture is designed to be used in the beamline instead of the single-hole aperture. Two image retrieval algorithms, the Wiener filter algorithm and the Richardson-Lucy algorithm, are evaluated by using analytical and Monte Carlo simulations. The simulation results reveal that the coded source imaging technique is suitable for the CPHS to partially solve the problem of low neutron flux.
Non-Uniform Contrast and Noise Correction for Coded Source Neutron Imaging
Santos-Villalobos, Hector J; Bingham, Philip R
2012-01-01
Since the first application of neutron radiography in the 1930s, the field of neutron radiography has matured enough to develop several applications. However, advances in the technology are far from concluded. In general, the resolution of scintillator-based detection systems is limited to the $10\\mu m$ range, and the relatively low neutron count rate of neutron sources compared to other illumination sources restricts time resolved measurement. One path toward improved resolution is the use of magnification; however, to date neutron optics are inefficient, expensive, and difficult to develop. There is a clear demand for cost-effective scintillator-based neutron imaging systems that achieve resolutions of $1 \\mu m$ or less. Such imaging system would dramatically extend the application of neutron imaging. For such purposes a coded source imaging system is under development. The current challenge is to reduce artifacts in the reconstructed coded source images. Artifacts are generated by non-uniform illumination of the source, gamma rays, dark current at the imaging sensor, and system noise from the reconstruction kernel. In this paper, we describe how to pre-process the coded signal to reduce noise and non-uniform illumination, and how to reconstruct the coded signal with three reconstruction methods correlation, maximum likelihood estimation, and algebraic reconstruction technique. We illustrates our results with experimental examples.
Monte Carlo Predictions of Prompt Fission Neutrons and Photons: a Code Comparison
NASA Astrophysics Data System (ADS)
Talou, P.; Kawano, T.; Stetcu, I.; Vogt, R.; Randrup, J.
2014-04-01
This paper reports on initial comparisons between the LANL CGMF and LBNL/LLNL FREYA codes, which both aim at computing prompt fission neutrons and gammas. While the methodologies used in both codes are somewhat similar, the detailed implementations and physical assumptions are different. We are investigating how some of these differences impact predictions.
Interpretation of TRIGA reactivity transients with RELAP5/PARCS coupled-code
Bandini, G.; Meloni, P.; Polidori, M.
2006-07-01
In the frame of future experiments to carried out upon TRIGA reactors, which aim to verify the real feasibility of the ADS (Accelerator Driven System) concept, it is essential to build a numerical tool able to simulate the dynamic behaviour of the reactor in subcritical configuration. This model developed to support the design of subcritical experiments and the safety analysis of the reactor, as a first step has to be assessed against the experimental data available for the critical reactor. To this purpose the thermal-hydraulic/ neutronic numerical model based on the RELAP5/PARCS coupled-code is been tested against the experimental reactivity transients conducted on the RC1-TRIGA reactor at the ENEA Casaccia Research Center in forecast of the TRADE (TRIGA Accelerator Driven Experiment) subcritical experience. The results of the calculations already performed show a qualitative good agreement with the experimental data and allow to address the future developments and improvements of the numerical model. (authors)
Current and anticipated uses of thermalhydraulic and neutronic codes at PSI
Aksan, S.N.; Zimmermann, M.A.; Yadigaroglu, G.
1997-07-01
The thermalhydraulic and/or neutronic codes in use at PSI mainly provide the capability to perform deterministic safety analysis for Swiss NPPs and also serve as analysis tools for experimental facilities for LWR and ALWR simulations. In relation to these applications, physical model development and improvements, and assessment of the codes are also essential components of the activities. In this paper, a brief overview is provided on the thermalhydraulic and/or neutronic codes used for safety analysis of LWRs, at PSI, and also of some experiences and applications with these codes. Based on these experiences, additional assessment needs are indicated, together with some model improvement needs. The future needs that could be used to specify both the development of a new code and also improvement of available codes are summarized.
Energy Science and Technology Software Center (ESTSC)
1991-05-01
Version 00 MORSE-CGA was developed to add the capability of modelling rectangular lattices for nuclear reactor cores or for multipartitioned structures. It thus enhances the capability of the MORSE code system. The MORSE code is a multipurpose neutron and gamma-ray transport Monte Carlo code. It has been designed as a tool for solving most shielding problems. Through the use of multigroup cross sections, the solution of neutron, gamma-ray, or coupled neutron-gamma-ray problems may be obtainedmore » in either the forward or adjoint mode. Time dependence for both shielding and criticality problems is provided. General three-dimensional geometry may be used with an albedo option available at any material surface. Isotropic or anisotropic scattering up to a P16 expansion of the angular distribution is allowed.« less
NASA Astrophysics Data System (ADS)
Wang, Guan-bo; Yang, Xin; Qian, Da-zhi; Li, Run-dong; Tang, Bin
2014-09-01
The T(D,n)4He and 7Li(p,n)7Be neutron sources have been used for decades in nuclear physics research, stellar nucleosynthesis research and neutron therapy research. In this work, the neutron characterization including neutron yield, spectra, and angular distribution for D-T and p-7Li sources have been studied with our new binary collision approximation (BCA) based direct collision coupling method. Distinguished from the traditional path integration method for getting the neutron weight, the new model establishes a relationship between the scattering cross section and the impact parameter, which allows the secondary neutron generation carrying out jointly with ions BCA tracking. The experimental measurements of neutron characterizations have been employed for these two reactions, and the new algorithm is validated.
Thermal neutron response of a boron-coated GEM detector via GEANT4 Monte Carlo code.
Jamil, M; Rhee, J T; Kim, H G; Ahmad, Farzana; Jeon, Y J
2014-10-22
In this work, we report the design configuration and the performance of the hybrid Gas Electron Multiplier (GEM) detector. In order to make the detector sensitive to thermal neutrons, the forward electrode of the GEM has been coated with the enriched boron-10 material, which works as a neutron converter. A total of 5×5cm(2) configuration of GEM has been used for thermal neutron studies. The response of the detector has been estimated via using GEANT4 MC code with two different physics lists. Using the QGSP_BIC_HP physics list, the neutron detection efficiency was determined to be about 3%, while with QGSP_BERT_HP physics list the efficiency was around 2.5%, at the incident thermal neutron energies of 25meV. The higher response of the detector proves that GEM-coated with boron converter improves the efficiency for thermal neutrons detection. PMID:25464183
A theory manual for multi-physics code coupling in LIME.
Belcourt, Noel; Bartlett, Roscoe Ainsworth; Pawlowski, Roger Patrick; Schmidt, Rodney Cannon; Hooper, Russell Warren
2011-03-01
The Lightweight Integrating Multi-physics Environment (LIME) is a software package for creating multi-physics simulation codes. Its primary application space is when computer codes are currently available to solve different parts of a multi-physics problem and now need to be coupled with other such codes. In this report we define a common domain language for discussing multi-physics coupling and describe the basic theory associated with multiphysics coupling algorithms that are to be supported in LIME. We provide an assessment of coupling techniques for both steady-state and time dependent coupled systems. Example couplings are also demonstrated.
Zhang, L.; Lanza, R.C.
1999-12-01
The authors have developed a near field coded aperture imaging system for use with fast neutron techniques as a tool for the detection of contraband and hidden explosives through nuclear elemental analysis. The technique relies on the prompt gamma rays produced by fast neutron interactions with the object being examined. The position of the nuclear elements is determined by the location of the gamma emitters. For existing fast neutron techniques, in Pulsed Fast Neutron Analysis (PFNA), neutrons are used with very low efficiency; in Fast Neutron Analysis (FNS), the sensitivity for detection of the signature gamma rays is very low. For the Coded Aperture Fast Neutron Analysis (CAFNA{reg{underscore}sign}) the authors have developed, the efficiency for both using the probing fast neutrons and detecting the prompt gamma rays is high. For a probed volume of n{sup 3} volume elements (voxels) in a cube of n resolution elements on a side, they can compare the sensitivity with other neutron probing techniques. As compared to PFNA, the improvement for neutron utilization is n{sup 2}, where the total number of voxels in the object being examined is n{sup 3}. Compared to FNA, the improvement for gamma-ray imaging is proportional to the total open area of the coded aperture plane; a typical value is n{sup 2}/2, where n{sup 2} is the number of total detector resolution elements or the number of pixels in an object layer. It should be noted that the actual signal to noise ratio of a system depends also on the nature and distribution of background events and this comparison may reduce somewhat the effective sensitivity of CAFNA. They have performed analysis, Monte Carlo simulations, and preliminary experiments using low and high energy gamma-ray sources. The results show that a high sensitivity 3-D contraband imaging and detection system can be realized by using CAFNA.
Rohatgi, U.S.; Cheng, H.S.; Khan, H.J.; Mallen, A.N.; Neymotin, L.Y.
1998-03-01
This document describes the major modifications and improvements made to the modeling of the RAMONA-3B/MOD0 code since 1981, when the code description and assessment report was completed. The new version of the code is RAMONA-4B. RAMONA-4B is a systems transient code for application to different versions of Boiling Water Reactors (BWR) such as the current BWR, the Advanced Boiling Water Reactor (ABWR), and the Simplified Boiling Water Reactor (SBWR). This code uses a three-dimensional neutron kinetics model coupled with a multichannel, non-equilibrium, drift-flux, two-phase flow formulation of the thermal hydraulics of the reactor vessel. The code is designed to analyze a wide spectrum of BWR core and system transients and instability issues. Chapter 1 is an overview of the code`s capabilities and limitations; Chapter 2 discusses the neutron kinetics modeling and the implementation of reactivity edits. Chapter 3 is an overview of the heat conduction calculations. Chapter 4 presents modifications to the thermal-hydraulics model of the vessel, recirculation loop, steam separators, boron transport, and SBWR specific components. Chapter 5 describes modeling of the plant control and safety systems. Chapter 6 presents and modeling of Balance of Plant (BOP). Chapter 7 describes the mechanistic containment model in the code. The content of this report is complementary to the RAMONA-3B code description and assessment document. 53 refs., 81 figs., 13 tabs.
Analyses of the MSLB benchmark V1000CT-2 by the coupled system code ATHLET-BIPR8KN
Nikonov, S. P.; Langenbuch, S.; Lizorkin, M. P.; Velkov, K.
2006-07-01
Within the activities of OECD/NEA is being initiated the second phase of the VVER-1000 Coolant Transient Benchmark (V1000CT-2). It considers the best estimate analyses of a Main Steam Line Break (MSLB) of a VVER-1000 NPP with two exercises. The analyses have been performed with the coupled system code ATHLET-BIPR8KN which enables to perform realistic simulation of three-dimensional neutron kinetics and thermal-hydraulic processes in VVER NPP. Results are presented and analysed for the two proposed scenarios. These results are supplemented by sensitivity studies varying the number of the thermo-hydraulic channels (THC) in the core and by comparisons with point kinetics calculations. This work is of considerable importance for the validation of the coupled system code ATHLET-BIPR8KN in case of asymmetric core inlet conditions. (authors)
A 2D/1D coupling neutron transport method based on the matrix MOC and NEM methods
Zhang, H.; Zheng, Y.; Wu, H.; Cao, L.
2013-07-01
A new 2D/1D coupling method based on the matrix MOC method (MMOC) and nodal expansion method (NEM) is proposed for solving the three-dimensional heterogeneous neutron transport problem. The MMOC method, used for radial two-dimensional calculation, constructs a response matrix between source and flux with only one sweep and then solves the linear system by using the restarted GMRES algorithm instead of the traditional trajectory sweeping process during within-group iteration for angular flux update. Long characteristics are generated by using the customization of commercial software AutoCAD. A one-dimensional diffusion calculation is carried out in the axial direction by employing the NEM method. The 2D and ID solutions are coupled through the transverse leakage items. The 3D CMFD method is used to ensure the global neutron balance and adjust the different convergence properties of the radial and axial solvers. A computational code is developed based on these theories. Two benchmarks are calculated to verify the coupling method and the code. It is observed that the corresponding numerical results agree well with references, which indicates that the new method is capable of solving the 3D heterogeneous neutron transport problem directly. (authors)
Image enhancement using MCNP5 code and MATLAB in neutron radiography.
Tharwat, Montaser; Mohamed, Nader; Mongy, T
2014-07-01
This work presents a method that can be used to enhance the neutron radiography (NR) image for objects with high scattering materials like hydrogen, carbon and other light materials. This method used Monte Carlo code, MCNP5, to simulate the NR process and get the flux distribution for each pixel of the image and determines the scattered neutron distribution that caused image blur, and then uses MATLAB to subtract this scattered neutron distribution from the initial image to improve its quality. This work was performed before the commissioning of digital NR system in Jan. 2013. The MATLAB enhancement method is quite a good technique in the case of static based film neutron radiography, while in neutron imaging (NI) technique, image enhancement and quantitative measurement were efficient by using ImageJ software. The enhanced image quality and quantitative measurements were presented in this work. PMID:24583508
Coupled Receiver/Decoders for Low-Rate Turbo Codes
NASA Technical Reports Server (NTRS)
Hamkins, Jon; Divsalar, Dariush
2005-01-01
been proposed for receiving weak single- channel phase-modulated radio signals bearing low-rate-turbo-coded binary data. Originally intended for use in receiving telemetry signals from distant spacecraft, the proposed receiver/ decoders may also provide enhanced reception in mobile radiotelephone systems. A radio signal of the type to which the proposal applies comprises a residual carrier signal and a phase-modulated data signal. The residual carrier signal is needed as a phase reference for demodulation as a prerequisite to decoding. Low-rate turbo codes afford high coding gains and thereby enable the extraction of data from arriving radio signals that might otherwise be too weak. In the case of a conventional receiver, if the signal-to-noise ratio (specifically, the symbol energy to one-sided noise power spectral density) of the arriving signal is below approximately 0 dB, then there may not be enough energy per symbol to enable the receiver to recover properly the carrier phase. One could solve the problem at the transmitter by diverting some power from the data signal to the residual carrier. A better solution . a coupled receiver/decoder according to the proposal . could reduce the needed amount of residual carrier power. In all that follows, it is to be understood that all processing would be digital and the incoming signals to be processed would be, more precisely, outputs of analog-to-digital converters that preprocess the residual carrier and data signals at a rate of multiple samples per symbol. The upper part of the figure depicts a conventional receiving system, in which the receiver and decoder are uncoupled, and which is also called a non-data-aided system because output data from the decoder are not used in the receiver to aid in recovering the carrier phase. The receiver tracks the carrier phase from the residual carrier signal and uses the carrier phase to wipe phase noise off the data signal. The receiver typically includes a phase-locked loop
Investigation of Coded Source Neutron Imaging at the North Carolina State University PULSTAR Reactor
Xiao, Ziyu; Mishra, Kaushal; Hawari, Ayman; Bingham, Philip R; Bilheux, Hassina Z; Tobin Jr, Kenneth William
2010-10-01
A neutron imaging facility is located on beam-tube #5 of the 1-MWth PULSTAR reactor at the North Carolina State University. An investigation has been initiated to explore the application of coded imaging techniques at the facility. Coded imaging uses a mosaic of pinholes to encode an aperture, thus generating an encoded image of the object at the detector. To reconstruct the image recorded by the detector, corresponding decoding patterns are used. The optimized design of coded masks is critical for the performance of this technique and will depend on the characteristics of the imaging beam. In this work, Monte Carlo (MCNP) simulations were utilized to explore the needed modifications to the PULSTAR thermal neutron beam to support coded imaging techniques. In addition, an assessment of coded mask design has been performed. The simulations indicated that a 12 inch single crystal sapphire filter is suited for such an application at the PULSTAR beam in terms of maximizing flux with good neutron-to-gamma ratio. Computational simulations demonstrate the feasibility of correlation reconstruction methods on neutron transmission imaging. A gadolinium aperture with thickness of 500 m was used to construct the mask using a 38 34 URA pattern. A test experiment using such a URA design has been conducted and the point spread function of the system has been measured.
Page, R.; Jones, J.R.
1997-07-01
Ensuring that safety analysis needs are met in the future is likely to lead to the development of new codes and the further development of existing codes. It is therefore advantageous to define standards for data interfaces and to develop software interfacing techniques which can readily accommodate changes when they are made. Defining interface standards is beneficial but is necessarily restricted in application if future requirements are not known in detail. Code interfacing methods are of particular relevance with the move towards automatic grid frequency response operation where the integration of plant dynamic, core follow and fault study calculation tools is considered advantageous. This paper describes the background and features of a new code TALINK (Transient Analysis code LINKage program) used to provide a flexible interface to link the RELAP5 thermal hydraulics code with the PANTHER neutron kinetics and the SIBDYM whole plant dynamic modelling codes used by Nuclear Electric. The complete package enables the codes to be executed in parallel and provides an integrated whole plant thermal-hydraulics and neutron kinetics model. In addition the paper discusses the capabilities and pedigree of the component codes used to form the integrated transient analysis package and the details of the calculation of a postulated Sizewell `B` Loss of offsite power fault transient.
Energy Science and Technology Software Center (ESTSC)
2008-02-29
Version 00 (1) Problems to be solved: MVP/GMVP II can solve eigenvalue and fixed-source problems. The multigroup code GMVP can solve forward and adjoint problems for neutron, photon and neutron-photon coupled transport. The continuous-energy code MVP can solve only the forward problems. Both codes can also perform time-dependent calculations. (2) Geometry description: MVP/GMVP employs combinatorial geometry to describe the calculation geometry. It describes spatial regions by the combination of the 3-dimensional objects (BODIes). Currently, themore » following objects (BODIes) can be used. - BODIes with linear surfaces : half space, parallelepiped, right parallelepiped, wedge, right hexagonal prism - BODIes with quadratic surface and linear surfaces : cylinder, sphere, truncated right cone, truncated elliptic cone, ellipsoid by rotation, general ellipsoid - Arbitrary quadratic surface and torus The rectangular and hexagonal lattice geometry can be used to describe the repeated geometry. Furthermore, the statistical geometry model is available to treat coated fuel particles or pebbles for high temperature reactors. (3) Particle sources: The various forms of energy-, angle-, space- and time-dependent distribution functions can be specified. See Abstract for more detail.« less
Energy Science and Technology Software Center (ESTSC)
1980-10-15
Version 00 PALLAS-2DCY-FX is a code for direct integration of the transport equation in two-dimensional (r,z) geometry. It solves the energy and angular-dependent Boltzmann transport equation with general anisotropic scattering in cylindrical geometry. Its principal applications are to neutron or gamma-ray transport problems in the forward mode. The code is particularly designed for and suited to the solution of deep penetration radiation transport problems with an external (fixed) source.
ACDOS2: an improved neutron-induced dose rate code
Lagache, J.C.
1981-06-01
To calculate the expected dose rate from fusion reactors as a function of geometry, composition, and time after shutdown a computer code, ACDOS2, was written, which utilizes up-to-date libraries of cross-sections and radioisotope decay data. ACDOS2 is in ANSI FORTRAN IV, in order to make it readily adaptable elsewhere.
Neutron Activation Analysis and Product Isotope Inventory Code System.
Energy Science and Technology Software Center (ESTSC)
1990-10-31
Version 00 NAC was designed to predict the neutron-induced gamma-ray radioactivity for a wide variety of composite materials. The NAC output includes the input data, a list of all reactions for each constituent element, and the end-of-irradiation disintegration rates for each reaction. NAC also compiles a product isotope inventory containing the isotope name, the disintegration rate, the gamma-ray source strength, and the absorbed dose rate at 1 meter from an unshielded point source. The inducedmore » activity is calculated as a function of irradiation and decay times; the effect of cyclic irradiation can also be calculated.« less
Girardi, E.; Guerin, P.; Dulla, S.; Nervo, M.; Ravetto, P.
2012-07-01
Quasi-Static (QS) methods are quite popular in the reactor physics community and they exhibit two main advantages. First, these methods overcome both the limits of the Point Kinetic (PK) approach and the issues of the computational effort related to the direct discretization of the time-dependent neutron transport equation. Second, QS methods can be implemented in such a way that they can be easily coupled to very different external spatial solvers. In this paper, the results of the coupling between the QS methods developed by Politecnico di Torino and the EDF R and D core code COCAGNE are presented. The goal of these activities is to evaluate the performances of QS methods (in term of computational cost and precision) with respect to the direct kinetic solver (e.g. {theta}-scheme) already available in COCAGNE. Additionally, they allow to perform an extensive cross-validation of different kinetic models (QS and direct methods). (authors)
Neutronics-thermalhydraulics coupling in a CANDU SCWR
NASA Astrophysics Data System (ADS)
Adouki, Pierre
In order to implement new nuclear technologies as a solution to the growing demand for energy, 10 countries agreed on a framework for international cooperation in 2002, to form the Generation IV International Forum (GIF). The goal of the GIF is to design the next generation of nuclear reactors that would be cost effective and would enhance safety. This forum has proposed several types of Generation IV reactors including the Supercritical Water-Cooled Reactor (SCWR). The SCWR comes in two main configurations: pressure vessel SCWR and pressure tube SCWR (PT-SCWR). In this study, the CANDU SCWR (a PT-SCWR) is considered. This reactor is oriented vertically and contains 336 channels with a length of 5 m. The target coolant inlet and outlet temperatures are 350 Celsius and 625 Celsius, respectively. The coolant flows downwards, and the reactor power is 2540 MWth. Various fuel designs have been considered in order not to exceed the linear element rating. However, the dependency between the core power and thermalhydraulics parameters results in the necessity to use a neutronics/thermalhydaulics coupling scheme to determine the core power and the thermalhydraulics parameters. The core power obtained has a power peaking factor of 1.4. The bundle power distribution for all channels has a peak at the third bundle from the inlet, but the value of this peak increases with the channel power. The heat-transfer coefficient and the specific-heat capacity have a peak at the same location in a channel, and this location shifts toward the inlet as the channel power increases. The exit coolant temperature increases with the channel power, while the exit coolant density and pressure decrease with the channel power. Also, higher channel powers lead to higher fuel and cladding temperatures. Moreover, as the coupling method is applied, the effective multiplication factor and the values of thermalhydaulics parameters oscillate as they converge.
Wilson, W B; Perry, R T; Charlton, W S; Parish, T A; Shores, E F
2005-01-01
SOURCES is a computer code that determines neutron production rates and spectra from (alpha,n) reactions, spontaneous fission and delayed neutron emission owing to the decay of radionuclides in homogeneous media, interface problems and three-region interface problems. The code is also capable of calculating the neutron production rates due to (alpha,n) reactions induced by a monoenergetic beam of alpha particles incident on a slab of target material. The (alpha,n) spectra are calculated using an assumed isotropic angular distribution in the centre-of-mass system with a library of 107 nuclide decay alpha-particle spectra, 24 sets of measured and/or evaluated (alpha,n) cross sections and product nuclide level branching fractions, and functional alpha particle stopping cross sections for Z < 106. Spontaneous fission sources and spectra are calculated with evaluated half-life, spontaneous fission branching and Watt spectrum parameters for 44 actinides. The delayed neutron spectra are taken from an evaluated library of 105 precursors. The code outputs the magnitude and spectra of the resultant neutron sources. It also provides an analysis of the contributions to that source by each nuclide in the problem. PMID:16381695
Neutron Star Structure in the Presence of Conformally Coupled Scalar Fields
NASA Technical Reports Server (NTRS)
Sultana, Joseph; Bose, Benjamin; Kazanas, Demosthenes
2014-01-01
Neutron star models are studied in the context of scalar-tensor theories of gravity in the presence of a conformally coupled scalar field, using two different numerical equations of state (EoS) representing different degrees of stiffness. In both cases we obtain a complete solution by matching the interior numerical solution of the coupled Einstein-scalar field hydrostatic equations, with an exact metric on the surface of the star. These are then used to find the effect of the scalar field and its coupling to geometry, on the neutron star structure, particularly the maximum neutron star mass and radius. We show that in the presence of a conformally coupled scalar field, neutron stars are less dense and have smaller masses and radii than their counterparts in the minimally coupled case, and the effect increases with the magnitude of the scalar field at the center of the star.
Neutron star structure in the presence of conformally coupled scalar fields
NASA Astrophysics Data System (ADS)
Sultana, Joseph; Bose, Benjamin; Kazanas, Demosthenes
2014-10-01
Neutron star models are studied in the context of scalar-tensor theories of gravity in the presence of a conformally coupled scalar field, using two different numerical equations of state (EoS) representing different degrees of stiffness. In both cases we obtain a complete solution by matching the interior numerical solution of the coupled Einstein-scalar field hydrostatic equations, with an exact metric on the surface of the star. These are then used to find the effect of the scalar field and its coupling to geometry, on the neutron star structure, particularly the maximum neutron star mass and radius. We show that in the presence of a conformally coupled scalar field, neutron stars are less dense and have smaller masses and radii than their counterparts in the minimally coupled case, and the effect increases with the magnitude of the scalar field at the center of the star.
Not Available
2006-09-01
This document describes the recent developments of the PEBBED code suite and its application to the PBMR-400 Coupled Code Benchmark. This report addresses an FY2006 Level 2 milestone under the NGNP Design and Evaluation Methods Work Package. The milestone states "Complete a report describing the results of the application of the integrated PEBBED code package to the PBMR-400 coupled code benchmark". The report describes the current state of the PEBBED code suite, provides an overview of the Benchmark problems to which it was applied, discusses the code developments achieved in the past year, and states some of the results attained. Results of the steady state problems generated by the PEBBED fuel management code compare favorably to the preliminary results generated by codes from other participating institutions and to similar non-Benchmark analyses. Partial transient analysis capability has been achieved through the acquisition of the NEM-THERMIX code from Penn State University. Phase I of the task has been achieved through the development of a self-consistent set of tools for generating cross sections for design and transient analysis and in the successful execution of the steady state benchmark exercises.
Santos-Villalobos, Hector J; Bingham, Philip R; Gregor, Jens
2013-01-01
The limitations in neutron flux and resolution (L/D) of current neutron imaging systems can be addressed with a Coded Source Imaging system with magnification (xCSI). More precisely, the multiple sources in an xCSI system can exceed the flux of a single pinhole system for several orders of magnitude, while maintaining a higher L/D with the small sources. Moreover, designing for an xCSI system reduces noise from neutron scattering, because the object is placed away from the detector to achieve magnification. However, xCSI systems are adversely affected by correlated noise such as non-uniform illumination of the neutron source, incorrect sampling of the coded radiograph, misalignment of the coded masks, mask transparency, and the imperfection of the system Point Spread Function (PSF). We argue that a model-based reconstruction algorithm can overcome these problems and describe the implementation of a Simultaneous Iterative Reconstruction Technique algorithm for coded sources. Design pitfalls that preclude a satisfactory reconstruction are documented.
Verification of a neutronic code for transient analysis in reactors with Hex-z geometry
Gonzalez-Pintor, S.; Verdu, G.; Ginestar, D.
2012-07-01
Due to the geometry of the fuel bundles, to simulate reactors such as VVER reactors it is necessary to develop methods that can deal with hexagonal prisms as basic elements of the spatial discretization. The main features of a code based on a high order finite element method for the spatial discretization of the neutron diffusion equation and an implicit difference method for the time discretization of this equation are presented and the performance of the code is tested solving the first exercise of the AER transient benchmark. The obtained results are compared with the reference results of the benchmark and with the results provided by PARCS code. (authors)
ABAREX -- A neutron spherical optical-statistical-model code -- A user`s manual
Smith, A.B.; Lawson, R.D.
1998-06-01
The contemporary version of the neutron spherical optical-statistical-model code ABAREX is summarized with the objective of providing detailed operational guidance for the user. The physical concepts involved are very briefly outlined. The code is described in some detail and a number of explicit examples are given. With this document one should very quickly become fluent with the use of ABAREX. While the code has operated on a number of computing systems, this version is specifically tailored for the VAX/VMS work station and/or the IBM-compatible personal computer.
Borehole parametric study for neutron induced capture gamma-ray spectrometry using the MCNP code.
Shahriari, M; Sohrabpour, M
2000-01-01
The MCNP Monte Carlo code has been used to simulate neutron transport from an Am-Be source into a granite formation surrounding a borehole. The effects of the moisture and the neutron poison on the thermal neutron flux distribution and the capture by the absorbing elements has been calculated. Thermal and nonthermal captures for certain absorbers having resonance structures in the epithermal and fast energy regions such as W and Si were performed. It is shown that for those absorbers having large resonances in the epithermal regions when they are present in dry formation or when accompanied by neutron poisons the resonance captures may be significant compared to the thermal captures. PMID:10670932
Reitsma, F.; Ivanov, K.; Downar, T.; De Haas, H.; Gougar, H. D.
2006-07-01
The Pebble Bed Modular Reactor (PBMR) is a High-Temperature Gas-cooled Reactor (HTGR) concept to be built in South Africa. As part of the verification and validation program the definition and execution of code-to-code benchmark exercises are important. The Nuclear Energy Agency (NEA) of the Organisation for Economic Cooperation and Development (OECD) has accepted, through the Nuclear Science Committee (NSC), the inclusion of the Pebble-Bed Modular Reactor (PBMR) coupled neutronics/thermal hydraulics transient benchmark problem in its program. The OECD benchmark defines steady-state and transients cases, including reactivity insertion transients. It makes use of a common set of cross sections (to eliminate uncertainties between different codes) and includes specific simplifications to the design to limit the need for participants to introduce approximations in their models. In this paper the detailed specification is explained, including the test cases to be calculated and the results required from participants. (authors)
Filippone, W.L.; Baker, R.S.
1990-12-31
The neutron transport equation is solved by a hybrid method that iteratively couples regions where deterministic (S{sub N}) and stochastic (Monte Carlo) methods are applied. Unlike previous hybrid methods, the Monte Carlo and S{sub N} regions are fully coupled in the sense that no assumption is made about geometrical separation or decoupling. The hybrid method provides a new means of solving problems involving both optically thick and optically thin regions that neither Monte Carlo nor S{sub N} is well suited for by themselves. The fully coupled Monte Carlo/S{sub N} technique consists of defining spatial and/or energy regions of a problem in which either a Monte Carlo calculation or an S{sub N} calculation is to be performed. The Monte Carlo region may comprise the entire spatial region for selected energy groups, or may consist of a rectangular area that is either completely or partially embedded in an arbitrary S{sub N} region. The Monte Carlo and S{sub N} regions are then connected through the common angular boundary fluxes, which are determined iteratively using the response matrix technique, and volumetric sources. The hybrid method has been implemented in the S{sub N} code TWODANT by adding special-purpose Monte Carlo subroutines to calculate the response matrices and volumetric sources, and linkage subrountines to carry out the interface flux iterations. The common angular boundary fluxes are included in the S{sub N} code as interior boundary sources, leaving the logic for the solution of the transport flux unchanged, while, with minor modifications, the diffusion synthetic accelerator remains effective in accelerating S{sub N} calculations. The special-purpose Monte Carlo routines used are essentially analog, with few variance reduction techniques employed. However, the routines have been successfully vectorized, with approximately a factor of five increase in speed over the non-vectorized version.
Ebert, D.
1997-07-01
This is a report on the CSNI Workshop on Transient Thermal-Hydraulic and Neutronic Codes Requirements held at Annapolis, Maryland, USA November 5-8, 1996. This experts` meeting consisted of 140 participants from 21 countries; 65 invited papers were presented. The meeting was divided into five areas: (1) current and prospective plans of thermal hydraulic codes development; (2) current and anticipated uses of thermal-hydraulic codes; (3) advances in modeling of thermal-hydraulic phenomena and associated additional experimental needs; (4) numerical methods in multi-phase flows; and (5) programming language, code architectures and user interfaces. The workshop consensus identified the following important action items to be addressed by the international community in order to maintain and improve the calculational capability: (a) preserve current code expertise and institutional memory, (b) preserve the ability to use the existing investment in plant transient analysis codes, (c) maintain essential experimental capabilities, (d) develop advanced measurement capabilities to support future code validation work, (e) integrate existing analytical capabilities so as to improve performance and reduce operating costs, (f) exploit the proven advances in code architecture, numerics, graphical user interfaces, and modularization in order to improve code performance and scrutibility, and (g) more effectively utilize user experience in modifying and improving the codes.
GPU-accelerated 3D neutron diffusion code based on finite difference method
Xu, Q.; Yu, G.; Wang, K.
2012-07-01
Finite difference method, as a traditional numerical solution to neutron diffusion equation, although considered simpler and more precise than the coarse mesh nodal methods, has a bottle neck to be widely applied caused by the huge memory and unendurable computation time it requires. In recent years, the concept of General-Purpose computation on GPUs has provided us with a powerful computational engine for scientific research. In this study, a GPU-Accelerated multi-group 3D neutron diffusion code based on finite difference method was developed. First, a clean-sheet neutron diffusion code (3DFD-CPU) was written in C++ on the CPU architecture, and later ported to GPUs under NVIDIA's CUDA platform (3DFD-GPU). The IAEA 3D PWR benchmark problem was calculated in the numerical test, where three different codes, including the original CPU-based sequential code, the HYPRE (High Performance Pre-conditioners)-based diffusion code and CITATION, were used as counterpoints to test the efficiency and accuracy of the GPU-based program. The results demonstrate both high efficiency and adequate accuracy of the GPU implementation for neutron diffusion equation. A speedup factor of about 46 times was obtained, using NVIDIA's Geforce GTX470 GPU card against a 2.50 GHz Intel Quad Q9300 CPU processor. Compared with the HYPRE-based code performing in parallel on an 8-core tower server, the speedup of about 2 still could be observed. More encouragingly, without any mathematical acceleration technology, the GPU implementation ran about 5 times faster than CITATION which was speeded up by using the SOR method and Chebyshev extrapolation technique. (authors)
A Numerical Model for Coupling of Neutron Diffusion and Thermomechanics in Fast Burst Reactors
Samet Y. Kadioglu; Dana A. Knoll; Cassiano De Oliveira
2008-11-01
We develop a numerical model for coupling of neutron diffusion adn termomechanics in order to stimulate transient behavior of a fast burst reactor. The problem involves solving a set of non-linear different equations which approximate neutron diffusion, temperature change, and material behavior. With this equation set we will model the transition from a supercritical to subcritical state and possible mechanical vibration.
Neutronics code VALE for two-dimensional triagonal (hexagonal) and three-dimensional geometries
Vondy, D.R.; Fowler, T.B.
1981-08-01
This report documents the computer code VALE designed to solve multigroup neutronics problems with the diffusion theory approximation to neutron transport for a triagonal arrangement of mesh points on planes in two- and three-dimensional geometry. This code parallels the VENTURE neutronics code in the local computation system, making exposure and fuel management capabilities available. It uses and generates interface data files adopted in the cooperative effort sponsored by Reactor Physics RRT Division of the US DOE. The programming in FORTRAN is straightforward, although data is transferred in blocks between auxiliary storage devices and main core, and direct access schemes are used. The size of problems which can be handled is essentially limited only by cost of calculation since the arrays are variably dimensioned. The memory requirement is held down while data transfer during iteration is increased only as necessary with problem size. There is provision for the more common boundary conditions including the repeating boundary, 180/sup 0/ rotational symmetry, and the rotational symmetry conditions for the 30/sup 0/, 60/sup 0/, and 120/sup 0/ triangular grids on planes. A variety of types of problems may be solved: the usual neutron flux eignevalue problem, or a direct criticality search on the buckling, on a reciprocal velocity absorber (prompt mode), or on nuclide concentrations. The adjoint problem and fixed source problem may be solved, as well as the dominating higher harmonic, or the importance problem for an arbitrary fixed source.
a New Method for Neutron Capture Therapy (nct) and Related Simulation by MCNP4C Code
NASA Astrophysics Data System (ADS)
Shirazi, Mousavi; Alireza, Seyed; Ali, Taheri
2010-01-01
Neutron capture therapy (NCT) is enumerated as one of the most important methods for treatment of some strong maladies among cancers in medical science thus is unavoidable controlling and protecting instances in use of this science. Among of treatment instances of this maladies with use of nuclear medical science is use of neutron therapy that is one of the most important and effective methods in treatment of cancers. But whereas fast neutrons have too destroyer effects and also sake of protection against additional absorbed energy (absorbed dose) by tissue during neutron therapy and also naught damaging to rest of healthy tissues, should be measured absorbed energy by tissue accurately, because destroyer effects of fast neutrons is almost quintuple more than gamma photons. In this article for neutron therapy act of male's liver has been simulated a system by the Monte Carlo method (MCNP4C code) and also with use of analytical method, thus absorbed dose by this tissue has been obtained for sources with different energies accurately and has been compared results of this two methods together.
Determination of the axial-vector weak coupling constant with ultracold neutrons.
Liu, J; Mendenhall, M P; Holley, A T; Back, H O; Bowles, T J; Broussard, L J; Carr, R; Clayton, S; Currie, S; Filippone, B W; García, A; Geltenbort, P; Hickerson, K P; Hoagland, J; Hogan, G E; Hona, B; Ito, T M; Liu, C-Y; Makela, M; Mammei, R R; Martin, J W; Melconian, D; Morris, C L; Pattie, R W; Pérez Galván, A; Pitt, M L; Plaster, B; Ramsey, J C; Rios, R; Russell, R; Saunders, A; Seestrom, S J; Sondheim, W E; Tatar, E; Vogelaar, R B; VornDick, B; Wrede, C; Yan, H; Young, A R
2010-10-29
A precise measurement of the neutron decay β asymmetry A₀ has been carried out using polarized ultracold neutrons from the pulsed spallation ultracold neutron source at the Los Alamos Neutron Science Center. Combining data obtained in 2008 and 2009, we report A₀ = -0.119 66±0.000 89{-0.001 40}{+0.001 23}, from which we determine the ratio of the axial-vector to vector weak coupling of the nucleon g{A}/g{V}=-1.275 90{-0.004 45}{+0.004 09}. PMID:21231098
Coupling external radiation transport code results to the GADRAS detector response function.
Mitchell, Dean J; Thoreson, Gregory G.; Horne, Steven M.
2014-01-01
Simulating gamma spectra is useful for analyzing special nuclear materials. Gamma spectra are influenced not only by the source and the detector, but also by the external, and potentially complex scattering environment. The scattering environment can make accurate representations of gamma spectra difficult to obtain. By coupling the Monte Carlo Nuclear Particle (MCNP) code with the Gamma Detector Response and Analysis Software (GADRAS) detector response function, gamma spectrum simulations can be computed with a high degree of fidelity even in the presence of a complex scattering environment. Traditionally, GADRAS represents the external scattering environment with empirically derived scattering parameters. By modeling the external scattering environment in MCNP and using the results as input for the GADRAS detector response function, gamma spectra can be obtained with a high degree of fidelity. This method was verified with experimental data obtained in an environment with a significant amount of scattering material. The experiment used both gamma-emitting sources and moderated and bare neutron-emitting sources. The sources were modeled using GADRAS and MCNP in the presence of the external scattering environment, producing accurate representations of the experimental data.
VENTURE/PC manual: A multidimensional multigroup neutron diffusion code system. Version 3
Shapiro, A.; Huria, H.C.; Cho, K.W.
1991-12-01
VENTURE/PC is a recompilation of part of the Oak Ridge BOLD VENTURE code system, which will operate on an IBM PC or compatible computer. Neutron diffusion theory solutions are obtained for multidimensional, multigroup problems. This manual contains information associated with operating the code system. The purpose of the various modules used in the code system, and the input for these modules are discussed. The PC code structure is also given. Version 2 included several enhancements not given in the original version of the code. In particular, flux iterations can be done in core rather than by reading and writing to disk, for problems which allow sufficient memory for such in-core iterations. This speeds up the iteration process. Version 3 does not include any of the special processors used in the previous versions. These special processors utilized formatted input for various elements of the code system. All such input data is now entered through the Input Processor, which produces standard interface files for the various modules in the code system. In addition, a Standard Interface File Handbook is included in the documentation which is distributed with the code, to assist in developing the input for the Input Processor.
VENTURE/PC manual: A multidimensional multigroup neutron diffusion code system
Shapiro, A.; Huria, H.C.; Cho, K.W. )
1991-12-01
VENTURE/PC is a recompilation of part of the Oak Ridge BOLD VENTURE code system, which will operate on an IBM PC or compatible computer. Neutron diffusion theory solutions are obtained for multidimensional, multigroup problems. This manual contains information associated with operating the code system. The purpose of the various modules used in the code system, and the input for these modules are discussed. The PC code structure is also given. Version 2 included several enhancements not given in the original version of the code. In particular, flux iterations can be done in core rather than by reading and writing to disk, for problems which allow sufficient memory for such in-core iterations. This speeds up the iteration process. Version 3 does not include any of the special processors used in the previous versions. These special processors utilized formatted input for various elements of the code system. All such input data is now entered through the Input Processor, which produces standard interface files for the various modules in the code system. In addition, a Standard Interface File Handbook is included in the documentation which is distributed with the code, to assist in developing the input for the Input Processor.
Crust-core coupling in rotating neutron stars
Glampedakis, Kostas; Andersson, Nils
2006-08-15
Motivated by their gravitational wave driven instability, we investigate the influence of the crust on r-mode oscillations in a neutron star. Using a simplistic model of an elastic neutron star crust with constant shear modulus, we carry out an analytic calculation with the main objective of deriving an expression for the slippage between the core and the crust. Our analytic estimates support previous numerical results and provide useful insights into the details of the problem.
DANDE: a linked code system for core neutronics/depletion analysis
LaBauve, R.J.; England, T.R.; George, D.C.; MacFarlane, R.E.; Wilson, W.B.
1985-06-01
This report describes DANDE - a modular neutronics, depletion code system for reactor analysis. It consists of nuclear data processing, core physics, and fuel depletion modules, and allows one to use diffusion and transport methods interchangeably in core neutronics calculations. This latter capability is especially important in the design of small modular cores. Additional unique features include the capability of updating the nuclear data file during a calculation; a detailed treatment of depletion, burnable poisons as well as fuel; and the ability to make geometric changes such as control rod repositioning and fuel relocation in the course of a calculation. The detailed treatment of reactor fuel burnup, fission-product creation and decay, as well as inventories of higher-order actinides is a necessity when predicting the behavior of reactor fuel under increased burn conditions. The operation of the code system is made clear in this report by following a sample problem.
Application of three-dimensional transport code to the analysis of the neutron streaming experiment
Chatani, K.; Slater, C.O.
1990-01-01
This paper summarized the calculational results of neutron streaming through a Clinch River Breeder Reactor (CRBR) Prototype coolant pipe chaseway. Particular emphasis is placed on results at bends in the chaseway. Calculations were performed with three three-dimensional codes: the discrete ordinates radiation transport code TORT and Monte Carlo radiation transport code MORSE, which were developed by Oak Ridge National Laboratory (ORNL), and the discrete ordinates code ENSEMBLE, which was developed in Japan. The purpose of the calculations is not only to compare the calculational results with the experimental results, but also to compare the results of TORT and MORSE with those of ENSEMBLE. In the TORT calculations, two types of difference methods, weighted-difference method was applied in ENSEMBLE calculation. Both TORT and ENSEMBLE produced nearly the same calculational results, but differed in the number of iterations required for converging each neutron group. Also, the two types of difference methods in the TORT calculations showed no appreciable variance in the number of iterations required. However, a noticeable disparity in the computer times and some variation in the calculational results did occur. The comparisons of the calculational results with the experimental results, showed for the epithermal neutron flux generally good agreement in the first and second legs and at the first bend where the two-dimensional modeling might be difficult. Results were fair to poor along the centerline of the first leg near the opening to the second leg because of discrete ordinates ray effects. Additionally, the agreement was good throughout the first and second legs for the thermal neutron region. Calculations with MORSE were made. These calculational results and comparisons are described also. 8 refs., 4 figs.
Simulations of inspiraling and merging double neutron stars using the Spectral Einstein Code
NASA Astrophysics Data System (ADS)
Haas, Roland; Ott, Christian D.; Szilagyi, Bela; Kaplan, Jeffrey D.; Lippuner, Jonas; Scheel, Mark A.; Barkett, Kevin; Muhlberger, Curran D.; Dietrich, Tim; Duez, Matthew D.; Foucart, Francois; Pfeiffer, Harald P.; Kidder, Lawrence E.; Teukolsky, Saul A.
2016-06-01
We present results on the inspiral, merger, and postmerger evolution of a neutron star-neutron star (NSNS) system. Our results are obtained using the hybrid pseudospectral-finite volume Spectral Einstein Code (SpEC). To test our numerical methods, we evolve an equal-mass system for ≈22 orbits before merger. This waveform is the longest waveform obtained from fully general-relativistic simulations for NSNSs to date. Such long (and accurate) numerical waveforms are required to further improve semianalytical models used in gravitational wave data analysis, for example, the effective one body models. We discuss in detail the improvements to SpEC's ability to simulate NSNS mergers, in particular mesh refined grids to better resolve the merger and postmerger phases. We provide a set of consistency checks and compare our results to NSNS merger simulations with the independent bam code. We find agreement between them, which increases confidence in results obtained with either code. This work paves the way for future studies using long waveforms and more complex microphysical descriptions of neutron star matter in SpEC.
Rohatgi, U.S.; Cheng, H.S.; Khan, H.J.; Mallen, A.N.; Neymotin, L.Y.
1998-03-01
This document is the User`s Manual for the Boiling Water Reactor (BWR), and Simplified Boiling Water Reactor (SBWR) systems transient code RAMONA-4B. The code uses a three-dimensional neutron-kinetics model coupled with a multichannel, nonequilibrium, drift-flux, phase-flow model of the thermal hydraulics of the reactor vessel. The code is designed to analyze a wide spectrum of BWR core and system transients. Chapter 1 gives an overview of the code`s capabilities and limitations; Chapter 2 describes the code`s structure, lists major subroutines, and discusses the computer requirements. Chapter 3 is on code, auxillary codes, and instructions for running RAMONA-4B on Sun SPARC and IBM Workstations. Chapter 4 contains component descriptions and detailed card-by-card input instructions. Chapter 5 provides samples of the tabulated output for the steady-state and transient calculations and discusses the plotting procedures for the steady-state and transient calculations. Three appendices contain important user and programmer information: lists of plot variables (Appendix A) listings of input deck for sample problem (Appendix B), and a description of the plotting program PAD (Appendix C). 24 refs., 18 figs., 11 tabs.
M. J. Russell
2006-06-01
This is an assessment of codes and standards applicable to a hydrogen production plant to be coupled to a nuclear reactor. The result of the assessment is a list of codes and standards that are expected to be applicable to the plant during its design and construction.
Ortiz-Rodriguez, J. M.; Reyes Alfaro, A.; Reyes Haro, A.; Solis Sanches, L. O.; Miranda, R. Castaneda; Cervantes Viramontes, J. M.; Vega-Carrillo, H. R.
2013-07-03
In this work the performance of two neutron spectrum unfolding codes based on iterative procedures and artificial neural networks is evaluated. The first one code based on traditional iterative procedures and called Neutron spectrometry and dosimetry from the Universidad Autonoma de Zacatecas (NSDUAZ) use the SPUNIT iterative algorithm and was designed to unfold neutron spectrum and calculate 15 dosimetric quantities and 7 IAEA survey meters. The main feature of this code is the automated selection of the initial guess spectrum trough a compendium of neutron spectrum compiled by the IAEA. The second one code known as Neutron spectrometry and dosimetry with artificial neural networks (NDSann) is a code designed using neural nets technology. The artificial intelligence approach of neural net does not solve mathematical equations. By using the knowledge stored at synaptic weights on a neural net properly trained, the code is capable to unfold neutron spectrum and to simultaneously calculate 15 dosimetric quantities, needing as entrance data, only the rate counts measured with a Bonner spheres system. Similarities of both NSDUAZ and NSDann codes are: they follow the same easy and intuitive user's philosophy and were designed in a graphical interface under the LabVIEW programming environment. Both codes unfold the neutron spectrum expressed in 60 energy bins, calculate 15 dosimetric quantities and generate a full report in HTML format. Differences of these codes are: NSDUAZ code was designed using classical iterative approaches and needs an initial guess spectrum in order to initiate the iterative procedure. In NSDUAZ, a programming routine was designed to calculate 7 IAEA instrument survey meters using the fluence-dose conversion coefficients. NSDann code use artificial neural networks for solving the ill-conditioned equation system of neutron spectrometry problem through synaptic weights of a properly trained neural network. Contrary to iterative procedures, in neural
NASA Astrophysics Data System (ADS)
Ortiz-Rodríguez, J. M.; Reyes Alfaro, A.; Reyes Haro, A.; Solís Sánches, L. O.; Miranda, R. Castañeda; Cervantes Viramontes, J. M.; Vega-Carrillo, H. R.
2013-07-01
In this work the performance of two neutron spectrum unfolding codes based on iterative procedures and artificial neural networks is evaluated. The first one code based on traditional iterative procedures and called Neutron spectrometry and dosimetry from the Universidad Autonoma de Zacatecas (NSDUAZ) use the SPUNIT iterative algorithm and was designed to unfold neutron spectrum and calculate 15 dosimetric quantities and 7 IAEA survey meters. The main feature of this code is the automated selection of the initial guess spectrum trough a compendium of neutron spectrum compiled by the IAEA. The second one code known as Neutron spectrometry and dosimetry with artificial neural networks (NDSann) is a code designed using neural nets technology. The artificial intelligence approach of neural net does not solve mathematical equations. By using the knowledge stored at synaptic weights on a neural net properly trained, the code is capable to unfold neutron spectrum and to simultaneously calculate 15 dosimetric quantities, needing as entrance data, only the rate counts measured with a Bonner spheres system. Similarities of both NSDUAZ and NSDann codes are: they follow the same easy and intuitive user's philosophy and were designed in a graphical interface under the LabVIEW programming environment. Both codes unfold the neutron spectrum expressed in 60 energy bins, calculate 15 dosimetric quantities and generate a full report in HTML format. Differences of these codes are: NSDUAZ code was designed using classical iterative approaches and needs an initial guess spectrum in order to initiate the iterative procedure. In NSDUAZ, a programming routine was designed to calculate 7 IAEA instrument survey meters using the fluence-dose conversion coefficients. NSDann code use artificial neural networks for solving the ill-conditioned equation system of neutron spectrometry problem through synaptic weights of a properly trained neural network. Contrary to iterative procedures, in neural
Energy Science and Technology Software Center (ESTSC)
1985-02-01
Version 00 TP1 is a transport theory code, developed to determine reactivity effects and kinetic parameters such as effective delayed neutron fractions and mean generation time by applying the usual perturbation formalism for one-dimensional geometry.
3D Neutron Transport PWR Full-core Calculation with RMC code
NASA Astrophysics Data System (ADS)
Qiu, Yishu; She, Ding; Fan, Xiao; Wang, Kan; Li, Zeguang; Liang, Jingang; Leroyer, Hadrien
2014-06-01
Nowadays, there are more and more interests in the use of Monte Carlo codes to calculate the detailed power density distributions in full-core reactors. With the Inspur TS1000 HPC Server of Tsinghua University, several calculations have been done based on the EDF 3D Neutron Transport PWR Full-core benchmark through large-scale parallelism. To investigate and compare the results of the deterministic method and Monte Carlo method, EDF R&D and Department of Engineering Physics of Tsinghua University are having a collaboration to make code to code verification. So in this paper, two codes are used. One is the code COCAGNE developed by the EDF R&D, a deterministic core code, and the other is the Monte Carlo code RMC developed by Department of Engineering Physics in Tsinghua University. First, the full-core model is described and a 26-group calculation was performed by these two codes using the same 26-group cross-section library provided by EDF R&D. Then the parallel and tally performance of RMC is discussed. RMC employs a novel algorithm which can cut down most of the communications. It can be seen clearly that the speedup ratio almost linearly increases with the nodes. Furthermore the cell-mapping method applied by RMC consumes little time to tally even millions of cells. The results of the codes COCAGNE and RMC are compared in three ways. The results of these two codes agree well with each other. It can be concluded that both COCAGNE and RMC are able to provide 3D-transport solutions associated with detailed power density distributions calculation in PWR full-core reactors. Finally, to investigate how many histories are needed to obtain a given standard deviation for a full 3D solution, the non-symmetrized condensed 2-group fluxes of RMC are discussed.
VizieR Online Data Catalog: Neutron-star matter unified EoS FORTRAN codes (Potekhin+, 2013)
NASA Astrophysics Data System (ADS)
Potekhin, A. Y.; Fantina, A. F.; Chamel, N.; Pearson, J. M.; Goriely, S.
2013-09-01
The FORTRAN codes implement the formulae for the unified EoSs of neutron-star matter BSk19, BSk20, and BSk21, supplemented with codes for masses and radii of static neutron stars, number fractions of particles, nuclear size and shape parameters, and effective masses of nucleons. Detailed descriptions of all subroutines are given by comment lines in the files (4 data files).
Development and preliminary verification of the 3D core neutronic code: COCO
Lu, H.; Mo, K.; Li, W.; Bai, N.; Li, J.
2012-07-01
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)
NASA Astrophysics Data System (ADS)
Nuttin, A.; Capellan, N.; David, S.; Doligez, X.; El Mhari, C.; Méplan, O.
2014-06-01
Safety analysis of innovative reactor designs requires three dimensional modeling to ensure a sufficiently realistic description, starting from steady state. Actual Monte Carlo (MC) neutron transport codes are suitable candidates to simulate large complex geometries, with eventual innovative fuel. But if local values such as power densities over small regions are needed, reliable results get more difficult to obtain within an acceptable computation time. In this scope, NEA has proposed a performance test of full PWR core calculations based on Monte Carlo neutron transport, which we have used to define an optimal detail level for convergence of steady state coupled neutronics. Coupling between MCNP for neutronics and the subchannel code COBRA for thermal-hydraulics has been performed using the C++ tool MURE, developed for about ten years at LPSC and IPNO. In parallel with this study and within the same MURE framework, a simplified code of nodal kinetics based on two-group and few-point diffusion equations has been developed and validated on a typical CANDU LOCA. Methods for the computation of necessary diffusion data have been defined and applied to NU (Nat. U) and Th fuel CANDU after assembly evolutions by MURE. Simplicity of CANDU LOCA model has made possible a comparison of these two fuel behaviours during such a transient.
Automated design of coupled RF cavities using 2-D and 3-D codes
Smith, Peter; Christiansen, D. W.; Greninger, P. T.; Spalek, G.
2001-01-01
Coupled RF cavities in the Accelerator Production of Tritium Project have been designed using a procedure in which a 2-D code (CCT) searches for a design that meets frequency and coupling requirements, while a 3-D code (HFSS) is used to obtain empirical factors used by CCT to characterize the coupling slot between cavities. Using assumed values of the empirical factors, CCT runs the Superfish code iteratively to solve for a trial cavity design that has a specified frequency and coupling. The frequency shifts and the coupling constant k of the slot are modeled in CCT using a perturbation theory, the results of which are adjusted using the empirical factors. Given a trial design, HFSS is run using periodic boundary conditions to obtain a mode spectrum. The mode spectrum is processed using the DISPER code to obtain values of the coupling and the frequencies with slots. These results are used to calculate a new set of empirical factors, which are fed back into CCT for another design iteration. Cold models have been fabricated and tested to validate the codes, and results will be presented.
Gleicher, Frederick N.; Williamson, Richard L.; Ortensi, Javier; Wang, Yaqi; Spencer, Benjamin W.; Novascone, Stephen R.; Hales, Jason D.; Martineau, Richard C.
2014-10-01
The MOOSE neutron transport application RATTLESNAKE was coupled to the fuels performance application BISON to provide a higher fidelity tool for fuel performance simulation. This project is motivated by the desire to couple a high fidelity core analysis program (based on the self-adjoint angular flux equations) to a high fidelity fuel performance program, both of which can simulate on unstructured meshes. RATTLESNAKE solves self-adjoint angular flux transport equation and provides a sub-pin level resolution of the multigroup neutron flux with resonance treatment during burnup or a fast transient. BISON solves the coupled thermomechanical equations for the fuel on a sub-millimeter scale. Both applications are able to solve their respective systems on aligned and unaligned unstructured finite element meshes. The power density and local burnup was transferred from RATTLESNAKE to BISON with the MOOSE Multiapp transfer system. Multiple depletion cases were run with one-way data transfer from RATTLESNAKE to BISON. The eigenvalues are shown to agree well with values obtained from the lattice physics code DRAGON. The one-way data transfer of power density is shown to agree with the power density obtained from an internal Lassman-style model in BISON.
Thermal neutron imaging through XRQA2 GAFCHROMIC films coupled with a cadmium radiator
NASA Astrophysics Data System (ADS)
Sacco, D.; Bedogni, R.; Bortot, D.; Palomba, M.; Pola, A.; Introini, M. V.; Lorenzoli, M.; Gentile, A.; Strigari, L.; Pressello, C.; Soriani, A.; Gómez-Ros, J. M.
2015-10-01
A simple and inexpensive method to perform passive thermal neutron imaging on large areas was developed on the basis of XRQA2 GAFCHROMIC films, commonly employed for quality assurance in radiology. To enhance their thermal neutron response, the sensitive face of film was coupled with a 1 mm thick cadmium radiator, forming a sandwich. By exchanging the order of Cd filter and sensitive film with respect to the incident neutron beam direction, two different configurations (beam-Cd-film and beam-film-Cd) were identified. These configurations were tested at thermal neutrons fluence values in the range 109-1010 cm-2, using the ex-core radial thermal neutron column of the ENEA Casaccia - TRIGA reactor. The results are presented in this work.
A new paradigm for local-global coupling in whole-core neutron transport.
Lewis, E.; Smith, M.; Palmiotti, G,; Nuclear Engineering Division; Northwestern Univ.; INL
2009-01-01
A new paradigm that increases the efficiency of whole-core neutron transport calculations without lattice homogenization is introduced. Quasi-reflected interface conditions are formulated to partially decouple periodic lattice effects from global flux gradients. The starting point is the finite subelement form of the variational nodal code VARIANT that eliminates fuel-coolant homogenization through the use of heterogeneous nodes. The interface spherical harmonics expansions that couple pin-cell-sized nodes are divided into low-order and high-order terms, and reflected interface conditions are applied to the high-order terms. Combined with an integral transport method within the node, the new approach dramatically reduces both the formation time and the dimensions of the nodal response matrices and leads to sharply reduced memory requirements and computational time. The method is applied to the two-dimensional C5G7 problem, an Organisation for Economic Co-operation and Development/Nuclear Energy Agency pressurized water reactor benchmark containing mixed oxide (MOX) and UO{sub 2} fuel assemblies, as well as to a three-dimensional MOX fuel assembly. Results indicate the new approach results in very little loss of accuracy relative to the corresponding full spherical harmonics expansions while reducing computational times by well over an order of magnitude.
Neutron beam characterization measurements at the Manuel Lujan Jr. neutron scattering center
Mocko, Michal; Muhrer, Guenter; Daemen, Luke L; Kelsey, Charles T; Duran, Michael A; Tovesson, Fredrik K
2010-01-01
We have measured the neutron beam characteristics of neutron moderators at the Manuel Lujan Jr. Neutron Scattering Center at LANSCE. The absolute thermal neutron flux, energy spectra and time emission spectra were measured for the high resolution and high intensity decoupled water, partially coupled liquid hydrogen and partially coupled water moderators. The results of our experimental study will provide an insight into aging of different target-moderator-reflector-shield components as well as new experimental data for benchmarking of neutron transport codes.
NASA Astrophysics Data System (ADS)
Nelson, Adam
Multi-group scattering moment matrices are critical to the solution of the multi-group form of the neutron transport equation, as they are responsible for describing the change in direction and energy of neutrons. These matrices, however, are difficult to correctly calculate from the measured nuclear data with both deterministic and stochastic methods. Calculating these parameters when using deterministic methods requires a set of assumptions which do not hold true in all conditions. These quantities can be calculated accurately with stochastic methods, however doing so is computationally expensive due to the poor efficiency of tallying scattering moment matrices. This work presents an improved method of obtaining multi-group scattering moment matrices from a Monte Carlo neutron transport code. This improved method of tallying the scattering moment matrices is based on recognizing that all of the outgoing particle information is known a priori and can be taken advantage of to increase the tallying efficiency (therefore reducing the uncertainty) of the stochastically integrated tallies. In this scheme, the complete outgoing probability distribution is tallied, supplying every one of the scattering moment matrices elements with its share of data. In addition to reducing the uncertainty, this method allows for the use of a track-length estimation process potentially offering even further improvement to the tallying efficiency. Unfortunately, to produce the needed distributions, the probability functions themselves must undergo an integration over the outgoing energy and scattering angle dimensions. This integration is too costly to perform during the Monte Carlo simulation itself and therefore must be performed in advance by way of a pre-processing code. The new method increases the information obtained from tally events and therefore has a significantly higher efficiency than the currently used techniques. The improved method has been implemented in a code system
NASA Astrophysics Data System (ADS)
Takeda, N.; Kudo, K.; Toyokawa, H.; Torii, T.; Hashimoto, M.; Sugita, T.; Dietze, G.; Yang, X.
1999-02-01
A Monte Carlo code Neutron RESPonse function for Gas counters (NRESPG) has been developed for the calculation of neutron response functions and efficiencies for neutron energies up to 20 MeV, which can be applied for 3He, H 2, or BF 3 gas proportional counters with or without moderator. This code can simulate the neutron behavior in a two-dimensional detector configuration and treat the thermal motion of a moderator atom which becomes important as the neutron energy becomes sufficiently low. Further, a more precise measured data was taken to simulate the position-dependent gas multiplication in the sensitive and insensitive gas region of a proportional counter. The NRESPG code has been applied for the calculation of response functions of 3He cylindrical proportional counters to determine neutron energy and neutron fluence in a monoenergetic calibration field. Thus, a remarkable discrepancy in the lower portion of the full-energy peak produced by the 3He(n,p)T reaction can be removed which results in a good agreement between simulations and experiments. The code has been also used for the simulation of the response of a McTaggart-type long counter consisting of a central cylindrical BF 3 counter surrounded by a polyethylene moderator. The results of the NRESPG simulations were compared with those obtained from MCNP calculations.
Search for strongly coupled Chameleon scalar field with neutron interferometry
NASA Astrophysics Data System (ADS)
Li, K.; Arif, M.; Cory, D.; Haun, R.; Heacock, B.; Huber, M.; Nsofini, J.; Pushin, D. A.; Saggu, P.; Sarenac, D.; Shahi, C.; Skavysh, V.; Snow, M.; Young, A.
2015-04-01
The dark energy proposed to explain the observed accelerated expansion of the universe is not understood. A chameleon scalar field proposed as a dark energy candidate can explain the accelerated expansion and evade all current gravity experimental bounds. It features an effective range of the chameleon scalar field that depends on the local mass density. Hence a perfect crystal neutron interferometer, that measures relative phase shift between two paths, is a prefect tool to search for the chameleon field. We are preparing a two-chamber helium gas cell for the neutron interferometer. We can lower the pressure in one cell so low that the chameleon field range expands into the cell and causes a measurable neutron phase shift while keeping the pressure difference constant. We expect to set a new upper limit of the Chameleon field by at least one order of magnitude. This work is supported by NSF Grant 1205977, DOE Grant DE-FG02-97ER41042, Canadian Excellence Research Chairs program, Natural Sciences and Engineering Research Council of Canada and Collaborative Research and Training Experience Program
MCNP: a general Monte Carlo code for neutron and photon transport
Forster, R.A.; Godfrey, T.N.K.
1985-01-01
MCNP is a very general Monte Carlo neutron photon transport code system with approximately 250 person years of Group X-6 code development invested. It is extremely portable, user-oriented, and a true production code as it is used about 60 Cray hours per month by about 150 Los Alamos users. It has as its data base the best cross-section evaluations available. MCNP contains state-of-the-art traditional and adaptive Monte Carlo techniques to be applied to the solution of an ever-increasing number of problems. Excellent user-oriented documentation is available for all facets of the MCNP code system. Many useful and important variants of MCNP exist for special applications. The Radiation Shielding Information Center (RSIC) in Oak Ridge, Tennessee is the contact point for worldwide MCNP code and documentation distribution. A much improved MCNP Version 3A will be available in the fall of 1985, along with new and improved documentation. Future directions in MCNP development will change the meaning of MCNP to Monte Carlo N Particle where N particle varieties will be transported.
The Transient 3-D Transport Coupled Code TORT-TD/ATTICA3D for High-Fidelity Pebble-Bed HTGR Analyses
NASA Astrophysics Data System (ADS)
Seubert, Armin; Sureda, Antonio; Lapins, Janis; Bader, Johannes; Laurien, Eckart
2012-01-01
This article describes the 3D discrete ordinates-based coupled code system TORT-TD/ATTICA3D that aims at steady state and transient analyses of pebble-bed high-temperature gas cooled reactors. In view of increasing computing power, the application of time-dependent neutron transport methods becomes feasible for best estimate evaluations of safety margins. The calculation capabilities of TORT-TD/ATTICA3D are presented along with the coupling approach, with focus on the time-dependent neutron transport features of TORT-TD. Results obtained for the OECD/NEA/NSC PBMR-400 benchmark demonstrate the transient capabilities of TORT-TD/ATTICA3D.
NASA Astrophysics Data System (ADS)
Karpov, A. V.; Rachkov, V. A.; Samarin, V. V.
2015-12-01
Background: Significant enhancement of sub-barrier fusion cross sections owing to neutron transfer with positive Q values was observed in many combinations of colliding nuclei. This degree of freedom has not yet been included into the rigorous quantum coupled-channels (QCC) approach. However, the empirical coupled-channels model with neutron rearrangement [Zagrebaev, Phys. Rev. C 67, 061601 (2003), 10.1103/PhysRevC.67.061601] has already been successfully used in several papers to reproduce and predict cross sections for sub-barrier fusion reactions of stable nuclei. Purpose: The objective of this study is to combine the QCC approach and the empirical model to account for additional channels of neutron rearrangement. Method: Coupling of relative motion to collective degrees of freedom (rotation of nuclei and/or their surface vibrations) are taken into account within the QCC approach. The probability of transfer of x neutrons with a given Q value is estimated semiclassically. Results: The proposed new model was successfully tested on a few combinations of fusing nuclei 40Ca+90,96Zr, 32S+96,90, and 60,64Ni+100Mo. The calculated fusion cross sections and barrier distribution functions agree well with experimental data. Conclusions: The model developed in this work confirms all the conclusions previously made within the empirical coupled-channels model with neutron rearrangement [see Rachkov et al., Phys. Rev. C 90, 014614 (2014), 10.1103/PhysRevC.90.014614]. Moreover, it has an advantage of a more reliable microscopic account for the coupling between relative motion and the collective degrees of freedom. The proposed model can also be used to reproduce the structure of the barrier distribution function. This is a step forward to a complete solution of the long-term problem of accounting for neutron transfer channels in the QCC model.
Bruyères-le-Châtel Neutron Evaluations of Actinides with the TALYS Code: The Fission Channel
NASA Astrophysics Data System (ADS)
Romain, P.; Morillon, B.; Duarte, H.
2016-01-01
For several years, various neutron evaluations of plutonium and uranium isotopes have been performed at Bruyères-le-Châtel (BRC), from 1 keV up to 30 MeV. Since only nuclear reaction models have been used to produce these evaluations, our approach was named the "Full Model" approach. Total, shape elastic and direct inelastic cross sections were obtained from the coupled channels model using a dispersive optical potential developed for actinides, with a large enough coupling scheme including the lowest octupolar band. All other cross sections were calculated using the Hauser-Feshbach theory (TALYS code) with a pre-equilibrium component above 8-10 MeV. In this paper, we focus our attention on the fission channel. More precisely, we will present the BRC contribution to fission modeling and the philosophy adopted in our "Full Model" approach. Performing evaluations with the "Full Model" approach implies the optimization of a large number of model parameters. With increasing neutron incident energy, many residual nuclei produced by nucleon emission also lead to fission. All available experimental data assigned to various fission mechanisms of the same nucleus were used to determine fission barrier parameters. For uranium isotopes, triple-humped fission barriers were required in order to reproduce accurately variations of the experimental fission cross sections. Our BRC fission modeling has shown that the effects of the class II or class III states located in the wells of the fission barrier sometimes provide an anti-resonant transmission rather than a resonant one. Consistent evaluations were produced for a large series of U and Pu isotopes. Resulting files were tested against integral data.
Verification of SMART Neutronics Design Methodology by the MCNAP Monte Carlo Code
Jong Sung Chung; Kyung Jin Shim; Chang Hyo Kim; Chungchan Lee; Sung Quun Zee
2000-11-12
SMART is a small advanced integral pressurized water reactor (PWR) of 330 MW(thermal) designed for both electricity generation and seawater desalinization. The CASMO-3/MASTER nuclear analysis system, a design-basis of Korean PWR plants, has been employed for the SMART core nuclear design and analysis because the fuel assembly (FA) characteristics and reactor operating conditions in temperature and pressure are similar to those of PWR plants. However, the SMART FAs are highly poisoned with more than 20 Al{sub 2}O{sub 3}-B{sub 4}C plus additional Gd{sub 2}O{sub 3}/UO{sub 2} BPRs each FA. The reactor is operated with control rods inserted. Therefore, the flux and power distribution may become more distorted than those of commercial PWR plants. In addition, SMART should produce power from room temperature to hot-power operating condition because it employs nuclear heating from room temperature. This demands reliable predictions of core criticality, shutdown margin, control rod worth, power distributions, and reactivity coefficients at both room temperature and hot operating condition, yet no such data are available to verify the CASMO-3/MASTER (hereafter MASTER) code system. In the absence of experimental verification data for the SMART neutronics design, the Monte Carlo depletion analysis program MCNAP is adopted as near-term alternatives for qualifying MASTER neutronics design calculations. The MCNAP is a personal computer-based continuous energy Monte Carlo neutronics analysis program written in C++ language. We established its qualification by presenting its prediction accuracy on measurements of Venus critical facilities and core neutronics analysis of a PWR plant in operation, and depletion characteristics of integral burnable absorber FAs of the current PWR. Here, we present a comparison of MASTER and MCNAP neutronics design calculations for SMART and establish the qualification of the MASTER system.
Energy Science and Technology Software Center (ESTSC)
1996-12-19
Version 03 The NJOY nuclear data processing system is a comprehensive computer code system for producing pointwise and multigroup cross sections and related quantities from ENDF/B evaluated nuclear data in the ENDF format, including the latest US library, ENDF/B-VI. The NJOY code works with neutrons, photons, and charged particles and produces libraries for a wide variety of particle transport and reactor analysis codes.
A portable, parallel, object-oriented Monte Carlo neutron transport code in C++
Lee, S.R.; Cummings, J.C.; Nolen, S.D. |
1997-05-01
We have developed a multi-group Monte Carlo neutron transport code using C++ and the Parallel Object-Oriented Methods and Applications (POOMA) class library. This transport code, called MC++, currently computes k and {alpha}-eigenvalues and is portable to and runs parallel on a wide variety of platforms, including MPPs, clustered SMPs, and individual workstations. It contains appropriate classes and abstractions for particle transport and, through the use of POOMA, for portable parallelism. Current capabilities of MC++ are discussed, along with physics and performance results on a variety of hardware, including all Accelerated Strategic Computing Initiative (ASCI) hardware. Current parallel performance indicates the ability to compute {alpha}-eigenvalues in seconds to minutes rather than hours to days. Future plans and the implementation of a general transport physics framework are also discussed.
Zhang Tiankui; Hu Huasi; Jia Qinggang; Zhang Fengna; Liu Zhihua; Hu Guang; Guo Wei; Chen Da; Li Zhenghong; Wu Yuelei
2012-11-15
Monte-Carlo simulation of neutron coded imaging based on encoding aperture for Z-pinch of large field-of-view with 5 mm radius has been investigated, and then the coded image has been obtained. Reconstruction method of source image based on genetic algorithms (GA) has been established. 'Residual watermark,' which emerges unavoidably in reconstructed image, while the peak normalization is employed in GA fitness calculation because of its statistical fluctuation amplification, has been discovered and studied. Residual watermark is primarily related to the shape and other parameters of the encoding aperture cross section. The properties and essential causes of the residual watermark were analyzed, while the identification on equivalent radius of aperture was provided. By using the equivalent radius, the reconstruction can also be accomplished without knowing the point spread function (PSF) of actual aperture. The reconstruction result is close to that by using PSF of the actual aperture.
Osiris: A Modern, High-Performance, Coupled, Multi-Physics Code For Nuclear Reactor Core Analysis
Procassini, R J; Chand, K K; Clouse, C J; Ferencz, R M; Grandy, J M; Henshaw, W D; Kramer, K J; Parsons, I D
2007-02-26
To meet the simulation needs of the GNEP program, LLNL is leveraging a suite of high-performance codes to be used in the development of a multi-physics tool for modeling nuclear reactor cores. The Osiris code project, which began last summer, is employing modern computational science techniques in the development of the individual physics modules and the coupling framework. Initial development is focused on coupling thermal-hydraulics and neutral-particle transport, while later phases of the project will add thermal-structural mechanics and isotope depletion. Osiris will be applicable to the design of existing and future reactor systems through the use of first-principles, coupled physics models with fine-scale spatial resolution in three dimensions and fine-scale particle-energy resolution. Our intent is to replace an existing set of legacy, serial codes which require significant approximations and assumptions, with an integrated, coupled code that permits the design of a reactor core using a first-principles physics approach on a wide range of computing platforms, including the world's most powerful parallel computers. A key research activity of this effort deals with the efficient and scalable coupling of physics modules which utilize rather disparate mesh topologies. Our approach allows each code module to use a mesh topology and resolution that is optimal for the physics being solved, and employs a mesh-mapping and data-transfer module to effect the coupling. Additional research is planned in the area of scalable, parallel thermal-hydraulics, high-spatial-accuracy depletion and coupled-physics simulation using Monte Carlo transport.
Brantley, P S
2006-09-27
We describe an asymptotic analysis of the coupled nonlinear system of equations describing time-dependent three-dimensional monoenergetic neutron transport and isotopic depletion and radioactive decay. The classic asymptotic diffusion scaling of Larsen and Keller [1], along with a consistent small scaling of the terms describing the radioactive decay of isotopes, is applied to this coupled nonlinear system of equations in a medium of specified initial isotopic composition. The analysis demonstrates that to leading order the neutron transport equation limits to the standard time-dependent neutron diffusion equation with macroscopic cross sections whose number densities are determined by the standard system of ordinary differential equations, the so-called Bateman equations, describing the temporal evolution of the nuclide number densities.
Incorporation of coupled nonequilibrium chemistry into a two-dimensional nozzle code (SEAGULL)
NASA Technical Reports Server (NTRS)
Ratliff, A. W.
1979-01-01
A two-dimensional multiple shock nozzle code (SEAGULL) was extended to include the effects of finite rate chemistry. The basic code that treats multiple shocks and contact surfaces was fully coupled with a generalized finite rate chemistry and vibrational energy exchange package. The modified code retains all of the original SEAGULL features plus the capability to treat chemical and vibrational nonequilibrium reactions. Any chemical and/or vibrational energy exchange mechanism can be handled as long as thermodynamic data and rate constants are available for all participating species.
Yang, W. S.; Lee, C. H.
2008-05-16
Under the fast reactor simulation program launched in April 2007, development of an advanced multigroup cross section generation code was initiated in July 2007, in conjunction with the development of the high-fidelity deterministic neutron transport code UNIC. The general objectives are to simplify the existing multi-step schemes and to improve the resolved and unresolved resonance treatments. Based on the review results of current methods and the fact that they have been applied successfully to fast critical experiment analyses and fast reactor designs for last three decades, the methodologies of the ETOE-2/MC{sup 2}-2/SDX code system were selected as the starting set of methodologies for multigroup cross section generation for fast reactor analysis. As the first step for coupling with the UNIC code and use in a parallel computing environment, the MC{sup 2}-2 code was updated by modernizing the memory structure and replacing old data management package subroutines and functions with FORTRAN 90 based routines. Various modifications were also made in the ETOE-2 and MC{sup 2}-2 codes to process the ENDF/B-VII.0 data properly. Using the updated ETOE-2/MC{sup 2}-2 code system, the ENDF/B-VII.0 data was successfully processed for major heavy and intermediate nuclides employed in sodium-cooled fast reactors. Initial verification tests of the MC{sup 2}-2 libraries generated from ENDF/B-VII.0 data were performed by inter-comparison of twenty-one group infinite dilute total cross sections obtained from MC{sup 2}-2, VIM, and NJOY. For almost all nuclides considered, MC{sup 2}-2 cross sections agreed very well with those from VIM and NJOY. Preliminary validation tests of the ENDF/B-VII.0 libraries of MC{sup 2}-2 were also performed using a set of sixteen fast critical benchmark problems. The deterministic results based on MC{sup 2}-2/TWODANT calculations were in good agreement with MCNP solutions within {approx}0.25% {Delta}{rho}, except a few small LANL fast assemblies
NASA Astrophysics Data System (ADS)
Tsokaros, Antonios; UryÅ«, Kōji; Rezzolla, Luciano
2015-05-01
We present the extension of our cocal—Compact Object CALculator—code to compute general-relativistic initial data for binary compact-star systems. In particular, we construct quasiequilibrium initial data for equal-mass binaries with spins that are either aligned or antialigned with the orbital angular momentum. The Isenberg-Wilson-Mathews formalism is adopted and the constraint equations are solved using the representation formula with a suitable choice of a Green's function. We validate the new code with solutions for equal-mass binaries and explore its capabilities for a wide range of compactnesses, from a white dwarf binary with compactness ˜1 0-4, up to a highly relativistic neutron-star binary with compactness ˜0.22 . We also present a comparison with corotating and irrotational quasiequilibrium sequences from the spectral code lorene [Taniguchi and Gourgoulhon, Phys. Rev. D 66, 104019 (2002)] and with different compactness, showing that the results from the two codes agree to a precision of the order of 0.05%. Finally, we present equilibria for spinning configurations with a nuclear-physics equation of state in a piecewise polytropic representation.
Santos-Villalobos, Hector J; Gregor, Jens; Bingham, Philip R
2014-01-01
At the present, neutron sources cannot be fabricated small and powerful enough in order to achieve high resolution radiography while maintaining an adequate flux. One solution is to employ computational imaging techniques such as a Magnified Coded Source Imaging (CSI) system. A coded-mask is placed between the neutron source and the object. The system resolution is increased by reducing the size of the mask holes and the flux is increased by increasing the size of the coded-mask and/or the number of holes. One limitation of such system is that the resolution of current state-of-the-art scintillator-based detectors caps around 50um. To overcome this challenge, the coded-mask and object are magnified by making the distance from the coded-mask to the object much smaller than the distance from object to detector. In previous work, we have shown via synthetic experiments that our least squares method outperforms other methods in image quality and reconstruction precision because of the modeling of the CSI system components. However, the validation experiments were limited to simplistic neutron sources. In this work, we aim to model the flux distribution of a real neutron source and incorporate such a model in our least squares computational system. We provide a full description of the methodology used to characterize the neutron source and validate the method with synthetic experiments.
RSAP - A Code for Display of Neutron Cross Section Data and SAMMY Fit Results
Sayer, R.O.
2001-02-02
RSAP is a computer code for display of neutron cross section data and selected SAMMY output. SAMMY is a multilevel R-matrix code for fitting neutron time-of-flight cross-section data using Bayes' method. RSAP, which runs on the Digital Unix Alpha platform, reads ORELA Data Files (ODF) created by SAMMY and uses graphics routines from the PLPLOT package. In addition, RSAP can read data and/or computed values from ASCII files with a format specified by the user. Plot output may be displayed in an X window, sent to a postscript file (rsap.ps), or sent to a color postscript file (rsap.psc). Thirteen plot types are supported, allowing the user to display cross section data, transmission data, errors, theory, Bayes fits, and residuals in various combinations. In this document the designations theory and Bayes refer to the initial and final theoretical cross sections, respectively, as evaluated by SAMMY. Special plot types include Bayes/Data, Theory--Data, and Bayes--Data. Output from two SAMMY runs may be compared by plotting the ratios Theory2/Theory1 and Bayes2/Bayes1 or by plotting the differences (Theory2-Theory1) and (Bayes2-Bayes1).
Validation of the MCNPX-PoliMi Code to Design a Fast-Neutron Multiplicity Counter
J. L. Dolan; A. C. Kaplan; M. Flaska; S. A. Pozzi; D. L. Chichester
2012-07-01
Many safeguards measurement systems used at nuclear facilities, both domestically and internationally, rely on He-3 detectors and well established mathematical equations to interpret coincidence and multiplicity-type measurements for verifying quantities of special nuclear material. Due to resource shortages alternatives to these existing He-3 based systems are being sought. Work is also underway to broaden the capabilities of these types of measurement systems in order to improve current multiplicity analysis techniques. As a part of a Material Protection, Accounting, and Control Technology (MPACT) project within the U.S. Department of Energy's Fuel Cycle Technology Program we are designing a fast-neutron multiplicity counter with organic liquid scintillators to quantify important quantities such as plutonium mass. We are also examining the potential benefits of using fast-neutron detectors for multiplicity analysis of advanced fuels in comparison with He-3 detectors and testing the performance of such designs. The designs are being developed and optimized using the MCNPX-PoliMi transport code to study detector response. In the full paper, we will discuss validation measurements used to justify the use of the MCNPX-PoliMi code paired with the MPPost multiplicity routine to design a fast neutron multiplicity counter with liquid scintillators. This multiplicity counter will be designed with the end goal of safeguarding advanced nuclear fuels. With improved timing qualities associated with liquid scintillation detectors, we can design a system that is less limited by nuclear materials of high activities. Initial testing of the designed system with nuclear fuels will take place at Idaho National Laboratory in a later stage of this collaboration.
Ortiz-Rodriguez, J. M.; Reyes Alfaro, A.; Reyes Haro, A.; Solis Sanches, L. O.; Miranda, R. Castaneda; Cervantes Viramontes, J. M.; Vega-Carrillo, H. R.
2013-07-03
In this work a neutron spectrum unfolding code, based on artificial intelligence technology is presented. The code called ''Neutron Spectrometry and Dosimetry with Artificial Neural Networks and two Bonner spheres'', (NSDann2BS), was designed in a graphical user interface under the LabVIEW programming environment. The main features of this code are to use an embedded artificial neural network architecture optimized with the ''Robust design of artificial neural networks methodology'' and to use two Bonner spheres as the only piece of information. In order to build the code here presented, once the net topology was optimized and properly trained, knowledge stored at synaptic weights was extracted and using a graphical framework build on the LabVIEW programming environment, the NSDann2BS code was designed. This code is friendly, intuitive and easy to use for the end user. The code is freely available upon request to authors. To demonstrate the use of the neural net embedded in the NSDann2BS code, the rate counts of {sup 252}Cf, {sup 241}AmBe and {sup 239}PuBe neutron sources measured with a Bonner spheres system.
NASA Astrophysics Data System (ADS)
Ortiz-Rodríguez, J. M.; Reyes Alfaro, A.; Reyes Haro, A.; Solís Sánches, L. O.; Miranda, R. Castañeda; Cervantes Viramontes, J. M.; Vega-Carrillo, H. R.
2013-07-01
In this work a neutron spectrum unfolding code, based on artificial intelligence technology is presented. The code called "Neutron Spectrometry and Dosimetry with Artificial Neural Networks and two Bonner spheres", (NSDann2BS), was designed in a graphical user interface under the LabVIEW programming environment. The main features of this code are to use an embedded artificial neural network architecture optimized with the "Robust design of artificial neural networks methodology" and to use two Bonner spheres as the only piece of information. In order to build the code here presented, once the net topology was optimized and properly trained, knowledge stored at synaptic weights was extracted and using a graphical framework build on the LabVIEW programming environment, the NSDann2BS code was designed. This code is friendly, intuitive and easy to use for the end user. The code is freely available upon request to authors. To demonstrate the use of the neural net embedded in the NSDann2BS code, the rate counts of 252Cf, 241AmBe and 239PuBe neutron sources measured with a Bonner spheres system.
Bandini, B.R. Los Alamos National Lab., NM
1990-05-01
No present light water reactor accident analysis code employs both high state of the art neutronics and thermal-hydraulics computational algorithms. Adding a modern three-dimensional neutron kinetics model to the present TRAC-PFI/MOD2 code would create a fully up to date pressurized water reactor accident evaluation code. After reviewing several options, it was decided that the Nodal Expansion Method would best provide the basis for this multidimensional transient neutronic analysis capability. Steady-state and transient versions of the Nodal Expansion Method were coded in both three-dimensional Cartesian and cylindrical geometries. In stand-alone form this method of solving the few group neutron diffusion equations was shown to yield efficient and accurate results for a variety of steady-state and transient benchmark problems. The Nodal Expansion Method was then incorporated into TRAC-PFl/MOD2. The combined NEM/TRAC code results agreed well with the EPRI-ARROTTA core-only transient analysis code when modelling a severe PWR control rod ejection accident.
Samarin, V. V.
2015-10-15
The fusion cross sections for the {sup 17,18}O+{sup 27}Al, {sup 18}O+{sup 58}Ni, and {sup 6}He+{sup 197}Au reactions were calculated by the coupled-channel method. The radial dependence of matrices that describe coupling to valence-neutron-rearrangement channels was determined with the aid of two-center wave functions. The coupling-strength parameters were evaluated on the basis of numerically solving the time-dependent Schrödinger equation. Satisfactory agreement with experimental data was obtained.
Daavittila, Antti; Haemaelaeinen, Anitta; Kyrki-Rajamaeki, Riitta
2003-05-15
All of the three exercises of the Organization for Economic Cooperation and Development/Nuclear Regulatory Commission pressurized water reactor main steam line break (PWR MSLB) benchmark were calculated at VTT, the Technical Research Centre of Finland. For the first exercise, the plant simulation with point-kinetic neutronics, the thermal-hydraulics code SMABRE was used. The second exercise was calculated with the three-dimensional reactor dynamics code TRAB-3D, and the third exercise with the combination TRAB-3D/SMABRE. VTT has over ten years' experience of coupling neutronic and thermal-hydraulic codes, but this benchmark was the first time these two codes, both developed at VTT, were coupled together. The coupled code system is fast and efficient; the total computation time of the 100-s transient in the third exercise was 16 min on a modern UNIX workstation. The results of all the exercises are similar to those of the other participants. In order to demonstrate the effect of secondary circuit modeling on the results, three different cases were calculated. In case 1 there is no phase separation in the steam lines and no flow reversal in the aspirator. In case 2 the flow reversal in the aspirator is allowed, but there is no phase separation in the steam lines. Finally, in case 3 the drift-flux model is used for the phase separation in the steam lines, but the aspirator flow reversal is not allowed. With these two modeling variations, it is possible to cover a remarkably broad range of results. The maximum power level reached after the reactor trip varies from 534 to 904 MW, the range of the time of the power maximum being close to 30 s. Compared to the total calculated transient time of 100 s, the effect of the secondary side modeling is extremely important.
Comparisons between the TOPLHA and the ALOHA codes on Lower Hybrid antenna coupling
NASA Astrophysics Data System (ADS)
Hillairet, J.; Meneghini, O.; Milanesio, D.; Voyer, D.; Goniche, M.
2009-11-01
Comparisons have been made between TOPLHA and ALOHA, two codes for Lower Hybrid antenna coupling to cold inhomogeneous plasmas, on the Tore-Supra C2 multi-junction and an ITER relevant Passive-Active Multi-junction (PAM). Calculation for several plasma scenarios of reflexion coefficients and launched spectra are in good agreement.
Modeling equal and unequal mass binary neutron star mergers using public codes
NASA Astrophysics Data System (ADS)
De Pietri, Roberto; Feo, Alessandra; Maione, Francesco; Löffler, Frank
2016-03-01
We present three-dimensional simulations of the dynamics of binary neutron star mergers from the late stage of the inspiral process up to ˜20 ms after the system has merged, either to form a hypermassive neutron star or a rotating black hole. We investigate five equal mass models of total gravitational mass 2.207, 2.373, 2.537, 2.697, and 2.854 M⊙, respectively, and four unequal mass models with MADM≃2.53 M⊙ and q ≃0.94 , 0.88, 0.83, and 0.77 (where q =M(1 )/M(2 ) is the mass ratio). We use a semirealistic equation of state, namely, the seven-segment piecewise polytropic SLyPP with a thermal component given by Γth=1.8 . We have also compared the resulting dynamics (for one model) using both the BSSN-NOK and CCZ4 methods for the evolution of the gravitational sector and also different reconstruction methods for the matter sector, namely, PPM, WENO, and MP5. Our results show agreement at high resolution, but superiority of BSSN-NOK supplemented by WENO reconstruction at lower resolutions. One of the important characteristics of the present investigation is that for the first time it has been done using only publicly available open source software: the Einstein Toolkit code, deployed for the dynamical evolution, and the LORENE code, for the generation of the initial models. All of the source code and parameters used for the simulations have been made publicly available. This not only makes it possible to rerun and reanalyze our data but also enables others to directly build upon this work for future research.
A Complex-Geometry Validation Experiment for Advanced Neutron Transport Codes
David W. Nigg; Anthony W. LaPorta; Joseph W. Nielsen; James Parry; Mark D. DeHart; Samuel E. Bays; William F. Skerjanc
2013-11-01
The Idaho National Laboratory (INL) has initiated a focused effort to upgrade legacy computational reactor physics software tools and protocols used for support of core fuel management and experiment management in the Advanced Test Reactor (ATR) and its companion critical facility (ATRC) at the INL.. This will be accomplished through the introduction of modern high-fidelity computational software and protocols, with appropriate new Verification and Validation (V&V) protocols, over the next 12-18 months. Stochastic and deterministic transport theory based reactor physics codes and nuclear data packages that support this effort include MCNP5[1], SCALE/KENO6[2], HELIOS[3], SCALE/NEWT[2], and ATTILA[4]. Furthermore, a capability for sensitivity analysis and uncertainty quantification based on the TSUNAMI[5] system has also been implemented. Finally, we are also evaluating the Serpent[6] and MC21[7] codes, as additional verification tools in the near term as well as for possible applications to full three-dimensional Monte Carlo based fuel management modeling in the longer term. On the experimental side, several new benchmark-quality code validation measurements based on neutron activation spectrometry have been conducted using the ATRC. Results for the first four experiments, focused on neutron spectrum measurements within the Northwest Large In-Pile Tube (NW LIPT) and in the core fuel elements surrounding the NW LIPT and the diametrically opposite Southeast IPT have been reported [8,9]. A fifth, very recent, experiment focused on detailed measurements of the element-to-element core power distribution is summarized here and examples of the use of the measured data for validation of corresponding MCNP5, HELIOS, NEWT, and Serpent computational models using modern least-square adjustment methods are provided.
Inlet-Compressor Analysis Performed Using Coupled Computational Fluid Dynamics Codes
NASA Technical Reports Server (NTRS)
Cole, Gary L.; Suresh, Ambady; Townsend, Scott
1999-01-01
A thorough understanding of dynamic interactions between inlets and compressors is extremely important to the design and development of propulsion control systems, particularly for supersonic aircraft such as the High-Speed Civil Transport (HSCT). Computational fluid dynamics (CFD) codes are routinely used to analyze individual propulsion components. By coupling the appropriate CFD component codes, it is possible to investigate inlet-compressor interactions. The objectives of this work were to gain a better understanding of inlet-compressor interaction physics, formulate a more realistic compressor-face boundary condition for time-accurate CFD simulations of inlets, and to take a first step toward the CFD simulation of an entire engine by coupling multidimensional component codes. This work was conducted at the NASA Lewis Research Center by a team of civil servants and support service contractors as part of the High Performance Computing and Communications Program (HPCCP).
Safety Related Investigations of the VVER-1000 Reactor Type by the Coupled Code System TRACE/PARCS
NASA Astrophysics Data System (ADS)
Jaeger, Wadim; Espinoza, Victor Hugo Sánchez; Lischke, Wolfgang
This study was performed at the Institute of Reactor Safety at the Forschungszentrum Karlsruhe. It is embedded in the ongoing investigations of the international code assessment and maintenance program (CAMP) for qualification and validation of system codes like TRACE(1) and PARCS(2). The chosen reactor type used to validate these two codes was the Russian designed VVER-1000 because the OECD/NEA VVER-1000 Coolant Transient Benchmark Phase 2(3) includes detailed information of the Bulgarian nuclear power plant (NPP) Kozloduy unit 6. The post-test investigations of a coolant mixing experiment have shown that the predicted parameters (coolant temperature, pressure drop, etc.) are in good agreement with the measured data. The coolant mixing pattern, especially in the downcomer, has been also reproduced quiet well by TRACE. The coupled code system TRACE/PARCS which was applied on a postulated main steam line break (MSLB) provided good results compared to reference values and the ones of other participants of the benchmark. The results show that the developed three-dimensional nodalization of the reactor pressure vessel (RPV) is appropriate to describe the coolant mixing phenomena in the downcomer and the lower plenum of a VVER-1000 reactor. This phenomenon is a key issue for investigations of MSLB transient where the thermal hydraulics and the core neutronics are strongly linked. The simulation of the RPV and core behavior for postulated transients using the validated 3D TRACE RPV model, taking into account boundary conditions at vessel in- and outlet, indicates that the results are physically sound and in good agreement to other participant's results.
An introduction to LIME 1.0 and its use in coupling codes for multiphysics simulations.
Belcourt, Noel; Pawlowski, Roger Patrick; Schmidt, Rodney Cannon; Hooper, Russell Warren
2011-11-01
LIME is a small software package for creating multiphysics simulation codes. The name was formed as an acronym denoting 'Lightweight Integrating Multiphysics Environment for coupling codes.' LIME is intended to be especially useful when separate computer codes (which may be written in any standard computer language) already exist to solve different parts of a multiphysics problem. LIME provides the key high-level software (written in C++), a well defined approach (with example templates), and interface requirements to enable the assembly of multiple physics codes into a single coupled-multiphysics simulation code. In this report we introduce important software design characteristics of LIME, describe key components of a typical multiphysics application that might be created using LIME, and provide basic examples of its use - including the customized software that must be written by a user. We also describe the types of modifications that may be needed to individual physics codes in order for them to be incorporated into a LIME-based multiphysics application.
A mono-dimensional nuclear fuel performance analysis code, PUMA, development from a coupled approach
Cheon, J. S.; Lee, B. O.; Lee, C. B.; Yacout, A. M.
2013-07-01
Multidimensional-multi-physical phenomena in nuclear fuels are treated as a set of mono-dimensional-coupled problems which encompass heat, displacement, fuel constituent redistribution, and fission gas release. Rather than uncoupling these coupled equations as in conventional fuel performance analysis codes, efforts are put into to obtain fully coupled solutions by relying on the recent advances of numerical analysis. Through this approach, a new SFR metal fuel performance analysis code, called PUMA (Performance of Uranium Metal fuel rod Analysis code) is under development. Although coupling between temperature and fuel constituent was made easily, the coupling between the mechanical equilibrium equation and a set of stiff kinetics equations for fission gas release is accomplished by introducing one-level Newton scheme through backward differentiation formula. Displacement equations from 1D finite element formulation of the mechanical equilibrium equation are solved simultaneously with stress equation, creep equation, swelling equation, and FGR equations. Calculations was made successfully such that the swelling and the hydrostatic pressure are interrelated each other. (authors)
Wangerin, K; Culbertson, C N; Jevremovic, T
2005-08-01
The goal of this study was to evaluate the COG Monte Carlo radiation transport code, developed and tested by Lawrence Livermore National Laboratory, for gadolinium neutron capture therapy (GdNCT) related modeling. The validity of COG NCT model has been established for this model, and here the calculation was extended to analyze the effect of various gadolinium concentrations on dose distribution and cell-kill effect of the GdNCT modality and to determine the optimum therapeutic conditions for treating brain cancers. The computational results were compared with the widely used MCNP code. The differences between the COG and MCNP predictions were generally small and suggest that the COG code can be applied to similar research problems in NCT. Results for this study also showed that a concentration of 100 ppm gadolinium in the tumor was most beneficial when using an epithermal neutron beam. PMID:16010124
Sayer, R.O.
2003-07-29
RSAP [1] is a computer code for display and manipulation of neutron cross section data and selected SAMMY output. SAMMY [2] is a multilevel R-matrix code for fitting neutron time-of-flight cross-section data using Bayes' method. This users' guide provides documentation for the recently updated RSAP code (version 6). The code has been ported to the Linux platform, and several new features have been added, including the capability to read cross section data from ASCII pointwise ENDF files as well as double-precision PLT output from SAMMY. A number of bugs have been found and corrected, and the input formats have been improved. Input items are parsed so that items may be separated by spaces or commas.
Wang, G. B.; Wang, K.; Liu, H. G.; Li, R. D.
2013-07-01
A Monte Carlo tool RSMC (Reaction Sequence Monte Carlo) was developed to simulate deuteron/triton transportation and reaction coupled problem. The 'Forced particle production' variance reduction technique was used to improve the simulation speed, which made the secondary product play a major role. The mono-energy 14 MeV fusion neutron source was employed as a validation. Then the thermal-to-fusion neutron convertor was studied with our tool. Moreover, an in-core conversion efficiency measurement experiment was performed with {sup 6}LiD and {sup 6}LiH converters. Threshold activation foils was used to indicate the fast and fusion neutron flux. Besides, two other pivotal parameters were calculated theoretically. Finally, the conversion efficiency of {sup 6}LiD is obtained as 1.97x10{sup -4}, which matches well with the theoretical result. (authors)
Code System for Calculating Alpha, N; Spontaneous Fission; and Delayed Neutron Sources and Spectra.
Energy Science and Technology Software Center (ESTSC)
2002-07-18
Version: 04 SOURCES4C is a code system that determines neutron production rates and spectra from (alpha,n) reactions, spontaneous fission, and delayed neutron emission due to radionuclide decay. In this release the three-region problem was modified to correct several bugs, and new documentation was added to the package. Details are available in the included LA-UR-02-1617 (2002) report. The code is capable of calculating (alpha,n) source rates and spectra in four types of problems: homogeneous media (i.e.,more » an intimate mixture of alpha-emitting source material and low-Z target material), two-region interface problems (i.e., a slab of alpha-emitting source material in contact with a slab of low-Z target material), three-region interface problems (i.e., a thin slab of low-Z target material sandwiched between alpha-emitting source material and low-Z target material), and (alpha,n) reactions induced by a monoenergetic beam of alpha-particles incident on a slab of target material. The process of creating a SOURCES input file (tape1) is streamlined with the Los Alamos SOURCES Tape1 Creator and Library Link (LASTCALL) Version 1. Intended to supplement the SOURCES manual, LASTCALL is a simple graphical user interface designed to minimize common errors made during input. The optional application, LASTCALL, consists of a single dialog window launched from an executable (lastcall.exe) on Windows-based personal computers.« less
Validation of the new code package APOLLO2.8 for accurate PWR neutronics calculations
Santamarina, A.; Bernard, D.; Blaise, P.; Leconte, P.; Palau, J. M.; Roque, B.; Vaglio, C.; Vidal, J. F.
2013-07-01
This paper summarizes the Qualification work performed to demonstrate the accuracy of the new APOLLO2.S/SHEM-MOC package based on JEFF3.1.1 nuclear data file for the prediction of PWR neutronics parameters. This experimental validation is based on PWR mock-up critical experiments performed in the EOLE/MINERVE zero-power reactors and on P.I. Es on spent fuel assemblies from the French PWRs. The Calculation-Experiment comparison for the main design parameters is presented: reactivity of UOX and MOX lattices, depletion calculation and fuel inventory, reactivity loss with burnup, pin-by-pin power maps, Doppler coefficient, Moderator Temperature Coefficient, Void coefficient, UO{sub 2}-Gd{sub 2}O{sub 3} poisoning worth, Efficiency of Ag-In-Cd and B4C control rods, Reflector Saving for both standard 2-cm baffle and GEN3 advanced thick SS reflector. From this qualification process, calculation biases and associated uncertainties are derived. This code package APOLLO2.8 is already implemented in the ARCADIA new AREVA calculation chain for core physics and is currently under implementation in the future neutronics package of the French utility Electricite de France. (authors)
A coupled deterministic/stochastic method for computing neutron capture therapy dose rates
NASA Astrophysics Data System (ADS)
Hubbard, Thomas Richard
Neutron capture therapy (NCT) is an experimental method of treating brain tumors and other cancers by: (1) injecting or infusing the patient with a tumor-seeking, neutron target-labeled drug; and (2) irradiating the patient in an intense epithermal neutron fluence. The nuclear reaction between the neutrons and the target nuclei (e.g. sp{10}B(n,alpha)sp7Lirbrack releases energy in the form of high-LET (i.e. energy deposited within the range of a cell diameter) reaction particles which selectively kill the tumor cell. The efficacy of NCT is partly dependent on the delivery of maximum thermal neutron fluence to the tumor and the minimization of radiation dose to healthy tissue. Since the filtered neutron source (e.g. research reactor) usually provides a broad energy spectrum of highly-penetrating neutron and gamma-photon radiation, detailed transport calculations are necessary in order to plan treatments that use optimal treatment facility configurations and patient positioning. Current computational methods for NCT use either discrete ordinates calculation or, more often, Monte Carlo simulation to predict neutron fluences in the vicinity of the tumor. These methods do not, however, accurately calculate the transport of radiation throughout the entire facility or the deposition of dose in all the various parts of the body due to shortcomings of using either method alone. A computational method, specifically designed for NCT problems, has been adapted from the MASH methodology and couples a forward discrete ordinates (Ssb{n}) calculation with an adjoint Monte Carlo run to predict the dose at any point within the patient. The transport from the source through the filter/collimator is performed with a forward DORT run, and this is then coupled to adjoint MORSE results at a selected coupling parallelepiped which surrounds human phantom. Another routine was written to allow the user to generate the MORSE models at various angles and positions within the treatment room. The
Observational constraints on neutron star crust-core coupling during glitches
NASA Astrophysics Data System (ADS)
Newton, W. G.; Berger, S.; Haskell, B.
2015-12-01
We demonstrate that observations of glitches in the Vela pulsar can be used to investigate the strength of the crust-core coupling in a neutron star and provide a powerful probe of the internal structure of neutron stars. We assume that glitch recovery is dominated by the torque exerted by the mutual friction-mediated recoupling of superfluid components of the core that were decoupled from the crust during the glitch. Then we use the observations of the recoveries from two recent glitches in the Vela pulsar to infer the fraction of the core that is coupled to the crust during the glitch. We then analyse whether crustal neutrons alone are sufficient to drive glitches in the Vela pulsar, taking into account crustal entrainment. We use two sets of neutron star equations of state (EOSs) which span crust and core consistently and cover a conservative range of the slope of the symmetry energy at saturation density 30 < L < 120 MeV. The two sets differ in the stiffness of the high density EOS. We find that for medium to stiff EOSs, observations imply >70 per cent of the moment of inertia of the core is coupled to the crust during the glitch, though for softer EOSs L ≈ 30 MeV as little as 5 per cent could be coupled. We find that only by extending the region where superfluid vortices are strongly pinned into the core by densities at least 0.016 fm-3 above the crust-core transition density does any EOS reproduce the observed glitch activity.
Luo, Qinlong; Martins, G. B.; Daghofer, M.; Yu, Rong; Yildirim, Yucel; Moreo, Adriana; Dagotto, Elbio R
2010-01-01
The results of neutron-scattering and angle-resolved photoemission experiments for the Fe-pnictide parent compounds, and their metallic nature, are shown to impose severe constraints on the range of values that can be considered realistic for the intraorbital Hubbard repulsion U and Hund coupling J in multiorbital Hubbard models treated in the mean-field approximation. Phase diagrams for three- and five-orbital models are here provided, and the physically realistic regime of couplings is highlighted, to guide future theoretical work into the proper region of parameters of Hubbard models. In addition, using the random phase approximation, the pairing tendencies in these realistic coupling regions are investigated. It is shown that the dominant spin-singlet pairing channels in these coupling regimes correspond to nodal superconductivity, with strong competition between several states that belong to different irreducible representations. This is compatible with experimental bulk measurements that have reported the existence of nodes in several Fe-pnictide compounds.
NASA Technical Reports Server (NTRS)
Armstrong, T. W.
1972-01-01
Several Monte Carlo radiation transport computer codes are used to predict quantities of interest in the fields of radiotherapy and radiobiology. The calculational methods are described and comparisions of calculated and experimental results are presented for dose distributions produced by protons, neutrons, and negatively charged pions. Comparisons of calculated and experimental cell survival probabilities are also presented.
ITS Version 6 : the integrated TIGER series of coupled electron/photon Monte Carlo transport codes.
Franke, Brian Claude; Kensek, Ronald Patrick; Laub, Thomas William
2008-04-01
ITS is a powerful and user-friendly software package permitting state-of-the-art Monte Carlo solution of lineartime-independent coupled electron/photon radiation transport problems, with or without the presence of macroscopic electric and magnetic fields of arbitrary spatial dependence. Our goal has been to simultaneously maximize operational simplicity and physical accuracy. Through a set of preprocessor directives, the user selects one of the many ITS codes. The ease with which the makefile system is applied combines with an input scheme based on order-independent descriptive keywords that makes maximum use of defaults and internal error checking to provide experimentalists and theorists alike with a method for the routine but rigorous solution of sophisticated radiation transport problems. Physical rigor is provided by employing accurate cross sections, sampling distributions, and physical models for describing the production and transport of the electron/photon cascade from 1.0 GeV down to 1.0 keV. The availability of source code permits the more sophisticated user to tailor the codes to specific applications and to extend the capabilities of the codes to more complex applications. Version 6, the latest version of ITS, contains (1) improvements to the ITS 5.0 codes, and (2) conversion to Fortran 90. The general user friendliness of the software has been enhanced through memory allocation to reduce the need for users to modify and recompile the code.
NASA Astrophysics Data System (ADS)
Chauvet, Yves
1985-07-01
This paper summarized two improvements of a real production code by using vectorization and multitasking techniques. After a short description of Monte Carlo algorithms employed in our neutron transport problems, we briefly describe the work we have done in order to get a vector code. Vectorization principles will be presented and measured performances on the CRAY 1S, CYBER 205 and CRAY X-MP compared in terms of vector lengths. The second part of this work is an adaptation to multitasking on the CRAY X-MP using exclusively standard multitasking tools available with FORTRAN under the COS 1.13 system. Two examples will be presented. The goal of the first one is to measure the overhead inherent to multitasking when tasks become too small and to define a granularity threshold that is to say a minimum size for a task. With the second example we propose a method that is very X-MP oriented in order to get the best speedup factor on such a computer. In conclusion we prove that Monte Carlo algorithms are very well suited to future vector and parallel computers.
A Visualization Code System for Gamma and Neutron Shielding Calculations, Version 2.0
Energy Science and Technology Software Center (ESTSC)
2008-08-01
EASYQAD, Version 2.0, is a standalone Windows XP or Windows 7 code system which facilitates gamma and neutron shielding calculations with user friendly graphical interfaces. It is used to analyze radiation shielding problems and includes: - 8 kinds of geometry types - Various flexible source options - Common material library - Various detector types The update contents of EASYQAD Version 2.0 are below: - Addition of starting option with P-code files - Addition of multi-sourcemore » calculation function - Expansion of source geometries - Addition of warning message - Modifications of EASYQAD program errors a. Coordination application problem in source division b. Source position error c. Rotation problem of source geometry d. Program running error in using more than six gamma energy distribution e. EASYQAD display problem of the right elliptic cylinder, ellipsoid and truncated right cone geometries Through intuitive windows and their interactions inside EASYQAD, the user can specify the dimensions of 3D-shapes, their material compositions, their densities, the type of radioactive sources, the locations of the sources, the type and positions of detectors. With the ease of using these sequences, shielding problems will become simpler and more clearly understandable to the analyzer. Furthermore, the error checking system can prevent users from making mistakes by automatically debugging the user inputs and giving modal dialog windows. The included AECL implementation of QAD-CGGP-A, Version 95.2 (C00645MNYCP00), is run from the user interface.« less
The TORT three-dimensional discrete ordinates neutron/photon transport code (TORT version 3)
Rhoades, W.A.; Simpson, D.B.
1997-10-01
TORT calculates the flux or fluence of neutrons and/or photons throughout three-dimensional systems due to particles incident upon the system`s external boundaries, due to fixed internal sources, or due to sources generated by interaction with the system materials. The transport process is represented by the Boltzman transport equation. The method of discrete ordinates is used to treat the directional variable, and a multigroup formulation treats the energy dependence. Anisotropic scattering is treated using a Legendre expansion. Various methods are used to treat spatial dependence, including nodal and characteristic procedures that have been especially adapted to resist numerical distortion. A method of body overlay assists in material zone specification, or the specification can be generated by an external code supplied by the user. Several special features are designed to concentrate machine resources where they are most needed. The directional quadrature and Legendre expansion can vary with energy group. A discontinuous mesh capability has been shown to reduce the size of large problems by a factor of roughly three in some cases. The emphasis in this code is a robust, adaptable application of time-tested methods, together with a few well-tested extensions.
Use of SUSA in Uncertainty and Sensitivity Analysis for INL VHTR Coupled Codes
Gerhard Strydom
2010-06-01
The need for a defendable and systematic Uncertainty and Sensitivity approach that conforms to the Code Scaling, Applicability, and Uncertainty (CSAU) process, and that could be used for a wide variety of software codes, was defined in 2008.The GRS (Gesellschaft für Anlagen und Reaktorsicherheit) company of Germany has developed one type of CSAU approach that is particularly well suited for legacy coupled core analysis codes, and a trial version of their commercial software product SUSA (Software for Uncertainty and Sensitivity Analyses) was acquired on May 12, 2010. This interim milestone report provides an overview of the current status of the implementation and testing of SUSA at the INL VHTR Project Office.
Burns, Kimberly A.
2009-08-01
The accurate and efficient simulation of coupled neutron-photon problems is necessary for several important radiation detection applications. Examples include the detection of nuclear threats concealed in cargo containers and prompt gamma neutron activation analysis for nondestructive determination of elemental composition of unknown samples.
Sartori, E.
1992-12-31
This paper presents a brief review of computer codes concerned with checking, plotting, processing and using of covariances of neutron cross-section data. It concentrates on those available from the computer code information centers of the United States and the OECD/Nuclear Energy Agency. Emphasis will be placed also on codes using covariances for specific applications such as uncertainty analysis, data adjustment and data consistency analysis. Recent evaluations contain neutron cross section covariance information for all isotopes of major importance for technological applications of nuclear energy. It is therefore important that the available software tools needed for taking advantage of this information are widely known as hey permit the determination of better safety margins and allow the optimization of more economic, I designs of nuclear energy systems.
Studies of coupled cavity LINAC (CCL) accelerating structures with 3-D codes
Spalek, G.
2000-08-01
The cw CCL being designed for the Accelerator Production of Tritium (APT) project accelerates protons from 96 MeV to 211 MeV. It consists of 99 segments each containing up to seven accelerating cavities. Segments are coupled by intersegment coupling cavities and grouped into supermodules. The design method needs to address not only basic cavity sizing for a given coupling and pi/2 mode frequency, but also the effects of high power densities on the cavity frequency, mechanical stresses, and the structure's stop band during operation. On the APT project, 3-D RF (Ansoft Corp.'s HFSS) and coupled RF/structural (Ansys Inc.'s ANSYS) codes are being used. to develop tools to address the above issues and guide cooling channel design. The code's predictions are being checked against available low power Aluminum models. Stop band behavior under power will be checked once the tools are extended to CCDTL structures that have been tested at high power. A summary of calculations made to date and agreement with measured results will be presented.
Barths, H.; Felsch, C.; Peters, N.
2008-11-15
The objective of this work is the development of a consistent mixing model for the two-way-coupling of a CFD code and a multi-zone code based on multiple zero-dimensional reactors. The two-way-coupling allows for a computationally efficient modeling of HCCI combustion. The physical domain in the CFD code is subdivided into multiple zones based on three phase variables (fuel mixture fraction, dilution, and total enthalpy). Those phase variables are sufficient for the description of the thermodynamic state of each zone, assuming that each zone is at the same pressure. Each zone in the CFD code is represented by a corresponding zone in the zero-dimensional code. The zero-dimensional code solves the chemistry for each zone, and the heat release is fed back into the CFD code. The difficulty in facing this kind of methodology is to keep the thermodynamic state of each zone consistent between the CFD code and the zero-dimensional code after the initialization of the zones in the multi-zone code has taken place. The thermodynamic state of each zone (and thereby the phase variables) will change in time due to mixing and source terms (e.g., vaporization of fuel, wall heat transfer). The focus of this work lies on a consistent description of the mixing between the zones in phase space in the zero-dimensional code, based on the solution of the CFD code. Two mixing models with different degrees of accuracy, complexity, and numerical effort are described. The most elaborate mixing model (and an appropriate treatment of the source terms) keeps the thermodynamic state of the zones in the CFD code and the zero-dimensional code identical. The models are applied to a test case of HCCI combustion in an engine. (author)
Barths, H.; Felsch, C.; Peters, N.
2009-01-15
The objective of this work is the development of a consistent mixing model for the two-way-coupling of a CFD code and a multi-zone code based on multiple zero-dimensional reactors. The two-way-coupling allows for a computationally efficient modeling of HCCI combustion. The physical domain in the CFD code is subdivided into multiple zones based on three phase variables (fuel mixture fraction, dilution, and total enthalpy). Those phase variables are sufficient for the description of the thermodynamic state of each zone, assuming that each zone is at the same pressure. Each zone in the CFD code is represented by a corresponding zone in the zero-dimensional code. The zero-dimensional code solves the chemistry for each zone, and the heat release is fed back into the CFD code. The difficulty in facing this kind of methodology is to keep the thermodynamic state of each zone consistent between the CFD code and the zero-dimensional code after the initialization of the zones in the multi-zone code has taken place. The thermodynamic state of each zone (and thereby the phase variables) will change in time due to mixing and source terms (e.g., vaporization of fuel, wall heat transfer). The focus of this work lies on a consistent description of the mixing between the zones in phase space in the zero-dimensional code, based on the solution of the CFD code. Two mixing models with different degrees of accuracy, complexity, and numerical effort are described. The most elaborate mixing model (and an appropriate treatment of the source terms) keeps the thermodynamic state of the zones in the CFD code and the zero-dimensional code identical. The models are applied to a test case of HCCI combustion in an engine. (author)
Chang H. Oh; Eung S. Kim; Mike Patterson
2010-06-01
Abstract – A tritium permeation analyses code (TPAC) was developed by Idaho National Laboratory for the purpose of analyzing tritium distributions in very high temperature reactor (VHTR) systems, including integrated hydrogen production systems. A MATLAB SIMULINK software package was used in developing the code. The TPAC is based on the mass balance equations of tritium-containing species and various forms of hydrogen coupled with a variety of tritium sources, sinks, and permeation models. In the TPAC, ternary fission and neutron reactions with 6Li, 7Li 10B, and 3He were taken into considerations as tritium sources. Purification and leakage models were implemented as main tritium sinks. Permeation of tritium and H2 through pipes, vessels, and heat exchangers were considered as main tritium transport paths. In addition, electroyzer and isotope exchange models were developed for analyzing hydrogen production systems, including high temperature electrolysis and sulfur-iodine processes.
MULTIDIMENSIONAL COUPLED PHOTON-ELECTRON TRANSPORT SIMULATIONS USING NEUTRAL PARTICLE SN CODES
Ilas, Dan; Williams, Mark L; Peplow, Douglas E.; Kirk, Bernadette Lugue
2008-01-01
During the past two years a study was underway at ORNL to assess the suitability of the popular SN neutral particle codes ANISN, DORT and TORT for coupled photon-electron calculations specific to external beam therapy of medical physics applications. The CEPXS-BFP code was used to generate the cross sections. The computational tests were performed on phantoms typical of those used in medical physics for external beam therapy, with materials simulated by water at different densities and the comparisons were made against Monte Carlo simulations that served as benchmarks. Although the results for one-dimensional calculations were encouraging, it appeared that the higher dimensional transport codes had fundamental difficulties in handling the electron transport. The results of two-dimensional simulations using the code DORT with an S16 fully symmetric quadrature set agree fairly with the reference Monte Carlo results but not well enough for clinical applications. While the photon fluxes are in better agreement (generally, within less than 5% from the reference), the discrepancy increases, sometimes very significantly, for the electron fluxes. The paper, however, focuses on the results obtained with the three-dimensional code TORT which had convergence difficulties for the electron groups. Numerical instabilities occurred in these groups. These instabilities were more pronounced with the degree of anisotropy of the problem.
The stability of tidally deformed neutron stars to three- and four-mode coupling
Venumadhav, Tejaswi; Zimmerman, Aaron; Hirata, Christopher M.
2014-01-20
It has recently been suggested that the tidal deformation of a neutron star excites daughter p- and g-modes to large amplitudes via a quasi-static instability. This would remove energy from the tidal bulge, resulting in dissipation and possibly affecting the phase evolution of inspiralling binary neutron stars and hence the extraction of binary parameters from gravitational wave observations. This instability appears to arise because of a large three-mode interaction among the tidal mode and high-order p- and g-modes of similar radial wavenumber. We show that additional four-mode interactions enter into the analysis at the same order as the three-mode terms previously considered. We compute these four-mode couplings by finding a volume-preserving coordinate transformation that relates the energy of a tidally deformed star to that of a radially perturbed spherical star. Using this method, we relate the four-mode coupling to three-mode couplings and show that there is a near-exact cancellation between the destabilizing effect of the three-mode interactions and the stabilizing effect of the four-mode interaction. We then show that the equilibrium tide is stable against the quasi-static decay into daughter p- and g-modes to leading order. The leading deviation from the quasi-static approximation due to orbital motion of the binary is considered; while it may slightly spoil the near-cancellation, any resulting instability timescale is at least of order the gravitational wave inspiral time. We conclude that the p-/g-mode coupling does not lead to a quasi-static instability, and does not impact the phase evolution of gravitational waves from binary neutron stars.
A new three-tier architecture design for multi-sphere neutron spectrometer with the FLUKA code
NASA Astrophysics Data System (ADS)
Huang, Hong; Yang, Jian-Bo; Tuo, Xian-Guo; Liu, Zhi; Wang, Qi-Biao; Wang, Xu
2016-07-01
The current commercially, available Bonner sphere neutron spectrometer (BSS) has high sensitivity to neutrons below 20 MeV, which causes it to be poorly placed to measure neutrons ranging from a few MeV to 100 MeV. The paper added moderator layers and the auxiliary material layer upon 3He proportional counters with FLUKA code, with a view to improve. The results showed that the responsive peaks to neutrons below 20 MeV gradually shift to higher energy region and decrease slightly with the increasing moderator thickness. On the contrary, the response for neutrons above 20 MeV was always very low until we embed auxiliary materials such as copper (Cu), lead (Pb), tungsten (W) into moderator layers. This paper chose the most suitable auxiliary material Pb to design a three-tier architecture multi-sphere neutron spectrometer (NBSS). Through calculating and comparing, the NBSS was advantageous in terms of response for 5-100 MeV and the highest response was 35.2 times the response of polyethylene (PE) ball with the same PE thickness.
Radiation Coupling with the FUN3D Unstructured-Grid CFD Code
NASA Technical Reports Server (NTRS)
Wood, William A.
2012-01-01
The HARA radiation code is fully-coupled to the FUN3D unstructured-grid CFD code for the purpose of simulating high-energy hypersonic flows. The radiation energy source terms and surface heat transfer, under the tangent slab approximation, are included within the fluid dynamic ow solver. The Fire II flight test, at the Mach-31 1643-second trajectory point, is used as a demonstration case. Comparisons are made with an existing structured-grid capability, the LAURA/HARA coupling. The radiative surface heat transfer rates from the present approach match the benchmark values within 6%. Although radiation coupling is the focus of the present work, convective surface heat transfer rates are also reported, and are seen to vary depending upon the choice of mesh connectivity and FUN3D ux reconstruction algorithm. On a tetrahedral-element mesh the convective heating matches the benchmark at the stagnation point, but under-predicts by 15% on the Fire II shoulder. Conversely, on a mixed-element mesh the convective heating over-predicts at the stagnation point by 20%, but matches the benchmark away from the stagnation region.
Integrated TIGER Series of Coupled Electron/Photon Monte Carlo Transport Codes System.
Energy Science and Technology Software Center (ESTSC)
2012-11-30
Version: 00 Distribution is restricted to US Government Agencies and Their Contractors Only. The Integrated Tiger Series (ITS) is a powerful and user-friendly software package permitting state-of-the-art Monte Carlo solution of linear time-independent coupled electron/photon radiation transport problems, with or without the presence of macroscopic electric and magnetic fields of arbitrary spatial dependence. The goal has been to simultaneously maximize operational simplicity and physical accuracy. Through a set of preprocessor directives, the user selects onemore » of the many ITS codes. The ease with which the makefile system is applied combines with an input scheme based on order-independent descriptive keywords that makes maximum use of defaults and internal error checking to provide experimentalists and theorists alike with a method for the routine but rigorous solution of sophisticated radiation transport problems. Physical rigor is provided by employing accurate cross sections, sampling distributions, and physical models for describing the production and transport of the electron/photon cascade from 1.0 GeV down to 1.0 keV. The availability of source code permits the more sophisticated user to tailor the codes to specific applications and to extend the capabilities of the codes to more complex applications. Version 6, the latest version of ITS, contains (1) improvements to the ITS 5.0 codes, and (2) conversion to Fortran 95. The general user friendliness of the software has been enhanced through memory allocation to reduce the need for users to modify and recompile the code.« less
Integrated TIGER Series of Coupled Electron/Photon Monte Carlo Transport Codes System.
VALDEZ, GREG D.
2012-11-30
Version: 00 Distribution is restricted to US Government Agencies and Their Contractors Only. The Integrated Tiger Series (ITS) is a powerful and user-friendly software package permitting state-of-the-art Monte Carlo solution of linear time-independent coupled electron/photon radiation transport problems, with or without the presence of macroscopic electric and magnetic fields of arbitrary spatial dependence. The goal has been to simultaneously maximize operational simplicity and physical accuracy. Through a set of preprocessor directives, the user selects one of the many ITS codes. The ease with which the makefile system is applied combines with an input scheme based on order-independent descriptive keywords that makes maximum use of defaults and internal error checking to provide experimentalists and theorists alike with a method for the routine but rigorous solution of sophisticated radiation transport problems. Physical rigor is provided by employing accurate cross sections, sampling distributions, and physical models for describing the production and transport of the electron/photon cascade from 1.0 GeV down to 1.0 keV. The availability of source code permits the more sophisticated user to tailor the codes to specific applications and to extend the capabilities of the codes to more complex applications. Version 6, the latest version of ITS, contains (1) improvements to the ITS 5.0 codes, and (2) conversion to Fortran 95. The general user friendliness of the software has been enhanced through memory allocation to reduce the need for users to modify and recompile the code.
Monte Carlo Code System for Calculation of Multiple Scattering of Neutrons in the Resonance Region.
Energy Science and Technology Software Center (ESTSC)
1983-01-25
Version 00 MCRTOF systematically calculates capture and scattering probabilities for neutrons incident on a material disk, with neutron cross sections calculated from the resonance parameters. Capture, front and rear face scattering, transmission, etc., probabilities are obtained from the average destinations of the incident neutrons.
Keney, G.S.
1981-08-01
A computer code has been written to calculate neutron induced activation of neutral-beam injector components and the corresponding dose rates as a function of geometry, component composition, and time after shutdown. The code, ACDOS1, was written in FORTRAN IV to calculate both activity and dose rates for up to 30 target nuclides and 50 neutron groups. Sufficient versatility has also been incorporated into the code to make it applicable to a variety of general activation problems due to neutrons of energy less than 20 MeV.
A multilevel method for coupling the neutron kinetics and heat transfer equations
Tamang, A.; Anistratov, D. Y.
2013-07-01
We present a computational method for adequate and efficient coupling the neutron transport equation with the precursor and heat transfer equations. It is based on the multilevel nonlinear quasidiffusion (QD) method for solving the multigroup transport equation. The system of equations includes the time-dependent high-order transport equation and time-dependent multigroup and effective one-group low-order QD equations. We also formulate a method applying the {alpha}-approximation for the time-dependent high-order transport equation. This approach enables one to avoid storing the angular flux from the previous time step. Numerical results for a model transient problem are presented. (authors)
Pedestal Fueling Simulations with a Coupled Kinetic-kinetic Plasma-neutral Transport Code
D.P. Stotler, C.S. Chang, S.H. Ku, J. Lang and G.Y. Park
2012-08-29
A Monte Carlo neutral transport routine, based on DEGAS2, has been coupled to the guiding center ion-electron-neutral neoclassical PIC code XGC0 to provide a realistic treatment of neutral atoms and molecules in the tokamak edge plasma. The DEGAS2 routine allows detailed atomic physics and plasma-material interaction processes to be incorporated into these simulations. The spatial pro le of the neutral particle source used in the DEGAS2 routine is determined from the uxes of XGC0 ions to the material surfaces. The kinetic-kinetic plasma-neutral transport capability is demonstrated with example pedestal fueling simulations.
Energy Conservation Tests of a Coupled Kinetic-kinetic Plasma-neutral Transport Code
Stotler, D. P.; Chang, C. S.; Ku, S. H.; Lang, J.; Park, G.
2012-08-29
A Monte Carlo neutral transport routine, based on DEGAS2, has been coupled to the guiding center ion-electron-neutral neoclassical PIC code XGC0 to provide a realistic treatment of neutral atoms and molecules in the tokamak edge plasma. The DEGAS2 routine allows detailed atomic physics and plasma-material interaction processes to be incorporated into these simulations. The spatial pro le of the neutral particle source used in the DEGAS2 routine is determined from the uxes of XGC0 ions to the material surfaces. The kinetic-kinetic plasma-neutral transport capability is demonstrated with example pedestal fueling simulations.
Chatani, K. )
1992-08-01
This report summarizes the calculational results from analyses of a Clinch River Breeder Reactor (CRBR) prototypic coolant pipe chaseway neutron streaming experiment Comparisons of calculated and measured results are presented, major emphasis being placed on results at bends in the chaseway. Calculations were performed with three three-dimensional radiation transport codes: the discrete ordinates code TORT and the Monte Carlo code MORSE, both developed by the Oak Ridge National Laboratory (ORNL), and the discrete ordinates code ENSEMBLE, developed by Japan. The calculated results from the three codes are compared (1) with previously-calculated DOT3.5 two-dimensional results, (2) among themselves, and (3) with measured results. Calculations with TORT used both the weighted-difference and nodal methods. Only the weighted-difference method was used in ENSEMBLE. When the calculated results were compared to measured results, it was found that calculation-to-experiment (C/E) ratios were good in the regions of the chaseway where two-dimensional modeling might be difficult and where there were no significant discrete ordinates ray effects. Excellent agreement was observed for responses dominated by thermal neutron contributions. MORSE-calculated results and comparisons are described also, and detailed results are presented in an appendix.
NASA Astrophysics Data System (ADS)
Pica, G.; Lovett, B. W.; Bhatt, R. N.; Schenkel, T.; Lyon, S. A.
2016-01-01
A scaled quantum computer with donor spins in silicon would benefit from a viable semiconductor framework and a strong inherent decoupling of the qubits from the noisy environment. Coupling neighboring spins via the natural exchange interaction according to current designs requires gate control structures with extremely small length scales. We present a silicon architecture where bismuth donors with long coherence times are coupled to electrons that can shuttle between adjacent quantum dots, thus relaxing the pitch requirements and allowing space between donors for classical control devices. An adiabatic SWAP operation within each donor/dot pair solves the scalability issues intrinsic to exchange-based two-qubit gates, as it does not rely on subnanometer precision in donor placement and is robust against noise in the control fields. We use this SWAP together with well established global microwave Rabi pulses and parallel electron shuttling to construct a surface code that needs minimal, feasible local control.
Wilson, W. B.; Perry, R. T.; Shores, E. F.; Charlton, W. S.; Parish, Theodore A.; Estes, G. P.; Brown, T. H.; Arthur, Edward D. ,; Bozoian, Michael; England, T. R.; Madland, D. G.; Stewart, J. E.
2002-01-01
SOURCES 4C is a computer code that determines neutron production rates and spectra from ({alpha},n) reactions, spontaneous fission, and delayed neutron emission due to radionuclide decay. The code is capable of calculating ({alpha},n) source rates and spectra in four types of problems: homogeneous media (i.e., an intimate mixture of a-emitting source material and low-Z target material), two-region interface problems (i.e., a slab of {alpha}-emitting source material in contact with a slab of low-Z target material), three-region interface problems (i.e., a thin slab of low-Z target material sandwiched between {alpha}-emitting source material and low-Z target material), and ({alpha},n) reactions induced by a monoenergetic beam of {alpha}-particles incident on a slab of target material. Spontaneous fission spectra are calculated with evaluated half-life, spontaneous fission branching, and Watt spectrum parameters for 44 actinides. The ({alpha},n) spectra are calculated using an assumed isotropic angular distribution in the center-of-mass system with a library of 107 nuclide decay {alpha}-particle spectra, 24 sets of measured and/or evaluated ({alpha},n) cross sections and product nuclide level branching fractions, and functional {alpha}-particle stopping cross sections for Z < 106. The delayed neutron spectra are taken from an evaluated library of 105 precursors. The code provides the magnitude and spectra, if desired, of the resultant neutron source in addition to an analysis of the'contributions by each nuclide in the problem. LASTCALL, a graphical user interface, is included in the code package.
NASA Astrophysics Data System (ADS)
Jiang, Li; Zhang, Guo-Feng
2016-08-01
The effects of nuclear field and spin-orbit interaction on dense coding and swap operation are studied in detail for both the antiferromagnetic (AFM) and ferromagnetic (FM) coupling cases. The conditions for a valid dense coding and under which swap operation is feasible are given.
MAC/GMC Code Enhanced for Coupled Electromagnetothermoelastic Analysis of Smart Composites
NASA Technical Reports Server (NTRS)
Bednarcyk, Brett A.; Arnold, Steven M.; Aboudi, Jacob
2002-01-01
Intelligent materials are those that exhibit coupling between their electromagnetic response and their thermomechanical response. This coupling allows smart materials to react mechanically (e.g., an induced displacement) to applied electrical or magnetic fields (for instance). These materials find many important applications in sensors, actuators, and transducers. Recently interest has arisen in the development of smart composites that are formed via the combination of two or more phases, one or more of which is a smart material. To design with and utilize smart composites, designers need theories that predict the coupled smart behavior of these materials from the electromagnetothermoelastic properties of the individual phases. The micromechanics model known as the generalized method of cells (GMC) has recently been extended to provide this important capability. This coupled electromagnetothermoelastic theory has recently been incorporated within NASA Glenn Research Center's Micromechanics Analysis Code with Generalized Method of Cells (MAC/GMC). This software package is user friendly and has many additional features that render it useful as a design and analysis tool for composite materials in general, and with its new capabilities, for smart composites as well.
Simulating Coupling Complexity in Space Plasmas: First Results from a new code
NASA Astrophysics Data System (ADS)
Kryukov, I.; Zank, G. P.; Pogorelov, N. V.; Raeder, J.; Ciardo, G.; Florinski, V. A.; Heerikhuisen, J.; Li, G.; Petrini, F.; Shematovich, V. I.; Winske, D.; Shaikh, D.; Webb, G. M.; Yee, H. M.
2005-12-01
The development of codes that embrace 'coupling complexity' via the self-consistent incorporation of multiple physical scales and multiple physical processes in models has been identified by the NRC Decadal Survey in Solar and Space Physics as a crucial necessary development in simulation/modeling technology for the coming decade. The National Science Foundation, through its Information Technology Research (ITR) Program, is supporting our efforts to develop a new class of computational code for plasmas and neutral gases that integrates multiple scales and multiple physical processes and descriptions. We are developing a highly modular, parallelized, scalable code that incorporates multiple scales by synthesizing 3 simulation technologies: 1) Computational fluid dynamics (hydrodynamics or magneto-hydrodynamics-MHD) for the large-scale plasma; 2) direct Monte Carlo simulation of atoms/neutral gas, and 3) transport code solvers to model highly energetic particle distributions. We are constructing the code so that a fourth simulation technology, hybrid simulations for microscale structures and particle distributions, can be incorporated in future work, but for the present, this aspect will be addressed at a test-particle level. This synthesis we will provide a computational tool that will advance our understanding of the physics of neutral and charged gases enormously. Besides making major advances in basic plasma physics and neutral gas problems, this project will address 3 Grand Challenge space physics problems that reflect our research interests: 1) To develop a temporal global heliospheric model which includes the interaction of solar and interstellar plasma with neutral populations (hydrogen, helium, etc., and dust), test-particle kinetic pickup ion acceleration at the termination shock, anomalous cosmic ray production, interaction with galactic cosmic rays, while incorporating the time variability of the solar wind and the solar cycle. 2) To develop a coronal
A Coupling Methodology for Mesoscale-informed Nuclear Fuel Performance Codes
Michael Tonks; Derek Gaston; Cody Permann; Paul Millett; Glen Hansen; Dieter Wolf
2010-10-01
This study proposes an approach for capturing the effect of microstructural evolution on reactor fuel performance by coupling a mesoscale irradiated microstructure model with a finite element fuel performance code. To achieve this, the macroscale system is solved in a parallel, fully coupled, fully-implicit manner using the preconditioned Jacobian-free Newton Krylov (JFNK) method. Within the JFNK solution algorithm, microstructure-influenced material parameters are calculated by the mesoscale model and passed back to the macroscale calculation. Due to the stochastic nature of the mesoscale model, a dynamic fitting technique is implemented to smooth roughness in the calculated material parameters. The proposed methodology is demonstrated on a simple model of a reactor fuel pellet. In the model, INL’s BISON fuel performance code calculates the steady-state temperature profile in a fuel pellet and the microstructure-influenced thermal conductivity is determined with a phase field model of irradiated microstructures. This simple multiscale model demonstrates good nonlinear convergence and near ideal parallel scalability. By capturing the formation of large mesoscale voids in the pellet interior, the multiscale model predicted the irradiation-induced reduction in the thermal conductivity commonly observed in reactors.
Phase-amplitude coupling supports phase coding in human ECoG
Watrous, Andrew J; Deuker, Lorena; Fell, Juergen; Axmacher, Nikolai
2015-01-01
Prior studies have shown that high-frequency activity (HFA) is modulated by the phase of low-frequency activity. This phenomenon of phase-amplitude coupling (PAC) is often interpreted as reflecting phase coding of neural representations, although evidence for this link is still lacking in humans. Here, we show that PAC indeed supports phase-dependent stimulus representations for categories. Six patients with medication-resistant epilepsy viewed images of faces, tools, houses, and scenes during simultaneous acquisition of intracranial recordings. Analyzing 167 electrodes, we observed PAC at 43% of electrodes. Further inspection of PAC revealed that category specific HFA modulations occurred at different phases and frequencies of the underlying low-frequency rhythm, permitting decoding of categorical information using the phase at which HFA events occurred. These results provide evidence for categorical phase-coded neural representations and are the first to show that PAC coincides with phase-dependent coding in the human brain. DOI: http://dx.doi.org/10.7554/eLife.07886.001 PMID:26308582
Cold Uniform Matter and Neutron Stars in the Quark-Meson-Coupling Model
J.R. Stone; P.A.M. Guichon; H.H. Matevosyan; A.W. Thomas
2007-08-01
A new density dependent effective baryon-baryon interaction has been recently derived from the quark-meson-coupling (QMC) model, offering impressive results in application to finite nuclei and dense baryon matter. This self-consistent, relativistic quark-level approach is used to construct the Equation of State (EoS) and to calculate key properties of high density matter and cold, slowly rotating neutron stars. The results include predictions for the maximum mass of neutron star models, together with the corresponding radius and central density, as well the properties of neutron stars with mass of order 1.4 M{sub {circle_dot}}. The cooling mechanism allowed by the QMC EoS is explored and the parameters relevant to slow rotation, namely the moment of inertia and the period of rotation investigated. The results of the calculation, which are found to be in good agreement with available observational data, are compared with the predictions of more traditional EoS, based on the A18+{delta}v+UIX* and modified Reid soft core potentials, the Skyrme SkM* interaction and two relativistic mean field (RMF) models for a hybrid stars including quark matter. The QMC EoS provides cold neutron star models with maximum mass 1.9-2.1 M{sub {circle_dot}}, with central density less than 6 times nuclear saturation density (n{sub 0} = 0.16 fm{sup -3}) and offers a consistent description of the stellar mass up to this density limit. In contrast with other models, QMC predicts no hyperon contribution at densities lower than 3n{sub 0}, for matter in {beta}-equilibrium. At higher densities, {Xi}{sup -,0} and {Lambda} hyperons are present. The absence of lighter {Sigma}{sup {+-},0} hyperons is understood as a consequence of antisymmetrization, together with the implementation of the color hyperfine interaction in the response of the quark bag to the nuclear scalar field.
Characterisation of the TRIUMF neutron facility using a Monte Carlo simulation code.
Monk, S D; Abram, T; Joyce, M J
2015-04-01
Here, the characterisation of the high-energy neutron field at TRIUMF (The Tri Universities Meson Facility, Vancouver, British Columbia) with Monte Carlo simulation software is described. The package used is MCNPX version 2.6.0, with the neutron fluence rate determined at three locations within the TRIUMF Thermal Neutron Facility (TNF), including the exit of the neutron channel where users of the facility can test devices that may be susceptible to the effects of this form of radiation. The facility is often used to roughly emulate the field likely to be encountered at high altitudes due to radiation of galactic origin and thus the simulated information is compared with the energy spectrum calculated to be due to neutron radiation of cosmic origin at typical aircraft altitudes. The calculated values were also compared with neutron flux measurements that were estimated using the activation of various foils by the staff of the facility, showing agreement within an order of magnitude. PMID:25342608
NASA Astrophysics Data System (ADS)
Feltus, Madeline Anne
1990-01-01
This thesis examines coupled time-dependent thermal -hydraulic (T/H) and neutronics solution methods for Pressurized Water Reactor (PWR) transient analysis. The degree of equivalence is evaluated between the typical quasi-static approach and a newly-developed iterative tandem method. Four specific PWR transients that exhibit a wide range of Reactor Coolant System (RCS) T/H response were investigated: (1) a Station Blackout Anticipated Transient Without Scram (ATWS), (2) a Loss of Feedwater ATWS, (3) a Total Loss of RCS Flow with Scram, and (4) a Main Steam Line Break (MSLB). Rather than using simplified RCS and core models, the theory and method in this thesis were applied practically by using realistic models for an actual four-loop Westinghouse PWR plant. The time-dependent STAR kinetics code, based on the QUANDRY Analytic Nodal Method, and the RETRAN and MCPWR T/H systems codes were used to develop a new, fully coupled, tandem STAR/MCPWRQ methodology that runs tandemly on an enhanced 386/387 IBM PC architecture. MCPWRQ uses externally calculated power input rather than point kinetics power level results. The tandem method was compared to quasi -static STAR and time-dependent STAR 2-D and 3-D kinetics results. The new STAR/MCPWRQ method uses RETRAN time-dependent T/H and point kinetics power input as a first estimate. STAR and MCPWRQ are used tandemly to couple STAR 3-D, time-dependent core power results with the MCPWRQ RCS T/H phenomena. This thesis shows that: (a) quasi-static and point kinetics methods are not able to describe severe PWR transient phenomena adequately; and (b) fully coupled, 3-D, time -dependent, tandem (or possibly parallel) analysis methods should be used for PWR reactor transients instead. By tandemly coupling the RCS response in terms of updated core inlet conditions with 3-D time-dependent core kinetics response, the core power response and T/H conditions are forced to be self-consistent during the entire transient. The transient analyses
NASA Astrophysics Data System (ADS)
Clarke, S. D.; Miller, E. C.; Flaska, M.; Pozzi, S. A.; Oberer, R. B.; Chiang, L. G.
2013-02-01
Neutron coincidence counting is widely used in nuclear safeguards. Simulations of these systems can be performed using Monte Carlo codes such as MCNPX to aid in calibration or measurement design. However, the MCNPX coincidence-counting routines treat particle histories individually, therefore the dead time of the acquisition electronics is not treated. The MCNPX-PoliMi code provides the ability to model detailed effects such as data-acquisition electronics and system dead times. A specialized post-processing code has been developed to interpret the collision-log file and determine the response of a 3He multiplicity counter. The MCNPX-PoliMi simulation provides the full neutron multiplicity distribution measured by the 3He tubes. This distribution is used to compute the singles, doubles, and triples rates which are the quantities used to determine 235U mass. MCNPX-PoliMi has previously been validated with passive multiplicity measurements. In this study, a detailed analysis of the measurement system operating in active mode is presented for uranium-oxide standards ranging from 0.5 to 4.0 kg with a Canberra JCC-51 active well coincidence counter. MCNPX-PoliMi calculations are also compared with MCNPX. The two codes agree to within 1% for the cases with negligible dead times. The simulations are validated with measurements performed at the Y-12 National Security Complex.
PHISICS multi-group transport neutronic capabilities for RELAP5
Epiney, A.; Rabiti, C.; Alfonsi, A.; Wang, Y.; Cogliati, J.; Strydom, G.
2012-07-01
PHISICS is a neutronic code system currently under development at INL. Its goal is to provide state of the art simulation capability to reactor designers. This paper reports on the effort of coupling this package to the thermal hydraulic system code RELAP5. This will enable full prismatic core and system modeling and the possibility to model coupled (thermal-hydraulics and neutronics) problems with more options for 3D neutron kinetics, compared to the existing diffusion theory neutron kinetics module in RELAP5 (NESTLE). The paper describes the capabilities of the coupling and illustrates them with a set of sample problems. (authors)
Coupling of Sph and Finite Element Codes for Multi-Layer Orbital Debris Shield Design
NASA Technical Reports Server (NTRS)
Fahrenthold, Eric P.
1997-01-01
Particle-based hydrodynamics models offer distinct advantages over Eulerian and Lagrangian hydrocodes in particular shock physics applications. Particle models are designed to avoid the mesh distortion and state variable diffusion problems which can hinder the effective use of Lagrangian and Eulerian codes respectively. However conventional particle-in-cell and smooth particle hydrodynamics methods employ particles which are actually moving interpolation points. A new particle-based modeling methodology, termed Hamiltonian particle hydrodynamics, was developed by Fahrenthold and Koo (1997) to provide an alternative, fully Lagrangian, energy-based approach to shock physics simulations. This alternative formulation avoids the tensile and boundary instabilities associated with standard smooth particle hydrodynamics formulations and the diffusive grid- to-particle mapping schemes characteristic of particle-in-cell methods. In the work described herein, the method of Fahrenthold and Koo has been extended, by coupling the aforementioned hydrodynamic particle model to a hexahedral finite element based description of the continuum dynamics. The resulting continuum model retains all of the features (including general contact-impact effects) of Hamiltonian particle hydrodynamics, while in addition accounting for tensile strength, plasticity, and damage effects important in the simulation of hypervelocity impact on orbital debris shielding. A three dimensional, vectorized, and autotasked implementation of the extended particle method described here has been coded for application to orbital debris shielding design. Source code for the pre-processor (PREP), analysis code (EXOS), post-processor (POST), and rezoner (ZONE), have been delivered separately, along with a User's Guide describing installation and application of the software.
Madland, D.G.; Arthur, E.D.; Estes, G.P.; Stewart, J.E.; Bozoian, M.; Perry, R.T.; Parish, T.A.; Brown, T.H.; England, T.R.; Wilson, W.B.; Charlton, W.S.
1999-09-01
SOURCES 4A is a computer code that determines neutron production rates and spectra from ({alpha},n) reactions, spontaneous fission, and delayed neutron emission due to the decay of radionuclides. The code is capable of calculating ({alpha},n) source rates and spectra in four types of problems: homogeneous media (i.e., a mixture of {alpha}-emitting source material and low-Z target material), two-region interface problems (i.e., a slab of {alpha}-emitting source material in contact with a slab of low-Z target material), three-region interface problems (i.e., a thin slab of low-Z target material sandwiched between {alpha}-emitting source material and low-Z target material), and ({alpha},n) reactions induced by a monoenergetic beam of {alpha}-particles incident on a slab of target material. Spontaneous fission spectra are calculated with evaluated half-life, spontaneous fission branching, and Watt spectrum parameters for 43 actinides. The ({alpha},n) spectra are calculated using an assumed isotropic angular distribution in the center-of-mass system with a library of 89 nuclide decay {alpha}-particle spectra, 24 sets of measured and/or evaluated ({alpha},n) cross sections and product nuclide level branching fractions, and functional {alpha}-particle stopping cross sections for Z < 106. The delayed neutron spectra are taken from an evaluated library of 105 precursors. The code outputs the magnitude and spectra of the resultant neutron source. It also provides an analysis of the contributions to that source by each nuclide in the problem.
Multigroup Time-Independent Neutron Transport Code System for Plane or Spherical Geometry.
Energy Science and Technology Software Center (ESTSC)
1986-12-01
Version 00 PALLAS-PL/SP solves multigroup time-independent one-dimensional neutron transport problems in plane or spherical geometry. The problems solved are subject to a variety of boundary conditions or a distributed source. General anisotropic scattering problems are treated for solving deep-penetration problems in which angle-dependent neutron spectra are calculated in detail.
Dawahra, S; Khattab, K; Saba, G
2015-10-01
A comparative study for fuel conversion from the HEU to LEU in the Miniature Neutron Source Reactor (MNSR) has been performed in this paper using the MCNP4C code. The neutron energy and lethargy flux spectra in the first inner and outer irradiation sites of the MNSR reactor for the existing HEU fuel (UAl4-Al, 90% enriched) and the potential LEU fuels (U3Si2-Al, U3Si-Al, U9Mo-Al, 19.75% enriched and UO2, 12.6% enriched) were investigated using the MCNP4C code. The neutron energy flux spectra for each group was calculated by dividing the neutron flux by the width of each energy group. The neutron flux spectra per unit lethargy was calculated by multiplying the neutron energy flux spectra for each energy group by the average energy of each group. The thermal neutron flux was calculated by summing the neutron fluxes from 0.0 to 0.625 eV, the fast neutron flux was calculated by summing the neutron fluxes from 0.5 MeV to 10 MeV for the existing HEU and potential LEU fuels. Good agreements have been noticed between the flux spectra for the potential LEU fuels and the existing HEU fuels with maximum relative differences less than 10% and 8% in the inner and outer irradiation sites. PMID:26142805
Wicaksono, D.; Zerkak, O.; Nikitin, K.; Ferroukhi, H.; Chawla, R.
2013-07-01
This paper reports refinement studies on the temporal coupling scheme and time-stepping management of TRACE/S3K, a dynamically coupled code version of the thermal-hydraulics system code TRACE and the 3D core simulator Simulate-3K. The studies were carried out for two test cases, namely a PWR rod ejection accident and the Peach Bottom 2 Turbine Trip Test 2. The solution of the coupled calculation, especially the power peak, proves to be very sensitive to the time-step size with the currently employed conventional operator-splitting. Furthermore, a very small time-step size is necessary to achieve decent accuracy. This degrades the trade-off between accuracy and performance. A simple and computationally cheap implementation of time-projection of power has been shown to be able to improve the convergence of the coupled calculation. This scheme is able to achieve a prescribed accuracy with a larger time-step size. (authors)
Code System to Calculate Neutron and Gamma-Ray Skyshine Doses Using the Integral Line-Beam Method.
Energy Science and Technology Software Center (ESTSC)
2000-11-16
Version 03 This package includes the SKYNEUT 1.1, SKYDOSE 2.3, MCSKY 2.3 and SKYCONES 1.1 codes plus the DLC-188/SKYDATA library to form a comprehensive system for calculating skyshine doses. See the author's web site for related information: http://athena.mne.ksu.edu/~jks/ SKYNEUT evaluates the neutron and neutron-induced secondary gamma-ray skyshine doses from an isotropic, point, neutron source collimated by three simple geometries: an open silo, a vertical black (perfectly absorbing) wall, and a rectangular building. The source maymore » emit monoenergetic neutrons or neutrons with an arbitrary multigroup spectrum of energies. SKYDOSE evaluates the gamma-ray skyshine dose from an isotropic, monoenergetic, point gamma-photon source collimated by three simple geometries: (1) a source in a silo, (2) a source behind an infinitely long, vertical, black wall, and (3) a source in a rectangular building. In all three geometries an optional overhead slab shield may be specified. MCSKY evaluates the gamma-ray skyshine dose from an isotropic, monoenergetic, point gamma-photon source collimated into either a vertical cone (i.e., silo geometry) or into a vertically oriented structure with an N-sided polygon cross section. An overhead laminate shield composed of two different materials is assumed, although shield thicknesses of zero may be specified to model an unshielded SKYSHINE source. SKYCONES evaluates the skyshine doses produced by a point neutron or gamma-photon source emitting, into the atmosphere, radiation that is collimated into an upward conical annulus between two arbitrary polar angles. The source is assumed to be axially (azimuthally) symmetric about a vertical axis through the source and can have an arbitrary polyenergetic spectrum. Nested contiguous annular cones can thus be used to represent the energy and polar-angle dependence of a skyshine source emitting radiation into the atmosphere.« less
Ui, Atsushi; Miyaji, Takamasa
2004-10-15
The best-estimate coupled three-dimensional (3-D) core and thermal-hydraulic code system TRAC-BF1/COS3D has been developed. COS3D, based on a modified one-group neutronic model, is a 3-D core simulator used for licensing analyses and core management of commercial boiling water reactor (BWR) plants in Japan. TRAC-BF1 is a plant simulator based on a two-fluid model. TRAC-BF1/COS3D is a coupled system of both codes, which are connected using a parallel computing tool. This code system was applied to the OECD/NRC BWR Turbine Trip Benchmark. Since the two-group cross-section tables are provided by the benchmark team, COS3D was modified to apply to this specification. Three best-estimate scenarios and four hypothetical scenarios were calculated using this code system. In the best-estimate scenario, the predicted core power with TRAC-BF1/COS3D is slightly underestimated compared with the measured data. The reason seems to be a slight difference in the core boundary conditions, that is, pressure changes and the core inlet flow distribution, because the peak in this analysis is sensitive to them. However, the results of this benchmark analysis show that TRAC-BF1/COS3D gives good precision for the prediction of the actual BWR transient behavior on the whole. Furthermore, the results with the modified one-group model and the two-group model were compared to verify the application of the modified one-group model to this benchmark. This comparison shows that the results of the modified one-group model are appropriate and sufficiently precise.
Effects of tensor couplings of ω and ρ mesons on 1S0 nucleon superfluidity in neutron star matter
NASA Astrophysics Data System (ADS)
Huang, Xiu-Lin; Xu, Yan; Liu, Cheng-Zhi; Liu, Guang-Zhou
2016-03-01
The 1S 0 nucleon superfluidity in neutron star matter was investigated in the framework of relativistic σ-ω-π-ρ model with the tensor couplings of ω and ρ mesons using the relativistic Hartree-Fock (RHF) approximation. It was found that the tensor couplings of ω and ρ mesons lead to a clear growth of the 1S 0 neutron pairing gap in the density range where there exists 1S 0 neutron superfluidity. The 1S 0 pairing gap of proton with the tensor couplings of ω and ρ mesons in the density range of ρB = 0.0-0.079fm-3 is lower and then in the density range of ρB = 0.079-0.383fm-3 higher than the corresponding value without the tensor couplings of ω and ρ mesons. Our results provide a basic to understand the influence of the tensor couplings of ω and ρ mesons on the cooling properties of neutron star.
Neutron scattering from coupled phonon-impurity modes in KCl1-c(KCN)c
NASA Astrophysics Data System (ADS)
Nicklow, R. M.; Crummett, W. P.; Mostoller, M.; Wood, R. F.
1980-09-01
The hybridization of host-lattice phonons with the internal-energy states of CN- impurities in KCl has been studied by inelastic neutron scattering as a function of temperature between 10 and 100 K for samples with impurity concentrations in the range c=0.4 to 6 at.%. A temperature- and concentration-dependent coupling between phonons with Eg symmetry and the Eg transitions of CN- ions is observed near a frequency of 0.5 THz, a value which is consistent with the energy-level spacings for CN- in KCl as deduced by Beyeler. However, an expected coupling of phonons and CN- transitions with T2g symmetry near the same frequency was not detected. A simple two-level model for the CN- impurity provides a rather good description of the data for the Eg coupled modes for c<2 at.%, but it deviates significantly for larger concentrations. Quasielastic scattering, which has a strong dependence on impurity concentration, temperature, and phonon wave vector, is also observed.
Simulations of the Dynamics of the Coupled Energetic and Relativistic Electrons Using VERB Code
NASA Astrophysics Data System (ADS)
Shprits, Y.; Kellerman, A. C.; Drozdov, A.
2015-12-01
Modeling and understanding of ring current and radiation belt coupled system has been a grand challenge since the beginning of the space age. In this study we show long term simulations with a 3D VERB code of modeling the radiation belts with boundary conditions derived from observations around geosynchronous orbit. We also present 4D VERB simulations that include convective transport, radial diffusion, pitch angle scattering and local acceleration. VERB simulations show that the lower energy inward transport is dominated by the convection and higher energy transport is dominated by the diffusive radial transport. We also show that at energies of 100s of keV a number of processes work simultaneously including convective transport, radial diffusion, local acceleration, loss to the loss cone and loss to the magnetopause. The results of the simulaiton of March 2013 storm are compared with Van Allen Probes observations.
Simulation of high-energy radiation belt electron fluxes using NARMAX-VERB coupled codes
Pakhotin, I P; Drozdov, A Y; Shprits, Y Y; Boynton, R J; Subbotin, D A; Balikhin, M A
2014-01-01
This study presents a fusion of data-driven and physics-driven methodologies of energetic electron flux forecasting in the outer radiation belt. Data-driven NARMAX (Nonlinear AutoRegressive Moving Averages with eXogenous inputs) model predictions for geosynchronous orbit fluxes have been used as an outer boundary condition to drive the physics-based Versatile Electron Radiation Belt (VERB) code, to simulate energetic electron fluxes in the outer radiation belt environment. The coupled system has been tested for three extended time periods totalling several weeks of observations. The time periods involved periods of quiet, moderate, and strong geomagnetic activity and captured a range of dynamics typical of the radiation belts. The model has successfully simulated energetic electron fluxes for various magnetospheric conditions. Physical mechanisms that may be responsible for the discrepancies between the model results and observations are discussed. PMID:26167432
Warthog: A MOOSE-Based Application for the Direct Code Coupling of BISON and PROTEUS
McCaskey, Alexander J.; Slattery, Stuart; Billings, Jay Jay
2015-09-01
The Nuclear Energy Advanced Modeling and Simulation (NEAMS) program from the Department of Energy s Office of Nuclear Energy provides a robust toolkit for the modeling and simulation of current and future advanced nuclear reactor designs. This toolkit provides these technologies organized across product lines: two divisions targeted at fuels and end-to-end reactor modeling, and a third for integration, coupling, and high-level workflow management. The Fuels Product Line and the Reactor Product line provide advanced computational technologies that serve each respective field well, however, their current lack of integration presents a major impediment to future improvements of simulation solution fidelity. There is a desire for the capability to mix and match tools across Product Lines in an effort to utilize the best from both to improve NEAMS modeling and simulation technologies. This report will detail a new effort to provide this Product Line interoperability through the development of a new application called Warthog. This application couples the BISON Fuel Performance application from the Fuels Product Line and the PROTEUS Core Neutronics application from the Reactors Product Line in an effort to utilize the best from all parts of the NEAMS toolkit and improve overall solution fidelity of nuclear fuel simulations. To acheive this, Warthog leverages as much prior work from the NEAMS program as possible, and in doing so, enables interoperability between the disparate MOOSE and SHARP frameworks, and the libMesh and MOAB mesh data formats. The remainder of this report will describe this work in full. We will begin with a detailed look at the individual NEAMS framework technologies used and developed in the various Product Lines, and the current status of their interoperability. We will then introduce the Warthog application: its overall architecture and the ways it leverages the best existing tools from accross the NEAMS toolkit to enable BISON-PROTEUS integration
Coupled Neutronics Thermal-Hydraulic Solution of a Full-Core PWR Using VERA-CS
Clarno, Kevin T; Palmtag, Scott; Davidson, Gregory G; Salko, Robert K; Evans, Thomas M; Turner, John A; Belcourt, Kenneth; Hooper, Russell; Schmidt, Rodney
2014-01-01
The Consortium for Advanced Simulation of Light Water Reactors (CASL) is developing a core simulator called VERA-CS to model operating PWR reactors with high resolution. This paper describes how the development of VERA-CS is being driven by a set of progression benchmark problems that specify the delivery of useful capability in discrete steps. As part of this development, this paper will describe the current capability of VERA-CS to perform a multiphysics simulation of an operating PWR at Hot Full Power (HFP) conditions using a set of existing computer codes coupled together in a novel method. Results for several single-assembly cases are shown that demonstrate coupling for different boron concentrations and power levels. Finally, high-resolution results are shown for a full-core PWR reactor modeled in quarter-symmetry.
Glushkov, Y.S.; Ponomarev-Stepnoi, N.N.; Kompanietz, G.V.; Gomin, Y.A.; Maiorov, L.V.; Lobynstev, V.A.; Polyakov, D.N.; Sapir, J.; Streetman, J.R.
1993-11-01
Topaz-II is a heterogeneous, epithermal reactor, fueled with highly enriched uranium-dioxide, cooled with NaK, and moderated with zirconium-hydride. The reactor core contains 37 single-cell thermionic fuel elements, and is surrounded by a radial beryllium reflector that contains 12 rotatable control drums with poison segments. For the physics analysis of TOPAZ II it is necessary to use the Monte Carlo method. The United States (US) and Russia used two different Monte Carlo codes, namely MCNP and MCU-2, respectively. The work described in this paper was aimed at comparing the codes and neutronic data used in the US and Russia for verification of Topaz-II nuclear safety. For this purpose, the US and Russia developed a joint benchmark model of the Topaz-II reactor. The American and Russian teams performed independent computations for a series of variants representing potential water immersion accidents. Comparison of the MCNP and MCU-2 codes showed somewhat different results both for the absolute values of k{sub eff} and for reactivity effects. Future calculations will be performed to obtain a detailed understanding of the reasons for such discrepancies. For these analyses it will be necessary for the US and Russian teams to exchange neutronic data on Topaz-II physics calculations.
NASA Technical Reports Server (NTRS)
Ghorai, S. K.
1983-01-01
The purpose of this project was to use a one-dimensional discrete coordinates transport code called ANISN in order to determine the energy-angle-spatial distribution of neutrons in a 6-feet cube rock box which houses a D-T neutron generator at its center. The project was two-fold. The first phase of the project involved adaptation of the ANISN code written for an IBM 360/75/91 computer to the UNIVAC system at JSC. The second phase of the project was to use the code with proper geometry, source function and rock material composition in order to determine the neutron flux distribution around the rock box when a 14.1 MeV neutron generator placed at its center is activated.
NASA Astrophysics Data System (ADS)
Ghorai, S. K.
1983-09-01
The purpose of this project was to use a one-dimensional discrete coordinates transport code called ANISN in order to determine the energy-angle-spatial distribution of neutrons in a 6-feet cube rock box which houses a D-T neutron generator at its center. The project was two-fold. The first phase of the project involved adaptation of the ANISN code written for an IBM 360/75/91 computer to the UNIVAC system at JSC. The second phase of the project was to use the code with proper geometry, source function and rock material composition in order to determine the neutron flux distribution around the rock box when a 14.1 MeV neutron generator placed at its center is activated.
Role of nuclear couplings in the inelastic excitation of weakly-bound neutron-rich nuclei
Dasso, C.H.; Lenzi, S.M.; Vitturi, A.
1996-12-31
Much effort is presently devoted to the study of nuclear systems far from the stability line. Particular emphasis has been placed in light systems such as {sup 11}Li, {sup 8}B and others, where the very small binding energy of the last particles causes their density distribution to extend considerably outside of the remaining nuclear core. Some of the properties associated with this feature are expected to characterize also heavier systems in the vicinity of the proton or neutron drip lines. It is by now well established that low-lying concentrations of multipole strength arise from pure configurations in which a peculiar matching between the wavelength of the continuum wavefunction of the particles and the range of the weakly-bound hole states occurs. To this end the authors consider the break-up of a weakly-bound system in a heavy-ion collision and focus attention in the inelastic excitation of the low-lying part of the continuum. They make use of the fact that previous investigations have shown that the multipole response in this region is not of a collective nature and describe their excited states as pure particle-hole configurations. Since the relevant parameter determining the strength distributions is the binding energy of the last bound orbital they find it most convenient to use single-particle wavefunctions generated by a sperical square-well potential with characteristic nuclear dimensions and whose depth has been adjusted to give rise to a situation in which the last occupied neutron orbital is loosely-bound. Spin-orbit couplings are, for the present purpose, ignored. The results of this investigation clearly indicate that nuclear couplings have the predominant role in causing projectile dissociation in many circumstances, even at bombarding energies remarkably below the Coulomb barrier.