SEARCH · Engineering Papers

Results for “3D gravity”

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

SourceResearch connectionFrom yearThrough year

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

At least 181 records · Page 10

All sources NASA NTRS DOE OSTI

NASA NTRS2006

Modeling the Observed QBO and Inter-Annual Variations of the Diurnal Tide in the Mesosphere

In the current version of the Numerical Spectral Model (NSM), the Quasi-biennial Oscillation (QBO) is generated primarily by small-scale gravity waves (GW) from Hines' Doppler Spread Parameterization (DSP). The model does not have topography, and the planetary waves are solely generated by instabilities. We discuss a 3D modeling study that describes the QBO extending from the stratosphere into the upper mesosphere, where the oscillation produces significant inter-annual variations in the diurnal tide. The numerical results are compared with temperature measurements from the SABER (TIMED) and MLS (UARS) instruments obtained by Huang et al. (2006). With a GW source that peaks at the Equator and is taken to be isotropic and independent of season, the NSM generates a QBO with variable periods around 26 months and zonal wind amplitudes of almost 25 m/s at 30 km. As reported earlier, the NSM reproduces the observed equinoctial maxima in the diurnal tide at altitudes around 95 km. The modeled QBO modulates the tide such that the seasonal amplitude maxima can vary from one year to another by as much as 30%. To shed light on the underlying mechanisms, the relative importance of the advection terms are discussed, and they are shown to be important in the stratosphere. At altitudes above 80 km, however, the QBO-related inter-annual variations of the tide are generated primarily by GW momentum deposition. In qualitative agreement with the SABER measurements, the model generates distinct zonal-mean QBO temperature variations in the stratosphere and mesosphere. In the stratosphere, the computed amplitudes are not much smaller than those observed, and the rate of downward propagation at the Equator is reproduced. The modeled temperature amplitudes in the mesosphere, however, are much smaller than those observed. The observed and computed temperature variations of the QBO peak at the Equator but extend with phase reversals to high latitudes, in contrast to the zonal winds that are confined to equatorial latitudes. Hemispherical asymmetries also appear in both the model results and the observations. The temperature amplitudes outside the equatorial region however tend to occur at lower latitudes in the model results. While there is qualitative agreement between the TIMED measurements and the model prediction, there are some areas of significant disagreement that require us to reexamine the present version of the NSM. The numerical results critically depend on the chosen parameters that determine the wave forcing, and there are a number of avenues to improve the performance of the model that had not been tuned to fit the observations. The GW spectrum and its latitude dependence in the troposphere are not well known, and numerical experiments are discussed that describe the related model response. While it appears that eastward propagating Kelvin waves and westward propagating Rossby gravity waves are not the primary source to generate the QBO, the GW forcing can seed the oscillation and act as a catalyst to enhance effectiveness of these planetary waves.

Mayr, Hans G.↗

DOE OSTI2024

High speed vibration compensation using magnetic fiducials via NIMBLE

Additive manufacturing (AM) technology, specifically 3D printing, holds great promise for in-orbit manufacturing. In-space printing can significantly reduce the mass, cost, and risk of long-term space exploration by enabling replacement parts to be made as needed and reducing dependence on Earth. However, printing in a zero-gravity environment poses challenges due to the absence of a rigid ground for the print platform, which can result in vibrational and rotational forces that may impact printing integrity. Here, to address this issue, this paper proposes a novel linear magnetic position tracking algorithm, named Navigation Integrating Magnets By Linear Estimation (NIMBLE), for dynamic vibration compensation during 3D printing of truss structures in space. Compared to the most commonly used nonlinear optimization method, the NIMBLE algorithm is more than two orders of magnitude faster. With only a single 3-axis magnet sensor and a small NdFeB magnet, the NIMBLE algorithm provides a simple and easily implemented tracking solution for in-orbit 3D printing.

3D printing in space↗

DOE OSTI2012

Astrophysical tests of gravity: a screening map of the nearby universe

