Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer Codes”

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

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

At least 1,153 records · Page 64

Substantiation of parameters of the geometric model of the research reactor core for the calculation using the Monte Carlo method

The choice of the spatial nodalization for the calculation of the power density and burnup distribution in a research reactor core with fuel assemblies of the IRT-3M and VVR-KN type using the program based on the Monte Carlo code is described. The influence of the spatial nodalization on the results of calculating basic neutronic characteristics and calculation time is investigated.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

On the image of AGS 3 He 2+ $\vec{n}$ 0 in the blue

This note addresses the transport of Helion spin $\vec{n}$ 0 vector, from its periodic orientation in the AGS to RHIC Blue ring injection kicker, via the AGS extraction system and the AtR line. The goal is to investigate optimal injection energy into RHIC, in the matter of Helion spin matching, in the hypothesis of equal warm and cold snake strengths in the AGS. The study uses recently computed OPERA 3-D field maps of the AGS cold snake, including possibility of independent solenoid and helix settings (as discussed in Tech. Note C-A/AP/485), together with the machinery of the AGS and AtR models developed in the stepwise ray-tracing code Zgoubi. Computing tools and methods employed are discussed as well, in order to facilitate possible further checks or investigations. They are however similar to those used in an earlier study regarding the image in RHIC Blue and Yellow of AGS $\vec{n}$ 0 via the AtR in the case of proton beam (Tech. Note C-A/AP/502), which can be referred to for additional details.

43 PARTICLE ACCELERATORS↗

Validation of the U.S. NRC coupled code system TRITON/TRACE/PARCS with the special power excursion reactor test III (SPERT III)

The Special Power Excursion Reactor Test III (SPERT III) was a series of reactivity insertion experiments conducted in the 1950's. This paper describes the validation of the U.S. NRC Coupled Code system TRITON/PARCS/TRACE to simulate reactivity insertion accidents (RIA) by using several of the SPERT III tests. The work here used the SPERT III E-core configuration tests in which the RIA was initiated by ejecting a control rod. The resulting super-prompt reactivity excursion and negative reactivity feedback produced the familiar bell shaped power increase and decrease. The energy deposition during such a power peak has important safety consequences and provides validation basis for core coupled multi-physics codes. The transients of five separate tests are used to benchmark the PARCS/TRACE coupled code. The models were thoroughly validated using the original experiment documentation. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Effect of IAEA patch for TRANSX2.15 - 25264

IAEA released two unofficial patches for TRANSX2.15 for the MATXS file of FENDL-2.0 in 1998, but the patches have not been known well. The first patch is for all MATXS files, while the second one is effective only for MATXS files with mixed data types, which are very rare. This paper demonstrates with simple calculations that the first patch is essential because it solves a problem of extremely large neutron fluxes caused by the original TRANSX2.15 code. We hope that the patch will be included in the official TRANSX2.15 release. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

AMPX: A Modern cross section processing system for generating nuclear data libraries - 14619

Nuclear data libraries are essential for nuclear analysis code systems like SCALE. AMPX has been continuously developed since the early 1970's and is used to generate these data libraries for SCALE. Recently, Oak Ridge National Laboratory (ORNL) initiated a major modernization effort to position AMPX for continued development and deployment on modern computing platforms. An important step toward modernization has been to merge AMPX under the SCALE continuous integration (CI) development framework. As a result, AMPX is now developed under the SCALE software quality assurance plan, providing increased confidence in the quality of data libraries developed and deployed by AMPX. This merge allowed us to develop novel methods to generate self-shielding factors that combine SCALE's resonance self-shielding capabilities with the AMPX library-generation capabilities. AMPX processes Evaluated Nuclear Data Format (ENDF/B) files. The ENDF community is working to develop a more modern format that will replace the existing ENDF/B format. Due to strict space and memory requirements when AMPX was first developed, the original AMPX processing and reading of ENDF data were tightly integrated. On newer computer systems, this is no longer necessary. Now there is a separation between reading and processing which allows for simultaneous support of current and new ENDF formats. As a result, AMPX is positioned to adopt any new ENDF format when available. For end users, the modernization effort has resulted in improved processing capabilities for generating nuclear data libraries. The objective of the paper is to document the latest AMPX cross-section processing capabilities and provide the current status of the AMPX modernization effort. (authors)

42 ENGINEERING↗

MARS-KS code validation activity through the atlas domestic standard problem

