Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Integrated Simulation”

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

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

At least 1,027 records · Page 57

Phonons reveal coupled cholesterol-lipid dynamics in ternary membranes

Experimental studies of collective dynamics in lipid bilayers have been challenging due to the energy resolution required to observe these low-energy phonon-like modes. However, inelastic X-ray scattering (IXS) measurements — a technique for probing vibrations in soft and biological materials — are now possible with sub-meV resolution, permitting direct observation of low energy, phonon-like modes in lipid membranes. Here, in this study, IXS measurements with sub-meV energy resolution reveal a low-energy optic-like phonon mode at roughly 3 meV in the liquid-ordered (L o ) and liquid-disordered (L d ) phases of a ternary lipid mixture. This mode is only observed experimentally at momentum transfers greater than 5 nm -1 in the L o system. A similar gapped mode is also observed in all-atom molecular dynamics (MD) simulations of the same mixture, indicating that the simulations accurately represnt the fast, collective dynamics in the L o phase. It’s optical nature and the Q range of the gap together suggest that the observed mode is due to the coupled motion of cholesterol-lipid pairs, separated by several hydrocarbon chains within the membrane plane. Analysis of the simulations provide molecular insight into the origin of the mode in transient, nanoscale substructures of hexagonally packed hydrocarbon chains. This nanoscale hexagonal packing was previously reported based on molecular dynamics simulations and later by NMR measurements. Here, however, the integration of IXS and MD simulations identifies a new signature of the L o substructure in the collective lipid dynamics, thanks to the recent confluence of IXS sensitivity and MD simulation capabilities.

36 MATERIALS SCIENCE↗

Physical properties of the inner shocks in hot, tilted black hole accretion flows

Simulations of hot, pressure-supported, tilted black hole accretion flows, in which the angular momentum of the flow is misaligned with the black hole spin axis, can exhibit two nonaxisymmetric shock structures in the inner regions of the flow. We analyze the strength and significance of these shock structures in simulations with tilt angles of 10° and 15°. By integrating fluid trajectories in the simulations through the shocks and tracking the variations of fluid quantities along these trajectories, we show that these shocks are strong, with substantial compression ratios, in contrast to earlier claims. However, they are only moderately relativistic. We also show that the two density enhancements resembling flow streams in their shape are in fact merely post-shock compressions, as fluid trajectories cut across, rather than flow along, them. The dissipation associated with the shocks is a substantial fraction (≅ 3%-12%) of the rest mass energy advected into the hole, and therefore comparable to the dissipation expected from turbulence. The shocks should therefore make order unity changes in the observed properties of black hole accretion flows that are tilted.

79 ASTRONOMY AND ASTROPHYSICS↗

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

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

54 ENVIRONMENTAL SCIENCES↗

Experimental Design and Laboratory Characterization of a Medium- and High-Efficiency Rooftop Unit for use in Building Energy Simulations

Commercial rooftop units that incorporate variable speed components are quickly entering into the market place; however, the design and modeling tools that calculate the energy benefits of these units cannot properly estimate performance. This is due to two factors: 1) There are insufficient measured performance data of real systems that can be entered into these simulations, and 2) building energy simulators are unable to properly interpret performance data. Therefore, publicly available datasets designed to encompass the entire performance map coupled with control sequence information are needed as inputs for these modelling programs to accurately estimate energy savings and promote variable speed RTUs. As part of the Verification and Uncertainty Characterization for Energy Simulation milestone funded by the Department of Energy, the National Renewable Energy Laboratory recorded and published data across the entire performance range for two commercially available RTUs; one with a two-stage compressor and variable speed supply fan, and another with a fully variable speed compressor, supply fan and condenser fan. This paper details 1) how the design of experiments was created as well as how the units were controlled, and 2) the operational limitations of the two units. The paper will also discuss methods for building energy simulation programs to properly interpret these data sets and integrate them into building energy simulation models.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Enforcing elemental mass and energy balances for reduced order models

