Search NASA⌕ Search

SEARCH · Search NASA

Results for “INTERNAL 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 505 records · Page 28

Kimberlina 1.2 CCUS Geophysical Models and Synthetic Data Sets

This synthetic multi-scale and multi-physics data set was produced in collaboration with teams at the Lawrence Berkeley National Laboratory, National Energy Technology Laboratory, Los Alamos National Laboratory, and Colorado School of Mines through the Science-informed Machine Learning for Accelerating Real-Time Decisions in Subsurface Applications (SMART) Initiative. Data are associated with the following publication: Alumbaugh, D., Gasperikova, E., Crandall, D., Commer, M., Feng, S., Harbert, W., Li, Y., Lin, Y., and Samarasinghe, S., “The Kimberlina Synthetic Geophysical Model and Data Set for CO2 Monitoring Investigations”, The Geoscience Data Journal, 2023, DOI: 10.1002/gdj3.191. The dataset uses the Kimberlina 1.2 CO2 reservoir flow model simulations based on a hypothetical CO2 storage site in California (Birkholzer et al., 2011; Wainwright et al., 2013). Geophysical properties models (P- and S-wave seismic velocities, saturated density, and electrical resistivity) were produced with an approach similar to that of Yang et al. (2019) and Gasperikova et al. (2022) for 100 Kimberlina 1.2 reservoir models. Links to individual resources are provided below: [CO2 Saturation Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-co2-saturation-models); Resistivity Models – [part 1](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-resistivity-models-part-1), [part 2](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-resistivity-models-part-2), and [part 3](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-resistivity-models-part-3); [Vp Velocity Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-vp-velocity-models); [Vs Velocity Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-vs-velocity-models); [Density Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-density-models). The 3D distributions of geophysical properties for the 33 time stamps of the SIM001 model were used to generate synthetic seismic, gravity, and electromagnetic (EM) responses for 33 times between zero and 200 years. Synthetic surface seismic data were generated using 2D and 3D finite-difference codes that simulate the acoustic wave equation (Moczo et al., 2007). 2D data were simulated for six point-pressure sources along a 2D line with 10 m receiver spacing and a time spacing of 0.0005 s. 3D simulations were completed for 25 surface pressure sources using a source separation of 1 km in both the x and y directions and a time spacing of 0.001 s. Links to individual resources are provided below: [2D velocity models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-2d-velocity-models) and [2D surface seismic data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-2d-surface-seismic-data). [3D velocity models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-velocity-models), and 3D seismic data [year0](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year0), [year1](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year1), [year2](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year2), [year5](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year5), [year10](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year10), [year15](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year15), [year20](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year20), [year25](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year25), [year30](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year30), [year35](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year35), [year40](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year40), [year45](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year45), [year49](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year49), [year50](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year50), [year51](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year51), [year52](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year52), [year55](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year55), [year60](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year60), [year65](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year65), [year70](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year70), [year75](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year75), [year80](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year80), [year85](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year85), [year90](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year90), [year95](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year95), [year100](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year100), [year110](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year110), [year120](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year120), [year130](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year130), [year140](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year140), [year150](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year150), [year175](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year175), [year200](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year200). The Python scripts to read these models and data are provided [here](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-python-scripts). EM simulations used a borehole-to-surface survey configuration, with the source located near the reservoir level and receivers on the surface using the code developed by Commer and Newman (2008). Pseudo-2D data for the source at [2500 m](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-pseudo-2d-csem-data-tz2500m) and [3025 m](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-pseudo-2d-csem-data-tz3025m), used a 2D inline receiver configuration to simulate a response over 3D resistivity models. The [3D data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-csem-data) contain electric fields generated by borehole sources at monitoring well locations and measured over a surface receiver grid. Vector gravity data, both on the surface and in boreholes, were simulated using a modeling code developed by Rim and Li (2015). The simulation scenarios were parallel to those used for the EM: [pseudo-2D data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-gravity-data) were calculated along the same lines and within the same boreholes, and [3D data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-gravity-data) were simulated over 3D models on the surface and in three monitoring wells. A series of [synthetic well logs](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-well-logs) of CO2 saturation, acoustic velocity, density, and induction resistivity in the injection well and three monitoring wells are also provided at 0, 1, 2, 5, 10, 15, and 20 years after the initiation of injection. These were constructed by combining the low-frequency trend of the geophysical models with the high-frequency variations of actual well logs collected in the Kimberlina 1 well that was drilled at the proposed site. Measurements of permeability and pore connectivity were made on cores of Vedder Sandstone, which forms the primary reservoir unit: [CT micro scans](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-ct-micro-scans-of-vedder-formation) and [Industrial CT Images](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-industrial-ct-images-vedder-formation). These measurements provide the range of scales in the otherwise synthetic data set to be as close to a real-world situation as possible. References: Birkholzer, J.T., Zhou, Q., Cortis, A. and Finsterle, S., 2011. A sensitivity study on regional pressure buildup from large-scale CO2 storage projects. Energy Procedia, 4, 4371-4378. Commer, M., and Newman, G.A., 2008. New advances in three-dimensional controlled-source electromagnetic inversion, Geophysical Journal International, 172, 513-535. Gasperikova, E., Appriou, D., Bonneville, A., Feng, Z., Huang, L., Gao, K., Yang, X., Daley, T., 2022, Sensitivity of geophysical techniques for monitoring secondary CO2 storage plumes, Int. J. Greenh. Gas Control, Volume 114, 103585, ISSN 1750-5836, https://doi.org/10.1016/j.ijggc.2022.103585. Moczo, P., J.O. Robertsson and L. Eisner, 2007, The finite-difference time-domain method for modeling of seismic wave propagation: Advances in geophysics, 48, 421-516. Rim, H., and Y. Li, 2015, Advantages of borehole vector gravity in density imaging, Geophysics, 80, G1-G13. Wainwright, H. M.; Finsterle, S.; Zhou, Q.; Birkholzer, J. T., 2013. Modeling the Performance of Large-Scale CO2 Storage Systems: A Comparison of Different Sensitivity Analysis Methods. International Journal of Greenhouse Gas Control, 17, 189205. https://doi.org/10.1016/j.ijggc.2013.05.007, DOI: 10.18141/1603331. Yang, X., Buscheck, T.A., Mansoor, K., Wang, Z., Gao, K., Huang, L., Appriou, D., and Carroll, S.A., 2019. Assessment of geophysical monitoring methods for detection of brine and CO2 leakage in drinking water aquifers, International Journal of Greenhouse Gas Control, 90, 102803, https://doi.org/10.1016/j.ijggc.2019.102803.

