Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time scale separation”

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,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

Variability of Springtime Transpacific Pollution Transport During 2000-2006: The INTEX-5 Mission in the Context of Previous Years

We analyze the transport of pollution across the Pacific during the NASA INTEX-B (Intercontinental Chemical Transport Experiment Part 8) campaign in spring 2006 and examine how this year compares to the time period for 2000 through 2006. In addition to aircraft measurements of carbon monoxide (CO) collected during INTEX-B, we include in this study multi-year satellite retrievals of CO from the Measurements of Pollution in the Troposphere (MOPITT) instrument and simulations from the chemistry transport model MOZART-4. Model tracers are used to examine the contributions of different source regions and source types to pollution levels over the Pacific. Additional modeling studies are performed to separate the impacts of inter-annual variability in meteorology and .dynamics from changes in source strength. interannual variability in the tropospheric CO burden over the Pacific and the US as estimated from the MOPITT data range up to 7% and a somewhat smaller estimate (5%) is derived from the model. When keeping the emissions in the model constant between years, the year-to-year changes are reduced (2%), but show that in addition to changes in emissions, variable meteorological conditions also impact transpacific pollution transport. We estimate that about 113 of the variability in the tropospheric CO loading over the contiguous US is explained by changes in emissions and about 213 by changes in meteorology and transport. Biomass burning sources are found to be a larger driver for inter-annual variability in the CO loading compared to fossil and biofuel sources or photochemical CO production even though their absolute contributions are smaller. Source contribution analysis shows that the aircraft sampling during INTEX-B was fairly representative of the larger scale region, but with a slight bias towards higher influence from Asian contributions.

Pfister, G. G.↗
NASA NTRS2016

CGRO/BATSE Data Support the New Paradigm For GRB Prompt Emission and the New L-i(nTh)-E-peak,i(nTh,rest) Relation

The paradigm for gamma-ray burst (GRB) prompt emission is changing. Since early in the Compton Gamma RayObservatory (CGRO) era, the empirical Band function has been considered a good description of the keV-MeV-gamma-ray prompt emission spectra despite the fact that its shape was very often inconsistent with the theoretical predictions, especially those expected in pure synchrotron emission scenarios. We have recently established a new observational model analyzing data of the NASA Fermi Gamma-ray Space Telescope. In this model, GRB prompt emission would be a combination of three main emission components: (i) a thermal-like component that we have interpreted so far as emission from the jet photosphere, (ii) a non-thermal component that we have interpreted so far as either synchrotron radiation from the propagating and accelerated charged particles within the jet or reprocessed jet photospheric emission, and (iii) an additional non-thermal (cutoff) power law (PL) extending from low to high energies in gamma-rays and most likely of inverse Compton origin. In this article we reanalyze some of the bright GRBs, namely GRBs 941017, 970111, and 990123, observed with the Burst And Transient Source Experiment (BATSE) on board CGRO with the new model. We conclude that BATSE data for these three GRBs are fully consistent with the recent results obtained with Fermi: some bright BATSE GRBs exhibit three separate components during the prompt phase with similar spectral parameters as those reported from Fermi data. In addition, the analysis of the BATSE GRBs with the new prompt emission model results in a relation between the time-resolved energy flux of the non-thermal component, F(in)(Th), and its corresponding nuFnu spectral peak energy,Epeak,inTh (i.e., FinThEpeak,inTh ), which has a similar index when fitted to a PL as the one initially derived from Fermi data. For GRBs with known redshifts (z) this results in a possible universal relation between the luminosity of the non-thermal component, LinTh, and its corresponding nuFnu spectral peak energy in the rest frame, Epeak,iNT,rest (i.e.,LinThEpeak,iNT,rest). We estimated the redshifts of GRBs 941017 and 970111 using GRB 990123 with z = 1.61as a reference. The estimated redshift for GRB 941017 is typical for long GRBs and the estimated redshift for GRB970111 is right in the range of the expected values for this burst.

acceleration of particles – black hole physics â↗
DOE OSTI2022

Mining Air for Fuels and Fine Chemicals

The Intergovernmental Panel on Climate Change’s 2021 report underscores that Direct Air Capture (DAC) is an essential component to limiting global warming to 2 °C, yet separating CO 2 from air and its subsequent concentration for conversion into carbon neutral fuels/chemicals or permanent storage remains energy intensive and costly. Further, most DAC technologies rely on repetitive capture and release cycles that reduces the amount of time the DAC materials are able to collect CO 2 and increases wear and tear. Moisture swing (MS) sorption within anion exchange materials utilizes differences in relative humidity to capture and concentrate CO 2 up to 500-fold with negligible energy inputs and low-cost sorbents. Molecular modeling predicts that membranes with a dry and wet side would actively transport CO 2 from its dry side to its wet side against a counter flow of H 2 O evaporating on the dry side, however, existing MS materials are too brittle to be processed into membranes to test this hypothesis. This project brought together experts in DAC materials and design, techno-economic analysis and project management at Arizona State University (ASU), polymer processing at University of Texas at Austin (UTA) and modeling gas transport and exchange at Norther Arizona University (NAU) to develop hollow fiber membranes that use energy from water evaporation to continuously pump CO 2 from air against a concentration gradient to the membranes interior and requiring significantly less energy than current DAC technologies. While this project did not meet its ultimate objectives of demonstrating an active CO 2 pumping membrane, a number of key accomplishments were made in developing low-cost, flexible anion exchange membranes (AEM) at 3.5 m 2 scale, characterizing the ionic and molecular transport within MS materials, and developing several analytical and mathematical models of the MS and pumping process leading to new fundamental knowledge about key rate limiting steps of CO 2 flux in anion exchange membranes at low water activity. Technoeconomic models show that if CO 2 fluxes ≥ 25 µmol CO 2 m -2 s -1 (membrane surface area) are achieved that cost ≤ $\$$100/tonne CO 2 captured, purified and compressed to pipeline quality is possible. The knowledge gained in this research will guide future research in developing new materials toward overcoming these rate limiting steps that in turn could enable transformative and disruptive DAC technologies needed for capturing gigatons of CO 2 per year needed to limiting global warming to 2 °C.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2018