The 2 nd Domestic Standard Problem (DSP-02) exercise using the ATLAS integral effect test data was executed to transfer the integral effect test data to domestic nuclear industries and to contribute to improving the safety analysis methodology for PWRs. A small break loss of coolant accident of a 6-inch break at the cold leg was determined as a target scenario by considering its technical importance and by incorporating interests from participants. Ten calculation results using MARS-KS code were collected, major prediction results were described qualitatively and code prediction accuracy was assessed quantitatively using the FFTBM. In addition, special code assessment activities were carried out to find out the area where the model improvement is required in the MARS-KS code. The lessons from this DSP-02 and recommendations to code developers are described in this paper. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Proton-induced fission of actinides at energies 26.5 and 62.9 MeV--Theoretical interpretation

Fission properties of proton-induced fission on {sup 232}Th, {sup 237}Np, {sup 238}U, {sup 239}Pu and {sup 241}Am targets, measured at the Louvain-la-Neuve cyclotron facility at proton energies of 26.5 and 62.9 MeV, are compared with the predictions of the state-of-the-art nuclear reaction code TALYS. The sensitivity of the calculations to the input parameters of the code and possible improvements are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Investigation on plasma generated during hypervelocity impact at different impact velocities and angles

A 3D Smoothed Particle Hydrodynamics code was developed to investigate plasma generation by considering a chemical reaction process in hypervelocity impacts of an aluminum projectile on an aluminum target. The chemical reaction process was described by the reaction rate based on the Arrhenius equation and used to calculate the plasma generation during the impact simulation. The predicted result was verified by empirical formulas and a new empirical formula was proposed based on the comparisons and analyses. The influence of the impact angle was discussed for different impact velocities. Then, the application of both the new and original empirical formulas for protection design from plasma generated by hypervelocity impact was discussed, which demonstrated that the code and model were useful in the prediction of hypervelocity impacts on spacecraft.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Single pin BWR benchmark problem for coupled Monte Carlo - Thermal hydraulics analysis

As part of the European NURISP research project, a single pin BWR benchmark problem was defined. The aim of this initiative is to test the coupling strategies between Monte Carlo and subchannel codes developed by different project participants. In this paper the results obtained by the Delft Univ. of Technology and Karlsruhe Inst. of Technology will be presented. The benchmark problem was simulated with the following coupled codes: TRIPOLI-SUBCHANFLOW, MCNP-FLICA, MCNP-SUBCHANFLOW, and KENO-SUBCHANFLOW. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

CATHARE thermal-hydraulic system code for HLM preliminary validation in natural convection tests

The innovative nuclear systems cooled by Heavy Liquid Metal (HLM) are the subject of an ongoing interest both in Europe and outside, evidenced by a number of projects in progress. In the frame of the European Framework Programmes have been evidenced the need to adopt a thermalhydraulic system code capable to treat lead and Lead-Bismuth Eutectic (LBE) systems, with a particular interest in developing a 'European' code. Considering this scenario, within a specific collaboration between ENEA and CEA, the CATHARE French system code has been modified extending its capabilities to simulate HLM systems. In the present paper, the state of the validation process of CATHARE-HLM is discussed. The activity aims to assess the capabilities and limitations of the code to simulate the behavior of integral facilities, in particular in natural circulation conditions. The experimental data come from NACIE LBE-cooled facility sited at the ENEA Brasimone laboratories. The results obtained show a good capability in reproducing the systems behavior, despite some uncertainties on the experimental measurements. Future improvements on the code are going to be planned within the collaboration ENEA/CEA. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Validation of continuous energy sensitivity coefficients calculations in the MORET 5 code - 14269

A continuous-energy sensitivity coefficients calculation capability has been developed in the MORET 5 Monte Carlo code recently. As discussed in this paper, it relies on the differential operator method, an adjoint weighting model (to avoid calculating the fission source derivatives) and a constrained sensitivity coefficient model for sensitivity profiles to distributions such as, for instance, the fission spectrum. For now, the coefficients are provided per medium, per isotope, per reaction on an energy binning defined by the user. The sensitivity coefficients are currently available for the following quantities: total, elastic scattering, total inelastic scattering, fission, total capture and radiative capture cross sections, average neutron produced by fission and fission spectrum. This paper deals with the validation of this new feature in the MORET code on several benchmarks and applications with various energy spectra and geometrical configurations. This validation is performed through the UACSA expert group and the Jezebel ICSBEP benchmarks, but also thanks to calculations performed in the frame of experimental program design in the STACY facility. Comparisons are made using the MCNP6.1 code and the SCALE6.1 package as references. Preliminary results show a pretty good agreement between the different software. (authors)