CCUS↗

Extreme compression of planetary gases: High-accuracy pressure-density measurements of hydrogen-helium mixtures above fourfold compression

Hydrogen (H 2 ) and helium (He), the most abundant elements in the universe, pose a unique challenge in measuring the equation of state of the mixture, owing to their differing physical properties. There remains a need for data with high enough precision to discriminate between existing equation of state (EOS) mix models in order to understand the internal structure of gas-giant planets. Here, we have measured the EOS of precompressed H 2 - He mixtures at conditions directly relevant to the planetary interiors using hypervelocity gas guns and Sandia’s Z machine with less than 10% uncertainty in density, enabling validation of mixture models. We precompressed 50:50 molar mixtures of H 2 -He to 0.1–0.2 GPa and directly measured particle velocity (in gas-gun experiments) and shock velocities (in Z-machine experiments). To complement the experimental efforts, we also computed the Hugoniots of precompressed H 2 -He mixtures using density-functional-theory-based molecular dynamics. Furthermore, we observe approximately 3- to 4.3-fold density compression at pressures up to 44 GPa.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Wind Turbine Sound Setbacks and Supply Curves: Ordinances and Extrapolated Trends, 110 Hub Height, 130 Rotor Diameter

This dataset provides a comprehensive set of wind turbine sound setbacks from every residential structure in the contiguous United States (CONUS). A sound setback is defined as the minimum required distance between a residential structure and a hypothetical turbine installation site to ensure that modeled sound levels received at the residence do not exceed local sound ordinances, which are commonly expressed in A-weighted decibels (dBA). Therefore, sound setbacks are a local spatial assessment combining multiple factors, including the sound pressure curve as a function of the observer location (distance and direction) relative to the turbine, local sound regulations, and the geographical distribution of residential structures. The dataset is organized into multiple scenario-based products, detailed as follows: 1. Existing and extrapolated sound setbacks. An existing scenario characterizes sound setbacks only in states or counties that have implemented sound regulations as of 2022. The extrapolated scenarios extend a constant sound threshold to counties that lack explicit sound regulations, with thresholds ranging from 35 to 60 dBA, in 5-dBA increments reflecting the variation observed in current sound ordinances. 2. Sound setbacks in directional and worst scenarios. The directional scenario accounts for the distance and orientation of residential structures relative to a hypothetical turbine location, utilizing the turbine's sound emissions in that specific direction. In contrast, the worst scenario takes loudest sound level at each distance step from the turbine, irrespective of directional considerations, which aligns with current industry practice. 3. Supply curves for Open and Reference Access scenarios. This dataset includes supply curves generated by the reV model, which integrates each of the above sound setbacks into both Open and Reference siting scenarios. In addition, two Open and Reference baselines scenarios were included which do not consider sound setbacks for comparative analysis. All sound setback data are stored in TIF files, with partial maps of the data provided in PNG format. The values in the sound setback raster range from 0 to 1, representing the fraction of developable land within a 90 meter by 90 meter pixel due to sound ordinances. A value of 0 indicates areas where wind energy development is prohibited, while a value of 1 signifies areas fully permissible. The wind turbine parameters used in the sound modeling are based on the land-based turbine from International Energy Agency (IEA), featuring a rated electrical power of 3.4 MW, a rotor diameter of 130 meters, and a hub height of 110 meters. The atmospheric conditions, including wind speed/direction, turbulence, air temperature, relative humidity, and air pressure, that drive the sound generation are obtained from the WIND Toolkit dataset.

