Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time dependent behavior”

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 757 records · Page 42

Parsed synthesis of pyocyanin via co-culture enables context-dependent intercellular redox communication

Abstract Background Microbial co-cultures and consortia are of interest in cell-based molecular production and even as “smart” therapeutics in that one can take advantage of division of labor and specialization to expand both the range of available functions and mechanisms for control. The development of tools that enable coordination and modulation of consortia will be crucial for future application of multi-population cultures. In particular, these systems would benefit from an expanded toolset that enables orthogonal inter-strain communication. Results We created a co-culture for the synthesis of a redox-active phenazine signaling molecule, pyocyanin (PYO), by dividing its synthesis into the generation of its intermediate, phenazine carboxylic acid (PCA) from the first strain, followed by consumption of PCA and generation of PYO in a second strain. Interestingly, both PCA and PYO can be used to actuate gene expression in cells engineered with the soxRS oxidative stress regulon, although importantly this signaling activity was found to depend on growth media. That is, like other signaling motifs in bacterial systems, the signaling activity is context dependent. We then used this co-culture’s phenazine signals in a tri-culture to modulate gene expression and production of three model products: quorum sensing molecule autoinducer-1 and two fluorescent marker proteins, eGFP and DsRed. We also showed how these redox-based signals could be intermingled with other quorum-sensing (QS) signals which are more commonly used in synthetic biology, to control complex behaviors. To provide control over product synthesis in the tri-cultures, we also showed how a QS-induced growth control module could guide metabolic flux in one population and at the same time guide overall tri-culture function. Specifically, we showed that phenazine signal recognition, enabled through the oxidative stress response regulon soxRS, was dependent on media composition such that signal propagation within our parsed synthetic system could guide different desired outcomes based on the prevailing environment. In doing so, we expanded the range of signaling molecules available for coordination and the modes by which they can be utilized to influence overall function of a multi-population culture. Conclusions Our results show that redox-based signaling can be intermingled with other quorum sensing signaling in ways that enable user-defined control of microbial consortia yielding various outcomes defined by culture medium. Further, we demonstrated the utility of our previously designed growth control module in influencing signal propagation and metabolic activity is unimpeded by orthogonal redox-based signaling. By exploring novel multi-modal strategies for guiding communication and consortia outcome, the concepts introduced here may prove to be useful for coordination of multiple populations within complex microbial systems.

Chun, Kayla↗

The dynamics of dense and dilute clusters of drops evaporating in large, coherent vortices

The behavior of evaporating clusters of drops embedded into large, coherent vortices is described using a formulation which is valid for both dense and dilute clusters. Drops and gas interact both dynamically and thermodynamically. Dynamic coupling occurs through a force on the drops due to drag resulting from a slip velocity between the two phases. The net interaction force on the gas with drops is due to a source thrust from evaporation plus drag on each drop. The drag coefficient accounts for blowing from the drop surface. Thermodynamic coupling is a result of drop heating and evaporation. Limitations due to drop proximity on heating and evaporation are taken into account. The vortical motion of the drops in the cluster results in the formation of a core region devoid of drops at the center of the vortex, and a shell region containing the drops and surrounding the inner core. Results are presented showing the dependence of the evaporation time, the final to initial volume ratio and the final to initial shell thickness ratio upon the initial air/fuel mass ratio and as a function of the initial tangential velocities, upon the initial Stokes number, initial drop radius and initial outer cluster radius. Differences in behavior between and control parameters of dense and dilute clusters are pointed out by these new results. It is found that for dense clusters the final to initial volume ratio and final to initial shell thickness scale with the initial Stokes number, a new result which must be validated experimentally.

Bellan, J.↗

Navigating Team Dynamics: Automated Detection of Micro-Behaviors Between Team Members Through Longitudinal Interaction Data

