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,441 records · Page 80
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

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↗
DOE OSTI2014

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↗
DOE OSTI2018

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↗
DOE OSTI2017

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↗
DOE OSTI2015

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 ↗
DOE OSTI2013

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↗
DOE OSTI2019

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↗
DOE OSTI2013

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↗
DOE OSTI2024

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↗
DOE OSTI2010

Recent Developments in Simulations of an Inverse Cyclotron for Intense Muon Beams

A number of recent developments have led to simulations of an inverse cyclotron for cooling intense muon beams for neutrino factories and muon colliders. Such a device could potentially act as a novel beam cooling mechanism for muons, and it would be significantly smaller and cheaper than other cooling channel designs. Realistic designs are still being explored, but the first simulations of particle tracking in the inverse cyclotron, with accumulation in the cyclotron core, have been done with electrostatic simulations in the particle-in-cell code VORPAL. We present an overview of the muon inverse cyclotron concept and recent simulation results.

43 PARTICLE ACCELERATORS↗
DOE OSTI2010

Quantum-capacity-approaching codes for the detected-jump channel

The quantum-channel capacity gives the ultimate limit for the rate at which quantum data can be reliably transmitted through a noisy quantum channel. Degradable quantum channels are among the few channels whose quantum capacities are known. Given the quantum capacity of a degradable channel, it remains challenging to find a practical coding scheme which approaches capacity. Here we discuss code designs for the detected-jump channel, a degradable channel with practical relevance describing the physics of spontaneous decay of atoms with detected photon emission. We show that this channel can be used to simulate a binary classical channel with both erasures and bit flips. The capacity of the simulated classical channel gives a lower bound on the quantum capacity of the detected-jump channel. When the jump probability is small, it almost equals the quantum capacity. Hence using a classical capacity-approaching code for the simulated classical channel yields a quantum code which approaches the quantum capacity of the detected-jump channel.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

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↗
DOE OSTI2011

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↗
DOE OSTI2015

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↗
DOE OSTI2017

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↗
DOE OSTI2012

Statistical uncertainty analysis applied to the DRAGONv4 code lattice calculations and based on JENDL-4 covariance data

In this paper, multi-group microscopic cross-section uncertainty is propagated through the DRAGON (Version 4) lattice code, in order to perform uncertainty analysis on k{infinity} and 2-group homogenized macroscopic cross-sections predictions. A statistical methodology is employed for such purposes, where cross-sections of certain isotopes of various elements belonging to the 172 groups DRAGLIB library format, are considered as normal random variables. This library is based on JENDL-4 data, because JENDL-4 contains the largest amount of isotopic covariance matrixes among the different major nuclear data libraries. The aim is to propagate multi-group nuclide uncertainty by running the DRAGONv4 code 500 times, and to assess the output uncertainty of a test case corresponding to a 17 x 17 PWR fuel assembly segment without poison. The chosen sampling strategy for the current study is Latin Hypercube Sampling (LHS). The quasi-random LHS allows a much better coverage of the input uncertainties than simple random sampling (SRS) because it densely stratifies across the range of each input probability distribution. Output uncertainty assessment is based on the tolerance limits concept, where the sample formed by the code calculations infers to cover 95% of the output population with at least a 95% of confidence. This analysis is the first attempt to propagate parameter uncertainties of modern multi-group libraries, which are used to feed advanced lattice codes that perform state of the art resonant self-shielding calculations such as DRAGONv4. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2011

Study of Heating and Fusion Power Production in ITER Discharges

ITER simulations, in which the temperatures, toroidal angular frequency and currents are evolved, are carried out using the PTRANSP code starting with initial profiles and boundary conditions obtained from TSC code studies. The dependence of heat deposition and current drive on ICRF frequency, number of poloidal modes, beam orientation, number of Monte Carlo particles and ECRH launch angles is studied in order to examine various possibilities and contingencies for ITER steady state and hybrid discharges. For the hybrid discharges, the fusion power production and fusion Q, computed using the Multi-Mode MMM v7.1 anomalous transport model, are compared with those predicted using the GLF23 model. The simulations of the hybrid scenario indicate that the fusion power production at 1000 sec will be approximately 500 MW corresponding to a fusion Q = 10.0. The discharge scenarios simulated aid in understanding the conditions for optimizing fusion power production and in examining measures of plasma performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Using SWIFT Data

SWIFT has been shown to be a promising new diagnostic tool with many applications. SWIFT diagnostics coupled with PDV and other tools represent an exciting new source of data. SIAS is the name for the methodology we use to couple our SWIFT data to calculations for maximum utilization. SIAS, formerly Schlieren Inverse Analysis Software -- now SWIFT Inverse Analysis Software, uses SWIFT data and a hydro code calculation repeated in a genetic algorithm loop to reconstruct pressure history of shock transducer material interface.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