SEARCH · Engineering Papers

Results for “Injection data analysis”

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.

SourceResearch connectionFrom yearThrough year

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

At least 703 records · Page 39

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

PyroCbs from Australia Fires and its Impact Using Satellite Observations from CrIS and TROPOMI and Reanalysis Data

Pyrocumulonimbus (pyroCb) clouds are thunder clouds created by intense heat from the Earth’s surface. They are formed similarly to cumulonimbus clouds, but the intense heat that results in the vigorous updraft comes from fire, either large wildfires or volcanic eruptions. Australia’s unprecedented fire disasters at the end of 2019 to early 2020 emitted huge amounts of carbon monoxide (CO) and fire aerosol particles to the atmosphere, particularly during the pyroCb outbreak that occurred in southeast Australia between 29 December 2019 and 4 January 2020. It was estimated that at least 18 pyroCbs were generated during this episode, and some of them injected ice, smoke, and biomass burning gases above the local tropopause. An unprecedented abundance of H2O and CO in the stratosphere, and the displacement of background ozone (O3) and N2O from rapid ascent of air from the troposphere and lower stratosphere were found from satellite observations. Some other studies also found that the fire emissions and their long-range transport resulted in stratospheric aerosol, temperature, and O3 anomalies after the 2020 Australian bushfires and altered the Antarctic ozone and vortex, posing great impact to local air quality and climate change. Further study on the atmospheric thermodynamic status of atmosphere associated with these pyroCbs, and the change of the cloud properties and trace gases during this unprecedented Australia fires will be made using a new single Field of View (SFOV) Sounder Atmospheric Products (SiFSAP). SiFSAP was developed by NASA using the Cross-track Infrared Sounder (CrIS) and Advanced Technology Microwave Sounder (ATMS) onboard SNPP and JPSS-1, and will soon be available to the public at NASA DAAC. Since this product has a spatial resolution of 15 km at nadir, which is better than most global weather and climate models and other current operational sounding products, a process-oriented analysis of the dynamic transport of CO and fire plumes during this unprecedented fire disasters will be made in this study. Based on a Principal Component Radiative Transfer Model (PCRTM) and an optimized estimation retrieval algorithm, a simultaneously retrieval is made using the whole spectral information measured by CrIS, and the derived SiFSAP include temperature, water vapor, trace gases (such as O3, CO2, CO, CH4 and N2O), cloud properties and surface properties. Use of ATMS together with CrIS allows SiFSAP to get accurate retrieval products under thick pyroCb conditions, and an algorithm to detect pyroCb based on the hyperspectral infrared sounder spectrum from CrIS will be developed and verified. In addition to SiFSAP sounding products, other products like CO, O3, NO2 from TROPOMI, O3 from OMPS will be used for retrospective analysis. The wind fields from the NASA’s Modern-Era Retrospective Analysis for Research and Applications Version-2 (MERRA-2) and ERA5 will be used to characterize the transport, and the SiFSAP temperature and water vapor profiles within and around pyroCbs will be compared with MERRA-2 and ERA5 products.

Xiaozhen (Shawn) Xiong↗

DOE OSTI2022

Time Sequence Machine Learning-Based Data Intrusion Detection for Smart Voltage Source Converter-Enabled Power Grid

Smart inverters of distributed energy resources can enable cloud computing, condition monitoring, result visualization, remote control, and peer-to-peer energy trading in advanced power systems. However, the advent of data injection attacks in the communication architecture can alter measurement characteristics of power grids and have devastating consequences. In this article, we propose a time sequence machine learning-based anomaly detection methodology for detecting cyber intrusion into control signal setpoints and dc voltage signal measurement bias of the voltage source converter (VSC) in wind generators. We first investigated the effects of four types of denial of service, tampering signal, and stealthy-type data intrusion attacks on smart VSCs and overall wind farms. We then proposed a novel time sequence machine learning-based intrusion detection framework that can be implemented to detect different cyberattacks in the VSCs. The performance of the proposed framework has been compared with that of autoencoder and clustering-based intrusion detection framework. The proposed framework was validated by using the IEEE 39 bus power system in the presence of four wind farms in different locations. Using several metrics for intrusion detection performance, we validated the effectiveness of the proposed framework.

42 ENGINEERING↗

NASA NTRS2005

Constraining Oceanic Sources and Sinks of CO2 by Inverse Modeling: A First Step Toward Assimilation of In Situ and Remote Data

