RESUMEN
A modified many-body dissipative particle dynamics (mDPD) model is rigorously calibrated to achieve realistic fluid-fluid/solid interphase properties and applied for mesoscale flow simulations to elucidate the transport mechanisms of heptane liquid and water, respectively, through pore networks formed by packed silica nanoparticles with a uniform diameter of 30 nm. Two million CPU core hours were used to complete the simulation studies. Results show reduction of permeability by 54-64% in heptane flow and by 88-91% in water flow, respectively, compared to the Kozeny-Carman equation. In these nanopores, a large portion of the fluids are in the near-wall regions and thus not mobile due to the confinement effect, resulting in reduced hydraulic conductivity. Moreover, intense oscillations in the calculated flow velocities also indicate the confinement effect that contests the external driven force to flow. The generic form of Darcy's law is considered valid for flow through homogeneous nanopore networks, while permeability depends collectively on pore size and surface wettability. This fluid-permeability dependency is unique to flow in nanopores. In addition, potential dependence of permeability on pore connectivity is observed when the porosity remains the same in different core specimens.
RESUMEN
The pressure-driven flow of long-chain hydrocarbons in nanosized pores is important in energy, environmental, biological, and pharmaceutical applications. This paper examines the flow of hexane, heptane, and decane in carbon nanotubes (CNTs) of pore diameters 1-8 nm using molecular dynamic simulations. Enhancement of water flow in CNTs in comparison to rates predicted by continuum models has been well established in the literature. Our work was intended to observe if molecular dynamic simulations of hydrocarbon flow in CNTs produced similar enhancements. We used the OPLS-AA force field to simulate the hydrocarbons and the CNTs. Our simulations predicted the bulk densities of the hydrocarbons to be within 3% of the literature values. Molecular sizes and shapes of the hydrocarbon molecules compared to the pore size create interesting density patterns for smaller sized CNTs. We observed moderate flow enhancements for all the hydrocarbons (1-100) flowing through small-sized CNTs. For very small CNTs the larger hydrocarbons were forced to flow in a cork-screw fashion. As a result of this flow orientation, the larger molecules flowed as effectively (similar enhancements) as the smaller hydrocarbons.
RESUMEN
Advances in imaging have made it possible to view nanometer and sub-nanometer structures that are either synthesized or that occur naturally. It is believed that fluid dynamic and thermodynamic behavior differ significantly at these scales from the bulk. From a materials perspective, it is important to be able to create complex structures at the nanometer scale, reproducibly, so that the fluid behavior may be studied. New advances in nanoscale-resolution 3D-printing offer opportunities to achieve this goal. In particular, additive manufacturing with two-photon polymerization allows creation of intricate structures. Using this technology, a creation of the first nano-3D-printed digital (shale) rock is reported. In this paper, focused ion beam-scanning electron microscopy (FIB-SEM) nano-tomography image dataset was used to reconstruct a high-resolution digital rock 3D model of a Marcellus Shale rock sample. Porosity of this 3D model has been characterized and its connected/effective pore system has been extracted and nano-3D-printed. The workflow of creating this novel nano-3D-printed digital rock 3D model is described in this paper.
RESUMEN
The phase behavior of hydrocarbon fluids confined in porous media has been reported to deviate significantly from that in the bulk environment due to the existence of sub-10 nm pores. Though experiments and simulations have measured the bubble/dew points and sorption isotherms of hydrocarbons confined in both natural and synthetic nanopores, the confinement effects in terms of the strength of fluid-pore interactions tuned by surface wettability and chemistry have received comparably less discussion. More importantly, the underlying physics of confinement-induced phenomena remain obfuscated. In this work, we studied the phase behavior and capillary condensation of n-hexane to understand the effects of confinement at the molecular level. To systematically investigate the pore effects, we constructed two types of wall confinements; one is a structureless virtual wall described by the Steele potential and the other one is an all-atom amorphous silica structure with surface modified by hydroxyl groups. Our numerical results demonstrated the importance of fluid-pore interaction, pore size, and pore morphology effects in mediating the pressure-volume-temperature (PVT) properties of hydrocarbons. The most remarkable finding of this work was that the saturation pressure predicted from the van der Waals-type adsorption isothermal loop could be elevated or suppressed relative to the bulk phase, as illustrated in the graphical abstract. As the surface energy (i.e., fluid-pore interaction) decreased, the isothermal vapor pressure increased, indicating a greater preference for the fluid to exist in the vapor state. Sufficient reduction of the fluid-pore interactions could even elevate the vapor pressure above that of the bulk fluid.
RESUMEN
Nanoconfined water plays a pivotal role in a vast number of fields ranging from biological and materials sciences to catalysis, nanofluidics and geochemistry. Here, we report the freezing and melting behavior of water (D2O) nanoconfined in architected silica-based matrices including Vycor glass and mesoporous silica SBA-15 and SBA-16 with pore diameters ranging between 4-15 nm, which are investigated using differential scanning calorimetry and 2H nuclear magnetic resonance spectroscopy. The results provide compelling evidence that the extreme dynamical heterogeneity of water molecules is preserved over distances as small as a few angstroms. Solidification progresses in a layer-by-layer fashion with a coexistence of liquid-like and solid-like dynamical fraction at all temperatures during the transition process. The previously reported fragile-to-strong dynamic transition in nanoconfined water is argued to be a direct consequence of the layer-by-layer solidification.
RESUMEN
Porosity and permeability are the key factors in assessing the hydrocarbon productivity of unconventional (shale) reservoirs, which are complex in nature due to their heterogeneous mineralogy and poorly connected nano- and micro-pore systems. Experimental efforts to measure these petrophysical properties posse many limitations, because they often take weeks to complete and are difficult to reproduce. Alternatively, numerical simulations can be conducted in digital rock 3D models reconstructed from image datasets acquired via e.g., nanoscale-resolution focused ion beam-scanning electron microscopy (FIB-SEM) nano-tomography. In this study, impact of reservoir confinement (stress) on porosity and permeability of shales was investigated using two digital rock 3D models, which represented nanoporous organic/mineral microstructure of the Marcellus Shale. Five stress scenarios were simulated for different depths (2,000-6,000 feet) within the production interval of a typical oil/gas reservoir within the Marcellus Shale play. Porosity and permeability of the pre- and post-compression digital rock 3D models were calculated and compared. A minimal effect of stress on porosity and permeability was observed in both 3D models. These results have direct implications in determining the oil-/gas-in-place and assessing the production potential of a shale reservoir under various stress conditions.
RESUMEN
Type I kerogen was isolated from Green River Shale and characterized using SEM, TGA, DSC, and nitrogen adsorption. The swelling behavior of this kerogen with decane was analyzed using traditional test-tube swelling experiment and Dynamic Light Scattering. The TGA and DSC were used to analyze the thermodynamic behavior of decane that was sorbed in the kerogen and show that kerogen suppresses the boiling point of decane due to the effect of confinement. However, the suppression is larger when oil (a multicomponent mixture) was used, possibly due to the combined effect of differential uptake of components by kerogen (kerogen prefers and sorbs polars and aromatics more than saturates, leading to splitting of oil into a sorbed and a free phase) and confinement in nano pores. Test-tube swelling, TGA, and DSC experiments were also performed on pyridine(polar-aromatic)-swelled kerogen. The combined and individual contributions from the two effects (the effect of confinement and differential uptake of hydrocarbon components) on properties of liquid in contact with kerogen, are studied in this work. Molecular Dynamics (MD) simulations revealed the variation in the swelling of type II kerogen in the presence of same amount of different liquids (differential swelling of kerogen).
RESUMEN
A parallel implementation of the three-dimensional Shan-and-Chen multicomponent, multiphase lattice Boltzmann method (LBM) was used to simulate the equilibrium distributions of two immiscible fluids in porous media. The simulations were successfully validated against cone-beam x-ray microtomographic data on the distribution of oil (decane), water, and air phases in a 5-mm cube of porous medium composed of packed quartz sand grains. The results confirm that LBM models allow for the straightforward incorporation of complex pore space geometry determined from x-ray microtomography measurements and that simulated wetting and nonwetting phase distributions are consistent with x-ray observations on both macroscopic and microscopic scales.
RESUMEN
Asphaltenic solids formed in the Rangely field in the course of a carbon dioxide flood and heptane insolubles in the oil from the same field were used in this study. Four different solvents were used to dissolve the asphaltenes. Near-infrared (NIR) spectroscopy was used to determine the onset of asphaltene precipitation by heptane titration. When the onset values were plotted versus asphaltene concentrations, distinct break points (called critical aggregation concentrations (CAC) in this paper) were observed. CACs for the field asphaltenes dissolved in toluene, trichloroethylene, tetrahydrofuran, and pyridine occurred at concentrations of 3.0, 3.7, 5.0, and 8.2 g/l, respectively. CACs are observed at similar concentrations as critical micelle concentrations (CMC) for the asphaltenes in the solvents employed and can be interpreted to be the points at which rates of asphaltene aggregations change. CMC values of asphaltenes determined from surface tension measurements (in pyridine and TCE) were slightly higher than the CAC values measured by NIR onset measurements. The CAC for heptane-insoluble asphaltenes in toluene was 3.1 g/l. Thermal gravimetric analysis (TGA) and elemental compositions of the two asphaltenes showed that the H/C ratio of the heptane-insoluble asphaltenes was higher and molecular weight (measured by vapor pressure osmometry) was lower.