Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Gradient information”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

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↗
DOE OSTI2024

Can Error Mitigation Improve Trainability of Noisy Variational Quantum Algorithms?

Variational Quantum Algorithms (VQAs) are often viewed as the best hope for near-term quantum advantage. However, recent studies have shown that noise can severely limit the trainability of VQAs, e.g., by exponentially flattening the cost landscape and suppressing the magnitudes of cost gradients. Error Mitigation (EM) shows promise in reducing the impact of noise on near-term devices. Thus, it is natural to ask whether EM can improve the trainability of VQAs. In this work, we first show that, for a broad class of EM strategies, exponential cost concentration cannot be resolved without committing exponential resources elsewhere. This class of strategies includes as special cases Zero Noise Extrapolation, Virtual Distillation, Probabilistic Error Cancellation, and Clifford Data Regression. Second, we perform analytical and numerical analysis of these EM protocols, and we find that some of them (e.g., Virtual Distillation) can make it harder to resolve cost function values compared to running no EM at all. As a positive result, we do find numerical evidence that Clifford Data Regression (CDR) can aid the training process in certain settings where cost concentration is not too severe. Our results show that care should be taken in applying EM protocols as they can either worsen or not improve trainability. On the other hand, our positive results for CDR highlight the possibility of engineering error mitigation methods to improve trainability.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2000

Gravitational Effects on Flow Instability and Transition in Low Density Jets

Recent experiments have shown that low-density gas jets injected into a high-density ambient gas undergo an instability mode, leading to highly-periodic oscillations in the flow-field for certain conditions. The transition from laminar to turbulent flow in these jets is abrupt, without the gradual change in scales. Even the fine scale turbulent structure repeats itself with extreme regularity from cycle to cycle. Similar observations were obtained in buoyancy-dominated and momentum-dominated jets characterized by the Richardson numbers, Ri = [gD(rho(sub a)-rho(sub j))/rho(sub j)U(sub j)(exp 2) ] where g is the gravitational acceleration, D is the jet diameter, rho(sub a) and rho(sub a) are, respectively, the free-stream and jet densities, and U(sub j) is the mean jet exit velocity. At high Richardson numbers, the instability is presumably caused by buoyancy since the flow-oscillation frequency (f) or the Strouhal number, St = [fD/U(sub j)] scales with Ri. In momentum-dominated jets, however, the Strouhal number of the oscillating flow is relatively independent of the Ri. In this case, a local absolute instability is predicted in the potential core of low-density jets with S [= rho(sub j)/rho(sub a)] < 0.7, which agrees qualitatively with experiments. Although the instability in gas jets of high Richardson numbers is attributed to buoyancy, direct physical evidence has not been acquired in experiments. If the instability is indeed caused by buoyancy, the near-field flow structure of the jet will change significantly when the buoyancy is removed, for example, in the microgravity environment. Thus, quantitative data on the spatial and temporal evolutions of the instability, length and time scale of the oscillating mode and its effects on the mean flow and breakdown of the potential core are needed in normal and microgravity to delineate gravitational effects in buoyant jets. In momentum dominated low-density jets, the instability is speculated to originate in the potential core. However, experiments have not succeeded in identifying the direct physical cause of the instability. For example, the theory predicts an oscillating mode for S<0.62 in the limit of zero momentum thickness, which contradicts with the experimental findings of Kyle and Sreenivasan. The analyses of momentum-dominated jets neglect buoyancy effects because of the small Richardson number. Although this assumption is appropriate in the potential core, the gravitational effects are important in the annular region surrounding the jet, where the density and velocity gradients are large. This reasoning provides basis for the hypothesis that the instability in low Richardosn number jets studied by Kyle and Sreenivasan and Monkewitz et al. is caused by buoyancy. The striking similarity in characteristics of the instability and virtually the identical conclusions reached by Subbarao and Cantwell in buoyant (Ri>0.5) helium jets on one hand and by Kyle and Sreenivasan in momentum-dominated (Ri<1x10(exp -3)) helium jets on the other support this hypothesis. However, quantitative experiments in normal and microgravity are necessary to obtain direct physical evidence of buoyancy effects on the flow instability and structure of momentum-dominated low-density jets. The primary objective of this new research project is to quantify how buoyancy affects the flow instability and structure in the near field of low-density jets. The flow will be described by the spatial and temporal evolutions of the instability, length and time scales of the oscillating mode, and the mean and fluctuating concentration fields. To meet this objective, concentration measurements will be obtained across the whole field using quantitative Rainbow Schlieren Deflectometry, providing spatial resolution of 0.1mm and temporal resolution of 0.017s to 1ms. The experimental effort will be supplemented with linear stability analysis of low-density jets by considering buoyancy. The first objective of this research is to investigate the effects of gravity on the flow instability and structure of low-density jets. The flow instability in these jets has been attributed to buoyancy. By removing buoyancy in our experiments, we seek to obtain the direct physical evidence of the instability mechanism. In the absence of the instability, the flow structure will undergo a significant change. We seek to quantify these changes by mapping the flow field (in terms of the concentration profiles) of these jets at non-buoyant conditions. Such information is presently lacking in the existing literature. The second objective of this research is to determine if the instability in momentum-driven, low-density jets is caused by buoyancy. At these conditions, the buoyancy effects are commonly ignored because of the small Richardson based on global parameters. By eliminating buoyancy in our experiments, globally as well as locally, we seek to examine the possibility that the instability mechanism in self-excited, buoyant or momentum-driven jets is the same. To meet this objective, we would quantify the jet flow in normal and microgravity, while systematically decreasing the Richardson number from buoyancy-driven to momentum driven flow regime. The third objective of this research is to perform a linear stability analysis of low-density gas jets by including the gravitational effects. The flow oscillations in these jets are attributed to an absolute instability, whereby the disturbance grows exponentially at the site to ultimately contaminate the entire flow field. We seek to study the characteristics of both convective and absolute instabilities and demarcate the boundary between them.

Agrawal A. K.↗
DOE OSTI2023

Ectomycorrhizal effects on decomposition are highly dependent on fungal traits, climate, and litter properties: A model-based assessment. Dataset.