Array↗

Effect of microstructural and environmental variables on ductility of austenitic stainless steels

Austenitic stainless steels are used extensively in harsh environments, including for high-pressure gaseous hydrogen service. However, the tensile ductility of this class of materials is very sensitive to materials and environmental variables. While tensile ductility is generally insufficient to qualify a material for hydrogen service, ductility is an effective tool to explore microstructural and environmental variables and their effects on hydrogen susceptibility, to inform understanding of the mechanisms of hydrogen effects in metals, and to provide insight to microstructural variables that may improve relative performance. In this study, hydrogen precharging was used to simulate high-pressure hydrogen environments to evaluate hydrogen effects on tensile properties. Several austenitic stainless steels were considered, including both metastable and stable alloys. Room temperature and subambient temperature tensile properties were evaluated with three different internal hydrogen contents for type 304L and 316L austenitic stainless steels and one hydrogen content for XM-11. Significant ductility loss was observed for both metastable and stable alloys, suggesting the stability of the austenitic phase is not sufficient to characterize the effects of hydrogen. Internal hydrogen does influence the character of deformation, which drives local damage accumulation and ultimately fracture for both metastable and stable alloys. While a quantitative description of hydrogen-assisted fracture in austenitic stainless steels remains elusive, these observations underscore the importance of the hydrogen-defect interactions and the accumulation of damage at deformation length scales.

36 MATERIALS SCIENCE↗

Heat transfer optimization of uo 2 -mo fuel using genetic algorithms

