Search NASA⌕ Search

SEARCH · Search NASA

Results for “gradient projection”

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 649 records · Page 36

Computation of graph hitting time moments; Chapel code implementation.

The project that developed this is unclassified, with the mandate to produce open source code. This code computes the hitting time moments of a graph using a linear algebra configuration. The goal of this work is to explore the performance capabilities of the Chapel programming language. Toward that end, we generate random adjacency matrices which represent a random graph. The code can also read in an adjacency matrix from a file. The main computation is the Conjugate Gradient method.SAND2020-12651 M. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Barrett, Richard↗

Program Verification for Extreme-Scale Applications (Final Scientific/Technical Report)

The project seeks to develop tools and techniques to help software developers verify the correctness and accuracy of their code. The target domain is scientific software of the kind widely used and developed in the Department of Energy research community, with a particular focus on "extreme scale" programs - those that are expected to involve possibly millions of parallel threads of execution. The report covers the University of Delaware contribution to the collaborative project. The project had a number of successful outcomes, especially regarding the development and extension of the CIVL software verification framework. CIVL is a verification tool for C or Fortran programs that use MPI, OpenMP, CUDA, and/or Pthreads for parallelization. CIVL went through vast improvements and extensions, and was successfully applied to a number of challenging codes. It found a subtle bug in the Devito PDE framework. It was able to verify the functional correctness of a conjugate gradient solver using a novel probabilistic technique with vanishingly small chance of error. CIVL was used very successfully in verification competitions, and the CIVL solutions to the competition challenges were presented and published. CIVL was also used successfully by other researchers on a computational chemistry kernel.

97 MATHEMATICS AND COMPUTING↗

Carbon content of soil fractions varies with season, rainfall, and soil fertility across a lowland tropical moist forest gradient

Tropical forests contain some of the largest soil carbon (C) stocks on Earth, making them broadly relevant to terrestrial-climate feedbacks, yet our understanding of how their soil organic C (SOC) fractions vary over space and time is limited. We studied effects of season, fertility, and mean annual precipitation (MAP) on the C contents of soil fractions across 14 lowland forests in Panama. We measured free-debris, occluded-debris, and mineral-associated SOC fractions, as well as soluble C associated with each fraction. We hypothesized that mineral-associated SOC would be greatest in infertile, strongly weathered soils with large reactive mineral contents. We also hypothesized that the debris SOC fractions would accumulate during the dry season, reflecting seasonal increases in litterfall. To address this, we compared soil fractions in wet and dry seasons from fertile and infertile soils across a range of 1809 – 2864 mm MAP. The C content (mg C / g soil) of all soil fractions varied with fertility and MAP: specifically, free-debris SOC was greatest in wet, high-fertility soils, and occluded-debris SOC was greater in high-fertility than low-fertility soils. The mineral-associated SOC fraction, which contained the majority of bulk soil C, showed increasing C content with greater MAP in infertile sites, presumably driving similar spatial patterns in the bulk soil. Only the free-debris SOC fraction showed strong seasonal variation, increasing in mass during the dry season. Nitrogen behaved similarly to C. In summary, soil C contents increased with MAP in infertile sites but not fertile sites, driven by the mineral-associated SOC fraction. The dry season had greater free-debris SOC, but this seasonal trend was not apparent in bulk soil C, likely because of the small size of the free-debris SOC fraction. Furthermore, changes in the quantity and seasonality of precipitation, which are projected for tropical forests, might shift spatial and temporal patterns of soil C storage, which would in turn influence forest-climate feedbacks for this C-rich biome.

54 ENVIRONMENTAL SCIENCES↗

Data for Stetten et al. (2025), "Biogeochemical controls on iron speciation and cycling across upland to shoreline gradients in freshwater and estuarine coastal soils (Lake Erie and Chesapeake Bay, United States)"