The primary goal of this collaborative project, lead by Nicholas Gruber at UCLA, is to better quantify the air-sea fluxes of pre-industrial and anthropogenic carbon dioxide CO2 for the global ocean with a newly developed inverse method. The method involves using ocean general circulation models (OGCMs) to characterize the subsurface distribution of tracers injected for different regions of the ocean surface. By linearly combining a series of these basis (or Greens) functions, one should be able to recreate any arbitrary ocean tracer field. In our case, we are using the global database of ocean dissolved inorganic carbon (DIC) from the WOCE/JGOFS global CO2 survey from the 1990s. The observed DIC field has been partitioned into pre-industrial and anthropogenic components using the empirical C* technique. We then solve for an optimal combination of the regional basis functions that best match the two sets observations in a least-squares sense with singular value decomposition. This in turn provides solutions for the regional air-sea CO2 fluxes (pre-industrial and anthropogenic) that are independent of more traditional estimates derived from surface water disequilibria in the partial pressure of CO2 and wind speed gas exchange parameterizations. We are apply the method to a suite of nine different OGCMs to test the robustness of the derived surface fluxes fields to differences in ocean circulation. Together with a collegue Keith Lindsay at NCAR, S. Doney's role in this project was to create the appropriate steady-state and temporally evolving basis functions from the NCAR CSM-1 ocean circulation model. Because of a slow start across the project in defining a common set of experimental protocols, we requested and were given a one year, no cost extension. Both the steady-state and temporally evolving basis functions were delivered to the UCLA group for further analysis. Additional information is included in the original extended abstract.

Doney, Scott↗

DOE OSTI2025

A comparative analysis of numerical approaches for the description of gas flow in clay-based repository systems: From a laboratory to a large-scale gas injection test

There is nowadays a consensus among many countries that geological disposal is a favourable solution for the long-term management. Although different host formations and different barrier systems are under consideration around the world, clay-based materials form an important component for waste isolation in most national programmes. Hence, a good comprehension of the effect of gas flow on the hydro-mechanical behaviour of clay-based soils is essential, both at laboratory and field scale. Task B under the international cooperative project DECOVALEX-2023 has recently shown that, after some enhancement, models can be employed to reproduce laboratory scale tests, even with different sample geometries 37 . However, further work is required to understand whether they can be applied to simulate a large-scale experiment. Up-scaling of models for the advective transport of gas through clay-based low permeable material presents a number of problems related to the difficulty in obtaining consistent hydrogeological parameters and constitutive relationships at both laboratory and field scale. Based on a unique dataset from a large-scale gas injection test (Lasgit) performed at the Äspö Hard Rock Laboratory (Sweden), Task B within DECOVALEX-2023 has explored the refinement of these numerical strategies applied to the simulation of gas flow. Work performed within the task reveals that codes do not need to be substantially modified from the laboratory models to reproduce full-scale tests: indeed, model parameters calibrated and validated at laboratory scale have been applied to predict field scale gas flow at Lasgit, including peak gas pressure and injected cumulative gas volume. By means of (1) the introduction of interfaces between blocks to reflect the experimental configuration and the (2) adjustment of some parameters (e.g., higher permeability), the updated models are able to represent most of the key features observed in the experimental data, even at a large scale.

Tamayo-Mas, E↗

DOE OSTI2024

Geologic Characterization of the South Georgia Rift Basin for Source Proximal CO2 Storage

The project Geologic Characterization of the South Georgia Rift Basin for Source Proximal CO2 Storage is one of 9 site characterization projects that were implemented as part of ARRA (American Recovery and Reinvestment Act). Data from this project was used to improve resolution of data in NATCARB in the area of study. Data related to this study has already been incorporated in NATCARB Atlas. The South Carolina Research Foundation and partners evaluated the feasibility of CCS in the Jurassic/ Triassic (J / TR) saline formations of the buried Mesozoic South Georgia Rift (SGR) Basin that extends from South Carolina into Georgia. The J / TR sequence, based on preliminary assessment of limited geologic and geophysical data, appears to have both the appropriate areal extent and multiple horizons to permanently and safely store CO2 The presence of several igneous rock layers within the sequence may potentially provide adequate seals to prevent upward CO2 migration into the Coastal Plain aquifer systems. Approximately 81 kilometers of 2-D seismic reflection data were collected by Bay Geophysical, Inc. to explore a portion of the SGR located in southern Georgia. The 81 kilometers were divided into two lines approximately 40.5 kilometers each, with Line 1 intersecting Georgia well GGS 3457. Line 2 intersects Line 1 at the southern portion of Line 1 to maximize the extent of coverage away from GGS-3457 (a deep well drilled in the 1980s for oil and gas exploration). This well had a set of usable logs, including gamma and neutron logs that provided promising results related to CO2 storage. Results showed sandstone with porosity values greater than 10 percent and a thickness of 120 meters. The design of the seismic shot was to extrapolate information away from the well and to better define the extent of the SGR and the necessary reservoir and caprock for a successful CO2 injection. A numerical simulation model of CO2 Injection and migration was developed based on the geology log for the GGS-3457 well. The simulation model was used to investigate the feasibility of injecting 30 million metric tons of CO2 into SGR J / TA sediments and integrity of the diabase layers as seals to prevent CO2 migration.

