Molecular Dynamics Study on the Growth of Structure I Methane Hydrate in Aqueous Solution of Sodium Chloride

Facilitation of Hydrate Dissociation and Structural Evolution by Major Marine Anions under Static Electric Fields

Electric fields have been proven to be capable of significantly affecting the equilibrium state of hydrates. In this study, the thermodynamic properties and structural changes of methane hydrate (MH) in various anion solutions in an electric field at 0.7 V/nm were investigated by molecular dynamics simulations. The presence of anions significantly enhances the instability of methane hydrates under electric fields, leading to a staged dissociation process. First, the anions coexist with MH to form a temporary metastable structure under the action of an electric field. Then, the migration of anions causes the dissociation of nearby hydrates and the formation of flow channels in the hydrate layer, which leads to the complete dissociation of MH after a period. The promotive effects of F-, Br-, I-, and Cl- ions were close, while SO42- was relatively weak. The anions are still in hydration shells in the MH phase, but the structure of the hydration shells differs slightly from that in solution (the coordination numbers of I- and SO42- ions increased). The migration resistances of multiple anions to cross the surface of the hydrate layer are similar. However, inside the hydrate phase, the anions with a larger radius have a higher migration resistance. It is difficult for SO42- ions to migrate inside the hydrate phase, and they tend to form a metastable structure on the hydrate surface. Combining our previous studies, SrCl2 solution has the best hydrate promotion under an electric field environment.

Molecular Dynamics Study of Clathrate-like Ordering of Water in Supersaturated Methane Solution at Low Pressure

Using molecular dynamics, the evolution of a metastable solution for “methane + water” was studied for concentrations of 3.36, 6.5, 9.45, 12.2, and 14.8 mol% methane at 270 K and 1 bar during 100 ns. We have found the intriguing behavior of the system containing over 10,000 water molecules: the formation of hydrate-like structures is observed at 6.5 and 9.45 mol% concentrations throughout the entire solution volume. This formation of “blobs” and the following amorphous hydrate were studied. The creation of a metastable methane solution through supersaturation is the key to triggering the collective process of hydrate formation under low pressure. Even the first stage (0–1 ns), before the first fluctuating cavities appear, is a collective process of H-bond network reorganization. The formation of fluctuation cavities appears before steady hydrate growth begins and is associated with a preceding uniform increase in the water molecule’s tetrahedrality. Later, the constantly presented hydrate cavities become the foundation for a few independent hydrate nucleation centers, this evolution is consistent with the labile cluster and local structure hypotheses. This new mechanism of hydrogen-bond network reorganization depends on the entropy of the cavity arrangement of the guest molecules in the hydrate lattice and leads to hydrate growth.

Promoting Effect of Common Marine Cations on Hydrate Dissociation and Structural Evolution under a Static Electric Field

Natural gas hydrate, a potential energy resource, is attracting worldwide attention. In this study, we propose a new method of hydrate dissociation which uses seawater and electrostatic fields (SE method) cooperatively. The hydrate molecular dissociation mechanism of gas hydrate is a key issue in studying the kinetic properties of gas hydrate using the SE method. Therefore, molecular dynamics simulations were used to investigate the thermodynamic properties and structural changes of methane hydrate (MH) in multiple kinds of salt solutions under an electrostatic field. The results show that the electric field can drive cations into the MH phase to form a series of random semiopen cages, which are essentially temporary and metastable. The variation in free energy indicates that it is more difficult for divalent cations to enter the hydrate phase than monovalent cations, meaning that the hydrate structures formed with divalent cations are more unstable. Then, the ion current occurred in the hydrate phase (called ion migration in this study), which greatly accelerated hydrate dissociation. In contrast, the promotion effect of cations with the same charge on MH dissociation is as follows: Sr²⁺ > K⁺ ≈ Na⁺ > Ca²⁺ ≈ Mg²⁺. In general, the presence of common marine cations enhanced the promotion effect of the electric field on gas hydrate dissociation.

Molecular Simulation of the Effects of Cyclic Organic Compounds on the Stability of Lccbm Hydrates