Coastal environments are dynamic interfaces that mediate carbon and nutrient exchanges between terrestrial landscapes and open waters, but it is unclear how biogeochemical reactions, in particular iron (Fe) redox transformations, affect the understanding and prediction of coastal ecosystem functions. This dataset includes measurements from two freshwater sites in the Western and Central basins of Lake Erie (Ohio, United States) and two estuarine sites in the Chesapeake Bay (Maryland, United States); the analytical results were reported by Stetten et al. (2025) in Science of the Total Environment. It was produced as part of the COMPASS-FME project, which seeks to advance a scalable, predictive understanding of the fundamental biogeochemical processes, ecological structure, and ecosystem dynamics that distinguish coastal terrestrial-aquatic interfaces from the purely terrestrial or aquatic systems to which they are coupled. The sites were sampled in November 2022 (CRC), December 2022 (MSM), February 2023 (GCW), and March 2023 (OWC); site codes follow those used by Pennington et al. (2025).The dataset consists of the following soil data:- Solid data (Fe concentration, etc.)- Porewater data (sulfate, sulfide, etc.)- Linear combination fitting results of X-ray absorption near edge structure (XANES) spectra; i.e., quantitative results of the oxidation state of Fe, indicated as a proportion of pure Fe(III) and Fe(II) model compounds- Linear combination fitting results of EXAFS (extended X-ray absorption fine structure) spectra, indicated as proportion of of Fe-model compounds (illite, smectite, etc.)Each data type has a single file in comma-separated value (CSV) format. No special software is required to read it.

54 ENVIRONMENTAL SCIENCES↗

NASA's B-57B Gust Gradient Program

The B-57B Gust Gradient Program is a joint effort of NASA Headquarters, Marshall Space Flight Center, Dryden Flight Research Facility, Langley Research Center, and Ames Research Center. The primary program goal is to measure spanwise variations of turbulent gusts across an airflow. To this end, the NASA B-57B aircraft was equipped with three component gust probes on each wing tip and on the nose. Early results of flights done in conjunction with the Joint Airport Weather Studies (JAWS) project are described.

Camp, D.↗

Enhanced heavy oil recovery by low salinity polymer flood combined with microgel treatment

“Heavy oil resources account for a large portion of the total oil reserves around the world. The target heavy oil reservoir is located on Alaska’s North Slope (ANS). Advantages of low-salinity HPAM polymer (LSP) over high-salinity polymer (HSP) were demonstrated. LSP could recover more oil with 40% less polymer consumption. No additional oil was recovered by HSP after LSP flood. The first-ever polymer flood pilot on ANS showed remarkable success regarding water cut reduction, oil production increase, delayed breakthrough, and projected oil recovery improvement. Polymer alone was insufficient to achieve satisfactory recovery as the reservoirs were highly heterogeneous. Microgels could improve the effectiveness of polymer flood by reducing water cut and increasing oil recovery. Favorable working conditions were identified. Microgel transport behavior was studied using superpermeable sandpacks (27-221 darcies) with multiple pressure sensors. The particle-to-pore matching size ratio significantly impacted the effectiveness of the gels. A threshold differential pressure (ΔPth) and critical pressure gradient (∇Pcr) were required to push the gels to penetrate and propagate through the channels. The ΔPth and ∇Pcr revealed the underlying mechanisms of selective penetration/placement behavior of microgels in heterogeneous reservoirs. Diagrams were developed to estimate the maximum propagation distance of the gels in channels in conceptual field applications. Sandwich-like channel models and methodologies were developed to comprehensively evaluate the effectiveness of gel materials. Gel retention in the channels was quantified. Results also indicated that the retained gels were dehydrated. Fluid diversion and sweep improvement after gel treatments were evaluated by tracer tests”--Abstract, page iv.

Zhao, Yang↗

Direct observation of oxygen vacancy ordering evolution in cerium oxide at varying concentrations via high-resolution transmission electron microscopy

