Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dynamic mechanical 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.

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

At least 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Statistical Analysis of Inter-Area Oscillations in the U.S. Eastern Interconnection: A 2017-2023 Perspective

Recent advancements and the accumulation of high-resolution, long-term phasor measurement unit (PMU) data have provided detailed insights into inter-area oscillations in power grids. This study conducts a comprehensive statistical analysis of inter-area oscillations within the United States Eastern Interconnection from 2017 to 2023. Utilizing data captured by the advanced wide-area Frequency Monitoring Network (FNET/GridEye), this investigation examines the occurrence patterns, dominant frequencies, damping ratios, and excitation mechanisms of these oscillations. Our analysis sheds light on the evolving statistical behaviors of inter-area oscillations, offering updated and critical information for grid operators and planners. The insights gained from this study can be instrumental in enhancing the operational resilience of the power network and guiding strategic developments in grid infrastructure to accommodate future challenges. Additionally, the study discusses emerging challenges associated with the modernization of the power grid, including increased renewable penetration, dynamic load variability, and cyber-physical vulnerabilities that complicate oscillation monitoring and control.

Inter-area oscillations↗
DOE OSTI2020

Persistent Protein Motions in a Rugged Energy Landscape Revealed by Normal Mode Ensemble Analysis

Proteins are allosteric machines that couple motions at distinct, often distant, sites to control biological function. Low-frequency structural vibrations are a mechanism of this long-distance connection and are often used computationally to predict correlations, but experimentally identifying the vibrations associated with specific motions has proved challenging. Spectroscopy is an ideal tool to explore these excitations, but measurements have been largely unable to identify important frequency bands. The result is at odds with some previous calculations and raises the question what methods could successfully characterize protein structural vibrations. Here we show the lack of spectral structure arises in part from the variations in protein structure as the protein samples the energy landscape. However, by averaging over the energy landscape as sampled using an aggregate 18.5 μs of all-atom molecular dynamics simulation of hen egg white lysozyme and normal-mode analyses, we find vibrations with large overlap with functional displacements are surprisingly concentrated in narrow frequency bands. These bands are not apparent in either the ensemble averaged vibrational density of states or isotropic absorption. However, in the case of the ensemble averaged anisotropic absorption, there is persistent spectral structure and overlap between this structure and the functional displacement frequency bands. We systematically lay out heuristics for calculating the spectra robustly, including the need for statistical sampling of the protein and inclusion of adequate water in the spectral calculation. The results show the congested spectrum of these complex molecules obscures important frequency bands associated with function and reveal a method to overcome this congestion by combining structurally sensitive spectroscopy with robust normal mode ensemble analysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Dynamics of Hydroxyl Anions Promotes Lithium Ion Conduction in Antiperovskite Li 2 OHCl

Li 2 OHCl is an exemplar of the antiperovskite family of ionic conductors, for which high ionic conductivities have been reported, but in which the atomic-level mechanism of ion migration is unclear. The stable phase is both crystallographically defective and disordered, having similar to 1/3 of the Li sites vacant, while the presence of the OH(- )anion introduces the possibility of rotational disorder that may be coupled to cation migration. In this work, complementary experimental and computational methods are applied to understand the relationship between the crystal chemistry and ionic conductivity in Li 2 OHCl , which undergoes an orthorhombic to cubic phase transition near 311 K (approximate to 38 degrees C) and coincides with the more than a factor of 10 change in ionic conductivity (from 1.2 x 10 -5 mS/cm at 37 degrees C to 1.4 x 10 -3 mS/cm at 39 degrees C). X ray and neutron experiments conducted over the temperature range 20-200 degrees C, including diffraction, quasi-elastic neutron scattering (QENS), the maximum entropy method (MEM) analysis, and ab initio molecular dynamics (AIMD) simulations, together show conclusively that the high lithium ion conductivity of cubic Li2OHCl is correlated to "paddlewheel" rotation of the dynamic OH - anion. The present results suggest that in antiperovskites and derivative structures a high cation vacancy concentration combined with the presence of disordered molecular anions can lead to high cation mobility.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Tetrahedral Tilting and Lithium‐Ion Transport in Halide Argyrodites Prepared by Rapid, Microwave‐Assisted Synthesis