Two genetic algorithm (GA) methods were applied to thermal finite element models to optimize the heat transfer efficacy of a UO 2 -Mo composite fuel pellet with typical pressurized water reactor fuel geometry. Mo additions to UO 2 have been shown to increase the thermal conductivity, thus reducing centerline temperatures and temperature gradients. Previous studies evaluated uniformly dispersed Mo or continuous Mo internal geometries (e.g., fins, plates, discs) that were selected using engineering intuition. The current study uses two different implementations of the same GA to optimize Mo placement and minimize the fuel temperature with the only constraint being a maximum 10% Mo volume fraction. One approach superimposed Mo line elements onto the monolithic UO 2 pellet model, and the other converted entire UO 2 volume elements to Mo. The former method generated 1D heat transfer connections between nodes, whereas the latter method allowed for the formation of 3D structures. Features of the optimal fuel design produced by the GAs included dispersed Mo near the centerline that shifted the peak fuel temperature outward by 0.6 mm, Mo chains in the high-heat-flux region in the mid-to-outer radial zone, and a large continuous structure that spanned the full radius and height of the pellet and accounted for 87.7 % of the total Mo in the pellet. Analysis of this design indicates that the optimal Mo configuration is a balance between creating continuous heat transfer pathways and optimally dispersing Mo to minimize the heat transfer distance through UO 2 . This architecture ultimately produced an effective thermal conductivity of 11.3 W/m·K under the assumed boundary conditions. This result is higher than any previous values from the literature. In conclusion, potential fabrication methods and challenges are discussed in addition to the implications on fuel performance.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Understanding How Ligand Functionalization Influences CO2 and N2 Adsorption in a Sodalite Metal–Organic Framework

In this work, a detailed study is conducted to understand how ligand substitution influences the CO2 and N2 adsorption properties of two highly crystalline sodalite metal–organic frameworks (MOFs) known as Cu–BTT (BTT–3 = 1,3,5-benzenetristetrazolate) and Cu–BTTri (BTTri–3 = 1,3,5-benzenetristriazolate). The enthalpy of adsorption and observed adsorption capacities at a given pressure are significantly lower for Cu–BTTri compared to its tetrazole counterpart, Cu–BTT. In situ X-ray and neutron diffraction, which allow visualization of the CO2 and N2 binding sites on the internal surface of Cu–BTTri, provide insights into understanding the subtle differences. As expected, slightly elongated distances between the open Cu2+ sites and surface-bound CO2 in Cu–BTTri can be explained by the fact that the triazolate ligand is a better electron donor than the tetrazolate. The more pronounced Jahn–Teller effect in Cu–BTTri leads to weaker guest binding. The results of the aforementioned structural analysis were complemented by the prediction of the binding energies at each CO2 and N2 adsorption site by density functional theory calculations. In addition, variable temperature in situ diffraction measurements shed light on the fine structural changes of the framework and CO2 occupancies at different adsorption sites as a function of temperature. Finally, simulated breakthrough curves obtained for both sodalite MOFs demonstrate the materials’ potential performance in dry postcombustion CO2 capture. The simulation, which considers both framework uptake capacity and selectivity, predicts better separation performance for Cu–BTT. The information obtained in this work highlights how ligand substitution can influence adsorption properties and hence provides further insights into the material optimization for important separations.

adsorption capacities↗

Modeling In-Reactor Temperature Gradients in Lithium Aluminate Ceramic Pellets

In reactor thermal performance of lithium aluminate ceramic pellets are calculated to determine operating temperatures if internal components shift during irradiation. A 2D conduction model was setup to model this scenario with additional physics, including thermal expansion of components to reactor temperatures, creep down of the cladding due to differential pressure between the inside of the rod and reactor coolant system, and net radiation exchange between components. The non-linear coupling of the governing equations presented convergence challenges due to the very different effects from a small temperature change on the radiation vs radial and circumferential conduction equations. Initial temperature guesses for the multidimensional solution were provided from conservative 1D models simulating bounding scenarios. Convergence was greatly improved by taking small steps towards the root solutions and averaging with initial guesses. The final method robustly converged all 132 points at an axial plane after a few dozen iterations. The temperature results show that the additional physics have a secondary effect as the dominant heat transfer mechanism is radial conduction across the gas gap followed by circumferential conduction around components. Circumferential conduction in the lithium aluminate pellets is enhanced by the relatively high thermal conductivity of this ceramic.