Under illumination by a 300 keV electron beam, oxygen vacancy ordering structures are induced within cerium oxide grains. Here, our high-resolution transmission electron microscopy (HRTEM) study, supported by image simulation, reveals the evolution of these structures as vacancy concentration increases. The observed fluorite-type superlattice structures are identified as CeO 1.825 , CeO 1.75 , Ce 2 O 3 , displaying a gradient in oxygen vacancy concentration moving away from the grain surface. Correspondingly, the structural sequence transitions from Ce 2 O 3 to CeO 1.75 and then to CeO 1.825 . Without the constraints of surrounding grains, fluorite-type Ce 2 O 3 nanocrystals show a preference for transformation into an A-type trigonal structure. Notably, at temperatures up to 200°C, only the perfect fluorite structure is observed. Structural models were validated through both [110] and [001] projections. Our findings further confirm lattice expansion associated with local oxygen vacancy enrichment, which can be compensated by the formation of stacking faults, where a {111} oxygen plane is lost at defect sites.

Cerium oxide↗

Wood Decomposition: Understanding Processes Regulating Carbon Transfer to Soil Carbon Pools Using FACE Wood at Multiple Scales

Dead wood comprises approximately 20% of forest biomass. While it is widely recognized for its role in ecosystem functions relative to habitat and nutrient cycling, there are considerable uncertainties associated with its role in the forest carbon (C) cycle. The principal uncertainty is whether dead wood contributes to the soil C pool. With little supporting evidence, the presumption is that C in wood is respired back to the atmosphere through decomposition processes. The inability to attribute forest soil C to its source components has been a barrier to addressing this important forest carbon cycle question. The U.S. Dept. of Energy’s Free Air CO 2 Enrichment (FACE) experiment produced wood that has a distinct δ 13 C signature that can be used to trace its fate through organic matter decomposition processes. Using two hardwood species (aspen, birch) from the Rhinelander FACE site and pine from the Duke FACE site, the FACE Wood Decomposition Experiment (FWDE) was established in 2011 to provide the foundation for addressing wood decomposition and its associated processes at the continental scale. Plots were installed on nine experimental forests representing a range in bioclimatic forest conditions within the continental U.S. Each identical plot has logs of the three species incubating on the soil surface and pine logs suspended in the air, representing standing dead wood. Abiotic conditions are monitored continuously, and soil and forest conditions characterized. Logs on the sites have been measured every two years to assess decay. This project was initiated in 2016 with specific objectives to: (a) determine whether dead wood C contributes to the soil C pool, (b) assess the fungal community composition mediating the decay, (c) assess the interactions of termites with microbial decay, and (d) develop a mechanistic model to simulate wood decomposition. The principal field sampling was conducted in 2017; correspondingly, the results reported here represent six years of decomposition. Wood decomposition varied from 13 to 75 % loss across the bioclimatic gradient. Sites in the southeastern U.S. coastal plain and piedmont had fastest decomposition rates for each of the species. In contrast sites in the Rocky Mountains had the slowest decay rates, which were associated with cold temperatures. Standing deadwood decomposes at a much slower rate than logs on the soil surface; a response due to differences in wood moisture content. Wood carbon was measured in the mineral soil across each of the sites; accordingly, this is the first study to affirm that a portion of the wood C does become part of the soil C pool. Averaged across the FWDE sites, approximately 1.4% of the soil C in the upper 10 cm of the mineral soil directly beneath the log was derived from the deadwood; within 50 cm of the log, wood C accounted for an average of 2.1%. These results demonstrate that the contribution of deadwood to the soil C pool is a soil volume surrounding the log; and it’s likely a larger volume than what we sampled. We also hypothesized that there would be more wood C in the mineral soil when termites were present. However, that hypothesis wasn’t supported by measurements on the Santee Experimental Forest termite exclusion (SEF-TX) plots under FACE logs that incubated for 6 years. Across all the sites, the majority of the wood C was transferred to the atmosphere. Approximately 7,000 fungal species were sequenced. The fungal communities were different among each of the sites suggesting local colonization; interestingly, within sites the three species also exhibited differences in fungal community composition. White rot was the dominant fungal form in pine and hardwood, which yields less lignin in the wood residue. Brown rot occurred in approximately 20% of the pine and aspen logs, and it did not occur in birch logs. The occurrence of brown rot was only weakly correlated to site factors. We conducted two additional experiments to assess the effects of termites on wood decomposition and the associated interactions with the fungal community. Pine and birch logs were incubated on the Santee Experimental Forest termite exclusion plots in upland and riparian settings. There was no difference in decay between areas where termites were able colonize and plots where they had been excluded. Highly variable foraging by termites was a primary contributor the result. Correspondingly, fungal communities did not vary on logs incubated among the two treatments, but they did reflect site and species differences. The second experiment was conducted on each of the nine FWDE sites where pine and birch wood blocks were incubated for two years with and without a wrapping designed to exclude termites and other arthropods. The effects of termites on wood decomposition could be confirmed on two of the four sites known to have subterranean termites. On the two southeastern sites, termites contributed approximately 8% to mass loss in pine blocks over two years. A major factor clouding the role of termites in wood decomposition is the inability to assess the time, duration, and intensity of the foraging with the logs. Another uncertainty is related to whether a piece of wood will be foraged. A process-based wood decomposition model was developed representing biological processes mediated by fungi, bacteria, termites, beetles, as well as physical and chemical processes. The Coarse Wood Decomposition Model (CWDM) was parameterized based on a thorough review of the literature and it was validated using data from the FWDE. In the initial setup, the model considers softwoods and hardwoods groups, as opposed to individual species. Sensitivity analyses incorporating climate data from 89 locations in North America demonstrated that CWDM was stable from the tropics to the boreal zone. Model calibration and validation showed it to perform very well. The CWDM is poised for additional testing with subsequent FWDE data as well as other data sets, and it is available for incorporation into forest biogeochemical models or DOEs Earth System Modeling framework.