CH₄ can be separated from low-concentration coal bed methane (LCCBM) by using the hydrate-based gas separation (HBGS) method. To study the contribution of different cyclic organic compounds to the separation of CH₄ in LCCBM, an LCCBM hydrate model was constructed. Based on the Monte Carlo and molecular dynamics theory, we simulated the effect of three cyclic organic compounds—cyclopentane (CP), cyclopentanone (CP-one), and cyclopentanol (CP-ol)—on the stability of the LCCBM hydrate at P = 2 MPa, various temperatures, and discussed the structural stability of the hydrate in depth in terms of final snapshots, radial distribution function, mean square displacement, diffusion coefficient, and potential energy change. The results showed that for the CH₄-N₂ LCCMM gas mixture, CP showed the best facilitation effect compared to the other two cyclic compounds by maintaining the stability of the LCCBM hydrate well at T = 293 K. The promotion effect of CP-one is between CP and CP-ol, and when the temperature increases to T = 293 K, the oxygen atoms in the water molecule can maintain the essential stability of the hydrate structure, although the orderliness decreases significantly. Moreover, the structure of the hydrate model containing CP-ol is destroyed at T = 293 K, and the eventual escape of CH₄ and N₂ molecules in solution occurs as bubbles. The research results are important for further exploration of the mechanism of action of cyclic promoter molecules with LCCBM hydrate molecules and promoter preferences.

Growth Rate and Morphology Study of Tetrahydrofuran Hydrate Single Crystals and the Effect of Salt

For an advanced gas hydrate management strategy in flow assurance as well as the application of a gas hydrate technology (e.g., desalination), it is important to understand the kinetics and...

Microscopic Molecular Insights into Hydrate Formation and Growth in Pure and Saline Water Environments

The growth dynamics of natural gas hydrates in saline water has been studied using copious experiments and spectroscopic observations; however, the microscopic evidences to the structural and molecular transformations that they have provided are poorly understood. In this view, we perform extensive molecular dynamics simulations to gain physical insights into the formation and growth mechanism of naturally occurring gas hydrates with a wide variation in the amount of methane (1:5 to 1:18 methane/water ratio) in pure and salt (0-5 wt %) water environments at 50 MPa and 260 K. A couple of new findings analyzed from the number of cages and F_4φ order parameter are as follows: (a) 1:6 (methane/water ratio) is an optimum ratio for the rapid growth of a properly ordered hydrate in pure water at which the hydrate growth retards with increasing salt concentration, (b) there is an inconsequential difference between methane hydrate dynamics in pure water and 0.8 and 1.5 wt % salt water at a ratio of 1:12 (methane/water), and (c) lower methane (1:18) and salt (0.8 wt %) concentrations promote hydrate growth. Besides, this study observes the structural coexistence of S-I and S-II methane hydrates as the large 5¹²6⁴ cages appear along with the small 5¹² and large 5¹²6² cages, in which the low methane concentration favors the S-II structure.

Molecular Dynamics Simulation Study on the Growth of Structure II Nitrogen Hydrate

Crystal growth of N₂ hydrate in a three-phase system consisting of N₂ hydrate, liquid water, and gaseous N₂ was performed by molecular dynamics simulation at 260 K. Pressure influence on hydrate growth was evaluated. The kinetic properties including the growth rates and cage occupancies of the newly formed hydrate and the diffusion coefficient and concentration of N₂ molecules in liquid phase were measured. The results showed that the growth of N₂ hydrate could be divided into two stages where N₂ molecules in gas phase had to dissolve in liquid phase and then form hydrate cages at the liquid-hydrate interface. The diffusion coefficient and concentration of N₂ in liquid phase increased linearly with increasing pressure. As the pressure rose from 50 to 100 MPa, the hydrate growth rate kept increasing from 0.11 to 0.62 cages·ns^-1·Å^-2 and then dropped down to around 0.40 cages·ns^-1·Å^-2 once the pressure surpassed 100 MPa. During the hydrate formation, the initial sII N₂ hydrate phase set in the system served as a template for the subsequent growth of N₂ hydrate so that no new crystal structure was found. Analysis on the cage occupancies revealed that the amount of cages occupied by two N₂ molecules increased evidently when the pressure was above 100 MPa, which slowed down the growth rate of hydrate cages. Additionally, a small fraction of defective cages including two N₂ molecules trapped in 5¹²6⁵ cages and three N₂ molecules trapped 5¹²6⁸ cages was observed during the hydrate growth.

Hydrophobic Hydration and the Effect of NaCl Salt in the Adsorption of Hydrocarbons and Surfactants on Clathrate Hydrates

Adsorption of functional molecules on the surface of hydrates is key in the understanding of hydrate inhibitors. We investigate the adsorption of a hydrocarbon chain, nonionic and ionic surfactants, and ions at the hydrate-aqueous interface. Our results suggest a strong connection between the water ordering around solutes in bulk and the affinity for the hydrates surface. We distinguish two types of water ordering around solutes: (i) hydrophobic hydration where water molecules form a hydrogen bond network similar to clathrate hydrates, and (ii) ionic hydration where water molecules align according to the polarity of an ionic group. The nonionic surfactant and the hydrocarbon chain induce hydrophobic hydration and are favorably adsorbed on the hydrate surface. Adsorption of ions and the ionic headgroups on the hydrate surface is not favorable because ionic hydration and the hydrogen bond structure of hydrates are incompatible. The nonionic surfactant is adsorbed by the headgroup and tail while adsorption of the ionic surfactants is not favorable through the head. Water ordering is analyzed using the hydrogen bond and tetrahedral density profiles as a function of the distance to the chemical groups. The adsorption of solutes is studied through the free energy profiles as a function of the distance to the hydrate surface. Salt lowers the melting temperature of hydrates, disrupts hydrophobic hydration, reduces the solubility of solutes in the aqueous solution, and increases the propensity of solutes to be adsorbed on hydrate surfaces. Our studies are performed by the unbiased and steered molecular dynamics simulations. The results are in line with experiments on the effect of salt and alkanes in hydrate antiagglomeration.

