Search NASA⌕ Search

SEARCH · Search NASA

Results for “controller stability”

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

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

At least 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Core plasma physics basis and its impacts on the FNSF

The FNSF core plasma physics is reportede on with detailed analysis to establish the basis for the reference configuration, R = 4.8 m, a = 1.2 m, Ip = 7.87 MA, and B T = 7.5 T, q 95 = 6, f BS = 0.52, β N < 2.7 established in Ref. [1]. Central solenoid (CS) and poloidal field (PF) coils are far from the plasma as in a power plant, and acceptable coil currents are determined for the rampup and flattop fiducial states. Time-dependent free-boundary plasma evolution simulations demonstrate that the FNSF plasma can be established, ramped up, and relaxed into flattop, including vertical stabilizers, internal feedback coils and feedback control on plasma current, position, and shape. A range of density (n o / = 1.3–1.5) and temperature (T o / = 2.2–2.7) profiles are examined, indicating that energy confinement of H 98 = 1.1–1.2 is required to provide 100% non-inductive plasma current in the FNSF. GLF23 theory based transport model predicted lower energy confinement of H 98 ~0.6–0.85. The EPED analysis shows that the pedestal temperature ranges from 4.0–4.7 keV for pedestal densities of 1.7–1.0 × 10 20 /m 3 . The n = 1 kink stability shows no-wall beta limits, using the pressure and current profiles associated with the transport and current drive sources, ranging from β N ~2.25–2.55 depending on l i . A conducting wall can extend these limits by 10–40% depending on li and wall location. At the lower beta’s of the reference plasma, a combination of 50 MW of NB, 30 MW of LH, 20 MW of ICRF, 20 MW of EC, and bootstrap current, are found to provide 100% of the plasma current with a stable current profile. Impacts on the FNSF of plasma physics are discussed and R&D challenges are recorded.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Erosional redistribution of topsoil controls soil nitrogen dynamics

In recent years, the role of soil erosion on terrestrial carbon sequestration had been the focus of a growing number of studies. However, relatively little attention has been paid so far to the role of erosion on the lateral distribution of soil nitrogen (N) and the role of geomorphic processes on soil N dynamics. We present primary data on the stock of nitrogen in soil and its rate of erosion at a relatively undisturbed, zero-order watershed in northern California. Erosion transports 0.26–0.47 g N m -2 year -1 from eroding slope positions (Summit and Slope), and about two-thirds of the eroded N enters depositional landform positions (Hollow and Plain). Our results show that depositional-position soil profiles contain up to 3 times more N than soil profiles in the eroding positions. More than 92% of all soil nitrogen was chemically bound to soil minerals in all the landform positions, compared to 2–4% each found in the free light and occluded light fractions. Nitrogen associated with the free light fraction in topsoil is particularly susceptible to loss by soil erosion. By comparison, soil N associated with the aggregate-protected occluded light fractions and the mineral-associated dense fractions is likely to be protected from gaseous and dissolved losses. On average, we found that soil N has mean residence time of 694 years in eroding landform positions, compared to 2951 years in depositional landform positions. Our results also show that microbial processing of organic matter exerts strong control on overall soil N storage and N stabilized through sorptive interactions with soil minerals only in poorly drained depositional landform positions. Finally, soil erosion exerts important control on stock, distribution, and long-term fate of soil N in dynamic landscapes.

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

High temperature shockwave stabilized single atoms

The stability of single-atom catalysts is critical for their practical applications. Although a high temperature can promote the bond formation between metal atoms and the substrate with an enhanced stability, it often causes atom agglomeration and is incompatible with many temperature-sensitive substrates. Here, we report using controllable high-temperature shockwaves to synthesize and stabilize single atoms at very high temperatures (1,500-2,000 K), achieved by a periodic on-off heating that features a short on state (55 ms) and a ten-times longer off state. The high temperature provides the activation energy for atom dispersion by forming thermodynamically favourable metal-defect bonds and the off-state critically ensures the overall stability, especially for the substrate. The resultant high-temperature single atoms exhibit a superior thermal stability as durable catalysts. In conclusion, the reported shockwave method is facile, ultrafast and universal (for example, Pt, Ru and Co single atoms, and carbon, C 3 N 4 and TiO 2 substrates), which opens a general route for single-atom manufacturing that is conventionally challenging.

