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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

SU-E-I-07: An Improved Technique for Scatter Correction in PET

Purpose: In positron emission tomography (PET), the single scatter simulation (SSS) algorithm is widely used for scatter estimation in clinical scans. However, bias usually occurs at the essential steps of scaling the computed SSS distribution to real scatter amounts by employing the scatter-only projection tail. The bias can be amplified when the scatter-only projection tail is too small, resulting in incorrect scatter correction. To this end, we propose a novel scatter calibration technique to accurately estimate the amount of scatter using pre-determined scatter fraction (SF) function instead of the employment of scatter-only tail information. Methods: As the SF depends on the radioactivity distribution and the attenuating material of the patient, an accurate theoretical relation cannot be devised. Instead, we constructed an empirical transformation function between SFs and average attenuation coefficients based on a serious of phantom studies with different sizes and materials. From the average attenuation coefficient, the predicted SFs were calculated using empirical transformation function. Hence, real scatter amount can be obtained by scaling the SSS distribution with the predicted SFs. The simulation was conducted using the SimSET. The Siemens Biograph™ 6 PET scanner was modeled in this study. The Software for Tomographic Image Reconstruction (STIR) was employed to estimate the scatter and reconstruct images. The EEC phantom was adopted to evaluate the performance of our proposed technique. Results: The scatter-corrected image of our method demonstrated improved image contrast over that of SSS. For our technique and SSS of the reconstructed images, the normalized standard deviation were 0.053 and 0.182, respectively; the root mean squared errors were 11.852 and 13.767, respectively. Conclusion: We have proposed an alternative method to calibrate SSS (C-SSS) to the absolute scatter amounts using SF. This method can avoid the bias caused by the insufficient tail information and therefore improve the accuracy of scatter estimation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

Reduced dose measurement of absolute myocardial blood flow using dynamic SPECT imaging in a porcine model

Purpose: Absolute myocardial blood flow (MBF) and myocardial flow reserve (MFR) measurements provide important additional information over traditional relative perfusion imaging. Recent advances in camera technology have made this possible with single-photon emission tomography (SPECT). Low dose protocols are desirable to reduce the patient radiation risk; however, increased noise may reduce the accuracy of MBF measurements. The authors studied the effect of reducing dose on the accuracy of dynamic SPECT MBF measurements. Methods: Nineteen 30–40 kg pigs were injected with 370 + 1110 MBq of Tc-99m sestamibi or tetrofosmin or 37 + 111 MBq of Tl-201 at rest + stress. Microspheres were injected simultaneously to measure MBF. The pigs were imaged in list-mode for 11 min starting at the time of injection using a Discovery NM 530c camera (GE Healthcare). Each list file was modified so that 3/4, 1/2, 1/4, 1/8, 1/16, and 1/32 of the original counts were included in the projections. Modified projections were reconstructed with CT-based attenuation correction and an energy window-based scatter correction and analyzed with FlowQuant kinetic modeling software using a 1-compartment model. A modified Renkin-Crone extraction function was used to convert the tracer uptake rate K1 to MBF values. The SPECT results were compared to those from microspheres. Results: Correlation between SPECT and microsphere MBF values for the full injected activity was r ≥ 0.75 for all 3 tracers and did not significantly degrade over all count levels. The mean MBF and MFR and the standard errors in the estimates were not significantly worse than the full-count data at 1/4-counts (Tc99m-tracers) and 1/2-counts (Tl-201). Conclusions: Dynamic SPECT measurement of MBF and MFR in pigs can be performed with 1/4 (Tc99m-tracers) or 1/2 (Tl-201) of the standard injected activity without significantly reducing accuracy and precision.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

SU-E-I-66: Radiomics and Image Registration Updates for the Computational Environment for Radiotherapy Research (CERR)