Astrophysical tests of modified gravity theories in the nearby universe have been emphasized recently by Hui 2009 and Jain 2011. A key element of such tests is the screening mechanism whereby general relativity is restored in massive halos or high density environments like the Milky Way. In chameleon theories of gravity, including all f(R) models, field dwarf galaxies may be unscreened and therefore feel an extra force, as opposed to screened galaxies. The first step to study differences between screened and unscreened galaxies is to create a 3D screening map. We use N-body simulations to test and calibrate simple approximations to determine the level of screening in galaxy catalogs. Sources of systematic errors in the screening map due to observational inaccuracies are modeled and their contamination is estimated. We then apply our methods to create a map out to 200 Mpc in the Sloan Digital Sky Survey footprint using data from the Sloan survey and other sources. In two companion papers this map will be used to carry out new tests of gravity using distance indicators and the disks of dwarf galaxies. We also make our screening map publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2002

The Influence of Tropospheric Processes in Modeling the Middle Atmosphere with Gravity Waves

Our Numerical Spectral Model (NSM) extends from the ground up into the thermosphere and has a vertical grid point resolution of about 0.5 km to resolve the interactions of gravity waves (GWs) described with Hines' Doppler Spread Parameterization (DSP). This model produces in the stratosphere and mesosphere the major features of QBO, SAO, tides, and planetary waves. The purpose of this paper is to discuss results from an initial study with our 3D model that shows how certain tropospheric processes can affect the dynamics of the middle atmosphere. Under the influence of tropospheric heating, and augmented by GW interactions, two distinct but related processes can be identified. (1) A meridional circulation develops in the stratosphere, with rising motions at low latitudes that are in magnitude comparable to the downward propagation of the QBO. As Dunkerton pointed out, a larger GW source is then required to reproduce the observed QBO, which tends to move us closer to the values recommended for the DSP. This has significant consequences for our model results that describe the upper mesosphere, considering the general importance of GWs for this region and in influencing planetary waves (e.g., 2-day wave) and tides in particular. (2) Tropospheric heating produces zonal jets near the tropopause that are related to latitudinal variations in pressure and reversing temperature variations (resembling the dynamical conditions near the mesopause), which in turn is conducive to generate baroclinic instability. Modeling results show that our ability to generate the QBO critically depends on the magnitude of the temperature reversal that is a measure of this instability. Planetary waves are generated in this process, which can apparently interfere with or augment the GW interactions. As originally demonstrated by Lindzen and Holton, the eastward propagating Kelvin waves and westward propagating Rossby gravity waves (generated by tropospheric convection) can in principle provide the acceleration to influence the QBO, and we were able to confirm this with our 3D model.

Mayr, H. G.↗

DOE OSTI2023

Secular outflows from 3D MHD hypermassive neutron star accretion disc systems

Magnetized hypermassive neutron stars (HMNSs) have been proposed as a way for neutron star mergers to produce high electron fraction, high-velocity ejecta, as required by kilonova models to explain the observed light curve of GW170817. The HMNS drives outflows through neutrino energy deposition and mechanical oscillations, and raises the electron fraction of outflows through neutrino interactions before collapsing to a black hole (BH). Here, in this study, we perform 3D numerical simulations of HMNS–torus systems in ideal magnetohydrodynamics, using a leakage/absorption scheme for neutrino transport, the nuclear APR equation of state, and Newtonian self-gravity, with a pseudo-Newtonian potential added after BH formation. Due to the uncertainty in the HMNS collapse time, we choose two different parametrized times to induce collapse. We also explore two initial magnetic field geometries in the torus, and evolve the systems until the outflows diminish significantly ($\sim\!\! 1\!\! - \!\!2\ \mathrm{s}$). We find bluer, faster outflows as compared to equivalent BH–torus systems, producing M ∼ 10 −3 M ⊙ of ejecta with Y e ≥ 0.25 and v ≥ 0.25c by the simulation end. Approximately half the outflows are launched in disc winds at times $t\lesssim 500 \ \mathrm{ms}$, with a broad distribution of electron fractions and velocities, depending on the initial condition. The remaining outflows are thermally driven, characterized by lower velocities and electron fractions. Nucleosynthesis with tracer particles shows patterns resembling solar abundances in all models. Although outflows from our simulations do not match those inferred from two-component modelling of the GW170817 kilonova, self-consistent multidimensional detailed kilonova models are required to determine whether our outflows can power the blue kilonova.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

The Northwest Geysers Enhanced Geothermal System Demonstration Project, The Geysers, California (Final Report)