How Can TOLNet Help to Better Understand Tropospheric Ozone? A Satellite Perspective

Potential sources of a priori ozone (O3) profiles for use in Tropospheric Emissions: Monitoring of Pollution (TEMPO) satellite tropospheric O3 retrievals are evaluated with observations from multiple Tropospheric Ozone Lidar Network (TOLNet) systems in North America. An O3 profile climatology (tropopause-based O3 climatology (TB-Clim), currently proposed for use in the TEMPO O3 retrieval algorithm) derived from ozonesonde observations and O3 profiles from three separate models (operational Goddard Earth Observing System (GEOS-5) Forward Processing (FP) product, reanalysis product from Modern-Era Retrospective analysis for Research and Applications version 2 (MERRA2), and the GEOS-Chem chemical transport model (CTM)) were: 1) evaluated with TOLNet measurements on various temporal scales (seasonally, daily, hourly) and 2) implemented as a priori information in theoretical TEMPO tropospheric O3 retrievals in order to determine how each a priori impacts the accuracy of retrieved tropospheric (0-10 km) and lowermost tropospheric (LMT, 0-2 km) O3 columns. We found that all sources of a priori O3 profiles evaluated in this study generally reproduced the vertical structure of summer-averaged observations. However, larger differences between the a priori profiles and lidar observations were observed when evaluating inter-daily and diurnal variability of tropospheric O3. The TB-Clim O3 profile climatology was unable to replicate observed inter-daily and diurnal variability of O3 while model products, in particular GEOS-Chem simulations, displayed more skill in reproducing these features. Due to the ability of models, primarily the CTM used in this study, on average to capture the inter-daily and diurnal variability of tropospheric and LMT O3 columns, using a priori profiles from CTM simulations resulted in TEMPO retrievals with the best statistical comparison with lidar observations. Furthermore, important from an air quality perspective, when high LMT O3 values were observed, using CTM a priori profiles resulted in TEMPO LMT O3 retrievals with the least bias. The application of time-specific (non-climatological) hourly/daily model predictions as the a priori profile in TEMPO O3 retrievals will be best suited when applying this data to study air quality or event-based processes as the standard retrieval algorithm will still need to use a climatology product. Follow-on studies to this work are currently being conducted to investigate the application of different CTM-predicted O3 climatology products in the standard TEMPO retrieval algorithm. Finally, similar methods to those used in this study can be easily applied by TEMPO data users to recalculate tropospheric O3 profiles provided from the standard retrieval using a different source of a priori.

Satellite↗
NASA NTRS1994

Effects of the intense solar activity of March/June 1991 observed in the outer heliosphere

The properties of the large-scale global merged interaction region (GMIR) generated by the intense solar events of March and June 1991 are studied using the available solar wind, interplanetary magnetic field, and energetic particle data from the observing network of Pioneer 10 and Voyagers 1 and 2 in the outer heliosphere. At heliocentric distances extending to 55 AU the delayed effects of this enhanced solar activity are observed in the form of large inceases in the solar wind velocity and interplanetary magnetic field and significant decreases in the galactic cosmic ray intensity. For low-energy ions (5-MeV protons) there was a single long-lived event extending over a period of some 6 months. Near the strongest interplanetary disturbances the H and He spectra are best represented by similar exponentials in momentum/nucleon (i.e., particle velocity at these at these energies). Over the rest of the event the characteristic momentum for He, (P(sub 0))(sub He) is generally approximately 0.66 for hydrogen. These spectra and the consistently low H/He ratio (25.3) at 2 MeV/nucleon closely resemble that observed in corrotating interaction regions events. Despite the strong north/south asymmetry in the solar activity, the interplanetary disturbances produced the same net decrease in the galactic cosmic ray intensity of ions greater than 70 MeV at the three widely separated spacecraft when the effects of the long-term recovery are taken into account. A comparison of the relative intensity of MeV ions at these three spacecraft suggest that the most intense solar events occurred on the back side of the Sun in time periods adjacent to the March and June episodes of solar activity. It is argued that this GMIR as a system is responsible for the low-frequency radio emission observed by the Voyager Plasma Wave experiment some 1.46 years after the onset of the March 1991 activity.

Mcdonald, F. B.↗
DOE OSTI2021

Modeling the Fundamental Radiation Chemistry of the Organic Diluent, and the Effect of Metal Ion Complexation on the Radiochemical Behavior of Active Compounds

All used nuclear fuel (UNF) reprocessing technologies must operate efficiently in the presence of an intense, multi-component (predominantly alpha, beta, and gamma) radiation field. Consequently, radiation-induced degradation of reprocessing systems is of concern, as it negatively impacts process performance over time due to the destruction of both active compounds (ligands, phase modifiers, holdback agents, etc.) and the formation of degradation products. Reprocessing solvent system radiolysis has been linked to changes in separation efficiency and physical properties of solvent mixtures, solvent-recycle longevity, crud formation, and other unexpected outcomes that impact the efficient recovery of valuable materials (e.g., the actinides) and the volume of hazardous radioactive waste for final disposal, i.e., in a geological repository. Consequently, a fundamental understanding of radiolytic processes and their effects on reprocessing solvent system performance is critical for: (i) the cost-effective development and innovation of separation technologies; (ii) the design and implementation of predictive radiation chemical models for process monitoring and lifetimes; and (iii) potentially the ability to exploit radiolytic phenomenon to our benefit, e.g., strategic radiolysis of active molecules to liberate specific degradation products that aid subsequent process stages. Despite extensive investigation into the radiolytic behavior of active solvent system compounds, little attention has been given to understanding (i) the radiation chemical behavior and modification of the organic diluent and (ii) the effect of metal ion complexation on the radiochemical behavior of active compounds.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