To simulate the effects of mycorrhizal fungi on soil organic matter cycling, we incorporated mycorrhizal processes into the Carbon, Organisms, Rhizosphere, and Protection in the Soil Environment (CORPSE) model to develop a new soil model Myco-CORPSE. The new model was calibrated and evaluated against soil measurements taken at temperate forests in New Hampshire (NH) and Georgia (GA). A series of scenario analysis were also conducted to explore the conditions under which ectomycorrhizal (ECM) N acquisition processes can induce different soil C accumulation in ECM systems compared to arbuscular (AM) systems.In this data package, we included:-The Python codes of the standard Myco-CORPSE model we developed: "Standard Myco_CORPSE python codes.zip". The main program is the "gradient_sim.py" which calculates the bulk soil microbes and CN content along a user defined gradient of clay, soil temperature, soil moisture and mycorrhizal dominance, and relies on two subprograms "CORPSE_deriv.py" and "CORPSE_integrate.py". "CORPSE_deriv.py" calculated the changes in all simulated soil stock within every time step and "CORPSE_integrate.py" integrate the changes in all simulated soil stock within simulated time period. The program "Plot.py" is used to plot the major outputs produced by the main program "gradient_sim.py".-The modified Python codes of Myco-CORPSE models with site-level environmental inputs (in NH and GA) used to conduct simulations in NH and GA sites: "NH_GA model simulations.zip". -The Python codes used to evaluate the Myco-CORPSE simulation outputs in NH and GA sites against site-level measurements: "Plot NH_GA simulation against measurements.zip". It includes both the evaluation Python code, the model outputs on NH and GA sites, and the measured soil properties in both sites.-The modified Python codes of Myco-CORPSE models "Scenario analysis_model simulations.zip" that is used to conduct scenario analysis of how different litter properties, mycorrhizal fungal traits, climate, and seasonal variation in temperature and vegetation phenology impact the mycorrhizal effects on soil CN properties. The sub file folder "Scenario analysis_litter traits" contains the codes for scenario analysis of different litter properties; The sub file folder "Scenario analysis_ECM types" contains the codes for scenario analysis of different ECM fungal traits; The sub file folder "Scenario analysis_climate&seasonality" contains the codes for scenario analysis of different climate and seasonalities;-"Scenario analysis_model results and plotting codes.zip" contains all the output files from the the scenario analysis of Myco-CORPSE model as described above and the plotting codes used to the generate the heatmaps and scatterplots shown in the manuscript "Ectomycorrhizal effects on decomposition are highly dependent on fungal traits, climate, and litter properties: A model-based assessment"The majority of the model outputs did not have specific geographic information or temporal coverage because the analysis we conducted are mainly hypothetical model simulations. We only provided geographic description, coordinates and temporal coverage for those soil measurements which we used for model evaluations (included in the "Plot NH_GA simulation against measurements.zip").

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Hawai‘i Play Fairway (Final Report)

Most of Hawai'i's geothermal resources are blind—their manifestations, such as hot springs and steam vents, do not appear on the ground surface because the heated water flows far below. With the exception of $K\bar{i}lauea East Rift Zone$, in most areas of Hawai'i, high lateral permeability in the first kilometer below ground surface prevents surface thermal features from developing. As a methodology for discovering these blind resources, Play Fairway Analysis (PFA) involves finding potential locations of blind hydrothermal systems and describing potential geothermal sources in rift-zone settings. Using the PFA to find Hawai'i's geothermal resources, the University of Hawai'i (UH) conducted the Hawai‘i Play Fairway Project, Hawai'i's first statewide geothermal resource assessment since 1985. Sponsored by the U.S. Department of Energy, the Hawai'i Play Fairway Project provided an updated resource assessment, a roadmap for additional exploration activities, and the identification of areas for further exploration. Benefitting from UH's core competency in earth sciences and experienced geothermal researchers, the project comprised three phases. During the first phase, the team identified, compiled, and ranked existing geologic, groundwater, and geophysical datasets relevant to subsurface heat, fluid and permeability. Using a Bayesian statistical approach, the team developed a statistical methodology to integrate these data into a resource probability map. The team evaluated the confidence in the probability value and considered development viability of areas with geothermal resources. With these analyses, the team identified 10 locations in the Hawaiian Islands for exploration activities. For the second phase, the team collected new groundwater data in 10 locations across the state and new geophysical data on $L\bar{a}na‘i, Maui$, and central Hawai'i Island and modeled topographically induced stress to better characterize subsurface permeability. Analyzing the subsurface stresses, the team evaluated the potential for fracture-induced permeability. The team inverted the MT and gravity data to produce 3D models of resistivity and density, respectively, on $L\bar{a}na‘i$, across $Haleakal\bar{a}'s$ SW rift (Maui), and surrounding Mauna Kea (Hawai‘i Island). The team developed and applied a new method for incorporating depth information about resistivity, density, and potential for fracture-induced permeability into the statistical method for computing resource probability in these three focus areas. The team incorporated the new groundwater results with the new geophysical results and the calculations of potential for fracture-induced permeability to produce updated maps of resource probability and confidence. Through combining data from the first and second phases, the team determined locations for further exploration during the third phase. For MT and gravity surveys, the team recommended $Kaua'i's$ $L\bar{i}hu'e$ $Basin$, the east rift of $Maui's$ $Haleakal\bar{a}$ volcano, and the southwest rift of Hawai'i Island's Mauna Loa volcano. The MT and gravity surveys aimed to enable improved confidence in the resource potential in these locations. For drilling deep groundwater well(s), the team recommended Southeast Mauna Kea and $L\bar{a}na's$ $P\bar{a}l\bar{a}wai$ $Basin$. During the third phase, further exploration involved drilling a groundwater well in $L\bar{a}na's$ $P\bar{a}l\bar{a}wai$ $Basin$ and performing more geophysical surveys. We deepened an existing water well proximal to our target area on $L\bar{a}na'i$ due to funding constraints that precluded us from spudding a new well that would exceed 1km depth. Drilling was preceded by a number of substantial elements including: writing an Environmental Assessment and the subsequent legal process, performance of deviation logging, lowering a camera down the well, coordinating site preparation with $P\bar{u}lama$ $L\bar{a}na'i$, shipping the UH-owned rig interisland, procuring supplies, and leading 3 community meetings on $L\bar{a}na'i$. Drilling occurred 24/7 the entire month of June 2019 over which time $L\bar{a}na'i$ $Well$ 10 was deepened from 427 m to 1057 m, with continuous core collected. We measured a roughly linear temperature gradient averaging 42°C/km and a maximum bottom hole temperature of 66°C. This gradient is more than twice the background for Hawai'i and within a range of gradients measured in this depth range for some exploration wells within KERZ. We consider these results encouraging for $L\bar{a}na'i's$ resource potential and recommend following with a slim hole within $L\bar{a}na'i's$ caldera (our target zone) to ~ 2 km. Further, the positive implications such results have for the island of O‘ahu are substantial - the shield stage of O'ahu's volcanoes ended 1-2 My earlier. However, O'ahu uses more electricity than the rest of the islands combined, and the utility recently called for 500-700MW of firm, dispatchable renewable electricity on O'ahu by 2033. In Phase 3, we also collected limited new encouraging groundwater data, and updated our thoughts on the probabilities of fluid and permeability at resource depths (Pr F = 1; Pr P = mostly unconstrained). Ultimately, we advocate for using our final probability of heat, and confidence in this probability, to drive the next phase of exploration. We contend further development of geothermal in Hawai‘i will enable the state to achieve its 100% renewable policy objective and Hawai'i to transition off of fossil fuels through geothermal discovery and development. The project not only produced a large amount of data and expanded the existing knowledge of Hawai'i's geothermal resources, but also produced publications, theses, presentations, core photos, datasets, media reports, television interviews, community events, and a blog. Students and new professionals benefitted from the project's hands-on research experiences and educational opportunities and earned awards and recognition.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