An Enhanced Geothermal System (EGS) Demonstration Project was completed in the northwest portion of The Geysers geothermal field in Northern California. The goal was to demonstrate the feasibility of stimulating a deep, high temperature reservoir (HTR) with measured temperatures up to 750 °F (400 °C). We note that while some of the traditional Hot Dry Rock EGS concepts involves one or a few repeated hydraulically pressurized stimulations, our concept involves continuous injection over a periods of years at relatively low injection flows and at low pressures. Important measures of success are a sustained steam production rate and decreasing noncondensible gas (NCG) concentrations within the EGS Demonstration area. The project was completed in three phases and this Final Report characterizes the EGS reservoir and documents the successful stimulation of the HTR by the continuous injection of treated wastewater under only a gravity pressure head.

15 GEOTHERMAL ENERGY↗

DOE OSTI2021

Protoplanetary Disk Rings as Sites for Planetesimal Formation

Axisymmetric dust rings are a ubiquitous feature of young protoplanetary disks. These rings are likely caused by pressure bumps in the gas profile; a small bump can induce a traffic-jam-like pattern in the dust density, while a large bump may halt radial dust drift entirely. The resulting increase in dust concentration may trigger planetesimal formation by the streaming instability (SI), as the SI itself requires some initial concentration of dust. Here we present the first 3D simulations of planetesimal formation in the presence of a pressure bump modeled specifically after those seen by Atacama Large Millimeter/submillimeter Array. We place a pressure bump at the center of a large 3D shearing box, along with an initial solid-to-gas ratio of Z = 0.01, and we include both particle back-reaction and particle self-gravity. We consider millimeter-sized and centimeter-sized particles separately. For simulations with centimeter-sized particles, we find that even a small pressure bump leads to the formation of planetesimals via the SI; a pressure bump does not need to fully halt radial particle drift for the SI to become efficient. Furthermore, pure gravitational collapse via concentration in pressure bumps (such as would occur at sufficiently high concentrations and without the SI) is not responsible for planetesimal formation. For millimeter-sized particles, we find tentative evidence that planetesimal formation does not occur. If this result is confirmed at higher resolution, it could put strong constraints on where planetesimals can form. Ultimately, our results show that for centimeter-sized particles planetesimal formation in pressure bumps is extremely robust.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRSDate not supplied

Three-dimensional in Situ Observations of Polycrystalline Microstructure Evolution During Directional Solidification of Transparent Alloys Aboard the ISS and Quantitative Comparison with Numerical Modeling

Cellular/dendritic microstructures formed during solidification have a crucial influence on the mechanical properties of a wide range of structural alloys. By minimizing the amount of gravity-induced convection in the liquid, directional solidification experiments using transparent organic alloys conducted in the DECLIC-DSI onboard the International Space Station have provided unique 3D in situ observations of the spatiotemporal evolution of the solid-liquid interface during the formation of cellular and dendritic microstructures under purely diffusive growth conditions in polycrystalline samples containing several grains with a small misorientation with respect to the temperature gradient. Those observations have made it possible to perform benchmark quantitative comparisons with the predictions of state-of-the-art phase-field simulations of microstructure formation in 3D on experimentally relevant length and time scales. This talk will report quantitative comparisons between microgravity experiments and phase-field simulations in succinonitrile-camphor alloys of two different compositions that shed new light on the role of subgrain boundaries in the spatiotemporal evolution of the primary cellular spacing and the selection of dendritic array structures.

Kaihua Ji↗

NASA NTRSDate not supplied

Three-Dimensional In Situ Observations of Polycrystalline Microstructure Evolution During Directional Solidification of Transparent Alloys Aboard the ISS and Quantitative Comparison with Numerical Modeling