Using radiochemistry to diagnose fuel-ablator mix in inertial confinement fusion studies at the National Ignition Facility: measurement of Tc/Mo isotopic ratios for the Pushered Single Shell campaign

The development of radiochemical measurement techniques as a diagnostic for fusion experiments at the National Ignition Facility enables a new method for assessing fuel-ablator mix and the impact of this mix on capsule performance. Diagnosing capsule mix in internal confinement fusion studies is difficult due to the small spatial scales (10s of µm) and short-time frames (100s of ps) over which the mix typically evolves in these experiments. For the Pushered Single Shell campaign, radiochemical measurements on debris collected from fusion experiments can be used to determine isotopic ratios of activation products, particularly 96g Tc/ 99 Mo and 95g Tc/ 99 Mo, to provide vital information on nuclear reactions in the burning plasma that can inform simulations that seek to understand the degree of capsule-fuel mix and the impact on the capsule performance. These radiochemical measurements have been conducted regularly since November 2023 providing data on a range of capsule designs and neutron yields. Data from eight NIF experiments is presented, the measured 96g Tc/ 99 Mo and 95g Tc/ 99 Mo range from (0.5–5) × 10 –4 to (0.3–3) × 10 –4 , respectively. The development of radiochemical diagnostics aids in understanding and optimizing the design of fusion experiments, providing unique and valuable insights into capsule behavior and directly measuring fuel-ablator mix.

and nuclear chemistry↗
NASA NTRS2006

Generating Aromatics From CO2 on Mars or Natural Gas on Earth

Methane to aromatics on Mars ( METAMARS ) is the name of a process originally intended as a means of converting Martian atmospheric carbon dioxide to aromatic hydrocarbons and oxygen, which would be used as propellants for spacecraft to return to Earth. The process has been demonstrated on Earth on a laboratory scale. A truncated version of the process could be used on Earth to convert natural gas to aromatic hydrocarbon liquids. The greater (relative to natural gas) density of aromatic hydrocarbon liquids makes it more economically feasible to ship them to distant markets. Hence, this process makes it feasible to exploit some reserves of natural gas that, heretofore, have been considered as being "stranded" too far from markets to be of economic value. In the full version of METAMARS, carbon dioxide is frozen out of the atmosphere and fed to a Sabatier reactor along with hydrogen (which, on Mars, would have been brought from Earth). In the Sabatier reactor, these feedstocks are converted to methane and water. The water is condensed and electrolyzed to oxygen (which is liquefied) and hydrogen (which is recycled to the Sabatier reactor). The methane is sent to an aromatization reactor, wherein, over a molybdenum-on-zeolite catalyst at a temperature 700 C, it is partially converted into aromatic hydrocarbons (specifically, benzene, toluene, and naphthalene) along with hydrogen. The aromatics are collected by freezing, while unreacted methane and hydrogen are separated by a membrane. Most of the hydrogen is recycled to the Sabatier reactor, while the methane and a small portion of the hydrogen are recycled to the aromatization reactor. The partial recycle of hydrogen to the aromatization reactor greatly increases the catalyst lifetime and eases its regeneration by preventing the formation of graphitic carbon, which could damage the catalyst. (Moreover, if graphitic carbon were allowed to form, it would be necessary to use oxygen to remove it.) Because the aromatics contain only one hydrogen atom per carbon atom, METAMARS produces four times as much propellant from a given amount of hydrogen as does a related process that includes the Sabatier reaction and electrolysis but not aromatization. In the terrestrial version of METAMARS, the Sabatier reactor and electrolyzer would be omitted, while the hydrogen/ methane membrane-separating membrane, the aromatization reactor, and the unreacted-gas-recycling subsystem would be retained. Natural gas would be fed directly to the aromatization reactor. Because natural gas consists of higher hydrocarbons in addition to methane, the aromatization subprocess should be more efficient than it is for methane alone.

Muscatello, Anthony C.↗
NASA NTRS2023

Automated Fiber Placement Through Thickness Defect Stacking Optimization

In its 2022 commercial market outlook, Boeing forecasted an 80% increase in the global fleet through 2041 compared to 2019 pre-pandemic levels. This sharp rise in demand will drive pressure onto airframe manufacturers to ramp up production and find more efficient ways to design and manufacture airplanes. Complicating this challenge is the industry’s recent transformation from traditional metal-based airframes towards hybrid composite-metal aircraft. While composites have been used in aviation for decades, aircraft manufacturers are still struggling to design and manufacture quality parts at a high rate. Automated Fiber Placement (AFP) is one of the main manufacturing techniques used to produce large-scale composite parts. After a design has been created, a manufacturing strategy has to be developed based on the working material, part geometry, and machine capabilities. This process planning stage is essential to the AFP workflow and currently requires a high level of manual input from an experienced process planner. In an effort to automate and optimize this stage, the Computer Aided Process Planning (CAPP) module was developed. CAPP assists process planners in identifying optimal starting point location and layup strategy for each ply of a laminate. This Ply-Level Optimization (PLO) phase operates on the quantification of ply quality through predictable geometry-based defects such as gaps, overlaps, angle deviation, and steering. As you move from PLO to Laminate-Level Optimization (LLO) the design space grows exponentially, emphasizing the need for automated optimization. The work presented in this thesis expands CAPP’s functionality by comparing the planned fiber paths through the thickness of the laminate to mitigate stacked area defects and achieve an optimal laminate-level manufacturing strategy. Within CAPP, predicted gap and overlap defects are imported from Vericut Composites Programming (VCP) and then discretized to streamline the through-thickness comparison. Two objective functions are used to score different combinations of ply layup strategies based on defect stacking both globally and locally. Four combinatorial optimization algorithms were coupled with these objective functions to investigate the laminate-level manufacturing strategy design space and converge on the optimal plan. These algorithms were evaluated based on accuracy and efficiency through virtual testing on a complex tool surface. A separate LLO approach was developed to achieve near-optimal laminates in significantly less time. The end result is a software package which greatly reduces the required input from process planners, shortening the design-build cycle time and improving part quality.