LLNL Kimberlina 1.2 NUFT Simulations June 2018 (v2)

This dataset contains the output 6,000, 3-dimensional reactive multi-phase flow and transport aquifer simulations of brine and CO2 leakage into a protective aquiver in California’s San Joaquin Valley and input data files detailing the geologic mesh, aquifer physical properties and CO2 and brine injection rates. This data set was generated as an ongoing effort with the US DOE National Risk Assessment Partnership (NRAP) to evaluate the effectiveness of monitoring techniques to detect brine and CO2 leakage from legacy wells into underground sources of drinking water overlaying a CO2 storage reservoir. Each simulation contains a unique set of input parameters, generated stochastically. The outputs consist of these upper three geologic layers (from top): the Etchegoin, Macoma-Chanac, Santa Margarita-McLure formations. These simulations span the several distances (1, 3 and 6 km or wells W31-0.2, W31-0.5 and W31-1.0, respectively) from the CO2 injector, initiated from bottom hole pressure and saturation to calculate wellbore leakage from the storage reservoir, with low and high regional groundwater gradients and wellbore leakage into 5 leaky nodes. The dataset includes 1,000 unique simulations for each distance, which each contain a unique aquifer heterogeneity, aquifer and caprock permeability, and two model generations are included with a high permeability (prod07) and hybrid permeability (prod09). The range of permeability distributions is listed in Table 1. Each model generation consists of 3,000 simulations. Included in the dataset are the leakage rates determined from 2D wellbore models which utilize the pressure and CO2 saturation from LBL's reservoir simulations, NUFT mesh files with distributed lithology, NUFT rocktab files which describe the material properties for the geologic layers and the NUFT input files and post-processed output 'ntab' files. Each ntab file contains spatial (rows) and temporal (columns) model output tables for each model cell, the locations (x,y,z) and dimensions for each cells (dx, dy, dz). Table 1. Permeability distribution ranges for prod07 and prod09 model generations Geologic Layer: Permeability Range (log10 m^2) prod07 prod09 Etchegoin -12.92 to -10.92 -13.70 to -11.44 Macoma-Chanac -12.72 to -10.72 -13.50 to -11.24 Santa Margarita-McLure -12.70 to -10.70 -13.48 to -11.22 The input files used to generate the model include which are included in the dataset are: Time series of CO2 leakage input into the model (ex: Q_brn.W31-0.2.sim1000.layers123.tab) Time series of CO2 leakage input into the model (ex: Q_CO2.W31-0.2.sim1000.layers123.tab) Physical properties of the aquifer materials detailing the aquifer porosity, solid density, partitioning coefficients, permeabilities and van-Genuchten parameters detailed in a NUFT rocktab file: (ex: sim1000.usnt.rocktab) Numerical mesh and geologic data assigned to each model cell detailed in a NUFT genmsh format (ex: sim1000.mesh_k16.prod07.trans.genmsh) The primary output parameters are: pH (use absolute value) Change in TDS (mg/kg) Change in Pressure (Pa) Change CO2 gas saturation (fraction range 0.0-1.0) for example, the directory /p/lscratchh/mansoor1/nrap/kimberlina/prod09/mainfiles/sim1000/W31- 0.2 contains: sim1000.W31-0.2.trans.pH.red.ntab sim1000.W31-0.2.no_bg.trans.TDS.red.ntab sim1000.W31-0.2.usnt.P.deltabg.red.ntab sim1000.W31-0.2.usnt.CO2_sat.deltabg.red.ntab Each row in the NTAB files consist of model output per numerical grid cell. Each output file contains 33 columns (variables), including the information of numerical records, geologic location and sizes and the simulated parameter values over time. The first 13 variables are about numerical records and relative geologic information for a simulation grid: 1. index: simulation index 2. i: the ith grid of x-axis 3. j: the ith grid of y-axis 4. k: the ith grid of z-axis 5. element_ref: element reference 6. nuft_ind: nuft index 7. x: grid location in the x axis direction 8. y: grid location in the y axis direction 9. z: grid location in the z axis direction 10. dx: grid length in the x axis direction 11. dy: grid length in the y axis direction 12. dz: grid length in the z axis direction 13. volume: volume of the simulation grid The remainder (14, 15, 16...) variables are the simulated parameter values over time, take Pressure as an example, are: 14. 0.0y: initial pressure per cell. 15. 10.0y: simulated pressure at the end of the 10th year. 16. 20.0y: simulated pressure at the end of the 20th year. ... (repeated for every 10 years until 200 years)... The model extends 10,000 m, 5,000 m and 1,411 m in the x,y and z dimensions, respectively. The mesh consists of 164,832 cells with mesh dimensions of 101 x 51 x 32 (nx, ny, nz), with cell dimensions ranging from 100 m laterally (along x and y-axis) and model layers are as designated in the z-axis: Layer 1: atmosphere (1e-30 m thick) Layer 2: upper caprock (10 m thick) Layers 3-13: Etchegoin (536.23 m thck) Layers 14-27: Macoma-Chanac (679.04 m thick) Layers 28-32: Santa Margarita-McLure (185.94 m thick) The wellbore is placed along node i=51, j=26, and extends vertically along 5 nodes from the top to the bottom of the model. Special instructions when extracting files: Each Gzip archive (ex: prod07.sim1000-sim00099.tar.gz) contains 100 simulations. Gzip archives should be transferred into base directories (ie. In Linux: mkdir prod07; mv prod07.*.tar.gz prod07/.) before extracting, or files will be overwritten. Each sub-simulation tree should have the following file structure pattern (using the linux 'tree' command): |-- prod07 | |-- sim0001 | |-- W31-0.2 | | |-- Q_brn.W31-0.2.sim0001.layers123.tab | | |-- Q_co2.W31-0.2.sim0001.layers123.tab | | |-- sim0001.W31-0.2.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.2.trans.pH.red.ntab | | |-- sim0001.W31-0.2.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.2.usnt.P.deltabg.red.ntab | |-- W31-0.5 | | |-- Q_brn.W31-0.5.sim0001.layers123.tab | | |-- Q_co2.W31-0.5.sim0001.layers123.tab | | |-- sim0001.W31-0.5.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.5.trans.pH.red.ntab | | |-- sim0001.W31-0.5.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.5.usnt.P.deltabg.red.ntab | |-- W31-1.0 | | |-- Q_brn.W31-1.0.sim0001.layers123.tab | | |-- Q_co2.W31-1.0.sim0001.layers123.tab | | |-- sim0001.W31-1.0.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-1.0.trans.pH.red.ntab | | |-- sim0001.W31-1.0.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-1.0.usnt.P.deltabg.red.ntab | |-- sim0001.mesh_k16.prod07.trans.genmsh Disclaimer This document was prepared as an account of work sponsored by an agency of the United States government. Neither the United States government nor Lawrence Livermore National Security, LLC, nor any of their employees makes any warranty, expressed or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States government or Lawrence Livermore National Security, LLC. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States government or Lawrence Livermore National Security, LLC, and shall not be used for advertising or product endorsement purposes. Lawrence Livermore National Laboratory is operated by Lawrence Livermore National Security, LLC, for the U.S. Department of Energy, National Nuclear Security Administration under Contract DE-AC52-07NA27344. This report was reviewed and released as LLNL-MI-753464.

