Search NASA⌕ Search

SEARCH · Search NASA

Results for “Spatial refinement”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Analysis of Hydrology and Interim Measure Performance in the Chromium Plume at Los Alamos National Laboratory - 20498

Hexavalent chromium, Cr(VI), is present in the regional aquifer beneath Los Alamos National Laboratory (LANL) at concentrations greater than the New Mexico groundwater standard (50 ppb). The source of the Cr is blowdown from a power plant where potassium dichromate was used as a corrosion inhibitor from 1956 to 1972. Blowdown comingled with effluent from another source and migrated approximately 3 miles down a canyon, primarily as surface flow, until reaching an alluvial infiltration zone. The infiltrating water moved vertically through a thick (approximately 900 ft), complex stratigraphic sequence of unsaturated tuffs, basalt, and unconsolidated alluvial deposits. Horizontal transport also took place within spatially limited perched groundwater horizons in the vadose zone (VZ). Breakthrough at the regional aquifer water table is thought to have occurred in multiple hydraulic windows or 'drip points' located approximately 3 miles from the initial release site at the power plant. The regional aquifer is highly heterogeneous and anisotropic, and displays unconfined behavior in an upper zone and confined behavior at depth. The Cr(VI) contamination has so far remained confined to the upper approximately 25 m of the aquifer. The ambient hydraulic gradient in the regional aquifer is toward the east/southeast, but the water table is extremely flat in the Cr plume area compared to surrounding areas. There is some local variability in flow direction, potentially related to effects of seasonal pumping from nearby Los Alamos county water-supply wells, which draw from the regional aquifer several hundred feet below the water table. An interim measure (IM) is currently operating at the site to achieve hydraulic control of plume migration using a series of extraction wells, ion-exchange treatment, and injection along the downgradient portion of the plume. The migration of Cr contamination in the aquifer has been previously evaluated using a calibrated probabilistic numerical model, which has also been used to design the spatial configuration of extraction and injection wells for the IM, and to identify locations for monitoring wells. The model is built using the Finite Element Heat and Mass (FEHM) transfer code, and the Model Analysis and Decision Support (MADS) software is used for analyses. To date, the portion of the IM system that has operated in a continuous manner has been very effective at hydraulic control. Modifications to the flow field associated with the influence of IM pumping and injection are evaluated here, with modeling constraints to be provided by recent observations of tracers injected into both monitoring wells and IM injection wells, as well as opportunistic observations of geochemical signatures of treated ion-exchange effluent injected into IM injection wells. The tracer and geochemical data provide insights into modified flow directions within the Cr plume area that cannot be otherwise obtained. Results of the modeling work are being used to inform adaptive management of IM operations and will be a valuable tool for development and refinement of long-term remediation strategies. (authors)

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

Combining Observations and Models: A Review of the CARDAMOM Framework for Data‐Constrained Terrestrial Ecosystem Modeling

The rapid increase in the volume and variety of terrestrial biosphere observations (i.e., remote sensing data and in situ measurements) offers a unique opportunity to derive ecological insights, refine process‐based models, and improve forecasting for decision support. However, despite their potential, ecological observations have primarily been used to benchmark process‐based models, as many past and current models lack the capability to directly integrate observations and their associated uncertainties for parameterization. In contrast, data assimilation frameworks such as the CARbon DAta MOdel fraMework (CARDAMOM) and its suite of process‐based models, known as the Data Assimilation Linked Ecosystem Carbon Model (DALEC), are specifically designed for model‐data fusion. This review, motivated by a recent CARDAMOM community workshop, examines the development and applications of CARDAMOM, with an emphasis on its role in advancing ecosystem process understanding. CARDAMOM employs a Bayesian approach, using a Markov Chain Monte Carlo algorithm to enable data‐driven calibration of DALEC parameters and initial states (i.e., carbon pool sizes) through observation operators. CARDAMOM's unique ability to retrieve localized model process parameters from diverse datasets—ranging from in situ measurements to global satellite observations—makes it a highly flexible tool for analyzing spatially variable ecosystem responses to environmental change. However, assimilating these data also presents challenges, including data quality issues that propagate into model skill, as well as trade‐offs between model complexity, parameter equifinality, and predictive performance. We discuss potential solutions to these challenges, such as reducing parameter equifinality by incorporating new observations. This review also offers community recommendations for incorporating emerging datasets, integrating machine learning techniques, strengthening collaboration with remote sensing, field, and modeling communities, and expanding CARDAMOM's relevance for localized ecosystem monitoring and decision‐making. CARDAMOM enables a deep, mechanistic understanding of terrestrial ecosystem dynamics that cannot be achieved through empirical analyses of observational datasets or weakly constrained models alone.