Abstract This study demonstrates a rapid, dry, microwave‐assisted (MW) synthesis method that enables preparation of halide argyrodites ( , , ) in less than 20 min. The structures and ion transport properties of the resulting materials are compared with those synthesized by conventional solid‐state synthesis methods. The microwave‐assisted method leads to increased site disorder and elevated Arrhenius prefactors (), which lead to an order of magnitude improvement in the 30 ionic conductivity of MW‐. X‐ray pair distribution function analysis (XPDF) reveals significant rotational disorder of the units, which is impacted by the synthesis method, choice of halide, and presence of / site disorder. These rotational displacements are strongly correlated with ion transport, specifically and entropy of migration (). Overall, this study demonstrates a rapid synthesis route for preparing high‐quality halide argyrodite solid‐state electrolytes in less than 20 min, and further unravels atomistic insights into the interplay of structural disorder, rotational dynamics, and ion transport mechanisms.

36 MATERIALS SCIENCE↗
DOE OSTI2025

LeaPP: Learning Pathways to Polymorphs through Machine Learning Analysis of Atomic Trajectories

Understanding the mechanisms underlying crystal nucleation and growth is crucial for many technological applications. Due to the short length and time scales involved, crystal nucleation is often studied using molecular simulations. Most existing approaches to extract the nucleation mechanism from simulations focus on the analysis of static snapshots of the configurations, potentially overlooking subtle local fluctuations and the history of the particles involved in the formation of solid nuclei. Here, in this work, we propose a novel methodology called LeaPP that categorizes nucleation trajectories based on the temporal information of their constituent particles. We leverage the time evolution of the local environment of the crystallizing particles to encapsulate the relationship between the structure and dynamics and distinguish between different evolving particle paths. Identification of the distinct particle paths further enables characterizing the nucleation trajectories into different pathways. Collectively, LeaPP provides a more nuanced understanding of nucleation through an unsupervised approach with lesser dependence on traditional order parameters. Furthermore, the pathways identified by LeaPP are predictive of the resulting polymorph. We demonstrate LeaPP on three different systems─Lennard-Jones-like particles, Ni 3 Al, and water on surfaces. The general methodology underlying LeaPP─considering the time evolution of the building blocks─applies to a wide range of self-assembly problems.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Contributions of extreme and non-extreme precipitation to California precipitation seasonality changes under warming

Changes in seasonality of extreme precipitation have important implications for flood and fire hazards and water resources. Here we explore the contributions of extreme and non-extreme precipitation in the sharpened wet season in California under warming and the underlying mechanisms. Modeling evidence shows that increased extreme precipitation due to enhanced extreme intensity and more extreme days dominates the wetter winter, while decreased non-extreme precipitation due to fewer wet days induces the dryer spring and fall. Moisture budget analysis indicates that increased moisture dominates the extreme precipitation increase in winter, while weakened circulation offsets the moisture increase, resulting in no changes in spring and fall. The sharpened seasonal cycle of extreme precipitation via both dynamic and thermodynamic effects is consistent with the seasonality changes in atmospheric river days and is mainly attributed to the seasonality changes in the number of extreme days rather than the intensity during extreme days.

Dong, Lu↗
DOE OSTI2025

Performance Analysis and Simulaion of the Hydraulic Scram System in TREAT Reactor

The Transient Reactor Test Facility (TREAT) at Idaho National Laboratory (INL) serves a vital role in nuclear fuel safety research, enabling transient experiments that simulate reactivity excursions and accident scenarios. Central to these operations is the transient control rod drive system (TCRDS), which drives rapid motion of the transient control rods such that TREAT can simulate rapid power changes typical of reactor accidents. The reliability and performance of this system are critical for protecting both fuel specimens and reactor infrastructure. This study presents the initial phase of a two-year investigation into the dynamics and reliability of the TREAT hydraulic TCRDS. Conducted in collaboration with INL, the research employs a combined computational and experimental approach to analyze the system's response time, pressure transients, and potential failure modes. Emphasis is placed on understanding how fluid characteristics influence the TCRDS’s ability to achieve both rapid power changes and mechanical stability. The TRDS and the skid that powers it will be analyzed throughout this investigation. Computational modeling using computational fluid dynamics (CFD) will simulate the hydraulic response under varying conditions. In parallel, experimental testing planned at INL will validate these models and capture key performance metrics. This paper outlines the system design, analytical framework, and modeling strategies that form the foundation for later testing. Ultimately, this work aims to support improvements to the TCRDS’s design and reliability, contributing to the broader goal of enhancing nuclear fuel safety and sustaining TREAT’s mission as a premier nuclear fuel test facility.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2021