aquifer↗
DOE OSTI2022

Data from: "Warming and the dependence of limber pine (Pinus flexilis) establishment on summer soil moisture within and above its current elevation range"

This data package contains data that were used for analysis in “Warming and the dependence of limber pine (Pinus flexilis) establishment on summer soil moisture within and above its current elevation range”, by Moyes et al. 2013. All data collection and field research were completed on Niwot Ridge, Colorado, USA.This data package contains nine comma-separated-values (.csv) files, one text (.txt) file, and two zipped seedling folders that were used for leaf area analysis. One zipped folder contains 468 .jpg photographs of seedlings, and the second contains 468 corresponding Image J-processed .jpg images that include silhouette leaf area values. .csv and .txt files can be opened using any compatible simple text-editor software such as TextEdit (Mac) and Notepad (Windows); .csv’s can also be opened using R and Microsoft Excel. Image files can be opened using Preview (Mac) and Photos (Windows). In addition, there are a total of 31 Microsoft Excel files: three .xlsx files, and 28 raw Li-Cor output .xls files. This data user’s guide is available in .pdf format, and can be opened using Adobe Acrobat Reader, or any other compatible file viewing software. Geospatial data showing field site locations are also included in the archive for use and reference. There are two geospatial formats in this archive: ESRI shapefiles (.shp) and keyhole markup-language (.kml) files. Both file types contain bounding box information, with the former being polygons, and the latter containing corner coordinates for each site. ESRI shapefiles can be opened using any geospatial software compatible with the file type (such as ESRI’s ArcGIS suite and QGIS), and .kml files are compatible with Google Earth and Google Maps.--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------Continued changes in climate are projected to alter the geographic distributions of plant species, in part by affecting where individuals can establish from seed. We tested the hypothesis that warming promotes uphill redistribution of subalpine tree populations by reducing cold limitation at high elevation and enhancing drought stress at low elevation. We seeded limber pine (Pinus flexilis) into plots with combinations of infrared heating and water addition treatments, at sites positioned in lower subalpine forest, the treeline ecotone, and alpine tundra. In 2010, first-year seedlings were assessed for physiological performance and survival over the snow-free growing season. Seedlings emerged in midsummer, about 5–8 weeks after snowmelt. Low temperature was not observed to limit seedling photosynthesis or respiration between emergence and October, and thus experimental warming did not appear to reduce cold limitation at high elevation. Instead, gas exchange and water potential from all sites indicated a prevailing effect of summer moisture stress on photosynthesis and carbon balance. Infrared heaters raised soil growing degree days (base 5 °C, p < 0.001) and August–September mean soil temperature (p < 0.001). Despite marked differences in vegetation cover and meteorological conditions across sites, volumetric soil moisture content (θ) at 5–10 cm below 0.16 and 0.08 m^3 m^(−3) consistently corresponded with moderate and severe indications of drought stress in midday stem water potential, stomatal conductance, photosynthesis, and respiration. Seedling survival was greater in watered plots than in heated plots (p = 0.01), and negatively related to soil growing degree days and duration of exposure to θ < 0.08 m^3 m^(−3) in a stepwise linear regression model (p < 0.0001). We concluded that seasonal moisture stress and high soil surface temperature imposed a strong limitation to limber pine seedling establishment across a broad elevation gradient, including at treeline, and that these limitations are likely to be enhanced by further climate warming.--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------Maintenance log:Amended file name error in Data User's Guide May 19 2022

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Ecohydrology of wetland plant communities along an estuarine to tidal river gradient

Abiotic filters that interact with wetland plant communities along tidal–fluvial gradients are highly dynamic, and understanding their quantitative thresholds and relationships to interspecific competition is important during an era of sea-level rise and watershed hydrologic change. Yet, landscape-scale studies of major coastal rivers from the river mouth to the head of tide, such as this study, remain rare. Here, we develop a new predictive framework for estuarine–tidal river research and management using a river-specific low-water datum and the wetland inundation indicator SEVg, the growing-season sum exceedance value of hourly surface-water depth. The distribution and variability of the wetland species pool (n = 203) on the 234 river kilometer (rkm) lower Columbia River and estuary floodplain are described for the first time. 4,940 quadrats at 50 marshes were surveyed (2005–2016). Throughout the estuarine–tidal river system, SEVg was well suited to describe the wetland inundation regime and its variability based on the combination of longitudinal river position and elevation. SEVg increased significantly landward. Two primary wetland inundation regimes were identified: the seaward-tidal, usually greatest during the winter months, and landward-fluvial, greatest during the growing season. Nearest the ocean, salinity is the abiotic factor limiting species richness and non-native species. Farther upriver, the daily, seasonal, and interannual variability of the wetting and drying cycle encourage disturbance-tolerant species and non-natives and limit the number of hydrophytes and total vegetative cover. Hence, the average between-year similarity of site-scale areal cover significantly decreased landward. Hierarchical cluster analysis indicated five vegetative groups and five ecohydrologic zones between rkm 0 and 234 were discriminated with 76 significant species–zone associations. All zones had unique indicator species. Species with high indicator values were Carex lyngbyei throughout the estuarine zones, and Eleocharis palustris, Sagittaria latifolia, and the invasive non-native Phalaris arundinacea in the upper estuarine and lower, middle, and upper tidal river zones (IV > 0.90). Competition from C. lyngbyei nearest the ocean and P. arundinacea in the tidal river was associated with reduced species richness when total cover was >65%. This framework of filters informs the design and prediction of future wetland plant communities on coastal river floodplains.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1996

Turbulent flame propagation in partially premixed flames