54 ENVIRONMENTAL SCIENCES↗

Development of Transformative Preparation Methods to Push up High Q&G Performance of FRIB Spare HWR Cryomodule Cavities

The FRIB accelerator project construction, a top priority of US nuclear science, was completed in January 2022, and is now moving to user operation. The stable and reliable operation of the accelerating cryomodules is essential in achieving/fulfilling DOE and user expectations. So far, FRIB cryomodules meet all FRIB specifications for cavity performance. However, during the lifetime of machine operation, degradation of cryomodule performance is possible, as reported in similar operating facilities (CEBAF, SNS). If cryomodule degradation is observed at FRIB, the under-performing cryomodule will require replacement/maintenance. In effort to manage operational reliability, FRIB plans to construct a 0.53 half-wave cryomodule to serve as an active spare. In a parallel effort, FRIB will also work toward increasing operational Q and gradient of spare cryomodule cavities to gain an overall performance margin to support future operational reliability. The current FRIB cavity designs have a potential to operate at gradients higher than 8 MV/m, but are currently limited by field emission (FE) and/or high field Q slope (HFQS); known issue in buffered chemical polished (BCP) treated cavities. The proposal looks to develop transformative surface preparation treatments to improve the operational gradient of spare cryomodules higher than 10 MV/m while maintaining high Q. Thus, increasing operational margin by 30 - 50%. With the goal to improve operational reliability set, the proposal will investigate multiple objectives as possible paths forward to achieve an overall increase in cavity performance and gain a better understanding of SRF limiting mechanisms. The proposal will study the application of different chemical surface treatments to 0.53 half-wave cavities, with the addition of low temperature bakes (LTB), and measure their effects on accelerating performance. Proposed chemical treatments to be explored in this proposal include conventional EP acid mixtures, as well as innovated EP and BCP acid mixtures designed to simplify processing paths in migrating FE and HFQS. The proposed transformative treatment wet N-doping also has the potential to replicate recent advancements in SRF technology relating to nitrogen doping and high Q operation without the requirement for an ultra-high vacuum annealing furnace; currently being developed at FNAL and JLAB. In parallel, high Q performance relating to flux trapping will be investigated with the installation of a second layer of magnetic shielding in the vertical test Dewar. The research objectives presented in the proposal, and their corresponding effects on cavity performance, will provide essential knowledge and future guidance to the SRF community and provide possible paths for future SRF based projects and applications.