Cellular/dendritic microstructures formed during solidification have a crucial influence on the mechanical properties of a wide range of structural alloys. By minimizing the amount of gravity-induced convection in the liquid, directional solidification experiments using transparent organic alloys conducted in the DECLIC-DSI onboard the International Space Station have provided unique 3D in situ observations of the spatiotemporal evolution of the solid-liquid interface during the formation of cellular and dendritic microstructures under purely diffusive growth conditions. Those observations have made it possible to perform benchmark quantitative comparisons with the predictions of state-of-the-art phase-field simulations of microstructure formation in 3D on experimentally relevant length and time scales. This talk will report quantitative comparisons between microgravity experiments and phase-field simulations in succinonitrile-camphor alloys of two different compositions, which shed new light on the selection of dendritic array structures, and on roles of macroscopic curvatures and subgrain boundaries between grains with a small misorientation with respect to the temperature gradient in the spatiotemporal evolution of the primary cellular spacing.

Kaihua Ji↗

DOE OSTI2016

A NEW MODEL FOR MIXING BY DOUBLE-DIFFUSIVE CONVECTION (SEMI-CONVECTION). III. THERMAL AND COMPOSITIONAL TRANSPORT THROUGH NON-LAYERED ODDC

Oscillatory double-diffusive convection (ODDC; also known as semi-convection) refers to a type of double-diffusive instability that occurs in regions of planetary and stellar interiors that have a destabilizing thermal stratification and a stabilizing mean molecular weight stratification. In this series of papers, we use an extensive suite of three-dimensional (3D) numerical simulations to quantify the transport of heat and chemical species by ODDC. Rosenblum et al. first showed that ODDC can either spontaneously form layers that significantly enhance the transport of heat and chemical species compared to microscopic transport or remain in a state dominated by large-scale gravity waves, in which there is a more modest enhancement of the turbulent transport rates. Subsequent studies in this series focused on identifying under what conditions layers form and quantifying transport through layered systems. Here we proceed to characterize transport through systems that are unstable to ODDC, but do not undergo spontaneous layer formation. We measure the thermal and compositional fluxes in non-layered ODDC from both two-dimensional (2D) and 3D numerical simulations, and show that 3D simulations are well approximated by similar simulations in a 2D domain. We find that the turbulent mixing rate in this regime is weak and can, to a first-level approximation, be neglected. We conclude by summarizing the findings of papers I through III into a single prescription for transport systems unstable to ODDC.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2019

Design of Experiment for the Measurement of Aerosol Droplet Size Distribution of Temperature-Controlled Thermally Atomized Printed Electronic Inks

In-Space Manufacturing (ISM) centers around NASA’s growing need and ability to produce space technologies on demand in space. As the future of long term presence in space and deep space exploration approach, fundamental questions of our dependence on earth resupply to Low Earth Orbit (LEO) remain unanswered. ISM is leading various effort to evaluate the feasibility of producing essential spares and redundant parts on demand to enable a sustainable space-based supply chain model for part resupply. Among the parts and systems being considered, Avionics form the neural network of modern day aircraft and space vehicles providing a wealth of information ranging from Guidance, Navigation, and Control (GN&C) systems to on board Environmental Control and Life Support (ECLS) systems. Recent advances in the use of Aerosol Jet Technology to print Avionics components ranging from electrical traces on a circuit board to complex transistors and sensors raise the possibility of using such technology to reproduce or recreate electronic parts on demand with the help of custom electronics 3D printers. The challenge herein lies within the ability of such printers to generate and deposit an aerosol of electronic material utilizing processes independent of or enhanced by gravity to ensure controllably identical or improved behavior of the aerosol in an International Space Station (ISS) laboratory and on the ground. The behavior as well as the hazards and properties associated with such aerosols in a microgravity environment must be understood well in order to merit a feasible approach to utilizing them for manufacturing in space. In this report, we outline the experimental setup of a modified conventional vaping device to be used as the ideal gravity independent thermal atomization mechanism to generate aerosol. Our objective is to identify the ideal mass, density, and volume of our aerosolized droplets of ink to conclude that there exist a threshold of aerosolized ink droplet sizes that are indeed independent of the effects of gravity and remain stable after atomization. We use a Malvern Spraytec® Spray Particle Size Analyzer to perform real-time laser diffraction measurement of our ink droplets during atomization. The droplet size between conductive ink, dielectric ink and vegetable glycerin have been measured and contrasted. Furthermore, the mechanism of thermal atomization versus traditional pneumatic and ultrasonic atomization for operation in microgravity have been explored.

Wu, Chih-Hao↗

DOE OSTI2016

