SEARCH · Search NASA

Results for “Divided differences”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,242 records · Page 69

All sources NASA NTRS DOE OSTI

DOE OSTI2017

Preliminary Evaluation of Removing Used Nuclear Fuel from Shutdown Sites

A preliminary evaluation of removing spent nuclear fuel (SNF) from 14 shutdown nuclear power plant sites was performed. At these shutdown sites the nuclear power reactors have been permanently shut down and the sites have been decommissioned or are undergoing decommissioning. The shutdown sites were Maine Yankee, Yankee Rowe, Connecticut Yankee, Humboldt Bay, Big Rock Point, Rancho Seco, Trojan, La Crosse, Zion, Crystal River, Kewaunee, San Onofre, Vermont Yankee, and Fort Calhoun. The evaluation was divided into four components: (1) characterization of the SNF and greater-than-Class C low-level radioactive waste (GTCC waste) inventory, (2) a description of the on-site infrastructure at the shutdown sites, (3) an evaluation of the near-site transportation infrastructure and transportation experience at the shutdown sites, and (4) an evaluation of the actions necessary to prepare for and remove SNF and GTCC waste. The primary sources for the inventory of SNF and GTCC waste were the U.S. Department of Energy (DOE) spent nuclear fuel inventory database, industry publications such as StoreFUEL, and government sources such as the U.S. Nuclear Regulatory Commission. The primary sources for information on the conditions of on-site infrastructure and near-site transportation infrastructure and experience included information collected during site visits, information provided by managers at the shutdown sites, Facility Interface Data Sheets compiled for DOE in 2005, Services Planning Documents prepared for DOE in 1993 and 1994, industry publications such as Radwaste Solutions, and Google Earth. State representatives, State Regional Group representatives, Tribal representatives, and a Federal Railroad Administration representatives have participated in eleven of the shutdown site visits. Every shutdown site was found to have at least one off-site transportation mode option for removing its SNF and GTCC waste; some have multiple options. Experience removing large components during reactor decommissioning provided an important source of information used to identify the transportation mode options for the sites. Especially important in conducting the evaluation were site visits, through which information was obtained that would not have been available otherwise. Extensive photographs taken during the site visits proved to be particularly useful in documenting the current conditions at or near the sites. Additional conclusions from this evaluation include: The 14 shutdown sites use designs from 4 different suppliers involving 11 different (horizontal and vertical) dry storage systems that would require the use of 9 different transportation cask designs to remove the SNF and GTCC waste from the shutdown sites. Although some changes to transportation certificates of compliance will be required, the SNF at the initial 9 shutdown sites (Maine Yankee, Yankee Rowe, Connecticut Yankee, Humboldt Bay, Big Rock Point, Rancho Seco, Trojan, La Crosse, and Zion) is in dual purpose dry storage canisters that can be transported, including a small amount of high-burnup fuel. Most sites indicated that 2-3 years of advance time would be required for its preparations before shipments could begin. Some sites could be ready in less time. As additional sites such as Palisades, Pilgrim, Three Mile Island, Oyster Creek, Indian Point, and Diablo Canyon shut down, these sites will be included in updates to the evaluation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS2018

An Evolutionary Computation System Design Concept for Developing Controlled Closed Ecosystems: An Intelligent Systems Approach to Foster Gravitational Ecosystem Research for Developing Sustainable Communities in Space and on Earth

An adjustably-autonomous intelligent systems approach for developing Closed Ecosystems (CESs) is presented, which includes a design concept and preliminary design details for the Controlled Closed-Ecosystem Development System (CCEDS) and the Orbiting Modular Artificial-Gravity Spacecraft (OMAGS). The paper is divided into three sections: CESs, the CCEDS Design Concept, and Orbiting Fractional-Gravity Closed Ecosystems OMAGS design concept. The first section briefly describes Closed EcoSystems (CESs), complex adaptive systems, biomes, microbial microbiomes, and their relevance for the study of astrobiology. This section also discusses initial efforts in the development of Closed Environment Life Support Systems (CELSSs) for sustainable communities in space and on Earth. This section concludes with a discussion of the bioregenerative life support system challenge of and the corresponding consequences due to the inverse relationship of the very small human biomass/non-human biomass ratio overall on the Earth with respect to the extremely large human biomass/non-human-biomass ratio found in cities and the International Space Station. The second section describes the CCEDS design concept, which consists of a population of controlled colonies of CES Modules (CESMs), each an integrated CES, continually generating data for an intelligent system that operates the CESs and their CESMs. A variety of CESM types and their use are briefly described. The CCEDS intelligent system uses an evolutionary computation algorithm described in this section to develop and optimize these CESs to increase their viability duration and the size of the animals they support with the ultimate goal to support populations of humans, both on Earth and in space. The CCEDS architecture, its five control subsystems, and its five evolutionary computation levels are also discussed. The section concludes with a discussion of several CCEDS design strategies. The third section summarizes the OMAGS design concept for a spacecraft with a payload consisting of CESs in an orbiting spacecraft centrifuge that operates for at least 5 years. The spacecraft concept is described including its 150cm-radius centrifuge with a 2 ton & 3,000 liter bioscience payload capacity for 24 CESMs. The centrifuge design has four physical levels for its CESMs, each level subject to a different fractional gravity level. This section presents the spacecraft benefits of being designed and operated such that the spacecraft and payload centrifuge wheel counter-rotate resulting in net zero angular momentum and zero gyroscopic forces. Artificial-gravity generation by centripetal acceleration is also discussed. This section concludes by showing the external specifications of the CESMs and their layout in the centrifuge, followed by discussing the multi-payload module rationale. In tandem, the CCEDS and OMAGS systems can be used to foster gravitational ecosystem research for developing sustainable communities in space and on Earth.

Dorais, Gregory A.↗

NASA NTRS2014

High-Energy Pulsar Models: Developments and New Questions

The past few years have seen a major advance in observational knowledge of high-energy (HE) pulsars. The Fermi Large Area Telescope (LAT) and AGILE have increased the number of known gamma-ray pulsars by an order of magnitude, its members being divided roughly equally among millisecond pulsars (MSPs), young radio-loud pulsars, and young radio-quiet pulsars. Many new and diverse emission characteristics are being measured, while radio and X-ray follow-up observations increase the pulsar detection rate and enrich our multiwavelength picture of these extreme sources. The wealth of new data has provided impetus for further development and improvement of existing theoretical pulsar models. Geometric light curve (LC) modelling has uncovered three broad classes into which HE pulsars fall: those where the radio profile leads, is aligned with, or lags the gamma-ray profile. For example, the original MSP and original black widow system are members of the second class, requiring co-located emission regions and thereby breaking with traditional notions of radio emission origin. These models imply narrow accelerator gaps in the outer magnetosphere, indicating copious pair production even in MSP magnetospheres that were previously thought to be pair-starved. The increased quality and variety of the LCs necessitate construction of ever more sophisticated models. We will review progress in global magnetosphere solutions which specify a finite conductivity on field lines above the stellar surface, filling the gap between the standard vacuum and force-free (FF; plasma-filled) models. The possibility of deriving phase-resolved spectra for the brightest pulsars, coupled with the fact that the HE pulsar population is sizable enough to allow sampling of various pulsar geometries, will enable much more stringent testing of future radiation models. Reproduction of the observed phase-resolved behavior of this disparate group will be one of the next frontiers in pulsar science, impacting on our understanding of particle acceleration, emission, and magnetosphere geometry. One may now also study evolutionary trends of the measured or inferred quantities, and probe pulsar visibility and population properties such as radiation beam sizes of different pulsar classes, as well as the distribution of spin-down power, gamma-ray luminosity, conversion efficiency, spectral index, and cutoff energy across the population. Lastly, the recent detection of very-high-energy (VHE) pulsations from the Crab pulsar generated quite a few ideas to explain this emission, leading to an extension of standard models and possibly even a bridge between the physics of pulsars and pulsar wind nebulae (PWNe).

gamma rays:↗

NASA NTRS1994

Catalogue of UBVRI photometry of T Tauri stars and analysis of the causes of their variability