The success in future long term space exploration missions will depend on the cooperation, coordination, and mutual understanding among the crew members. Micro-behaviors are momentary, subtle linguistic and paralinguistic indicators of thinking and feeling toward another member of the team (Cortina et al., 2001; Smith & Griffiths, 2022) that can significantly impact team dynamics and influence the overall team performance (Paromita & Chaspari, 2024). Due to their interactive nature, micro-behaviors have a sender (i.e., the team member expressing the micro-behavior) and a target (the team member impacted by the micro-behavior). Detection of these behaviors can assist in avoiding possible conflict among crew members and promoting the overall team success. Our prior research focused on an initial proof of concept of machine learning (ML) models and natural language processing (NLP) techniques that were used for automatically detect micro-behaviors among crew members of the US National Aeronautics and Space Administration’s (NASA) Human Exploration Research Analog (HERA) Campaigns 4 and 5 missions (Paromita et al., 2023). Results underscored the importance of incorporating contextual information in the ML models in the form of sentiment analysis, type of task, and dyadic interaction among team members. Here, we expand the scope of our prior work in two ways. First, we assess ML/NLP methods on new behavioral annotations coded using an adapted version of Smith & Griffins (2022) theoretical framework in terms of Violation (i.e., presence of valenced behavior, uplifting/positive or discouraging/negative), Intensity (i.e., force of behavior in terms of how uplifting or discouraging is the behavior), and Intent (i.e., motive of the behavior in terms of whether it was deliberate or unintentional). Second, we expand the design of the ML model to preserve information about the role of each team member within the occurrence of the micro-behavior (in contrast to the previous model that only considered the sender and the target without determining the team member role). This allows to consider all team members' contributions in the conversation and model long-term dependencies in the dialogue. Our experiments for this study are conducted on data from 5 teams of the NASA HERA C4 (NASA grant NNX16AQ48G (PI: Bell)). Conversations were extracted from the 1.5 hour Team Interaction Battery (TIB) task that occurred 5 times in-mission per crew. This resulted in a total of 13,058 conversational turns (i.e., 17.8% uplifting, 3.3% discouraging, 75.76% neutral, 3.14% nulls). Our findings with the revised behavioral coding and ML/NLP models indicate a 43.66% macro F1-score (i.e., 38.29% precision (P), 50.8% recall (R)) for a dialog state-tracking model that includes information from the sender only, and a 40.9% F1-score (i.e., 38.7% P, 43.36% R) for the same model that includes information from both the sender and the target of the micro-behavior. These are significantly higher compared to simple random forest models that classify behaviors strictly based on speech content and do not consider iterative team dynamics, achieving a 36.07% F1-score (i.e., 39.04% R, 33.53% P). Our findings demonstrate potential ways to leverage large conversational datasets to better capture complex team dynamics. We will discuss future directions including proposed models that can incorporate additional mission days and tasks beyond the TIB for objectively quantifying team behavior at high temporal resolution in space exploration missions.

Projna Paromita↗

Adoption of AI in the Utility T&D Sector: Use Cases, Consequence, Assessment and Benefits

Digital transformation and utilization of artificial intelligence (AI) in the electric grid are fundamentally changing the industry’s approach to common problems and enabling a broader paradigm shift in grid planning and operations. The change in approach is circularly both enabling and driving modernization, with load growth and reliable management of data center and AI infrastructure shifting away from planning approaches with relatively predictable behaviors and toward a mix of consumer and industrial choices that surpass human cognitive abilities to process. This movement has potential to condition humans to not understand the system on which the AI depends, while requiring it for development of the necessary infrastructure. Approaches which would address most likely grid conditions and events, such as faults, aging of equipment, and weather, now must also account for large loads which shift not based upon weather or time of day, but the computational load. Quantifying computational load is independent of the traditional grid forecasting variables, where a data center’s aggregate load is determined by user and AI system behavior and decoupled from normal grid planning and operations. AI is both the cause and solution for these challenges, with new grid planning tools integrating massive amounts of decisions into frameworks.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Progress Towards a Predictive Eagle Behavior and Risk Modeling Framework: Overview and Recent Validation Efforts

