Search NASA⌕ Search

SEARCH · Search NASA

Results for “surface pressure”

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

Surface Chemistry of Liquid Bismuth under Oxygen and Water Vapor Studied by Ambient Pressure X-ray Photoelectron Spectroscopy

Bismuth (Bi) has a fairly low melting point of 544 K making it practical as a liquid metal medium in a number of applications. Under ambient atmospheric conditions the surface of solid Bi oxidizes. While the solid Bi interface has been fairly well characterized with surface science studies upon exposure to oxygen, to date little is known about the molecular level reactivity of the liquid interface. Here, using ambient pressure X-ray photoelectron spectroscopy the liquid-gas interfacial chemistry of Bi was examined upon exposure to oxygen gas and water vapor at 550 K up to a maximum pressure of 0.5 Torr. Water vapor remained unreactive towards the liquid interface over the entire pressure range. Oxygen remined unreactive up to 10 –4 Torr, whereas above this pressure oxidation was observed forming Bi 2 O 3 . The oxidation exposure onset was > 5 x 10 5 Langmuirs, significantly higher than what is required for solid Bi interfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Pressurized Gas-Driven Elemental Redistribution Enables Ultrastable PtNi Catalysts for Heavy-Duty Vehicles

Gas-driven element redistribution, characterized by the preferential enrichment of one element at the surface relative to the bulk, is frequently observed in multicomponent alloys. Using L1 0 -ordered PtNi as a model system, we reveal that gas pressure plays a critical role in governing adsorption-driven surface composition during annealing in reducing gases: low pressure favors Pt surface segregation, while high pressure facilitates Ni surface enrichment. In this study, we developed a high-pressure nitriding (HPN) strategy that modulates the surface structure and composition of PtNi catalysts. The resulting HPN-PtNi exhibits enhanced performance and durability in membrane electrode assemblies for heavy-duty fuel cell applications, maintaining a high current density of 1.19 A cm −2 at 0.7 V after 90,000 voltage cycles. Through a combination of experimental and theoretical analyses, we reveal that the HPN process forms additional stabilizing Ni−N bonds and induces elemental redistribution with Ni surface enrichment and a Ni-deficient Pt subsurface. These modifications alter the atomic coordination environment of the ordered PtNi phase. This work presents a generalizable strategy to design robust and high-performing Pt-based catalysts by controlling gas-pressure-driven elemental redistribution and dopant incorporation.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2024

Impact of breech geometry and propellant flow on the release of large pellets for the ITER disruption mitigation system

Studies have been performed on the release mechanism for large pellets using high pressure gas in a shattered pellet injector. Typically, pellets are dislodged from the cryogenic surface and accelerated down a barrel using high pressure gas delivered by a fast-acting propellant valve. The pellets impact an angled surface which shatters the pellet into many small fragments before entering the plasma. This technique was initially demonstrated on DIII-D (Commaux et al 2016 Nucl. Fusion 56 046007) and is now deployed on JET, KSTAR, ASDEX-Upgrade, and other tokamaks around the world in support of ITER's disruption mitigation system design and physics basis. The large hydrogen, 28.5 mm diameter, 2 length-to-diameter ratio, pellets foreseen for ITER SPI operation have low material strength and low heat of sublimation, which cause the pellets to be fragile and highly reactive to the impact of warm propellant gas. Due to the size of the pellets, significantly more propellant gas is required to dislodge and accelerate them. This creates a potentially significant propellant gas removal issue as 2–6 bar-L of gas is expected to be required for release and speed control. The research presented in this paper is an in-depth exploration of the parameters that are keys to reliable pellet release and speed control. Computational fluid dynamics (CFD) modeling of propellant flows through various breech designs was conducted to determine the force generated on the back surface of a pellet. These simulations assumed the use of the ORNL designed flyer plate valve. CFD modeling combined with experimental measurements provide adequate insight to determine a path to an optimal valve and breech design for ITER SPI pellet release and speed control while minimizing propellant gas usage.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Failure in Confined Brazilian Tests on Sandstone