Development of economically feasible gasification and carbon capture, utilization and storage (CCUS) technologies requires a variety of software tools to optimize the designs of not only the key devices involved (e., g., gasifier, CO{sub 2} adsorber) but also the entire power generation system. High-fidelity models such as Computational Fluid Dynamics (CFD) models are capable of accurately simulating the detailed flow dynamics, heat transfer, and chemistry inside the key devices. However, the integration of CFD models within steady-state process simulators, and subsequent optimization of the integrated system, still presents significant challenges due to the scale differences in both time and length, as well the high computational cost. A reduced order model (ROM) generated from a high-fidelity model can serve as a bridge between the models of different scales. While high-fidelity models are built upon the principles of mass, momentum, and energy conservations, ROMs are usually developed based on regression-type equations and hence their predictions may violate the mass and energy conservation laws. A high-fidelity model may also have the mass and energy balance problem if it is not tightly converged. Conservations of mass and energy are important when a ROM is integrated to a flowsheet for the process simulation of the entire chemical or power generation system, especially when recycle streams are connected to the modeled device. As a part of the Carbon Capture Simulation Initiative (CCSI) project supported by the U.S. Department of Energy, we developed a software framework for generating ROMs from CFD simulations and integrating them with Process Modeling Environments (PMEs) for system-wide optimization. This paper presents a method to correct the results of a high-fidelity model or a ROM such that the elemental mass and energy are conserved perfectly. Correction factors for the flow rates of individual species in the product streams are solved using a minimization algorithm based on Lagrangian multiplier method. Enthalpies of product streams are also modified to enforce the energy balance. The approach is illustrated for two ROMs, one based on a CFD model of an entrained-flow gasifier and the other based on the CFD model of a multiphase CO{sub 2} adsorber.

01 COAL, LIGNITE, AND PEAT↗

Failure Diagnosis for the Holdup Tank System via ISFA

This paper discusses the use of the integrated system failure analysis (ISFA) technique for fault diagnosis for the holdup tank system. ISFA is a simulation-based, qualitative and integrated approach used to study fault propagation in systems containing both hardware and software subsystems. The holdup tank system consists of a tank containing a fluid whose level is controlled by an inlet valve and an outlet valve. We introduce the component and functional models of the system, quantify the main parameters and simulate possible failure-propagation paths based on the fault propagation approach, ISFA. The results show that most component failures in the holdup tank system can be identified clearly and that ISFA is viable as a technique for fault diagnosis. Since ISFA is a qualitative technique that can be used in the very early stages of system design, this case study provides indications that it can be used early to study design aspects that relate to robustness and fault tolerance. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Simulating a solar parabolic trough collector plant used for industrial process heat using an optimized operating scheme that utilizes flexible heat integration

This study focuses on the modeling and simulation of a novel design and operation of a solar industrial process heat (SIPH) plant that uses a parabolic trough collector system for generating process heat. The SIPH plant incorporates flexible heat integration (FHI) by having two options for heat sinks as well as a flexible collection temperature for the heat transfer fluid. Leveraging the degrees of freedom created by FHI allows an optimized operating scheme to be created which maximizes the solar heat delivered to the industrial processes. Yearly results for a 27 MW-t plant located in Salt Lake City, UT are simulated for a base case and a case utilizing optimized FHI. For the main process temperature case presented, FHI increases the total solar share by 13.0% and reduces the levelized cost of heat of the SIPH plant by 10.4% relative to the base case. Furthermore, three case study days are presented which highlight how FHI maximizes solar share, specifically during mornings, afternoons, and days with intermittent solar conditions. An environmental analysis is also performed showing that the SIPH plant can reduce harmful emissions by 15.4% for the base case and 17.4% for the FHI case compared to a plant that only utilizes a natural gas furnace for process heat. Overall, the work thoroughly explores an effective operating scheme for a SIPH plant, helping to make it more economical and environmentally beneficial.

14 SOLAR ENERGY↗

Ice and water droplets on graphite: A comparison of quantum and classical simulations

Ice and water droplets on graphite have been studied by quantum path integral and classical molecular dynamics simulations. The point-charge q-TIP4P/F potential was used to model the interaction between flexible water molecules, while the water-graphite interaction was described by a Lennard-Jones potential previously used to reproduce the macroscopic contact angle of water droplets on graphite. Several energetic and structural properties of water droplets with sizes between 10{sup 2} and 10{sup 3} molecules were analyzed in a temperature interval of 50–350 K. The vibrational density of states of crystalline and amorphous ice drops was correlated to the one of ice Ih to assess the influence of the droplet interface and molecular disorder on the vibrational properties. The average distance of covalent OH bonds is found 0.01 Å larger in the quantum limit than in the classical one. The OO distances are elongated by 0.03 Å in the quantum simulations at 50 K. Bond distance fluctuations are large as a consequence of the zero-point vibrations. The analysis of the H-bond network shows that the liquid droplet is more structured in the classical limit than in the quantum case. The average kinetic and potential energy of the ice and water droplets on graphite has been compared with the values of ice Ih and liquid water as a function of temperature. The droplet kinetic energy shows a temperature dependence similar to the one of liquid water, without apparent discontinuity at temperatures where the droplet is solid. However, the droplet potential energy becomes significantly larger than the one of ice or water at the same temperature. In the quantum limit, the ice droplet is more expanded than in a classical description. Liquid droplets display identical density profiles and liquid-vapor interfaces in the quantum and classical limits. The value of the contact angle is not influenced by quantum effects. Contact angles of droplets decrease as the size of the water droplet increases which implies a positive sign of the line tension of the droplet.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Enhancing Modeling and Simulation for Effective Protection Strategies