42 ENGINEERING↗
DOE OSTI2023

Structural water in amorphous carbonate minerals: ab initio molecular dynamics simulations of X-ray pair distribution experiments

Water is known to play a controlling role in directing mineralization pathways and stabilizing metastable amorphous intermediates in hydrous carbonate mineral MCO 3 ·nH 2 O systems, where M 2+ is a divalent metal cation. Despite this recognition, the nature of the controls on crystallization is poorly understood, largely owing to the difficulty in characterizing the dynamically disordered structures of amorphous intermediates at the atomic scale. Here, we present a series of atomistic models, derived from ab initio molecular dynamics simulation, across a range of experimentally relevant cations (M = Ca, Mg, Sr) and hydration levels (0 ≤ n ≤ 2). Theoretical simulations of the dependence of the X-ray pair distribution function on the hydration level n show good agreement with available experimental data and thus provide further evidence for a lack of significant nanoscale structure in amorphous carbonates. Upon dehydration, the metal coordination number does not change significantly, but the relative extent of water dissociation increases, indicating that a thermodynamic driving force exists for water dissociation to accompany dehydration. Mg strongly favors monodentate conformation of carbonate ligands and shows a marked preference to exchange monodentate carbonate O for water O upon hydration, whereas Ca and Sr exchange mono- and bidentate carbonate ligands with comparable frequency. Finally, water forms an extensive hydrogen bond network among both water and carbonate groups that exhibits frequent proton transfers for all three cations considered suggesting that proton mobility is likely predominantly due to water dissociation and proton transfer reactions rather than molecular water diffusion.

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

A gimbal platform stabilization for topographic applications

The aim of this work is the stabilization of a Gimbal platform for optical sensors acquisitions in topographic applications using mobile vehicles. The stabilization of the line of sight (LOS) consists in tracking the command velocity in presence of nonlinear noise due to the external environment. The hardware architecture is characterized by an Ardupilot platform that allows the control of both the mobile device and the Gimbal. Here we developed a new approach to stabilize the Gimbal platform, which is based on neural network. For the control system, we considered a plant that represents the transfer function of the servo system control model for an inertial stabilized Gimbal platform. The transductor used in the feed-back line control is characterized by the Rate Gyro transfer function installed onboard of Ardupilot. For the simulation and investigation of the system performance, we used the Simulink tool of Matlab. Results show that the hardware/software approach is efficient, reliable and cheap for direct photogrammetry, as well as for general purpose applications using mobile vehicles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Robust H{sub ∞} control for singular time-delay systems with saturating actuators via static output feedback

The problem of H{sub ∞} robust exponential stability and static output feedback stabilization of uncertain singular time-delay systems with saturating actuators are investigated. First, the delay-dependent condition is given to guarantee that the singular time-delay systems are not only exponential admissible, but also satisfies H{sub ∞} performance γ. Moreover, the estimation of stability region and the design method of H{sub ∞} static output feedback controller are given by solving a convex optimization problem. Some simulation examples are provided to demonstrate the merit of the proposed design methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Making the Grid "Smart" Through "Smart" Microgrids: Real-Time Power Management of Microgrids with Multiple Distributed Generation Sources Using Intelligent Control