Strength of rocks in the confined tension region, where the minimum principal stress is tensile, has only infrequently been measured and is not well understood. Quasi-static confined Brazilian tests under a range of confining stresses (2.76 to 27.58 MPa) where used to determine the strength of sandstone in the confined tension region. The test results indicate that the strength in the confined tension region was a strong function of the intermediate principal stress: increasing the intermediate principal stress significantly increased the strength of the sandstone. The strength data were well fit by the Mogi–Coulomb criterion, which accounts for the intermediate principal stress. Unconfined Brazilian strength data were not well fit to the Mogi–Coulomb criterion derived from the confined Brazilian test data, consistent with a transition from tensile to shear processes dominating failure with increasing confining pressure. Observations of post-failure fracture surfaces reveal more indication of shear processes with increasing confining pressure. Numerical simulations from combined finite-discrete element method are compared to the experimental results and reflect similar conditions for failure compared to the experimental tests in the confined tension region.

42 ENGINEERING↗
DOE OSTI2021

Model, data, and code for paper "Modeling of streamflow in a 30-kilometer-long reach spanning 5 years using OpenFOAM 5.x"

The data package includes data, model, and code that support the analyses and conclusions in the paper titled “modeling of streamflow in a 30-kilometer-long reach spanning 5 years using OpenFOAM 5.x”. The primary goal of this paper is to demonstrate that key streamflow properties such as water depth, flow velocity, and dynamic pressure in a natural river at 30-kilometer scale over 5 years can be reliably and efficiently modeled using the computational framework presented in this paper. To support the paper, various data types from remote sensing, field observations, and computational models are used. Specific details are described as follows. Firstly, the river bathymetry data was obtained from a Light Detection and Ranging (LiDAR) survey. This data is then converted to a triangulated surface format, STL, for mesh generation in OpenFOAM. The STL data can be found in Model_Setups/BaseCase_2013To2015/constant/triSurface. The OpenFOAM mesh generated using this STL file can be found in constant/polyMesh. Other model setups, boundary and initial conditions can be found in /system and /0.org under folder BaseCase_2013To2015. A similar data structure can also be found in BaseCase_2018To2019 for the simulations during 2018 and 2019. Secondly, the OpenFOAM simulations need the upstream discharge and water depth information at the upstream boundary to drive the model. These data are generated from a one-dimensional hydraulic model and the data can be found under the folder Model_Setups /1D model Mass1 data. The mass1_65.csv and mass1_191.csv files include the results of the 1D model at the model inlet and outlet, respectively. The Matlab source code Mass1ToOFBC20182019.m is used to convert these data into OpenFOAM boundary condition setups.With the above OpenFOAM model, it can generate data for water surface elevation, flow velocity, and dynamic pressure. In this paper, the water surface elevation was measured at 7 locations during different periods between 2011 and 2019. The exact survey locations (see Fig1_SurveyLocations.txt) can be found in folder Fig_1. The variation of water stage over time at the 7 locations can be found in folder /Observation_WSE. The data type include .txt, .csv, .xlsx, and .mat. The .mat data can be loaded by Matlab.We also measured the flow velocities at 12 cross-sections along the river. At each cross-section, we recorded the x, y locations, depth, three velocity components u,v,w. These data are saved to a Matlab format which can be found under folder /Observation_Velocity and /Fig_1. The relative locations of velocity survey locations to the river bathymetry can be found in Figure 1c.The water stage data at the 7 locations from OpenFOAM, 1D, and 2D hydraulic models are also provided to evaluate the long-term performance of 3D models vs 1D/2D models. The water stage data for the 7 locations from OpenFOAM have been saved to .mat format and can be found in /OpenFOAM_WSE. The water stage data from the 1D model are saved in .csv format and can be found in /Mass1_WSE. The water stage from the 2D model is saved as .mat format and can be found in / Mass2_WSEIn addition, the OpenFOAM model outputs the information of hydrostatic and hydrodynamic pressure. They are saved as .mat format under folder /Fig_11/2013_1. As the files are too large, we only uploaded the data for January 2013. The area of different ratio of dynamic pressure to static pressure for all simulation range, i.e., 2013-2015, are saved to .mat format. They can be found in /Fig_11/PA. Further, the data of wall clock time versus the solution time of the OpenFOAM modeling are also saved to .mat format under folder /Fig_13/LogsMat. In summary, the data package contains seven data types, including .txt, .csv, .xlsx, .dat, .stl, .m, and .mat. The former 4 types can be directly open using a text editor or Microsoft Office. The .mat format needs to be read by Matlab. The Matlab source code .m files need to be run with Matlab. The OpenFOAM setups can be visualized in ParaView. The .stl file can be opened in ParaView or Blender. The data in subfolders Fig_1 to Fig_10 and Fig_12 are copied from the aforementioned data folders to generate specific figures for the paper. A readME.txt file is included in each subfolder to further describe how the data in each folder are generated and used to support the paper.Please use the data package's DOI to cite the data package. Please contact yunxiang.chen@pnnl.gov if you need more data related to the paper.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Probing the surface chemistry for reverse water gas shift reaction on Pt(1 1 1) using ambient pressure X-ray photoelectron spectroscopy

