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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Implications for Post-processing Nucleosynthesis of Core-collapse Supernova Models with Lagrangian Particles

We investigate core-collapse supernova (CCSN) nucleosynthesis with self-consistent, axisymmetric (2D) simulations performed using the neutrino hydrodynamics code Chimera. Computational costs have traditionally constrained the evolution of the nuclear composition within multidimensional CCSN models to, at best, a 14-species α-network capable of tracking only (α,γ) reactions from {sup 4}He to {sup 60}Zn. Such a simplified network limits the ability to accurately evolve detailed composition and neutronization or calculate the nuclear energy generation rate. Lagrangian tracer particles are commonly used to extend the nuclear network evolution by incorporating more realistic networks into post-processing nucleosynthesis calculations. However, limitations such as poor spatial resolution of the tracer particles; inconsistent thermodynamic evolution, including misestimation of expansion timescales; and uncertain determination of the multidimensional mass cut at the end of the simulation impose uncertainties inherent to this approach. We present a detailed analysis of the impact of such uncertainties for four self-consistent axisymmetric CCSN models initiated from solar-metallicity, nonrotating progenitors of 12, 15, 20, and 25 M{sub ⊙} and evolved with the smaller α-network to more than 1 s after the launch of an explosion.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Interrelationships between information and energy using knowledge management tools

Edward Fredkin was an enthusiastic advocate of information-based theoretical physics, who, in the early 1980s, proposed a new theory of physics based on the idea that the universe is ultimately composed of software. According to Fredkin, reality should be considered as being composed not of particles, matter and forces or energy but of bits of data or information modified according to computational rules. Fredkin went on to demonstrate that, while energy is necessary for storing and retrieving information, it can be arbitrarily reduced in order to carry out any particular instance of information processing, and this operation does not have a lower bound. This implies that it is information rather than matter or energy that should be considered at the ultimate fundamental constituent of reality. This possibility had already been suggested by other scientists. Norbert Wiener heralded a fundamental shift from energy to information and suggested that the universe was founded essentially on the transformation of information, not energy. However, Konrad Zuse was the first, back in 1967, to defend the idea that a digital computer is computing the universe. Richard P. Feynman showed this possibility in a similar light in his reflections on how information related to matter and energy. Other pioneering research on the theory of digital physics was published by Kantor in 1977 and more recently by Stephen Wolfram in 2002, who thereby joined the host of voices upholding that it is patterns of information, not matter and energy, that constitute the cornerstones of reality. In this paper, we introduce the use of knowledge management tools for the purpose of analysing this topic.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Integrating Energy-Efficient Computing with Computational Research to Accelerate Energy Technology

NREL's computational sciences center hosts the largest high performance computing (HPC) capabilities dedicated to energy research while functioning as a living laboratory for energy-efficient computing. NREL's HPC capabilities support the research needs of the Department of Energy's Office of Energy Efficiency and Renewable Energy (EERE). In ten years of operation, HPC use in EERE-sponsored research has grown by a factor of 30, including work in electricity generation, energy efficiency, transportation, and energy system modeling. This paper analyzes this research portfolio, providing examples of individual use cases. The paper documents NREL's history of operating one of the world's most energy-efficient data centers while examining pathways to reduce economic and environmental impact beyond reduction of Power Usage Efficiency (PUE). This paper concludes by examining the unique opportunities created for accelerating improvements in data center efficiency created by combining an HPC system dedicated to energy research and a research program in energy-efficient computing.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

SU-E-T-61: A Practical Process for Fabricating Passive Scatter Proton Beam Modulation Compensation Filters Using 3D Printing

Purpose: The purpose of this project was to devise a practical fabrication process for passive scatter proton beam compensation filters (CF) that is competitive in time, cost and effort using 3D printing. Methods: DICOM compensator filter files for a proton beam were generated by our Eclipse (Varian, Inc.) treatment planning system. The compensator thickness specifications were extracted with in-house software written in Matlab (MathWorks, Inc.) code and written to a text file which could be read by the Rhinoceros 5, computer-aided design (CAD) package (Robert McNeel and Associates), which subsequently generated a smoothed model in a STereoLithographic also known as a Standard Tesselation Language file (STL). The model in the STL file was subsequently refined using Netfabb software and then converted to printing instructions using Cura. version 15.02.1. for our 3D printer. The Airwolf3D, model HD2x, fused filament fabrication (FFF) 3D printer (Airwolf3D.com) was used for our fabrication system with a print speed of 150mm per second. It can print in over 22 different plastic filament materials in a build volume of 11” x 8” x 12”. We choose ABS plastic to print the 3D model of the imprint for our CFs. Results: Prints of the CF could be performed at a print speed of 70mm per second. The time to print the 3D topology for the CF for the 14 cm diameter snout of our Mevion 250 proton accelerator was less than 3 hours. The printed model is intended to subsequently be used as a mold to imprint a molten wax cylindrical to form the compensation after cooling. The whole process should be performed for a typical 3 beam treatment plan within a day. Conclusion: Use of 3D printing is practical and can be used to print a 3D model of a CF within a few hours.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