In this Project we collaborated with two DOE National Laboratories, Pacific Northwest National Lab (PNNL) and Lawrence Berkeley National Lab (LBL). Dr. Hammerstrom of PNNL initially supported our project and was on the graduate committee of one of the Ph.D. students (graduated in 2014) who was supported by this project. He is also a committee member of a current graduate student of the PI who was supported by this project in the last two years (August 2014-July 2016). The graduate student is now supported be the Electrical and Computer Engineering (ECE) Department at Montana State University (MSU). Dr. Chris Marney of LBL provided actual load data, and the software WEBOPT developed at LBL for microgrid (MG) design for our project. NEC-Labs America, a private industry, also supported our project, providing expert support and modest financial support. We also used the software “HOMER,” originally developed at the National Renewable Energy Laboratory (NREL) and the most recent version made available to us by HOMER Energy, Inc., for MG (hybrid energy system) unit sizing. We compared the findings from WebOpt and HOMER and designed appropriately sized hybrid systems for our case studies. The objective of the project was to investigate real-time power management strategies for MGs using intelligent control, considering maximum feasible energy sustainability, reliability and efficiency while, minimizing cost and undesired environmental impact (emissions). Through analytic and simulation studies, we evaluated the suitability of several heuristic and artificial-intelligence (AI)-based optimization techniques that had potential for real-time MG power management, including genetic algorithms (GA), ant colony optimization (ACO), particle swarm optimization (PSO), and multi-agent systems (MAS), which is based on the negotiation of smart software-based agents. We found that PSO and MAS, in particular, distributed MAS, were more efficient and better suited for our work. We investigated the following: • Intelligent load control - demand response (DR) - for frequency stabilization in islanded MGs (partially supported by PNNL). • The impact of high penetration of solar photovoltaic (PV)-generated power at the distribution level (partially supported by PNNL). • The application of AI approaches to renewable (wind, PV) power forecasting (proposed by the reviewers of our proposal). • Application of AI approaches and DR for real-time MG power management (partially supported by NEC Labs-America) • Application of DR in dealing with the variability of wind power • Real-time MG power management using DR and storage (partially supported by NEC Labs-America) • Application of DR in enhancing the performance of load-frequency controller • MAS-based whole-sale and retail power market design for smart grid A

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Stability of multi-population traffic flows

Traffic waves, known also as stop-and-go waves or phantom jams, appear naturally as traffic instabilities, also in confined environments as a ring-road. A multi-population traffic is studied on a ring-road, comprised of drivers with stable and unstable behavior. There exists a critical penetration rate of stable vehicles above which the system is stable, and under which the system is unstable. In the latter case, stop-and-go waves appear, provided enough cars are on the road. The critical penetration rate is explicitly computable, and, in reasonable situations, a small minority of aggressive drivers is enough to destabilize an otherwise very stable flow. This is a source of instability that a single population model would not be able to explain. Also, the multi-population system can be stable below the critical penetration rate if the number of cars is sufficiently small. Instability emerges as the number of cars increases, even if the traffic density remains the same (i.e., number of cars and road size increase similarly). This shows that small experiments could lead to deducing imprecise stability conditions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

A new passivity-based nonlinear causal control technique for wave energy converters with finite stroke

A new technique is proposed for the design of nonlinear causal feedback controllers for wave energy converters, which explicitly protect against stroke saturation. Here, the technique is an improvement of a technique originally proposed. The primary challenge associated with the problem is the design of nonlinear controllers for which global closed-loop stability can be assured with probability 1 for stationary stochastic response, without also unduly sacrificing power generation performance. The proposed technique consists of three steps: 1) designing a linear feedback controller using multi-objective optimization techniques; 2) augmenting this design with an extra input channel that adheres to a closed-loop passivity condition; and 3) designing an outer, nonlinear passive feedback loop that controls this augmented input in such a way as to ensure stroke limits are maintained. The primary contribution of this paper is a fundamental improvement to the second of these steps, for which the previous technique was not guaranteed to return a valid solution. The proposed technique, by contrast, is shown to always produce a nonlinear controller with the desired properties, under very general modeling assumptions. The technique is illustrated on a simple heaving-buoy wave energy converter example.

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

Pointing stabilization of a 1 Hz high-power laser via machine learning

Abstract High-power lasers are vital for particle acceleration, imaging, fusion and materials processing, requiring precise control and high-energy delivery. Laser plasma accelerators (LPAs) demand laser positional stability at focus to ensure consistent electron beams in applications such as X-ray free-electron lasers and high-energy colliders. Achieving this stability is especially challenging for the low-repetition-rate lasers in current LPAs. We present a machine learning method that predicts and corrects laser pointing instabilities in real-time using a high-frequency pilot beam. By preemptively adjusting a correction mirror, this approach overcomes traditional feedback limits. Demonstrated on the BELLA petawatt laser operating at the terawatt level (30 mJ amplification), our method achieved root mean square pointing stabilization of 0.34 and 0.59 $\unicode{x3bc} \mathrm{rad}$ in the x and y directions, reducing jitter by 65% and 47%, respectively. This is the first successful application of predictive control for shot-to-shot stabilization in low-repetition-rate laser systems, paving the way for full-energy petawatt lasers and transformative advances across science, industry and security.

Amodio, Alessio↗
DOE OSTI2024

Slug-flow synthesis of NCMA: Effect of substitution of cobalt with aluminum on the electrochemical performance of Ni-rich cathode for lithium-ion battery