Using ambient pressure XPS (APXPS), we explored carbon dioxide (CO 2 ) adsorption and CO 2 hydrogenation on Pt(1 1 1) single crystal surface to observe the activation of CO 2 and the subsequent reaction mechanism. In pure CO 2 , we observed CO adsorbates and adsorbed oxygen on Pt(1 1 1) derived from CO 2 dissociation at room temperature. The introduction of H 2 (at a pressure ratio of 1:1 (H 2 :CO 2 )) increased the production of CO across all temperatures by facilitating the removal of surface oxygen. As a consequence, the surface could expose sites that could then be utilized for producing CO. Under these conditions, the reverse water–gas shift (RWGS) reaction was observed starting at 300 °C. At higher H 2 partial pressure (10:1 (H 2 :CO 2 )), the RWGS reaction initiated at a lower temperature of 200 °C and continued to enhance the conversion of CO 2 with increasing temperatures. Our results revealed that CO 2 was activated on a clean Pt(1 1 1) surface through the dissociation mechanism to form adsorbed CO and O at room temperature and at elevated temperatures. Introducing H 2 facilitated the RWGS as adsorbed oxygen was consumed continuously to form H 2 O, and adsorbed CO desorbed from the surface at elevated temperatures. This work clearly provides direct experimental evidence for the surface chemistry of CO 2 dissociation and demonstrates how hydrogen impacts the RWGS reaction on a platinum surface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Revealing CO 2 dissociation pathways at vicinal copper (997) interfaces

Size- and shape-tailored copper (Cu) nanocrystals can offer vicinal planes for facile carbon dioxide (CO 2 ) activation. Despite extensive reactivity benchmarks, a correlation between CO 2 conversion and morphology structure has not yet been established at vicinal Cu interfaces. Herein, ambient pressure scanning tunneling microscopy reveals step-broken Cu nanocluster evolutions on the Cu(997) surface under 1 mbar CO 2 (g). The CO 2 dissociation reaction produces carbon monoxide (CO) adsorbate and atomic oxygen (O) at Cu step-edges, inducing complicated restructuring of the Cu atoms to compensate for increased surface chemical potential energy at ambient pressure. The CO molecules bound at under-coordinated Cu atoms contribute to the reversible Cu clustering with the pressure gap effect, whereas the dissociated oxygen leads to irreversible Cu faceting geometries. Synchrotron-based ambient pressure X-ray photoelectron spectroscopy identifies the chemical binding energy changes in CO-Cu complexes, which proves the characterized real-space evidence for the step-broken Cu nanoclusters under CO(g) environments. Our in situ surface observations provide a more realistic insight into Cu nanocatalyst designs for efficient CO 2 conversion to renewable energy sources during C 1 chemical reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Surface Interactions and Nanoconfinement of Methane and Methane plus CO2 Revealed by High-Pressure Magic Angle Spinning NMR Spectroscopy and Molecular Dynamics

