SEARCH · Engineering Papers

Results for “gas diffusion layers”

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 235 records · Page 13

All sources NASA NTRS DOE OSTI

NASA NTRS1972

Mercury contamination study for flight system safety

The effects and prevention of possible mercury pollution from the failure of solar electric propulsion spacecraft using mercury propellant were studied from tankage loading of post launch trajector injection. During preflight operations and initial flight mode there is little danger of mercury pollution if proper safety precautions are taken. Any spillage on the loading, mating, transportation, or launch pad areas is obvious and can be removed by vacuum cleaning soil and chemical fixing. Mercury spilled on Cape Kennedy ground soil will be chemically complexed and retained by the sandstone subsoil. A cover layer of sand or gravel on spilled mercury which has settled to the bottom of a water body adjacent to the system operation will control and eliminate the formation of toxic organic mercurials. Mercury released into the earth's atmosphere through leakage of a fireball will be diffused to low concentration levels. However, gas phase reactions of mercury with ozone could cause a local ozone depletion and result in serious ecological hazards.

Gorzynski, C. S., Jr.↗

DOE OSTI2023

Simultaneous measurements of kinetic and scalar energy spectrum time evolution in the Richtmyer–Meshkov instability upon reshock

The Richtmyer–Meshkov instability (Richtmyer, Commun. Pure Appl. Maths , vol. 13, issue 2, 1960, pp. 297–319; Meshkov, Fluid Dyn. , vol. 4, issue 5, 1972, pp. 101–104) of a twice-shocked gas interface is studied using both high spatial resolution single-shot (SS) and lower spatial resolution, time-resolved, high-speed (HS) simultaneous planar laser-induced fluorescence and particle image velocimetry in the Wisconsin Shock Tube Laboratory's vertical shock tube. The initial condition (IC) is a shear layer with broadband diffuse perturbations at the interface between a helium–acetone mixture and argon. This IC is accelerated by a shock of nominal strength Mach number $M = 1.75$ , and then accelerated again by the transmitted shock that reflects off the end wall of the tube. An ensemble of experiments is analysed after reshock while the interface mixing width grows linearly with time. The kinetic and scalar energy spectra and the terms of their evolution equation are calculated and compared between SS and HS experiments. The inertial range scaling of the scalar power spectrum is found to follow Gibson's relation (Gibson, Phys. Fluids , vol. 11, issue 11, 1968, pp. 2316–2327) as a function of Schmidt number when the effective turbulent Schmidt number is used in place of the material Schmidt number that controls equilibrium scaling. Further, the spatially integrated scalar flux follows similar behaviour observed for the kinetic energy in large eddy simulation studies by Zeng et al. ( Phys. Fluids , vol. 30, issue 6, 2018, 064106) while the spatially varying scalar flux exhibits back scatter along the centre of the mixing layer and forward energy transfer in the spike and bubble regions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1975

Flow visualization of discrete hole film cooling for gas turbine applications

Film injection from discrete holes in a three row staggered array with 5-diameter spacing is studied for three different hole angles: (1) normal, (2) slanted 30 deg to the surface in the direction of the mainstream, and (3) slanted 30 deg to the surface and 45 deg laterally to the mainstream. The boundary layer thickness-to-hole diameter ratio and Reynolds number are typical of gas turbine film cooling applications. Two different injection locations are studied to evaluate the effect of boundary layer thickness on film penetration and mixing. Detailed streaklines showing the turbulent motion of the injected air are obtained by photographing very small neutrally buoyant helium filled 'soap' bubbles which follow the flow field. Unlike smoke, which diffuses rapidly in the high turbulent mixing region associated with discrete hole blowing, the bubble streaklines passing downstream injection locations are clearly identifiable and can be traced back to their origin. Visualization of surface temperature patterns obtained from infrared photographs of a similar film cooled surface are also included.

Colladay, R. S.↗

DOE OSTI2024

Materials Science of the Interstitial Doping Process