Bayesian inference↗
NASA NTRS1998

A 2-D Interface Element for Coupled Analysis of Independently Modeled 3-D Finite Element Subdomains

Over the past few years, the development of the interface technology has provided an analysis framework for embedding detailed finite element models within finite element models which are less refined. This development has enabled the use of cascading substructure domains without the constraint of coincident nodes along substructure boundaries. The approach used for the interface element is based on an alternate variational principle often used in deriving hybrid finite elements. The resulting system of equations exhibits a high degree of sparsity but gives rise to a non-positive definite system which causes difficulties with many of the equation solvers in general-purpose finite element codes. Hence the global system of equations is generally solved using, a decomposition procedure with pivoting. The research reported to-date for the interface element includes the one-dimensional line interface element and two-dimensional surface interface element. Several large-scale simulations, including geometrically nonlinear problems, have been reported using the one-dimensional interface element technology; however, only limited applications are available for the surface interface element. In the applications reported to-date, the geometry of the interfaced domains exactly match each other even though the spatial discretization within each domain may be different. As such, the spatial modeling of each domain, the interface elements and the assembled system is still laborious. The present research is focused on developing a rapid modeling procedure based on a parametric interface representation of independently defined subdomains which are also independently discretized.

Kandil, Osama A.↗
NASA NTRS2004

[Analysis of Multiplatform CO (Carbon Monoxide) Measurements During Trace-P Mission]

Carbon monoxide is considered mission critical (TRACE-P NRA) because it is one of the gases involved in controlling the oxidizing power of the atmosphere and, as a tracer gas, is valuable in interpreting mission data sets. Carbon monoxide exhibits interannual differences, suggesting relatively short-term imbalances in sources and sinks. Sources of CO are dominated by fossil fuel combustion, biomass burning, and the photochemical oxidation of CH4 and nonmethane hydrocarbons while reaction with OH is believed to be the major sink for atmospheric CO, with additional losses due to soil uptake. Uncertainties in the magnitude and distribution of both sources and sinks remain fairly large however, and additional data are required to refine the global budget. Seasonal changes and a northern hemispheric latitudinal gradient have been described for a variety of Pacific basin sites through long-term monitoring of surface background levels. Latitudinal variations have also recently been described at upper tropospheric altitudes over a multi-year period by. TRACE-P will provide an aircraft survey of CO over the northern Pacific in the northern spring when CO concentrations are at their seasonal maximum in the northern hemisphere (NH) and at their seasonal minimum in the southern hemisphere (SH). Previous GTE missions, Le., PEM West-B and PEM Tropics-B, ground-based, and satellite observations (MAPS, April 1994) give us a general picture of the distribution of CO over the northern Pacific during this season. Based on these measurements, background CO levels over remote ocean areas are anticipated to be in the range of 110 - 180 ppbv, while those closer to the Asian continent may rise as high as 600 ppbv. These measurements also reveal high spatial variability (both horizontal and vertical) as well as temporal variations in CO over the area planned for the TRACE-P mission. This variability is a result of multiple CO sources, the meteorological complexity of transport processes, and the photochemical aging of air masses. The influence of biomass burning in the southern Pacific should be relatively small since the mission coincides with the southern tropical wet season when agricultural burning is at its seasonal low. The proposed CO measurements taken during TRACE-P should therefore largely be a function of the impact of various NH sources, primarily Asian and predominantly fossil fuel combustion and biomass burning. These processes are also major sources of many other atmospheric pollutants, consequently making accurate and precise CO measurements is one of the highest TRACE-P priorities [TRACE-P NRA]. The TRACE-P mission emphasizes the dual objectives of assessing the magnitude of the transport of chemically and radiatively important gases such as CO from Asia to the western Pacific, and determining how emissions change and are modified during this transport.

Pougatchev, Nikita S.↗
NASA NTRS2019

Acoustic Phased Array Quantification of Quiet Technology Demonstrator 3 Advanced Inlet Liner Noise Component