2-D seismic↗

DOE OSTI2020

Scale and Regionality of Nonelectric Markets for U.S. Nuclear Light Water Reactors

This study assesses existing and potential industries that could conceivably be directly coupled to existing nuclear reactors. The goal is to identify the scale, location, and accessibility of the candidate industrial-product markets, as well as process feedstocks that are available near the plants to establish new industries. For example, CO 2 as a feedstock can be combined with H 2 to produce formic acid (FA), transportation fuels, and lubricants. These new plants can be entirely supported with the heat and electricity provided by a nearby NPP. The potential demand for nonelectric industrial products was assessed by documenting current and possible growth of nonelectricity product markets considered. This assessment used DOE- and industry-supported tools, data, and projections to capture regional industrial market opportunities. Electricity-capacity markets that reward large and reliable generators, such as NPPs, were considered because the electricity market will likely continue to be an important revenue source to NPPs. The key is to balance the needs of energy customers so as to optimize revenue for the affiliated energy customers or partners. In most cases, flexible plant energy delivery and power generation for the grid will require either energy storage or a stock of intermediate products to sustain the industrial customers when the NPP dispatches electricity to the grid. A diverse mix of temperate regions with operating NPPs around the U.S.—representing a variety of operating markets, local generation mix, and seasonal climates—were chosen for this market study. Both current and future market opportunities for candidate industrial-product markets surrounding these NPPs were studied. Figure 2 illustrates the regions chosen for this study. The success of developing nonelectric industrial-product markets as alternative revenue-generating sources for LWRs depends, not only on demand from growing existing markets, such as petroleum refining and NH 3 production, but also on the development of new markets such as light-duty (LD) and heavy-duty (HD) hydrogen FCEVs, synfuels, chemical production, biofuels, metal refining, injection of hydrogen into NG pipelines for gas power-generating units, FA, polymers, and close-coupled industrial heat applications, all of which can significantly increase demand relative to current levels while decarbonizing energy sectors. This study also presents a sample analysis of the economics of hydrogen production in an area of Minnesota, considering the capital and operating costs of a hydrogen plant as well as the local market demand for hydrogen. It includes some assumptions on electricity-grid pricing, showing how hydrogen could be integrated with an NPP and be competitive with the incumbent hydrogen-production process, steam methane reforming (SMR). The objectives of this study include: Provide U.S. NPP operators a robust sampling of the market demand location, scale, and accessibility (including storage and transportation) of the wide variety of industrial-product choices that can be produced using nuclear thermal energy and electricity proximate to a subset of U.S. NPPs to inform the industry of the potential opportunity; show examples and trade-off analyses of how U.S. LWR operators can access these markets, including storage and transportation of industrial products to their intended markets; and present a general analysis example for one industrial product (hydrogen) in one region (Minnesota area), including production, storage, and transportation, to show how nuclear-hybrid integrated energy systems (IESs) could access local markets and improve the profitability of an NPP.

03 NATURAL GAS↗

DOE OSTI2022

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↗

NASA NTRS2004

Characterization of Composite Fan Case Resins

The majority of commercial turbine engines that power today s aircraft use a large fan driven by the engine core to generate thrust which dramatically increases the engine s efficiency. However, if one of these fan blades fails during flight, it becomes high energy shrapnel, potentially impacting the engine or puncturing the aircraft itself and thus risking the lives of passengers. To solve this problem, the fan case must be capable of containing a fan blade should it break off during flight. Currently, all commercial fan cases are made of either just a thick metal barrier or a thinner metal wall surrounded by Kevlar-an ultra strong fiber that elastically catches the blade. My summer 2004 project was to characterize the resins for a composite fan case that will be lighter and more efficient than the current metal. The composite fan case is created by braiding carbon fibers and injecting a polymer resin into the braid. The resin holds the fibers together, so at first using the strongest polymer appears to logically lead to the strongest fan case. Unfortunately, the stronger polymers are too viscous when melted. This makes the manufacturing process more difficult because the polymer does not flow as freely through the braid, and the final product is less dense. With all of this in mind, it is important to remember that the strength of the polymer is still imperative; the case must still contain blades with high impact energy. The research identified which polymer had the right balance of properties, including ease of fabrication, toughness, and ability to transfer the load to the carbon fibers. Resin deformation was studied to better understand the composite response during high speed impact. My role in this research was the testing of polymers using dynamic mechanical analysis and tensile, compression, and torsion testing. Dynamic mechanical analysis examines the response of materials under cyclic loading. Two techniques were used for dynamic mechanical analysis. The ARES Instrument analyzed the material through torsion. The second machine, TA Instruments apparatus, applied a bending force to the specimen. These experiments were used to explore the effects of temperature and strain rate on the stiffness and strength of the resins. The two different types of loading allowed us to verify our results. An axial-torsional load frame, manufactured by MTS Systems, Inc., was used to conduct the tensile, compression, and torsional testing. These tests were used to determine the stress-strain curves for the resins. The elastic and plastic deformation data was provided to another team member for characterization of high fidelity material property predictions. This information was useful in having a better understanding of the polymers so that the fan cases could be as sturdy as possible. Deformation studies are the foundation for the computational modeling that provides the structural design of a composite engine case as well as detailed analysis of the blade impact event.