Adsorption of Kinetic Hydrate Inhibitors on Growing Surfaces: A Molecular Dynamics Study

We investigate the mechanism of a typical kinetic hydrate inhibitor (KHI), polyvinylcaprolactam (PVCap), which has been applied to prevent hydrate plugs from forming in gas pipe lines, using molecular dynamics simulations of crystal growth of ethylene oxide hydrate. Water-soluble ethylene oxide is chosen as a guest species to avoid problems associated with the presence of the gas phase in the simulation cell such as slow crystal growth. A PVCap dodecamer adsorbs irreversibly on the hydrate surface which grows at supercooling of 3 K when the hydrophobic part of two pendent groups are trapped in open cages at the surface. The amide hydrogen bonds make no contribution to the adsorption. PVCap can adsorb on various crystallographic planes of sI hydrate. This is in contrast to antifreeze proteins, each of which prefers a specific plane of ice. The trapped PVCap gives rise to necessarily the concave surface of the hydrate. The crystal growth rate decreases with increasing surface curvature, indicating that the inhibition by PVCap is explained by the Gibbs-Thomson effect.

The molecular mechanism of the inhibition effects of PVCaps on the growth of sI hydrate: an unstable adsorption mechanism

Adsorption and non-binding-hydrate sites produce synergistic effects that lead to unstable adsorption and inhibition effects (recurring destruction of hydrate).

A molecular dynamic study on the dissociation mechanism of SI methane hydrate in inorganic salt aqueous solutions

Gas hydrate is not only a potential energy resource, but also almost the biggest challenge in oil/gas flow assurance. Inorganic salts such as NaCl, KCl and CaCl₂ are widely used as the thermodynamic inhibitor to reduce the risk caused by hydrate formation. However, the inhibition mechanism is still unclear. Therefore, molecular dynamic (MD) simulation was performed to study the dissociation of structure I (SI) methane hydrate in existence of inorganic salt aqueous solution on a micro-scale. The simulation results showed that, the dissociation became stagnant due to the presence of liquid film formed by the decomposed water molecules, and more inorganic ions could shorten the stagnation-time. The diffusion coefficients of ions and water molecules were the largest in KCl system. The structures of ion/H₂O and H₂O/H₂O were the most compact in hydrate/NaCl system. The ionic ability to decompose hydrate cells followed the sequence of: Ca²⁺>2K⁺>2Cl^->2Na⁺.

Effect of Greenhouse Gases Dissolved in Seawater

A molecular dynamics simulation has been performed on the greenhouse gases carbon dioxide and methane dissolved in a sodium chloride aqueous solution, as a simple model of seawater. A carbon dioxide molecule is also treated as a hydrogen carbonate ion. The structure, coordination number, diffusion coefficient, shear viscosity, specific heat, and thermal conductivity of the solutions have been discussed. The anomalous behaviors of these properties, especially the negative pressure dependence of thermal conductivity, have been observed in the higher-pressure region.

Effects of thermodynamic inhibitors on the dissociation of methane hydrate: a molecular dynamics study

We investigate the effects of methanol and NaCl, which are known as thermodynamic hydrate inhibitors, on the dissociation kinetics of methane hydrate in aqueous solutions by using molecular dynamics simulations. It is shown that the dissociation rate is not constant but changes with time. The dissociation rate in the initial stage is increased by methanol whereas it is decreased by NaCl. This difference arises from the opposite effects of the two thermodynamic inhibitors on the hydration free energy of methane. The dissociation rate of methane hydrate is increased by the formation of methane bubbles in the aqueous phase because the bubbles absorb surrounding methane molecules. It is found that both methanol and NaCl facilitate the bubble formation. However, their mechanisms are completely different from each other. The presence of ions enhances the hydrophobic interactions between methane molecules. In addition, the ions in the solution cause a highly non-uniform distribution of dissolved methane molecules. These two effects result in the easy formation of bubbles in the NaCl solution. In contrast, methanol assists the bubble formation because of its amphiphilic character.