Particle accelerators are an increasingly important tool for frontier science. Growing initial and operating costs are a significant barrier for upgrades and for new machines. While everything matters, the major cost contributor is the SRF cavities and their ancillary facilities (e.g., cryoplant). Accordingly, the accelerator science community devotes much R&D effort to improving their energy efficiency (increased Q o ) and gradient (E acc ). While improved gradient is at the forefront for certain machines (ILC), improved quality factor has broader impact, benefitting all SRF applications. An important opportunity for accelerator science and technology to move forward arose in the course of building the LCLS-II, the second generation Linac Coherent Light Source at SLAC. At more or less the same time, researchers at Fermilab discovered that introducing a small amount of nitrogen to the niobium surface could improve the mid-range quality factor as much as three-fold. A firm resolution of how nitrogen confers its benefit attracts much current research interest. The key elements of the “nitrogen doping” process were vacuum bake at 800 °C, brief exposure to mTorr of nitrogen at several hundred degrees followed by electropolish (EP) to remove several microns from the surface to eliminate unwanted nitrides: While the process as a whole was novel, the comprising unit operations are familiar to the accelerator community. It was judged reasonable to adopt it as a cost-reduction technology for LCLS-II. Researchers carried out a program of varying process parameters and measuring performance in single-cell cavities, leading to a consensus stable protocol for the project. The transition to vendor fabrication and multiple niobium sources has presented unforeseen challenges of performance variation evidently not connected to anything that could be incorporated in a purchase specification. Moving beyond the high temperature N-doping process above, researchers reported a simplified process consisting entirely of tens of hours anneal in a N atmosphere at low temperatures (~120°C – 160°C) after 800 °C UHV bake. These processes typically yield a few-nm doped layer while the high temperature process yields at least a few to many micron doped layer. Even more recently, oxygen has been used to dope instead of nitrogen which leaves a few-µm O-alloyed layer upon vacuum annealing for 300 °C for ~3 hours. The investigation of these materials is just beginning, but the process simplification they may offer is surely attractive. The very low quantity of material that appears to be significant in the “infusion” process indicates the need for very careful control of gas species available for diffusion into the surface during low temperature treatment, both for process control and research to characterize the underlying dynamics. Oxygen alloying offers the further opportunity to utilize the decomposition of the surface native oxide as the dopant source. It is necessary to understand and (thus) manage this process. We have been supported by the Department of Energy Offices of High Energy Physics and Nuclear Physics to pursue this goal.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Enhancing reactant selectivity for Ni/Mg reforming catalysts using silicalite-1 shells: A modeling study

The ability of a zeolite shell to enhance the selective conversion of hydrocarbons through diffusional limitations was investigated using a multi-scale model of a fixed-bed reactor. The impact of shell thickness and molecule/pore size on the catalytic performance of silicalite-1 zeolite encapsulated nickel catalyst pellets for steam reforming of C1-C7 hydrocarbons is reported. A reaction–diffusion model using kinetic expressions established in literature was employed. The model was verified through comparison with reported experimental results for steam reforming data over a temperature range of 748 – 1113 K and pressure of 1 – 10 bar. Comparisons are also made against experimental data for steam reforming in the presence of a zeolite shell. Evaluation of the Weisz-Prater criterion for both the core and encapsulated catalyst confirmed mass transfer limitation induced by the utilization of a zeolite shell. The model was used to suggest an optimal thickness that balances diffusional limitations imposed by the zeolite layer on methane versus that of the heavier hydrocarbons. The optimum thickness varied as a function of hydrocarbon size and shape which determined the diffusion rates. For toluene and heptane, a 50 nm thick shell was sufficient to wholly prevent reaction. Hydrocarbons like propane and butane required a shell 7.5 and 5 μm thick. Increasing the gas-hourly-space-velocity from 10,000 to 60,000 h –1 caused a decrease in the optimum shell thickness. Furthermore, this approach can be modified for application to other mixed hydrocarbon systems to predict optimal catalyst design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Amorphous fluorinated copolymer gas separation membranes

Membranes having a permselective active layer of a copolymerized perfluorinated monomer and an non-fluorinated alkylvinylester monomer demonstrate superior selective permeability performance for separating gas mixtures compared to membranes of exclusively perfluorinated polymers. Preferred active layer compositions are copolymers of perfluoro-2,2-dimethyl-1,3 dioxole (PDD) copolymerized with an alkylvinyl ester such as vinyl acetate, and vinyl pivalate, and with alkylvinyl esters that are substantially hydrolyzed to provide copolymerized vinyl alcohol functionality. The membranes can have a thin, high diffusion rate, “gutter layer” of a fluorinated polymer highly permeable to nitrogen positioned between the active layer and a porous support layer. A novel copolymer effective in selectively permeable membranes is a copolymer of PDD and an alkylvinyl ester compound having the formula H2C═CHOC(O)R1 in which R1 is a linear or branched alkyl group of from 2 to 5 carbon atoms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2005