AFP↗
NASA NTRS2002

Major Effects of Nonmetallic Inclusions on the Fatigue Life of Disk Superalloy Demonstrated

The fatigue properties of modern powder metallurgy disk alloys can vary because of the different steps of materials and component processing and machining. Among these variables, the effects of nonmetallic inclusions introduced during the powder atomization and handling processes have been shown to significantly degrade low-cycle fatigue life. The levels of inclusion contamination have, therefore, been reduced to less than 1 part per million in state-of-the-art nickel disk powder-processing facilities. Yet the large quantities of compressor and turbine disks weighing from 100 to over 1000 lb have enough total volume and surface area for these rare inclusions to still be present and limit fatigue life. The objective of this study was to investigate the effects on fatigue life of these inclusions, as part of the Crack Resistant Disk Materials task within the Ultra Safe Propulsion Project. Inclusions were carefully introduced at elevated levels in a nickel-base disk superalloy, U720, produced using powder metallurgy processing. Multiple strain-controlled fatigue tests were then performed on extracted test specimens at 650 C. Analyses were performed to compare the low-cycle fatigue lives and failure initiation sites as functions of inclusion content and fatigue conditions. Powder of the nickel-base superalloy U720 was atomized in argon at Special Metals Corporation, Inc., using production-scale high-cleanliness powder-processing facilities and handling practices. The powder was then passed through a 270-mesh screen. One portion of this powder was set aside for subsequent consolidation without introduced inclusions. Two other portions of this powder were seeded with alumina inclusions. Small, polycrystalline soft (Type 2) inclusions of about 50 mm diameter were carefully prepared and blended into one powder lot, and larger hard (Type 1) inclusions of about 150 mm mean diameter were introduced into the other seeded portion of powder. All three portions of powder were then sealed in separate containers, hot isostatically pressurized, extruded, forged into subscale disks, and heat treated. Low-cycle-fatigue specimens were then extracted, machined, and tested. Fatigue tests were performed at 650 C in closed-loop servohydraulic testing machines using induction heating and axial extensometers. All tests were continued to failure, and fractographic evaluations were performed on all specimens to determine the crack initiation sites. A large majority of the failures in specimens with introduced inclusions occurred at cracks initiating from inclusions at the specimen surface, as shown for each type of inclusion in the following bar chart. The inclusions significantly reduced fatigue life from unseeded material levels, as shown in the bar chart. These effects were found to depend on the strain range, strain ratio, and inclusion size. Tests at lower strain ranges and higher strain ratios resulted in larger effects of inclusions on life. Inclusion effects on life were thereby maximized in tests at the lowest strain range of 0.6 percent and the most positive strain ratio of 0.5. Under these conditions, small Type 2 inclusions reduced life substantially-- about 20 times, whereas large Type 1 inclusions dramatically reduced life 100 times. These results clearly demonstrate that it is essential to include the effects of inclusions for realistic predictions of disk fatigue life. Important issues, including temperature dependence, crack initiation versus propagation, surface treatments, realistic disk features and machining, and realistic disk spin testing will be addressed to accurately model inclusion effects on disk fatigue life. Fatigue life varied from well over 105 cycles for no inclusions to a little over 103 cycles for 100-micrometer inclusions. A single crack initiating at a surface-connected seeded inclusion caused failure in each case.

Gabb, Timothy P.↗
DOE OSTI2021

Specific absorbed fractions for a revised series of the UF/NCI pediatric reference phantoms: internal photon sources

Assessment of radiation absorbed dose to internal organs of the body from the intake of radionuclides, or in the medical setting through the injection of radiopharmaceuticals, is generally performed based upon reference biokinetic models or patient imaging data, respectively. Biokinetic models estimate the time course of activity localized to source organs. The time-integration of these organ activity profiles are then scaled by the radionuclide S-value, which defines the absorbed dose to a target tissue per nuclear transformation in various source tissues. S-values are computed using established nuclear decay information (particle energies and yields), and a parameter termed the specific absorbed fraction (SAF). The SAF is the ratio of the absorbed fraction—fraction of particle energy emitted in the source tissue that is deposited in the target tissue—and the target organ mass. While values of the SAF may be computed using patient-specific or individual-specific anatomic models, they have been more widely available through the use of computational reference phantoms. In this study, we report on an extensive series of photon SAFs computed in a revised series of the University of Florida and the National Cancer Institute pediatric reference phantoms which have been modified to conform to the specifications embodied in the ICRP reference adult phantoms of Publication 110 (e.g. organs modeled, organ ID numbers, blood contribution to elemental compositions). Following phantom anatomical revisions, photon radiation transport simulations were performed using MCNPX v2.7 in each of the ten phantoms of the series—male and female newborn, 1 year old, 5 year old, 10 year old, and 15 year old—for 60 different tissues serving as source and/or target regions. A total of 25 photon energies were considered from 10 keV to 10 MeV along a logarithm energy grid. Detailed analyses were conducted of the relative statistical errors in the Monte Carlo target tissue energy deposition tallies at low photon energies and over all energies for source–target combinations at large intra-organ separation distances. Based on these analyses, various data smoothing algorithms were employed, including multi-point weighted data smoothing, and log–log interpolation at low energies (1 keV and 5 keV) using limiting SAF values based upon target organ mass to bound the interpolation interval. The final dataset is provided in a series of ten electronic supplemental files in MS Excel format. Finally, the results of this study were further used as the basis for assessing the radiative component of internal electron source SAFs as described in our companion paper (Schwarz et al 2021) for this same pediatric phantom series.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2021

