Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Pore size distribution”

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.

360 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Clear as mud redefined: Tunable transparent mineral scaffolds for visualizing microbial processes below ground

Microbes inhabiting complex porous microenvironments in sediments and aquifers catalyze reactions that are critical to global biogeochemical cycles and ecosystem health. However, the opacity and complexity of porous sediment and rock matrices have considerably hindered the study of microbial processes occurring within these habitats. Here, we generated microbially compatible, optically transparent mineral scaffolds to visualize and investigate microbial colonization and activities occurring in these environments, in laboratory settings and in situ. Using inexpensive synthetic cryolite mineral, we produced optically transparent scaffolds mimicking the complex 3D structure of sediments and rocks by adapting a suspension-based, freeze-casting technique commonly used in materials science. Fine-tuning of parameters, such as freezing rate and choice of solvent, provided full control of pore size and architecture. The combined effects of scaffold porosity and structure on the movement of microbe-sized particles, tested using velocity tracking of fluorescent beads, showed diverse yet reproducible behaviors. The scaffolds we produced are compatible with epifluorescence microscopy, allowing the fluorescence-based identification of colonizing microbes by DNA-based staining and fluorescence in situ hybridization (FISH) to depths of 100 µm. Additionally, Raman spectroscopy analysis indicates minimal background signal in regions used for measuring deuterium and 13 C enrichment in microorganisms, highlighting the potential to directly couple D 2 O or 13 C stable isotope probing and Raman-FISH for quantifying microbial activity at the single-cell level. To demonstrate the relevance of cryolite scaffolds for environmental field studies, we visualized their colonization by diverse microorganisms within rhizosphere sediments of a coastal seagrass plant using epifluorescence microscopy. The tool presented here enables highly resolved, spatially explicit, and multimodal investigations into the distribution, activities, and interactions of underground microbes typically obscured within opaque geological materials until now.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Sequentially coupled flow and geomechanical simulation with a discrete fracture model for analyzing fracturing fluid recovery and distribution in fractured ultra-low permeability gas reservoirs

More accurate characterization and prediction of the in-situ distribution of fracturing fluid in fractured reservoirs are needed for enhancing well productivity. In this study, an implicit-sequentially coupled flow/geomechanics simulator incorporating an efficient discrete fracture model is developed to model fluid distribution and recovery performance of ultra-low permeability gas reservoirs. The finite-volume and finite-element methods are used for space discretization of the flow and geomechanics equations, respectively, while the backward Euler method is employed for time discretization. The flow and geomechanics equations are solved sequentially based on fixed-stress splitting. An efficient discrete-fracture model is used to explicitly model the fractured system. Flexible unstructured gridding is employed to model arbitrarily-oriented fractures. The interrelations among pore volume, permeability and geomechanical conditions are considered dynamically using two-way coupled flow and geomechanics computations. The geometry of fracture (networks) due to hydraulic fracturing has significant impacts on the fracturing fluid recovery efficiency and ensuing fluid distribution. Under the same injection volume, the fracturing fluid recovery is higher when the fracture geometry is planar. Fluid recovery is relatively lower whenever natural fractures are activated during fracturing treatments; flowback time is also shortened when complex fracture network with enlarged fracture interface is present. Fracturing fluid in hydraulic fractures may leak off into the natural fractures and subsequently imbibes into the surrounding matrix due to capillarity effects. The fracturing fluid recovery and in-situ fluid distribution are sensitive to the shut-in duration and fracture closure behavior. This study analyzes the coupled flow-geomechanical responses of fractured gas reservoirs during the post-fracturing periods. Understanding the fate of the fracturing fluid can provide insights on, to some extent, the stimulated fracture volume, size of the water invasion zone, and efficiency of the fracturing design. The simulation predictions can also provide more accurate initial reservoir conditions (e.g. distributions of different phases and pressure) for long-term well performance estimation.

42 ENGINEERING↗
DOE OSTI2022

Fatigue performance of laser powder bed fusion hydride-dehydride Ti-6Al-4V powder