Linear Stability Analysis of Gravitational Effects on a Low-Density Gas Jet Injected into a High-Density Medium

The objective of this study was to determine the effects of buoyancy on the absolute instability of low-density gas jets injected into high-density gas mediums. Most of the existing analyses of low-density gas jets injected into a high-density ambient have been carried out neglecting effects of gravity. In order to investigate the influence of gravity on the near-injector development of the flow, a spatio-temporal stability analysis of a low-density round jet injected into a high-density ambient gas was performed. The flow was assumed to be isothermal and locally parallel; viscous and diffusive effects were ignored. The variables were represented as the sum of the mean value and a normal-mode small disturbance. An ordinary differential equation governing the amplitude of the pressure disturbance was derived. The velocity and density profiles in the shear layer, and the Froude number (signifying the effects of gravity) were the three important parameters in this equation. Together with the boundary conditions, an eigenvalue problem was formulated. Assuming that the velocity and density profiles in the shear layer to be represented by hyperbolic tangent functions, the eigenvalue problem was solved for various values of Froude number. The Briggs-Bers criterion was combined with the spatio-temporal stability analysis to determine the nature of the absolute instability of the jet whether absolutely or convectively unstable. The roles of the density ratio, Froude number, Schmidt number, and the lateral shift between the density and velocity profiles on the absolute instability of the jet were determined. Comparisons of the results with previous experimental studies show good agreement when the effects of these variables are combined together. Thus, the combination of these variables determines how absolutely unstable the jet will be.

Lawson, Anthony L.↗

NASA NTRS1987

Composite battery separator

A composite battery separator comprises a support element (10) having an open pore structure such as a ribbed lattice and at least one liquid permeable sheet (20,22) to distribute the compressive force evenly onto the surfaces of the layers (24, 26) of negative active material and positive active material. In a non-flooded battery cell the compressible, porous material (18), such as a glass mat which absorbs the electrolyte, is compressed into a major portion of the pores or openings (16) in the support element. The unfilled pores in the material (18) form a gas diffusion path as the channels (41) formed between adjacent ribs in the lattice element (30,36). Facing two lattice elements (30, 31) with acute angled cross-ribs (34, 38) facing each other prevents the elements from interlocking and distorting a porous, separator (42) disposed between the lattice elements.

Edwards, Dean B.↗

NASA NTRS2015

Remote Sensing Measurements of the CO2 Mixing Ratio in the Planetary Boundary Layer Using Cloud Slicing with Airborne Lidar

We have measured the CO2 volume mixing ratio (VMR) within the planetary boundary layer (PBL) using cloud slicing with an airborne pulsed integrated path differential absorption (IPDA) lidar from flight altitudes of up to 13 km. During a flight over Iowa in summer 2011, simultaneous measurement of the optical range and CO2 absorption to clouds and the ground were made using time-resolved detection of pulse echoes from each scattering surface.We determined the CO2 absorption in the PBL by differencing the two lidar-measured absorption line shapes, one to a broken shallow cumulus cloud layer located at the top of the PBL and the other to the ground. Solving for the CO2 VMR in the PBL and that of the free troposphere, we measured a ≈15 ppm (4%) drawdown in the PBL. Both CO2 VMRs were within ≈3 ppm of in situ CO2 profile measurements. We have also demonstrated cloud slicing using scatter from thin, diffuse cirrus clouds and cumulus clouds, which allowed solving for the CO2 VMR for three vertical layers. The technique and retrieval algorithm are applicable to a space-based lidar instrument as well as to lidar IPDA measurements of other trace gases. Thus, lidar cloud slicing also offers promise toward space-based remote sensing of vertical trace gas profiles in the atmosphere using a variety of clouds.

VMR↗

DOE OSTI2024

"Source Term Modeling for Advanced Gas Micro-Reactors"