43 PARTICLE ACCELERATORS↗

Topography-Induced TKE Budget Behavior Over an Amazon Forest

Data from three different heights (35, 50 and 81 m) of one of the ATTO Project towers in the Amazon forest were used to calculate the TKE (turbulence kinetic energy) budget and some other statistics within the RSL (roughness sublayer). The statistical analyses were carried out for unstable and stable cases. The vertical transport and the vertical advection terms do not explain the imbalances found in the TKE budget, highlighting the impor- tance of horizontal terms in complex terrain. Indeed, patterns independent of the time of day were found for vertical transport, horizontal and vertical velocity skewness, the mean vertical component velocity, and momentum and vertical TKE flux. Here, the highest values of mean vertical velocity were observed from the direction where there is a valley near the tower, which may indicate a topographical effect. Monin-Obukhov Similarity Theory (MOST) predictions for the mean velocity gradient, strictly applicable only in the inertial sublayer (ISL), gave reasonable predictions at 35 and 50 m, but dimensionless dissipation rates displayed large scatter. At 81 m, there were clear deviations from MOST for the two functions, disclosing the effect of topography. In addition, a neutral LES (large-eddy simu- lation) study showed good agreement with tower data for the same general wind direction.

ATTO project↗

Thermal analysis of heat storage canisters for a solar dynamic, space power system

A thermal analysis was performed of a thermal energy storage canister of a type suggested for use in a solar receiver for an orbiting Brayton cycle power system. Energy storage for the eclipse portion of the cycle is provided by the latent heat of a eutectic mixture of LiF and CaF2 contained in the canister. The chief motivation for the study is the prediction of vapor void effects on temperature profiles and the identification of possible differences between ground test data and projected behavior in microgravity. The first phase of this study is based on a two-dimensional, cylindrical coordinates model using an interim procedure for describing void behavor in 1-g and microgravity. The thermal analysis includes the effects of solidification front behavior, conduction in liquid/solid salt and canister materials, void growth and shrinkage, radiant heat transfer across the void, and convection in the melt due to Marangoni-induced flow and, in 1-g, flow due to density gradients. A number of significant differences between 1-g and o-g behavior were found. This resulted from differences in void location relative to the maximum heat flux and a significantly smaller effective conductance in 0-g due to the absence of gravity-induced convection.

Wichner, R. P.↗

Study of Self-Interaction Errors in Density Functional Calculations of Magnetic Exchange Coupling Constants Using Three Self-Interaction Correction Methods

We examine the role of self-interaction errors (SIE) removal on the evaluation of magnetic exchange coupling constants. In particular we analyze the effect of scaling down the self-interaction-correction (SIC) for three non-empirical density functional approximations (DFAs) namely, the local spin density approximation, the Perdew-Burke-Ernzerhof generalized gradient approximation, and recent SCAN family of meta-GGA functionals. To this end, we employ three one-electron SIC methods: Perdew-Zunger [Perdew, J. P.; Zunger, A. Phys. Rev. B, 1981, 23, 5048] SIC, the orbitalwise scaled SIC method [Vydrov, O. A. et al., J. Chem. Phys. 2006, 124, 094108], and the recent local scaling method [Zope, R. R. et al., J. Chem. Phys. 2019, 151, 214108]. We compute the magnetic exchange coupling constants using the spin projection and non projection approaches for sets of molecules composed of dinuclear and polynuclear H-He models, organic radical molecules, and chlorocuprate, and compare these results against accurate theories and experiment. Our results show that for the systems that mainly consist of single electron regions, PZSIC performs well but for more complex organic systems and the chlorcuprates, an overcorrecting tendency of PZSIC combined with the DFAs utilized in this work is more pronounced, and in such cases LSIC with kinetic energy density ratio performs better than PZSIC. Furthermore, analysis of the results in terms of SIC corrections to the density and to the total energy shows that both density and energy correction are required to obtain an improved prediction of magnetic exchange couplings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Sedimentary Geothermal Resources in Nevada, Utah, Colorado, and Texas