This report was created by Sandia National Laboratories (SNL) to document the principles and methodology of performance data collection and integration with modeling and simulation tools to better facilitate the performance evaluation of physical protection systems (PPS). Results and conclusions from the use of modeling and simulation tools are only as good as the data employed by the tools when conducting analysis. Acquiring the performance testing data necessary to ensure effective evaluation can be a complex and sometimes daunting process. It is the desire of the organization to provide guidance that eases the burdens associated with pursuit of these objectives. This document draws heavily upon longstanding principles of systems engineering that have been developed and employed by SNL in the discipline of security since the 1970s. The scope of this document is constrained to the testing of system components and integration of data that is applicable within the context of PPS performance analysis using two tools that have been developed by SNL, PathTrace© and Scribe3D©. For guidance related to testing and evaluation more broadly, the manuals and reports referenced by this document can be consulted.

42 ENGINEERING↗

ALTRIOS - Advanced Locomotive Technology and Rail Infrastructure Optimization System

The Advanced Locomotive Technology and Rail Infrastructure Optimization System (ALTRIOS) is a unique, fully integrated, open-source software tool to evaluate strategies for deploying advanced locomotive technologies and associated infrastructure for cost-effective decarbonization. ALTRIOS simulates energy conversion and storage dynamics, locomotive and train dynamics, meet-pass planning (detailed train timetabling), and freight-demand-driven train scheduling in a Pareto optimization. Because new locomotives represent a significant long-term capital investment and new technologies must be thoroughly demonstrated before deployment, this tool provides guidance on the risk/reward tradeoffs of different technology rollout strategies. An open, integrated simulation tool will be invaluable for identifying future research needs and making decisions on technology development, routes, and train selection. ALTRIOS development is a collaboration between NREL, the University of Illinois Urbana-Champaign (UIUC) RailTEC, Southwest Research Institute (SwRI), and BNSF.

ADVANCED PROPULSION SYSTEMS,MATHEMATICS AND COMPUT↗

SYSTEMATICS IN METALLICITY GRADIENT MEASUREMENTS. I. ANGULAR RESOLUTION, SIGNAL TO NOISE, AND ANNULAR BINNING

With the rapid progress in metallicity gradient studies at high redshift, it is imperative that we thoroughly understand the systematics in these measurements. This work investigates how the [N II]/H{alpha}-ratio-based metallicity gradients change with angular resolution, signal to noise (S/N), and annular binning parameters. Two approaches are used: (1) we downgrade the high angular resolution integral-field data of a gravitationally lensed galaxy and re-derive the metallicity gradients at different angular resolution; (2) we simulate high-redshift integral field spectroscopy observations under different angular resolution and S/N conditions using a local galaxy with a known gradient. We find that the measured metallicity gradient changes systematically with angular resolution and annular binning. Seeing-limited observations produce significantly flatter gradients than higher angular resolution observations. There is a critical angular resolution limit beyond which the measured metallicity gradient is substantially different to the intrinsic gradient. This critical angular resolution depends on the intrinsic gradient of the galaxy and is {<=}0.''02 for our simulated galaxy. We show that seeing-limited high-redshift metallicity gradients are likely to be strongly affected by resolution-driven gradient flattening. Annular binning with a small number of annuli produces a more flattened gradient than the intrinsic gradient due to weak line smearing. For three-annulus bins, a minimum S/N of {approx}5 on the [N II] line is required for the faintest annulus to constrain the gradients with meaningful errors.

79 ASTRONOMY AND ASTROPHYSICS↗

ALTRIOS (Advanced Locomotive Technology and Rail Infrastructure Optimization System) [SWR-22-54]