MO-DE-207B-11: Reliability of PET/CT Radiomics Features in Functional and Morphological Components of NSCLC Lesions: A Repeatability Analysis in a Prospective Multicenter Cohort

Purpose: The goal of this study was to evaluate the repeatability of radiomics features (intensity, shape and heterogeneity) in both PET and low-dose CT components of test-retest FDG-PET/CT images in a prospective multicenter cohort of 74 NSCLC patients from ACRIN 6678 and a similar Merck trial. Methods: Seventy-four patients with stage III-IV NCSLC were prospectively included. The primary tumor and up to 3 additional lesions per patient were analyzed. The Fuzzy Locally Adaptive Bayesian algorithm was used to automatically delineate metabolically active volume (MAV) in PET. The 3D SlicerTM software was exploited to delineate anatomical volumes (AV) in CT. Ten intensity first-order features, as well as 26 textural features and four 3D shape descriptors were calculated from tumour volumes in both modalities. The repeatability of each metric was assessed by Bland-Altman analysis. Results: One hundred and five lesions (primary tumors and nodal or distant metastases) were delineated and characterized. The MAV and AV determination had a repeatability of −1.4±11.0% and −1.2±18.7% respectively. Several shape and heterogeneity features were found to be highly or moderately repeatable (e.g., sphericity, co-occurrence entropy or intensity size-zone matrix zone percentage), whereas others were confirmed as unreliable with much higher variability (more than twice that of the corresponding volume determination). Conclusion: Our results in this large multicenter cohort with more than 100 measurements confirm the PET findings in previous studies (with <30 lesions). In addition, our study is the first to explore the repeatability of radiomics features in the low-dose CT component of PET/CT acquisitions (previous studies considered dosimetry CT, CE-CT or CBCT). Several features were identified as reliable in both PET and CT components and could be used to build prognostic models. This work has received a French government support granted to the CominLabs excellence laboratory and managed by the National Research Agency in the “Investing for the Future” program under reference ANR-10-LABX-07-01, and support from the city of Brest.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

A PARALLEL MONTE CARLO CODE FOR SIMULATING COLLISIONAL N-BODY SYSTEMS

We present a new parallel code for computing the dynamical evolution of collisional N-body systems with up to N {approx} 10{sup 7} particles. Our code is based on the Henon Monte Carlo method for solving the Fokker-Planck equation, and makes assumptions of spherical symmetry and dynamical equilibrium. The principal algorithmic developments involve optimizing data structures and the introduction of a parallel random number generation scheme as well as a parallel sorting algorithm required to find nearest neighbors for interactions and to compute the gravitational potential. The new algorithms we introduce along with our choice of decomposition scheme minimize communication costs and ensure optimal distribution of data and workload among the processing units. Our implementation uses the Message Passing Interface library for communication, which makes it portable to many different supercomputing architectures. We validate the code by calculating the evolution of clusters with initial Plummer distribution functions up to core collapse with the number of stars, N, spanning three orders of magnitude from 10{sup 5} to 10{sup 7}. We find that our results are in good agreement with self-similar core-collapse solutions, and the core-collapse times generally agree with expectations from the literature. Also, we observe good total energy conservation, within {approx}< 0.04% throughout all simulations. We analyze the performance of the code, and demonstrate near-linear scaling of the runtime with the number of processors up to 64 processors for N = 10{sup 5}, 128 for N = 10{sup 6} and 256 for N = 10{sup 7}. The runtime reaches saturation with the addition of processors beyond these limits, which is a characteristic of the parallel sorting algorithm. The resulting maximum speedups we achieve are approximately 60 Multiplication-Sign , 100 Multiplication-Sign , and 220 Multiplication-Sign , respectively.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Assessment of Mechanical Draft Cooling Tower Thermal Emissions from Visual Images of Plumes

Using a one-dimensional code, we computed the power (enthalpy discharge rate) of a twelve-cell mechanical draft cooling tower (MDCT) using over two hundred visible condensed water vapor plume volume measurements derived from images, weather data, and tower operating conditions. The plume images were simultaneously captured by multiple stationary digital cameras surrounding the cooling tower. An analysis technique combining structure from motion (SfM), a neural-network-based image segmentation algorithm, and space carving was used to quantify the volumes. Afterwards, the power output was computed using novel techniques in the one-dimensional code that included cooling tower exhaust plume adjacency effects implemented with a modified version of the entrainment function, weather data averaged from eleven stations, and fan operations at the times when plume volumes were measured. The model was then compared with the averaged observed power output, and it validated well with an average error ranging from 6 to 12%, depending on the meteorological data used in the simulations. This methodology can possibly determine power plant fuel consumption rates by applying visible imagery.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2017