Maintaining the safety of the public, environment, and operating personnel is the most important factor in designing, operating, maintaining, and decommissioning nuclear reactors. In recent years, there has been a growing interest in the development of micro-reactors employing TRi-structural ISOtropic (TRISO)-coated particle fuel. In gas reactors, TRISO fuel plays an important role in the safety case for high temperature reactors because of the fission product retention properties of the fuel. This ability enables the use of a functional containment strategy for the reactor where multiple barriers are used to prevent fission product release to the environment. Part of the safety analysis of these advanced reactors is the assessment of radionuclide releases under normal and accident conditions through the multiple credited safety barriers. Using conservative assumptions, a mechanistic analysis can be performed to quantify these releases that combines the probabilistic assessment of failure with analytic solutions to radionuclide transport equations. Source term modeling for TRISO fuel has been performed for previous reactor designs; however, these models are outdated, in many cases proprietary, and need updates to be applied to the current state of TRISO fuel technology and alternative gas reactor core configurations [1]. Currently, the only publicly available source term assessment for gas reactors is an expert-based Monte Carlo simulation based on the effectiveness of the fuel kernel, coating layers, and graphite block in a modular high temperature gas reactor [2]. Thus, there is a need to develop a simple, versatile, and mechanistic model of fission product release and transport in gas reactor cores that could be applied to a variety of reactors through user inputs and reactor-specific radionuclide inventories. The release is calculated by the diffusion of the key safety important fission products through the kernel, silicon carbide (SiC), graphite for both intact and defective TRISO particles based on fuel and graphite temperatures in the reactor under normal operation. These releases from the fuel enter the coolant where they can plate-out on cooler surfaces. A clean-up model is included for designs with a coolant purification system to remove fission gases. This initial distribution of fission products in the reactor serves as an initial condition for potential releases under postulated accident conditions. The model then can calculate the fission product release for any transient temperature profile and fission product releases can then be used to assess radiological dose to the workers and the public using conventional dose tools. Data on the diffusion of fission products is based on historic German TRISO experiments and the more current Department of Energy (DOE) Advanced Gas Reactor (AGR) TRISO fuel development program. The model is coded in python with inputs and outputs in excel spreadsheets, as well as python plotting utilities to aid in the interpretation of the results. References: [1] INL, NGNP Mechanistic Source Term White Paper, INL-10-17997, July 2010. [2] David A. Petti, Richard R. Hobbins, Peter Lowry, Hans Gougar, “Representative Source Terms and The Influence of Reactor Attributes on Functional Containment in Modular High Temperature Gas-cooled Reactors,” Nuclear Technology, Vol. 184, p. 181-197, Nov. 2013.

07 ISOTOPE AND RADIATION SOURCES↗

DOE OSTI2023

Gas nitriding behavior of refractory metals and implications for multi-principal element alloy design

Multi-principal element alloys (MPEAs) comprise a large, flexible compositional space that enables tuning of their chemistry, structure, and properties. Here, to facilitate the development of nitriding-based surface-enhancement strategies that harness a broad compositional space, this study examined the gas nitriding behavior of Hf, Mo, Nb, Ta, Ti, and Zr as a function of time, temperature (750 and 1000 °C), and nitriding potential (i.e. ammonia-to-hydrogen ratio). These metals were selected because they have a strong driving force to form nitrides, and appear in many promising refractory MPEA compositions. The nitriding temperatures were selected based on the phase transformation temperature of Ti and Zr, and the nitriding potentials were chosen such that all elements are expected to form nitrides. Mass gain measurements indicate that all six elements follow parabolic kinetics. The microstructure observations and quantitative microchemical analysis show formation of dense and well-adhered compound layers for Mo, Nb, and Ta. Thick diffusion zones appear in Hf, Ta, Ti, and Zr, and diffusion coefficients were fit to the composition profiles. Partial delamination of the compound layer occurred for Ti and Zr. Peak hardness values above 30 GPa are obtained in the dense compound layers, and the solute hardening of the underlying alloy is correlated with the nitrogen content. The results provide insight into the dynamics of nitride compound formation relative to interstitial dissolution of nitrogen, and are discussed in the context of MPEA composition and processing design.

36 MATERIALS SCIENCE↗

DOE OSTI2020

In Situ Oxidation Studies of High-Entropy Alloy Nanoparticles