This presentation summarizes progress to date of the U.S. Department of Energy project, "Development of a computational framework for modeling golden eagles (Aquila chrysaetos) near wind farms," which focuses on stochastic behavioral modeling of soaring raptors across landscape, facility, and turbine spatiotemporal scales. This publicly available, open-source modeling framework includes behavioral models based on three different underlying principles: energy minimization at landscape scale, behavioral heuristics at landscape-facility scale, and data-driven behaviors at the facility-micro-scale. We will briefly overview the key advancements in the behavioral modeling state of the art, which leverages multiple high-resolution telemetry data sources combined with high-fidelity atmospheric flow modeling insights. We then present preliminary results from a validation study in Altamont, California. This new study involves a novel application of the Stochastic Soaring Raptor Simulator (SSRS), in a new geographic locale, to understand facility scale eagle movement patterns over time scales representative of a wind project's lifetime. For this desktop analysis (that does not depend on any high-performance computing resources), SSRS simultaneously considers a variety of wind conditions and eagle approach vectors toward a project site of interest. This work demonstrates the integration of publicly available landscape-scale atmospheric datasets, our recently improved engineering updraft models (see presentation from Thedin et al.), and our energy minimization behavioral models within the SSRS framework. While we only present results from a single behavioral model, the integration of these three modeling components forms the foundation for our more sophisticated behavioral models (see presentations from Brandes et al., Sandhu et al.) that are under active development. Results are presented in the form of presence maps, which may be applied to estimate risk to wildlife, augment ground survey data, inform wind-plant operations, or incorporated into wind-plant designs.

agent-based modeling↗

Analysis of time-dependent strain heterogeneity in shales using X-ray microscopy and digital volume correlation

In this paper, we report the results of a 4-dimensional (3 spatial dimensions plus time) in-situ creep experiment on a milli-core shale plug from Barnett formation in Texas. An outstanding challenge in shale mechanics is connecting observable characteristics to predictable mechanical and hydraulic behavior. We employ novel imaging and analysis techniques to statistically link sample mineralogy and structure to observed deformation. In particular, we focus on mechanisms of porosity reduction during transient response to step changes in stress, analogous to the aseismic deformation that occurs after an episode of fracturing. Our main goal in this study is to take the first steps to visualize and capture the small-scale deformation mechanisms involved in creep of shales, using imaging techniques. To do so, we conducted a uniaxial creep experiment inside an Xradia Versa 520 micro-CT system for 12 h at room temperature under a constant uniaxial stress of 25 MPa. We employed a Digital Volume Correlation (DVC) method to measure micro-scale deforspmation of the sample by comparing the resulting volumetric images with a high-resolution scan of the sample obtained prior to the creep experiment. The results indicate that the time-dependent deformation in the compliant clay plus kerogen matrix constituents facilitates rigid grain rotation and compaction of intergranular porosity at the boundary between rigid grains and matrix.

58 GEOSCIENCES↗

Macroscopic instructions vs microscopic operations in quantum circuits

In many experiments on microscopic quantum systems, it is implicitly assumed that when a macroscopic procedure or “instruction” is repeated many times – perhaps in different contexts – each application results in the same microscopic quantum operation. But in practice, the microscopic effect of a single macroscopic instruction can easily depend on its context. If undetected, this can lead to unexpected behavior and unreliable results. Here, we design and analyze several tests to detect context-dependence. They are based on invariants of matrix products, and while they can be as data intensive as quantum process tomography, they do not require tomographic reconstruction, and are insensitive to imperfect knowledge about the experiments. We also construct a measure of how unitary (reversible) an operation is, and show how to estimate the volume of physical states accessible by a quantum operation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ensembles of Metastable States Govern Heterogeneous Catalysis on Dynamic Interfaces

