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,243 records · Page 69

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↗

Development of INSPCT-S for inspection of spent fuel pool

In this paper, we discuss an accurate and fast software tool (INSPCT-S, Inspection of Nuclear Spent fuel-Pool Calculation Tool, version Spreadsheet) developed for calculation of the response of fission chambers placed in a spent fuel pool, such as Atucha-I. INSPCT-S is developed for identification of suspicious regions of the pool that may have missing or substitute assemblies. INSPCT-S uses a hybrid algorithm based on the adjoint function methodology. The neutron source is comprised of spontaneous fission, ({alpha}, n) interactions, and subcritical multiplication. The former is evaluated using the ORIGEN-ARP code, and the latter is obtained with the fission matrix (FM) formulation. The FM coefficients are determined using the MCNP Monte Carlo code, and the importance function is determined using the PENTRAN 3-D parallel Sn code. Three databases for the neutron source, FM elements, and adjoint flux are prepared as functions of different parameters including burnup, cooling time, enrichment, and pool lattice size. INSPCT-S uses the aforementioned databases and systems of equations to calculate detector responses, which are subsequently compared with normalized experimental data. If this comparison is not satisfied, INSPCT-S utilizes color coding to identify the suspicious regions of a spent fuel pool. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

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↗

Human Liver Epithelium Response to HCoV-229E Infection Epigenomics (ACS-DP4)

The purpose of this experiment was to evaluate how wild-type Human coronavirus strain 229E (HCoV-229E) infection alters chromatin accessibility in infected cells only. Sample data was obtained for mock and infected (standard and UV-inactivated) immortalized human liver cells (HuH-7) and collected 24 hrs. post infection. Samples were processed using assay for transposase-accessible chromatin using high-throughput sequencing (ATAC-Seq) and generated bar coded library samples were evaluated for RNA sequencing (RNA-Seq) expression analysis. Processed ATAC-Seq datasets are openly accessible from the download button and contain secondary processed RNA-Seq results files and supporting metadata materials. Data download includes a sample naming key, infection titer metadata, normalized counts, and relevant computational source code information supporting data transparency and reuse.

59 BASIC BIOLOGICAL SCIENCES↗

Analysis of Spectral Lines in Large Databases of Synthetic Spectra for Massive Stars

In this paper, we describe a program that identifies in the optical spectrum the main parameters of a spectral line, namely the initial and final wavelengths, and the line depth. Moreover, using numerical calculations, it identifies and removes adjacent lines. Next, the program calculates the equivalent width and the FWHM. The software was tested in a sample of 300 lines in two databases of synthetic spectra generated by the CMFGEN and PoWR codes, and 300 lines in observed spectra from the IACOB database, showing a Gaussian distribution of relative errors, from which it is inferred that 80% of the measured lines have errors less than 17% and only 5% of the lines have errors greater than 26%. The program was also run on the entire database of 45,000 CMFGEN and 202 POWR synthetic spectra, generating a library of H i, He i, and He ii lines necessary to feed the FITspec code for the derivation of stellar parameters: effective temperature, surface gravity, and luminosity.

47 OTHER INSTRUMENTATION↗

Impedance Measurement of an Antenna with Hydrogen Plasma Driven by 13.56 MHz-rf for J-PARC H{sup -} ion source

In order to satisfy the Japan Proton Accelerator Research Complex (J-PARC) second stage requirements of an H{sup -} ion beam current of 50 mA and a life-time of 2000 hours, the development of a 13.56 MHz-rf-driven H{sup -} ion source was started by using an antenna developed at the Spallation Neutron Source (SNS). As the first step, the impedance of an antenna with hydrogen plasma, which is one of the most important parameters, was measured. For an rf-power of 0 kW, 1.25 kW or 7.45 kW, it was measured as j84.19+5.107 {Omega}, j80.22+11.82 {Omega} or j58.24+21.16 {Omega}, respectively. Its real-part had a rather large value and rather large rf-power dependence unexpectedly. By using the measured impedance, the circuit to produce rf-plasma in the plasma chamber on a -50 kV potential with a 13.56 MHz-rf-source on the ground potential was designed by the circuit simulation code LTSpice IV.

43 PARTICLE ACCELERATORS↗

Hadron polarizability data analysis: GoAT

The A2 Collaboration at the Institute for Nuclear Physics in Mainz, Germany, is working towards determining the polarizabilities of hadrons from nonperturbative quantum chromodynamics through Compton scattering experiments at low energies. The asymmetry observables are directly related to the scalar and spin polarizabilities of the hadrons. Online analysis software, which will give real-time feedback on asymmetries, efficiencies, energies, and angle distributions, has been developed. The new software is a big improvement over the existing online code and will greatly develop the quality of the acquired data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

CHEMICAL EVOLUTION LIBRARY FOR GALAXY FORMATION SIMULATION

We have developed a software library for chemical evolution simulations of galaxy formation under the simple stellar population (SSP) approximation. In this library, all of the necessary components concerning chemical evolution, such as initial mass functions, stellar lifetimes, yields from Type II and Type Ia supernovae, asymptotic giant branch stars, and neutron star mergers, are compiled from the literature. Various models are pre-implemented in this library so that users can choose their favorite combination of models. Subroutines of this library return released energy and masses of individual elements depending on a given event type. Since the redistribution manner of these quantities depends on the implementation of users’ simulation codes, this library leaves it up to the simulation code. As demonstrations, we carry out both one-zone, closed-box simulations and 3D simulations of a collapsing gas and dark matter system using this library. In these simulations, we can easily compare the impact of individual models on the chemical evolution of galaxies, just by changing the control flags and parameters of the library. Since this library only deals with the part of chemical evolution under the SSP approximation, any simulation codes that use the SSP approximation—namely, particle-base and mesh codes, as well as semianalytical models—can use it. This library is named “CELib” after the term “Chemical Evolution Library” and is made available to the community.

79 ASTRONOMY AND ASTROPHYSICS↗