A computer-based catalogue of UBVRI photoelectric photometry of T Tauri stars and their earlier type analogs has been compiled. It presently includes over 10 000 entries on 80 stars and will be updated on a regular basis; it is available on Internet. The catalogue is used to analyze the sometimes bizarre light variations of pre-main-sequence stars on time scales of days to months in an attempt to illuminate the nature and causes of the phenomenon. It is useful in discussing their light variations to divide the stars into three groups according to their spectra. These are: weak T Tauri stars (WTTS; spectral class later than K0 and W(sub H-alpha less than 10 A), classical T Tauri stars (CTTS; spectral class later than K0 and W(sub H-alpha) greater than 10 A), and early type T Tauri stars (ETTS; spectral class of K0 or earlier). Three distinct types of variability are displayed by stars in the catalogue. Type I variations are periodic in VRI and undoubtedly caused by rotational modulation of a star with an asymmetric distribution of cool spots on its surface. Irregular flare activity is sometimes seen on such stars in U and B. Type I variations are easiest to see on WTTS but are clearly present on CTTS and ETTS as well. Type II variations are caused by hot 'spots' or zones and, it is argued, result from changes in the excess or 'veiling' continuum commonly attributed to an accretion boundary layer or impact zone of a magnetically channeled accretion flow. This type of variation is seen predominantly or solely in CTTS. A subcategory, designated Type IIp, consists of stars which display periodic variations caused by hot spots. Whereas cool spots may last for hundreds or thousands of rotations, hot spots appear to come and go on a much shorter time scale. This suggests that both unsteady accretion and rotation of the star contribute to Type II variations. It is shown that a third type of variation exists among ETTS, including stars as early as A type. UX Ori is a typical example and we call these Type III variables or UXors. Their distinguishing characteristic is that they can display very large amplitudes (exceeding 2.8 mag in V) while showing little or no evidence for a veiling continuum or any substantial change in their photospheric spectra. If Type III variations are caused by changes in accretion luminosity, then boundary layers or impact zones in ETTS must be much different from CTTS which, of course, is possible since mass accretion rates are probably much higher. However, the leading hypothesis for explaining Type III variations is variable obscuration by circumstellar dust. It is argued that the putative dust clumps causing such variations cannot be confined to a disk; otherwise UXors would be rare. Perhaps magnetic effects are involved in levitating accreting dust out of the plane, as has been suggested for CTTS, or perhaps we are witnessing continuing infall of clumps from placental clouds. A third possibility is that dust may be condensing in an outflow.

Herbst, William↗

DOE OSTI2024

DECOVALEX-2023: Task E Final Report

This is the Task E final report for DECOVALEX-2023. Task E is focused on understanding thermal, two-phase hydrological, and mechanical (TH2M) processes, especially related to predicting brine migration in the excavation damaged zone around a heated excavation in salt. Salt is attractive as a disposal medium for radioactive waste because it is self-healing and is essentially impermeable and essentially non-porous in the far field (away from excavations). Investigation of the short-term (days to years) near-field (centimeters to tens of meters) behavior of salt is important for radioactive waste disposal because this early period strongly controls the amount of brine in a salt repository. Brine leads to corrosion of waste forms and waste packages, and possible dissolution of radionuclides with brine transport being a potential transport vector to the accessible environment. The main test case used in Task E is the ongoing Brine Availability Test in Salt (BATS) heater test located underground at the Waste Isolation Pilot Plant (WIPP) near Carlsbad, New Mexico, USA. The Task was divided into a series of Steps. Step 0 was an introduction to processes in salt, that included matching historical unheated brine inflow data from boreholes at WIPP and matching temperature observations during BATS heater test 1a. Step 1 included validation of models against a thermo-poroelastic analytical solution relevant to heated boreholes in salt, and two-phase flow around an excavation in salt. Step 2 required all the individual components covered in steps 0 and 1 to come together to match observed brine inflow behavior during the BATS 1a heater test. There were a range of approaches from the teams, from mechanistic to prescriptive. Given the uncertainties in the problem, some teams used one- or two-dimensional models of the processes, while other teams included more geometrical complexity in three-dimensional models. The key learning points from Task E have been: • Heat conduction through salt typically requires non-linear thermal conductivity (as a function of temperature), but most models do a good job matching observations, given appropriate adjustments to the applied power and some thermocouple locations. • Thermal pressurization requires coupled thermal-hydrological-mechanical (THM) responses that consider the thermal expansion of the fluid and solid phases. • Initialization of two-phase flow models around a borehole or excavation in salt are more realistically represented as “wetting up”, rather than “drying down” (i.e., the initial state after excavation is mostly dry, rather than mostly wet). • The BATS 1a heater test includes a significant release of brine after the end of heating, which requires a large increase in permeability to recreate. Task E has been a great learning experience for all the teams involved, and feedback from the modeling teams has led to changes in the design of follow-on BATS experiments, which are now ongoing underground at WIPP. There was a balance throughout the task between freedom to model phenomena how each team saw fit, and prescriptiveness in problem design to bring the modeling teams closer together to allow attribution of smaller differences between models to different modeling choices. The modeling approaches seem to go through two phases: an early phase of discovery or testing, and a later phase of refinement and improvement. In future modeling efforts, different field data could be used (e.g., BATS 2) and more time should be included in the processes for teams to make multiple model refinement or even significant changes to their conceptual model or setup, based on lessons learned from the modeling exercise.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2023

DEEPEN 3D PFA Favorability Models and 2D Favorability Maps at Newberry Volcano

DEEPEN stands for DE-risking Exploration of geothermal Plays in magmatic ENvironments. Part of the DEEPEN project involved developing and testing a methodology for a 3D play fairway analysis (PFA) for multiple play types (conventional hydrothermal, superhot EGS, and supercritical). This was tested using new and existing geoscientific exploration datasets at Newberry Volcano. This GDR submission includes images, data, and models related to the 3D favorability and uncertainty models and the 2D favorability and uncertainty maps. The DEEPEN PFA Methodology is based on the method proposed by Poux et al. (2020), which uses the Leapfrog Geothermal software with the Edge extension to conduct PFA in 3D. This method uses all available data to build a 3D geodata model which can be broken down into smaller blocks and analyzed with advanced geostatistical methods. Each data set is imported into a 3D model in Leapfrog and divided into smaller blocks. Conditional queries can then be used to assign each block an index value which conditionally ranks each block's favorability, from 0-5 with 5 being most favorable, for each model (e.g., lithologic, seismic, magnetic, structural). The values between 0-5 assigned to each block are referred to as index values. The final step of the process is to combine all the index models to create a favorability index. This involves multiplying each index model by a given weight and then summing the resulting values. The DEEPEN PFA Methodology follows this approach, but split up by the specific geologic components of each play type. These components are defined as follows for each magmatic play type: 1. Conventional hydrothermal plays in magmatic environments: Heat, fluid, and permeability 2. Superhot EGS plays: Heat, thermal insulation, and producibility (the ability to create and sustain fractures suitable for and EGS reservoir) 3. Supercritical plays: Heat, supercritical fluid, pressure seal, and producibility (the proper permeability and pressure conditions to allow production of supercritical fluid) More information on these components and their development can be found in Kolker et al., 2022. For the purposes of subsurface imaging, it is easier to detect a permeable fluid-filled reservoir than it is to detect separate fluid and permeability components. Therefore, in this analysis, we combine fluid and permeability for conventional hydrothermal plays, and supercritical fluid and producibility for supercritical plays. More information on this process is described in the following sections. We also project the 3D favorability volumes onto 2D surfaces for simplified joint interpretation, and we incorporate an uncertainty component. Uncertainty was modeled using the best approach for the dataset in question, for the datasets where we had enough information to do so. Identifying which subsurface parameters are the least resolved can help qualify current PFA results and focus future efforts in data collection. Where possible, the resulting uncertainty models/indices were weighted using the same weights applied to the respective datasets, and summed, following the PFA methodology above, but for uncertainty. There are two different versions of the Leapfrog model and associated favorability models: - v1.0: The first release in June 2023 - v2.1: The second release, with improvements made to the earthquake catalog (included additional identified events, removed duplicate events), to the temperature model (fixed a deep BHT), and to the index models (updated the seismicity-heat source index models for supercritical and EGS, and the resistivity-insulation index models for all three play types). Also uses the jet color map rather than the magma color map for improved interpretability. - v2.1.1: Updated to include v2.0 uncertainty results (see below for uncertainty model versions) There are two different versions of the associated uncertainty models: - v1.0: The first release in June 2023 - v2.0: The se...

15 GEOTHERMAL ENERGY↗

NASA NTRS2011

Electrically-Active Convection and Tropical Cyclogenesis in the Atlantic and East Pacific

It has been hypothesized that deep, intense convective-scale hot towers may aid the process of tropical cyclogenesis and intensification through dynamic and thermodynamic feedbacks on the larger meso-to-synoptic scale circulation. In this study, we make use of NCEP Reanalysis data and Tropical Rainfall Measurement Mission (TRMM) lightning imaging sensor (LIS), precipitation radar (PR), and microwave imager (TMI) data to investigate the role that widespread and/or intense lightning-producing convection (i.e., electrically-hot towers) present in African easterly waves (AEWs) may play in tropical cyclogenesis over the Atlantic, Caribbean, and East Pacific regions. NCEP Reanalysis 700 hPa meridional winds for the months of June to November for the years 2001-2009 were analyzed for the domain of 5 deg. N-20 deg. N and 130 deg. W-20 deg. E in order to partition individual AEWs into northerly, southerly, trough, and ridge phases. Subsequently, information from National Hurricane Center (NHC) storm reports was used to divide the waves into developing and non-developing waves. In addition, information from the NHC reports was used to further divide the developing waves into those waves that spawned storms that only developed to tropical storm strength and those that spawned storms that reached hurricane strength. The developing waves were also divided by the region in which they developed. To assess the evolution of convection associated with the AEWs as they propagated across our analysis domain, the full 130 deg. W-20 deg. E domain was divided into five longitude bands, and waves were analyzed for each band. To help determine the gross nature of the smaller convective scale, composites were created of all developing and non-developing waves as a function of AEW wave phase over the full analysis domain and each longitude band by compositing TRMM PR, TMI, LIS, and IR brightness temperature data extracted from the NASA global-merged infrared brightness temperature dataset. Finally, similar composites were created using various NCEP variables to assess the nature of the larger scale environment and circulation. Results suggest a clear distinction between developing and non-developing waves as developing waves near their region of development in terms of the intensity of convection (indicated by lightning flash rate), coverage of cold cloudiness (indicated by the percentage of a 2.5 deg. by 2.5 deg. box covered by IR brightness temperatures less than 240 deg. K and 210 deg. K), and large-scale variables, such as midlevel moisture and upper-level upward motion. For example, waves that developed in the East Pacific longitude band (i.e., 130 deg. W 95 deg. W) were observed in that band to have a flash rate of 56.4 flashes per day, a coverage by brightness temperatures less than 240 deg. K equal to 15.9%, a coverage by brightness temperatures less than 210 deg. K equal to 2.2%, a 700-hPa specific humidity anomaly of 0.4 g per kilogram, and a 300-hPa omega value of -0.04 Pascals per second in the trough phase compared to the non-developing wave trough values of 22.1 flashes per day, a coverage by brightness temperatures less than 240 deg. K equal to 8.1%, a coverage by brightness temperatures less than 210 deg. K equal to 0.9%, a 700-hPa specific humidity anomaly of -0.3 g per kilogram, and a 300-hPa omega value of -0.01 Pascals per second. Further analysis is being conducted to determine if the aforementioned behavior is observed for developing waves farther from their region of development and to determine any significant differences between waves that spawned storms that reached tropical storm strength and those that spawned storms that reached hurricane strength.

Leppert, Kenneth D., II↗

NASA NTRS2012

NASA Advanced Explorations Systems: Concepts for Logistics to Living

The NASA Advanced Exploration Systems (AES) Logistics Reduction and Repurposing (LRR) project strives to enable a largely mission-independent cradle-to-grave-to-cradle approach to minimize logistics contributions to total mission architecture mass. The goals are to engineer logistics materials, common crew consumables, and container configurations to meet the following five basic goals: 1. Minimize intrinsic logistics mass and improve ground logistics flexibility. 2. Allow logistics components to be directly repurposed for on-orbit non-logistics functions (e.g., crew cabin outfitting) thereby indirectly reducing mass/volume. 3. Compact and process logistics that have not been directly repurposed to generate useful on-orbit components and/or compounds (e.g., radiation shielding, propellant, other usable chemical constituents). 4. Enable long-term stable storage and disposal of logistics end products that cannot be reused or repurposed (e.g., compaction for volume reduction, odor control, and maintenance of crew cabin hygienic conditions). 5. Allow vehicles in different mission phases to share logistics resources. This paper addresses the work being done to meet the second goal, the direct repurposing of logistics components to meet other on-orbit needs, through a strategy termed Logistics to Living (L2L). L2L has several areas but can be defined as repurposing or converting logistical items (bags, containers, foam, components, etc.) into useful crew items or life support augmentation on-orbit after they have provided their primary logistics function. The intent is that by repurposing items, dedicated crew items do not have to be launched and overall launch mass is decreased. For non-LEO missions, the vehicle interior volume will be relatively fixed so L2L will enable this volume to be used more effectively through reuse and rearrangement of logistical components. Past work in the area of L2L has already conceptually developed several potential technologies [Howe, Howard 2010]. Several of the L2L concepts that have shown the most potential in the past are based on NASA cargo transfer bags (CTBs) or their equivalents which are currently used to transfer cargo to and from the ISS. A high percentage of all logistics supplies are packaging mass and for a 6-month mission a crew of four might need over 100 CTBs. These CTBs are used for on-orbit transfer and storage but eventually becomes waste after use since down mass is very limited. The work being done in L2L also considering innovative interior habitat construction that integrate the CTBs into the walls of future habitats. The direct integration could provide multiple functions: launch packaging, stowage, radiation protection, water processing, life support augmentation, as well as structure. Reuse of these CTBs would reduce the amount of waste generated and also significantly reduce future up mass requirements for exploration missions. Also discussed here is the L2L water wall , an innovative reuse of an unfolded CTB as a passive water treatment system utilizing forward osmosis. The bags have been modified to have an inner membrane liner that allows them to purify wastewater. They may also provide a structural water-wall element that can be used to provide radiation protection and as a structural divider. Integration of the components into vehicle/habitat architecture and consideration of operations concepts and human factors will be discussed. In the future these bags could be designed to treat wastewater, concentrated brines, and solid wastes, and to dewater solid wastes and produce a bio-stabilized construction element. This paper will describe the follow-on work done in design, fabrication and demonstrations of various L2L concepts, including advanced CTBs for reuse/repurposing, internal outfitting studies and the CTB-based forward osmosis water wall.

Shull, Sarah A.↗

DOE OSTI2024

DECOVALEX-2023: Task D Final Report

Task D of DECOVALEX-2023 is focused on the simulation of the coupled thermal hydraulic-mechanical (THM) behaviour in the full-scale engineered barrier system (EBS). The Horonobe EBS experiment is the demonstration of the full-scale EBS in the underground research laboratory (URL) (performed by JAEA in the Horonobe URL in Japan). Task D consisted of the three steps, a preliminary step (Step 0), simulation of the laboratory tests (Step 1) and simulation of the in-situ full-scale EBS experiment (Step 2). Since the Horonobe EBS experiment demonstrates the vertical emplacement option of the EBS, the experiment gallery is also backfilled with the backfill material. Therefore, interaction between the EBS and the backfill material can also be demonstrated, such as deformation (change of density) of the buffer material. The underground water in the Horonobe URL is saline. This fact adds chemical processes to THM behaviour. For example, mechanical properties (such as swelling pressure of the buffer material and backfill material) and hydraulic properties (such as permeability of the buffer material and backfill material) change depending on the water chemistry. Task D was therefore a challenging Task focused on not only the relatively simple THM behaviour but also complex THM behaviour including chemical processes. Six research teams (BGR, CAS, JAEA, KAERI, SNL and Taipower) participated the Task D. BGR, CAS, JAEA, KAERI and Taipower research teams selected a THM approach, while the SNL research team selected a TH approach. Step 1 involved the simulation of laboratory test results and was important to check the numerical codes developed by the research teams. Step 1 was divided into four sub steps. The simulation results through the Step 1 identified the parameters for simulation of the Step 2. Basic parameters of the materials (buffer material, backfill material, rock mass, concrete, sand) were provided by JAEA. Special parameters which research team needed were identified by back analysis of Step 1. Most notably the mechanical behaviour of swelling and displacement depended on the applied model (elastic model or elastoplastic model). Parameters such as Young’s modulus were found to need smaller values than characterised in the fundamental laboratory test results (Step 1-1, 1-2) for the elastic model. Although laboratory experiments are usually simple, test results contained some error. For example, if the saturation level is 100 % or higher, it should be considered an error. This situation was presented in the Step 1-3. A possible reason is that the buffer material is a mixture of bentonite and silica sand. When a specimen is cut to measure volume or weight, sand grains will affect the measurement data. In Step 2, boundary conditions such as temperature on the surface of the simulated overpack, heater power of the electrical heaters installed in the simulated overpack, injection pressure and inflow rate of the test water, were applied. The outer boundary conditions can be selected using measured data (injection pressure and inflow rate of the test water that is controlled by the injection systems installed in the sand layer around the buffer material and in the boundary between backfill material and concrete support). Since such measured data has some noise, research teams developed their own simplified boundary conditions. Inner boundary conditions can be selected using measured data as heater power and temperature on the surface of the simulated overpack. These data also contain some noise, so research teams developed their own simplified developed boundary conditions. Task D validated various approaches thorough the simulation of the in-situ full scale EBS system including backfill of the gallery: variations in the coupling processes (THM or THC), analysis codes, and boundary conditions. Temperature distribution in the buffer material was simulated well by all research teams. This means thermal behaviour is not sensitive to the simulation approaches. Although the water content distribution on the outside of the buffer material was well simulated by all research teams, the simulation results differ from the measured values inside the buffer material (at the centre and inside, near the simulated overpack). The buffer material is made from tap water, but in the in-situ experiment, saline groundwater infiltrates the buffer material. Therefore, the selection of the hydraulic parameters of the buffer material greatly affects the simulation results of the re saturation behaviour of the buffer material. In the Horonobe EBS experiment, measured values suitable for validating the simulation results were not obtained near the simulated overpack. When simulating the pressure and deformation of the buffer material, the measurement data is easily affected by the installation conditions of the measurement sensors, so verifying the measurement data itself remains an issue. Mechanical simulation results differ depending on whether they are considered as elastic or elastoplastic phenomena. The accuracy of measured in-situ data can be assessed by detailed analysis comparing sampling specimen analysis and measured data. The Horonobe EBS experiment is scheduled to be dismantled in the future (FY2026 and 2027). This detailed dismantling investigation will finally confirm the measured data.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

NASA NTRS2003

Unified Models of Turbulence and Nonlinear Wave Evolution in the Extended Solar Corona and Solar Wind

The PI (Cranmer) and Co-I (A. van Ballegooijen) made significant progress toward the goal of building a "unified model" of the dominant physical processes responsible for the acceleration of the solar wind. The approach outlined in the original proposal comprised two complementary pieces: (1) to further investigate individual physical processes under realistic coronal and solar wind conditions, and (2) to extract the dominant physical effects from simulations and apply them to a one-dimensional and time-independent model of plasma heating and acceleration. The accomplishments in the report period are thus divided into these two categories: 1a. Focused Study of Kinetic MHD Turbulence. We have developed a model of magnetohydrodynamic (MHD) turbulence in the extended solar corona that contains the effects of collisionless dissipation and anisotropic particle heating. A turbulent cascade is one possible way of generating small-scale fluctuations (easy to dissipate/heat) from a pre-existing population of low-frequency Alfven waves (difficult to dissipate/heat). We modeled the cascade as a combination of advection and diffusion in wavenumber space. The dominant spectral transfer occurs in the direction perpendicular to the background magnetic field. As expected from earlier models, this leads to a highly anisotropic fluctuation spectrum with a rapidly decaying tail in the parallel wavenumber direction. The wave power that decays to high enough frequencies to become ion cyclotron resonant depends on the relative strengths of advection and diffusion in the cascade. For the most realistic values of these parameters, though, there is insufficient power to heat protons and heavy ions. The dominant oblique waves undergo Landau damping, which implies strong parallel electron heating. We thus investigated the nonlinear evolution of the electron velocity distributions (VDFs) into parallel beams and discrete phase-space holes (similar to those seen in the terrestrial magnetosphere) which are an alternate means of heating protons via stochastic interactions similar to particle-particle collisions. 1b. Focused Study of the Multi-Mode Detailed Balance Formalism. The PI began to explore the feasibility of using the "weak turbulence," or detailed-balance theory of Tsytovich, Melrose, and others to encompass the relevant physics of the solar wind. This study did not go far, however, because if the "strong" MHD turbulence discussed above is a dominant player in the wind's acceleration region, this formalism is inherently not applicable to the corona. We will continue to study the various published approaches to the weak turbulence formalism, especially with an eye on ways to parameterize nonlinear wave reflection rates. 2. Building the Unified Model Code Architecture. We have begun developing the computational model of a time-steady open flux tube in the extended corona. The model will be "unified" in the sense that it will include (simultaneously for the first time) as many of the various proposed physical processes as possible, all on equal footing. To retain this generality, we have formulated the problem in two interconnected parts: a completely kinetic model for the particles, using the Monte Carlo approach, and a finite-difference approach for the self-consistent fluctuation spectra. The two codes are run sequentially and iteratively until complete consistency is achieved. The current version of the Monte Carlo code incorporates gravity, the zero-current electric field, magnetic mirroring, and collisions. The fluctuation code incorporates WKJ3 wave action conservation and the cascade/dissipation processes discussed above. The codes are being run for various test problems with known solutions. Planned additions to the codes include prescriptions for nonlinear wave steepening, kinetic velocity-space diffusion, and multi-mode coupling (including reflection and refraction).

Cranmer, Steven R.↗

DOE OSTI2020

In situ monitoring of the stopped muon flux at Mu2e

The specific object of Mu2e search is the neutrino-less coherent $\mu^-\rightarrow e^-$ conversion in the field of an aluminum nucleus and the signal is a monoenergetic electron of energy $ \approx 104.97$ MeV. This process is forbidden in the Standard Model but allowed in many of its extensions: with minimal changes to include neutrino masses and oscillation, the Branching Ratio of this process, or the similarly interesting $\mu\rightarrow e \gamma$ decay, is expected to be of the order of $\mathcal{O}(10^{-54})$. Values like these are below any currently achievable experimental sensitivity, and the observation of a signal would be unambiguous evidence of New Physics. The upper limit on the muon conversion was set by SINDRUM II at $7\times10^{-13}$ (90\% C.L.) and the goal of the Mu2e collaboration is an improvement of 4 order of magnitudes.\\ The Mu2e experiment can be conceptually divided into three stages: interaction of the primary proton pulse with the tungsten target a nd production of $\pi$ and $\mu$; collection and transport of the produced particles down to the aluminum stopping target; interaction with the stopping target and measurement of the output particles. The measurement of the stopped muon flux is of cardinal importance for the normalization of the Mu2e results. The baseline design includes two detectors developed to measure this flux, counting $\gamma$ emitted by the muonic atoms: HPGe and LaBr$_3$(Ce). The studies performed by the Mu2e Collaboration show that these two systems can be reliably used to determine the overall normalization.\\ The number of stopped muons is proportional to the number of protons on target, which itself depends on the extraction system for the proton beam. The Mu2e beam delivery system will use a resonant extraction to create the proton pulses and this method is characterized by intensity fluctuations on the time scale of milliseconds. These fluctuations, aside from the proportional fluctuation of t he number of muons, have a non-trivial impact on the overall! performance of the apparatus. For example, higher intensity proton bunches translate into higher veto rates by the Cosmic Rays Veto system, which is triggered by secondary particles generated by the proton interaction. Another example is the effect of the instantaneous luminosity on the reconstruction efficiency: higher intensity translates into higher detector occupancy, which in turn reduces the reconstruction efficiency.\\ The two cited detectors encounter limitations when trying to monitor the flux at the millisecond timescales: HPGe is by construction a slow detector while the rate of emission of the $\gamma (1809\textrm{ keV})$ the LaBr$_3$(Ce) will measure is too low. As of today, no system allows to monitor these fluctuations and the goal of this Thesis is to try filling this gap. The method we developed to monitor the fluctuations of the stopped muon flux, and by extension, the fluctuation in the proton pulses intensity is based on counting the number of muons capt ured in the stopping target by counting the number of protons produced in the process of the nuclear muon capture. This is one of the possible processes a stopped muon can undergo and it can lead to the ejection of charged particles.\\ The reconstruction algorithms were developed by the Mu2e collaboration with the explicit (and sometimes implicit) aim to reconstruct tracks of electrons in an energy range of a conversion electron. The protons we are interested in behave quite differently from these electrons and the cardinal task of this study has been to tailor the reconstruction routines to these particles. This has been done by analyzing simulated events consisting of a single proton to optimize the procedure and understand the features of these tracks. A major difference with respect to the electrons is the non-relativistic nature of $100\div300$ MeV$/c$ protons: accounting for the low velocity is both a challenge and a unique signature of the tracks.\\ The use of a pulse d beam and the need to allow produced $\pi$ to decay and $\m! u$ to int! eract in the stopping target forces the definition of a Mu2e \textit{event}: an event is comprised by everything which happens in a $1.7\ \mu$s window between the incidence two consecutive proton pulses. This means that the reconstruction will be performed in an environment crowded by many different particles. The next step of this study has been to estimate the performance of the method in a fully-fledged Mu2e simulated event and to study what kind of selections would improve it.\\ Our results are satisfactory since the number of reconstructed protons is significant (few per event), and show that a monitor on the timescale of milliseconds is indeed possible.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRSDate not supplied

Monitoring Airspace Complexity and Determining Contributing Factors

The national airspace has evolved over many years to accommodate increased traffic demand [1] while simultaneously maintaining one of the safest forms of transportation [2], [3]. One of the reasons for this success is the ability of the system and the operators to adapt and accommodate to situations that routinely disrupt optimal operations. These situations may include: adverse weather, delays, early arrivals, equipment outages, and other factors that are outside the operators ability to control. These factors can lead to states where automation is unable to properly handle these issues and therefore air traffic controllers and pilots have to intervene, ultimately increasing communication between operators resulting in higher workload. As controller workload increases to handle sub-optimal operating conditions this can be viewed as an increase in complexity. The reasoning for this is because humans are now required to make tactical decisions in response to external factors, resulting in a departure from the strategic plan where operations would be more efficiently managed. Human operators control airspace complexity under rigid regulations that are constantly changing. The airspace is divided into sectors and the number of aircraft assigned to each controller is limited for safe handling. There has been past work that devised airspace complexity metrics in commercial aviation and related these metrics to controller workload (e.g., [4],[5]). The upper bounds on the system load are pre-determined. Such bounds on complexity make for a safe system, but the system cannot scale and adapt to autonomous, dense, and heterogeneous traffic, including the many types of Unmanned Aerial Vehicles (UAVs) envisioned to be added to the operations. We hypothesize that, as traffic density and heterogeneity grow, and other key metrics change, there will be phase transitions at which the way traffic should be managed changes significantly [6]. We offer a method for in-time detection of contributing factors that lead to phase transitions, characterized by increased complexity. To the best of our knowledge, there is no tool similar to our proposed effort that identifies such contributing factors or precursor patterns. To define the scope we are proposing to measure complexity from the viewpoint of the Terminal Radar Approach Control Facilities (TRACON) controller’s perspective. In particular we are analyzing arrivals into KSFO. With safety as the top concern for airspace operators, it is important to recognize that as density and heterogeneity grow, the focus of the system will change. Times of the day when the airspace has low density and heterogeneity, the flights will follow more efficient paths where the aircraft move on established routes that are more or less directly to the destination. However, when density and heterogeneity increases, the system will begin changing focus to avoiding conflicts and collisions and route the flights in a more flexible way. Higher flexibility requires more communication and coordination between controllers and pilots which the current automation is unable to handle. This paper proposes a novel approach that monitors airspace complexity at multiple scales, uses a Machine Learning-based tool that predicts when operations will transition to a regime of greater complexity, and identifies actions that can reduce the complexity while still maintaining efficient and safe operations. We demonstrate our proposed approach using data from multiple complementary sources. This includes, but is not limited to: historical aircraft surveillance data from NASA’s Sherlock Data Warehouse [7], METAR weather data, and airport configuration data from Aviation System Performance Metrics (ASPM). The surveillance data flight paths are sampled at a variable sample rate — increasing as the aircraft approaches the airport. This is due to how Sherlock manages flight track stitching between different radar facilities which have different sampling rates. The weather and performance data are logged at defined intervals throughout the day at a courser refresh rate. In addition to the logged data and metrics, we leverage pre-defined Standard Terminal Arrival Routes (STARs) procedures to characterize the path of each flight. Each flight files for one of these routes in the flight plan well before entering the terminal airspace, and approximately follows the route until it leaves the STAR, typically on the final fix of a runway transition. However, most flights do not always fly the full STAR procedure to completion [8], but the majority do adhere to the fixes within the common route of the procedure. Our approach leverages fixes in the common route of each of the STARs to build a reference path to the airport. This allows us to characterize the flight paths in what we are defining as the “maneuvering area” (the airspace between the STAR and before the flight is lined up on the runway’s final approach) to determine how off nominal the flights are to calculate its complexity score. Determining airspace complexity is a concept that does not have a concrete answer. In designing this metric, we consider what increases the workload for the air traffic controllers. Consequently more specialized vectoring maneuvers results in higher workload. Accordingly, we start with a theory: each flight has a direct path it takes from the STAR’s common route to the final approach’s outer marker fix for the flight’s landing runway. It is important to note that the direct path is only used as a reference. If the majority of the flights have a large consistent offset as compared to other routes it does not necessarily mean that those flights have higher complexity. We are merely building a distribution based on this direct path for that particular STAR and runway pair to determine the normal mode of operations for that route. Flights that are in the upper tail of these distributions will result in higher complexity scores and flights that fly in the median will represent the normal mode of operations and therefore will have lower complexity scores. Since flights following each STAR route take different paths to the airport, we have a different distribution for each STAR route and therefore can model these distributions to compute a complexity score from their respective normalized distributions. To evaluate the effectiveness of our proposed airspace complexity metric we will compare against an established approach based on trajectory clustering [9]. This unsupervised learning technique consists of the following steps: (1) identify the general maneuvering areas (waypoints) by performing $\kappa$-means or DBSCAN clustering on locations where aircraft frequently turn based on the surveillance radar track data, (2) map flight trajectories onto sequences of waypoints, and (3) cluster the sequences based on their common subsequences. From a high-level perspective, this baseline model learns nominal operations in the airspace through the sequence of waypoints that are representative of where aircraft change direction and defines deviations from the nominal operations as “complex.” Therefore, more deviations from the nominal operations correspond to higher complexity values. For our validation, we re-implemented this technique and tune model hyper-parameters to correctly detect waypoints for the arrival traffic into the San Francisco bay area. We will compute the complexity measure over a one-year period using our proposed technique as well as the baseline. Our validation will be based on each technique’s ability to detect a set of undesirable outcomes (e.g., go-arounds, holding patterns, average time in the airspace, etc.). Since our current complexity metric is derived from the offset from the direct reference path, it’s important to understand what causes these offsets. In many of the flights with high offset distance, flights performing holding patterns and S turns can be observed. These maneuvering tactics are utilized to add distance between the aircraft and the destination runway to prevent multiple flights from having conflicting arrival times. In order to predict a rise in complexity (or the precursor to complexity), it’s necessary to be able to identify these potential conflicts (which in turn, result in higher offsets). To do this, we define a “representative flight” for each STAR route and runway pair. This flight is approximately the path the flight would take if there was a clear path with no other flights in the airspace — including the time remaining to the airport. We first identify the flights for a given STAR runway pair using the offset to the reference path distributions that fall between the 44-55 percentiles. This yields the flights that conform to the most normal mode of operation. Each of these flights is partitioned based on the percent complete from the entry point into the maneuvering areas from 0\% – 100\% complete. Then for each percent “bin”, we take the median value of the flight’s latitude/longitude coordinates, airspeed, and (non causal) time remaining to the airport to construct a lookup table for each percent complete bin on a given route. As a flight enters the maneuvering area, we can find the estimated arrival time of a flight to the airport by finding the closest point to the representative path’s percent complete bin (relative to the flight’s current position at any snapshot in the airspace) and therefore retrieve the corresponding remaining time left on the “representative path”. We assume that the flight will follow the representative path to completion when deriving these estimates. We can then compare these estimated arrival times against other flights for the same snapshot in time to identify potential conflicts. If more flights are estimated to arrive within a tolerance window than there are runways available, then we have a potential conflict. We can use this derived measure along with other factors expected to add disruption to the operation such as weather and runway configuration changes as an input to machine learning tools to detect precursors that increases in our complexity measure. This novel method will assist in uncovering insights into the contributing factors that lead to increased complexity that may allow for in-time responses to avoid reaching a high complexity state in the airspace.

complexity↗

NASA NTRS1999

Instrumentation Working Group Summary

The Instrumentation Working Group compiled a summary of measurement techniques applicable to gas turbine engine aerosol precursors and particulates. An assessment was made of the limits, accuracy, applicability, and technology readiness of the various techniques. Despite advances made in emissions characterization of aircraft engines, uncertainties still exist in the mechanisms by which aerosols and particulates are produced in the near-field engine exhaust. To adequately assess current understanding of the formation of sulfuric acid aerosols in the exhaust plumes of gas turbine engines, measurements are required to determine the degree and importance of sulfur oxidation in the turbine and at the engine exit. Ideally, concentrations of all sulfur species would be acquired, with emphasis on SO2 and SO3. Numerous options exist for extractive and non-extractive measurement of SO2 at the engine exit, most of which are well developed. SO2 measurements should be performed first to place an upper bound on the percentage of SO2 oxidation. If extractive and non-extractive techniques indicate that a large amount of the fuel sulfur is not detected as SO2, then efforts are needed to improve techniques for SO3 measurements. Additional work will be required to account for the fuel sulfur in the engine exhaust. Chemical Ionization Mass Spectrometry (CI-MS) measurements need to be pursued, although a careful assessment needs to be made of the sampling line impact on the extracted sample composition. Efforts should also be placed on implementing non-intrusive techniques and extending their capabilities by maximizing exhaust coverage for line-of-sight measurements, as well as development of 2-D techniques, where feasible. Recommendations were made to continue engine exit and combustor measurements of particulates. Particulate measurements should include particle size distribution, mass fraction, hydration properties, and volatile fraction. However, methods to ensure that unaltered samples are obtained need to be developed. Particulate speciation was also assigned a high priority for quantifying the fractions of carbon soot, PAH, refractory materials, metals, sulfates, and nitrates. High priority was also placed on performing a comparison of particle sizing instruments. Concern was expressed by the workshop attendees who routinely make particulate measurements about the variation in number density measured during in-flight tests by different instruments. In some cases, measurements performed by different groups of researchers during the same flight tests showed an order of magnitude variation. Second priority was assigned to measuring concentrations of odd hydrogen and oxidizing species. Since OH, HO2, H2O2, and O are extremely reactive, non-extractive measurements are recommended. A combination of absorption and fluorescence is anticipated to be effective for OH measurements in the combustor and at the engine exit. Extractive measurements of HO2 have been made in the stratosphere, where the ambient level of OH is relatively low. Use of techniques that convert HO2 to OH for combustor and engine exit measurements needs to be evaluated, since the ratio of HO2/OH may be 1% or less at both the combustor and engine exit. CI-MS might be a viable option for H2O2, subject to sampling line conversion issues. However, H2O2 is a low priority oxidizing species in the combustor and at the engine exit. Two candidates for atomic oxygen measurements are Resonance Enhanced Multi-Photon Ionization (REMPI) and Laser-Induced Fluorescence (LIF). Particulate measurement by simultaneous extractive and non-extractive techniques was given equal priority to the oxidizer measurements. Concern was expressed over the ability of typical ground test sampling lines to deliver an unaltered sample to a remotely located instrument. It was suggested that the sampling probe and line losses be checked out by attempting measurements using an optical or non-extractive technique immediately upstream of the sampling probe. This is a possible application for Laser Induced Incandescence (LII) as a check on the volume fraction of soot. Optical measurements of size distribution are not well developed for ultrafine particles less than about 20 nm in diameter, so a non-extractive technique for particulate size distribution cannot be recommended without further development. Carbon dioxide measurements need to be made to complement other extractive measurement techniques. CO2 measurements enable conversion of other species concentrations to emission indices. Carbon monoxide, which acts as a sink for oxidizing species, should be measured using non-extractive techniques. CO can be rapidly converted to CO2 in extractive probes, and a comparison between extractive and non-extractive measurements should be performed. Development of non-extractive techniques would help to assess the degree of CO conversion, and might be needed to improve the concentration measurement accuracy. Measurements of NO(x) will continue to be critical due to the role of NO and NO2 in atmospheric chemistry, and their influence on atmospheric ozone. Time-resolved measurements of temperature, velocity, and species concentrations were included on the list of desired measurement. Thermocouples are typically adequate for engine exit measurements. PIV and LDV are well established for obtaining velocity profiles. The techniques are listed in the accompanying table; are divided into extractive and non-extractive techniques. Efforts were made to include a measurement uncertainty for each technique. An assessment of the technology readiness was included.

Zaller, Michelle↗

DOE OSTI2021

The Role of Quantum Science Concepts in Enhancing Sensing and Imaging Technologies: Applications for Biology: Proceedings of a Workshop

Quantum concepts hold the potential to enable significant advances in sensing and imaging technologies that could be vital to the study of biological systems. The workshop Quantum Science Concepts in Enhancing Sensing and Imaging Technologies: Applications for Biology, held online March 8–10, 2021, was organized to examine the research and development needs to advance biological applications of quantum technology. Hosted by the National Academies of Sciences, Engineering, and Medicine, the event brought together experts working on state-of-the-art, quantum-enabled technologies and scientists who are interested in applying these technologies to biological systems. Through talks, panels, and discussions, the workshop facilitated a better understanding of the current and future biological applications of quantum-enabled technologies in fields such as microbiology, molecular biology, cell biology, plant science, mycology, and many others. The workshop was organized around three main themes. The first, quantum in biology, examined quantum concepts that are hypothesized to be important for life processes and that researchers are working to observe through biological imaging and sensing. The second, quantum for biology, addressed ways to use quantum concepts to enhance technologies for biological imaging and sensing. The third, biology for quantum, offered a wider discussion of how the frontiers of biological imaging and sensing could enable future study using quantum concepts, tools, or technologies. Throughout the workshop, participants identified a wide range of emerging approaches and opportunities at the intersection of quantum physics and biological sensing and imaging. During the workshop, there were some differences in how each speaker defined the term quantum. During one of the panels, Prem Kumar offered thoughts on what phenomena are classical versus quantum, explaining that techniques get progressively more quantum as you move from just having superposition to having superposition with measurement and entanglement. Another explanation from Clarice Aiello delineates the definition into several levels. This includes a base level of “quantum-ness,” which reflects that all matter is made of atoms, and when these particles are isolated they behave based on quantum mechanical principles. A second level is related to quantum coherence, where a single quantum object might be found in a coherent superposition state. A final level, which she described as the quantum-entangled level, involves multiple quantum systems which are entangled among themselves. Overall, the workshop touched on concepts such as superposition, entanglement, and squeezing and their potential implications for communication, computing, and simulation, in addition to the workshop’s main focal area, biological sensing and imaging. At the opening of the workshop, Thorsten Ritz of the University of California, Irvine, identified two questions at the heart of quantum biology: Is the machinery of life quantum mechanical, and can quantum mechanics be used to study the machinery of life in new ways? Participants highlighted systems in which researchers have explored these questions, from the vast array of molecular interactions involved in biological processes such as photosynthesis, to the mechanics involved in cellular functions such as differentiation and aggregation, to the role of oscillating magnetic fields in flight orientation among birds. Sensing and imaging technologies are crucial to biological research; these technologies could both enhance the study of quantum effects and be enhanced by quantum concepts. A critical challenge in biological research is to develop imaging and sensing tools that do not damage or interfere with the often fragile and fleeting systems being studied. Attendees discussed a variety of established and emerging technologies that could enhance noninvasive biological imaging, including single- and two-photon spectroscopy, single-molecule spectroscopy, quantum illumination, ghost imaging, and cryo-electron microscopy. One example came from Marlan Scully who gave a keynote address on the first day of the workshop. Scully emphasized the use of different laser technologies, which exhibit coherence and other quantum properties, in moving toward real-world biological applications, such as the detection of SARS-CoV-2. Both tools and theory will play an important role in advancing quantum biology research and applications. Several participants suggested theorists and experimentalists should work in tandem to understand and model biological processes. While physics often reduces systems to their simplest forms for fundamental insights, participants also noted the value of observing and understanding biological systems in all their “messiness,” capturing both the inner workings of biological systems and the complex interactions that occur within and between organisms. In discussions among participants, several attendees stressed the need to match emerging tools with the right scientific questions. Rather than developing quantum technologies as “a hammer looking for a nail,” participants emphasized a focus on exploring the problems these technologies are best suited to address. For example, it is important to consider the size of the phenomenon being studied, the timescales that are important in answering the scientific question, and other relevant considerations. Every tool along the spectrum from classical to quantum involves its own set of trade-offs. For example, Ted Laurence of the Lawrence Livermore National Laboratory said that quantum measurements, such as single photon counting, fluorescence transitions, and lasers, can take longer to produce the same results as classical measurements. These quantum measurements, however, do not require calibration and can enable new research questions to be answered. Prem Kumar, Northwestern University, noted that, despite their promise, quantum approaches should not be used simply for the sake of using quantum, especially in situations where classical approaches better meet the needs of the researcher. Understanding and applying quantum concepts could enable advances in a wide range of application areas including energy, synthetic biology, medicine, and sustainability. For example, Michelle O’Malley, University of California, Santa Barbara, described how improved noninvasive imaging approaches could help capture the complex interactions and functions involved in the breakdown of organic matter by microbial communities and lead to new technologies for capturing valuable products from plant waste. Several other participants discussed needs in tracking the movement of metabolites and molecules in microbial communities for insights into nutrient cycling in environments such as soil. Margaret Ahmad, Sorbonne University, discussed potential opportunities to leverage the magnetic properties of cryptochromes to advance new treatment approaches for diseases such as COVID-19 and cancer. Looking toward the future development of the field, participants discussed challenges to advancing quantum biology that arise from disciplinary disconnects between physicists and biologists. The siloing of academic research disciplines represents a significant barrier to progress. Disconnects in terminology, motivations and priorities, and structural barriers to collaborative work underscore the need for concerted efforts to bridge these divides. Attendees and speakers offered suggestions for resolving these divergences, establishing a shared language, moving the field forward, and fostering meaningful feedback between disciplines. Overall, participants stressed a need for balance, open communication, collaboration, unity, and clear dialogue on trade-offs between quantum and classical approaches. Keiko Torii, The University of Texas at Austin, said that the best collaborations happen when the project provides mutual advantages that can show off everyone’s talents, each team finds the work interesting, and partners develop a camaraderie to pursue new knowledge. To enable near- and long-term opportunities in this space, participants suggested that exploratory, high-risk funding could improve existing instrumentation to explore quantum enhancement collaboratively. They also emphasized the need for collaboration between quantum physicists and sensing/imaging scientists, which could be advanced through a dedicated quantum biology investigator program. People, even more than technology, will be crucial to the future of quantum biology. Participants explored training, education, and workforce needs to further develop this burgeoning field and cultivate the next generation of scientists. While many programs are still in their nascent stages, participants highlighted examples of approaches and programs being developed at various types of institutions to engage students and professional scientists in quantum biology research. Several participants stressed the need for an inclusive approach, spanning disciplines as well as communities to foster a diverse field fueled by the intellectual contributions of a wide range of people, including historically under-resourced schools and students. To increase awareness and excitement about quantum physics and related areas of biology, attendees suggested capitalizing on the “buzz” around quantum. Several participants emphasized the need to start early, introducing students to quantum concepts and their appealing “weirdness” in K–12 education. Engaging students early—before they become entrenched in traditional disciplinary siloes as typically happens in graduate school—could help to galvanize interest in the area and foster a generation of scientists with the interdisciplinary mindset and skills needed to advance this interdisciplinary field. While these efforts could be advanced at many levels and across multiple sectors, several participants suggested a national quantum biology center could be a valuable hub to coordinate and support quantum biology education and workforce development across academia, industry, and government.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS2021

Pluto’s Sputnik Planitia: Composition of geological units from infrared spectroscopy

We have compared spectroscopic data of Sputnik Planitia on Pluto, as acquired by New Horizons’ Linear Etalon Imaging Spectral Array (LEISA) instrument, to the geomorphology as mapped by White et al. (2017) using visible and panchromatic imaging acquired by the LOng-Range Reconnaissance Imager (LORRI) and the Multi-spectral Visible Imaging Camera (MVIC). We have focused on 13 of the geologic units identified by White et al. (2017), which include the plains and mountain units contained within the Sputnik basin. We divided the map of Sputnik Planitia into 15 provinces, each containing one or more geologic units, and we use LEISA to calculate the average spectra of the units inside the 15 provinces. Hapke-based modeling was then applied to the average spectra of the units to infer their surface composition, and to determine if the composition resulting from the modeling of LEISA spectra reflects the geomorphologic analyses of LORRI data, and if areas classified as being the same geologically, but which are geographically separated, share a similar composition. We investigated the spatial distribution of the most abundant ices on Pluto’s surface - CH4, N2, CO, H2O, and a non-ice component presumed to be a macromolecular carbon-rich material, termed a tholin, that imparts a positive spectral slope in the visible spectral region and a negative spectral slope longward of ~1.1 μm. Because the exact nature of the non-ice component is still debated and because the negative spectral slope of the available tholins in the near infrared does not perfectly match the Pluto data, for spectral modeling purposes we reference it generically as the negative spectral slope endmember (NSS endmember). We created maps of variations in the integrated band depth (from LEISA data) and areal mass fraction (from the modeling) of the components. The analysis of correlations between the occurrences of the endmembers in the geologic units led to the observation of an anomalous suppression of the strong CH4 absorption bands in units with compositions that are dominated by H2O ice and the NSS endmember. Exploring the mutual variation of the CH4 and N2 integrated band depths with the abundance of crystalline H2O and NSS endmember revealed that the NSS endmember is primarily responsible for the suppression of CH4 absorptions in mountainous units located along the western edge of Sputnik Planitia. Our spectroscopic analyses have provided additional insight into the geological processes that have shaped Sputnik Planitia. A general increase in volatile abundance from the north to the south of Sputnik Planitia is observed. Such an increase first observed and interpreted by Protopapa et al., 2017 and later confirmed by climate modeling (Bertrand et al., 2018) is expressed geomorphologically in the form of preferential deposition of N2 ice in the upland and mountainous regions bordering the plains of southern Sputnik Planitia. Relatively high amounts of pure CH4 are seen at the southern Tenzing Montes, which are a natural site for CH4 deposition owing to their great elevation and the lower insolation they are presently receiving. The NSS endmember correlates the existence of tholins within certain units, mostly those coating the low-latitude mountain ranges that are co-latitudinal with the tholin-covered Cthulhu Macula. The spectral analysis has also revealed compositional differences between the handful of occurrences of northern non-cellular plains and the surrounding cellular plains, all of which are located within the portion of Sputnik Planitia that is presently experiencing net sublimation of volatiles, and which do not therefore exhibit a surface layer of bright, freshly-deposited N2 ice. The compositional differences between the cellular and non-cellular plains here hint at the effectiveness of convection in entraining and trapping tholins within the body of the cellular plains, while preventing the spread of such tholins to abutting non-cellular plains.

Pluto's Sputnik Planitia↗

NASA NTRS2000

Phase-Shifting Liquid Crystal Point-Diffraction Interferometry

Microgravity fluid physics experiments frequently measure concentration and temperature. Interferometers such as the Twyman Green illustrated have performed full-field measurement of these quantities. As with most such devices, this interferometer uses a reference path that is not common with the path through the test section. Recombination of the test and reference wavefronts produces interference fringes. Unfortunately, in order to obtain stable fringes, the alignment of both the test and reference paths must be maintained to within a fraction of the wavelength of the light being used for the measurement. Otherwise, the fringes will shift and may disappear. Because these interferometers are extremely sensitive to bumping, jarring and transmitted vibration, they are typically mounted on optical isolation tables. Schlieren deflectometers or the more recent Shack-Hartmann wavefront sensors also measure concentration and temperature in laboratory fluid flows. Ray optics describe the operation of both devices. In a schlieren system, an expanded, collimated beam passes through a test section where refractive index gradients deflect rays. A lens focuses the beam to a filter placed in the rear focal plane of the decollimating lens. In a quantitative color schlieren system, gradients in the index of refraction appear as colors in the field of view due to the action of the color filter. Since sensitivity is a function of the focal length of the decollimating lens, these systems are rather long and filter fabrication and calibration is rather difficult. A Shack-Hartmann wavefront sensor is an array of small lenslets. Typical diameters are on the order of a few hundred microns. Since these lenslets divide the test section into resolution elements, the spatial resolution can be no smaller than an individual lenslet. Such a device was recently used to perform high-speed tomography of heated air exiting a 1.27 cm diameter nozzle. While these wavefront sensors are very compact, the limited spatial resolution and the methods required for data reduction suggest that a more useful instrument needs to be developed. The category of interferometers known as common path interferometers can eliminate much of the vibration sensitivity associated with traditional interferometry as described above. In these devices, division of the amplitude of the wavefront following the test section produces the reference beam. Examples of these instruments include shearing and point diffraction interferometers. In the latter case, shown schematically, a lens focuses light passing through the test section onto a small diffracting object. Such objects are typically either a circle of material on a high quality glass plate or a small sphere in a glass cell. The size of the focused spot is several times larger than the object so that the light not intercepted by the diffracting object forms the test beam while the diffracted light generates a spherical reference beam. While this configuration is mechanically stable, phase shifting one beam with respect to the other is difficult due to the common path. Phase shifting enables extremely accurate measurements of the phase of the interferogram using only gray scale intensity measurements and is the de facto standard of industry. Mercer and Creath 2 demonstrated phase shifting in a point diffraction interferometer using a spherical spacer in a liquid crystal cell as the diffracting object. By changing the voltage across the cell, they were able to shift the phase of the undiffracted beam relative to the reference beam generated by diffraction from the sphere. While they applied this technology to fluid measurements, the device shifted phase so slowly that it was not useful for studying transient phenomena. We have identified several technical problems that precluded operation of the device at video frame rates and intend to solve them to produce a phase-shifting liquid crystal point-diffraction interferometer operating at video frame rates. The first task is to produce high contrast fringes. Since the diffracted beam is much weaker than the transmitted beam, interferograms have poor contrast unless a dye is added to the liquid crystal to reduce the intensity of the undiffracted light. Dyes previously used were not rigorously characterized and suffered from hysteresis in both the initial alignment state of the device and the electro-optic switching characteristics. Hence, our initial effort will identify and characterize dyes that do not suffer from these difficulties and are readily soluble in the liquid crystal host. Since the ultimate goal of this research is to produce interferometers capable of phase shifting at video frame rates, we will quantify the difference in switching times between ferroelectric and nematic liquid crystals. While we have more experience with nematic crystals, they typically switch more slowly than ferroelectric cells. As part of that effort, we will investigate the difference in the modulation of the interferograms as a function of the type of liquid crystal in the cell. Because the temporal switching response of a liquid crystal cell is directly related its thickness, we intend to explore techniques required to produce cells that are as thin as possible. However, the cells must still produce a total phase shift of two pi radians.

Griffin, DeVon W.↗

NASA NTRS2014

Dividing the Concentrator Target From the Genesis Mission

The Genesis spacecraft, launched in 2001, traveled to a Lagrangian point between the Earth and Sun to collect particles from the solar wind and return them to Earth. However, during the return of the spacecraft in 2004, the parachute failed to open during descent, and the Genesis spacecraft crashed into the Utah desert. Many of the solar wind collectors were broken into smaller pieces, and the field team rapidly collected the capsule and collector pieces for later assessment. On each of the next few days, the team discovered that various collectors had survived intact, including three of four concentrator targets. Within a month, the team had imaged more than 10,000 fragments and packed them for transport to the Astromaterials Acquisition and Curation Office within the ARES Directorate at JSC. Currently, the Genesis samples are curated along with the other extraterrestrial sample collections within ARES. Although they were broken and dirty, the Genesis solar wind collectors still offered the science community the opportunity to better understand our Sun and the solar system as a whole. One of the more highly prized concentrator collectors survived the crash almost completely intact. The Genesis Concentrator was designed to concentrate the solar wind by a factor of at least 20 so that solar oxygen and nitrogen isotopes could be measured. One of these materials was the Diamond-on-Silicon (DoS) concentrator target. Unfortunately, the DoS concentrator broke on impact. Nevertheless, the scientific value of the DoS concentrator target was high. The Genesis Allocation Committee received a request for approximately 1 cm(sup 2) of the DoS specimen taken near the focal point of the concentrator for the analysis of solar wind nitrogen isotopes. The largest fragment, Genesis sample 60000, was designated for this allocation and needed to be precisely cut. The requirement was to subdivide the designated sample in a manner that prevented contamination of the sample and minimized the risk of losing or breaking the precious requested sample fragment. The Genesis curator determined that the use of laser scribing techniques to "cut" a precise line and subsequently cleave the sample (in a controlled break of the sample along that line) was the best method for accomplishing the sample subdivision. However, there were risks, including excess heating of the sample, that could cause some of the implanted solar wind to be lost via thermal diffusion. Accidentally breaking the sample during the handling and cleaving process was an additional risk. Early in fiscal year 2013, to address this delicate, complicated task, the ARES Directorate assembled its top scientists to develop a cutting plan that would ensure success when applied to the actual concentrator target wafer; i.e., to produce an approximately 1 cm(sup 2) piece from the requested area of the wafer. The team, subsequently referred to as the JSC Genesis Tiger Team, spent months researching and testing parameters and techniques related to scribing, cleaving, transporting, handling, and holding (i.e., mounting) the specimen. The investigation required considerable "thinking outside the box," and many, many trials using nonflight wafer analogs. After all preliminary testing, the following method was adopted as the final cutting plan. It was used in two final end-to-end practice runs before being used on the actual flight target wafer. The wafer was oriented on the laser cutting stage with the 100 and 010 directions of the wafer parallel to the corresponding X and Y directions of the cutting stage. The laser was programed to scribe 31 lines of the appropriate length along the Y stage direction. The programed scribe lines were separated by 5 micron in the X direction. The laser parameters were set as follows: (1) The laser power was 0.5 watts; (2) each line consisted of 50 passes, with the Z position being advanced 5 micron per pass; and (3) 30 s would elapse before the next line was scribed to allow for wafer cool down from any possible heating via the laser. The ablated material that "stuck" in the "scribe-cut" was removed from the "cut" using an ultrasonic micro-tool. After all the ablated silicon was removed from the wafer, the wafer was repositioned in exactly the same orientation on the laser stage. The laser was focused using the bottom of the wafer channel, and the 31-line scribing pattern described above was reprogrammed using the Z position of the groove bottom as the starting Z value instead of the top wafer surface, which was used previously. Upon completion of the second set of scribes, the ultrasonic micro-tool was again used to clean out the cut. The wafer was remounted on the stage in exactly the same orientation as before. The laser was again focused on the bottom of the groove. This time, however, the laser was.programed to scribe only one line down the exact center of the channel. The final scribe line consisted of 100 passes with a Z advance of 5 micron per pass and with the laser power set at 0.5 watts. As mentioned above, the final cutting plan was practiced in two end-to-end trials using non-flight, triangular-shaped silicon wafers similar in size and orientation to the actual DOS 60000 target sample. The actual scribing of the triangular-shaped wafers required scribing two lines and cleaving (i.e. scribe-cleave, then scribe-cleave) to obtain the piece requested for allocation. Early in December 2012, after many months of experiments and practicing and perfecting the techniques and procedures, the team successfully subdivided the Genesis DoS 60000 target sample, one of the most scientifically important samples from the Genesis mission (figure 2). On December 17, 2012, the allocated piece of concentrator target sample was delivered to the requesting principal investigator.The cutting plan developed for the subdivision of this sample will be used as the model for subdividing future requested Genesis flight wafers (appropriately modified for different wafer types).

Lauer, H. V., Jr.↗

DOE OSTI2023

Multi-Sourced Collaboration for the Production and Refining of Rare Elements and Critical Metals (Final Technical Report)

The project objective was to develop a feasible and cost-effective method for recovering rare earth elements (REEs) and critical materials (CMs) from coal and coal byproducts, resulting in high-purity individually separated REEs and CMs. The targeted REEs included Y, Pr, Nd, Gd, Dy, and Sm, with a purity of over 99.5%, while the CMs included Co, Mn, Ga, Sr, Li, Ni, Zn, and Ge, with a purity of over 90%. The project aimed to design a prototype facility capable of producing 1-3 tonnes/day of high-purity REO mixes. The work was divided into four designated circuits: 1) REE extraction and concentration, 2) REE separation and purification, 3) RE metal production, and 4) CM production. To achieve these goals, the project involved 11 tasks, including technology reviews, research, process flow diagram development, mass balance estimation, and preliminary technical-economic analysis. The project team included researchers from the University of Kentucky, University of Alabama and Virginia Tech as well as process specialists from Argonne National Laboratory. MP Materials provided technical support regarding rare earth markets and processing while Alliance Coal performed resource assessment. The project included a market analysis for Nd/Pr, Tb, Dy, Gd, Y, Co, Mn, Li, Sr, Ga, Ni, Zn, and Ge. These analyses provided insights into the supply and demand trends as well as historic and future projections of market price relative to purity requirements for these elements. Two coal resources were selected for the project: the West Kentucky No. 13 (Baker) Seam and an undisclosed lignite resource in the Illinois coal basin. The estimated quantities of REEs in these resources were calculated based on production samples and drilling data. It was estimated that there is adequate supply for an operation producing one metric ton daily of higher purity mixed rare earth oxides (MREO) for approximately 20 years at a site located in western Kentucky. In Circuit 1, project data was obtained from a pilot heap leach and REE concentration facility. It was concluded that the existing circuit, which generated a MREO concentrate, two types of CM mixed products, and Li- and Sr-containing waters, would be suitable feed for circuits 2-4. Data from the first-of-its-kind coal coarse refuse heap leach pilot pad played a crucial role in estimating reliable elemental concentrations of the pregnant leaching solution (PLS). The average total REE concentration in the PLS was found to be 28.6 ppm. In Circuit 2, several concepts were explored including a novel process referred to as solvent-assisted chromatography (SAC). This concept involved a novel columnar reactor that incorporated multiple mixer/settlers, thereby enabling the operation of counter-flowing aqueous and organic phases. Unfortunately, due to project time constraints, a complete fundamental modeling analysis could not be completed to fully evaluate the technology. Molten salt electrowinning was considered as an alternative for circuit 3 following circuit 2 purification circuit utilizing the novel SAC process. A mass and energy balance of Nd reduction to metal in a fluoride containing molten salt electrolyte was conducted. Comparisons were made with the current state of Asian molten salt electrorefining, and potential improvements in siphoning rare earth metals (REM) from the reactor were presented. A cost estimate was performed for the production of 1 tonne per day, which yielded a total of $2.29 million for the nine electrowinning (EW) cells required. The selected option for circuits 2 and 3 was a plasma distillation process, which initially separates rare earth elements (REEs) from other elements. This is followed by selective electrowinning in various ionic liquids. The selection was made on the basis of thermodynamic modeling and experimental data previously published by a project partner. The combination offers an innovative approach to integrated refining and RE metal production. For Circuit 4, an extensive literature review was conducted for the processing of the CMs. The ultimate decision was to utilize a combined plasma and ionic liquid process as well to produce individual high-purity concentrates of Zn, Ni, Co, Mn, and Mg. A separate flowsheet for Li and Sr was recommended, which would yield carbonates of these elements. Due to the lack of suitable experimental data at this time, a process recommendation could not be provided but several methods have been proposed for consideration. Lastly, a techno-economic analysis (TEA) was conducted to assess the effectiveness of the proposed process for further investigation. The TEA results revealed a capital expense (CapEx) of $737 million and an annual operational expense (OpEx) of $220 million. Due to the selected elements, the hypothetical heap leach pad can produce 1 metric tonne per day of REO equivalent, but a conscious decision was made to only treat targeted REEs, resulting in the production of 0.4 metric tonne of REM. An estimated annual revenue of $90.87 million was projected based on standard market pricing information provided by the funding agency. During the TEA, ten different modules were evaluated for costing purposes. The precipitation circuit was identified as the largest single operational expense, followed by the Mg/Mn process due to the amount of treated metal. In terms of capital expenditures, the heap leach process incurred the highest cost, followed by the Mg/Mn process. The scalability of the plasma process is a crucial consideration since the reactors cannot be scaled beyond the largest demonstrated size due to their reliance on surface area of the slag and vapor phase. The purity estimate for the REEs are generally 98%±2% to produce a metal. The purity level being lower than the project objective was due to the lack of specific experimental data needed to tighten the tolerance of the estimates. Based on literature and previous experience, the CMs are estimated as follows; Ga (95%+, metal), Sr (95%+, carbonate), Li (95%+, carbonate), Ni (98%±2%, metal), Zn (95%+, metal sponge), Ge (95%+, metal), Co (98%±2%, metal), and Mn (98%±2%, metal).

01 COAL, LIGNITE, AND PEAT↗