Heterogeneous catalysis is at the heart of chemical industry. Being able to tune and design efficient catalysts for processes of interest is of an utmost importance, and for this, the molecular-level understanding of heterogeneous catalysts is the first step, and indeed a prime focus of the modern catalysis research. For a long time, a single most thermodynamically stable structure of the catalytic interface attained in reaction conditions had been envisioned as the reactive phase. However, some catalytic interfaces continue to undergo structural dynamics in the steady state, triggered by high temperatures, pressures, and binding and changing reagents. Among particularly dynamic interfaces are such widely-used catalysts as crystalline and amorphous surfaced supporting (sub-)nano metallic clusters. Recently, it became clear that this dynamic fluxionality causes the supported clusters to populate many distinct structural and stoichiometric states in catalytic conditions. Hence, the catalytic interface should be viewed as an evolving statistical ensemble of many (not one) structures. As such, every member in the ensemble contributes to the properties of the catalyst differently, and in proportion to its probability to be populated. This new notion flips the established paradigm and calls for new theory, modeling approaches, operando measurements, and updated design strategies. The statistical ensemble nature of surface-supported sub-nano cluster catalysts can be exemplified by oxide-supported and adsorbate-covered Pt, Pd, Cu, CuPd clusters, catalytic toward oxidative and non-oxidative dehydrogenation. They have access to a variety of 3D and quasi-2D shapes. The compositions of their thermal ensembles are dependent on the cluster size, leading to size-specific catalytic activities and the famous “every atom counts” phenomenon. The support and adsorbates affect catalyst structures, and state of the reacting species causes the ensemble to change in every reaction intermediate. The most stable member of the ensemble dominates the thermodynamic properties of the corresponding intermediate, whereas the kinetics can be determined by more active but less populated metastable catalyst states, and that suggests that many earlier studies might have overlooked the actual active sites. Both effects depend on the relative time-scales of catalyst restructuring and reaction dynamics. The catalyst may routinely operate off-equilibrium. Ensemble phenomena lead to surprising exceptions from established rules of catalysis, such as scaling relations, and Arrhenius behavior. Catalyst deactivation is also an ensemble property, and its extent of mitigation can be predicted through the new paradigm. These findings were enabled by advances in theory, such as global optimization and subsequent utilization of multiple local minima, and pathways sampling, as well as operando catalyst characterization. The fact that the per-site and per-species resolution is needed for the description and predicting of catalyst properties gives theory the central role in catalysis research, as most experiments provide ensemble-average information and cannot detect the crucial minority species that may be responsible for the catalytic activity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ultrafast Nanoimaging of Carrier Funneling in Composition-Graded Semiconductor Nanowires

Recent advances in bandgap engineering of low-dimensional semiconductors have enabled high-efficiency carrier transport in miniaturized electronic and optoelectronic devices. The physical properties and functionalities of these materials are governed by complex carrier dynamics coupled with multiple transport mechanisms in tailored band structures. Here, we report ultrafast nanoimaging of carrier funneling and recombination in composition-grade CdSxSe1-x nanowires using pump-probe near-field nanoscopy. Leveraging the high resolution of this technique in both space and time, we resolve nanoscale local carrier dynamics along composition-graded nanowires, revealing the local variation of composition-dependent carrier mobilities and lifetimes that significantly differ from their uniform composition counterparts. Furthermore, we demonstrate a length-dependent behavior wherein shorter nanowires exhibit enhanced funneling effects, accelerating carrier transport by up to 33%. Our findings provide direct visualization of nanoscale carrier transport while supporting an effective approach for investigating complex carrier interactions in inhomogeneous semiconductor nanostructures, with implications for optimizing next-generation optoelectronic devices.

Yang, Rundi↗

Deciphering helix assembly in the heliconical nematic phase via tender resonant X-ray scattering

Being a link between a uniaxial nematic and a chiral nematic, the twist bend nematic (N TB ) has been an intriguing topic for research over the last decade as a key to understanding chirality generation. Accurate description and understanding of resonance effects in helical structures provide crucial knowledge on the phase behavior beyond positional ordering. We examined the manifestation of resonance effects in N TB via tender resonant X-ray scattering (TReXS) at the sulfur K-edge. For the first time we demonstrate quantitatively that the energy dependence of the scattering peak in the N TB phase follows the energy dependence of the complex refractive indices measured by X-ray absorption. Furthermore, taking advantage of molecular energy calculation, we decipher the distinct helical pitch variation trend for symmetric/asymmetric dimers from the perspective of hybridization and reveal the strong effect of substituting oxygen for sulfur atoms on the helical pitch and the thermal helix extension.

36 MATERIALS SCIENCE↗

Elucidating plasma dynamics in Hasegawa–Wakatani turbulence by information geometry