Carstens, Nathan A.↗

ECAR-5067 Computational Evaluation on Effect of Edge-rounding in Circular to Annular Flow Transitions within a Hexagonal Duct

In this ECAR, the computational fluid dynamics (CFD) simulation of the duct with a cross-section varying from circle to annulus was performed to evaluate the effect of edge-rounding of the internal edges. The edge of circular channel and the edge of inner annular channel were rounded with the radius of curvature ranged from 0 mm to 20 mm. The CFD results show that the pressure loss of the duct can be reduced by increasing the radius of the rounded edge. Case-A results showed that the pressure drop can be significantly reduced by the edge-rounding, e.g., the pressure drop was reduced approximately 56% by the radius of curvature of 5.0mm, and reduced 70% by the radius of curvature of 20.0 mm although its effectiveness would be attenuated as the radius of curvature of edge increases. Results from Case-B showed the impact of rounding the inlet edge was practically negligible compared to the outlet edge. This ECAR could provide a guideline to determine the optimum inlet module design of the VTR fuel rod bundle.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Internal combustion engine as a chemical reactor to produce synthesis gas from hydrocarbon feeds

An internal combustion engine is operated at fuel-rich conditions by adjusting one or more operating parameters such as, for example, a throttle, an ignition timing, a load coupled to the engine, a fuel pressure, power to a supercharger, and power to a preheater to maintain a specified engine speed and a temperature of an exhaust gas. Operating the engine under these conditions allows the engine to function as a reformer producing a synthesis gas comprising hydrogen and carbon monoxide.

Carpenter, III, John Reeves↗

High-Temperature Particle Heat Exchanger for sCO2 Power Cycles [Award 30342]

This report describes the design, development, and testing of a prototype 100 kWt particle-to-supercritical CO 2 (sCO2) heat exchanger. An analytic hierarchy process was implemented to compare and evaluate alternative heat-exchanger designs (fluidized bed, shell-and-plate moving packed bed, and shell-and-tube moving packed bed) that could meet the high pressure (≥ 20 MPa) and high temperature (≥ 700 °C) operational requirements associated with sCO2 power cycles. Cost, heat-transfer coefficient, structural reliability, manufacturability, parasitics and heat losses, scalability, compatibility, erosion and corrosion, transient operation, and inspection ease were considered in the evaluation. A 100 kWt shell-and-plate design was selected for construction and integration with Sandia’s falling particle receiver system that heats the particles using concentrated sunlight. Sandia worked with industry to design and construct the moving packed-bed shell-and-plate heat exchanger. Tests were performed to evaluate its performance using both electrical heating and concentrated sunlight to heat the particles. Overall heat transfer coefficients at off-design conditions (reduced operating temperatures and only three stainless steel banks in the counter-crossflow heat exchanger) were measured to be approximately ~25 - 70 W/m 2 -K, significantly lower than simulated values of >100 W/m 2 -K. Tests using the falling particle receiver to heat the particles with concentrated sunlight yielded overall heat transfer coefficients of ~35 – 80 W/m 2 -K with four banks (including a nickel-alloy bank above the three stainless steel banks). The overall heat transfer coefficient was observed to decrease with increasing particle inlet temperatures, which contrasted the results of simulations that showed an increase in heat transfer coefficient with temperature due to increased effective particle-bed thermal conductivity from radiation. The likely cause of the discrepancy was particle-flow maldistributions and funnel flow within the heat exchanger caused by internal ledges and cross-bracing, which could have been exacerbated by increased particle-wall friction at higher temperatures. Additional heat loss at higher temperatures may also contribute to a lower overall heat-transfer coefficient. Design challenges including pressure drop, particle and sCO2 flow maldistribution, and reduced heat transfer coefficient are discussed with approaches for mitigation in future designs. Lessons learned regarding instrumentation, performance characterization, and operation of particle components and sCO2 flow loops are also discussed. Finally, a 200 MWt commercial-scale shell-and-plate heat-exchanger design based on the concepts investigated in this report is proposed.