This study explores the fundamental, molecular- to microscopic-level behavior of methane gas confined into nanoporous silica proxies with different pore diameters and surface-to-volume (S/V) ratios. Surfaces and pore walls of nanoporous silica matrices are decorated with hydroxyl (-OH) groups, resembling natural heterogeneity. High-pressure MAS NMR was utilized to characterize the interactions between methane and the engineered nanoporous silica proxies under various temperature and pressure regimes. There was a change in the chemical shift position of confined methane slightly in the mixtures with nanoporous silica up to 393 K, as shown by high-pressure 13C-NMR. The 13C-NMR chemical shift of methane was changed by pressure, explained by the densification of methane inside the nanoporous silica materials. The influence of pore diameter and S/V of the nanoporous silica materials on the behaviors and dynamics of methane were studied. The presence of CO2 in mixtures of silica and methane needs analysis with caution because CO2 in a supercritical state and gaseous CO2 change the original structure of nanoporous silica and change surface area and pore volume. According to simulation, the picosecond scale dynamics of methane confined in larger pores of amorphous silica is faster. In the 4 nm pore, the diffusivity obtained from MD simulations in the pore with a higher S/V ratio is slower due to the trapping of methane molecules in adsorbed layers close to the corrugated pore surface. In contrast, relaxation measured with NMR for smaller pores (higher S/V) exhibits larger T1, indicating slower relaxation.

03 NATURAL GAS↗
DOE OSTI2023

A review of in situ/operando studies of heterogeneous catalytic hydrogenation of CO 2 to methanol

Repurposing CO 2 into chemicals, one of the Carbon Dioxide Removal (CDR) strategies, still faces significant challenges in conversion and energy efficiency due to the lack of effective catalysts and processes. Fundamental understanding through in situ/operando investigations of the reaction mechanisms and catalyst structures is pivotal for developing the efficient catalysts. This paper reviews the past and recent in situ/operando studies of methanol synthesis from heterogeneous CO 2 hydrogenation over a few typical catalysts including Cu-based, oxide-based, and noble-metal-based catalysts. With the development of high-pressure reactors, in situ/operando IR, X-ray absorption spectroscopy (XAS), X-ray diffraction (XRD), neutron diffraction and imaging have been used to reveal the surface intermediates and structures of working catalysts under CO 2 hydrogenation conditions. On the one hand, the combined operando techniques shed light on working mechanisms for some catalytic systems. On the other hand, due to the complexity of selective CO 2 hydrogenation reaction and limited sensitivity of current accessible operando techniques to the surface structure of catalysts, it is still murky about the exact nature of the active sites at the surface/interface and how different sites promote different reaction paths. Furthermore, it is concluded that new methodologies for differentiating signals from different sites, and the development of surface-sensitive techniques for high pressure reactions are needed for providing structural descriptors of highly active and selective CO 2 conversion catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Convection in Scaled Turbine Internal Cooling Passages With Additive Manufacturing Roughness