Purpose: To present new tools in CERR for Radiomics, image registration and other software updates and additions. Methods: Radiomics: CERR supports generating 3-D texture metrics based on gray scale co-occurance. Two new ways to calculate texture features were added: (1) Local Texture Averaging: Local texture is calculated around a voxel within the userdefined bounding box. The final texture metrics are the average of local textures for all the voxels. This is useful to detect any local texture patterns within an image. (2) Image Smoothing: A convolution ball of user-defined radius is rolled over an image to smooth out artifacts. The texture metrics are then computed on the smooth image. Image Registration: (1) Support was added to import deformation vector fields as well as non-deformable transformation matrices generated by vendor software and stored in standard DICOM format. (2) Support was added to use image within masks while computing image deformations. CT to MR registration is supported. This registration uses morphological edge information within the images to guide the deformation process. In addition to these features, other noteworthy additions to CERR include (1) Irregularly shaped ROI: This is done by taking intersection between infinitely extended irregular polygons drawn on any of the two views. Such an ROI is more conformal and useful in avoiding any unwanted parts of images that cannot be avoided with the conventional cubic box. The ROI is useful to generate Radiomics metrics. (2) Ability to insert RTDOSE in DICOM format to existing CERR plans. (3) Ability to import multi-frame PET-CT and SPECT-CT while maintaining spatial registration between the two modalities. (4) Ability to compile CERR on Unix-like systems. Results: The new features and updates are available via https://www.github.com/adityaapte/cerr . Conclusion: Features added to CERR increase its utility in Radiomics, Image-Registration and Outcomes modeling.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2016

Critical and Subcritical 0-Power Experiment at Rensselaer (CaSPER)

This report discusses the 0-power experiment at Rensselaer Polytechnic Institute (CaSPER). K eff simulation results, list-mode multiplication results, and related work are included. The aim of the work is subcritical measurements for code and nuclear data validation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

Numerical Simulation and Model Analyses of In-situ Biogeochemical Remediation of Chromium at the LANL - 18306