Understanding Hepatitis B Virus Dynamics and the Antiviral Effect of Interferon Alpha Treatment in Humanized Chimeric Mice

Whereas the mode of action of lamivudine (LAM) against hepatitis B virus (HBV) is well established, the inhibition mechanism(s) of interferon alpha (IFN-α) is less completely defined. To advance our understanding, we, for this study, mathematically modeled HBV kinetics during 14-day pegylated IFN-α-2a (pegIFN), LAM, or pegIFN-plus-LAM (pegIFN+LAM) treatment of 39 chronically HBV-infected humanized uPA/SCID chimeric mice. Serum HBV DNA and intracellular HBV DNA were measured frequently. We developed a multicompartmental mathematical model and simultaneously fit it to the serum and intracellular HBV DNA data. Unexpectedly, even in the absence of an adaptive immune response, a biphasic decline in serum HBV DNA and intracellular HBV DNA was observed in response to all treatments. Kinetic analysis and modeling indicate that the first phase represents inhibition of intracellular HBV DNA synthesis and secretion, which was similar under all treatments with an overall mean efficacy of 98%. In contrast, there were distinct differences in HBV decline during the second phase, which was accounted for in the model by a time-dependent inhibition of intracellular HBV DNA synthesis, with the steepest decline observed during pegIFN+LAM treatment (1.28/day) and the slowest (0.1/day) during pegIFN monotherapy. Reminiscent of observations in patients treated with pegIFN and/or LAM, a biphasic HBV decline was observed in treated humanized mice in the absence of an adaptive immune response. Interestingly, combination treatment did not increase the initial inhibition of HBV production but rather enhanced second-phase decline, providing insight into the dynamics of HBV treatment response and the mode of action of IFN-α against HBV.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Tailoring electrophotonic capabilities of atomically thin GeS through controlled organometallic intercalation

The unique structure of van der Waals materials facilitates a robust mechanism for precise control of physical properties. In this study, we present a comprehensive analysis based on the intercalation of organometallics to modulate the optoelectronic behavior of two-dimensional germanium sulfide (GeS). Advanced computational exploration reveals significant and tunable features in the intercalated material. Additionally, the weak chemical interactions between organometallics and GeS support the electric-field-mediated drift and charge–discharge processes in intercalants. Controlling the concentration of organometallics in this manner enables the dynamic emergence of novel characteristics post-intercalation. These include flatbands near the Fermi level, significant enhancement of carrier mobility, and a magnetic ground state that is atypical for pristine GeS. Our findings demonstrate that organometallic intercalation offers a powerful strategy for tailoring the optoelectronic and magnetic characteristics of GeS, paving the way for harnessing emerging features for applications in next-generation devices.

Physics↗
DOE OSTI2022

Quantum advantage for differential equation analysis

Quantum algorithms for differential equation solving, data processing, and machine learning potentially offer an exponential speedup over all known classical algorithms. However, there also exist obstacles to obtaining this potential speedup in useful problem instances. The essential obstacle for quantum differential equation solving is that outputting useful information may require difficult postprocessing, and the essential obstacle for quantum data processing and machine learning is that inputting the data is a difficult task just by itself. In this study, we demonstrate that, when combined, these difficulties solve one another. We show how the output of quantum differential equation solving can serve as the input for quantum data processing and machine learning, allowing dynamical analysis in terms of principal components, power spectra, and wavelet decompositions. To illustrate this, we consider continuous-time Markov processes on epidemiological and social networks. These quantum algorithms provide an exponential advantage over existing classical Monte Carlo methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Assessing individual movement, habitat use, and behavior of non-breeding marine birds in relation to prey availability in the US Atlantic