Turbulent premixed flame propagation is essential in many practical devices. In the past, fundamental and modeling studies of propagating flames have generally focused on turbulent flame propagation in mixtures of homogeneous composition, i.e. a mixture where the fuel-oxidizer mass ratio, or equivalence ratio, is uniform. This situation corresponds to the ideal case of perfect premixing between fuel and oxidizer. In practical situations, however, deviations from this ideal case occur frequently. In stratified reciprocating engines, fuel injection and large-scale flow motions are fine-tuned to create a mean gradient of equivalence ratio in the combustion chamber which provides additional control on combustion performance. In aircraft engines, combustion occurs with fuel and secondary air injected at various locations resulting in a nonuniform equivalence ratio. In both examples, mean values of the equivalence ratio can exhibit strong spatial and temporal variations. These variations in mixture composition are particularly significant in engines that use direct fuel injection into the combustion chamber. In this case, the liquid fuel does not always completely vaporize and mix before combustion occurs, resulting in persistent rich and lean pockets into which the turbulent flame propagates. From a practical point of view, there are several basic and important issues regarding partially premixed combustion that need to be resolved. Two such issues are how reactant composition inhomogeneities affect the laminar and turbulent flame speeds, and how the burnt gas temperature varies as a function of these inhomogeneities. Knowledge of the flame speed is critical in optimizing combustion performance, and the minimization of pollutant emissions relies heavily on the temperature in the burnt gases. Another application of partially premixed combustion is found in the field of active control of turbulent combustion. One possible technique of active control consists of pulsating the fuel flow rate and thereby modulating the equivalence ratio (Bloxsidge et al. 1987). Models of partially premixed combustion would be extremely useful in addressing all these questions related to practical systems. Unfortunately, the lack of a fundamental understanding regarding partially premixed combustion has resulted in an absence of models which accurately capture the complex nature of these flames. Previous work on partially premixed combustion has focused primarily on laminar triple flames. Triple flames correspond to an extreme case where fuel and oxidizer are initially totally separated (Veynante et al. 1994 and Ruetsch et al. 1995). These flames have a nontrivial propagation speed and are believed to be a key element in the stabilization process of jet diffusion flames. Different theories have also been proposed in the literature to describe a turbulent flame propagating in a mixture with variable equivalence ratio (Muller et al. 1994), but few validations are available. The objective of the present study is to provide basic information on the effects of partial premixing in turbulent combustion. In the following, we use direct numerical simulations to study laminar and turbulent flame propagation with variable equivalence ratio.

Poinsot, T.↗
DOE OSTI2022

Microstructural Characterization of the Second High Fluence Baffle-Former Bolt Retrieved from a Westinghouse Two-loop Downflow Type PWR

As one of the pressurized water reactor (PWR) internal components, baffle-former bolts (BFBs) are subjected to significant mechanical stress and neutron irradiation from the reactor core during the plant operation. Over the long operation period, these conditions lead to potential degradation and reduced load-carrying capacity of the bolts. In support of evaluating long-term operational performance of materials used in core internal components, the Oak Ridge National Laboratory (ORNL), through the Department of Energy (DOE), Light Water Reactor Sustainability (LWRS) Program, Materials Research Pathway (MRP) has harvested two high fluence BFBs from a commercial Westinghouse two-loop downflow type PWR. The two bolts of interest, i.e. bolts # 4412 and 4416, were withdrawn from service in 2011 as part of a preventative replacement plan. No identification of cracking or potential damage was found for these bolts during their removal in 2011. However, the bolts required a lower torque for removal from the baffle structure than the original torque specified during installation. Irradiation displacement damage levels in the bolts range from 15 to 41 displacements per atom. The goal of this project is to perform detailed microstructural and mechanical property characterization of BFBs following in-service exposures. The information from these bolts will be integral to the LWRS program initiatives in evaluating end of life microstructure and properties. Furthermore, valuable data will be obtained that can be incorporated into model predictions of long-term irradiation behavior and compared to results obtained in high flux experimental reactor conditions. In this report, we present our latest study in FY22 on microstructural characterizations of the second high fluence baffle-former bolt, i.e., bolt # 4412. Analytical electron microscopy and atom probe tomography characterization were performed. The radiation-induced defects in the material add to the large wealth of knowledge for neutron-induced defects in 304/316 grades of stainless steels, specifically for radiation-induced precipitation after high fluence commercial PWR irradiation. The main findings are summarized as follows: 1) The cavity size was considerably larger in the bolt thread section than in the bolt head, with the bolt thread section having a bimodal distribution of cavities greater than ~6 nm in diameter and less than ~3 nm in diameter. The bolt head only had the small-sized cavities. In addition, there was a denuded zone of large cavities near grain boundaries in the thread section of the bolt. 2) Radiation-induced precipitation in the BFB #4412 was highly complex, with the volume fraction, size, and number density of Ni/Si and Cu-rich precipitates depending strongly on the radiation temperature/dose. In many cases, co-precipitates of adjoined clusters were found with Ni/Si-rich precipitates sandwiched between Cu-rich clusters and Mo/Cr/P-rich clusters. 3) Solute segregation out of solution was highest for most solutes in the thread section of the bolt #4412 with the exception of Cu, which experienced more separation out of solution into Curich clusters in the bolt head section. This highlights the difference in the mechanisms for precipitation of Ni/Si clusters, which have the Ni 3 Si phase composition, and precipitation of Cu-rich clusters. 4) There appear to be multiple simultaneous influences that affect the microstructural variation along the length of the bolt that overcomes the ~2X difference in irradiation dose between the bolt head and the bolt thread. The irradiation temperature, thermal/fast neutron ratio variation, potential strain gradient, and exposure to PWR coolant water that each section of the bolt sees may have more influence on the microstructural evolution than the total irradiation dose. The bolt thread and shank, with higher temperature, higher relative fast neutron flux, higher strain, and exposure to coolant but lower dose, underwent more enhanced cavity formation, precipitation, and solute segregation than the bolt head section.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS2001

Proton Transport

The transport of protons across membranes is an essential process for both bioenergetics of modern cells and the origins of cellular life. All living systems make use of proton gradients across cell walls to convert environmental energy into a high-energy chemical compound, adenosine triphosphate (ATP), synthesized from adenosine diphosphate. ATP, in turn, is used as a source of energy to drive many cellular reactions. The ubiquity of this process in biology suggests that even the earliest cellular systems were relying on proton gradient for harvesting environmental energy needed to support their survival and growth. In contemporary cells, proton transfer is assisted by large, complex proteins embedded in membranes. The issue addressed in this Study was: how the same process can be accomplished with the aid of similar but much simpler molecules that could have existed in the protobiological milieu? The model system used in the study contained a bilayer membrane made of phospholipid, dimyristoylphosphatidylcholine (DMPC) which is a good model of the biological membranes forming cellular boundaries. Both sides of the bilayer were surrounded by water which simulated the environment inside and outside the cell. Embedded in the membrane was a fragment of the Influenza-A M$_2$ protein and enough sodium counterions to maintain system neutrality. This protein has been shown to exhibit remarkably high rates of proton transport and, therefore, is an excellent model to study the formation of proton gradients across membranes. The Influenza M$_2$ protein is 97 amino acids in length, but a fragment 25 amino acids long. which contains a transmembrane domain of 19 amino acids flanked by three amino acids on each side. is sufficient to transport protons. Four identical protein fragments, each folded into a helix, aggregate to form small channels spanning the membrane. Protons are conducted through a narrow pore in the middle of the channel in response to applied voltage. This channel is large enough to contain water molecules. and is normally filled with water. In analogy to the mechanism of proton transfer in some other channels, it has been postulated that protons are translocated along the network of water molecules filling the pore of the channel. This mechanism, however, must involve an additional important step because the channel contains four histidine amino acid residues, one from each of the helices, which are sufficiently large to occlude the pore and interrupt the water network. The histidine residues ensure channel selectivity by blocking transport of small ions, such as sodium or potassium. They have been also implicated in gating protons due to the ability of each histidine to become positively charged by accepting an additional proton. Two mechanisms of gating have been proposed. In one mechanism, all four histidines acquire an additional proton and, due to repulsion between their positive charges, move away from one another, thus opening the channel. The alternative mechanism relies of the ability of protons to move between different atoms in a molecule (tautomerization). Thus, a proton is captured on one side of the gate while another proton is released from the opposite side, and the molecule returns to the initial state through tautomerization. The simulations were designed to test these two mechanisms. Large-scale, atomic-level molecular dynamics simulations of the channel with the histidine residues in different protonation states revealed that all intermediate states of the system involved in the tautomerization mechanism are structurally stable and the arrangement of water molecules in the channel is conducive to the proton transport. In contrast, in the four-protonated state, postulated to exist in the gate-opening mechanism, the electrostatic repulsion between the histidine residues appears to be so large that the channel loses its structural integrity and one helix moves away from the remaining three. Additional information is contained within the original extended abstract.

