Search NASA⌕ Search

SEARCH · Search NASA

Results for “C Code”

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,189 records · Page 66

SU-F-T-308: Mobius FX Evaluation and Comparison Against a Commercial 4D Detector Array for VMAT Plan QA

Purpose: To evaluate the use of MobiusFX as a pre-treatment verification IMRT QA tool and compare it with a commercial 4D detector array for VMAT plan QA. Methods: 15 VMAT plan QA of different treatment sites were delivered and measured by traditional means with the 4D detector array ArcCheck (Sun Nuclear corporation) and at the same time measurement in linac treatment logs (Varian Dynalogs files) were analyzed from the same delivery with MobiusFX software (Mobius Medical Systems). VMAT plan QAs created in Eclipse treatment planning system (Varian) in a TrueBeam linac machine (Varian) were delivered and analyzed with the gamma analysis routine from SNPA software (Sun Nuclear corporation). Results: Comparable results in terms of the gamma analysis with 99.06% average gamma passing with 3%,3mm passing rate is observed in the comparison among MobiusFX, ArcCheck measurements, and the Treatment Planning System dose calculated. When going to a stricter criterion (1%,1mm) larger discrepancies are observed in different regions of the measurements with an average gamma of 66.24% between MobiusFX and ArcCheck. Conclusion: This work indicates the potential for using MobiusFX as a routine pre-treatment patient specific IMRT method for quality assurance purposes and its advantages as a phantom-less method which reduce the time for IMRT QA measurement. MobiusFX is capable of produce similar results of those by traditional methods used for patient specific pre-treatment verification VMAT QA. Even the gamma results comparing to the TPS are similar the analysis of both methods show that the errors being identified by each method are found in different regions. Traditional methods like ArcCheck are sensitive to setup errors and dose difference errors coming from the linac output. On the other hand linac log files analysis record different errors in the VMAT QA associated with the MLCs and gantry motion that by traditional methods cannot be detected.

60 APPLIED LIFE SCIENCES↗

SU-F-T-351: Establishing a Workflow for IMRT Pre-Treatment Reviews for NRG-GY006 Clinical Trial