Acoustic phased array flyover noise measurements were acquired as part of the Boeing 737 MAX-7 NASA Advanced Inlet Liner segment of the Quiet Technology Demonstrator 3 (QTD3) flight test program. This paper reports on the processes used for separating and quantifying the engine inlet, exhaust and airframe noise source components and provides sample phased array-based comparisons of the component noise source levels associated with the inlet liner treatment configurations. Full scale flyover noise testing of NASA advanced inlet liners was conducted as part of the Quiet Technology Demonstrator 3 flight test program in July and August of 2018. Details on the inlet designs and testing are provided in the companion paper of Reference 1. The present paper provides supplemental details relating to the acoustic phased array portion of the analyses provided in Ref. 1. In brief, the test article was a Boeing 737MAX-7 aircraft with a modified right hand (starboard side) engine inlet, which consisted of either a production inlet liner, a NASA designed inlet liner or a simulated hard wall configuration (accomplished by applying speed tape over the inlet acoustic treatment areas). In all three configurations, the engine forward fan case acoustic panel was replaced with a unperforated (hardwall) panel. No other modifications to any other acoustic treatment areas were made. The left hand (port side) engine was a production engine and was flown at idle thrust for all measurements in order to isolate the effects of the inlet liners to the right hand engine. As described in Ref. 1, the NASA inlet treatment consists of laterally cut slots (cut perpendicular to the flow direction) which are designed to reduce excrescence drag while maintaining or exceeding the liner acoustic noise reduction capabilities. The NASA inlet liner consists of a Multi-Degree of Freedom (MDOF) design with two breathable septum layers inserted into each honeycomb cell [1]. The aircraft noise measurements were acquired for both takeoff (flaps 1 setting, gear up) and approach (flaps 30 gear up and gear down) configurations. The inlet and flight test configurations are summarized in Table 1. Table 1: Inlet Treatment and Flight Configurations Inlet Forward Fan Case Aircraft Production Hardwall Flaps 1, gear up; flaps 30 gear up; flaps 30 gear down NASA Hardwall Flaps 1, gear up; flaps 30 gear up; flaps 30 gear down Hardwall Hardwall Flaps 1, gear up; flaps 30 gear up; flaps 30 gear down III.Test Description and Hardware The flight testing was conducted at the Grant County airport in Moses Lake, WA, between 27 July and 6 August 2018. The noise measurement instrumentation included 8 flush dish microphones arranged in a noise certification configuration as well as an 840 microphone phased array. The flush dish microphones were used to quantify the levels and differences in levels between the various inlet treatments. The phased array was used to separate and quantify the narrowband (tonal) and broadband noise component levels from the engine inlet/exhaust and from the airframe. Phased array extraction of the broadband component was critical to this study because it allowed for the separation of the inlet component from the total airplane level noise even when it was significantly below the total level. Figure 1 provides an overview of the phased array microphone layout as well as a detailed image of an individual phased array microphone mounted in a plate holder (the microphone sensor is the dot in the center of the plate). The ground plane ensemble array microphones (referred to as “ensemble array” in this paper) were mounted in plates with “flower petal” edges designed to minimize edge scattering effects. Fig. 1 Flyover test microphone layout. The phased array configuration was the result of a progressive development of concepts originally implemented in Ref. 2 and refined over the following years, consisting namely of multiple multi-arm logarithmic spiral subarrays designed to cover overlapping frequency ranges and optimized for various aircraft emission angles. For the present case, the signals from all 840 microphones were acquired on a single system. The 840 microphones were parsed into 11 primary subarray sets spanning from smallest to largest aperture size and labeled accordingly as “a, b, …, k”, where “a” corresponds to the smallest fielded subarray and “k” corresponds to the largest aperture subarray. The apertures ranged from approximately 10 ft to 427 ft in size (in the flight direction) with the subarrays consisting of between 215 and 312 microphones. Figure 2 shows three such subarrays, k, h and a. As done in Ref. 2, microphones were shared between subarrays in order to reduce total channel count. Fig. 2 Sample subarray sizes (20° from overhead – refer to Figure 3a discussion). In addition to the above, each of the 11 primary subarray sets consisted of four subarrays optimized to provide near equivalent array spatial resolution in both the flight and lateral directions within 30 degrees of overhead (i.e., airplane directly above the center of the array), namely, at angles  of 0, ±10, ±20 and ±30 degrees relative to overhead where angle  is defined as shown in Figure 3a. This allowed for optimized aircraft noise measurements from 60 to 120 degree emission angle.6 An example of this pletharray design is shown in Figure 3b for the k subarray. When the aircraft is at overhead, the microphones indicated by the blue markers are used for beamforming. When the aircraft is at angles ±10 degrees from overhead, both the blue and red colored microphones are used, and so on for the ±20 and ±30 degree aircraft locations. See Ref. 3 for extensive details on pletharray design for aeroacoustic phased array testing. 6 In the discussions that follow, emission angle values are used. These are the angles at the time sound is emitted relative to the engine axis and are calculated based on flight path angle, body aircraft body angle with respect to the relative wind direction, and engine axis angle relative to aircraft body angle.

Brusniak, Leon↗
NASA NTRS1994

Design of a composite filter realizable on practical spatial light modulators