We report hydride-dehydride (HDH) Ti-6Al-4V alloy with particle size distribution of 50–120 µm is laser powder bed fusion (L-PBF) processed using optimum processing parameters and a near-fully dense structure with a density of 99.9 % is achieved. Microstructural observations and phase analyses indicate formation of columnar β grains with acicular α/α' phases in as-built condition. The roughness of the as-fabricated samples is significant with an average roughness of R a = 15.71 ± 3.96 µm and a root mean square roughness of R rms = 108.4 ± 24.9 µm, however, both values are reduced to R a = 0.19 ± 0.04 µm and R rms = 4.9 ± 0.6 µm after mechanical grinding. Mechanical tests are carried out on as-fabricated specimens followed by stress relief treatment. All samples are tested to failure in fatigue, under fully-reversed tension-compression conditions of R = –1. The as-built samples failed from the surface with crack initiation mainly at micro-notches, whereas after mechanically grinding, crack initiation changed to subsurface defects such as pores. Minimizing surface roughness by mechanically grinding eliminates surface micro-notches which improves fatigue strength in the high cycle fatigue region. Fatigue notch factor calculations showed that the effect of surface roughness was significantly lower when HDH powder is used compared to standard spherical powder. X-ray diffraction analysis revealed an in-plane compressive stress, micro-strain and grain refinement on the surface of the mechanically ground samples. Fractography observations (macroscale) revealed a fully brittle fracture in the first stage of crack growth with a transition to a dominantly ductile fracture in the third stage of crack growth. On the other hand, at the micro scale, even the brittle fracture regions showed evidence of ductile fracture within the α' martensite laths.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Understanding Platinum Ionomer Interface Properties of Polymer Electrolyte Fuel Cells

A well-designed cathode catalyst layer with optimal ionomer distribution is critical to minimizing amount of Platinum (Pt) content in polymer electrolyte fuel cells (PEFCs). The impact of Pt loading, ionomer content and carbon support types on the catalyst/ionomer interface were also investigated at dry and wet conditions. Higher Pt loadings resulted in higher double layer capacity (C dl ) and similar electrochemical surface area (ECSA) due to well dispersed ionic phase material. Higher ionomer content resulted in higher ionic conductivity but also showed similar SO 3 − group coverage. High surface area (HSA) carbon support had larger ECSA and C dl at both dry and wet conditions, as less agglomerated Pt was well dispersed in the meso-pores of the support. Lower SO 3 − group coverages were observed for HSA carbon than for Vulcan carbon due to Pt particles being buried within the porous HSA carbon support. The effect of cell conditioning and voltage recovery on the PEFC cathode catalyst layer was shown to have minimal impact on SO 3 − group coverage despite a decrease in C dl and ECSA due to the size increase of Pt particles. At dry condition, a significant increases in SO 3 − group coverage were observed for all MEAs due to higher adsorptivity of ionomer in dry conditions.

25 ENERGY STORAGE↗
NASA NTRS2019

Cleanroom Contamination Identification Method Development

During fabrication, assembly, and testing of spacecraft and flight hardware it is vital to avoid contaminants that can cause degradation and could result in significant failure. Yet, there is no existing contamination monitoring method that provides the identity of airborne particles in a cleanroom facility. Knowing the particle identities, would allow scientists and engineers to determine the source of the contaminants and prevent setbacks before they occur or cause damage. Current cleanliness monitoring methods include airborne particle counters (APCs), fallout filters, and visual inspections. Particle counts from APCs are the primary metric used to define a cleanroom class and hence its level of cleanliness, but do not provide identification nor can they differentiate between large and small sizes of particles. In addition, using fallout filters is not a proactive, timely, or representative approach to cleanroom contamination monitoring because these samples are only retrieved after 30 days and are placed away from spacecraft processing to avoid interference with operations. In contrast, the forced air sampling method can collect a sample within an hour at any location required and provide results in less than a day. This system uses a cassette and filter sample medium to capture airborne particles which are then taken to a scanning electron microscope with energy dispersive spectroscopy (SEM/EDS) to identify and size the captured particles. Development of forced air sampling into an established laboratory capability will allow for fast sampling and routine identification of unknown contamination sources within the cleanroom. The test method development required market research for an air sampling cassette that increases sample collection efficiency and a filter with low enough background contamination to allow differentiation between a blank (control) and the collected sample. It was determined that a conductive black cassette and a polycarbonate filter were the best options. Conductive black cassettes, in comparison to the standard styrene, are manufactured using polypropylene filled with carbon. This makes the cassette conductive and minimizes the tendency of particles to stick to the wall of the cassette due to electrostatic force. In previous trials a mixed cellulose ester (MCE) filter was used to capture the contaminants, however the rougher surface of the filter contributed to entrapment of the particles within the filter structure and made it harder to identify the particles. In comparison, track etched polycarbonate filters have random cylindrical pores and a smooth surface which contributes to uniform sample distribution on the surface of the filter. Future work includes: testing the system using control samples to determine the efficiency and suitability of the medium, performing sample collection in various environments to establish ideal operating parameters and analyzing contaminant particles using SEM/EDS and assistant characterization techniques. Once fully developed, employing the forced air sampling method will help to prevent damage to spacecraft, avoid schedule delays, and allow for mission success.