Purpose: To establish a workflow for NRG-GY006 IMRT pre-treatment reviews, incorporating advanced radiotherapy technologies being evaluated as part of the clinical trial. Methods: Pre-Treatment reviews are required for every IMRT case as part of NRG-GY006 (a randomized phase II trial of radiation therapy and cisplatin alone or in combination with intravenous triapine in women with newly diagnosed bulky stage I B2, stage II, IIIB, or IVA cancer of the uterine cervix or stage II-IVA vaginal cancer. The pretreatment review process includes structures review and generating an active bone marrow(ABM)- to be used as an avoidance structure during IMRT optimization- and evaluating initial IMRT plan quality using knowledgeengineering based planning (KBP). Institutions will initially submit their simulation CT scan, structures file and PET/CT to IROC QA center for generating ABM. The ABM will be returned to the institution for use in planning. Institutions will then submit an initial IMRT plan for review and will receive information back following implementation of a KBP algorithm, for use in re-optimization, before submitting the final IMRT used for treatment. Results: ABM structure is generated using MIM vista software (Version 6.5, MIM corporation, Inc.). Here, the planning CT and the diagnostic PET/CT are fused and a sub threshold structure is auto segmented above the mean value of the SUV of the bone marrow. The generated ABM were compared with those generated with other software system (e.g. Velocity, Varian) and Dice coefficient (reflects the overlap of structures) ranged between 80 – 90% was achieved. A KBP model was built in Varian Eclipse TPS using the RapidPlan KBP software to perform plan quality assurance. Conclusion: The workflow for IMRT pretreatment reviews has been established. It represents a major improvement of NRG Oncology clinical trial quality assurance and incorporates the latest radiotherapy technologies as part of NCI clinical trials. This project was supported by grants U24CA180803 (IROC), UG1CA189867 (NCORP), U10CA180868 (NRG Oncology Operations), U10CA180822 (NRG Oncology SDMC) from the National Cancer Institute (NCI) and PA CURE grant.

60 APPLIED LIFE SCIENCES↗

SU-F-T-443: Quantification of Dosimetric Effects of Dental Metallic Implant On VMAT Plans

Purpose: To evaluate the dosimetric impact of metallic implant that correlates with the size of targets and metallic implants and distance in between on volumetric-modulated arc therapy (VMAT) plans for head and neck (H&N) cancer patients with dental metallic implant. Methods: CT images of H&N cancer patients with dental metallic implant were used. Target volumes with different sizes and locations were contoured. Metal artifact regions excluding surrounding critical organs were outlined and assigned with CT numbers close to water (0HU). VMAT plans with half-arc, one-full-arc and two-full-arcs were constructed and same plans were applied to structure sets with and without CT number assignment of metal artifact regions and compared. D95% was utilized to investigate PTV dose coverage and SNC Patient− Software was used for the analysis of dose distribution difference slice by slice. Results: For different targets sizes, variation of PTV dose coverage (Delta-D95%) with and without CT number replacement reduced with larger target volume for all half-arc, one-arc and two-arc VMAT plans even though there were no clinically significant differences. Additionally, there were no significant variations of the maximum percent difference (max.%diff) of dose distribution. With regard to the target location, Delta-D95% and max. %diff dropped with increasing distance between target and metallic implant. Furthermore, half-arc plans showed greater impact than one-arc plans, and two-arc plans had smallest influence for PTV dose coverage and dose distribution. Conclusion: The target size has less correlation of doseimetric impact than the target location relative to metallic implants. Plans with more arcs alleviate the dosimetric effect of metal artifact because of less contribution to the target dose from beams going through the regions with metallic artifacts. Incorrect CT number causes inaccurate dose distribution, therefore appropriately overwriting metallic artifact regions with reasonable CT numbers is recommended. More patient data are collected and under further analysis.

60 APPLIED LIFE SCIENCES↗

Approaches to reducing photon dose calculation errors near metal implants

Purpose: Dose calculation errors near metal implants are caused by limitations of the dose calculation algorithm in modeling tissue/metal interface effects as well as density assignment errors caused by imaging artifacts. The purpose of this study was to investigate two strategies for reducing dose calculation errors near metal implants: implementation of metal-based energy deposition kernels in the convolution/superposition (C/S) dose calculation method and use of metal artifact reduction methods for computed tomography (CT) imaging. Methods: Both error reduction strategies were investigated using a simple geometric slab phantom with a rectangular metal insert (composed of titanium or Cerrobend), as well as two anthropomorphic phantoms (one with spinal hardware and one with dental fillings), designed to mimic relevant clinical scenarios. To assess the dosimetric impact of metal kernels, the authors implemented titanium and silver kernels in a commercial collapsed cone C/S algorithm. To assess the impact of CT metal artifact reduction methods, the authors performed dose calculations using baseline imaging techniques (uncorrected 120 kVp imaging) and three commercial metal artifact reduction methods: Philips Healthcare’s O-MAR, GE Healthcare’s monochromatic gemstone spectral imaging (GSI) using dual-energy CT, and GSI with metal artifact reduction software (MARS) applied. For the simple geometric phantom, radiochromic film was used to measure dose upstream and downstream of metal inserts. For the anthropomorphic phantoms, ion chambers and radiochromic film were used to quantify the benefit of the error reduction strategies. Results: Metal kernels did not universally improve accuracy but rather resulted in better accuracy upstream of metal implants and decreased accuracy directly downstream. For the clinical cases (spinal hardware and dental fillings), metal kernels had very little impact on the dose calculation accuracy (<1.0%). Of the commercial CT artifact reduction methods investigated, the authors found that O-MAR was the most consistent method, resulting in either improved dose calculation accuracy (dental case) or little impact on calculation accuracy (spine case). GSI was unsuccessful at reducing the severe artifacts caused by dental fillings and had very little impact on calculation accuracy. GSI with MARS on the other hand gave mixed results, sometimes introducing metal distortion and increasing calculation errors (titanium rectangular implant and titanium spinal hardware) but other times very successfully reducing artifacts (Cerrobend rectangular implant and dental fillings). Conclusions: Though successful at improving dose calculation accuracy upstream of metal implants, metal kernels were not found to substantially improve accuracy for clinical cases. Of the commercial artifact reduction methods investigated, O-MAR was found to be the most consistent candidate for all-purpose CT simulation imaging. The MARS algorithm for GSI should be used with caution for titanium implants, larger implants, and implants located near heterogeneities as it can distort the size and shape of implants and increase calculation errors.

60 APPLIED LIFE SCIENCES↗

Data for Responsiveness of Miscanthus and Switchgrass Yields to Stand Age and Nitrogen Fertilization: A Meta-regression Analysis

The compiled datasets include plot level observations of energy crops (miscanthus and switchgrass) from recent experimental field trials in the US including dry biomass yield, location, state, region, harvest year, growing season degree days (GDD), winter season heating degree days (HDD), growing season cumulative precipitation, annual nitrogen application rate, age of the pant when harvested, National Commodity Crop Productivity Index (NCCPI) values, and cultivar type (switchgrass) from various published and unpublished sources. The stata codes include estimation procedures for four different specifications, i.e., Model A includes deterministic effect without interaction terms; Model B includes deterministic effect with interaction terms (N2, age2, N × age, GDD2, precip2, N × NCCPI); Model C includes deterministic effect with interaction terms, study, and location random effect; Model D includes deterministic effect with interaction terms, harvest year augmented study, and location random effect.

Age↗

Synthesis and Crystal Structure of the Zintl Phases NaSrSb, NaBaSb and NaEuSb

This work details the synthesis and the crystal structures of the ternary compounds NaSrSb, NaBaSb and NaEuSb. They are isostructural and adopt the hexagonal ZrNiAl-type structure (space group P6¯2m; Pearson code hP9). The structure determination in all three cases was performed using single-crystal X-ray diffraction methods. The structure features isolated Sb3– anions arranged in layers stacked along the crystallographic c-axis. In the interstices, alkali and alkaline-earth metal cations are found in tetrahedral and square pyramidal coordination environments, respectively. The formal partitioning of the valence electrons adheres to the valence rules, i.e., Na+Sr2+Sb3–, Na+Ba2+Sb3– and Na+Eu2+Sb3– can be considered as Zintl phases with intrinsic semiconductor behavior. Electronic band structure calculations conducted for NaBaSb are consistent with this notion and show a direct gap of approx. 0.9 eV. Additionally, the calculations hint at possible inverted Dirac cones, a feature that is reminiscent of topological quantum materials.

36 MATERIALS SCIENCE↗

Competing Use Analysis for Offshore Renewables and Subsea Cables

Demand for abundant and diverse resources in the oceans is growing, necessitating marine spatial planning. To inform development of Marine Hydrokinetic (MHK) and Offshore Wind (OSW) resources, the Department of Energy (DOE) has asked NREL to identify the competing uses areas between promising MHK/OSW sites and submarine power and telecommunications cables. The first step in this work is to identify and quantify the overlap between the MHK/OSW resource availability and existing cable routes. The analysis is done in terms of resource area because the task of quantifying actual impacts on available resource is a non-trivial undertaking that involves subjective decisions of identifying resource opportunities. Quantifying overlap in-terms of resource area?on the other hand?is significantly more straight forward, and useful to marine spatial planners. This dataset of polygons describes recommended setback areas from submarine cables for marine renewable energy Development. Two archive files are provided: 1. NREL Cables Setback-Polygons.zip contains the following GeoJSON files: 1.a. cables_100ft.geojson: The Code of Federal Regulations (CFR 585.301) specifies that the legal right of way for submarine cables is 100 ft (~30 meters) to either side of the cable (i.e., 200 ft wide). 1.b. cables_2z.geojson: The International Cable Protection Committee (ICPC) of the North American Submarine Cable Association (NASCA) recommends setback distances for new facilities as the maximum of 500 m or twice the bottom depth (2z), per ICPC Recommendation 13 No. 2 (Communications Security, Reliability and Interoperability Council IV 2014). For depths 250 m, 2 * depth is to be used. 1.c. cables_3z.geojson: The ICPC of NASCA recommends setback distances for new cables as the maximum of 500 m or three times the bottom depth (3z), per ICPC Recommendation 2 No. 10 (Communications Security, Reliability and Interoperability Council IV 2014). So for depths 167 m, 3 * depth is to be used. 2. NREL Submarine Cable Analysis Code Data Presentations.zip contains all files (code, data, README, etc) from release 0.1 of this Github repository for possibly regenerating the product as submarine cable configurations change: https://github.com/ecoquants/nrel-cables. For more details, please see the link to the full report for which this data product is a supplemental output: Submarine Cable Analysis for U.S. Marine Renewable Energy Development by Ben Best and Levi Kilcher National Renewable Energy Laboratory

16 TIDAL AND WAVE POWER↗

Lecture Notes on Criticality Safety Validation Using MCNP & Whisper

Training classes for nuclear criticality safety, MCNP documentation. The need for, and problems surrounding, validation of computer codes and data area considered first. Then some background for MCNP & Whisper is given--best practices for Monte Carlo criticality calculations, neutron spectra, S(α,β) thermal neutron scattering data, nuclear data sensitivities, covariance data, and correlation coefficients. Whisper is computational software designed to assist the nuclear criticality safety analyst with validation studies with the Monte Carlo radiation transport package MCNP. Whisper's methodology (benchmark selection – C k 's, weights; extreme value theory – bias, bias uncertainty; MOS for nuclear data uncertainty – GLLS) and usage are discussed.