SOLAR PROMINENCES EMBEDDED IN FLUX ROPES: MORPHOLOGICAL FEATURES AND DYNAMICS FROM 3D MHD SIMULATIONS

The temporal evolution of a solar prominence inserted in a three-dimensional magnetic flux rope is investigated numerically. Using the model of Titov and Démoulin under the regime of weak twist, the cold and dense prominence counteracts gravity by modifying the initially force-free magnetic configuration. In some cases a quasi-stationary situation is achieved after the relaxation phase, characterized by the excitation of standing vertical oscillations. These oscillations show a strong attenuation with time produced by the mechanism of continuum damping due to the inhomogeneous transition between the prominence and solar corona. The characteristic period of the vertical oscillations does not depend strongly on the twist of the flux rope. Nonlinearity is responsible for triggering the Kelvin–Helmholtz instability associated with the vertical oscillations and that eventually produces horizontal structures. Contrary to other configurations in which the longitudinal axis of the prominence is permeated by a perpendicular magnetic field, like in unsheared arcades, the orientation of the prominence along the flux rope axis prevents the development of Rayleigh–Taylor instabilities and therefore the appearance of vertical structuring along this axis.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

Characterizing the pressure smoothing scale of the intergalactic medium

The thermal state of the intergalactic medium (IGM) at z < 6 constrains the nature and timing of cosmic reionization events, but its inference from the Lyα forest is degenerate with the 3D structure of the IGM on ∼100 kpc scales, where, analogous to the classical Jeans argument, the pressure of the T ≃ 10{sup 4} K gas supports it against gravity. We simulate the IGM using smoothed particle hydrodynamics, and find that, at z < 6, the gas density power spectrum does not exhibit the expected filtering scale cutoff, because dense gas in collapsed halos dominates the small-scale power masking pressure smoothing effects. We introduce a new statistic, the real-space Lyα flux, F{sub real}, which naturally suppresses dense gas, and is thus robust against the poorly understood physics of galaxy formation, revealing pressure smoothing in the diffuse IGM. The F{sub real} power spectrum is accurately described by a simple fitting function with cutoff at λ{sub F}, allowing us to rigorously quantify the pressure smoothing scale for the first time: we find λ{sub F} = 79 kpc (comoving) at z = 3 for our fiducial thermal model. This statistic has the added advantage that it directly relates to observations of correlated Lyα forest absorption in close quasar pairs, recently proposed as a method to measure the pressure smoothing scale. Our results enable one to quantify the pressure smoothing scale in simulations, and ask meaningful questions about its dependence on reionization and thermal history. Accordingly, the standard description of the IGM in terms of the amplitude T{sub 0} and slope γ of the temperature–density relation T=T{sub 0}(ρ/ ρ-bar ){sup γ−1} should be augmented with a third pressure smoothing scale parameter λ{sub F}.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Kinetic Simulations of the Kruskal–Schwarzschild Instability in Accelerating Striped Outflows: Dynamics and Energy Dissipation

Astrophysical relativistic outflows are launched as Poynting-flux dominated, yet the mechanism governing efficient magnetic dissipation, which powers the observed emission, is still poorly understood. We study magnetic energy dissipation in relativistic "striped" jets, which host current sheets separating magnetically dominated regions with opposite field polarity. The effective gravity force g in the rest frame of accelerating jets drives the Kruskal–Schwarzschild instability (KSI), a magnetic analog of the Rayleigh–Taylor instability. By means of 2D and 3D particle-in-cell simulations, we study the linear and nonlinear evolution of the KSI. The linear stage is well described by linear stability analysis. The nonlinear stages of the KSI generate thin (skin-depth-thick) current layers, with length comparable to the dominant KSI wavelength. There, the relativistic drift-kink mode and the tearing mode drive efficient magnetic dissipation. The dissipation rate can be cast as an increase in the effective width Δ eff of the dissipative region, which follows dΔ eff /dt ≃ 0.05 $\sqrt{Δ_{eff} g}$. Our results have important implications for the location of the dissipation region in gamma-ray burst and active galactic nuclei jets.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1994

Neonatal rat heart cells cultured in simulated microgravity