Dvoracek, Charlene M.↗

DOE OSTI2024

NORTH DAKOTA CARBONSAFE PHASE III: SITE CHARACTERIZATION AND PERMITTING OF GEOLOGIC STORAGE OF CARBON DIOXIDE

The Energy & Environmental Research Center (EERC), in partnership with Minnkota Power Cooperative Inc. (Minnkota), SLB, and Computer Modelling Group Ltd. (CMG), supported wide-scale deployment of carbon capture and storage (CCS) as part of the U.S. Department of Energy (DOE) National Energy Technology Laboratory Carbon Storage Assurance Facility Enterprise (CarbonSAFE) Initiative Phase III. This phase included the acquisition, analysis, and development of information to fully characterize two storage complexes to demonstrate viable storage resources for commercial volumes of CO2 (defined by DOE as a minimum of 50 million tonnes [MMt] of CO2 within a 30-year period) (National Energy Technology Laboratory, 2024). Phase III also involved the preparation, submission, and approval of North Dakota underground injection control (UIC) Class VI storage facility permits (SFPs)—required precursors to applications for Class VI injection well permits. The presumed viability of commercial-scale CCS, situated adjacent to Minnkota’s Milton R. Young Station (MRYS), is validated by Minnkota’s continued pursuit of Project Tundra—an initiative to build the world’s largest lignite-based CCS project in central North Dakota (www.projecttundrand.com). Project Tundra comprises two scopes of work, Tundra Capture (installation of postcombustion CO2 capture at MRYS) and Tundra SGS (secure geologic storage). The efforts of North Dakota CarbonSAFE Phase III, Site Characterization and Permitting, supported Tundra SGS. Extensive site-specific characterization activities included a successful multimeasurement geophysical approach and drilling a stratigraphic test well (J-ROC 1, subsequently renamed Liberty-1) adjacent to MRYS. Core collection and analyses, downhole testing and fluid sampling, and geophysical logging were performed on J-ROC 1 and on a nearby stratigraphic test well (known as J-LOC 1), which was drilled, cored, and tested under a complementary project funded by the North Dakota Lignite Research Program. The injection tests performed on J-LOC 1 positively impacted the CarbonSAFE project, resulting in fewer proposed injection wells and significant construction, operations, and monitoring cost savings. The characterization data collected and analyses performed were integrated into geologic models, and successive numerical simulations were run to determine CO2 plume extent and subsurface pressure buildup associated with the planned CO2 injection rate of nearly 4 MMt per year. The latter doubles the CarbonSAFE Initiative goal with an estimated 100 MMt of CO2 stored in 20 years. Application of the U.S. Environmental Protection Agency’s (EPA’s) method for estimating the Class VI Rule area of review (AOR) to the overpressurized Broom Creek Formation inspired an alternative method of calculation, called risk-based AOR delineation. This peer-reviewed method was applied for the first time during the storage facility-permitting process. The two SFP applications submitted in 2021 successfully resulted in North Dakota Industrial Commission (NDIC) orders in 2022 authorizing the creation of the storage facility areas and amalgamation of pore space as well as establishing financial responsibility requirements. After approval of the SFPs, Minnkota filed in 2022 applications for permits to reenter the J-ROC 1 well and to drill two new wells—all with the intended purpose to become Class VI injection wells. To establish eligibility under the Internal Revenue Code for Section 45Q tax incentives, a monitoring, reporting, and verification (MRV) plan was prepared and submitted by Minnkota to EPA in 2022, resulting in the first such plan approved in North Dakota. Also in 2022, under the National Environmental Policy Act (NEPA), Minnkota prepared and submitted an environmental information volume (EIV) describing the proposed CCS project and associated potential environmental impacts. Based on the EIV, DOE determined that the proposed construction project required an environmental assessment, and Minnkota submitted the first draft in 2023 and a revised draft in 2024. Both submissions were followed by a public comment period. Subsequently, DOE issued a finding of no significant impact (FONSI) on September 13, 2024. A successful outreach program, strongly based in the production, presentation, and dissemination of informational material, fostered an environment to aid stakeholders in making informed decisions regarding the planned project. Opportunities for public input were provided at various steps along the way, including at county planning and zoning meetings, before and during the SFP administrative hearing, and during environmental assessment public comment periods. In addition, land/pore space owners and mineral owners had various points of contact, including granting access rights, securing pore space leasing, and mineral owner notifications. Based upon the successful storage facility permitting issued by NDIC, approval of the MRV plan by EPA, and receipt of a FONSI under the NEPA, Minnkota is continuing its pursuit of Project Tundra. In December 2023, the Office of Clean Energy Demonstrations under its Carbon Capture Demonstrations Projects Program announced funding for the capture system (Office of Clean Energy Demonstrations, 2023) and a proposal for CarbonSAFE Phase IV: Construction funding was submitted in March 2024 for the storage project. A go/no-go decision to proceed with construction and operations in the Broom Creek Formation is anticipated in 2024. References National Energy Technology Laboratory, CarbonSafe Initiative, https://netl.doe.gov/carbon-management/carbon-storage/carbonsafe (accessed August 2024). Office of Clean Energy Demonstrations, 2023, OCED selects three projects in CA, ND, and TX to reduce harmful carbon pollution, create new economic opportunities, and advance carbon reducing technologies, December, www.energy.gov/oced/articles/oced-selects-three-projects-ca-nd-and-tx-reduce-harmful-carbon-pollution-create-new (accessed August 2024).