Optimizing Performance of Combustion Chemistry Solvers on Intel's Many Integrated Core (MIC) Architectures

This work investigates novel algorithm designs and optimization techniques for restructuring chemistry integrators in zero and multidimensional combustion solvers, which can then be effectively used on the emerging generation of Intel's Many Integrated Core/Xeon Phi processors. These processors offer increased computing performance via large number of lightweight cores at relatively lower clock speeds compared to traditional processors (e.g. Intel Sandybridge/Ivybridge) used in current supercomputers. This style of processor can be productively used for chemistry integrators that form a costly part of computational combustion codes, in spite of their relatively lower clock speeds. Performance commensurate with traditional processors is achieved here through the combination of careful memory layout, exposing multiple levels of fine grain parallelism and through extensive use of vendor supported libraries (Cilk Plus and Math Kernel Libraries). Important optimization techniques for efficient memory usage and vectorization have been identified and quantified. These optimizations resulted in a factor of ~ 3 speed-up using Intel 2013 compiler and ~ 1.5 using Intel 2017 compiler for large chemical mechanisms compared to the unoptimized version on the Intel Xeon Phi. The strategies, especially with respect to memory usage and vectorization, should also be beneficial for general purpose computational fluid dynamics codes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Development of Parallel Computing Framework to Enhance Radiation Transport Code Capabilities for Rare Isotope Beam Facility Design

A parallel computing framework has been developed to use with general-purpose radiation transport codes. The framework was implemented as a C++ module that uses MPI for message passing. It is intended to be used with older radiation transport codes implemented in Fortran77, Fortran 90 or C. The module is significantly independent of radiation transport codes it can be used with, and is connected to the codes by means of a number of interface functions. The framework was developed and tested in conjunction with the MARS15 code. It is possible to use it with other codes such as PHITS, FLUKA and MCNP after certain adjustments. Besides the parallel computing functionality, the framework offers a checkpoint facility that allows restarting calculations with a saved checkpoint file. The checkpoint facility can be used in single process calculations as well as in the parallel regime. The framework corrects some of the known problems with the scheduling and load balancing found in the original implementations of the parallel computing functionality in MARS15 and PHITS. The framework can be used efficiently on homogeneous systems and networks of workstations, where the interference from the other users is possible.

43 PARTICLE ACCELERATORS↗
DOE OSTI2019

Code Verification of a Warm Electron Diode Using the EMPIRE Plasma Simulation Code.

This report documents recent code verification exercises for a warm electron diode problem using the EMPIRE plasma code. This computationally expensive test was performed three times, including two different code versions and two different time integration algorithms, and the resulting code responses were analyzed for convergence to the analytical solution and orders-of-convergence using the StREEQ numerical error estimation tool. Significant deviations from the exact solution and expected orders-of-convergence, as well as changes in code behavior over time and due to choice of time integration algorithm were observed, illustrating the need to fix fundamental code issues as well as additional code verification testing before the code can be relied upon to accurately solve critical problems within its application space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

A method to synchronize signals from multiple patient monitoring devices through a single input channel for inclusion in list-mode acquisitions

Purpose: This technical note documents a method that the authors developed for combining a signal to synchronize a patient-monitoring device with a second physiological signal for inclusion into list-mode acquisition. Our specific application requires synchronizing an external patient motion-tracking system with a medical imaging system by multiplexing the tracking input with the ECG input. The authors believe that their methodology can be adapted for use in a variety of medical imaging modalities including single photon emission computed tomography (SPECT) and positron emission tomography (PET). Methods: The authors insert a unique pulse sequence into a single physiological input channel. This sequence is then recorded in the list-mode acquisition along with the R-wave pulse used for ECG gating. The specific form of our pulse sequence allows for recognition of the time point being synchronized even when portions of the pulse sequence are lost due to collisions with R-wave pulses. This was achieved by altering our software used in binning the list-mode data to recognize even a portion of our pulse sequence. Limitations on heart rates at which our pulse sequence could be reliably detected were investigated by simulating the mixing of the two signals as a function of heart rate and time point during the cardiac cycle at which our pulse sequence is mixed with the cardiac signal. Results: The authors have successfully achieved accurate temporal synchronization of our motion-tracking system with acquisition of SPECT projections used in 17 recent clinical research cases. In our simulation analysis the authors determined that synchronization to enable compensation for body and respiratory motion could be achieved for heart rates up to 125 beats-per-minute (bpm). Conclusions: Synchronization of list-mode acquisition with external patient monitoring devices such as those employed in motion-tracking can reliably be achieved using a simple method that can be implemented using minimal external hardware and software modification through a single input channel, while still recording cardiac gating signals.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