Although high-entropy alloys (HEAs) have shown tremendous potential for elevated temperature, anticorrosion, and catalysis applications, little is known on how HEA materials behave under complex service environments. Herein, we studied the high-temperature oxidation behavior of Fe 0.28 Co 0.21 Ni 0.20 Cu 0.08 Pt 0.23 HEA nanoparticles (NPs) in an atmospheric pressure dry air environment by in situ gas-cell transmission electron microscopy. It is found that the oxidation of HEA NPs is governed by Kirkendall effects with logarithmic oxidation rates rather than parabolic as predicted by Wagner's theory. Further, the HEA NPs are found to oxidize at a significantly slower rate compared to monometallic NPs. The outward diffusion of transition metals and formation of disordered oxide layer are observed in real time and confirmed through analytical energy dispersive spectroscopy, and electron energy loss spectroscopy characterizations. Localized ordered lattices are identified in the oxide, suggesting the formation of Fe 2 O 3 , CoO, NiO, and CuO crystallites in an overall disordered matrix. Hybrid Monte Carlo and molecular dynamics simulations based on first-principles energies and forces support these findings and show that the oxidation drives surface segregation of Fe, Co, Ni, and Cu, while Pt stays in the core region. The present work offers key insights into how HEA NPs behave under high-temperature oxidizing environment and sheds light on future design of highly stable alloys under complex service conditions.

36 MATERIALS SCIENCE↗

NASA NTRS1998

SiC Recession Due to SiO2 Scale Volatility Under Combustion Conditions: Thermodynamics and Gaseous Diffusion Model - Part 2

In combustion environments, volatilization of SiO2 to Si-O-H(g) species is a critical issue. Available thermochemical data for Si-O-H(g) species were used to calculate boundary layer controlled fluxes from SiO2. Calculated fluxes were compared to volatilization rates Of SiO2 scales grown on SiC which were measured in Part 1 of this paper. Calculated volatilization rates were also compared to those measured in synthetic combustion gas furnace tests. Probable vapor species were identified in both fuel-lean and fuel-rich combustion environments based on the observed pressure, temperature and velocity dependencies as well as the magnitude of the volatility rate. Water vapor is responsible for the degradation of SiO2 in the fuel-lean environment. Silica volatility in fuel-lean combustion environments is attributed primarily to the formation of Si(OH)4(g) with a small contribution of SiO(OH)2(g).

Opila, Elizabeth J.↗

DOE OSTI2025

Multiscale, mechanistic modeling of cesium transport in silicon carbide for TRISO fuel performance prediction

Understanding cesium (Cs) transport in TRistructural ISOtropic (TRISO) particle fuel is crucial for predicting fission product release in high-temperature reactors. However, current challenges include significant scatter in diffusivity data and unexplained temperature-dependent diffusion regimes in the silicon carbide layer. This study addresses these challenges by developing a multiscale, mechanistic Cs transport model integrating atomistic simulations and phase field modeling. Our model quantifies temperature and grain size effects on Cs diffusivity, attributing experimentally observed regimes to a transition from bulk-dominated diffusivity at high temperatures to grain boundary-dominated diffusivity at lower temperatures. The model, validated against diffusion measurements and advanced gas reactor (AGR)-1 and AGR-2 post-irradiation fission product release data, enhances the predictive capability of the BISON fuel performance code. This study advances our understanding of Cs release from TRISO particles and its dependence on temperature and silicon carbide grain size, with implications for the safety and efficiency of high-temperature nuclear reactors.

BISON↗

NASA NTRS1984

Analytic modeling of a spray diffusion flame

A detailed model for a spray diffusion flame is described. The model is based on the boundary layer form of the equations of motion, with droplet transport accounted for using a discretized droplet size distribution function. Interphase transport of mass and energy are accounted for, with a flame-sheet model used to describe the combustion process on a droplet scale. Near dynamic equilibrium is assumed for the description of droplet transport; droplets can diffuse relative to the gas phase. Gas-phase mixing is accounted for using a two-equation turbulence model; buoyancy effects are included, with a temperature fluctuation equation used to account for buoyancy effects on turbulence structure. Thermal radiation from gas-phase CO2 and H2O is included. Gas-phase chemical kinetics are modeled using a 20-reaction, 10-species version of the advanced quasi-global chemical kinetics formulation. Results are compared with data for a vaporizing Freon spray and a pentane spray flame. It is shown that the computational approach provides a reasonably valid picture of the overall development of a spray diffusion flame, and, furthermore, provides a useful tool for the parametric examination of the spray combustion process.

Harsha, P. T.↗

NASA NTRS1990

Characteristics of the Martian atmosphere surface layer