Specific absorbed fractions for a revised series of the UF/NCI pediatric reference phantoms: internal electron sources

Assessment of radiation absorbed dose to internal organs of the body from the intake of radionuclides, or in the medical setting through the injection of radiopharmaceuticals, is generally performed based upon reference biokinetic models or patient imaging data, respectively. Biokinetic models estimate the time course of activity localized to source organs. The time-integration of these organ activity profiles are then scaled by the radionuclide S-value, which defines the absorbed dose to a target tissue per nuclear transformation in various source tissues. S-values are computed using established nuclear decay information (particle energies and yields), and a parameter termed the specific absorbed fraction (SAF). The SAF is the ratio of the absorbed fraction—fraction of particle energy emitted in the source tissue that is deposited in the target tissue—and the target organ mass. While values of the SAF may be computed using patient-specific or individual-specific anatomic models, they have been more widely available through the use of computational reference phantoms. In this study, we report on an extensive series of photon SAFs computed in a revised series of the University of Florida and the National Cancer Institute pediatric reference phantoms which have been modified to conform to the specifications embodied in the ICRP reference adult phantoms of Publication 110 (e.g. organs modeled, organ ID numbers, blood contribution to elemental compositions). Following phantom anatomical revisions, photon radiation transport simulations were performed using MCNPX v2.7 in each of the ten phantoms of the series—male and female newborn, 1 year old, 5 year old, 10 year old, and 15 year old—for 60 different tissues serving as source and/or target regions. A total of 25 photon energies were considered from 10 keV to 10 MeV along a logarithm energy grid. Detailed analyses were conducted of the relative statistical errors in the Monte Carlo target tissue energy deposition tallies at low photon energies and over all energies for source–target combinations at large intra-organ separation distances. Based on these analyses, various data smoothing algorithms were employed, including multi-point weighted data smoothing, and log–log interpolation at low energies (1 keV and 5 keV) using limiting SAF values based upon target organ mass to bound the interpolation interval. The final dataset is provided in a series of ten electronic supplemental files in MS Excel format. Finally, the results of this study were further used as the basis for assessing the radiative component of internal electron source SAFs as described in our companion paper (Schwarz et al 2021) for this same pediatric phantom series.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Developing And Scaling an OpenFOAM Model to Study Turbulent Flow in a HFIR Coolant Channel

Improving the understanding of how computational fluid dynamics (CFD) direct numerical simulations (DNS) of flows in the High Flux Isotope Reactor (HFIR) perform when run in parallel using the high performance computing (HPC) platform Summit at the Oak Ridge Leadership Computing Facility (OLCF) is of particular importance to boost the computational tools used to support HFIR conversion to low enriched fuel (LEU). Evaluation of scaling performance was driven by the increasing importance of graphics processing unit (GPU) usage in HPC, which is becoming the standard for modern supercomputers such as Summit. The desired results are to obtain a strong positive correlation between the computational resources dedicated to a problem and the relative speed-up of the simulation in comparison to a benchmark. This capability will allow substantially improvement in HFIR flow analytical capabilities, specifically when predicting turbulence properties at high Reynolds numbers. The study leverages previous simulation results performed with code PHASTA (finite element) on HPC platforms Cori (NERSC) and Theta (ALCF) [1] with computing options provided in the computing platform OpenFOAM (finite volume) at OLCF. Transitioning from PHASTA to OpenFOAM will (1) eliminate dependence on third-party software for mesh generation and manipulation, (2) reduce resource needs by employing modern architectures, and (3) build expertise for future modeling of HFIR-specific problems like heat transfer in involute geometry, entrance effects, flow structure in channel corners, and so on—all important issues when defining the available thermal margins in the transition to LEU. CPUs and GPUs differ significantly in their architecture and utilization, as discussed in the literature [2]. The most important differences are in the approach to computations and their memory. A single GPU contains a large quantity of cores, enabling it to perform with a much higher throughput than a CPU, but execution requires a different approach. GPU codes execute instructions using the Single-Instruction Multiple-Thread (SIMT) approach in which a single instruction is used for groups of threads called warps. A warp typically consists of 32 threads which must execute the same set of instructions, although on separate threads. Alternately, a CPU has far fewer cores that are much more flexible in their operation, excelling at quickly performing more complex serial computations. This is why GPUs have greater throughput when properly utilized. The second important difference is seen when comparing their memory spaces. Limited memory allocations and CPU–GPU communications cause a significant bottleneck in GPU-accelerated programs. Further study was required to properly take advantage of GPU resources. A comprehensive analysis of code performance and the model-specific features of turbulence constitutes the core of this work. In this study, a DNS simulation of HFIR channel turbulence was performed with the finite volume CFD code OpenFOAM v2112 and CUDA v11.0 on Red Hat Enterprise Linux v8.2. The OpenFOAM installation had AMGx integrated to enable GPU acceleration and utilizes the PETSc4FOAM library. The computational resources and the problem size were scaled on CPU and CPU + GPU architectures to gain a better understanding of the performance of a DNS problem on modern computing hardware. The study aimed to analyze the scaling of the code exclusively on CPUs and then to examine the scaling of the codes with GPU acceleration enabled. Scaling studies included CPU and GPU acceleration on a mesh of varying resolution to analyze the impact of problem size relative to computational resources. In the course of preparing the GPU configuration on Summit, mainly using the AMGX solvers, difficulties were encountered stemming from constant changes resulting from extensive ongoing development activities and the changing environment. This resulted in the inability to complete the GPU portion of the work. The code was compiled and tested, but production runs to assess acceleration were not performed because the used discretional compute time allocation expired as year-end approached. The Summit HPC platform is scheduled for decommissioning in 2024, making it unattractive for future use with Nvidia-based GPUs. Therefore, the work will be moved onto NERSC machines in FY24. An application was prepared and submitted, and sufficient node-hours were awarded to continue the research in the next calendar year. This report summarizes work performed thus far, which mostly focused on CPU OpenFOAM computing.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS2016