14 SOLAR ENERGY↗

All‐Solid‐State Batteries With Mechanically Stable Interfaces Consisting of a Zero‐Strain Cation‐Disordered Rocksalt Cathode

Interface stabilization is critical to the development of working all‐solid‐state batteries. Rigid cathode/solid electrolyte interfaces often disintegrate due to anisotropic volume change of cathode‐active materials, resulting in irreversible capacity loss. Herein, we demonstrate that Li 1.211 Mo 0.467 Cr 0.3 O 2 (LMCO), a pioneering cation‐disordered rocksalt oxide (DRX) cathode that has intrinsically small volume change upon lithium intercalation, can be integrated with a thiophosphate‐based solid electrolyte for all‐solid‐state batteries. Interface stability of the all‐solid LMCO cell was investigated by electrochemical impedance spectroscopy, X‐ray micro‐computed tomography, and electron microscopy. Since LMCO was initially synthesized as a layered phase exhibiting a large volume change, interface disintegration can be observable in the early cycles. As layered LMCO phase‐transformed into DRX LMCO in subsequent cycles, reintegration of the interfaces occurs within a pressurized cell as a result of its zero‐stain behavior. Consequently, the DRX LMCO cathode maintains interface integrity, and thus electrical wiring, over an extended number of cycles, leading to improved capacity retention with small internal cell resistance.

DRX↗

Dynamic fracture of glass fiber-reinforced ductile polymer matrix composites and loading rate effect

Here, the dynamic fracture of S-2 glass fiber-reinforced polymer matrix composites (FRPMCs) was investigated in this study. The matrix ductility was improved by a recently developed network topology modification technique via mixing partially reacted substructures (mPRS). The composite material was manufactured and characterized by micro-CT scanning and scanning electron microscopy (SEM). Dynamic single-edge notched bending (d-SENB) experiments were performed on the composites by using a modified split-Hopkinson pressure bar. Each specimen’s fracture process was visualized by ultrafast X-ray imaging. Such in-situ radiography enabled identifying the damage initiation below 50-micron scale and inspecting its propagation through the internal structures of opaque composites, thereby accurately quantifying the composites’ mechanical properties. Furthermore, the identical d-SENB experiments were designed and the digital image correlation (DIC) was employed to monitor the stress wave propagation on the composite specimens. The force and deflection measurements were modified and correlated to the physical damage processes. Besides, quasi-static SENB experiments were conducted to identify the loading rate effects on the composites’ fracture behaviors. The force and deflection history, bending stiffness, energy dissipation, and fracture toughness at different loading rates were quantified and compared. Finally, post-fracture analysis by micro-CT scanning and SEM provided physical observations on the variation of the fracture morphology by different loading rates.

42 ENGINEERING↗

Efficient and robust phase-split computations in the internal energy, volume, and moles ( UVN ) space

Phase-split computations in an isolated system, which in general may be defined as one where the total internal energy (U), volume (V), and the number of moles (N = N 1 , N 2 ,..., N n ) of the components are fixed at some specified set of values, involves the determination of the temperature (T), pressure (P), plus the amount and composition of the various phases that constitute the system. A simpler, but analogous problem is one where T, P, and N are specified instead. In TPN space, one may first perform a stability analysis to determine whether the system is stable, meaning whether at equilibrium it will split up into multiple phases or remain in single phase. If the single-phase state is unstable, the stability analysis reliably provides a good set of initial guesses in the subsequent phase-split computations. In UVN space, however, we demonstrate that the stability analysis (which is the main subject of our earlier study [1]) may not in general provide good enough initial guesses; we offer alternative strategies for setting up good initial guesses. Furthermore, we show that a combination of successive substitution iteration (SSI) and Newton's method---two iterative methods that are prevalent in the literature in TPN space---facilitates a robust and efficient algorithm for phase-split computations in isolated systems. Finally, this combination has so far not been applied in UVN space.