The objectives of this project were to (1) perform a literature review of sedimentary geothermal resources, (2) identify data sources and develop data-collection methodologies that characterize selected resources, (3) screen sedimentary basins and formations for sedimentary geothermal potential, and (4) evaluate the technical feasibility of one or more selected locations. Numerous publications have characterized geothermal resources within sedimentary basins. A literature search reviewed publications describing resources located in Colorado, Louisiana, Nevada, Texas, Utah, and Wyoming. The most attractive resources have high temperature gradients, low drilling costs, and reservoir permeabilities greater than 10 millidarcies (mD). Prospects in Colorado, Nevada, Texas, and Utah exhibit attractive characteristics and were chosen for further analysis. Sedimentary resources in Nevada and Utah are most attractive, followed by tested resources in Texas and untested resources in Colorado. The identified resources in Wyoming and Louisiana had lower geothermal gradients and were not evaluated. Reservoir modeling and techno-economic analysis were performed at Marys River Basin–North in Nevada. Geothermal energy production at this location is expected to have a levelized cost of energy (LCOE) ranging between 10 and 20 cents/kWh. Additional work may result in lower LCOE estimates at this location and at other attractive prospects in these three regions. The Great Basin carbonate and alluvial aquifer system of eastern Nevada and western Utah includes a lower carbonate aquifer unit, which has the potential for hosting both conduction- and convection-dominated geothermal systems. Mapping has identified lateral thickness variability expressed as a roughly 120–150-km-wide central corridor, which hosts the thickest and most continuous formations and extends from near Las Vegas north to the Idaho border. Purchase and analysis of privately held legacy seismic data could potentially compensate for the lack of sufficient data documenting measured depth to the lower carbonate aquifer unit. Multiple orogenies, extension episodes, and intrusive events have deformed and displaced the target formations of the lower carbonate aquifer unit. This structural complexity and potential for dipping reservoirs emphasizes the need for detailed geological, geophysical, and reservoir modeling. Heat flow within three Colorado sedimentary basins reviewed as part of this study was calculated in targeted studies by the Colorado Geologic Survey and Colorado School of Mines. These calculations are based on bottom-hole temperature data sets with significant limitations and some variability but produce values consistently higher than the global continental average of 65 mW/m2 for all three basins. Heat flow in the Raton Basin is the highest; however, permeability measurements from specific sedimentary formations with high heat flow have not been obtained. Promising formations for sedimentary geothermal systems were found in all three regions studied—Nevada-Utah, Colorado, and Texas. The next steps for developing sedimentary geothermal resources vary due to differences in available data and resource uncertainties. Additional scopes of work are recommended for identified basins in these three regions.

15 GEOTHERMAL ENERGY↗

Dual laser optical system and method for studying fluid flow

A dual laser optical system and method is disclosed for visualization of phenomena in transport substances which induce refractive index gradients such as fluid flow and pressure and temperature gradients in fluids and gases. Two images representing mutually perpendicular components of refractive index gradients may be viewed simultaneously on screen. Two lasers having wave lengths in the visible range but separated by about 1000 angstroms are utilized to provide beams which are collimated into a beam containing components of the different wave lengths. The collimated beam is passed through a test volume of the transparent substance. The collimated beam is then separated into components of the different wave lengths and focused onto a pair of knife edges arranged mutually perpendicular to produce and project images onto the screen.

Owen, R. B.↗

Large Coded Aperture Mask for Spaceflight Hard X-ray Images