In vitro characteristics of cardiac cells cultured in simulated microgravity are reported. Tissue culture methods performed at unit gravity constrain cells to propagate, differentiate, and interact in a two dimensional (2D) plane. Neonatal rat cardiac cells in 2D culture organize predominantly as bundles of cardiomyocytes with the intervening areas filled by non-myocyte cell types. Such cardiac cell cultures respond predictably to the addition of exogenous compounds, and in many ways they represent an excellent in vitro model system. The gravity-induced 2D organization of the cells, however, does not accurately reflect the distribution of cells in the intact tissue. We have begun characterizations of a three-dimensional (3D) culturing system designed to mimic microgravity. The NASA designed High-Aspect-Ratio-Vessel (HARV) bioreactors provide a low shear environment which allows cells to be cultured in static suspension. HARV-3D cultures were prepared on microcarrier beads and compared to control-2D cultures using a combination of microscopic and biochemical techniques. Both systems were uniformly inoculated and medium exchanged at standard intervals. Cells in control cultures adhered to the polystyrene surface of the tissue culture dishes and exhibited typical 2D organization. Cells in cultured in HARV's adhered to microcarrier beads, the beads aggregated into defined clusters containing 8 to 15 beads per cluster, and the clusters exhibited distinct 3D layers: myocytes and fibroblasts appeared attached to the surfaces of beads and were overlaid by an outer cell type. In addition, cultures prepared in HARV's using alternative support matrices also displayed morphological formations not seen in control cultures. Generally, the cells prepared in HARV and control cultures were similar, however, the dramatic alterations in 3D organization recommend the HARV as an ideal vessel for the generation of tissue-like organizations of cardiac cells in simulated microgravity.

Akins, Robert E.↗

NASA NTRS1997

Neonatal rat heart cells cultured in simulated microgravity

In vitro characteristics of cardiac cells cultured in simulated microgravity are reported. Tissue culture methods performed at unit gravity constrain cells to propagate, differentiate, and interact in a two-dimensional (2D) plane. Neonatal rat cardiac cells in 2D culture organize predominantly as bundles of cardiomyocytes with the intervening areas filled by nonmyocyte cell types. Such cardiac cell cultures respond predictably to the addition of exogenous compounds, and in many ways they represent an excellent in vitro model system. The gravity-induced 2D organization of the cells, however, does not accurately reflect the distribution of cells in the intact tissue. We have begun characterizations of a three-dimensional (3D) culturing system designed to mimic microgravity. The NASA-designed High-Aspect Ratio Vessel (HARV) bioreactors provide a low shear environment that allows cells to be cultured in static suspension. HARV-3D cultures were prepared on microcarrier beads and compared to control-2D cultures using a combination of microscopic and biochemical techniques. Both systems were uniformly inoculated and medium exchanged at standard intervals. Cells in control cultures adhered to the polystyrene surface of the tissue culture dishes and exhibited typical 2D organization. Cells cultured in HARVs adhered to microcarrier beads, the beads aggregated into defined clusters containing 8 to 15 beads per cluster, and the clusters exhibited distinct 3D layers: myocytes and fibroblasts appeared attached to the surfaces of beads and were overlaid by an outer cell type. In addition, cultures prepared in HARVs using alternative support matrices also displayed morphological formations not seen in control cultures. Generally, the cells prepared in HARV and control cultures were similar; however, the dramatic alterations in 3D organization recommend the HARV as an ideal vessel for the generation of tissuelike organization of cardiac cells in vitro.

Non-NASA Center↗

NASA NTRS2002

Dynamics of Sheared Granular Materials