Hernandez Melendez, Jailyn M.↗
DOE OSTI2022

Mineralogical, magnetic and geochemical data constrain the pathways and extent of weathering of mineralized sedimentary rocks

The oxidative weathering of sulfidic rock can profoundly impact watersheds through the resulting export of acidity and metals. Weathering leaves a record of mineral transformation, particularly involving minor redox-sensitive phases, that can inform the development of conceptual and quantitative models. In sulfidic sedimentary rocks, however, variations in depositional history, diagenesis and mineralization can change or overprint the distributions of these trace minerals, complicating the interpretation of weathering signatures. Here we show that a combination of bulk mineralogical and geochemical techniques, micrometer-resolution X-ray fluorescence microprobe analysis and rock magnetic measurements, applied to drill core samples and single weathered fractures, can provide data that enable the development of a geochemically consistent weathering model. This work focused on one watershed in the Upper Colorado River Basin sitting within the Mesaverde Formation, a sedimentary sandstone bedrock with disseminated sulfide minerals, including pyrite and sphalerite, that were introduced during diagenesis and subsequent magmatic-hydrothermal mineralization. Combined analytical methods revealed the pathways of iron (Fe), carbonate and silicate mineral weathering and showed how pH controls element retention or release from the actively weathering fractured sandstone. Drill core logging, whole rock X-ray diffraction, and geochemical measurements document the progression from unweathered rock at depth to weathered rock at the surface. X-ray microprobe analyses of a 1-cm size weathering profile along a fracture surface are consistent with the mobilization of Fe(II) and Fe(III) into acidic pore water from the dissolution of primary pyrite, Fe-sphalerite, chlorite, and minor siderite and pyrrhotite. These reactions are followed by the precipitation of secondary minerals such as of goethite and jarosite, a Fe-(oxyhydr)oxide and hydrous Fe(III) sulfate, respectively. Microscale analyses also helped explain the weathering reactions responsible for the mineralogical transformations observed in the top and most weathered section of the drill core. For example, dissolution of feldspar and chlorite neutralizes the acidity generated by Fe and sulfide mineral oxidation, oversaturating the solution in both Fe-oxides. The combination of X-ray spectromicroscopy and magnetic measurements show that the Fe(III) product is goethite, mainly present either as a coating on fracture surfaces in the actively weathering region of the core or more homogeneously contained within the unconsolidated regolith at the top of the core. Low-temperature magnetic data reveal the presence of ferromagnetic Fe-sulfide pyrrhotite that, although it occurs at trace concentrations, could provide a qualitative proxy for unweathered sulfide minerals because the loss of pyrrhotite is associated with the onset of oxidative weathering. Pyrrhotite loss and goethite formation are detectable through room-temperature magnetic coercivity changes, suggesting that rock magnetic measurements can determine weathering intensity in rock samples at many scales. In conclusion, this work contributes evidence that the weathering of sulfidic sedimentary rocks follows a geochemical pattern in which the abundance of sulfide minerals controls the generation of acidity and dissolved elements, and the pH-dependent mobility of these elements controls their export to the ground- and surface-water.

58 GEOSCIENCES↗
DOE OSTI2023

Early self-organization of fission gas bubble superlattice formation in neutron-irradiated monolithic U-10Mo fuels

Self-organization of defect superlattices in far-from-equilibrium systems presents a promising way to mitigate swelling concerns in nuclear materials. The gas bubble superlattice (GBS) is a highly ordered, three-dimensional complex defect structure that can retain fission gasses in Uranium-Molybdenum (U-Mo) fuels. Transmission electron microscopy (TEM) investigation of monolithic U-10Mo fuel irradiated to 1.15 × 10 21 fissions/cm 3 and 1.30 × 10 21 fissions/cm 3 revealed that early-stage ordering preferentially occurs at the grain boundaries (GB) and that the critical bubble size for complete ordering is ~3 nm. Once formed at the GB, the GBS extends towards the grain interior; however, the spread in distance from the GB varies likely depending on the type and strength of the GB sink. TEM results also showed a possible correlation between the growth and evolution of the intragranular disordered bubbles and large dislocation networks. The fission product distribution in and outside of the GBS was also investigated confirming the presence of xenon in the GBS, as well as other fission products including cesium, barium, lanthanum, and cerium.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

A New Model for Simulating the Imbibition of a Wetting-Phase Fluid in a Matrix-Fracture Dual Connectivity System