SU-E-J-262: Variability in Texture Analysis of Gynecological Tumors in the Context of An 18F-FDG PET Adaptive Protocol

Purpose: This study examines the effect on texture analysis due to variable reconstruction of PET images in the context of an adaptive FDG PET protocol for node positive gynecologic cancer patients. By measuring variability in texture features from baseline and intra-treatment PET-CT, we can isolate unreliable texture features due to large variation. Methods: A subset of seven patients with node positive gynecological cancers visible on PET was selected for this study. Prescribed dose varied between 45–50.4Gy, with a 55–70Gy boost to the PET positive nodes. A baseline and intratreatment (between 30–36Gy) PET-CT were obtained on a Siemens Biograph mCT. Each clinical PET image set was reconstructed 6 times using a TrueX+TOF algorithm with varying iterations and Gaussian filter. Baseline and intra-treatment primary GTVs were segmented using PET Edge (MIM Software Inc., Cleveland, OH), a semi-automatic gradient-based algorithm, on the clinical PET and transferred to the other reconstructed sets. Using an in-house MATLAB program, four 3D texture matrices describing relationships between voxel intensities in the GTV were generated: co-occurrence, run length, size zone, and neighborhood difference. From these, 39 textural features characterizing texture were calculated in addition to SUV histogram features. The percent variability among parameters was first calculated. Each reconstructed texture feature from baseline and intra-treatment per patient was normalized to the clinical baseline scan and compared using the Wilcoxon signed-rank test in order to isolate variations due to reconstruction parameters. Results: For the baseline scans, 13 texture features showed a mean range greater than 10%. For the intra scans, 28 texture features showed a mean range greater than 10%. Comparing baseline to intra scans, 25 texture features showed p <0.05. Conclusion: Variability due to different reconstruction parameters increased with treatment, however, the majority of texture features showed significant changes during treatment independent of reconstruction effects.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

WE-G-209-03: PET

Digital radiography, CT, PET, and MR are complicated imaging modalities which are composed of many hardware and software components. These components work together in a highly coordinated chain of events with the intent to produce high quality images. Acquisition, processing and reconstruction of data must occur in a precise way for optimum image quality to be achieved. Any error or unexpected event in the entire process can produce unwanted pixel intensities in the final images which may contribute to visible image artifacts. The diagnostic imaging physicist is uniquely qualified to investigate and contribute to resolution of image artifacts. This course will teach the participant to identify common artifacts found clinically in digital radiography, CT, PET, and MR, to determine the causes of artifacts, and to make recommendations for how to resolve artifacts. Learning Objectives: Identify common artifacts found clinically in digital radiography, CT, PET and MR. Determine causes of various clinical artifacts from digital radiography, CT, PET and MR. Describe how to resolve various clinical artifacts from digital radiography, CT, PET and MR.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

WE-G-209-04: MRI

Digital radiography, CT, PET, and MR are complicated imaging modalities which are composed of many hardware and software components. These components work together in a highly coordinated chain of events with the intent to produce high quality images. Acquisition, processing and reconstruction of data must occur in a precise way for optimum image quality to be achieved. Any error or unexpected event in the entire process can produce unwanted pixel intensities in the final images which may contribute to visible image artifacts. The diagnostic imaging physicist is uniquely qualified to investigate and contribute to resolution of image artifacts. This course will teach the participant to identify common artifacts found clinically in digital radiography, CT, PET, and MR, to determine the causes of artifacts, and to make recommendations for how to resolve artifacts. Learning Objectives: Identify common artifacts found clinically in digital radiography, CT, PET and MR. Determine causes of various clinical artifacts from digital radiography, CT, PET and MR. Describe how to resolve various clinical artifacts from digital radiography, CT, PET and MR.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Group Ring Ideals Related to Reed–Muller Codes

Reed–Muller codes are one of the most well-studied families of codes; however, there are till open problems regarding their structure. Recently a new ring-theoretic approach has emerged that provides a rather intuitive construction of these codes. This approach is centered around the notion of basic Reed–Muller codes. It is known that basic Reed–Muller codes ℳ{sub π}(m, k) over a prime field are powers of the radical RS of a corresponding group algebra and over a nonprime field there are no such equalities, except for trivial ones. In this paper, we consider the ideals ℜ{sub S}ℳ{sub π}(m, k) that arise while studying the inclusions of the basic codes and radical powers.

97 MATHEMATICS AND COMPUTING↗