Environmental management problems associated with contamination of groundwater resources with heavy metals (such as chromium) typically require complex biogeochemical model simulations and analysis of in situ remediation. Frequently, the in situ remediation is the only viable treatment option. In these cases, a three-dimensional reactive transport model of relevant biogeochemical processes is invaluable. Recently, we have developed a novel unique biogeochemical simulator called CHROTRAN (http://chrotran.lanl.gov; Hansen et al. 2017, Vesselinov et al. 2018a) capable to represent the full three-dimensional dynamics of processes associated with in situ biogeochemical remediation of heavy metals in heterogeneous aquifers. CHROTRAN accounts for spatial and temporal transients of (1) the heavy metal to be remediated, (2) introduced amendments (e.g., through well injection), and (3) biomass growth and decay. The introduced amendments can include (1) an electron donor (e.g., molasses), (2) a nontoxic-conservative bio-inhibitor (e.g., ethanol), and (3) a biocide (e.g., dithionite). In addition, direct abiotic reduction by donor-metal interaction as well as donor-driven biomass growth and bio-reduction are modeled in CHROTRAN. Other critical processes, such as donor sorption, bio-fouling and biomass death are also modeled. CHROTRAN can simulate the groundwater flow dynamics introduced due a series of injection and extraction wells and account for reduction of aquifer permeability near the wells due bio clogging. CHROTRAN is based on the existing PFLOTRAN (http://pflotran.org; Lichtner et al, 2017a) code framework, taking advantage of the reaction interface allowing for implementation of complex model features and leverage other facets of PFLOTRAN, such as multi-component geochemistry and high-performance computing. Here, we present numerical simulations and model analyses related to the Los Alamos National Laboratory (LANL) chromium site. The analyses are performed by coupling of CHROTRAN with MADS (Model Analysis and Decision Support; http://mads.lanl.gov; Vesselinov et al., 2018b). We present analyses related to uncertainty quantification of model predictions, sensitivity of model parameters and remedial decision support. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Open Chemistry, JupyterLab, REST, and quantum chemistry

Quantum chemistry must evolve if it wants to fully leverage the benefits of the internet age, where the worldwide web offers a vast tapestry of tools that enable users to communicate and interact with complex data at the speed and convenience of a button press. The Open Chemistry project has developed an open-source framework that offers an end-to-end solution for producing, sharing, and visualizing quantum chemical data interactively on the web using an array of modern tools and approaches. These tools build on some of the best open-source community projects such as Jupyter for interactive online notebooks, coupled with 3D accelerated visualization, state-of-the-art computational chemistry codes including NWChem and Psi4, and emerging machine learning and data mining tools such as ChemML and ANI. They offer flexible formats to import and export data, along with approaches to compare computational and experimental data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

Three dosimetry models of lipoma arborescens treated by {sup 90}Y synovectomy

Purpose: Lipoma arborescens (LA) is a benign intra-articular lipomatous proliferation of the synovial membrane. This extremely rare condition has previously been treated by intra-articular{sup 90}Y radiosynoviorthesis but dosimetry literature on this form of radionuclide therapy is nonexistent. The authors detail methodology for successful treatment of LA and provide for the first time estimates of radiation dosimetry. The authors also analyze the biodistribution of the radiopharmaceutical over the course of the patient's treatment through sequential imaging. Methods: A patient with bilateral LA underwent intracavity injection of{sup 90}Y citrate colloid to the right and left knee joint spaces (181 and 198 MBq, respectively). SPECT/CT datasets were acquired over 9 days to quantify the biodistribution and kinetics of the radiopharmaceutical. Radiation dosimetry was performed using the MIRD schema (through OLINDA software), a custom voxel-based method, and a direct Monte Carlo calculation (OEDIPE). Results: Follow-up MRI showed marked reduction in LA size in both knees. Mean absorbed doses to the LA were 21.2 ± 0.8 and 42.9 ± 2.3 Gy using OLINDA, 8.1 ± 0.3 and 16.7 ± 0.5 Gy using voxel based methodology, and 8.2 ± 0.3 and 15.7 ± 0.5 Gy for OEDIPE in the right and left LA, respectively. Distribution of the radiopharmaceutical within the joint space alters over the imaging period, with less than 1% of the remaining activity having moved posteriorly in the knee cavity. No uptake was detected outside of the joint space after assessment with whole-body scintigraphy. Conclusions: An activity of approximately 185 MBq successfully relieved clinical symptoms of LA. There was good correlation between direct Monte Carlo and voxel based techniques, but OLINDA was shown to overestimate the absorbed dose to the tumor. Accurate dosimetry may help select an activity more tailored to the specific size and location of the LA.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2016

Extending the Coyote emulator to dark energy models with standard w {sub 0}- w {sub a} parametrization of the equation of state

We discuss an extension of the Coyote emulator to predict non-linear matter power spectra of dark energy (DE) models with a scale factor dependent equation of state of the form w = w {sub 0}+(1- a ) w {sub a} . The extension is based on the mapping rule between non-linear spectra of DE models with constant equation of state and those with time varying one originally introduced in ref. [40]. Using a series of N-body simulations we show that the spectral equivalence is accurate to sub-percent level across the same range of modes and redshift covered by the Coyote suite. Thus, the extended emulator provides a very efficient and accurate tool to predict non-linear power spectra for DE models with w {sub 0}- w {sub a} parametrization. According to the same criteria we have developed a numerical code that we have implemented in a dedicated module for the CAMB code, that can be used in combination with the Coyote Emulator in likelihood analyses of non-linear matter power spectrum measurements. All codes can be found at https://github.com/luciano-casarini/pkequal.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Model Exploration of an Information-Based Healthcare Intervention Using Parallelization and Active Learning

This paper describes the application of a large-scale active learning method to characterize the parameter space of a computational agent-based model developed to investigate the impact of CommunityRx, a clinical information-based health intervention that provides patients with personalized information about local community resources to meet basic and self-care needs. Additionally, the diffusion of information about community resources and their use is modeled via networked interactions and their subsequent effect on agents' use of community resources across an urban population. A random forest model is iteratively fitted to model evaluations to characterize the model parameter space with respect to observed empirical data. We demonstrate the feasibility of using high-performance computing and active learning model exploration techniques to characterize large parameter spaces; by partitioning the parameter space into potentially viable and non-viable regions, we rule out regions of space where simulation output is implausible to observed empirical data. We argue that such methods are necessary to enable model exploration in complex computational models that incorporate increasingly available micro-level behavior data. We provide public access to the model and high-performance computing experimentation code.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Plant-wide simulation model for transient studies on the XE-100 - 18464

This paper gives an overview of the plant-wide simulation model that is used by X-energy LLC for developing their plant design and control philosophy, as well as assessing plant upset conditions. The simulation model was created in Flownex{sup R} thermal-fluid system simulation software, which solves mass, momentum and energy equations for compressible, incompressible and homogeneous two-phase flow. The plant model encompasses all major components, from reactor to condenser, including main control loop elements. The paper pays particular attention to the validation of the two-phase flow framework used for the secondary steam cycle. Two steady-state validation cases are done using operational data from the THTR. Validation of transient operation is done by simulating density wave oscillations, a type of system-level boiling phenomenon. This is regarded as a thorough test of transient capability because of the challenges this phenomenon poses to several aspects of the calculation and to the consistency of the code as a whole. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

ColDICE: A parallel Vlasov–Poisson solver using moving adaptive simplicial tessellation

Resolving numerically Vlasov–Poisson equations for initially cold systems can be reduced to following the evolution of a three-dimensional sheet evolving in six-dimensional phase-space. We describe a public parallel numerical algorithm consisting in representing the phase-space sheet with a conforming, self-adaptive simplicial tessellation of which the vertices follow the Lagrangian equations of motion. The algorithm is implemented both in six- and four-dimensional phase-space. Refinement of the tessellation mesh is performed using the bisection method and a local representation of the phase-space sheet at second order relying on additional tracers created when needed at runtime. In order to preserve in the best way the Hamiltonian nature of the system, refinement is anisotropic and constrained by measurements of local Poincaré invariants. Resolution of Poisson equation is performed using the fast Fourier method on a regular rectangular grid, similarly to particle in cells codes. To compute the density projected onto this grid, the intersection of the tessellation and the grid is calculated using the method of Franklin and Kankanhalli [65–67] generalised to linear order. As preliminary tests of the code, we study in four dimensional phase-space the evolution of an initially small patch in a chaotic potential and the cosmological collapse of a fluctuation composed of two sinusoidal waves. We also perform a “warm” dark matter simulation in six-dimensional phase-space that we use to check the parallel scaling of the code.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Time Dependent Data Mining in RAVEN

RAVEN is a generic software framework to perform parametric and probabilistic analysis based on the response of complex system codes. The goal of this type of analyses is to understand the response of such systems in particular with respect their probabilistic behavior, to understand their predictability and drivers or lack of thereof. Data mining capabilities are the cornerstones to perform such deep learning of system responses. For this reason static data mining capabilities were added last fiscal year (FY 15). In real applications, when dealing with complex multi-scale, multi-physics systems it seems natural that, during transients, the relevance of the different scales, and physics, would evolve over time. For these reasons the data mining capabilities have been extended allowing their application over time. In this writing it is reported a description of the new RAVEN capabilities implemented with several simple analytical tests to explain their application and highlight the proper implementation. The report concludes with the application of those newly implemented capabilities to the analysis of a simulation performed with the Bison code.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

COG11.1 description, new features, and development activities - 14521

The latest COG code version, COG11.1, is now available. New features include: hybrid criticality detector variance reduction (CritDetVR); delayed fission gammas (DFG); nuclear resonance fluorescence (NRF); radiation sources (RadSrc); an updated COG dictionary (COGLEX); and additional nuclear data libraries. Ongoing software quality assurance and code development activities are also described. (authors)

42 ENGINEERING↗
DOE OSTI2015

Calculation of the dynamic first electronic hyperpolarizability β(−ω{sub σ}; ω{sub 1}, ω{sub 2}) of periodic systems. Theory, validation, and application to multi-layer MoS{sub 2}

We describe our implementation of a fully analytical scheme, based on the 2n + 1 rule, for computing the coupled perturbed Hartree Fock and Kohn-Sham dynamic first hyperpolarizability tensor β(−ω{sub σ}; ω{sub 1}, ω{sub 2}) of periodic 1D (polymer), 2D (slab), and 3D (crystal) systems in the CRYSTAL code [R. Dovesi et al., Int. J. Quantum Chem. 114, 1287 (2014)], which utilizes local Gaussian type basis sets. The dc-Pockels (dc-P) and second harmonic generation (SHG) tensors are included as special cases. It is verified that (i) symmetry requirements are satisfied; (ii) using LiF as an example, the infinite periodic polymer result agrees with extrapolated finite oligomer calculations and, likewise, for the build-up to a 2D slab and a 3D crystal; (iii) the values converge to the static case for low frequencies; and (iv) the Bishop-deKee dispersion formulas relating dc-P, SHG, and general processes are reproduced through quartic terms. Preliminary SHG calculations on multi-layer MoS{sub 2} satisfactorily reproduce experimental data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

GEO3D - Three-Dimensional Computer Model of a Ground Source Heat Pump System

This file is the setup file for GEO3D, a computer program written by Jim Menart to simulate vertical wells in conjunction with a heat pump for ground source heat pump (GSHP) systems. This is a very detailed three-dimensional computer model. This program produces detailed heat transfer and temperature field information for a vertical GSHP system.

15 GEOTHERMAL ENERGY↗
DOE OSTI2014

Compartment modeling of dynamic brain PET—The impact of scatter corrections on parameter errors

Purpose: The aim of this study was to investigate the effect of scatter and its correction on kinetic parameters in dynamic brain positron emission tomography (PET) tumor imaging. The 2-tissue compartment model was used, and two different reconstruction methods and two scatter correction (SC) schemes were investigated. Methods: The GATE Monte Carlo (MC) software was used to perform 2 × 15 full PET scan simulations of a voxelized head phantom with inserted tumor regions. The two sets of kinetic parameters of all tissues were chosen to represent the 2-tissue compartment model for the tracer 3′-deoxy-3′-({sup 18}F)fluorothymidine (FLT), and were denoted FLT{sub 1} and FLT{sub 2}. PET data were reconstructed with both 3D filtered back-projection with reprojection (3DRP) and 3D ordered-subset expectation maximization (OSEM). Images including true coincidences with attenuation correction (AC) and true+scattered coincidences with AC and with and without one of two applied SC schemes were reconstructed. Kinetic parameters were estimated by weighted nonlinear least squares fitting of image derived time–activity curves. Calculated parameters were compared to the true input to the MC simulations. Results: The relative parameter biases for scatter-eliminated data were 15%, 16%, 4%, 30%, 9%, and 7% (FLT{sub 1}) and 13%, 6%, 1%, 46%, 12%, and 8% (FLT{sub 2}) for K{sub 1}, k{sub 2}, k{sub 3}, k{sub 4}, V{sub a}, and K{sub i}, respectively. As expected, SC was essential for most parameters since omitting it increased biases by 10 percentage points on average. SC was not found necessary for the estimation of K{sub i} and k{sub 3}, however. There was no significant difference in parameter biases between the two investigated SC schemes or from parameter biases from scatter-eliminated PET data. Furthermore, neither 3DRP nor OSEM yielded the smallest parameter biases consistently although there was a slight favor for 3DRP which produced less biased k{sub 3} and K{sub i} estimates while OSEM resulted in a less biased V{sub a}. The uncertainty in OSEM parameters was about 26% (FLT{sub 1}) and 12% (FLT{sub 2}) larger than for 3DRP although identical postfilters were applied. Conclusions: SC was important for good parameter estimations. Both investigated SC schemes performed equally well on average and properly corrected for the scattered radiation, without introducing further bias. Furthermore, 3DRP was slightly favorable over OSEM in terms of kinetic parameter biases and SDs.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-J-196: A Prototype System for Portal Imaging for Intensity Modulated Neutron Therapy

Purpose: Fast neutron therapy is offered at the University of Washington Medical Center for treatment of selected cancers. The hardware and control systems of the UW Clinical Neutron Therapy System are undergoing upgrades to enable delivery of IMNT. To clinically implement IMNT, dose verification tools need to be developed. We propose a portal imaging system that relies on the creation of positron emitting isotopes ({sup 11}C and {sup 15}O) through (n, 2n) reactions with a PMMA plate placed below the patient. After field delivery, the plate is retrieved from the vault and imaged using a reader that detects the annihilation photons. The pattern of activity produced in the plate provides information to reconstruct the neutron fluence map that can be compared to fluence maps from Monte Carlo (MCNP) simulations to verify treatment delivery. We have previously performed Monte Carlo simulations of the portal imaging system (GATE simulations) and the beam line (MCNP simulations). In this work, initial measurements using a prototype system are presented. Methods: Custom electronics were developed for BGO detectors read out with photomultiplier tubes (previous generation PET detectors from a CTI ECAT 953 scanner). Two detectors were placed in coincidence, with a detector separation of 2 cm. Custom software was developed to create the crystal look up tables and perform a limited angle planar reconstruction with a stochastic normalization. To test the initial capabilities of the system, PMMA squares were irradiated with neutrons at a depth of 1.5 cm and read out using the prototype system. Doses ranging from 10–200 cGy were delivered. Results: Using the prototype system, dose differences in the therapeutic range could be determined. Conclusion: The prototype portal imaging system is capable of detecting neutron doses as low as 10–50 cGy and shows great promise as a patient QA tool for IMNT.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Release of MCNP6.2 and WHISPER-1.1 - Guidance for NCS users

In 2017, a new version of the MCNP6 Monte Carlo code package was released: MCNP6.2. This is the first version of MCNP to include the WHISPER package to support sensitivity-uncertainty based nuclear criticality safety validation. The coding changes to MCNP6.2 physics are relatively insignificant. Corrections to the S(α,β) thermal scattering numerics are generally negligible relative to problem statistics (or, in rare cases, prevent aborts). Similarly, the changes to adjoint-weighting for computing kinetics parameters may result in small differences, generally negligible compared to problem statistics. The change to the MCNP6.2 geometry treatment to correctly handle coincident surfaces in problems with universe/fill features will produce different round-off in the geometry tracking. This will produce differences in results relative to previous versions, but those differences should be small relative to problem statistics, and are not a concern. An immediate benefit of using MCNP6.2 (rather than MCNP6.1) is that the new version is typically 1.5-2 times faster.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