The imbibition experiment is an effective approach for measuring petrophysical properties of porous media, with many such experiments performed over the past decade. Quite some empirical, analytical, and numerical models have been developed to simulate spontaneous imbibition of the wetting phase fluid into porous media, but limitations still exist. In previous studies, the imbibition process has been considered to give a piston-like displacement or the porous medium modeled as multiply-sized pores linked with bonds; both approaches fail to yield comprehensive results due to their neglect of the presence of irregular fractures or nonuniform flow paths through the matrix. By building a numerical model for simulating laboratory-scale experimental data, we performed imbibition tests on several fractured Barnett Shale samples having fractures either parallel ( P ) or transverse ( T ) to the bedding plane and used MATLAB to build a new numerical model by combining the imbibition process in fractures and the matrix using concepts from percolation theory. The experimental data show that the rocks with P -direction fractures have a more steady increase of imbibition rates than the case of T -direction one. As the shale matrix with low pore connectivity hampers the upward water movement, the imbibition rate of shales with T -direction fractures will decrease suddenly after the bottom layer in contact with water is saturated during the initial period. This wetting phase movement (WPM) model can simulate 3D porous media with 2D fractures. The rate of imbibition by fractured porous media is associated with physical parameters such as porosity and fracture distribution (e.g., the number and angle of fractures). Using Monte Carlo methods, we examined fracture parameters and predicted elapsed time and cumulative water imbibition, for the Barnett Shale samples. The results show that the rate of imbibed water mass is sensitive to the number of fractures directly connected to water source, and the connectivity between two neighboring grid cells is a key parameter for the wetting-front progression. The findings of this study can help to better understand the imbibition process with multiple influencing processes and factors in fractured-matrix rocks. Although the experiments, data simulation, and prediction results are based only on Barnett Shale samples, the model is readily applicable to imbibition tests of other fractured rocks to show the spatial and temporal behavior during a dynamic imbibition process that are not easily captured experimentally.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Quantifying Properties for a Mechanistic, Predictive Understanding of Aqueous Impact on Aging of Medium and Low Voltage AC and DC Cabling in Nuclear Power Plants