Peck, Wesley↗

DOE OSTI2021

Skin factor and potential formation damage from chemical and mechanical processes in a naturally fractured carbonate aquifer with implications to CO 2 sequestration

Here, in this study, we investigate formation damage due to acidization and water injection tests into the naturally fractured carbonate Middle Duperow Formation at Kevin Dome, Montana, potentially diminishing the chance of a future successful Geological Carbon Sequestration (GCS) project. Multiple well-test analytical models, correlated with core description and lithology data, are used to determine flow behavior and communication between the water injection interval and surrounding formations. An improved three-dimensional (3D) geologic model with dual-continuum matrix and fracture properties is constructed based on most recent seismic, core, and water sample measurements. Brine injection is simulated to verify the interpretation from the analytical models, followed by CO 2 injection simulation. Geochemical calculations are performed to understand the in-situ processes that led to formation damage. Our findings suggest: (1) there are two possible scenarios that could lead to a positive total effective skin factor and permeability decline: partial penetration and formation damage; (2) analytical models indicate a positive total skin factor, contradicting results of a previous study suggesting that the well was mildly stimulated; (3) numerical simulation supports the formation damage hypothesis by matching the pressure buildup observed during the latter two brine injection tests; (4) several mechanical and chemical processes may have occurred during injection to clog the matrix/fracture system: anhydrite fines migration and/or calcite precipitation. We then make preventative suggestions for future GCS projects into carbonate reservoirs and remediation recommendations for GCS operation at the Kevin Dome.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1996

Transneuronal pathways to the vestibulocerebellum

The alpha-herpes virus (pseudorabies, PRV) was used to observe central nervous system (CNS) pathways associated with the vestibulocerebellar system. Retrograde transneuronal migration of alpha-herpes virions from specific lobules of the gerbil and rat vestibulo-cerebellar cortex was detected immunohistochemically. Using a time series analysis, progression of infection along polyneuronal cerebellar afferent pathways was examined. Pressure injections of > 20 nanoliters of a 10(8) plaque forming units (pfu) per ml solution of virus were sufficient to initiate an infectious locus which resulted in labeled neurons in the inferior olivary subnuclei, vestibular nuclei, and their afferent cell groups in a progressive temporal fashion and in growing complexity with increasing incubation time. We show that climbing fibers and some other cerebellar afferent fibers transported the virus retrogradely from the cerebellum within 24 hours. One to three days after cerebellar infection discrete cell groups were labeled and appropriate laterality within crossed projections was preserved. Subsequent nuclei labeled with PRV after infection of the flocculus/paraflocculus, or nodulus/uvula, included the following: vestibular (e.g., z) and inferior olivary nuclei (e.g., dorsal cap), accessory oculomotor (e.g., Darkschewitsch n.) and accessory optic related nuclei, (e.g., the nucleus of the optic tract, and the medial terminal nucleus); noradrenergic, raphe, and reticular cell groups (e.g., locus coeruleus, dorsal raphe, raphe pontis, and the lateral reticular tract); other vestibulocerebellum sites, the periaqueductal gray, substantia nigra, hippocampus, thalamus and hypothalamus, amygdala, septal nuclei, and the frontal, cingulate, entorhinal, perirhinal, and insular cortices. However, there were differences in the resulting labeling between infection in either region. Double-labeling experiments revealed that vestibular efferent neurons are located adjacent to, but are not included among, flocculus-projecting supragenual neurons. PRV transport from the vestibular labyrinth and cervical muscles also resulted in CNS infections. Virus propagation in situ provides specific connectivity information based on the functional transport across synapses. The findings support and extend anatomical data regarding vestibulo-olivo-cerebellar pathways.