Hybrid optical correlator systems use two spatial light modulators (SLM's), one at the input plane and the other at the filter plane. Currently available SLM's such as the deformable mirror device (DMD) and liquid crystal television (LCTV) SLM's exhibit arbitrarily constrained operating characteristics. The pattern recognition filters designed with the assumption that the SLM's have ideal operating characteristic may not behave as expected when implemented on the DMD or LCTV SLM's. Therefore it is necessary to incorporate the SLM constraints in the design of the filters. In this report, an iterative method is developed for the design of an unconstrained minimum average correlation energy (MACE) filter. Then using this algorithm a new approach for the design of a SLM constrained distortion invariant filter in the presence of input SLM is developed. Two different optimization algorithms are used to maximize the objective function during filter synthesis, one based on the simplex method and the other based on the Hooke and Jeeves method. Also, the simulated annealing based filter design algorithm proposed by Khan and Rajan is refined and improved. The performance of the filter is evaluated in terms of its recognition/discrimination capabilities using computer simulations and the results are compared with a simulated annealing optimization based MACE filter. The filters are designed for different LCTV SLM's operating characteristics and the correlation responses are compared. The distortion tolerance and the false class image discrimination qualities of the filter are comparable to those of the simulated annealing based filter but the new filter design takes about 1/6 of the computer time taken by the simulated annealing filter design.

Rajan, P. K.↗
DOE OSTI2021

Ordering effects on deformation substructures and strain hardening behavior of a CrCoNi based medium entropy alloy

Here, a CrCoNi based medium entropy alloy with small additions of Ti, Al and Nb (denoted as (CrCoNi) 93 Al 4 Ti 2 Nb) in the as-quenched condition, exhibits tensile properties comparable to those of the equiatomic CrCoNi alloy at room temperature. Dark field transmission electron microscopy (TEM), atomic resolution high-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM) together with atom probe tomography (APT) show that spatially-localized long range ordering (LRO) L1 2 domains exist in this alloy. The evolution of deformation substructure with plastic deformation in this alloy was characterized using electron backscatter diffraction (EBSD), electron channeling contrast imaging (ECCI) and STEM based techniques including the recently developed weak beam dark field STEM imaging. Plastic deformation occurs by the slip of a/2<110>dislocations, which are narrowly dissociated into Shockley partial dislocations on {111} slip planes. Their dissociation distances in the (CrCoNi) 93 Al 4 Ti 2 Nb alloy are much smaller than the widths of the corresponding partials in the equiatomic CrCoNi alloy due to one or more of the minor alloying elements (Al, Ti, Nb). Dislocation slip in this alloy has a pronounced planar character. The leading dislocations in slip bands glide as pairs due to the existence of LRO domains. Multipoles were formed through the slip of dislocations with opposite signs on adjacent {111} slip planes. Those multipoles serve as building blocks for the formation of subgrain structures consisting of fine slip bands. The distances between slip bands were continuously refined during plastic deformation and dynamic refinement of slip bands plays a crucial role in strain hardening. The effects of LRO domains on planar dislocation slip, the deactivation of deformation twinning and strain hardening of this alloy are discussed.

36 MATERIALS SCIENCE↗
NASA NTRS2016

Soil Moisture Active Passive Mission L4_SM Data Product Assessment (Version 2 Validated Release)