42 ENGINEERING↗

HopBox: An image analysis pipeline to characterize hop cone morphology

Abstract Hop cone morphology can influence picking and drying ability, and color can impact consumer preference and may be indicative of quality. However, these characteristics are not generally evaluated in hop breeding programs due to the tedious nature of trait quantification and the extensive variation among cones within a genotype. We developed the HopBox, which is a simply constructed light box with a camera mount, and a publicly available image processing pipeline that identifies hop cones within color‐corrected images, reads a QR code within the image, and outputs data on hop cone length, width, area, perimeter, openness, weight, color, and density. The trained model was applied to images of 500 cones each from 15 replicated advanced hop genotypes from the USDA‐ARS breeding program in Prosser, Washington. Analysis of variance revealed significant ( p < 0.001) differences between genotypes for all traits measured, enabling breeders to discriminate between genotypes for selection purposes. Broad sense heritability for all traits ranged from 0.23 to 0.59. A random sampling of hop cones from the complete dataset revealed that imaging only 5–10 cones adequately captured genotypic variation and provided acceptable rank correlations ( r s > 0.75); however, increasing the sample size to 30 provided optimal precision. Instructions for constructing a HopBox and the code for the analysis pipeline are publicly available online and have wide applicability for hop breeding and research.

Altendorf, Kayla R.↗