To address the gap in knowledge in understanding degradation in relevant service conditions, this project aimed to develop a mechanistic, predictive model of medium and low voltage cable failure based on the primary environmental degradation parameters of aqueous immersion time, temperature, and the oxidation extent. To do so, we took a two-fold approach toward evaluating degradation as related to the aqueous condition. First, we evaluated the chemical, mechanical, and electrical properties of polymers that comprise the cabling insulation under varied aqueous conditions at different temperatures. Secondly, we evaluated the performance of insulation materials under similar accelerated aging conditions in addition to the dielectric breakdown of the cabling under submersion conditions. Thermal oxidation and immersion of LDPE thin films were done using a Parr vessel and the extent of oxidation was monitored using the carbonyl index from ATR-FTIR spectrum, a ratio quantifying the carbonyl functionality produced on the LDPE film surface. Increasing the temperature, oxidation time, and oxygen pressure increased the water-vapor permeability and the carbonyl content. Thermogravimetric analysis of the films confirmed that there is indeed an initial weight loss of the smaller molecular weight molecules and volatiles. At higher temperatures, there is a secondary mass loss. Up to 70°C, the mass loss curves were similar until 80°C, when the aged films decreased at a much greater rate with respect to increasing temperature. This was seen in the 10, 50, 90% weight loss temperatures. The mechanical properties of sheet PP and HDPE as well as injection molded LDPE, HDPE, HDPP dog bones were plotted with respect to the time spent in the Parr vessel during thermal oxidation. For the sheet PP and HDPE, the UTS and the modulus of elasticity decreased while the elongation increased. The injection-molded LDPE, HDPE, and HDPP dog bones experienced a similar trend where the changes in the properties were within experimental uncertainty. This analysis of LDPE in dry and immersive oxygen-rich environments is an important baseline for future experiments looking at other degradation mechanisms. A predictive aging model for low-density polyethylene insulative cable housings was developed. This model revealed a combination of physical and chemical processes which lead to an increase in permeability and a decrease in strength of the insulator. In creating this predictive model, a gap of knowledge in the literature was revealed in two areas: understanding how the crystallinity of a polymer changes with time and visualizing the predicted pores formed through the insulator which leads to failure. Developing and using an adapted ATR-FTIR method, crystallinity was monitored and compared to bulk crystallinity found via DSC. Inhomogeneity in crystallinity changes were observed but the FTIR method but was found to be less reliable. Pore visualization was found utilizing electrical impedance spectroscopy saturated aged polyethylene films with synthesized citrate-capped gold nanoparticles, and chloroauric acid precursor. results indicate a decreasing impedance of the polyethylene films caused by the transport of ions through the film. SEM imaging was then performed for elemental analysis of the films and counter electrodes for the presence of ions. Chloride and gold ions were detected; however, no nanoparticles were found. This gives an estimated pore size of at least 0.3 nm through the aged film. Cyclic submergence of polyethylene and polypropylene in aqueous solutions of copper sulfate and Harrison’s solution were examined. It was discovered that aging in these mixed conditions and at 90°C did cause a small, yet significant increase to tensile strength for both unaged polyethylene and polypropylene. The PE and PP in cycled and submerged conditions did not have tensile strengths significantly different from dry-aged after 16 weeks. For both solutions, it was determined that capacitance increases with both water tree depth and cable temperature and is relatively independent of changing water tree AR. As for resistance, there is no apparent change between depth percentages of 10 and 80 for both solutions. This is because the water tree has not yet entered the conductor and thus there is no shorting current flowing through the tree yet. Between 80 and 100% water tree depth, it is observed that for both solutions there is a rapid decrease in resistance because the tree is now impacting the conductor shield and allowing shorting current to flow through it. The relationship between resistance and temperature was found to be inversed. As temperature increases, resistance decreases for both water tree solutions with distilled water having the most apparent difference across the range of temperatures simulated. With regards to water tree AR, resistance was found to be relatively independent at depth percentages of 10, 50, and 70, but there was some relationship at 100% depth for both solutions. At this depth, as water tree AR increased, resistance was found to also increase. Through the examination of voltage and electric field distribution plots, it was confirmed for EPR cables that increasing water tree depth leads to increasing distortion. The increase in shorting current at 100% depth for aqueous copper sulfate compared to distilled water was also visualized. Lastly, the relationship between shorting current and temperature for narrow water trees was visualized and validated using the distribution plots. This study has led to a better understanding of the effect of cable temperature, water tree solution, and geometry of the tree region on cable degradation. Furthermore, this study has documented simulation results of water tree degradation in EPR MV cables, which had been previously lacking in information. Future work in this area will add frequency into the mix and examine what effect it has on the rate of cable degradation as a result of water treeing. To compliment the work being carried by UMD on different polymer chemistries, harvested medium voltage (MV) cables were selected for accelerated aging study at ORNL. These medium voltage cables are representative of the vintages and materials that are currently in use at existing nuclear power plants (NPPs). These cables were immersed in water at 90°C and energized for a period of two years at elevated voltage. Partial discharge was measured during this period of time to track potential degradation. Unfortunately, due to the absence of failure, the extent of integration between the UMD techniques and the harvested MV insulation was limited. However, based on the findings and from the UMD research in the previous sections on polyethylene and polypropylene, follow-on characterization for harvested MV cable insulation should focus on sample preparation to take advantage of the UMD techniques to better understand mechanistic degradation in MV cable insulations. This would include: 1.) Utilization of solutions and impedance measurement insulation to track permeability in harvested and accelerated aging insulation samples, 2.) Adaptation of gold nanoparticle synthesis for electrical impedance spectroscopy with supporting SEM characterization to study pore formation in harvested and accelerated aged insulation samples, & 3.) ATR-FTIR crystallinity with multiple temperature dependent harvested MV cable insulation under different air and submerged accelerated aging. In addition, low voltage dielectric spectroscopy and high voltage dissipation factor, tan δ, could be implemented as a condition monitoring tool for harvested MV cable insulation in submerged environments.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Improved Understanding of Zircaloy Hydrogen Pickup Mechanisms in BWRs. NSUF Report (Final Report)