Resource availability is a key factor driving marine bird movements and distributions, but direct information on prey availability is difficult to obtain at relevant scales. We present novel methods for describing multi-scale trophic associations, combining movement analyses of marine birds with estimates of forage fish surface aggregations from digital aerial survey data and species occupancy from bottom trawl survey data. We analyzed satellite telemetry data from northern gannets Morus bassanus , red-throated loons Gavia stellata , and long-tailed ducks Clangula hyemalis in the US Atlantic during the non-breeding period. Using discrete-time hidden Markov models to distinguish area-restricted (i.e. putative foraging) from transit movements, we examined how environmental factors influence movement, and how forage fish species distributions and surface aggregations influence habitat use by gannets and loons that have greater dietary reliance. Our results suggest that chlorophyll a concentration significantly affected movement behavior across species, highlighting the importance of higher-productivity areas around estuaries during colder months when regional productivity is low. Though variable across species and seasons, spatial cross-correlation analysis revealed that herring species (Family Clupeidae), including Atlantic menhaden Brevoortia tyrannus , may be important resources; it also showed positive spatial correlations with forage fish aggregations. This suggests that prey patch dynamics and factors driving aggregation formation may be as important as species composition. However, spatial patterns were generally low (<0.3), suggesting a mismatch in spatiotemporal resolution, exemplifying the challenges in quantifying trophic relationships in marine systems. Disentangling predator-prey relationships is critical to understanding the mechanisms driving marine bird behavior in rapidly changing marine systems.

Environmental Sciences & Ecology↗
DOE OSTI2020

Changes in DNA Methylation in Response to 6-Benzylaminopurine Affect Allele-Specific Gene Expression in Populus Tomentosa

Cytokinins play important roles in the growth and development of plants. Physiological and photosynthetic characteristics are common indicators to measure the growth and development in plants. However, few reports have described the molecular mechanisms of physiological and photosynthetic changes in response to cytokinin, particularly in woody plants. DNA methylation is an essential epigenetic modification that dynamically regulates gene expression in response to the external environment. In this study, we examined genome-wide DNA methylation variation and transcriptional variation in poplar (Populus tomentosa) after short-term treatment with the synthetic cytokinin 6-benzylaminopurine (6-BA). We identified 460 significantly differentially methylated regions (DMRs) in response to 6-BA treatment. Transcriptome analysis showed that 339 protein-coding genes, 262 long non-coding RNAs (lncRNAs), and 15,793 24-nt small interfering RNAs (siRNAs) were differentially expressed under 6-BA treatment. Among these, 79% were differentially expressed between alleles in P. tomentosa, and 102,819 allele-specific expression (ASE) loci in 19,200 genes were detected showing differences in ASE levels after 6-BA treatment. Combined DNA methylation and gene expression analysis demonstrated that DNA methylation plays an important role in regulating allele-specific gene expression. To further investigate the relationship between these 6-BA-responsive genes and phenotypic variation, we performed SNP analysis of 460 6-BA-responsive DMRs via re-sequencing using a natural population of P. tomentosa, and we identified 206 SNPs that were significantly associated with growth and wood properties. Association analysis indicated that 53% of loci with allele-specific expression had primarily dominant effects on poplar traits. Our comprehensive analyses of P. tomentosa DNA methylation and the regulation of allele-specific gene expression suggest that DNA methylation is an important regulator of imbalanced expression between allelic loci.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Structure and thermodynamics of water adsorption in NU-1500-Cr