The Viability of Trajectory Analysis for Diagnosing Dynamical and Chemical Influences on Ozone Concentrations in the UTLS

The viability of trajectory analysis for diagnosing the interplay between chemistry and dynamics is investigated by comparing ozone mixing ratios modelled using air-parcel pathways to values observed along flight tracks during ATTREX (Airborne Tropical TRopopause EXperiment). Trajectories are initiated at the locations of ozone observations and tracked backward in time to their sources at termini of backward trajectories. The modelled values of ozone utilize 3-dimensional analysis fields from WACCM (Whole Atmosphere Community Climate Model) (a chemical-climate model with dynamical fields nudged towards MERRA (Modern-Era Retrospective Analysis and Research Applications) reanalysis) and ERA-interim (product of ECMWF - the European Centre for Medium-Range Weather Forecasts) to determine source mixing ratios with chemical production and loss terms derived from the ozone chemistry used in WACCM. A statistical base of modelled ozone is constructed with 6 trajectory platforms (adiabatic, diabatic, and kinematic forced by ERA-interim and MERRA), two chemical models (WACCM chemistry and no chemistry), and 4 trajectory lengths (5, 10, 20, and 30 days). Linear regression is employed to separate systematic errors from random errors and to characterize the impact of source mixing ratios, path length, vertical motion, and chemistry on modelled ozone errors. Errors in the analysis ozone fields are large, if not dominant, contributors to model error. Random errors are particularly large for point-by-point comparisons, however averaging over 800 km (75 minutes) flight segments substantially reduces random error and exposes systematic errors. Of the two analysis ozone data sets, WACCM, which incorporates detailed chemistry, provides the smaller systematic errors while ERA-interim, which has crude chemistry but assimilates observational data, has the smaller random errors. Of the different trajectory platforms, adiabatic calculations produce the smaller random errors (irrespective of the use of chemistry) but both vertical motion and chemistry are required to optimally reduce systematic errors. These results suggest that meaningful analysis of dynamical and chemical interactions that control ozone mixing ratios are viable on spatial scales larger than a few reanalysis grid spaces, that errors in the analyzed ozone data sets are large but not prohibitively so, and that vertical velocities and heating rates from reanalysis data, while problematic, contain useful information [on the ozone concentrations in the UTLS (Upper Troposphere/Lower Stratosphere)].

ozone↗
NASA NTRSDate not supplied

Accuracy of Center of Pressure Determination via Motion Capture

BACKGROUND: This study was conducted to support the stability assessment for tasks in lunar gravity and exercises on a Vibration Isolation and Stabilization (VIS) system in microgravity based on the dynamic feasibility criterion of whether the calculated position of the center of pressure (COP) falls within the base of support (BOS) which outlines the subject’s feet. Motion capture data combined with biomechanical modeling and simulation allows the forces and moments between the human and the VIS platform to be computed and the position of the COP as well as the location and shape of the BOS to be determined. The goal of this study was to assess the accuracy of the COP trajectory calculated using motion capture-based data. METHODS AND RESULTS: To obtain the dynamic quantities from which COP is calculated, motion capture data is first collected in the 1g lab environment by recording the trajectories of passive retroreflective markers placed on a subject during exercise or performance of a given task. The OpenSim [1] inverse kinematics (IK) tool is used to fit a scaled subject model to recorded marker trajectories while minimizing marker error to obtain joint angles. Then, a custom OpenSim plugin [2] is used to determine the subject’s time-varying moment of inertia and its time derivative, center of mass (CM) position, velocity, and acceleration, as well as the angular momentum and its time derivative relative to the subject’s CM. Some of these quantities are not needed for modeling tasks performed on a stationary lunar surface but, due to the moving exercise platform, are needed to model VIS response to the subject’s motion. Hand positions, used in calculating a cable force if present, are recorded as well. These quantities are used to calculate the total force (F ⃗^((plate) )) and moment (M ⃗^((plate) )) exerted by the lunar surface or the VIS plate on the subject’s shoe soles. COP is then calculated from the following equations: r_x^((cop) )= M_z^((plate) )/F_y^((plate) ) and r_z^((cop) )= 〖-M〗_x^((plate) )/F_y^((plate) ), where the y axis is normal to the surface. COP accuracy for feasibility assessments is then determined by whether it falls within the BOS, which is also computed by the plugin. To study the accuracy of COP calculated from motion capture, we first investigated whether COP remained within the BOS, as it must, for exercises performed in the 1g lab environment. Standard exercises such as back squat and deadlift were analyzed, as well as more explosive exercises including hang clean and press. Cases in which the COP exited the BOS indicated that COP accuracy required further investigation. In this study, an exercise device with cables was used, so cable force modeling accuracy should also be considered. In a separate study, we collected motion capture and force plate data for twenty-seven motions not involving an exercise device. About a third were genuine countermeasures exercises (e.g., hang clean and press), some were relevant for lunar tasks (e.g., object pick up), and the rest were of a “unit test” nature (e.g., swaying back and forth or side to side). Motion capture-based COP positions were compared with force plate measured COP. We found that while force plate measured COP remained within the BOS, motion capture-based COP was observed to briefly exit the BOS on occasion. Techniques to mitigate IK artifacts and filtering of calculated data could be used to improve the agreement of calculated and measured results, resolving excursions from the BOS within this dataset. The mean error between calculated and measured COP was found to be less than 6 mm. Additionally, we derived and investigated equations for the COP in terms of the cable force, cable location, as well as the human CM position, acceleration, and angular momentum with respect to the CM, and analyzed them for sensitivity to errors in individual quantities. Several were found, but the most significant one was that when the vertical force on the feet approaches zero, indicating a near-detachment or ‘jump off’ condition, errors are amplified. This is consistent with the observation that in the absence of pressure, the concept of the center of pressure would become meaningless.