The 2.6 square meter coded aperture mask is a vital part of the Burst Alert Telescope on the Swift mission. A random, but known pattern of more than 50,000 lead tiles, each 5 mm square, was bonded to a large honeycomb panel which projects a shadow on the detector array during a gamma ray burst. A two-year development process was necessary to explore ideas, apply techniques, and finalize procedures to meet the strict requirements for the coded aperture mask. Challenges included finding a honeycomb substrate with minimal gamma ray attenuation, selecting an adhesive with adequate bond strength to hold the tiles in place but soft enough to allow the tiles to expand and contract without distorting the panel under large temperature gradients, and eliminating excess adhesive from all untiled areas. The largest challenge was to find an efficient way to bond the > 50,000 lead tiles to the panel with positional tolerances measured in microns. In order to generate the desired bondline, adhesive was applied and allowed to cure to each tile. The pre-cured tiles were located in a tool to maintain positional accuracy, wet adhesive was applied to the panel, and it was lowered to the tile surface with synchronized actuators. Using this procedure, the entire tile pattern was transferred to the large honeycomb panel in a single bond. The pressure for the bond was achieved by enclosing the entire system in a vacuum bag. Thermal vacuum and acoustic tests validated this approach. This paper discusses the methods, materials, and techniques used to fabricate this very large and unique coded aperture mask for the Swift mission.

Vigneau, Danielle N.↗

Coupled Gravity and Elevation Measurement of Ice Sheet Mass Change

During June 2003, we measured surface gravity at six locations about a glaciological measurement site located on the South-central Greenland Ice. We operated a GPS unit for 90 minutes at each site -the unit was operated simultaneously with a base station unit in Sondrestrom Fjord so as to enable differential, post-processing of the data. We installed an aluminum, accumulation-rate-pole at each site. The base section of the pole also served as the mount for the GPS antenna. Two gravimeters were used simultaneously at each site. Measurements were repeated at each site with at time lapse of at least 50 minutes. We measured snow physical properties in two shallow pits The same measurement sites were occupied in 1981 and all were part of a hexagonal network of geodetic and glaciological measurements established by The Ohio State University in 1980. Additional gravity observations were acquired at three of the sites in 1993 and 1995. Gravity data were collected in conjunction with Doppler satellite measurements of position and elevation in 1981 and global positioning system measurements subsequently. The use of satellite navigation techniques permitted reoccupation of the same sites in each year to within a few 10 s of meters or better. After detrending the gravity data, making adjustments for tides and removing the residual effects of local spatial gradients in gravity, we observe an average secular decrease in gravity of about 0.01 milligal/year, but with tenths of milligal variations about the mean trend. The trend is consistent with a nearly linear increase in surface elevation of between 7 to 10 c d y r (depending on location) as measured by repeated airborne laser altimeter, surface Doppler satellite and GPS elevation measurements. Differences between the residual gravity anomalies after free air correction may be attributable to local mass changes. This project is a collaboration between the Byrd Polar Research Center of the Ohio State University and the Arctic Technology Center of the Danish Technical University.

Jezek, K. C.↗

Failure Analysis for Molten Salt Thermal Energy Storage Tanks for In-Service CSP Plants