Abstract Additive manufacturing processes, such as direct metal laser sintering (DMLS), enable the creation of novel turbine cooling internal passages and systems. However, the DMLS method produces a significant and unique surface roughness. Previous work in scaled passages analyzed pressure losses and friction factors associated with the rough surfaces, as well as investigated the velocity profiles and turbulent flow characteristics within the passage. In this study, the heat transfer characteristics of scaled additively manufactured surfaces were measured using infrared (IR) thermography. Roughness panels were CNC machined from plates of aluminum 6061 to create near isothermal roughness elements when heated. Fluid resistance differences between the aluminum roughness panels and roughness panels constructed from ABS plastic using the same roughness patterns from McClain et al. (2020, “Flow in a Simulated Turbine Blade Cooling Channel With Spatially Varying Roughness Caused by Additive Manufacturing Orientation,” ASME Turbo Expo 2020, Turbomachinery Technical Conference and Exposition, Virtual Conference, Sept. 21–25, GT2020-16069) were investigated. Finally, the overall thermal performance enhancements and friction losses were assessed through the calculation of surface averaged “global thermal performance” ratios. The global thermal performance characterizations indicate results in-line with those found for traditional commercial roughness and slightly below traditional internal passage convection enhancement methods such as swirl chambers, dimples, and ribs. The passages investigated in this study do not include compressibility effects or the long-wavelength artifacts and channel geometric deviations observed by Wildgoose et al. (2020, “Impact of Additive Manufacturing on Internal Cooling Channels with Varying Diameters and Build Directions,” ASME Turbo Expo 2020, Turbomachinery Technical Conference and Exposition, Virtual Conference, Sept. 21–25, GT2020-15049). However, the results of this study indicate that, based on the roughness augmentation alone, artificial convective cooling enhancers such as turbulators or dimples may still be required for additively manufactured turbine component cooling.

Engineering↗
DOE OSTI2021

Numerical study to predict the effect of surface roughness on the thermal and hydraulic performance of additively manufactured heat exchangers

Heat exchangers are used in many industries to transfer thermal energy between two or more objects with an associated thermal and hydraulic performance indicated by thermal effectiveness and pressure drop respectively. Many investigations have aimed to enhance the thermal effectiveness of heat exchangers by introducing turbulent flow inducers that increase the heat transfer rate at the expense of increasing pressure drop. Flow inducers like surface roughness can be controlled with advanced manufacturing methods like additive manufacturing for heat exchanger production. Selecting and manufacturing an optimal surface roughness for heat exchangers requires a more comprehensive approach that accounts for the surface roughness effect on the thermal and hydraulic performance of a sample. Since fabricating samples with additive manufacturing to a desired surface roughness can be an extensive experimental process, this study will fundamentally investigate the effect of surface roughness on the thermal and velocity boundary layers for various surface roughness profiles with physics-based computational modeling. To accurately model the surface roughness effect on thermohydraulic performance, correlations like roughness ratio, friction factor, and Nusselt number are investigated. This study expects to indicate correlations for different flow regimes which can successfully predict the effect of surface roughness on thermohydraulic performance. Results generated from modeling surface roughness effect with unique correlations are then compared to experimental data obtained from literature for validation. The results of this study will provide insight to understand the thermal fluid physics associated with surface roughness and design guidelines for heat exchangers produced with additive manufacturing.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Structure of pressure-gradient-driven current singularity in ideal magnetohydrodynamic equilibrium

Abstract Singular currents typically appear on rational surfaces in non-axisymmetric ideal magnetohydrodynamic (MHD) equilibria with a continuum of nested flux surfaces and a continuous rotational transition. These currents have two components: a surface current (Dirac δ -function in flux surface labeling) that prevents the formation of magnetic islands, and an algebraically divergent Pfirsch–Schlüter current density when a pressure gradient is present across the rational surface. On flux surfaces adjacent to the rational surface, the traditional treatment gives the Pfirsch–Schlüter current density scaling as J ∼ 1 / Δ ι , where Δ ι is the difference of the rotational transform relative to the rational surface. If the distance s between flux surfaces is proportional to Δ ι , the scaling relation J ∼ 1 / Δ ι ∼ 1 / s will lead to a paradox that the Pfirsch–Schlüter current is not integrable. In this work, we investigate this issue by considering the pressure-gradient-driven singular current in the Hahm–Kulsrud–Taylor problem, which is a prototype for singular currents arising from resonant magnetic perturbations. We show that not only the Pfirsch–Schlüter current density but also the diamagnetic current density are divergent as ∼ 1 / Δ ι . However, due to the formation of a Dirac δ -function current sheet at the rational surface, the neighboring flux surfaces are strongly packed with s ∼ ( Δ ι ) 2 . Consequently, the singular current density J ∼ 1 / s , making the total current finite, thus resolving the paradox. Furthermore, the strong packing of flux surfaces causes a steepening of the pressure gradient near the rational surface, with ∇ p ∼ d p / d s ∼ 1 / s . In general non-axisymmetric MHD equilibrium, contrary to Grad’s conjecture that the pressure profile is flat around densely distributed rational surfaces, our result suggests a pressure profile that densely steepens around them.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Probing FeSi, a d -electron topological Kondo insulator candidate, with magnetic field, pressure, and microwaves