NASA Discipline Neuroscience↗

NASA NTRS1995

Two cis elements collaborate to spatially repress transcription from a sea urchin promoter

The expression pattern of many territory-specific genes in metazoan embryos is maintained by an active process of negative spatial regulation. However, the mechanism of this strategy of gene regulation is not well understood in any system. Here we show that reporter constructs containing regulatory sequence for the SM30-alpha gene of Stronglyocentrotus purpuratus are expressed in a pattern congruent with that of the endogenous SM30 gene(s), largely as a result of active transcriptional repression in cell lineages in which the gene is not normally expressed. Chloramphenicol acetyl transferase assays of deletion constructs from the 2600-bp upstream region showed that repressive elements were present in the region from -1628 to -300. In situ hybridization analysis showed that the spatial fidelity of expression was severely compromised when the region from -1628 to -300 was deleted. Two highly repetitive sequence motifs, (G/A/C)CCCCT and (T/C)(T/A/C)CTTTT(T/A/C), are present in the -1628 to -300 region. Representatives of these elements were analyzed by gel mobility shift experiments and were found to interact specifically with protein in crude nuclear extracts. When oligonucleotides containing either sequence element were co-injected with a correctly regulated reporter as potential competitors, the reporter was expressed in inappropriate cells. When composite oligonucleotides, containing both sequence elements, were fused to a misregulated reporter, the expression of the reporter in inappropriate cells was suppressed. Comparison of composite oligonucleotides with oligonucleotides containing single constituent elements show that both sequence elements are required for effective spatial regulation. Thus, both individual elements are required, but only a composite element containing both elements is sufficient to function as a tissue-specific repressive element.

Non-NASA Center↗

NASA NTRSDate not supplied

A Comparison of Stratospheric Smoke Events Seen by SAGE II (1984-2005) and SAGE III (2017-2021)

The ongoing Stratospheric Aerosol and Gas Experiment aboard the International Space Station (SAGE III/ISS) has mostly occurred during a period of relatively low aerosol loading. At the same time, it has almost observed a number of small to moderate injections of volcanic material in the stratosphere from the Aoba/Ambae eruption of 2018 to that of Karymsky in 2022. In addition, two of the largest known smoke intrusions into the stratosphere have occurred during the SAGE III/ISS mission (the British Columbia fire of August 2017 and the Australian Wildfire Event of 2019/2020). At least two other much smaller smoke events have also been observed in the SAGE III data (California Creek in September 2020 and McKay Creek in June 2021). Particularly with the two large events, it is tempting to speculate that these events are increasing in intensity and frequency possibly due to ongoing climate change. Realistically, change is difficult to assess since instrument data sets that can be used for quantifying historical events can be quite different in what and where they provide measurements. This complicates comparing one event to another potentially separated by many years. An exception to this issues are the measurements by SAGE III/ISS and those of its predecessor SAGE II (1984-2005). Both use the same measurement strategy, measure aerosol extinction coefficient at 1020 nm, and have similar spatial/temporal coverage. Both instruments are known to observe smoke events in the stratosphere. However, both have the same low observation density characteristic of solar occultation such that near-global coverage takes almost a month and characterizing effectively instant events like smoke and volcanic events take a matter of weeks depending on the vagaries of where the observations are being made at a critical juncture unlike instruments with high temporal/spatial resolution like limb-scatter instruments such as OMPS and OSIRIS. The goal of this work will be to use the same analysis approach to both the SAGE II and SAGE III data and compare both the frequency and intensity of smoke events observed by these instruments. While we do not expect to be able to make conclusions about change between the periods, the analysis is a first step in interpreting changes in the contribution of smoke to the stratospheric aerosol budget between the time frames encompassed by these instruments.

L. W. Thomason↗

DOE OSTI2023

The Cost and Benefit of Enhancing Cybersecurity for Hybrid AC/DC Grids