42 ENGINEERING↗

SU-E-T-103: Development and Implementation of Web Based Quality Control Software

Purpose: Historically many radiation medicine programs have maintained their Quality Control (QC) test results in paper records or Microsoft Excel worksheets. Both these approaches represent significant logistical challenges, and are not predisposed to data review and approval. It has been our group's aim to develop and implement web based software designed not just to record and store QC data in a centralized database, but to provide scheduling and data review tools to help manage a radiation therapy clinics Equipment Quality control program. Methods: The software was written in the Python programming language using the Django web framework. In order to promote collaboration and validation from other centres the code was made open source and is freely available to the public via an online source code repository. The code was written to provide a common user interface for data entry, formalize the review and approval process, and offer automated data trending and process control analysis of test results. Results: As of February 2014, our installation of QAtrack+ has 180 tests defined in its database and has collected ∼22 000 test results, all of which have been reviewed and approved by a physicist via QATrack+'s review tools. These results include records for quality control of Elekta accelerators, CT simulators, our brachytherapy programme, TomoTherapy and Cyberknife units. Currently at least 5 other centres are known to be running QAtrack+ clinically, forming the start of an international user community. Conclusion: QAtrack+ has proven to be an effective tool for collecting radiation therapy QC data, allowing for rapid review and trending of data for a wide variety of treatment units. As free and open source software, all source code, documentation and a bug tracker are available to the public at https://bitbucket.org/tohccmedphys/qatrackplus/.

07 ISOTOPE AND RADIATION SOURCES↗

Qualification of Numerical and Semi-Analytical Modeling Codes for Thermal Design of HLW Repository - 18455

Disposal in a geological repository is a viable option for safely isolating high-level radioactive waste. Thermal processes induced by decay heat from the radioactive waste may impact the integrity of the engineered barrier system, and may lead to fluid flow and transport; and the associated migration of radionuclides. This study looks at an analysis of thermal conduction for the emplacement of high-level waste in a repository located in bedded salt. The simulations used analytical and numerical methods to study conductive heat transport in a repository. This study concentrated on heat flow in a central waste package with contributions from adjacent waste packages and drifts. Of particular importance to this work is comparison of predictions of peak temperature at the drift wall and at the waste package surface using different solution methods. Simulation results of analytical and numerical modeling methods are compared. The results show that predictions of the codes were comparable under identical initial and boundary conditions. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Verification of a Multiphysics Code with Method of Manufactured Solutions

To push the modeling and simulation to a new level in simulating realistic problems, a common practice is to couple different single-physics codes together to simulate problems involving multiple physics. Although the method of manufactured solutions (MMS) has been shown to be a powerful tool in code verification, there has been little work on applying MMS to a multi-physics code system. Additionally, due to the scarcity of multi-physics validation data, the significance of code verification becomes more prominent as an important part in software quality assurance. With this, the application of MMS to a coupled code system, e.g. VERA-CS, has become a crucial step, which helps identify bugs affecting the order of accuracy, thus increasing the credibility and reliability of the code. In this work, we investigate the application of MMS to verify a multi-physics code system and establish a framework for this application, and we demonstrate this with a neutronics code coupled with a thermal hydraulics (TH) code. The convection module is not yet added to this work, so by TH, only conduction is considered.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

SU-E-T-557: Monte Carlo Modeling of Philips RT-250 Orthovoltage Unit for Beam Spectrum Modulation