Recently, evidence for a conducting surface state (CSS) below 19 K was reported for the correlated d -electron small gap semiconductor FeSi. In the work reported herein, the CSS and the bulk phase of FeSi were probed via electrical resistivity ρ measurements as a function of temperature T , magnetic field B to 60 T, and pressure P to 7.6 GPa, and by means of a magnetic field-modulated microwave spectroscopy (MFMMS) technique. The properties of FeSi were also compared with those of the Kondo insulator SmB 6 to address the question of whether FeSi is a d -electron analogue of an f -electron Kondo insulator and, in addition, a “topological Kondo insulator” (TKI). The overall behavior of the magnetoresistance of FeSi at temperatures above and below the onset temperature T S = 19 K of the CSS is similar to that of SmB 6 . The two energy gaps, inferred from the ρ( T ) data in the semiconducting regime, increase with pressure up to about 7 GPa, followed by a drop which coincides with a sharp suppression of T S . Several studies of ρ( T ) under pressure on SmB 6 reveal behavior similar to that of FeSi in which the two energy gaps vanish at a critical pressure near the pressure at which T S vanishes, although the energy gaps in SmB 6 initially decrease with pressure, whereas in FeSi they increase with pressure. The MFMMS measurements showed a sharp feature at T S ≈ 19 K for FeSi, which could be due to ferromagnetic ordering of the CSS. However, no such feature was observed at T S ≈ 4.5 K for SmB 6 .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Surface Treatment Strategies for Mitigating Gas Hydrate & Asphaltene Formation, Growth, and Deposition in Flowloops

Abstract Numerous solids including gas hydrates, waxes, and asphaltenes have the potential to form in the production lines of gas and oil fields. This creates a highly non-ideal scenario as the accumulation of said species leads to flow assurance issues, especially with long-term processes like deposition. Since an ever-increasing amount of material is deposited in place at the pipe surface, production stoppage or active mitigation efforts become inevitable. The latter production issues result in increased safety risks and operational expenditures. Therefore, a cost-effective, passive deposition mitigation technology, such as a pipeline coating or surface treatment is especially appealing. The ability to address multiple pipeline flow assurance issues simultaneously without actively disrupting production would represent a dramatic step forward in this area. This study is part of a long-term ongoing effort that evaluates the performance and application of an omniphobic surface treatment for solids deposition prevention in industrially relevant systems. In particular, this specific work concentrates on the efficacy and robustness of the treatment under fully flowing conditions. The apparatuses utilized for this include two flowloops: a lab-scale, high-pressure flowloop for gas hydrate and surface treatment durability studies, and a bench-scale, atmospheric pressure loop for crude oil and asphaltene experiments. Film growth in high-pressure flowloop tests corroborated previous reports of delayed gas hydrate nucleation observed in rocking cells. Without the aid of the memory effect, treated oil-dominated experiments never experienced hydrate formation, spending upwards of a week in the hydrate stability zone (at the subcooled/fluid test conditions). Subsequent tests which utilized the memory effect then revealed that the hydrate formation rate reduced in the presence of the surface treatment compared to a bare stainless-steel surface. This testing was part of a larger set of trials conducted in the flowloop, which lasted about one year. The surface treatment durability under flowing conditions was evaluated during this time. Even after experiencing ∼4000 operating hours and 2 full pressure cycles, no evidence of delamination or damage was detected. Finally, as part of an extension to previous work, corroded surface asphaltene deposition experiments were performed in a bench-top flowloop. Treated experiments displayed an order of magnitude reduction in both total oil (all fractions of crude oil) and asphaltene fraction deposited.