Feasibility and performance of the staged Z-pinch: A one-dimensional study with FLASH and MACH2

Z-pinch platforms constitute a promising pathway to fusion energy research. Here, we present a one-dimensional numerical study of the staged Z-pinch (SZP) concept using the FLASH and MACH2 codes. We discuss the verification of the codes using two analytical benchmarks that include Z-pinch-relevant physics, building confidence on the codes' ability to model such experiments. Then, FLASH is used to simulate two different SZP configurations: a xenon gas-puff liner (SZP1*) and a silver solid liner (SZP2). The SZP2 results are compared against previously published MACH2 results, and a new code-to-code comparison on SZP1* is presented. Using an ideal equation of state and analytical transport coefficients, FLASH yields a fuel convergence ratio (CR) of approximately 39 and a mass-averaged fuel ion temperature slightly below 1 keV for the SZP2 scheme, significantly lower than the full-physics MACH2 prediction. For the new SZP1* configuration, full-physics FLASH simulations furnish large and inherently unstable CRs (>300) but achieve fuel ion temperatures of many kilo-electron volts. While MACH2 also predicts high temperatures, the fuel stagnates at a smaller CR. The integrated code-to-code comparison reveals how magnetic insulation, heat conduction, and radiation transport affect platform performance and the feasibility of the SZP concept.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Tutorial on Bayesian analysis of linear shock compression data

Gas gun and other shock compression experiments often produce shock wave velocity measurements that are linearly associated with particle velocity. Traditionally, this empirical relationship is quantified with a single Hugoniot curve that is estimated using least squares regression. However, for downstream modeling and simulation tasks, it is often more useful to have multiple Hugoniot curves in the pressure–volume plane that are consistent with the data. We employ Bayesian uncertainty quantification methods as a framework for propagating measurement uncertainty through to model parameters and predictions. Specifically, this Tutorial shows how to sample multiple Hugoniot curves in the pressure–volume plane that are consistent with the shock wave-particle velocity measurements in a two-step Bayesian approach. First, we obtain an analytical expression for the posterior distribution of the linear model parameters using Bayesian linear regression. Second, we propagate samples from the posterior distribution through the Rankine–Hugoniot equations to yield Hugoniot curves in the pressure–volume plane. The procedure is demonstrated with publicly available data on argon, copper, and nickel, and compared against bootstrapping and linear regression. The Bayesian procedure is shown to be interpretable, computationally inexpensive, and less sensitive than an alternative bootstrapping approach to the removal of the point in the copper dataset that has the largest particle velocity. As a Tutorial on Bayesian methodology for the shock compression community, we provide several derivations and explanations that make this paper self-contained, and make all code and data available at github.com/llnl/BALSCD.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Metrics and extrapolation of resonant magnetic perturbation thresholds for ELM suppression