Using an instrumented autoclave operating at 288°C and 1100 psi, electrical impedance spectroscopy (EIS) measurements were taken from three neutron irradiated fuel channel samples (Zircaloy-2 and Zircaloy-4). The three high priority fuel channel samples were tested at 288°C using a simulated PWR water with the AC frequency range extending down through 10 -4 Hz, with corresponding DC measurements to validate the resistance measurements. EIS measurements were also taken at room temperature using a simulated PWR water from four neutron irradiated fuel channel samples (Zircaloy-2 and Zircaloy-4). and four water rod samples (Zircaloy-2). EIS measurements were also performed on unirradiated Zircaloy-2 tubes that were previously oxidized in autoclaves at various temperatures for different time durations to obtain oxide layers of different thicknesses. These experiments proved to be more challenging than anticipated because of the small size of the samples, which required development of fixtures to allow valid testing of the materials, some of which are radioactive, under pressure and temperature. Equivalent circuits were used to model the data. Some samples resulted in Nyquist plots where the arcs were too close to each other resulting in the inability to separate the different oxide morphological thicknesses. Microstructural characterization in the form of optical microscopy, scanning electron microscopy and (scanning) transmission electron microscopy was performed to document the microstructural features of the oxidized Zr-2 and Zr-4 samples (oxide and metal/oxide interface) and to interpret the EIS results. At the start of this work, it was assumed that despite the laminar cracks and porosity visible in both optical and SEM, the oxide scale is dense near the metal, effectively preventing the water from contacting the metal directly. While it is generally accepted that the large laminar cracks pores allow the water to percolate into the oxide film, it is unknown just how deep it can go, and how interconnected is the porosity, which will dictate how much the porosity is connected electrically to the water within the oxide scale. The EIS results demonstrated that there is in fact a distribution of oxide thicknesses over the various samples, and within each material condition, that are communicating to the liquid. This is derived from the fact that Nyquist plots have an arc radius whose center point is often below the Z’ axis, implying that instead of a single oxide thickness, these oxides have a distribution of thicknesses. This distribution of thicknesses produces arcs that are not clearly separable in most of the samples.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Quantifying CO 2 Plume Stabilization at Carbon Storage Projects, North Dakota, USA

This study presents an approach for quantifying when injected carbon dioxide (CO 2 ) stabilizes pursuant to carbon capture and storage (CCS) project permitting and site closure requirements. The distribution of mobile-phase CO 2 (CO 2 plume) will evolve within the storage reservoir during and after injection through both physical and chemical trapping mechanisms. CCS policies generally agree that the CO 2 plume’s migratory behavior in post-injection should demonstrate nonendangerment to the environment but do not provide specific guidance on how to meet the definition of plume stabilization, generating some uncertainty for operators. Plume stability herein means the CO 2 plume 1) changes size minimally and predictably in the storage reservoir such that it will not cross key boundaries identified in the permit and 2) does not pose a threat to human health, underground sources of drinking water (USDWs), and the environment because of lateral migration to areas where leakage pathways may exist. Published literature on plume metrics was reviewed to determine which metric(s) may be most appropriate for determining CO 2 plume stability. A technical approach that defines plume stabilization by estimating the rate of change in the geographic footprint of the CO 2 plume with respect to time was developed and illustrated using a case study from North Dakota, USA, as a proposed solution for CCS operators to apply at the project permitting stage. Any prospective CCS operator may benefit from using the same approach to inform the selection of pore space lease and monitoring areas and develop post-injection site care plans.

03 NATURAL GAS↗
NASA NTRS2001

Application of Computer Tomography for Life Detection

Perhaps one of the most fundamentally difficult challenges facing those who would search for life is that of scale determination. Spatial scales of life on Earth range over more than 15 orders of magnitude in mass and volume, and more than 8 orders of magnitude in 2 dimensional space. If the distribution of life is sparse in comparison to the background on which it is found, then the choice of the right scale is critical to finding that life. But how does one identify the proper scale? To put this in other words, how does one recognize the "haystacks" in which the needles (biosignatures and evidence of life) might be most profitably searched for? The problem is further exacerbated when conditions get extreme because much of the life moves from the clement surface environment into the pores and more clement environments inside of rocks, minerals and soils. Once encased in their lithic homes, these microbes become nearly impossible to study by standard techniques because of the opacity of the rocks. It is this problem that we propose to address in the work proposed here. Computer Tomography (CT) has been a very valuable tool in medicine, where the best resolution available has typically been of the order of about 0.5 mm. However, to adapt the approach for life detection of microbial endoliths, the resolution needs to be moved to the micrometer and even submicrometer levels. Thus for the studies proposed here, we begin with a commercially available instrument that can yield resolution of approximately 10 micrometers. The rational for this is twofold: first, this is the "state of the art" in laboratory instruments; and second, that while the usual size of a microbial cell is about 1 micron, microorganisms tend to live in communities that usually exceed the 10 micrometer size range. The resolution also depends on the sample size itself, so having a small lab instrument into which small samples can be placed will be beneficial to the resolution. We have now used several different CT systems, beginning with the medical scanners (Arcadia CT group) for the detection of layered communities in sandstone rocks from Antarctica. Even this crude instrument was able to point to the areas of the rock that were dominated by microbial populations - this provides the critical first information that says, "Go back and look at these sites with other methods." We showed that without sample preparation or destruction it was possible to gain knowledge as to the presence of density differences suggestive of life.

Tsapin, A.↗
DOE OSTI2021

Influence of Al location on formation of silver clusters in mordenite