During the post-launch SMAP calibration and validation (Cal/Val) phase there are two objectives for each science data product team: 1) calibrate, verify, and improve the performance of the science algorithm, and 2) validate the accuracy of the science data product as specified in the science requirements and according to the Cal/Val schedule. This report provides an assessment of the SMAP Level 4 Surface and Root Zone Soil Moisture Passive (L4_SM) product specifically for the product's public Version 2 validated release scheduled for 29 April 2016. The assessment of the Version 2 L4_SM data product includes comparisons of SMAP L4_SM soil moisture estimates with in situ soil moisture observations from core validation sites and sparse networks. The assessment further includes a global evaluation of the internal diagnostics from the ensemble-based data assimilation system that is used to generate the L4_SM product. This evaluation focuses on the statistics of the observation-minus-forecast (O-F) residuals and the analysis increments. Together, the core validation site comparisons and the statistics of the assimilation diagnostics are considered primary validation methodologies for the L4_SM product. Comparisons against in situ measurements from regional-scale sparse networks are considered a secondary validation methodology because such in situ measurements are subject to up-scaling errors from the point-scale to the grid cell scale of the data product. Based on the limited set of core validation sites, the wide geographic range of the sparse network sites, and the global assessment of the assimilation diagnostics, the assessment presented here meets the criteria established by the Committee on Earth Observing Satellites for Stage 2 validation and supports the validated release of the data. An analysis of the time average surface and root zone soil moisture shows that the global pattern of arid and humid regions are captured by the L4_SM estimates. Results from the core validation site comparisons indicate that "Version 2" of the L4_SM data product meets the self-imposed L4_SM accuracy requirement, which is formulated in terms of the ubRMSE: the RMSE (Root Mean Square Error) after removal of the long-term mean difference. The overall ubRMSE of the 3-hourly L4_SM surface soil moisture at the 9 km scale is 0.035 cubic meters per cubic meter requirement. The corresponding ubRMSE for L4_SM root zone soil moisture is 0.024 cubic meters per cubic meter requirement. Both of these metrics are comfortably below the 0.04 cubic meters per cubic meter requirement. The L4_SM estimates are an improvement over estimates from a model-only SMAP Nature Run version 4 (NRv4), which demonstrates the beneficial impact of the SMAP brightness temperature data. L4_SM surface soil moisture estimates are consistently more skillful than NRv4 estimates, although not by a statistically significant margin. The lack of statistical significance is not surprising given the limited data record available to date. Root zone soil moisture estimates from L4_SM and NRv4 have similar skill. Results from comparisons of the L4_SM product to in situ measurements from nearly 400 sparse network sites corroborate the core validation site results. The instantaneous soil moisture and soil temperature analysis increments are within a reasonable range and result in spatially smooth soil moisture analyses. The O-F residuals exhibit only small biases on the order of 1-3 degrees Kelvin between the (re-scaled) SMAP brightness temperature observations and the L4_SM model forecast, which indicates that the assimilation system is largely unbiased. The spatially averaged time series standard deviation of the O-F residuals is 5.9 degrees Kelvin, which reduces to 4.0 degrees Kelvin for the observation-minus-analysis (O-A) residuals, reflecting the impact of the SMAP observations on the L4_SM system. Averaged globally, the time series standard deviation of the normalized O-F residuals is close to unity, which would suggest that the magnitude of the modeled errors approximately reflects that of the actual errors. The assessment report also notes several limitations of the "Version 2" L4_SM data product and science algorithm calibration that will be addressed in future releases. Regionally, the time series standard deviation of the normalized O-F residuals deviates considerably from unity, which indicates that the L4_SM assimilation algorithm either over- or under-estimates the actual errors that are present in the system. Planned improvements include revised land model parameters, revised error parameters for the land model and the assimilated SMAP observations, and revised surface meteorological forcing data for the operational period and underlying climatological data. Moreover, a refined analysis of the impact of SMAP observations will be facilitated by the construction of additional variants of the model-only reference data. Nevertheless, the “Version 2” validated release of the L4_SM product is sufficiently mature and of adequate quality for distribution to and use by the larger science and application communities.

SMAP L4_SM↗
NASA NTRS2018

Spatial and Seasonal Variations in C3Hx Hydrocarbon Abundance in Titan’s Stratosphere from Cassini CIRS Observations

Of the C3Hx hydrocarbons, propane (C3H8) and propyne (methylacetylene, CH3C2H) were first detected in Titan’s atmosphere during the Voyager 1 flyby in 1980. Propene (propylene, C3H6) was first detected in 2013 with data from the Composite InfraRed Spectrometer (CIRS) instrument on Cassini. We present the first measured abundance profiles of propene on Titan from radiative transfer modeling, and compare our measurements to predictions derived from several photochemical models. Near the equator, propene is observed to have a peak abundance of 10 ppbv at a pressure of 0.2 mbar. Several photochemical models predict the amount at this pressure to be in the range 0.3 - 1 ppbv and also show a local minimum near 0.2 mbar which we do not see in our measurements. We also see that propene follows a different latitudinal trend than the other C3 molecules. While propane and propyne concentrate near the winter pole, transported via a global convective cell, propene is most abundant above the equator. We retrieve vertical abundances profiles between 125 km and 375 km for these gases for latitude averages between 60◦S to 20◦S, 20◦S to 20◦N, and 20◦N to 60◦N over two time periods, 2004 through 2009 representing Titan’s atmosphere before the 2009 equinox, and 2012 through 2015 representing time after the equinox. Additionally, using newly corrected line data, we determined an updated upper limit for allene (propadiene, CH2CCH2), the isomer of propyne). The measurements we present will further constrain photochemical models, by refining reaction rates and the transport of these gases throughout Titan’s atmosphere.

Nicholas A. Lombardo↗
NASA NTRS2021

Casual Rerouting of AERONET Sun/Sky Photometers: Toward a New Network of Ground Measurements Dedicated to the Monitoring of Surface Properties?