Pohorille, Andrew↗
NASA NTRS2005

Depth-Penetrating Luminescence Thermography of Thermal- Barrier Coatings

A thermographic method has been developed for measuring temperatures at predetermined depths within dielectric material layers . especially thermal-barrier coatings. This method will help satisfy a need for noncontact measurement of through-the-thickness temperature gradients for evaluating the effectiveness of thermal- barrier coatings designed to prevent overheating of turbine blades, combustor liners, and other engine parts. Heretofore, thermography has been limited to measurement of surface and near-surface temperatures. In the thermographic method that is the immediate predecessor of the present method, a thermographic phosphor is applied to the outer surface of a thermal barrier coating, luminescence in the phosphor is excited by illuminating the phosphor at a suitable wavelength, and either the time dependence of the intensity of luminescence or the intensities of luminescence spectral lines is measured. Then an emissivity-independent surface-temperature value is computed by use of either the known temperature dependence of the luminescence decay time or the known temperature dependence of ratios between intensities of selected luminescence spectral lines. Until now, depth-penetrating measurements have not been possible because light of the wavelengths needed to excite phosphors could not penetrate thermal-barrier coating materials to useful depths. In the present method as in the method described above, one exploits the temperature dependence of luminescence decay time. In this case, the phosphor is incorporated into the thermal-barrier coat at the depth at which temperature is to be measured. To be suitable for use in this method, a phosphor must (1) exhibit a temperature dependence of luminescence decay time in the desired range, (2) be thermochemically compatible with the thermal-barrier coating, and (3) exhibit at least a minor excitation spectral peak and an emission spectral peak, both peaks being at wavelengths at which the thermal-barrier coating is transparent or at least translucent. Conventional thermographic phosphors are not suitable because they are most efficiently excited by ultraviolet light, which does not penetrate thermal- barrier coating materials. (Typical thermal-barrier coating materials include or consist of various formulations of yttria-stabilized zirconia.) Only a small fraction of phosphor candidates have significant excitation at wavelengths long enough (>500 nm) for sufficient penetration of thermal-barrier coatings. One suitable phosphor material . yttria doped with europium (Y2O3:Eu) . has a minor excitation peak at 532 nm and an emission peak at 611 nm. In experiments, this material was incorporated beneath a 100-microns-thick thermal-barrier coating and subjected to excitation and measurement by the luminescence-decay-time technique. These experiments were found to yield reliable temperature values up to 1,100 C. At the time of reporting the information for this article, a search for suitable phosphors other than (Y2O3:Eu) was continuing.

Eldridge, Jeffrey↗
NASA NTRS1999

Simulating the Formation of Lunar Floor-Fracture Craters Using Elastoviscoplastic Relaxation