Formation of zeolite supported Ag 0 clusters depends on a combination of thermodynamically stable atomic configurations, charge balance considerations, and mobility of species on the surface and within pores. Periodic density functional theory (DFT) calculations were performed to evaluate how the location of Al in the mordenite (MOR) framework and humidity control Ag 0 nanocluster formation. Four Al framework sites were studied (T1-T4) and the Al positions in the framework were identified by the shifts in the differential Al…Al pair distribution function (PDF). Furthermore, structural information about the Ag 0 nanoclusters, such as dangling bonds, can be identified by Ag…Ag PDF data. For Ag 0 formation in vacuum MOR structures with a Si:Al ratio of 5:1 with Al in the T1 position resulted in the most framework flexibility and the lowest Ag 0 nanocluster charge, indicating the best result for formation of charge neutral nanoclusters. When water is present, Al in the T3 and T4 positions results in the formation of the smallest average Ag 0 nanoclusters plus greater expansion of the O-T-O bond angle than in vacuum, indicating easier diffusion of the Ag 0 nanoclusters to the surface. Here, the presence of Al in 4-membered rings and in pairs indicates favorable MOR structures for formation of single Ag atoms, despite the existence of synthesis challenges. Therefore, Al in the T2 position is the least favorable for Ag 0 nanocluster formation in both vacuum and in the presence of water. Al in the T1, T3, and T4 positions provides beneficial effects through framework flexibility and changes in nanocluster size or charge that can be leveraged for design of zeolites for formation of metallic nanoclusters.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Microstructure‐Dependent Sodium Storage Mechanisms in Hard Carbon Anodes

Sustainable energy storage is essential to support the transition to renewables and meet the increasing demand for energy. Sodium‐ion batteries (NIBs) are attractive for grid‐scale energy storage due to the abundance and low cost of sodium, sustainability of other battery components, and electrochemical performance. Hard carbon (HC) is a leading anode material for NIBs, but its complex microstructure complicates the understanding of sodium storage mechanisms. Using X‐ray total scattering and density functional theory calculations, this study clarifies how HC's microstructural variations influence sodium storage across the slope (high potential) and plateau (low potential) regions of the potential capacity curve. In the slope region, sodium initially adsorbs at high‐binding energy defect sites and subsequently intercalates between graphene layers, adsorbing at low‐binding energy defect sites, correlating with different slopes observed during initial sodiation. Initial irreversibility arises from sodium trapping at surface defects and solid electrolyte interface formation. In the plateau region, sodium simultaneously intercalates and fills pores, influenced by pore size, interlayer spacing, and defect concentration. HCs with larger pore sizes form larger sodium clusters. In conclusion, the proposed mechanism underscores the role of microstructure engineering in enhancing HC performance and advancing NIBs for grid‐scale energy storage.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Investigation of the Effect of Injected CO 2 on the Morrow B Sandstone through Laboratory Batch Reaction Experiments: Implications for CO 2 Sequestration in the Farnsworth Unit, Northern Texas, USA

About one million tons of CO 2 have been injected into the Farnsworth unit to date. The target reservoir for CO 2 injection is the Morrow B Sandstone, which is primarily made of quartz with lesser amounts of albite, calcite, chlorite, and clay minerals. The impact of CO 2 injection on the mineralogy, porosity, and pore water composition of the Morrow B Sandstone is a major concern. Although numerical modeling studies suggest that porosity changes will be minimal, significant alterations to mineralogy and pore water composition are expected. Given the implications for CO 2 storage effectiveness and risk assessment, it is crucial to verify the accuracy of theoretical model predictions through laboratory experiments. To this end, batch reaction experiments were conducted to model conditions near an injection well in the Morrow B Sandstone and at locations further away, where the CO 2 has been diluted by formation water. The laboratory experiments involved submerging thin sections of both coarse- and fine-grained facies of the Morrow B Sandstone in formation water samples with varying levels of CO 2 . The experiments were conducted at the reservoir temperature of 75 °C. Two experimental runs were conducted, one lasting for 61 days and the other for 72 days. The initial fluid composition used in the second run was the same as in the first. The mineralogy changes in the thin sections were analyzed using SEM and the Tescan Integrated Mineral Analyzer (TIMA), while changes in the composition of the formation water were determined using ICP-AES. During each experiment, a thin layer of white fine-grained particles consisting mainly of dolomite and silica formed on the surface of the thin sections, leading to significant reductions in Ca, Mg, and Sr in the formation water. This outcome is consistent with numerical model predictions that dolomite would be the primary mineral that would react with injected CO 2 and that silica would be oversaturated in the formation water. Changes in mineral abundance in the thin sections themselves were much less systematic than in the theoretical modeling experiments, perhaps reflecting heterogeneities in the mineral grain size surface area to volume ratios and mineral distributions in the thin sections not considered in the numerical models.