03 NATURAL GAS↗
DOE OSTI2020

A spatial-statistical investigation of surface expressions associated with cyclic steaming in the Midway-Sunset Oil Field, California

In the Midway Sunset Oil Field in Central California, operators inject steam into the shallow diatomite formation to enhance heavy oil recovery through imbibition, wettability alteration, and viscosity reduction, among other mechanisms. The injected steam, however, does not always remain in the reservoir or return through the wells. In two zones in the study area, the steam comes out at the surface, creating sinkholes, seeps, and steam outlets. These phenomena, called “surface expressions,” pose safety and environmental hazards. Even though these surface expressions are a widespread problem in Central California, they are not well documented and understood. Possible causes of the surface expressions include: high injection pressure, structurally controlled flow patterns, leakage of steam through old improperly abandoned wells, high injection volumes, or flow along naturally occurring faults, among other possible factors. This work examines attributes of the zones with surface expressions in order to determine factors that may contribute to their occurrence. Spatial statistical analysis using logistic regression, random forests, and classification trees is used to explore the relationship between the surface expressions and geological and production-related attributes. The results point to a significant spatial correlation between the surface expressions and two predictors: concentration of plugged wells and geologic seal thickness. The results guide follow-up studies to further investigate the role of well abandonment and seal thickness in the occurrence of surface expressions.

02 PETROLEUM↗
DOE OSTI2025

An Assessment of the Error Due to Computing Waste Isolation Pilot Plant Porosity Using the Porosity Response Surface Approach

The Waste Isolation Pilot Plant Performance Assessment (WIPP PA) must predict the likelihood that radionuclides will escape into the biosphere via mechanisms that depend on geohydraulic flow. Ideally, one would predict the geohydraulic flow using coupled geohydraulic and geomechanical simulations, but such coupled simulations are not computationally tractable. Instead, Sandia has historically used a look-up table of porosities for a given fluid pressure and time, called the porosity response surface, but this approach can introduce porosity errors because it largely ignores the porosity’s dependence on the past fluid pressure history. This report discusses efforts to quantify these porosity errors for both the legacy and new porosity response surfaces. Six hundred different fluid pressure histories were fed through the legacy/new geomechanical model and the legacy/new porosity response surface to generate six hundred porosity error histories. The error associated with the legacy porosity surface was substantial, while the error associated with the new porosity surface was typically small, except when fluid pressures exceeded the lithostatic pressure at the repository. In response to the errors at high pressures, a preliminary study of the WIPP PA’s sensitivity to these porosity errors was conducted. The study found that reducing the porosity errors at high pressures negligibly affected predictions of radionuclide releases. Finally, an initial machine-learned model for porosity was developed. This ML model significantly reduced the porosity error at high pressures, but sizable errors remained, so more development is necessary before coupling an ML model to the geohydraulic model.

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

Atomic-Scale Investigation of Oxidation at the Black Phosphorus Surface

Black phosphorus (BP) exhibits extraordinary electronic properties that are desirable for a wide variety of electronic and optoelectronic applications. However, applications of BP are hindered by its rapid degradation in ambient conditions. Despite significant advances that have been made in understanding the degradation mechanism, no consensus has yet been reached on how BP oxidation occurs at the atomic scale as experimental studies have been mostly restricted to averaged effects of degradation over a micron- to millimeter-sized region. Here, BP oxidation is investigated using scanning tunneling microscopy/spectroscopy (STM/S). Introducing O 2 gas to the BP surface in ultrahigh vacuum at a pressure of 10 –5 mbar for 1 min creates two new types of defects on the surface. We identify these defects as dangling atomic oxygen and phosphorus multivacancies using density functional theory simulations. In addition to the structural changes to the surface, the electronic structure is also drastically altered by the introduction of oxygen. The 300 meV band gap of BP is lifted due to dosing. This change in the electronic structure is reversible through STM tip manipulation. As a result, these are the first experimental results showing the atomic-scale oxidation of BP, an important step toward understanding the degradation process

2D materials↗