Lunar floor-fracture craters formed during the height of mare basalt emplacement. Due to a general temporal and spatial relation with the maria, these craters, numbering some 200, may be diagnostic of the thermal structure of the crust during this time. As the name suggests, these craters exhibit brittle failure, generally limited to the central floor region. That, and a shallower depth than fresh lunar craters, has led to two main theories as to their formation: laccolith emplacement under the crater and viscous relaxation. The implications of each model for the state of the Moon's crust during this time are quite different, so the viability of each model must be checked. Laccolith emplacement has been treated elsewhere. However, previous attempts to study the relaxation of the craters have assumed only a uniform, Newtonian viscous response of the near surface to the topographic driving forces, and simply postulated that the fractures resulted from tensile stresses associated with floor uplift. Here, we use a more sophisticated rheological model that includes not only non-Newtonian viscous behavior (i.e., the viscosity is stress-dependent), but also incorporates elastic behavior and a plastic component to the rheology to directly simulate the formation of the floor fractures. The results of our simulations show that while elastoviscoplastic relaxation is potentially viable for larger floor-fracture craters, it is not viable for craters with diameters < or = 60 km, the size of the majority of floor-fracture craters. We employ the finite element method, a numerical technique well suited for boundary-value problems, via the commercially available MARC software package. To test the viability of topographic relaxation, our goal is to prepare the simulations as to maximize the amount of relaxation. We take advantage of the natural axisymmetry of craters, simulating one radial plane. Initial shapes are based on data for fresh craters from Pike. To simplify implementation, a fourth order polynomial is used for the basin, while a third order inverse function is used for the rim. This form closely approximates the long-wavelength behavior of complex craters, while ignoring higher-frequency topography, save the rim. This approximation is appropriate because crater relaxation is strongly controlled by long-wavelength topography. Loading is accomplished assuming a uniform gravity field (1.62 m/s-square) and a uniform density of 2900 kg/cubic m. The initial stress state is set to be hydrostatic, with an additional pressure term to account for any overlying topography. This additional pressure term is tapered exponentially with depth. While the simulations quickly settle on a preferred stress state, and while the final solution is fairly insensitive to the choice of the e-folding depth of the taper, selecting an e-folding depth close to the diameter of the crater sets the initial stress state near the preferred state. We assume a diuranally averaged surface temperature of -20 C, and allow temperature to increase with depth at a rate of 50 K//km. Assuming a thermal conductivity of 2 W/in/K, this gradient translates to a heat flow of 100mW/square m, an extremely high value for the Moon. Temperature, of course, will not increase without bound. To maximize relaxation, we allow our temperature profile to increase linearly until it reaches the solidus (assumed to be 1200C) at a depth of 24.4 km, at which point it is kept constant. The presence of melt will drop the bulk viscosity; however, we have no rheological control for partial melts. Therefore, we make no attempt to simulate this situation. Elastoviscoplastic rheological model. In general, geologic materials can behave in three main ways: elastically, viscously (via solid-state creep), and brittly (plasticity is a continuum approach to simulate this phenomenon). We combine these three deformation mechanisms in an extended Maxwell solid, where the total strain can be broken down into a simple summation of the elastic, creep, and plastic strains. In relaxation phenomena in general, the system takes advantage of any means possible to eliminate deviatoric stresses by relaxing away the topography. Previous analyses have only modeled the viscous response. Comparatively, the elastic response in our model can augment the relaxation, to a point. This effect decreases as the elastic response becomes stiffer; indeed, in the limit of infinite elastic Young's modulus (and with no plasticity), the solution converges on the purely viscous solution. Igneous rocks common to the lunar near-surface have Young's modulii in the range of 10-100 GPa. To maximize relaxation, we use a Young's modulus of 10 GPa. (There is negligible sensitivity to the other elastic modulus, the Poisson's ratio; we use 0.25.) For the viscous response, we use a flow law for steady-state creep in thoroughly dried Columbia diabase, because the high plagioclase (about 70 vol%) and orthopyroxene (about 17 vol%) content is similar to the composition of the lunar highland crust as described by remote sensing and sample studies: noritic anorthosite. This flow law is highly non-Newtonian, i.e., the viscosity is highly stress dependent. That, and the variability with temperature, stands in strong contrast to previous examinations of lunar floor-fracture crater relaxation. To model discrete, brittle faulting, we assume "Byerlee's rule," a standard geodynamical technique. We implement this "rule" with an-angle of internal friction of about 40 deg, and a higher-than-normal cohesion of about 3.2 MPa (to approximate the breaking of unfractured rock). The actual behavior of geologic materials is more complex than in our rheological model, so the uncertainties in the plasticity do not represent the state-of-the-art error. Additional information is contained in the original.

Dombard, A. J.↗
NASA NTRS1996

A technique for determining cloud free vs cloud contaminated pixels in satellite imagery

Since the first earth orbiting satellite sent pictures of the earth back to them, atmospheric scientists have been focused on the possibilities of using that information as both a forecasting tool and as a meteorological research tool. With the latest generation of Geostationary Operational Environmental Satellites (GOES) now entering service, that view of the earth yields views at a frequency and resolution never before available. These satellites have imagers with a five band multi-spectral capability with high spatial resolution. In addition, the sounder has eighteen thermal infrared (IR) channels plus one low-resolution visible band. With a resolution as small as one kilometer, GOES provides scientists with a powerful eye on the atmosphere. Menzel and Purdom (1994) detail both the imager and sounder capability as well as other systems on the GOES satellites. Immediately apparent in the visible channel are the patterns of clouds swirling over both oceans and continents. These clouds range in size from huge planetary systems covering thousands of kilometers to puffy fair weather cumulus clouds on the order of half a kilometer in size. With the IR sensors temperature patterns are observed. High clouds appear very cold, while low stratus field show temperatures near that of the surface. The surface, in turn, generally appears warmer than the clouds. It would seem then a simple manner to determine cloud and surface temperature from the imagery, but such is not the case. While most of the atmospheric constituents are well mixed and homogeneous, water vapor is not. The water molecule, because of its unique structure and vibration modes, affects the transmittance of the atmosphere most notably in the infrared regions. There are regions of the IR spectrum where water vapor acts as a strong absorber, and at others it is nearly transparent. The transparent wavelengths are called windows, and one such window occurs at 11.2 microns. Adjacent to this window at 12.7 microns which is strongly absorbed by water vapor. These two wavelengths form what is known as a split window, the utility of which was used. Using the linearized form of the radiative transfer equation, they were able to use the split window to determine the amount of water vapor present in the atmosphere. Jedlovec developed the physical split-window (PSW) technique which determines the integrated water content (IWC). The PSW method using Visible Infrared Spin Scan Radiometer (VISSR) Atmospheric Sounder (VAS) found on the older versions of the GOES satellites was used. Recently, Jedlovec and colleagues have been attempting to apply the PSW method using full disk IR imagery obtained by the new generation of GOES satellites. IWC is essential for improved analysis and prediction of convective storms which have been observed to develop in regions of both strong and rapidly evolving moisture gradients. It has also been used in the prediction of clouds and precipitation.

Wohlman, Richard A.↗
NASA NTRS2003

The Galileo Probe: How it Has Changed Our Understanding of Jupiter

The Galileo Mission to Jupiter, which arrived in December of 1995, provided the first study by an orbiter, and the first in-situ sampling via an entry probe, of an outer planet atmosphere. The rationale for an entry probe is that, even from an orbiter, remote sensing of the jovian atmosphere could not adequately retrieve the information desired. This paper provides a current summary of the most significant aspects of the data returned from the Galileo entry probe. As a result of the probe measurements, there has been a reassessment of our understanding of outer planet formation and evolution of the solar system. The primary scientific objective of the Galileo probe was to determine the composition of the jovian atmosphere, which from remote sensing remained either very uncertain, or completely unknown, with respect to several key elements. The probe found that the global He mass fraction is. significantly above the value reported from the Voyager Jupiter flybys but is slightly below the protosolar value, implying that there has been some settling of He to the deep jovian interior. The probe He measurements have also led to a reevaluation of the Voyager He mass fraction for Saturn, which is now determined to be much closer to that of Jupiter. The elements C, N, S, Ar, Kr, Xe were all found to have global abundances approximately 3 times their respective solar abundances. This result has raised a number of fundamental issues with regard to properties of planetesimals and the solar nebula at the time of giant planet formation. Ne, on the other hand, was found to be highly depleted, probably as the result of it being carried along with helium as helium settles towards the deep interior. The global abundance of O was not obtained by the probe because of the influence of local processes at the probe entry site (PES), processes which depleted condensible species, in this case H2O, well below condensation levels. Other condensible species, namely NH3 and H2S, were similarly affected but attained their deep equilibrium mixing ratios before the maximum depth sampled by the probe. Processes that might be capable of producing such effects on the condensibles are still under investigation. Measured isotopic ratios of noble gases and other heavy elements are solar, and (D + (Sup 3)He)/H is the same to within measurement uncertainties as in the local interstellar medium. No thick clouds were detected, and in particular no significant water cloud, but the PES location clearly affected the probe measurements of clouds. In fact, the probe data must be understood in the context of the location of the PES, which was within what is termed a 5 micron hot spot, a local clearing in the clouds that is bright near the 5 microns spectral region. The thermal structure at the PES was determined from approximately 1000 km above the 1 bar pressure level (10(exp -9 bars)) to 132 km 1 bar (22bars). The probe showed the atmosphere to have a generally sub-adiabatic temperature gradient (static stability) of = 0.1 K/km to as deep as the probe made measurements. In the upper atmosphere the probe derived a maximum positive vertical temperature gradient of approximately 5 K/km, and maximum temperature of = 900 K. The energy sources producing the warm upper atmosphere have yet to be completely identified. At first glance, Doppler tracking of the probe indicates that the long observed cloud level zonal winds extend to levels at least as deep as the probe made measurements. Zonal wind increases from = 80 m/s at pressures less than a bar to about 180 m/s near 5 bars, and remains approximately constant with depth thereafter. However, there is a question as to whether the winds measured from probe tracking are representative of the general wind field, or are considerably influenced by localized winds associated with the PES.

Young, Richard E.↗
NASA NTRS2011

Spatially Resolved, In Situ Carbon Isotope Analysis of Archean Organic Matter

Spatiotemporal variability in the carbon isotope composition of sedimentary organic matter (OM) preserves information about the evolution of the biosphere and of the exogenic carbon cycle as a whole. Primary compositions, and imprints of the post-depositional processes that obscure them, exist at the scale of individual sedimentary grains (mm to micron). Secondary ion mass spectrometry (SIMS) (1) enables analysis at these scales and in petrographic context, (2) permits morphological and compositional characterization of the analyte and associated minerals prior to isotopic analysis, and (3) reveals patterns of variability homogenized by bulk techniques. Here we present new methods for in situ organic carbon isotope analysis with sub-permil precision and spatial resolution to 1 micron using SIMS, as well as new data acquired from a suite of Archean rocks. Three analytical protocols were developed for the CAMECA ims1280 at WiscSIMS to analyze domains of varying size and carbon concentration. Average reproducibility (at 2SD) using a 6 micron spot size with two Faraday cup detectors was 0.4 %, and 0.8 % for analyses using 1 micron and 3 micron spot sizes with a Faraday cup (for C-12) and an electron multiplier (for C-13). Eight coals, two ambers, a shungite, and a graphite were evaluated for micron-scale isotopic heterogeneity, and LCNN anthracite (delta C-13 = -23.56 +/- 0.1 %, 2SD) was chosen as the working standard. Correlation between instrumental bias and H/C was observed and calibrated for each analytical session using organic materials with H/C between 0.1 and 1.5 (atomic), allowing a correction based upon a C-13H/C-13 measurement included in every analysis. Matrix effects of variable C/SiO2 were evaluated by measuring mm to sub-micron graphite domains in quartzite from Bogala mine, Sri Lanka. Apparent instrumental bias and C-12 count rate are correlated in this case, but this may be related to a crystal orientation effect in graphite. Analyses of amorphous Archean OM suggest that instrumental bias is consistent for 12C count rates as low as 10% relative to anthracite. Samples from the ABDP-9 (n=3; Mount McRae Shale, approximately 2.5 Ga), RHDH2a (n=2; Carrawine Dolomite and Jeerinah Fm, approximately 2.6 Ga), WRL1 (n=3; Wittenoom Fm, Marra Mamba Iron Formation, and Jeerinah Fm, approximately 2.6 Ga), and SV1 (n=1; Tumbiana Fm, approximately 2.7 Ga) drill cores, each previously analyzed for bulk organic carbon isotope composition, yielded 100 new, in situ data from Neoarchean sedimentary OM. In these samples, delta C-13 varies between -53.1 and -28.3 % and offsets between in situ and bulk compositions range from -8.3 to 18.8%. In some cases, isotopic composition and mode of occurrence (e.g. morphology and mineral associations) are statistically correlated, enabling the identification of distinct reservoirs of OM. Our results support previous evidence for gradients of oxidation with depth in Neoarchean environments driven by photosynthesis and methane metabolism. The relevance of these findings to questions of bio- and syngenicity as well as the alteration history of previously reported Archean OM will be discussed.

Williford, Kenneth H.↗
DOE OSTI2020

Progression to Compatibility Evaluations in Flowing Molten Salts

Molten salt compatibility with structural alloys has been identified as a key issue for the development of Generation 3 concentrating solar power (CSP) systems with thermal storage. To accelerate this evaluation to pumped systems, the goal of this project was to conduct thermal convection loop (TCL) experiments with a peak temperature of ≥700°C and a typical temperature gradient of ~100°C. The experiments indicated that conventional ~16wt.%Cr Ni-based alloys are compatible up to 700°C with purified (i.e. low O) or dried (low H₂O) industrial-sourced Mg-K-Na chloride salt with <10 μm/yr loss. This two-year project was conducted based on the experimental and mechanistic understanding developed more than 60 years ago at Oak Ridge National Laboratory (ORNL). The first TCL experiment met the <15 μm/yr corrosion metric for this project with specimens of Ni-based alloy 600 exposed at 580°-700°C for 1000 h and post-exposure room temperature tensile tests showed minimal degradation. The second TCL experiment successfully deployed an electrochemical sensor from Argonne National Laboratory and had a peak temperature of 750°C but only ran for ~110 h due to a furnace failure. Both experiments used highly purified industrial-sourced salt with an O content of ~3 μg O/g salt and a Mg addition of 0.04 wt.%. In the second year, the Chloride Collective developed a more economical drying procedure such that the O content was much higher (>20,000 μg O/g salt). A third TCL experiment was conducted with a sensor and a peak temperature of 700°C using dried salt from the same industrial source and increased NaCl content (~20 wt.%). In addition to a 0.05% Mg addition to the salt, a Mg coupon was added in the coldest part of the loop which dissolved during the experiment. Again, small mass changes were noted for specimens of alloys 600 and C276 but the values were slightly higher than those measured in purified salt with Mg. A thin non-continuous and non-adherent oxide layer was deposited on most specimens containing Mg, Si and Al but Cr depletion also was observed. Both years included facilities qualification crucible experiments and then capsule experiments to confirm a baseline isothermal reaction rate. The first year capsule experiment led to the conclusion that the two-stage ORNL purification process using NH₄Cl and CCl₄ left the salt with a high Cl potential and a Mg addition (~0.05wt.%) was needed to lower the potential. The second year capsule experiments at 600° and 700°C found little difference in depth of attack for 0-0.25%Mg additions. In general, the complex reactions where salt can be trapped in the porous surface layer of metal indicated that mass change is an unreliable metric and average depth of Cr depletion is a better metric for assessing the extent of attack. These results have created a new baseline that is contrary to the recent published literature for chloride salts where mass losses have been reported that can be extrapolated to very significant annual metal loss rates. Chloride salt corrosion can be controlled and the results also indicate that salt purification to low O levels may not be necessary. However, additional TCL experiments are needed at different times and temperatures to generate important engineering information such as temperature dependent corrosion rates and reaction rate laws for extrapolation to long-term behavior and isolate the effect of salt additives and impurities.

36 MATERIALS SCIENCE↗
DOE OSTI2025

GEPA: Reflective Prompt Evolution Can Outperform Reinforcement Learning

Large language models (LLMs) are increasingly adapted to downstream tasks via reinforcement learning (RL) methods like Group Relative Policy Optimization (GRPO), which often require thousands of rollouts to learn new tasks. We argue that the interpretable nature of language often provides a much richer learning medium for LLMs, compared to policy gradients derived from sparse, scalar rewards. To test this, we introduce GEPA (Genetic-Pareto), a prompt optimizer that thoroughly incorporates natural language reflection to learn high-level rules from trial and error. Given any AI system containing one or more LLM prompts, GEPA samples trajectories (e.g., reasoning, tool calls, and tool outputs) and reflects on them in natural language to diagnose problems, propose and test prompt updates, and combine complementary lessons from the Pareto frontier of its own attempts. As a result of GEPA's design, it can often turn even just a few rollouts into a large quality gain. Across six tasks, GEPA outperforms GRPO by 6% on average and by up to 20%, while using up to 35x fewer rollouts. GEPA also outperforms the leading prompt optimizer, MIPROv2, by over 10% (e.g., +12% accuracy on AIME-2025), and demonstrates promising results as an inference-time search strategy for code optimization. We release our code at https://github.com/gepa-ai/gepa.

97 MATHEMATICS AND COMPUTING↗