58 GEOSCIENCES↗
DOE OSTI2024

Impact of Pressure-Dependent Interfacial Tension and Contact Angle on Capillary Heterogeneity Trapping of CO2 in Storage Aquifers

Summary Carbon dioxide (CO2) capillary trapping increases the total amount of CO2 that can be effectively immobilized in storage aquifers. This trapping, manifesting itself as accumulated CO2 columns at a continuum scale, is because of capillary threshold effects that occur below low-permeability barriers. Considering that capillary pressure is dictated by heterogeneous pore throat size, the trapped CO2 column height and associated CO2 saturation will vary spatially within a storage aquifer. This variation will be influenced by two pressure-dependent interfacial parameters—CO2/brine interfacial tension (IFT) and CO2/brine/rock contact angle. Our objective is to understand how the pressure dependence of these two parameters affects the heterogeneity of capillary trapped CO2 at a continuum scale. Our conceptual model is a 1D two-zone system with the upper zone being a flow barrier (low permeability) and the lower zone being a flow path (high permeability). The inputs to this model include microfacies-dependent capillary pressure vs. saturation curves and permeability values. The input capillary pressure curves were collected in the literature that represents carbonate microfacies (e.g., dolograinstone) in a prevalent formation in the Permian Basin. We then used the Leverett j-function to scale the capillary pressure curve for the two zones that are assigned with the same or different microfacies. During scaling, we considered the influence of pressure on both the IFT and contact angle of CO2/brine/dolomite systems. We varied the zone permeability contrast ratio from 2 to 50. We then assumed capillary gravity equilibriums and calculated the CO2 saturation buildup corresponding to various trapped CO2 column heights. The CO2 saturation buildup is defined as the CO2 saturation in the lower layer minus that in the upper one. We found that the saturation buildup can be doubled when varying pressure in a storage aquifer, after considering pressure-dependent IFT and contact angles. Thus, assuming these two parameters to be constant across such aquifers would cause large errors in the quantification of capillary trapping of CO2. The whole study demonstrates the importance of considering pressure-dependent interfacial properties in predicting the vertical distribution of capillary trapped CO2. It has important implications in developing a better understanding of leakage risks and consequent storage safety.

Engineering↗
DOE OSTI2025

Statistical evaluation of microscale stress conditions leading to void nucleation in the weak shock regime

Here, we investigate the heterogeneity of the stress state driven by anisotropic deformation response at the single crystal level through five statistical volume element (SVE) calculations of polycrystalline BCC tantalum. This work focuses on grain boundaries as a prominent material defect type prone to void nucleation based upon experimental observations of predominantly intergranular void nucleation in this material. The SVEs are constructed to be statistically representative of larger volumes of material and are meshed such that mean and standard deviation of grain size and orientation information is reconstructed. The computational meshes feature hexahedral (brick) elements and smooth conformal grain boundaries where significant stress concentration is known to occur, a tail effect of interest in the extreme events process of dynamic ductile damage. An existing micromechanical crystallographic plasticity model shown to capture the single crystal behavior of BCC tantalum well is used to perform the polycrystal calculations. The model includes representation of the non-Schmid effect of non-planar screw dislocation kinetics in tantalum. A three-dimensional stress state time profile predicted by damage modeling of a flyer plate impact experiment is applied as boundary conditions to each SVE. Resulting grain boundary stress state statistics are strongly non-Gaussian. Significant structural evolution is observed within the compressive hold before unloading into tension in the stress profile. Strong angular dependence of grain boundary traction magnitude with shock direction is observed. Non-Schmid effects continue to suggest their influence on propensity of microstructural defect types to nucleate voids. A general void nucleation criterion is proposed using probability theory. The general framework is specified to polycrystalline BCC tantalum in the weak shock regime to include the SVE calculations and literature molecular dynamics calculations of grain boundary void nucleation strength. Probability density functions (PDFs) are used to describe the interaction between the local stress state heterogeneity and the distributed grain boundary void nucleation strength state. A causation entropy maximization procedure removes the requirement for ad hoc selection of a PDF functional form and provides a rigorous procedure for data-based PDF determination. The resulting physically informed PDF describes the spatial appearance frequency of nucleated voids as a function of applied macroscale pressure. Lower length scale physics are thus packaged in a precise and computationally efficient way to provide computational plasticity insight to macroscale dynamic ductile damage models.

36 MATERIALS SCIENCE↗