02 PETROLEUM↗

Entropy generation from hydrodynamic mixing in inertial confinement fusion indirect-drive targets

The increase in entropy from the physical mixing of two adjacent materials in inertial confinement fusion (ICF) implosions and gas-filled hohlraums is analytically assessed. An idealized model of entropy generation from the mixing of identical ideal-gas particles across a material interface in the presence of pressure and temperature gradients is applied. Physically, mix-driven entropy generation refers to the work done by the gases in expanding into a larger common volume from atomic mixing under the condition of no internal energy change, or work needed to restore the initial unmixed state. Furthermore, the effect of a mix-generated entropy increase is analytically shown to lead to less compression of the composite ICF fluid under adiabatic conditions. The amount of entropy generation is estimated to be ~10 J for a Rayleigh–Taylor-induced micrometer-scale annular mixing layer between the solid deuterium–tritium fuel and (undoped) high-density carbon pusher of an imploding capsule at the National Ignition Facility (NIF). This level of entropy generation is consistent with lower-than-expected fuel compressions measured on the NIF [Hurricane et al., Phys. Plasmas 26, 052704 (2019)]. The degree of entropy increase from mixing of high-Z hohlraum wall material and low-Z, moderate- to high-density gas fills is estimated to lead to ~100 kJ of heat generation for NIF-scale experiments [Moody et al., Phys. Plasmas 21, 056317 (2014)]. This value represents a significant fraction of the inferred missing x-ray drive energy based on observed delays in capsule implosion times compared with mainline simulations [Jones et al., Phys. Plasmas 19, 056315 (2012)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Instrumentation for the In-Core Real-Time Mechanical Testing of Structural Materials (INCREASE) Project

Idaho National Laboratory (INL), in collaboration with the Electric Power Research Institute (EPRI), the Nuclear Regulatory Commission (NRC), the French Atomic and Alternative Energies Commission (CEA), the Joint Research Center (JRC), the Nuclear Research and Consultancy Group (NRG), and the Research Center Rez (CVR), started a Joint Experimental Program (JEEP) project that operates within the Nuclear Energy Agency’s Framework for Irradiation Experiments (FIDES II) program in order to develop capabilities for the in-core real-time mechanical testing of structural materials. This effort will focus on designing a shared capsule capable of housing a variety of in-core mechanical testing instrumentation allowing enhanced experiments for the material science community. The outcome of the project would be high-priority, stress relaxation data for stainless-steel-based materials provided by EPRI and CEA. Stress relaxation is a major phenomenon that contributes to material degradation in nuclear reactor components. Currently, nuclear material stress relaxation is assessed both before and after irradiation, using complex and costly post-irradiation examination (PIE) activities. In-situ data would support the development of precision modeling and simulation of this degradation phenomena and would provide validation and benchmarking for existing models using the PIE data. As part of the U.S. Department of Energy Advanced Sensor and Instrumentation (ASI) program, INL has fabricated and tested out-of-core mechanical test instrumentation. This instrumentation was designed for easy adaptation to the irradiation capsule proposed under this JEEP, and can be deployed to measure real-time stress relaxation under pressurized-water reactor (PWR) conditions. This initial effort will partially serve to replace the testing capabilities lost because of shutting down the Halden Boiling Water Reactor (HBWR). The project would provide these capabilities to the international community via a shared capsule design that is easily adaptable to additional material test reactors. The design features will incorporate expansion to PWR and non-light-water reactor (LWR) environments that will be developed in future work. The capsule and instrumentation will be demonstrated in the Massachusetts Institute of Technology Reactor (MITR) for phase I and the Petten High Flux Reactor (HFR) for phase II irradiations to deliver real-time stress relaxation data on high priority stainless steel structural materials.

36 MATERIALS SCIENCE↗

A Pilot-Scale Evaluation of Natural Gas-Based Foam at Elevated Pressure and Temperature Conditions

Abstract Nitrogen (N2) and Carbon Dioxide (CO2) foams have been used as hydraulic fracturing fluids for several decades to reduce water usage and minimize damage in water-sensitive reservoirs. These foam treatments require gases to be liquefied and transported to site. An alternative approach would be to use natural gas (NG) that is readily available from nearby wells, pipelines, and processing facilities as the internal, gaseous phase to create a NG-based foam. Hydraulic fracturing with NG foam is a relatively inexpensive option, makes use of an abundant and often wasted resource, and may even provide production benefits in certain reservoirs. As part of an ongoing development project sponsored by the Department of Energy (DOE), the surface process to create NG foam is being developed and the properties of NG foam are being explored. This paper presents recent results from a rigorous pilot-scale demonstration of NG foam over a range of operating scenarios relevant to surface and bottomhole conditions with a variety of base-fluid mixtures. The Pilot-scale Foam Test Facility (PFTF) used in these investigations is first described. The PFTF is capable of generating foamed fluids at pressures up to 7,500 psig and at temperatures in excess of 300°F. Then, results from several investigations aimed at proving NG foam at conditions relevant to the field are presented. NG foam was characterized using rheology measurements and flow visualization techniques. Experiments were performed to investigate the texture and stability of NG foam generated by two different mixing methods: one using a custom designed tee to match mixing velocities in the field where the gas phase is jetted into the aqueous stream, and another to ensure comprehensive mixing for laboratory analysis. Parametric studies were conducted to explore the effects of flow rate, foam quality, and temperature on the stability of NG foam. Moreover, different fluid preparations were used to investigate the effect of base fluid and additive concentrations on the stability of NG foams. Additional laboratory studies that investigated foam stability with produced water and multicomponent NG mixtures are also reported. The NG foams explored in these investigations exhibited typical, shear-thinning behavior observed in rheological studies of N2- and CO2-based foams. The measured viscosity and observed stability indicate that NG foams are well suited for fracturing applications. Like other foams, NG foam exhibits sensitivity to operating temperature characterized by a decrease in apparent viscosity as temperature increases. Rapid foam breakdown was observed at significantly elevated temperatures exceeding 290°F. In addition to fluid characterization, these investigations also yielded several key lessons that should be applied to future field demonstrations of NG foam.

Beck, Griffin↗

New Dimension in Ab Initio Electronic Structure Theory: Temperature, Pressure, and Chemical Potential

Ab initio electronic structure theory has transformed gas-phase molecular science with its predictive ability. In the attempt to bring such predictive ability to macroscopic systems and condensed matter, the theory must integrate quantum mechanics with statistical thermodynamics, so that thermodynamic functions such as free energy, internal energy, entropy, and chemical potentials are computed as functions of temperature in a systematically converging series of approximations. Here, a general, versatile strategy of elevating ab initio electronic structure theory to nonzero temperatures is introduced and discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗

Ferroelectricity in hafnia controlled via surface electrochemical state

Ferroelectricity in binary oxides including hafnia and zirconia has riveted the attention of the scientific community due to the highly unconventional physical mechanisms and the potential for the integration of these materials into semiconductor workflows. Over the last decade, it has been argued that behaviours such as wake-up phenomena and an extreme sensitivity to electrode and processing conditions suggest that ferroelectricity in these materials is strongly influenced by other factors, including electrochemical boundary conditions and strain. Here we argue that the properties of these materials emerge due to the interplay between the bulk competition between ferroelectric and structural instabilities, similar to that in classical antiferroelectrics, coupled with non-local screening mediated by the finite density of states at surfaces and internal interfaces. Via the decoupling of electrochemical and electrostatic controls, realized via environmental and ultra-high vacuum piezoresponse force microscopy, we show that these materials demonstrate a rich spectrum of ferroic behaviours including partial-pressure-induced and temperature-induced transitions between ferroelectric and antiferroelectric behaviours. In conclusion, these behaviours are consistent with an antiferroionic model and suggest strategies for hafnia-based device optimization.

36 MATERIALS SCIENCE↗