This paper presents an innovative method for observing vegetation health at a very high spatial resolution (~5 × 5 cm) and low cost by upgrading an existing Aerosol RObotic NETwork (AERONET) ground station dedicated to the observation of aerosols in the atmosphere. This study evaluates the capability of a sun/sky photometer to perform additional surface reflectance observations. The ground station of Toulouse, France, which belongs to the AERONET sun/sky photometer network, is used for this feasibility study. The experiment was conducted for a 5-year period (between 2016 and 2020). The sun/sky photometer was mounted on a metallic structure at a height of 2.5 m, and the acquisition software was adapted to add a periodical (every hour) ground-observation scenario with the sun/sky photometer observing the surface instead of being inactive. Evaluation is performed by using a classical metric characterizing the vegetation health: the normalized difference vegetation index (NDVI), using as reference the satellite NDVI derived from a Sentinel-2 (S2) sensor at 10 × 10 m resolution. Comparison for the 5-year period showed good agreement between the S2 and sun/sky photometer NDVIs (i.e., bias = 0.004, RMSD = 0.082, and R = 0.882 for a mean value of S2A NDVI around 0.6). Discrepancies could have been due to spatial-representativeness issues (of the ground measurement compared to S2), the differences between spectral bands, and the quality of the atmospheric correction applied on S2 data (accuracy of the sun/sky photometer instrument was better than 0.1%). However, the accuracy of the atmospheric correction applied on S2 data in this station appeared to be of good quality, and no dependence on the presence of aerosols was observed. This first analysis of the potential of the CIMEL CE318 sun/sky photometer to monitor the surface is encouraging. Further analyses need to be carried out to estimate the potential in different AERONET stations. The occasional rerouting of AERONET stations could lead to a complementary network of surface reflectance observations. This would require an update of the software, and eventual adaptations of the measurement platforms to the station environments. The additional cost, based on the existing AERONET network, would be quite limited. These new surface measurements would be interesting for measurements of vegetation health (monitoring of NDVI, and also of other vegetation indices such as the leaf area and chlorophyll indices), for validation and calibration exercise purposes, and possibly to refine various scientific algorithms (i.e., algorithms dedicated to cloud detection or the AERONET aerosol retrieval algorithm itself). CIMEL is ready to include the ground scenario used in this study in all new sun/sky photometers.

Dominique Carrer↗
DOE OSTI2023

Characterization of Rare-Earth Elements in Lignite Coal of the Williston Basin: Past Efforts and Ongoing Work

Rare-earth elements (REEs) have been a subject area of high interest for their unique properties. REEs are crucial materials used in an incredible array of consumer goods, energy system components, and military defense applications. While the United States has one operating REE mine, the product is sent overseas for refining into usable metals making the United States 100% import-reliant on these critical materials. This has led the Federal Government to declare the REE market an issue of national security. The Energy and Environmental Research Center (EERC) under funding provided by the Department of Energy (DOE) has undertaken efforts to determine if the lignite coal found in the Williston Basin has the potential to be an ore body containing sufficient quantities of REEs and Critical Minerals (CM) for extraction and processing. In one such effort, the EERC collected over 400 samples from the Williston Basin lignite coal seams including outcrops as well as active mines. Those efforts have been followed up with ongoing work under the U.S. DOE’s Carbon Ore, Rare Earth and Critical Minerals Initiative (CORE-CM) currently ongoing in the Williston Basin as well as other basins in the United States. This ongoing work in the Williston Basin characterizing REEs in coal has focused on the collection of new sampling and analysis of REEs in coal and building upon previous characterization work. This information is being used to understand the spatial distribution of REEs through mapping and 3D modeling. The goals of these efforts are to better understand the mechanisms for distribution of the REEs in lignite coal as well as their concentration and determine knowledge gaps in characterization and to start to build the database required to understand the potential resource in the Williston Basin.

Feole, Ian K.↗
NASA NTRS1998

Theoretical Accuracy of Global Snow-Cover Mapping Using Satellite Data in the Earth Observing System (EOS) Era

Following the launch of the Earth Observing System first morning (EOS-AM1) satellite, daily, global snow-cover mapping will be performed automatically at a spatial resolution of 500 m, cloud-cover permitting, using Moderate Resolution Imaging Spectroradiometer (MODIS) data. A technique to calculate theoretical accuracy of the MODIS-derived snow maps is presented. Field studies demonstrate that under cloud-free conditions when snow cover is complete, snow-mapping errors are small (less than 1%) in all land covers studied except forests where errors are greater and more variable. The theoretical accuracy of MODIS snow-cover maps is largely determined by percent forest cover north of the snowline. Using the 17-class International Geosphere-Biosphere Program (IGBP) land-cover maps of North America and Eurasia, the Northern Hemisphere is classified into seven land-cover classes and water. Snow-mapping errors estimated for each of the seven land-cover classes are extrapolated to the entire Northern Hemisphere for areas north of the average continental snowline for each month. Average monthly errors for the Northern Hemisphere are expected to range from 5 - 10%, and the theoretical accuracy of the future global snow-cover maps is 92% or higher. Error estimates will be refined after the first full year that MODIS data are available.