Nickle-rich Li[Ni 1-x-y Co x Mn y ]O 2 (x, y≤0.1) (NCM) layered materials are known as promising cathode materials for next-generation lithium-ion batteries and electric vehicles owing to their high-reversible capacity and operating voltage of up to 3.6 vs Li/Li + . However, issues, such as irreversible phase transition, cation mixing, microcrack formation, thermal and structural stability of the material prevent its widespread adoption. Although, cation doping is a well-known technique to enhance the electrochemical performance of the NCM-based cathode material, the performance of the material is very sensitive to the doping amount. In this study, three Al-doped quaternary Ni-rich cathode materials Li[Ni 0.85 Co (0.1-x) Mn 0.05 Al x ]O 2 (where, x=0–0.04) (NCMA) are synthesized through three-phase slug-flow based continuous manufacturing process followed by high temperature calcination to study the effect of Aldoping on the performance of the cathode material while reducing Co. The slug flow-based production platform has several advantages, like particle size uniformity, high production rate, and homogeneity in elemental distribution. It is found that with an increase in Al content, the specific capacity decreases but the cyclic stability and rate capability increases. Furthermore, optimum Al-doping not only compensates for the adverse effect of low Co by decreasing the extent of cation mixing, but it also minimizes the electrode polarization and cracking of the particles.

25 ENERGY STORAGE↗
DOE OSTI2022

Using Al 3+ to Tailor Graphene Oxide Nanochannels: Impact on Membrane Stability and Permeability

Graphene oxide (GO) membranes, which form from the lamination of GO sheets, attract much attention due to their unique nanochannels. There is much interest in controlling the nanochannel structures and improving the aqueous stability of GO membranes so they can be effectively used in separation and filtration applications. This study employed a simple yet effective method of introducing trivalent aluminum cations to a GO sheet solution through the oxidation of aluminum foil, which modifies the nanochannels in the self-assembled GO membrane by increasing the inter-sheet distance while decreasing intra-sheet spacing. The Al 3+ modification resulted in an increase in membrane stability in water, methanol, ethanol, and propanol, yet decreased membrane permeability to water and propanol. These changes were attributed to strong interactions between Al 3+ and the membrane oxygenated functional groups, which resulted in an increase in membrane hydrophobicity and a decrease in the intra-sheet spacing as supported by surface tension, contact angle, atomic force microscopy, and X-ray photoelectron spectroscopy measurements. Our approach for forming Al 3+ modified GO membranes provides a method for improving the aqueous stability and tailoring the permeation selectivity of GO membranes, which have the potential to be implemented in vapor separation and fuel purification applications.

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

NSLS-II storage ring insertion device and front-end commissioning and operation

The National Synchrotron Light Source II (NSLS-II) is a state of the art 3 GeV third generation light source at Brookhaven National Laboratory. During spring/ summer of 2014, the storage ring was commissioned up to 50 mA without insertion devices. In the fall of 2014, we began commissioning of the project beamlines, which included seven insertion devices on six ID ports. Beamlines IXS, HXN, CSX-1, CSX-2, CHX, SRX, and XPD-1 consist of elliptically polarized undulator (EPU), damping wigglers (DW) and in-vacuum undulators (IVU) covering from VUV to hard x-ray range. In this paper, experience with commissioning and operation is discussed. We focus on reaching storage ring performance with IDs, including injection, design emittance, compensation of orbit distortions caused by ID residual field, source point stability, beam alignment and tools for control, monitoring and protection of the ring chambers from ID radiation.

43 PARTICLE ACCELERATORS↗
DOE OSTI2015

Reliability of unstable periodic orbit based control strategies in biological systems