Elements of various terrestrial boundary layer models are extended to Mars in order to estimate sensible heat, latent heat, and momentum fluxes within the Martian atmospheric surface ('constant flux') layer. The atmospheric surface layer consists of an interfacial sublayer immediately adjacent to the ground and an overlying fully turbulent surface sublayer where wind-shear production of turbulence dominates buoyancy production. Within the interfacial sublayer, sensible and latent heat are transported by non-steady molecular diffusion into small-scale eddies which intermittently burst through this zone. Both the thickness of the interfacial sublayer and the characteristics of the turbulent eddies penetrating through it depend on whether airflow is aerodynamically smooth or aerodynamically rough, as determined by the Roughness Reynold's number. Within the overlying surface sublayer, similarity theory can be used to express the mean vertical windspeed, temperature, and water vapor profiles in terms of a single parameter, the Monin-Obukhov stability parameter. To estimate the molecular viscosity and thermal conductivity of a CO2-H2O gas mixture under Martian conditions, parameterizations were developed using data from the TPRC Data Series and the first-order Chapman-Cowling expressions; the required collision integrals were approximated using the Lenard-Jones potential. Parameterizations for specific heat and binary diffusivity were also determined. The Brutsart model for sensible and latent heat transport within the interfacial sublayer for both aerodynamically smooth and rough airflow was experimentally tested under similar conditions, validating its application to Martian conditions. For the surface sublayer, the definition of the Monin-Obukhov length was modified to properly account for the buoyancy forces arising from water vapor gradients in the Martian atmospheric boundary layer. It was found that under most Martian conditions, the interfacial and surface sublayers offer roughly comparable resistance to sensible heat and water vapor transport and are thus both important in determining the associated fluxes.

Clow, G. D.↗

NASA NTRS1996

NASA flow fields analysis

The objectives of the present research are to improve design capabilities for low thrust rocket engines through understanding the detailed mixing and combustion processes in a representative combustor. Of particular interest is a small gaseous hydrogen-oxygen thruster which is considered as a coordinated part of an on-going experimental program at NASA LERC. Detailed computational modeling involves the solution of both the two- and three-dimensional Navier Stokes equations, coupled with chemical reactions and the species diffusion equations. Computations of interest include both steady state and time-accurate flowfields and are obtained by means of LU approximate factorization in time and flux split upwinding differencing in space. The emphasis in the research is focused on using numerical analysis to understand detailed combustor flowfields, including the shear layer dynamics created between fuel film cooling and the core gas in the vicinity on the nearby combustor wall; the integrity and effectiveness of the coolant film; and three-dimensional fuel and oxidizer jet injection/mixing/combustion characteristics in the primary combustor along with their joint impacts on global engine performance.

Merkle, Charles L.↗

DOE OSTI2022

Surface Interactions and Nanoconfinement of Methane and Methane plus CO2 Revealed by High-Pressure Magic Angle Spinning NMR Spectroscopy and Molecular Dynamics

This study explores the fundamental, molecular- to microscopic-level behavior of methane gas confined into nanoporous silica proxies with different pore diameters and surface-to-volume (S/V) ratios. Surfaces and pore walls of nanoporous silica matrices are decorated with hydroxyl (-OH) groups, resembling natural heterogeneity. High-pressure MAS NMR was utilized to characterize the interactions between methane and the engineered nanoporous silica proxies under various temperature and pressure regimes. There was a change in the chemical shift position of confined methane slightly in the mixtures with nanoporous silica up to 393 K, as shown by high-pressure 13C-NMR. The 13C-NMR chemical shift of methane was changed by pressure, explained by the densification of methane inside the nanoporous silica materials. The influence of pore diameter and S/V of the nanoporous silica materials on the behaviors and dynamics of methane were studied. The presence of CO2 in mixtures of silica and methane needs analysis with caution because CO2 in a supercritical state and gaseous CO2 change the original structure of nanoporous silica and change surface area and pore volume. According to simulation, the picosecond scale dynamics of methane confined in larger pores of amorphous silica is faster. In the 4 nm pore, the diffusivity obtained from MD simulations in the pore with a higher S/V ratio is slower due to the trapping of methane molecules in adsorbed layers close to the corrugated pore surface. In contrast, relaxation measured with NMR for smaller pores (higher S/V) exhibits larger T1, indicating slower relaxation.

03 NATURAL GAS↗