The Advanced Locomotive Technology and Rail Infrastructure Optimization System (ALTRIOS) is a unique, fully integrated, open-source software tool to evaluate strategies for deploying advanced locomotive technologies and associated infrastructure for cost-effective decarbonization. ALTRIOS simulates freight-demand driven train scheduling, mainline meet-pass planning, locomotive dynamics, train dynamics, energy conversion efficiencies, and energy storage dynamics of line-haul train operations. Because new locomotives represent a significant long-term capital investment and new technologies must be thoroughly demonstrated before deployment, this tool provides guidance on the risk/reward tradeoffs of different technology rollout strategies. An open, integrated simulation tool is invaluable for identifying future research needs and making decisions on technology development, routes, and train selection. ALTRIOS was developed as part of a collaborative effort by a team comprising the National Renewable Energy Laboratory (NREL), University of Illinois Urbana-Champaign (UIUC), Southwest Research Institute (SwRI), and BNSF Railway. Python Package: https://pypi.org/project/altrios/ Rust Crate: https://crates.io/crates/altrios-core

Baker, Chad↗

Geomechanical Framework for Secure CO 2 Storage in Fractured Reservoirs and Caprocks for Sedimentary Basins in theMidwest United States

This report presents final technical results for the project Geomechanical Framework for Secure CO 2 Storage in Fractured Reservoirs and Caprocks for Sedimentary Basins in the Midwest United States (DE-FE0023330). The project was a three-year effort consisting of seven technical tasks focused on defining geomechanical factors for CO 2 storage applications in deep saline rock formations in Ohio and the Midwest United States, because geomechancial issues have been identified as a significant risk factor for large-scale CO 2 storage applications. A basin-scale stress-strain analysis was completed to describe the geomechanical setting for rock formations of Ordovician-Cambrian age in Ohio and adjacent areas of the Midwest United States in relation to geologic CO 2 storage applications. The tectonic setting, stress orientation-magnitude, and geomechanical and petrophysical parameters for CO 2 storage zones and caprocks in the region were cataloged. Ten geophysical image logs were analyzed for natural fractures, borehole breakouts, and drilling-induced fractures. The logs indicated mostly less than 10 fractures per 100 vertical feet in the borehole, with mostly N65E principal stress orientation through the section. Geophysical image logs and other logs were obtained for three wells located near the sites where specific models were developed for geomechanical simulations: Arches site in Boone County, Kentucky; Northern Appalachian Basin site in Chautauqua County, New York; and E-Central Appalachian Basin site in Tuscarawas County, Ohio. For these three wells, 9,700 feet of image logs were processed and interpreted to provide a systematic review of the distribution within each well of natural fractures, wellbore breakouts, faults, and drilling induced fractures. There were many borehole breakouts and drilling-induced tensile fractures but few natural fractures. Concentrated fractures were present at the Rome-basal sandstone and basal sandstone-Precambrian contacts at the Arches and East-Central Appalachian Basin sites. Geophysical logs were utilized to develop local-scale geologic models by determining geomechanical and petrophysical parameters within the geologic formations. These data were ported to coupled fluid-flow and reservoir geomechanics multi-phase CO 2 injection simulations. The models were developed to emphasize the geomechanical layers within the CO 2 storage zones and caprocks. A series of simulations were completed for each site to evaluate whether commercial-scale CO 2 could be safely injected into each site, given site-specific geologic and geomechanical controls. This involved analyzing the simulation results for the integrity of the caprock, intermediate, and reservoir zones, as well quantifying the areal uplift at the surface. Simulation results were also examined to ensure that the stress-stress perturbations were isolated within the subsurface, and that there was only limited upward migration of the CO 2 . Simulations showed capacity to inject more than 10 million metric tons of CO 2 in a single well at the Arches and East Central Appalachian Basin sites without excessive geomechanical risks. Low-permeability rock layers at the Northern Appalachian Basin study area well resulted in very low CO 2 injection capacity. Fracture models developed for the sites suggests that the sites have sparse fracture network in the deeper Cambrian rocks. However, there were indicators in image logs of a moderate fracture matrix in the Rose Run Sandstone at the Northern Appalachian Basin site. Dual permeability fracture matrix simulations suggest the much higher injection rates may be feasible in the fractured interval. Guidance was developed for geomechanical site characterization in the areas of geophysical logging, rock core testing, well testing, and site monitoring. The guidance demonstrates that there is a suitable array of options for addressing geomechanical issues at CO 2 storage sites. Finally, a review of Marcellus and Utica-Point Pleasant shale gas wells and CO 2 storage intervals indicates that these items are vertically separated, except for the Oriskany sandstone and Marcellus wells in southwest Pennsylvania and northern West Virginia. Together, project results present a more realistic portrayal of geomechanical risk factors related to CO 2 storage for existing and future coal-fired power plants in Ohio.

58 GEOSCIENCES↗

Optimal Realization of Yang–Baxter Gate on Quantum Computers