Presence of recurrent and statistically significant unstable periodic orbits (UPOs) in time series obtained from biological systems is now routinely used as evidence for low dimensional chaos. Extracting accurate dynamical information from the detected UPO trajectories is vital for successful control strategies that either aim to stabilize the system near the fixed point or steer the system away from the periodic orbits. A hybrid UPO detection method from return maps that combines topological recurrence criterion, matrix fit algorithm, and stringent criterion for fixed point location gives accurate and statistically significant UPOs even in the presence of significant noise. Geometry of the return map, frequency of UPOs visiting the same trajectory, length of the data set, strength of the noise, and degree of nonstationarity affect the efficacy of the proposed method. Results suggest that establishing determinism from unambiguous UPO detection is often possible in short data sets with significant noise, but derived dynamical properties are rarely accurate and adequate for controlling the dynamics around these UPOs. A repeat chaos control experiment on epileptic hippocampal slices through more stringent control strategy and adaptive UPO tracking is reinterpreted in this context through simulation of similar control experiments on an analogous but stochastic computer model of epileptic brain slices. Reproduction of equivalent results suggests that far more stringent criteria are needed for linking apparent success of control in such experiments with possible determinism in the underlying dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Multiple-time-stepping generalized hybrid Monte Carlo methods

Performance of the generalized shadow hybrid Monte Carlo (GSHMC) method [1], which proved to be superior in sampling efficiency over its predecessors [2–4], molecular dynamics and hybrid Monte Carlo, can be further improved by combining it with multi-time-stepping (MTS) and mollification of slow forces. We demonstrate that the comparatively simple modifications of the method not only lead to better performance of GSHMC itself but also allow for beating the best performed methods, which use the similar force splitting schemes. In addition we show that the same ideas can be successfully applied to the conventional generalized hybrid Monte Carlo method (GHMC). The resulting methods, MTS-GHMC and MTS-GSHMC, provide accurate reproduction of thermodynamic and dynamical properties, exact temperature control during simulation and computational robustness and efficiency. MTS-GHMC uses a generalized momentum update to achieve weak stochastic stabilization to the molecular dynamics (MD) integrator. MTS-GSHMC adds the use of a shadow (modified) Hamiltonian to filter the MD trajectories in the HMC scheme. We introduce a new shadow Hamiltonian formulation adapted to force-splitting methods. The use of such Hamiltonians improves the acceptance rate of trajectories and has a strong impact on the sampling efficiency of the method. Both methods were implemented in the open-source MD package ProtoMol and were tested on a water and a protein systems. Results were compared to those obtained using a Langevin Molly (LM) method [5] on the same systems. The test results demonstrate the superiority of the new methods over LM in terms of stability, accuracy and sampling efficiency. This suggests that putting the MTS approach in the framework of hybrid Monte Carlo and using the natural stochasticity offered by the generalized hybrid Monte Carlo lead to improving stability of MTS and allow for achieving larger step sizes in the simulation of complex systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

The Inner Shell Spectroscopy beamline at NSLS-II: a facility for in situ and operando X-ray absorption spectroscopy for materials research

The Inner Shell Spectroscopy (ISS) beamline on the 8-ID station at the National Synchrotron Light Source II (NSLS-II), Upton, NY, USA, is a high-throughput X-ray absorption spectroscopy beamline designed for in situ, operando, and time-resolved material characterization using high monochromatic flux and scanning speed. This contribution discusses the technical specifications of the beamline in terms of optics, heat load management, monochromator motion control, and data acquisition and processing. Results of the beamline tests demonstrating the quality of the data obtainable on the instrument, possible energy scanning speeds, as well as long-term beamline stability are shown. The ability to directly control the monochromator trajectory to define the acquisition time for each spectral region is highlighted. Examples of studies performed on the beamline are presented. The paper is concluded with a brief outlook for future developments.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Manipulating Conjugated Polymer Backbone Dynamics through Controlled Thermal Cleavage of Alkyl Side Chains

The morphological stability of an organic photovoltaic (OPV) device is greatly affected by the dynamics of donors and acceptors occurring near the device's operational temperature. These dynamics can be quantified by the glass transition temperature (T g ) of conjugated polymers (CPs). Because flexible side chains possess much faster dynamics, the cleavage of the alkyl side chains will reduce chain dynamics, leading to a higher T g . In this work, the T g s for CPs are systematically studied with controlled side chain cleavage. Isothermal annealing of polythiophenes featuring thermally cleavable side chains at 140 °C, is found to remove more than 95% of alkyl side chains in 24 h, and raise the backbone T g from 23 to 75 °C. Coarse grain molecular dynamics simulations are used to understand the T g dependence on side chain cleavage. X-ray scattering indicates that the relative degree of crystallization remains constantduring isothermal annealing process. The effective conjugation length is not influenced by thermal cleavage; however, the density of chromophore is doubled after the complete removal of alkyl side chains. Here, the combined effect of enhancing T g and conserving crystalline structures during the thermal cleavage process can provide a pathway to improving the stability of optoelectronic properties in future OPV devices.

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