This large database study of resonant magnetic perturbation (RMP) edge localized mode (ELM) suppression thresholds in the AUG, DIII-D, EAST, and KSTAR tokamaks details the key strengths and weaknesses of RMP metrics. The RMP ELM suppression database used for this work contains plasma information at the time of transition from ELMing to ELM suppressed states where a clear experimental threshold is identified. The experimental threshold distributions are compared for five metrics: (1) the island overlap width, (2) pedestal top Chirikov overlap, (3) peeling edge displacement, (4) pedestal top resonant drive, and (5) edge dominant mode overlap. The distributions, the regularity of the dependence on RMP coil currents, and the sensitivities of a given metric to equilibrium reconstruction details are compared. The overlap metric proves to be a good compromise between including the appropriate plasma response physics and maintaining a numerical robustness. This quantity does not exhibit clear power-law scalings for projection, but machine learning can assist in predicting thresholds within the existing parameter ranges and providing uncertainty quantification of those predictions. Two new first-principles models, one utilizing a threshold from the non-linear Modified Rutherford equation evaluated at the pedestal top and one utilizing the SLAYER code to calculate the linear tearing threshold from torque balance, offer possible paths to extrapolation beyond the existing database parameter space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulations of future particle accelerators: issues and mitigations

The ever increasing demands placed upon machine performance have resulted in the need for more comprehensive particle accelerator modeling. Computer simulations are key to the success of particle accelerators. Many aspects of particle accelerators rely on computer modeling at some point, sometimes requiring complex simulation tools and massively parallel supercomputing. Examples include the modeling of beams at extreme intensities and densities (toward the quantum degeneracy limit), and with ultra-fine control (down to the level of individual particles). In the future, adaptively tuned models might also be relied upon to provide beam measurements beyond the resolution of existing diagnostics. Much time and effort has been put into creating accelerator software tools, some of which are highly successful. However, there are also shortcomings such as the general inability of existing software to be easily modified to meet changing simulation needs. In this paper possible mitigating strategies are discussed for issues faced by the accelerator community as it endeavors to produce better and more comprehensive modeling tools. This includes lack of coordination between code developers, lack of standards to make codes portable and/or reusable, lack of documentation, among others.

43 PARTICLE ACCELERATORS↗

Modelling type 1 quasar colours in the era of Rubin and Euclid

We construct a parametric SED model which is able to reproduce the average observed SDSS–UKIDSS–WISE quasar colours to within one-tenth of a magnitude across a wide range of redshift (0 < z < 5) and luminosity (–22 > M i > –29). This model is shown to provide accurate predictions for the colours of known quasars which are less luminous than those used to calibrate the model parameters, and also those at higher redshifts z > 5. Using a single parameter, the model encapsulates an up-to-date understanding of the intra-population variance in the rest-frame ultraviolet and optical emission lines of luminous quasars. At fixed redshift, there are systematic changes in the average quasar colours with apparent i-band magnitude, which we find to be well explained by the contribution from the host galaxy and our parametrization of the emission-line properties. By including redshift as an additional free parameter, the model could be used to provide photometric redshifts for individual objects. For the population as a whole we find that the average emission line and host-galaxy contributions can be well described by simple functions of luminosity which account for the observed changes in the average quasar colours across 18.1 < i AB < 21.5. Finally, we use these trends to provide predictions for quasar colours at the luminosities and redshifts which will be probed by the Rubin Observatory LSST and ESA-Euclid wide survey. The model code is applicable to a wide range of upcoming photometric and spectroscopic surveys, and is made publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗

Anomalous Absorption by the Two-Plasmon Decay Instability

Radiation-hydrodynamic simulations of directly driven fusion experiments at the Omega Laser Facility predict absorption accurately when targets are driven at low overlapped laser intensity. Discrepancies appear at increased intensity, however, with higher-than-expected laser absorption on target. Strong correlations with signatures of the two-plasmon decay (TPD) instability—including half-harmonic and hard-x-ray emission—indicate that TPD is responsible for this anomalous absorption. Scattered light data suggest that up to ≈ 30 % of the laser power reaching quarter-critical density can be absorbed locally when the TPD threshold is exceeded. A scaling of absorption versus TPD threshold parameter was empirically determined and validated using the laser–plasma simulation environment code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mutually unbiased bases for time-bin qudits

We introduce a method for generation and detection of mutually unbiased bases for time-bin qudits employing electro-optic phase modulator – coded fiber Bragg grating pairs. Our approach uses one spatial mode and can switch rapidly between bases.

Lukens, Joseph↗

Enhanced Visualization Capabilities in OpenMC: Report of Activities under a DOE GAIN Voucher

This work was performed under a Cooperative Research and Development Agreement (CRADA) through the U.S. Department of Energy’s Gateway for Accelerated Innovations in Nuclear (GAIN) program in partnership with Oklo Inc., the industry participant. The broad goal of this project was to improve the visualization capabilities of OpenMC related to model development and simulation results. OpenMC is a Monte Carlo particle transport code that is developed and maintained by a broad community of contributors, of which Argonne National Laboratory plays a leading role.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