C A Bell↗
DOE OSTI2022

NRAP-Open-IAM: NRAP Open Source Integrated Assessment Model

Note: This is the last version (a2.6.1) of NRAP-Open-IAM released during NRAP Phase II in 2022. The latest version of NRAP-Open-IAM is available here: https://edx.netl.doe.gov/dataset/phase-iii-nrap-open-iam NRAP-Open-IAM is an open-source software product that enables quantification of containment effectiveness and leakage risk at storage sites in the context of system uncertainties and variability. NRAP-Open-IAM represents the next-generation in a line of systems-based computational models developed for quantitative geological carbon storage (GCS) risk assessment. The model comprises a set of reduced-order and analytical models of various components of the GCS system, potential leakage pathways, receptors of concern including impact to groundwater resources and the atmosphere, a framework to support stochastic simulation, time stepping, uncertainty quantification, other analytical functionality for scenario and risk-performance evaluation, and a basic graphical user interface to support scenario development, data input simulation definition, and basic post-processing and results display. As the NRAP Open-IAM functionality continues to evolve, we continue to add to its capability to develop quantitative, probabilistic, and time-dependent profiles of the evolution of risk at a GCS site and evaluate the influence of uncertain parameters on uncertainty in predicted risk. It can be used to quantify the dynamics of reservoir saturation plume and pressure-affected area, for evaluation of the area of potential groundwater impact (i.e., Area of Review) and monitoring requirements to support cost and regulatory analysis, and for consideration of different post-injection site care and closure scenarios. This submission contains the current version of NRAP-Open-IAM available for evaluation and testing. To use the NRAP-Open-IAM, download the source code (https://edx.netl.doe.gov/dataset/08f396c8-bc5f-44ad-a028-6e98a6ea6d70/resource/4c24a3da-3b40-4ffe-9892-c807ae9f8760) then open the NRAP-Open-IAM user's guide (https://edx.netl.doe.gov/dataset/08f396c8-bc5f-44ad-a028-6e98a6ea6d70/resource/8b27335a-343c-4836-b8a3-3ad0bdc9e669) to read more about the tool. Installation instructions for Windows, Mac, and Linux can be found in the "installers" folder of the extracted NRAP-Open-IAM folder and describe setup of environment (e.g., Python libraries) needed for proper work of the tool. Test of installation can be done by running "python openiam_setup_tests.py" in the "setup" folder. The installation test also runs a test suite to see if the NRAP-Open-IAM has been installed correctly. To run the test suite separately, run "python iam_test.py" in the "test" folder. User's guide: https://edx.netl.doe.gov/dataset/08f396c8-bc5f-44ad-a028-6e98a6ea6d70/resource/8b27335a-343c-4836-b8a3-3ad0bdc9e669 Developer's guide: https://edx.netl.doe.gov/dataset/08f396c8-bc5f-44ad-a028-6e98a6ea6d70/resource/3bc6ee7d-609d-4eb6-80ba-fa6130ee0313 Reservoir simulation data used in some examples distributed with NRAP-Open-IAM: - Kimberlina: https://edx.netl.doe.gov/dataset/08f396c8-bc5f-44ad-a028-6e98a6ea6d70/resource/eb62cece-61b2-4037-9b6d-32407dde2ab8 - Kimberlina (compartmentalized): https://edx.netl.doe.gov/dataset/08f396c8-bc5f-44ad-a028-6e98a6ea6d70/resource/366f9530-3b32-4b84-affe-ab2df1d9a8b5 - FutureGen 2.0: https://edx.netl.doe.gov/dataset/futuregen-2-0-1008-simulation-reservoir-lookup-table NRAP-Open-IAM GitLab repository: https://gitlab.com/NRAP/OpenIAM Related publications: - Bacon, D., Yonkofski, C., Brown, C., Demirkanli, D. and Whiting, J., 2019. Risk-based post injection site care and monitoring for commercial-scale carbon storage: Reevaluation of the FutureGen 2.0 site using NRAP-Open-IAM and DREAM. International Journal of Greenhouse Gas Control 90: 102784. - Bacon, D. Demirkanli, D., and White, S., 2020. Probabilistic risk-based Area of Review (AoR) determination for a deep-saline carbon storage site. International Journal of Greenhouse Gas Control 102: 103153. - Harp, D., Oldenburg, C., and Pawar, R., 2019. A metric for evaluating conformance robustness during geologic CO2 sequestration operations. International Journal of Greenhouse Gas Control 85: 100-108. - Lackey, G., Vasylkivska, V., Huerta, N., King, S., and Dilmore, R., 2019. Managing well leakage risks at a geologic carbon storage site with many wells, International Journal of Greenhouse Gas Control, 88 :182-194. - Vasylkivska, V., Dilmore, R., Lackey, G., Zhang, Y., King, S., Bacon, D., Chen, B., Mansoor, K., and Harp, D., 2021. NRAP-Open-IAM: A flexible open-source integrated assessment model for geologic carbon storage risk assessment and management, Environmental Modelling & Software, 143: 105114. Presentations: - Chen, B., Harp, D., and Pawar, R., A data assimilation approach (ES-MDA) coupling with NRAP-Open-IAM for quantifying uncertainty reduction in geological CO2 sequestration. AGUFM 2019: T44A-02. - Chen, B., and Harp, D., Improving risk analysis precision for geologic CO2 sequestration by quantifying the uncertainty reduction before and after acquiring monitoring data. 14th Greenhouse Gas Control Technologies Conference, Melbourne, Australia, 2018, pp. 21-26. - Harp, D., National Risk Assessment Partnership Task 2: Containment Assurance. No. LA-UR-19-28654, Los Alamos National Laboratory (LANL), Los Alamos, NM (United States), 2019. - Vasylkivska, V., King, S., Bacon, D., Harp, D., Chen, B., Mansoor, K., Onishi, T., Yang, Y., Zhang, Y., and Keating, E., NRAP-Open-IAM: An open-source integrated assessment model, poster, Mastering the Subsurface Through Technology Innovation, Partnerships and Collaboration: Carbon Storage and Oil and Natural Gas Technologies Review Meeting, Pittsburgh, PA, August 13-16, 2018. - Vasylkivska, V., Lackey, G., King, S., Wentworth, A., Huerta, N., Creason, C., DiGiulio, J., Yang, Y., and Dilmore, R., Long-term risk analysis of a geologic CO2 storage project during the post-injection period, SIAM Conference on Computational Science and Engineering, Spokane, WA, February 25-March 1, 2019. - Vasylkivska, V., Overview of the NRAP-Open-IAM tool for carbon storage (beta release), 2019 Annual NRAP Tool Users Meeting, Pittsburgh, PA, August 27, 2019. - Vasylkivska, V., Bacon, D., Chen, B., Dilmore, R., Harp, D., King, S., Lackey, G., Lindner, E., Liu, G., Mansoor, K. and Zhang, Y., NRAP-Open-IAM: A new, open-source code for integrated assessment of geologic carbon storage containment effectiveness and leakage risk, poster, American Geophysical Union Fall Meeting 2020 (virtual meeting), December 2020. - Vasylkivska, V., NRAP open-source integrated assessment model and relevant application, oral presentation, NRAP workshop "NRAP Tools for Geologic Carbon Storage Risk-Based Decision Making" held in conjunction with Groundwater Protection Council (GWPC) 2021 Annual Forum (virtual meeting), Salt Lake City, UT, September 2021. - Vasylkivska, V., NRAP-Open-IAM: open-source integrated assessment model, digital poster/demonstration, software demonstration session, 2022 Carbon Management Project Review Meeting, August 16, 2022