Mechanisms Leading to Co-Existence of Gas Hydrate in Ocean Sediments [Part 2 of 2]

In this project we have sought to explain the co-existence of gas and hydrate phases in sediments within the gas hydrate stability zone. We have focused on the gas/brine interface at the scale of individual grains in the sediment. The capillary forces associated with a gas/brine interface play a dominant role in many processes that occur in the pores of sediments and sedimentary rocks. The mechanical forces associated with the same interface can lead to fracture initiation and propagation in hydrate-bearing sediments. Thus the unifying theme of the research reported here is that pore scale phenomena are key to understanding large scale phenomena in hydrate-bearing sediments whenever a free gas phase is present. Our analysis of pore-scale phenomena in this project has delineated three regimes that govern processes in which the gas phase pressure is increasing: fracturing, capillary fingering and viscous fingering. These regimes are characterized by different morphology of the region invaded by the gas. On the other hand when the gas phase pressure is decreasing, the corresponding regimes are capillary fingering and compaction. In this project, we studied all these regimes except compaction. Many processes of interest in hydrate-bearing sediments can be better understood when placed in the context of the appropriate regime. For example, hydrate formation in sub-permafrost sediments falls in the capillary fingering regime, whereas gas invasion into ocean sediments is likely to fall into the fracturing regime. Our research provides insight into the mechanisms by which gas reservoirs are converted to hydrate as the base of the gas hydrate stability zone descends through the reservoir. If the reservoir was no longer being charged, then variation in grain size distribution within the reservoir explain hydrate saturation profiles such as that at Mt. Elbert, where sand-rich intervals containing little hydrate are interspersed between intervals containing large hydrate saturations. Large volumes (of order one pore volume) of gaseous and aqueous phases must be transported into the gas hydrate stability zone. The driver for this transport is the pressure sink induced by a reduction in occupied pore volume that accompanies the formation of hydrate from gas and water. Pore-scale imbibition models and bed-scale multiphase flow models indicate that the rate-limiting step in converting gas to hydrate is the supply of water to the hydrate stability zone. Moreover, the water supply rate is controlled by capillarity-driven flux for conditions typical of the Alaska North Slope. A meter-scale laboratory experiment confirms that significant volumes of fluid phases move into the hydrate stability zone and that capillarity is essential for the water flux. The model shows that without capillarity-driven flux, large saturations of hydrate cannot form. The observations of thick zones of large saturation at Mallik and Mt Elbert thus suggest that the primary control on these systems is the rate of transport of gaseous and aqueous phases, driven by the pressure sink at the base of the gas hydrate stability zone. A key finding of our project is the elucidation of ?capillary fracturing? as a dominant gas transport mechanism in low-permeability media. We initially investigate this phenomenon by means of grain-scale simulations in which we extended a discrete element mechanics code (PFC, by Itasca) to incorporate the dynamics of first singlephase and then multiphase flow. A reductionist model on a square lattice allows us to determine some of the fundamental dependencies of the mode of gas invasion (capillary fingering, viscous fingering, and fracturing) on the parameters of the system. We then show that the morphology of the gas-invaded region exerts a fundamental control on the fabric of methane hydrate formation, and on the overpressures caused by methane hydrate dissociation. We demonstrate the existence of the different invasion regimes by means of controlled laboratory experiments in a radial cell. We collapse the behavior in the form of a phase diagram fully characterized by two dimensionless groups: a modified capillary number and a ?fracturing number? that reflects the balance between the pressure forces that act to open conduits in the granular pack, and frictional forces that resist it. We use all this small-scale knowledge to propose simple mechanistic models of gas migration and hydrate formation at the geologic bed scale. We propose that methane transport in lake and oceanic sediments is controlled by dynamic conduits, which dilate and release gas as the falling hydrostatic pressure reduces the effective stress below the tensile strength of the sediments. We test our model against a four-month record of hydrostatic load and methane flux in Upper Mystic Lake, Mass., USA, and show that it captures the complex episodicity of methane ebullition. Our quantitative conceptualization opens the door to integrated modeling of methane transport to constrain global methane release from lakes and other methane-rich sediment systems, and to assess its climate feedbacks.

58 GEOSCIENCES↗