This work focuses on the properties of sheared granular materials near the jamming transition. The project currently involves two aspects. The first of these is an experiment that is a prototype for a planned ISS (International Space Station) flight. The second is discrete element simulations (DES) that can give insight into the behavior one might expect in a reduced-g environment. The experimental arrangement consists of an annular channel that contains the granular material. One surface, say the upper surface, rotates so as to shear the material contained in the annulus. The lower surface controls the mean density/mean stress on the sample through an actuator or other control system. A novel feature under development is the ability to 'thermalize' the layer, i.e. create a larger amount of random motion in the material, by using the actuating system to provide vibrations as well control the mean volume of the annulus. The stress states of the system are determined by transducers on the non-rotating wall. These measure both shear and normal components of the stress on different size scales. Here, the idea is to characterize the system as the density varies through values spanning dense almost solid to relatively mobile granular states. This transition regime encompasses the regime usually thought of as the glass transition, and/or the jamming transition. Motivation for this experiment springs from ideas of a granular glass transition, a related jamming transition, and from recent experiments. In particular, we note recent experiments carried out by our group to characterize this type of transition and also to demonstrate/ characterize fluctuations in slowly sheared systems. These experiments give key insights into what one might expect in near-zero g. In particular, they show that the compressibility of granular systems diverges at a transition or critical point. It is this divergence, coupled to gravity, that makes it extremely difficult if not impossible to characterize the transition region in an earth-bound experiment. In the DE modeling, we analyze dynamics of a sheared granular system in Couette geometry in two (2D) and three (3D) space dimensions. Here, the idea is to both better understand what we might encounter in a reduced-g environment, and at a deeper level to deduce the physics of sheared systems in a density regime that has not been addressed by past experiments or simulations. One aspect of the simulations addresses sheared 2D system in zero-g environment. For low volume fractions, the expected dynamics of this type of system is relatively well understood. However, as the volume fraction is increased, the system undergoes a phase transition, as explained above. The DES concentrate on the evolution of the system as the solid volume fraction is slowly increased, and in particular on the behavior of very dense systems. For these configurations, the simulations show that polydispersity of the sheared particles is a crucial factor that determines the system response. Figures 1 and 2 below, that present the total force on each grain, show that even relatively small (10 %) nonuniformity of the size of the grains (expected in typical experiments) may lead to significant modifications of the system properties, such as velocity profiles, temperature, force propagation, and formation shear bands. The simulations are extended in a few other directions, in order to provide additional insight to the experimental system analyzed above. In one direction, both gravity, and driving due to vibrations are included. These simulations allow for predictions on the driving regime that is required in the experiments in order to analyze the jamming transition. Furthermore, direct comparison of experiments and DES will allow for verification of the modeling assumptions. We have also extended our modeling efforts to 3D. The (preliminary) results of these simulations of an annular system in zero-g environment will conclude the presentation.

Kondic, Lou↗

DOE OSTI2016

Linking pore-scale and basin-scale effects on diffusive methane transport in hydrate bearing environments through multi-scale reservoir simulations

We explore the gas hydrate-generating capacity of diffusive methane transport induced by solubility gradients due to pore size contrasts in lithologically heterogeneous marine sediments. Through the use of 1D, 2D, and 3D reactive transport simulations, we investigate scale-dependent processes in diffusion-dominated gas hydrate systems. These simulations all track a sand body, or series of sands, surrounded by clays as they are buried through the gas hydrate stability zone. Methane is sourced by microbial methanogenesis in the clays surrounding the sand layers. In 1D, simulations performed in a Lagrangian reference frame demonstrate that gas hydrate in thin sands (3.6 m thick) can occur in high saturations (upward of 70%) at the edges of sand bodies within the upper 400 meters below the seafloor. Diffusion of methane toward the center of the sand layer depends on the concentration gradient within the sand: broader sand pore size distributions with smaller median pore sizes enhance diffusive action toward the sand’s center. Incorporating downhole log- and laboratory-derived sand pore size distributions, gas hydrate saturations in the center of the sand can reach 20% of the hydrate saturations at the sand’s edges. Furthermore, we show that hydrate-free zones exist immediately above and below the sand and are approximately 5 m thick, depending on the sand-clay solubility contrast. A moving reference frame is also adopted in 2D, and the angle of gravity is rotated relative to the grid system to simulate a dipping sand layer. This is important to minimize diffusive edge effects or numerical diffusion that might be associated with a dipping sand in an Eulerian grid system oriented orthogonal to gravity. Two-dimensional simulations demonstrate the tendency for gas hydrate to accumulate downdip in a sand body because of greater methane transport at depth due to larger sand-clay solubility contrasts. In 3D, basin-scale simulations illuminate how convergent sand layers in a multilayered system can compete for diffusion from clays between them, resulting in relatively low hydrate saturations. All simulations suggest that when hydrate present in clays dissociates with burial, the additional dissolved methane is soaked up by nearby sands preserving high hydrate saturations.

Nole, Michael↗