AoR↗
NASA NTRS2006

Low- to Mid-Latitude X-Ray Emission from Jupiter

The Chandra X-ray Observatory (CXO) observed Jupiter during the period 2003 February 24-26 for approx.40 hours (4 Jupiter rotations), using both the spectroscopy array of the Advanced CCD Imaging Spectrometer (ACIS-S) and the imaging array of the High-Resolution Camera (HRC-I). Two ACIS-S exposures, each approx.8.5 hr long, were separated by an HRC-I exposure of approx.20 hr. The low- to mid-latitude non-auroral disk X-ray emission is much more spatially uniform than the auroral emission. However, the low- to mid-latitude X-ray count rate shows a small but statistically significant hour angle dependence, and is higher in regions of relatively low surface magnetic field strength, confirming ROSAT results. In addition, the spectrum from the low surface field region shows an enhancement in the energy band 1.14- 1.38 keV, perhaps partly due to line emission from that region. Correlation of surface magnetic field strength with count rate is not found for the 2000 December HRC-I data, at a time when solar activity was high. The low- to mid-latitude disk X-ray count rate observed by the HRC-I in the 2003 February observation is about 50% of that observed in 2000 December, roughly consistent with a decrease in the solar activity index (F10.7 cm flux) by a similar amount over the same time period. The low- to mid-latitude X-ray emission does not show any oscillations similar to the -45 minute oscillations sometimes seen from the northern auroral zone. The temporal variation in Jupiter's non-auroral X-ray emission exhibits similarities to variations in solar X-ray flux observed by GOES and TIMED/SEE. The two ACIS-S 0.3-2 keV low- to mid-latitude X-ray spectra are harder than the auroral spectrum, and are different from each other at energies above 0.7 keV, showing variability in Jupiter s non-auroral X-ray emission on a time scale of a day. The 0.3-2.0 keV X-ray power emitted at low- to mid-latitudes is 0.21 GW and 0.39 GW for the first and second ACIS-S exposures, respectively. We suggest that X-ray emission from Jupiter's disk may be largely generated by solar X-rays resonantly and fluorescently scattered in its upper atmosphere, especially at times of high incident solar X-ray flux. However, the correlation of higher count rate with low surface magnetic-field strength indicates the presence of some secondary component, possibly ion precipitation from radiation belts closer to the planet than elsewhere at low- to mid-latitudes.

Bhardwaj, Anil↗