Hall, D. K.↗
DOE OSTI2022

A weighted state redistribution algorithm for embedded boundary grids

State redistribution is an algorithm that stabilizes cut cells for embedded boundary grid methods. This work extends the earlier algorithm in several important ways. First, state redistribution is extended to three spatial dimensions. Second, we discuss several algorithmic changes and improvements motivated by the more complicated cut cell geometries that can occur in higher dimensions. In particular, we introduce a weighted version with less dissipation in an easily generalizable framework. Third, we demonstrate that state redistribution can also stabilize a solution update that includes both advective and diffusive contributions. Notably, the stabilization algorithm is shown to be effective for incompressible as well as compressible reacting flows. Finally, we discuss the implementation of the algorithm for several exascale-ready simulation codes based on AMReX, demonstrating ease of use in combination with domain decomposition, hybrid parallelism and complex physics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A Computational Framework to design 3D stiffness gradient acoustic metamaterials for impedance matching

Acoustic waves play a crucial role in various applications, including medical imaging, non-destructive testing, and sonar systems. One of the significant challenges in these applications is impedance matching, which is essential for minimizing reflections and maximizing the transfer of acoustic energy between different media. Acoustic metamaterials offer a promising solution to this challenge. In addition to impedance control, gradient stiffness can enhance structural efficiency and enable spatial control of wave propagation, making it a valuable feature in acoustic metamaterial design. In this pa- per, we present our developed computational method to design 3D stiffness gradient acoustic metamaterials for impedance matching. The key steps in our approach include generating initial designs using a periodic covariance function to provide unit cells that are both periodic on the boundaries and randomly formed inside the unit cell. Furthermore, we integrated manufacturing constraints into the design process, ensuring that the structures are interconnected for fabrication. We propose two computational optimization algorithms: GenUnit, based on a non-dominated sorting genetic algorithm (NSGA-II), and MLMatch, which leverages differentiable machine learning. The two approaches are not separate contributions but complementary com- ponents of a unified framework. GenUnit requires no training data and directly interfaces with physics-based simulations, making it highly accurate but slower for large-scale exploration. In contrast, MLMatch is data-hungry during training but, once trained, enables near-instantaneous inference and broad design-space coverage. Together, they form a hybrid strategy: ML- Match rapidly explores the global design space, and GenUnit provides local refinement with high-fidelity accuracy. This balance between training cost, inference time, and precision is the motivation for including both methods in the same study. We applied this dual-algorithm framework to generate two metallic-based metamaterial designs that match the acoustic impedance of water while exhibiting a controlled gradient in stiffness (from stiff to soft). The stiffness gradient is particularly advantageous in applications where one side of the structure must interface with soft or sensitive surfaces, such as human tissue or delicate components. Here, this work paves the way for improved materials in various acoustic applications, particularly in ultrasound devices, by providing better impedance.

Metamaterial↗
DOE OSTI2021

A New Martian Crater Chronology: Implications for Jezero Crater

Crater chronologies are a fundamental tool to assess the relative and absolute ages of planetary surfaces when direct radiometric dating is not available. Martian crater chronologies are derived from lunar crater spatial densities on terrains with known radiometric ages, and thus they critically depend on the Moon-to-Mars extrapolation. This extrapolation requires knowledge of the time evolution of the impact flux, including contributions from various impactor populations, factors that are not trivially connected to the dynamical evolution of the early Solar System. In this paper, we will present a new Martian crater chronology based on current dynamical models, and consider the main sources of uncertainties (e.g., impactor size–frequency distribution; dynamical models with late and early instabilities, etc.). The resulting “envelope” of Martian crater chronologies significantly differs from previous chronologies. The new Martian crater chronology is discussed using two interesting applications: Jezero crater’s dark terrain (relevant to the NASA Mars 2020 mission) and the southern heavily cratered highlands. Our results indicate that Jezero’s dark terrain may have formed ∼3.1 Ga, i.e., up to 0.5 Gyr older than previously thought. In addition, available crater chronologies (including our own) overestimate the number of craters larger than 150 km on the southern highlands, suggesting either that large craters have been efficiently erased over Martian history or that dynamical models need further refinement. Further, our chronology constrains the age of Isidis basin to be 4.05–4.2 Ga and that of the Borealis basin to be 4.35–4.40 Ga; these are predictions that can be tested with future sample and return missions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Boosting background suppression in the NEXT experiment through Richardson-Lucy deconvolution