Purpose: To develop and validate a Monte Carlo (MC) model of a Phillips RT-250 orthovoltage unit to test various beam spectrum modulation strategies for in vitro/vivo studies. A model of this type would enable the production of unconventional beams from a typical orthovoltage unit for novel therapeutic applications such as gold nanoparticle-aided radiotherapy. Methods: The MCNP5 code system was used to create a MC model of the head of RT-250 and a 30 × 30 × 30 cm{sup 3} water phantom. For the x-ray machine head, the current model includes the vacuum region, beryllium window, collimators, inherent filters and exterior steel housing. For increased computational efficiency, the primary x-ray spectrum from the target was calculated from a well-validated analytical software package. Calculated percentage-depth-dose (PDD) values and photon spectra were validated against experimental data from film and Compton-scatter spectrum measurements. Results: The model was validated for three common settings of the machine namely, 250 kVp (0.25 mm Cu), 125 kVp (2 mm Al), and 75 kVp (2 mm Al). The MC results for the PDD curves were compared with film measurements and showed good agreement for all depths with a maximum difference of 4 % around dmax and under 2.5 % for all other depths. The primary photon spectra were also measured and compared with the MC results showing reasonable agreement between the two, validating the input spectra and the final spectra as predicted by the current MC model. Conclusion: The current MC model accurately predicted the dosimetric and spectral characteristics of each beam from the RT-250 orthovoltage unit, demonstrating its applicability and reliability for beam spectrum modulation tasks. It accomplished this without the need to model the bremsstrahlung xray production from the target, while significantly improved on computational efficiency by at least two orders of magnitude. Supported by DOD/PCRP grant W81XWH-12-1-0198.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

On the image in RHIC of AGS spin vector via the AtR

This note re-visits the transport of the spin nn0 vector, from its periodic orientation in the AGS, to the downstream end of the Blue and Yellow RHIC ring injection kickers via the AGS extraction system and the AtR line. The goal is to verify the optimal injection energy into RHIC, in matter of spin matching. This is done using newly computed OPERA 3-D field maps of the AGS cold snake, including a new setting of the cold snake solenoid as discussed in a companion Tech. Note (C-A/AP/485, July 2013), together with the machinery of the AGS and AtR models developed in the stepwise ray-tracing code Zgoubi. Computing tools and methods employed in these investigations are discussed as well, in order to facilitate possible further checks or investigations.

43 PARTICLE ACCELERATORS↗

PPT-GPU: Scalable GPU Performance Modeling

Performance modeling is a challenging problem due to the complexities of hardware architectures. In this paper, we present PPT-GPU, a scalable and accurate simulation framework that enables GPU code developers and architects to predict the performance of applications in a fast, and accurate manner on different GPU architectures. PPT-GPU is part of the open source project, Performance Prediction Toolkit (PPT) developed at the Los Alamos National Laboratory. We extend the old GPU model in PPT that predict the runtimes of computational physics codes to offer better prediction accuracy, for which, we add models for different memory hierarchies found in GPUs and latencies for different instructions. To further show the utility of PPT-GPU, we compare our model against real GPU device(s) and the widely used cycle-accurate simulator, GPGPU-Sim using different workloads from RODINIA and Parboil benchmarks. The results indicate that the predicted performance of PPT-GPU is within a 10% error compared to the real device(s). In addition, PPT-GPU is highly scalable, where it is up to 450x faster than GPGPU-Sim with more accurate results.

97 MATHEMATICS AND COMPUTING↗

ABSORPTION EFFICIENCIES OF FORSTERITE. I. DISCRETE DIPOLE APPROXIMATION EXPLORATIONS IN GRAIN SHAPE AND SIZE

We compute the absorption efficiency (Q{sub abs}) of forsterite using the discrete dipole approximation in order to identify and describe what characteristics of crystal grain shape and size are important to the shape, peak location, and relative strength of spectral features in the 8-40 {mu}m wavelength range. Using the DDSCAT code, we compute Q{sub abs} for non-spherical polyhedral grain shapes with a{sub eff} = 0.1 {mu}m. The shape characteristics identified are (1) elongation/reduction along one of three crystallographic axes; (2) asymmetry, such that all three crystallographic axes are of different lengths; and (3) the presence of crystalline faces that are not parallel to a specific crystallographic axis, e.g., non-rectangular prisms and (di)pyramids. Elongation/reduction dominates the locations and shapes of spectral features near 10, 11, 16, 23.5, 27, and 33.5 {mu}m, while asymmetry and tips are secondary shape effects. Increasing grain sizes (0.1-1.0 {mu}m) shifts the 10 and 11 {mu}m features systematically toward longer wavelengths and relative to the 11 {mu}m feature increases the strengths and slightly broadens the longer wavelength features. Seven spectral shape classes are established for crystallographic a-, b-, and c-axes and include columnar and platelet shapes plus non-elongated or equant grain shapes. The spectral shape classes and the effects of grain size have practical application in identifying or excluding columnar, platelet, or equant forsterite grain shapes in astrophysical environs. Identification of the shape characteristics of forsterite from 8 to 40 {mu}m spectra provides a potential means to probe the temperatures at which forsterite formed.

79 ASTRONOMY AND ASTROPHYSICS↗