Quantum computers provide a promising method to study the dynamics of many-body systems beyond classical simulation. On the other hand, the analytical methods developed and results obtained from the integrable systems provide deep insights on the many-body system. Quantum simulation of the integrable system not only provides a valid benchmark for quantum computers but is also the first step in studying integrable-breaking systems. The building block for the simulation of an integrable system is the Yang–Baxter gate. It is vital to know how to optimally realize the Yang–Baxter gates on quantum computers. Based on the geometric picture of the Yang–Baxter gates, the optimal realizations of two types of Yang–Baxter gates with a minimal number of controlled NOT (CNOT) or gates are presented. It is also shown how to systematically realize the Yang–Baxter gates via the pulse control. The different realizations on IBM quantum computers are tested and compared. It is found that the pulse realizations of the Yang–Baxter gates always have a higher gate fidelity compared to the optimal CNOT or realizations. On the basis of the above optimal realizations, the simulation of the Yang–Baxter equation on quantum computers is demonstrated. Finally, these results provide a guideline and standard for further experimental studies based on the Yang–Baxter gate.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

GPU-accelerated dislocation dynamics using subcycling time-integration

Discrete dislocation dynamics (DDD) simulations are being increasingly employed to investigate metal plasticity at the mesoscale. However, in spite of its ability to access beyond the length and time limits of atomistic methods, the DDD model is still limited by its high computational cost, with ranges of achievable strains too low and strain rates too high by several orders of magnitude compared with typical experimental conditions. By combining the efficiency of the recently developed subcycling time-integrator with the highly parallel architecture of graphics processing unit (GPU) devices, we developed a DDD model that provides significant acceleration compared to existing implementations. Our GPU-accelerated implementation enables large-scale DDD simulations that can reach relevant levels of strain using a moderate amount of computational resources.

36 MATERIALS SCIENCE↗

The GENGA code: gravitational encounters in N-body simulations with GPU acceleration

We describe an open source GPU implementation of a hybrid symplectic N-body integrator, GENGA (Gravitational ENcounters with Gpu Acceleration), designed to integrate planet and planetesimal dynamics in the late stage of planet formation and stability analyses of planetary systems. GENGA uses a hybrid symplectic integrator to handle close encounters with very good energy conservation, which is essential in long-term planetary system integration. We extended the second-order hybrid integration scheme to higher orders. The GENGA code supports three simulation modes: integration of up to 2048 massive bodies, integration with up to a million test particles, or parallel integration of a large number of individual planetary systems. We compare the results of GENGA to Mercury and pkdgrav2 in terms of energy conservation and performance and find that the energy conservation of GENGA is comparable to Mercury and around two orders of magnitude better than pkdgrav2. GENGA runs up to 30 times faster than Mercury and up to 8 times faster than pkdgrav2. GENGA is written in CUDA C and runs on all NVIDIA GPUs with a computing capability of at least 2.0.

79 ASTRONOMY AND ASTROPHYSICS↗

Magnetic Designs of New First Target Beamline Magnets for the ORNL SNS Upgrade

The Spallation Neutron Source (SNS) at ORNL currently is being upgraded from 1.0 GeV to 1.3 GeV. Several water-cooled magnets should be upgraded to transport 30% of higher beam energy. New chicane, injection/extraction septum, and Lambertson magnets were designed. Designing the magnets was a challenging task because the new magnets required good combined integrated field quality and needed to occupy the old magnets space but with about 20% greater integrated magnetic field. Additional strong requirements applied to the magnets fringe field do not disturb the circulating beam. The special field profiles had to be provided in foil areas between magnets. The analysis described here was based on OPERA3D simulations. A special technique was used for analyzing the integrated field harmonics. Initially, the particle track was simulated, and integrated field components were calculated along this track for the reference radius, which were used for the harmonics analysis. In addition, 3D field maps were provided for beam optics simulations. The final beamline analysis confirmed good beam transmission and low losses.

3d field maps↗

Subsurface H 2 Storage: A Williston Basin Commercial-Scale Resource Study

Poster for the 2024 NETL Resource Sustainability Project Review Meeting, Pittsburgh, Pennsylvania, April 2-4, 2024. This poster presents a commercial‑scale assessment of subsurface hydrogen storage potential in the North Dakota portion of the Williston Basin, evaluating saline formations, depleted oil and gas reservoirs, and salt formations. The study integrates laboratory characterization, reservoir simulation, and basinwide analysis to assess storage capacity, injectivity, recovery, and risks related to geochemical, microbial, and wellbore interactions. Results support the feasibility of large‑volume, secure hydrogen storage and provide a framework to guide future hydrogen commercialization and infrastructure development.

08 HYDROGEN↗