Thermal Energy Storage (TES) is a fundamental component in concentrating solar power (CSP) plants to increase the plant's dispatchability, capacity factor, while reducing the levelized cost of electricity. In central receivers CSP plants, nitrate molten salts have been used for several years for operation temperatures of up to 565 degrees C. Despite many efforts to advance nitrate salt to higher operation temperatures (even considering a replacement with molten chloride salts) to achieve higher energy conversion efficiencies, the 565 degrees C temperature is currently considered the state-of-the art. Although molten salt tanks have been broadly deployed in commercial CSP plants worldwide, several failures have been reported in these tanks after a few months or years of operation, causing significant economic loss and mistrust in CSP technologies. Most of these failures are associated with the infancy of the technology and multiple issues related to tank design, fabrication, commissioning, and aggressive operation. A technical standard dedicated to the design and fabrication of molten nitrate TES tanks does not exist today. Current in-service molten salt tanks have been generally designed based on the American Petroleum Institute's (API) 650 and ASME Section II standards. The API 650 code provides guidelines for dimensions and fabrication for oil storage tanks up to 260 degrees C. The ASME standard provides allowable stress values for various materials at a range of temperatures and conditions. Both standards seem to be limited for molten salt TES tanks where high temperatures, thermal cycling, and transient conditions are expected. In 2020, NREL released the Concentrating Solar Power Best Practices Study (NREL/TP-5500-75763) that summarized multiple issues in CSP plants, along with potential alternatives and recommendations to address those issues based on information collected from participants representing about 80% of operating CSP plants in the world. One of the recommendations from this study was the development of accurate and validated models to evaluate the plant's transient operation, capable of capturing the effect of short-term clouds and operator response, while being flexible in being adapted to various spatial and temporal resource data. The "Failure Analysis for Molten Salt Thermal Energy Tanks for In-Service CSP Plants" project was inspired on this recommendation and was focused on (1) the development and validation of a physics-based model for a representative, commercial-scale molten salt tank, (2) performing simulations to evaluate the behavior of the tank as a function of typical plant operation conditions, (3) understanding tank failures mechanisms, (4) determining the residual stress and distortion in the tank floor after welding fabrication and evaluating their impact in the stresses developed in the tank during operation, (5) assessing the impact of key operation parameters on the temperature and stress distribution, (6) conduct a preliminary evaluation of design features to reduce stress and improve tank's reliability, and (7) estimate the tank's service life based on the stress developed under diverse operation scenarios. From the analysis conducted in the project and presented in this report, it was found that maximum stresses surpassing the yield strength point of the stainless steel (SS) 347H are developed on the tank floor near the perimeter. These large stresses are strongly influenced by the initial residual stresses and distortion of the tank floor after welding fabrication. During operation, large stresses are developed in the tank floor at high operation temperatures with large salt inventory levels during transient operation. High stresses are also related to elevated temperature gradients in the tank floor that could be attributed to insufficient mixing within the salt inflow and the salt inventory. Based on the analysis, creep is the predominant failure mechanism. However, the large stress levels could favor the plastic deformation into buckles, and crack formation due to stress relaxation cracking during cycle operation. A lifetime below 3 years was estimated for the typical plant operation conditions and a specific initial residual stress and deformation distribution of the tank floor. The estimated life agrees with the service time to failure reported in several commercial molten salt tanks. Desing and operation guidelines can be extracted from the analysis presented in this report, which could be adopted by tank manufacturers and CSP operators to advance toward an ultimate solution for tank failures by reducing residual and operational stresses to achieve a tank service life of more than 30 years. Addressing failures in molten salt TES tanks is fundamental for the CSP industry's survivability, but it is also important for other industrial and power generation applications using this technology, including nuclear and concentrating solar thermal.

14 SOLAR ENERGY↗

Thermal Characterization and Exhumation of Northwest San Juan Basin Area, NM

The San Juan Basin (SJB), located in southwestern Colorado northwestern New Mexico, containing Paleozoic through middle Cenozoic strata, formed as a partitioned basin approximately 80 Ma in response to the Laramide orogeny. The SJB is a commercially mature, petroleum-producing basin, and currently being explored for CCUS and geothermal resources. Modern heat flow within the SJB is spatially variable with higher geothermal gradients in its northern and eastern portions. The temporal history of the basin thermal history is essential for understanding capacity for carbon storage and geothermal exploration. We present thermochronometric analyses of two surface and two subsurface samples to constrain spatial-temporal thermal evolution of the northwestern San Juan Basin. Apatite (U-Th)/He thermochronometric analyses were conducted on Cliff House and Kirtland Formation outcrops and subsurface samples of the Ojo Alamo and Pictured Cliffs Formations from a DOE funded, CCUS project pilot-well. In addition, vitrinite reflectance data and 1D basin modelling constrain possible time-temperature pathways of the samples in an area where no previous studies have constrained the uplift/exhumation of the Hogback Monocline using thermochronometry. Modeling suggests that elevated subsurface temperatures developed simultaneously with regional late Oligiocene volcanism in the San Juan Volcanic field and timing of maximum burial in the SJB. Late Miocene to Pliocene cooling/exhumation through the apatite (U-Th)/He partial retention zone support stratigraphic evidence of the Colorado Plateau uplift, possibly regional epeirogenic uplift driven by mantle processes.

02 PETROLEUM↗