Abstract Metal-organic frameworks (MOFs) are a class of materials with diverse chemical and structural properties, and have been shown to effectively adsorb various types of guest molecules. The mechanism of water adsorption in NU-1500-Cr, a high-performance atmospheric water harvesting MOF, is investigated using a combination of molecular dynamics simulations and infrared spectroscopy. Calculations of thermodynamic and dynamical properties of water as a function of relative humidity allow for following the adsorption process from the initial hydration stage to complete filling of the MOF pores. Initial hydration begins at the water molecules that saturate the open Cr 3+ sites of the framework, which is then followed by the formation of water chains that extend along the channels connecting the hexagonal pores of the framework. Water present in these channels gradually coalesces and fills the hexagonal pores sequentially after the channels are completely hydrated. The development of hydrogen-bond networks inside the MOF pores as a function of relative humidity is characterized at the molecular level using experimental and computational infrared spectroscopy. A detailed analysis of the OH-stretch vibrational band indicates that the low-frequency tail stems from strongly polarized hydrogen-bonded water molecules, suggesting the presence of some structural disorder in the experimental samples. Strategies for designing efficient water harvesting MOFs are also proposed based on the mechanism of water adsorption in NU-1500-Cr.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Dominant balance-based adaptive mesh refinement for incompressible fluid flows

This work introduces a novel adaptive mesh refinement (AMR) method that utilizes dominant balance analysis (DBA) for efficient and accurate grid adaptation in computational fluid dynamics (CFD) simulations. The proposed method leverages a Gaussian mixture model (GMM) to classify grid cells into active and passive regions based on the dominant physical interactions within the equation space. By modeling truncation error probabilistically from discretized terms, the method identifies regions of high interaction where numerical accuracy is most sensitive to resolution. Unlike traditional AMR strategies, this approach does not rely on heuristic-based sensors or user-defined thresholds, providing a fully automated and problem-independent framework for AMR. Applied to the incompressible Navier-Stokes equations for steady and unsteady flow past a cylinder, the DBA-based AMR method achieves comparable accuracy to high-resolution grids while reducing computational costs by up to 70 %. The validation highlights the method’s effectiveness in capturing complex flow features while minimizing grid cells, directing computational resources toward regions with the most critical dynamics. This modular and scalable strategy is adaptable to a wide range of applications, presenting a promising tool for efficient high-fidelity simulations in CFD and other multiphysics domains.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Optimum power analysis of a self-reactive wave energy point absorber with mechanically-driven power take-offs

This paper studies a self-reactive ocean wave energy converter (WEC) where energy is absorbed through the relative motion between a floating buoy on the ocean surface and a submerged body. Two types of direct-drive power take off (PTO) systems are examined. One adopts a ball screw system in inverse driving mode with mechanical motion rectifier (MMR) that converts the oscillating bi-directional input heave motion into unidirectional rotation to drive an electromagnetic generator, and another one uses a ball screw system to directly drive the generator (called non-MMR). Dynamic models for both types of PTOs are established and integrated with the overall WEC system model in both the time and frequency domains. Numerical simulation is used to investigate the dynamic performance of both PTOs. The influences of the PTO inerter, which is from the rotational inertia of the flywheel and generator, and the WEC drag damping coefficient are discussed and explored. The analytical closed-form solution for the optimum condition of the non-MMR system is derived based on the criterion of maximum power absorption, and numerical simulation yields the optimum condition for the MMR system. Overall, the results under regular and irregular waves show that the PTO inerter in the MMR system can improve the power absorption for small wave periods and maintain the same performance for large wave periods. The PTO inerter in the non-MMR system shifts the peak frequency. The WEC drag damping negatively influences the power absorption in both PTO systems, especially at large wave periods.

16 TIDAL AND WAVE POWER↗
DOE OSTI2023

Investigations of nonlinear polarization transfer between obliquely intersecting beams

Electromagnetic waves propagating through plasma can interact nonlinearly through a variety of different mechanisms. The excitation of a plasma beat wave (ions or electrons) can create a refractive index modulation that changes the dispersion of the interacting beams. Alternatively, high-intensity beams can enter the regime where relativistic nonlinearities influence the propagation dynamics. In recent studies [ Opt. Express 29 , 1162 ( 2021 ) OPEXFF 1094-4087 10.1364/OE.413064 ], it was proposed that two beams propagating along the same axis can exchange their polarization state due to nonlinear interaction. Here we present a numerical analysis of two laser beams intersecting in a nonlinear medium at varying angles. Polarization transfer is observed as predicted by analytical theory for a range of angles. For small angles, it is found that filamentation of the interacting beams becomes important. Analytical estimates of the filamentation threshold are presented, and good agreement is found with the simulation data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A Novel Experimental Approach to Understand the Transport of Nanodrugs