As critical interfaces of AC grids and DC grids inside a hybrid AC/DC grid, the voltage-sourced-converter (VSC) has been demonstrated to be vulnerable to false data injection (FDI) cyber-attacks. As a result, the cyber-attack-induced AC grid frequency deviations and DC grid voltage deviations threaten the secure operation. Here, to enhance cybersecurity in a not only feasible but also cost-effective manner, this paper proposes a cost-benefit-based cyber-defense strategy for a hybrid AC/DC grid. First, this paper establishes a spatial-temporal dual cyber-attack evaluation model, in which the cyber-attack-induced frequency and voltage deviations are modelled in both a spatially and temporally dual manner. Then, the proposed cost-benefit-based cyber-defense strategy is modelled as a VSC commitment problem to achieve the trade-off between maximizing the cyber-defense benefits and minimizing the cyber-defense costs. The VSC commitment problem is then mathematically convexified into a mixed-integer second-order cone programming (MISOCP) problem, which could be efficiently solved in an event-triggered manner against unfolding cyber-attack events. Simulation results on a test hybrid AC/DC grid verified the feasibility and the cost-effectiveness of the proposed cyber-defense strategy.

cost-benefit analysis↗

NASA NTRS1992

Dynamic interactions between hypersonic vehicle aerodynamics and propulsion system performance

Described here is the development of a flexible simulation model for scramjet hypersonic propulsion systems. The primary goal is determination of sensitivity of the thrust vector and other system parameters to angle of attack changes of the vehicle. Such information is crucial in design and analysis of control system performance for hypersonic vehicles. The code is also intended to be a key element in carrying out dynamic interaction studies involving the influence of vehicle vibrations on propulsion system/control system coupling and flight stability. Simple models are employed to represent the various processes comprising the propulsion system. A method of characteristics (MOC) approach is used to solve the forebody and external nozzle flow fields. This results in a very fast computational algorithm capable of carrying out the vast number of simulation computations needed in guidance, stability, and control studies. The three-dimensional fore- and aft body (nozzle) geometry is characterized by the centerline profiles as represented by a series of coordinate points and body cross-section curvature. The engine module geometry is represented by an adjustable vertical grid to accommodate variations of the field parameters throughout the inlet and combustor. The scramjet inlet is modeled as a two-dimensional supersonic flow containing adjustable sidewall wedges and multiple fuel injection struts. The inlet geometry including the sidewall wedge angles, the number of injection struts, their sweepback relative to the vehicle reference line, and strut cross-section are user selectable. Combustion is currently represented by a Rayleigh line calculation including corrections for variable gas properties; improved models are being developed for this important element of the propulsion flow field. The program generates (1) variation of thrust magnitude and direction with angle of attack, (2) pitching moment and line of action of the thrust vector, (3) pressure and temperature distributions throughout the system, and (4) performance parameters such as thrust coefficient, specific impulse, mass flow rates, and equivalence ratio. Preliminary results are in good agreement with available performance data for systems resembling the NASP vehicle configuration.

Flandro, G. A.↗

DOE OSTI2023

Development and evaluation of a skeletal mechanism for EHN additized gasoline mixtures in large Eddy simulations of HCCI combustion

Advanced Low Temperature Combustion modes, such as the Sandia proposed Additive-Mixing Fuel Injection (AMFI), can unlock significant potential to boost fuel conversion efficiency and ultimately improve the energy conversion of internal combustion engines. This is a novel improved combustion process that is enabled by supplying small (<5%) variable amounts of autoignition improver to the fuel to enhance the engine operation and control. Common, diesel-fuel ignition-quality enhancing additive, 2-ethylexyl nitrate (EHN), is doped into gasoline to enable Sandia LTGC + AMFI combustion. This manuscript focuses on the development of a reduced sub-mechanism for EHN chemical kinetics at engine relevant conditions that is implemented into a skeletal mechanism for chemical kinetic studies of gasoline surrogate fuels. The mechanism validation utilized zero-dimensional numerical simulations and comparison to shock tube ignition-delay data of pure and EHN-doped n-heptane. Additional validation is presented with Homogeneous Charge Compression-Ignition (HCCI) engine data of pure and EHN-doped research-grade E10 gasoline. Then, the mechanism was deployed in a 3-D computational fluid dynamics (CFD) using Large Eddy Simulations (LES) to model the HCCI engine experiments of 0.4% vol EHN additized E10 gasoline at several equivalence ratios. Simulations showed a very good performance of the mechanism, and the model accurately reproduced (a) the ignition point, (b) combustion phasing, (c) combustion duration, and (d) the peak of the heat release rates of the engine experiments. The results show that EHN promotes Low-Temperature Heat Release, ultimately driving the gasoline to autoignite at thermodynamic conditions where the fuel would not otherwise ignite. Overall, this work demonstrates a viable reduced chemical-kinetic mechanism for EHN and shows that it can be combined with a skeletal gasoline mechanism for CFD-LES analysis of well-mixed LTGC that matches well with experimental results. The CFD-LES analysis also shows the spatial distribution of EHN-fuel interactions that control the autoignition throughout the combustion chamber.