In this work, the impact of adiabatic electrons on drift-wave turbulence, modeled by the Hasegawa–Wakatani equations, is studied using information length. Information length is a novel theoretical method for measuring distances between statistical states represented by different probability distribution functions (PDFs) along the path of a system and represents the total number of statistically different states that a system evolves through in time. Specifically, the time-dependent PDFs of turbulent fluctuations for a given adiabatic index $A$ are computed. The changes in fluctuation statistics are then quantified in time by using information length. The numerical results provide time traces exhibiting intermittent plasma dynamics, and such behavior is identified by a rapid change in the information length. The effects of $A$ are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

In Situ High-Temperature Ultrafast Electron Diffraction through Integrated Furnace and MEMS Platforms

Temperature fundamentally governs phase stability, defect evolution, and transport behavior in materials. Despite its central role, direct measurements of structural evolution at elevated temperatures on ultrafast timescales have remained limited. Here, we report the design, integration, and validation of 2 complementary in situ heating platforms that substantially extend the thermal operating range of ultrafast electron diffraction (UED). A compact furnace-type heating stage enables stable diffraction measurements from room temperature to 800 K with ±0.1 K stability under ultrahigh vacuum, achieved through multi-sensor feedback control, dual air-cooling channels, and a thermally isolated motion stage. In parallel, a microelectromechanical system (MEMS)-based heating platform provides rapid thermal response and access to extreme temperatures ≥1,373 K with ±0.1 K stability over hundreds-micrometer regions while supporting simultaneous electrical biasing for electrothermal coupling studies. Absolute temperature calibration is established using diffraction-based thermometry via aluminum lattice expansion and independently validated through in situ melting of bismuth thin films. UED measurements further reveal pronounced temperature-dependent nonequilibrium lattice dynamics in bismuth, including modifications to electron–phonon coupling and Debye–Waller behavior, as well as enhanced ultrafast diffuse scattering in aluminum at elevated temperatures. Together, these developments establish a practical framework for quantitative, time-resolved studies of temperature-driven kinetics and nonequilibrium structural dynamics under extreme thermal environments.