Next-generation neutrinoless double beta decay experiments aim for half-life sensitivities of ~ 10 27 yr, requiring suppressing backgrounds to < 1 count/tonne/yr. For this, any extra background rejection handle, beyond excellent energy resolution and the use of extremely radiopure materials, is of utmost importance. The NEXT experiment exploits differences in the spatial ionization patterns of double beta decay and single-electron events to discriminate signal from background. While the former display two Bragg peak dense ionization regions at the opposite ends of the track, the latter typically have only one such feature. Thus, comparing the energies at the track extremes provides an additional rejection tool. The unique combination of the topology-based background discrimination and excellent energy resolution (1% FWHM at the Q-value of the decay) is the distinguishing feature of NEXT. Previous studies demonstrated a topological background rejection factor of ~ 5 when reconstructing electron-positron pairs in the 208Tl 1.6 MeV double escape peak (with Compton events as background), recorded in the NEXT-White demonstrator at the Laboratorio Subterráneo de Canfranc, with 72% signal efficiency. This was recently improved through the use of a deep convolutional neural network to yield a background rejection factor of ~ 10 with 65% signal efficiency. Here, we present a new reconstruction method, based on the Richardson-Lucy deconvolution algorithm, which allows reversing the blurring induced by electron diffusion and electroluminescence light production in the NEXT TPC. The new method yields highly refined 3D images of reconstructed events, and, as a result, significantly improves the topological background discrimination. When applied to real-data 1.6 MeV e - e + pairs, it leads to a background rejection factor of 27 at 57% signal efficiency.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS2018

Spatial and Seasonal Variations in C3Hx Hydrocarbon Abundance in Titan's Stratosphere from Cassini CIRS Observations

Of the C3Hx hydrocarbons, propane (C3H8) and propyne (methylacetylene, CH3C2H) were first detected in Titan’s atmosphere during the Voyager 1 flyby in 1980. Propene (propylene, C3H6) was first detected in 2013 with data from the Composite InfraRed Spectrometer (CIRS) instrument on Cassini. We present the first measured abundance profiles of propene on Titan from radiative transfer modeling, and compare our measurements to predictions derived from several photochemical models. Near the equator, propene is observed to have a peak abundance of 10 ppbv at a pressure of 0.2 mbar. Several photochemical models predict the amount at this pressure to be in the range 0.3–1 ppbv and also show a local minimum near 0.2 mbar which we do not see in our measurements. We also see that propene follows a different latitudinal trend than the other C3 molecules. While propane and propyne concentrate near the winter pole, transported via a global convective cell, propene is most abundant above the equator. We retrieve vertical abundances profiles between 125 km and 375 km for these gases for latitude averages between 60°S–20°S, 20°S–20°N, and 20°N–60°N over two time periods, 2004 through 2009 representing Titan’s atmosphere before the 2009 equinox, and 2012 through 2015 representing time after the equinox. Additionally, using newly corrected line data, we determined an updated upper limit for allene (propadiene, CH2CCH2, the isomer of propyne). We claim a 3-σ upper limit mixing ratio of 2.5 × 10(exp −9) within 30° of the equator. The measurements we present will further constrain photochemical models by refining reaction rates and the transport of these gases throughout Titan’s atmosphere.

Lombardo, Nicholas A.↗
NASA NTRS1994

Experiments with the Mesoscale Atmospheric Simulation System (MASS) using the synthetic relative humidity

This study is intended to examine the impact of the synthetic relative humidity on the model simulation of mesoscale convective storm environment. The synthetic relative humidity is derived from the National Weather Services surface observations, and non-conventional sources including aircraft, radar, and satellite observations. The latter sources provide the mesoscale data of very high spatial and temporal resolution. The synthetic humidity data is used to complement the National Weather Services rawinsonde observations. It is believed that a realistic representation of initial moisture field in a mesoscale model is critical for the model simulation of thunderstorm development, and the formation of non-convective clouds as well as their effects on the surface energy budget. The impact will be investigated based on a real-data case study using the mesoscale atmospheric simulation system developed by Mesoscale Environmental Simulations Operations, Inc. The mesoscale atmospheric simulation system consists of objective analysis and initialization codes, and the coarse-mesh and fine-mesh dynamic prediction models. Both models are a three dimensional, primitive equation model containing the essential moist physics for simulating and forecasting mesoscale convective processes in the atmosphere. The modeling system is currently implemented at the Applied Meteorology Unit, Kennedy Space Center. Two procedures involving the synthetic relative humidity to define the model initial moisture fields are considered. It is proposed to perform several short-range (approximately 6 hours) comparative coarse-mesh simulation experiments with and without the synthetic data. They are aimed at revealing the model sensitivities should allow us both to refine the specification of the observational requirements, and to develop more accurate and efficient objective analysis schemes. The goal is to advance the MASS (Mesoscal Atmospheric Simulation System) modeling expertise so that the model output can provide reliable guidance for thunderstorm forecasting.

Chang, Chia-Bo↗