Guleria, Gaurav↗

NASA NTRS1996

An Acoustic Emission and Acousto-Ultrasonic Analysis of Impact Damaged Composite Pressure Vessels

The research presented herein summarizes the development of acoustic emission (AE) and acousto-ultrasonic (AU) techniques for the nondestructive evaluation of filament wound composite pressure vessels. Vessels fabricated from both graphite and kevlar fibers with an epoxy matrix were examined prior to hydroburst using AU and during hydroburst using AE. A dead weight drop apparatus featuring both blunt and sharp impactor tips was utilized to produce a single known energy 'damage' level in each of the vessels so that the degree to which the effects of impact damage could be measured. The damage levels ranged from barely visible to obvious fiber breakage and delamination. Independent neural network burst pressure prediction models were developed from a sample of each fiber/resin material system. Here, the cumulative AE amplitude distribution data collected from low level proof test (25% of the expected burst for undamaged vessels) were used to measure the effects of the impact on the residual burst pressure of the vessels. The results of the AE/neural network model for the inert propellant filled graphite/epoxy vessels 'IM7/3501-6, IM7/977-2 and IM7/8553-45' demonstrated that burst pressures can be predicted from low level AE proof test data, yielding an average error of 5.0%. The trained network for the IM7/977-2 class vessels was also able to predict the expected burst pressure of taller vessels (three times longer hoop region length) constructed of the same material and using the same manufacturing technique, with an average error of 4.9%. To a lesser extent, the burst pressure prediction models could also measure the effects of impact damage to the kevlar/epoxy 'Kevlar 49/ DPL862' vessels. Here though, due to the higher attenuation of the material, an insufficient amount of AE amplitude information was collected to generate robust network models. Although, the worst case trial errors were less than 6%, when additional blind predictions were attempted, errors as high as 50% were produced. An acousto-ultrasonic robotic evaluation system (AURES) was developed for mapping the effects of damage on filament wound pressure vessels prior to hydroproof testing. The AURES injects a single broadband ultrasonic pulse into each vessel at preprogrammed positions and records the effects of the interaction of that pulse on the material volume with a broadband receiver. A stress wave factor in the form of the energy associated with the 750 to 1000 kHz and 1000 to 1250 kHz frequency bands were used to map the potential failure sites for each vessel. The energy map associated with the graphite/epoxy vessels was found to decrease in the region of the impact damage. The kevlar vessels showed the opposite trend, with the energy values increasing around the damage/failure sites.

Walker, James L.↗

DOE OSTI2024

A deep learning-based workflow for fast prediction of 3D state variables in geological carbon storage: A dimension reduction approach

Deep learning (DL) models are extensively used as surrogate models for high-fidelity simulations of multiphase fluid flow in porous media at large scales, enabling fast forecasts of the spatial–temporal evolution of three-dimensional (3D) state variables in geological carbon storage (GCS). However, training these models in high-dimensional space remains computationally demanding and prone to overfitting because of limited training data. This paper presents a novel workflow to address these challenges by integrating dimension reduction (DR) methods. Here, the proposed workflow employed pre-trained DR models to extract the latent variables of geological models and state variables and utilized the multi-layer perceptron (MLP) for constructing mapping functions between the input and output variables in latent spaces. Subsequently, the pre-trained reconstruction models converted the MLP-predicted latent state variables to their original high-dimensional form. Furthermore, we proposed a novel strategy for the DR and reconstruction of 3D saturation fields to account for the unique data characteristics of sparsity, nonuniformity, and discontinuity. The proposed strategy applied PCA and inverse PCA for 2D average saturation fields and developed a DL-based 3D reconstruction model, leveraging three 2D average saturation fields as input to produce a 3D saturation field as output. The pre-training of DR and reconstruction models and training of MLP models were conducted on 84 Gulf of Mexico (GoM) simulations and evaluated on 12 testing simulations. Each simulation contained 720 monthly time steps, with the first 360 months as the injection period and the rest as the post-injection period. The proposed workflow, incorporating DR and DL models, accurately predicts the normalized 3D pressure fields, achieving mean square error (MSE) of 2.92 × 10 -7 compared to the ground truth obtained from a full-physics simulator. Furthermore, the proposed strategy outperformed PCA and convolutional autoencoder (CAE) models on 3D saturation fields, resulting in minor workflow prediction errors with an MSE of 2.93 × 10 -5 . The results suggest the proposed workflow provides sufficient predictive fidelity across temporal and spatial scales, and enables a speedup of 160 times compared to the full-physics simulator, facilitating improved decision-making and risk assessment for large-scale GCS management in real-time scenarios.

3D reconstruction model↗