Nanoparticle-based drugs offer attractive advantages like targeted delivery to the diseased site and size and shape-controlled properties. Therefore, understanding the particulate flow of the nanodrugs is important for effective delivery, accurate prediction of required dosage, and developing efficient drug delivery platforms for nanodrugs. In this study, the transport of nanodrugs including flow velocity and deposition is investigated using three model metal oxide nanodrugs of different sizes including iron oxide, zinc oxide, and combined Cu-Zn-Fe oxide synthesized via a modified polyol approach. The hydrodynamic size, size, morphology, chemical composition, crystal phase, and surface functional groups of the water-soluble nanodrugs were characterized via dynamic light scattering, transmission electron microscopy, scanning electron microscopy-energy dispersive X-ray, X-ray diffraction, and fourier transform infrared spectroscopy, respectively. Two different biomimetic flow channels with customized surfaces are developed via 3D printing to experimentally monitor the velocity and deposition of the different nanodrugs. A diffusion dominated mechanism of flow is seen in size ranges 92 nm to 110 nm of the nanodrugs, from the experimental velocity and mass loss profiles. The flow velocity analysis also shows that the transport of nanodrugs is controlled by sedimentation processes in the larger size ranges of 110–302 nm. However, the combined overview from experimental mass loss and velocity trends indicates presence of both diffusive and sedimentation forces in the 110–302 nm size ranges. It is also discovered that the nanodrugs with higher positive surface charges are transported faster through the two test channels, which also leads to lower deposition of these nanodrugs on the walls of the flow channels. The results from this study will be valuable in realizing reliable and cost-effective in vitro experimental approaches that can support in vivo methods to predict the flow of new nanodrugs.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Pore-Scale Dynamics of Liquid CO2–Water Displacement in 2D Axisymmetric Porous Micromodels Under Strong Drainage and Weak Imbibition Conditions: High-Speed μPIV Measurements

Resolving pore-scale transient flow dynamics is crucial to understanding the physics underlying multiphase flow in porous media and informing large-scale predictive models. Surface properties of the porous matrix play an important role in controlling such physics, yet interfacial mechanisms remain poorly understood, in part due to a lack of direct observations. This study reports on an experimental investigation of the pore-scale flow dynamics of liquid CO 2 and water in two-dimensional (2D) circular porous micromodels with different surface characteristics employing high-speed microscopic particle image velocimetry (μPIV). The design of the micromodel minimized side boundary effects due to the limited size of the domain. The high-speed μPIV technique resolved the spatial and temporal dynamics of multiphase flow of CO 2 and water under reservoir-relevant conditions, for both drainage and imbibition scenarios. When CO 2 displaced water in a hydrophilic micromodel (i.e., drainage), unstable capillary fingering occurred and the pore flow was dominated by successive pore-scale burst events (i.e., Haines jumps). When the same experiment was repeated in a nearly neutral wetting micromodel (i.e., weak imbibition), flow instability and fluctuations were virtually eliminated, leading to a more compact displacement pattern. Energy balance analysis indicates that the conversion efficiency between surface energy and external work is less than 30%, and that kinetic energy is a disproportionately smaller contributor to the energy budget. This is true even during a Haines jump event, which induces velocities typically two orders of magnitude higher than the bulk velocity. These novel measurements further enabled direct observations of the meniscus displacement, revealing a significant alteration of the pore filling mechanisms during drainage and imbibition. While the former typically featured burst events, which often occur only at one of the several throats connecting a pore, the latter is typically dominated by a cooperative filling mechanism involving simultaneous invasion of a pore from multiple throats. This cooperative filling mechanism leads to merging of two interfaces and releases surface energy, causing instantaneous high-speed events that are similar, yet fundamentally different from, burst events. Finally, pore-scale velocity fields were statistically analyzed to provide a quantitative measure of the role of capillary effects in these pore flows.

Li, Yaofa↗