Bai, Qianqian [Chinese Academy of Sciences (CAS), ↗

Transition to hard turbulence in thermal convection at infinite Prandtl number

Direct numerical simulations of two-dimensional high Rayleigh (Ra) number, base-heated thermal convection in large aspect-ratio boxes are presented for infinite Prandtl number fluids, as applied to the earth's mantle. A transition is characaterized in the flow structures in the neighborhood of Ra between 10 to the 7th and 10 to the 8th. These high Ra flows consist of large-scale cells with strong intermittent, boundary-layer instabilities. For Ra exceeding 10 to the 7th it is found that the heat-transfer mechanism changes from one characterized by mushroom-like plumes to one consisting of disconnected ascending instabilities, which do not carry with them all the thermal anomaly from the bottom boundary layer. Plume-plume collisions become much more prominent in high Ra situations and have a tendency of generating a pulse-like behavior in the fixed plume. This type of instability represents a distinct mode of heat transfer in the hard turbulent regime. Predictions of this model can be used to address certain issues concerning the mode of time-dependent convection in the earth's mantle.

Hansen, Ulrich↗

A new finite element formulation for computational fluid dynamics. IX - Fourier analysis of space-time Galerkin/least-squares algorithms

A Fourier stability and accuracy analysis of the space-time Galerkin/least-squares method as applied to a time-dependent advective-diffusive model problem is presented. Two time discretizations are studied: a constant-in-time approximation and a linear-in-time approximation. Corresponding space-time predictor multi-corrector algorithms are also derived and studied. The behavior of the space-time algorithms is compared to algorithms based on semidiscrete formulations.

Shakib, Farzin↗

Using Indirect Turbulence Measurements for Real-Time Parameter Estimation in Turbulent Air

The use of indirect turbulence measurements for real-time estimation of parameters in a linear longitudinal dynamics model in atmospheric turbulence was studied. It is shown that measuring the atmospheric turbulence makes it possible to treat the turbulence as a measured explanatory variable in the parameter estimation problem. Commercial off-the-shelf sensors were researched and evaluated, then compared to air data booms. Sources of colored noise in the explanatory variables resulting from typical turbulence measurement techniques were identified and studied. A major source of colored noise in the explanatory variables was identified as frequency dependent upwash and time delay. The resulting upwash and time delay corrections were analyzed and compared to previous time shift dynamic modeling research. Simulation data as well as flight test data in atmospheric turbulence were used to verify the time delay behavior. Recommendations are given for follow on flight research and instrumentation.

Martos, Borja↗

Cloud Computing Methods for Near Rectilinear Halo Orbit Trajectory Design

Complicated mission design problems require innovative computational solutions. As spacecraft depart from a proposed Gateway in a Near Rectilinear Halo Orbit (NRHO), recontact analysis is required to avoid risk of collision and ensure safe operations. Escape dynamics from NRHOs are governed by multiple gravitational bodies, yielding a trajectory design space that is exhaustively large. This paper summarizes the recontact analysis for departure from the NRHO and describes how the Deep Space Trajectory Explorer (DSTE) trajectory design software incorporates high performance cloud computing to compute and visualize the orbit design space. Recent focus on exploration missions to cislunar space has kindled accelerated interest in multibody orbit solutions. Trajectory analysis in the presence of multiple gravity fields is complex, and innovative computational tools are needed to simplify complicated design spaces, to generate large quantities of data quickly, and to visualize the output for user accessibility. The Gateway mission is a prime example. The Gateway1 is proposed as a human outpost in deep space. The current baseline orbit for the Gateway is a Near Rectilinear Halo Orbit (NRHO) near the Moon.2 The NRHO exists in a regime that experiences the gravitational effects of the Earth and the Moon simultaneously, complicating orbit analysis. The mission design process benefits greatly from updated computational tools for multibody missions like the Gateway. As an example, consider the problem of assessing the risk of collision in an NRHO. As a staging location to missions to the lunar surface and beyond the Earth-Moon system, the Gateway will experience spacecraft and other objects regularly arriving and departing. Departing objects potentially include spent logistics modules, visiting crew vehicles, debris objects, wastewater particles, and cubesats. Each departure is governed by the dynamics of the Gateway orbit and the surrounding dynamical environment. Over time, any unmaintained object in such an orbit eventually departs due to the small instabilities associated with the NRHOs. A separation maneuver speeds the departure from the NRHO, but the effects of the maneuver on the spacecraft behavior depend on the location, magnitude, and direction of the burn. Escape dynamics from the NRHO with regard to these maneuver options open up an enormous potential trajectory design space where subtle changes in input can produce dramatically large changes in the results. Any departing object must avoid recontacting the Gateway as it leaves the lunar vicinity, and a recontact analysis thus involves a significant number of computations and extensive output data. To explore the dynamics of this extensive design space, the Deep Space Trajectory Explorer3 (DSTE) trajectory design software incorporates new High Performance Computing (HPC) services and novel interactive visualizations. This paper details the HPC and cloud infrastructure techniques that are implemented in the DSTE, applying the new capabilities to analysis of recontact risk with the Gateway in NRHO. NEAR RECTILINEAR HALO ORBITS The Gateway is planned to fly in a lunar NRHO as its baseline orbit. The NRHO families of orbits are subsets of the larger halo families, which originate from planar orbits near the L1 and L2 libration points; the Earth-Moon L2 halo family appears in Figure 1. Each halo orbit is perfectly periodic in the Circular Restricted 3-Body Problem (CR3BP) and becomes a quasi-periodic orbit in a higher fidelity ephemeris force model. The NRHOs are defined as those members of the halo family with bounded stability properties;2 they pass near the Moon at perilune and are nearly polar. Families exist with apolunes located both above the lunar north pole and above the lunar south pole; the Gateway is planned to reside in a southern L2 NRHO in a 9:2 resonance with the lunar synodic period. The 9:2 NRHO is characterized by a period of about 6.5 days, a perilune radius of about 3,500 km, and an apolune radius of about 71,000 km; it is strongly affected by the gravity of both the Earth and the Moon simultaneously. This NRHO offers extended communications with assets on the south pole of the Moon,4 as well as low-cost orbit maintenance and attitude control,5 favorable eclipse avoidance properties,6 and inexpensive transfers from Earth and to other destinations.5,7 The NRHO portion of the southern L2 halo family is highlighted in black in Figure 1, and the 9:2 NRHO appears in blue.

Phillips, Sean M.↗

Unveiling the Potential of MeshGraphNets for Predicting Subsurface Evolution in Carbon Storage Projects

This is the conference paper accompanying an oral presentation “Unveiling the Potential of MeshGraphNets for Predicting Subsurface Evolution in Carbon Storage Projects” at the 17th International Conference on Greenhouse Gas Control Technologies GHGT-17 held in Calgary, Canada, October 20-24 , 2024. Carbon capture and storage (CCS) technology is critical for mitigating climate change but requires effective subsurface reservoir management to ensure safe containment of injected CO2. Accurate predictions of reservoir pressure and saturation are essential for assessing long-term CCS performance. Traditional numerical simulations, while effective, are computationally intensive, time-consuming, and constrained by data discretization. Previous work has shown the effectiveness of MeshGraphNets (MGN), a graph-based machine learning framework, as an innovative alternative for predicting reservoir behavior. MGN leverages graph neural networks (GNNs) and mesh representations to model complex geological formations, offering superior adaptability across different discretizations and reservoir configurations. Classic MGN implementations utilize an autoregressive technique to predict future behavior based on current predictions, but this technique is hampered by error accumulation over time. To enhance the model accuracy in time-series predictions, this study implemented a multi-step rollout strategy that integrates autoregressive predictions during training to stabilize prediction of saturation over time. Using the Illinois Basin – Decatur Project (IBDP) dataset, comprising 100 simulations of CO2 injection, pressure, and saturation changes, the framework demonstrated its ability to learn spatial dependencies and temporal dynamics. With inputs including permeabilities, porosities, and injection rates, MGN accurately predicted CO2 plume evolution over time, even with limited training data. Moreover, the addition of a multi-step rollout procedure during training improved the ability of MGN to predict stably over time by ~15%. This research positions MGN, enhanced with multi-step rollout capabilities, as a robust and efficient tool for CCS applications. It advances the field by enabling precise, computationally efficient predictions of reservoir behavior, providing a foundation for the broader adoption of machine learning frameworks in CCS and other geoscience domains.

Holcomb, Paul↗

Gravity Effects in Diffusive Coarsening of Bubble Lattices: von Neumann's Law

von Neumann modelled the evolution of two-dimensional soap froths as a purely diffusive phenomenon; the area growth of a given cell was found to depend only on the geometry of the bubble lattice. In the model, hexagons are stable, pentagons shrink and heptagons grow. The simplest equivalent to the area growth law is / approximately t(sub beta). The result depends on assuming (1) an incompressible gas; (2) bubble walls which meet at 120 deg and (3) constant wall thickness and curvature. Each assumption is borne out in experiments except the last one: bubble wall thickness between connecting cells varies in unit gravity because of gravity drainage. The bottom part of the soap membrane is thickened, the top part is thinned, such that gas diffusion across the membrane shows a complex dependence on gravity. As a result, experimental tests of von Neumann's law have been influenced by effects of gravity; fluid behavior along cell borders can give non-uniform wall thicknesses and thus alter the effective area and gas diffusion rates between adjacent bubbles. For area plotted as a function of time, Glazier (J.A. Glazier, S.P. Gross, and I. Stavans, Phys. Rev. A. 36, 306 (1987); J. Stavans, J.A, Glazier, Phys. Rev. Lett. 62, 1318 (1989).) suggest that in some cases their failure to observe von Neumann's predicted growth exponent ((sup beta)theor(sup =1; beta)exp(sup =0.70 + 0.10)) may have been the result of such "fluid drainage onto the lower glass plate". Additional experiments which varied plate spacing gave different beta exponents in a fashion consistent with this suggestion. During preliminary long duration experiments (approximately 100 h) aboard Spacelab-J, a low-gravity test of froth coarsening has examined (1) power law scaling of von Neumann's law (beta values) in the appropriate diffusive limits; (2) new bubble lattice dynamics such as greater fluid wetting behavior on froth membranes in low gravity; and (3) explicit relations for the gravity dependence of the second moment (or disorder parameter) governing the geometric spread in cell-sidedness around the mean of perfect hexagonal filling. By reducing the gravity-induced distortion in lattice wall thickness, the diffusion-limited regime of bubble coarsening becomes available for performing critical tests of network dynamics.

Noever, David A.↗