Thermal conductivity of methane hydrate from experiment and molecular simulation

A single-sided transient plane source technique has been used to determine the thermal conductivity and thermal diffusivity of a compacted methane hydrate sample over the temperature range of 261.5-277.4 K and at gas-phase pressures ranging from 3.8 to 14.2 MPa. The average thermal conductivity, 0.68 +/- 0.01 W/(m K), and thermal diffusivity, 2.04 x 10(-7) +/- 0.04 x 10(-7) m2/s, values are, respectively, higher and lower than previously reported values. Equilibrium molecular dynamics (MD) simulations of methane hydrate have also been performed in the NPT ensemble to estimate the thermal conductivity for methane compositions ranging from 80 to 100% of the maximum theoretical occupation, at 276 K and at pressures ranging from 0.1 to 100 MPa. Calculations were performed with three rigid potential models for water, namely, SPC/E, TIP4P-Ew, and TIP4P-FQ, the last of which includes the effects of polarizability. The thermal conductivities predicted from MD simulations were in reasonable agreement with experimental results, ranging from about 0.52 to 0.77 W/(m K) for the different potential models with the polarizable water model giving the best agreement with experiments. The MD simulation method was validated by comparing calculated and experimental thermal conductivity values for ice and liquid water. The simulations were in reasonable agreement with experimental data. The simulations predict a slight increase in the thermal conductivity with decreasing methane occupation of the hydrate cages. The thermal conductivity was found to be essentially independent of pressure in both simulations and experiments. Our experimental and simulation thermal conductivity results provide data to help predict gas hydrate stability in sediments for the purposes of production or estimating methane release into the environment due to gradual warming.     

Measurement and modeling of the effective thermal conductivity for porous methane hydrate samples

The effective thermal conductivities of gas-saturated porous methane hydrates were measured by a single-sided transient plane source (TPS) technique and simulated by a generalized fractal model of porous media that based on self-similarity.The density of porous hydrate,measured by the volume of the sample in the experimental system,was used to evaluate the porosity of methane hydrate samples.The fractal model was based on Sierpinski carpet,a thermal-electrical analogy technique and one-dimensional heat flow assumption.Both the experimental and computational results show the effective thermal conductivity of methane hydrate decreases with the porosity increase.The porosity of 0.3 can reduce the thermal conductivity of the methane hydrate by 25%.By analysis of the experimental data and the simulative result,the optimized thermal conductivity of the zero-porosity methane hydrate is about 0.7 W m-1K-1.     

Molecular dynamics simulations of the mechanisms of thermal conduction in methane hydrates

The thermal conductivity of methane hydrate is an important physical parameter affecting the processes of methane hydrate exploration,mining,gas hydrate storage and transportation as well as other applications.Equilibrium molecular dynamics simulations and the Green-Kubo method have been employed for systems from fully occupied to vacant occupied sI methane hydrate in order to estimate their thermal conductivity.The estimations were carried out at temperatures from 203.15 to 263.15 K and at pressures from 3 to 100 MPa.Potential models selected for water were TIP4P,TIP4P-Ew,TIP4P/2005,TIP4P-FQ and TIP4P/Ice.The effects of varying the ratio of the host and guest molecules and the external thermobaric conditions on the thermal conductivity of methane hydrate were studied.The results indicated that the thermal conductivity of methane hydrate is essentially determined by the cage framework which constitutes the hydrate lattice and the cage framework has only slightly higher thermal conductivity in the presence of the guest molecules.Inclusion of more guest molecules in the cage improves the thermal conductivity of methane hydrate.It is also revealed that the thermal conductivity of the sI hydrate shows a similar variation with temperature.Pressure also has an effect on the thermal conductivity,particularly at higher pressures.As the pressure increases,slightly higher thermal conductivities result.Changes in density have little impact on the thermal conductivity of methane hydrate.     

Estimation of ultra-stability of methane hydrate at 1 atm by thermal conductivity measurement

Thermal conductivity of methane hydrate was measured in hydrate dissociation self-preservation zone by means of the transient plane source (TPS) technique developed by Gustafsson. The sample was formed from 99.9% (volume ratio) methane gas with 280 ppm sodium dodecyl sulfate (SDS) solution under 6.6 MPa and 273.15 K. The methane hydrate sample was taken out of the cell and moved into a low temperature chamber when the conversion ratio of water was more than 90%. In order to measure the thermal conductivity, the sample was compacted into two columnar parts by compact tool at 268.15 K. The measurements are carried out in the temperature ranging from 263.15 K to 271.15 K at atmospheric pressure. Additionally, the relationship between thermal conductivity and time is also investigated at 263.15 K and 268.15 K, respectively. In 24 h, thermal conductivity increases only 5.45% at 268.15 K, but thermal conductivity increases 196.29% at 263.15 K. Methane hydrates exhibit only minimal decomposition at 1 atm and the temperature ranging from 263.15 K to 271.15 K. At 1 atm and 268.15 K, the total gas that evolved after 24 h was amounted to less than 0.71% of the originally stored gas, and this ultra-stability was maintained if the test was lasted for more than two hundreds hours before terminating.     

甲烷水合物导热机理的分子动力学模拟

甲烷水合物导热系数是甲烷水合物勘探、开采、储运以及其他应用过程中一个十分重要的物理参数.我们采用平衡分子动力学（EMD）方法Green-Kubo理论计算温度203.15～263.15K、压力范围3～100MPa、晶穴占有率为0～1的sI甲烷水合物的导热系数,采用的水分子模型包括TIP4P、TIP4P-Ew、TIP4P-FQ、TIP4P/2005、TIP4P/Ice.研究了主客体分子、外界温压条件等对甲烷水合物导热性能的影响.研究结果显示甲烷水合物的低导热性能由主体分子构建的sI笼型结构决定,而客体分子进入笼型结构后,使得笼型结构导热性能增强,同时进入笼型结构的客体分子越多,甲烷水合物导热性能越强.研究结果还显示在高温区域（T〉TDebye/3）内不同温度作用下,所有sI水合物具有相似的导热规律.压力对导热系数有一定影响,尤其是在较高压力条件下,压力越高,导热系数越大.而在不同温度和不同压力作用过程中,密度的改变对导热系数的增大或减小几乎没有影响.     

Molecular dynamics simulations of the thermal conductivity of methane hydrate

Nonequilibrium molecular dynamics simulations with the nonpolarizable SPC/E (Berendsen et al., J. Phys. Chem. 1987, 91, 6269) and the polarizable COS/G2 (Yu and van Gunsteren, J. Chem. Phys. 2004, 121, 9549) force fields have been employed to calculate the thermal conductivity and other associated properties of methane hydrate over a temperature range from 30 to 260 K. The calculated results are compared to experimental data over this same range. The values of the thermal conductivity calculated with the COS/G2 model are closer to the experimental values than are those calculated with the nonpolarizable SPC/E model. The calculations match the temperature trend in the experimental data at temperatures below 50 K; however, they exhibit a slight decrease in thermal conductivity at higher temperatures in comparison to an opposite trend in the experimental data. The calculated thermal conductivity values are found to be relatively insensitive to the occupancy of the cages except at low (T相似文献   

Experimental and modeling study on decomposition kinetics of methane hydrates in different media

The decomposition kinetic behaviors of methane hydrates formed in 5 cm3 porous wet activated carbon were studied experimentally in a closed system in the temperature range of 275.8-264.4 K. The decomposition rates of methane hydrates formed from 5 cm3 of pure free water and an aqueous solution of 650 g x m(-3) sodium dodecyl sulfate (SDS) were also measured for comparison. The decomposition rates of methane hydrates in seven different cases were compared. The results showed that the methane hydrates dissociate more rapidly in porous activated carbon than in free systems. A mathematical model was developed for describing the decomposition kinetic behavior of methane hydrates below ice point based on an ice-shielding mechanism in which a porous ice layer was assumed to be formed during the decomposition of hydrate, and the diffusion of methane molecules through it was assumed to be one of the control steps. The parameters of the model were determined by correlating the decomposition rate data, and the activation energies were further determined with respect to three different media. The model was found to well describe the decomposition kinetic behavior of methane hydrate in different media.     

Comparison of the water change characteristics between the formation and dissociation of methane hydrate and the freezing and thawing of ice in sand

Hydrate formation and dissociation processes are always accompanied by water migration in porous media, which is similar to the ice. In our study, a novel pF-meter sensor which could detect the changes of water content inside sand was first applied to hydrate formation and dissociation processes. It also can study the water change characteristics in the core scale of a partially saturated silica sand sample and compare the differences of water changes between the processes of formation and dissociation of methane hydrate and freezing and thawing of ice. The experimental results showed that the water changes in the processes of formation and dissociation of methane hydrate were basically similar to that of the freezing and thawing of ice in sand. When methane hydrate or ice was formed, water changes showed the decrease in water content on the whole and the pF values rose following the formation processes. However, there were very obvious differences between the ice thawing and hydrate dissociation.     

活性炭中甲烷水合物的分解动力学

在封闭体系内,在初始分解压力0.1 MPa,温度范围276～265 K之间,测定了 五组甲烷水合物在活性炭中的解动力学数据。分析了甲烷水合物在活性炭中分解的 物理过程,提出了以微分方程表达的宏观分解动力学模型。使用单步积分的吉尔（ Gear)方法解得微分方程的数值解,结合单纯形最优化方法拟合模型参数,模型计 算值与实验值符合良好。     

Texture change of ice on anomalously preserved methane clathrate hydrate

We used a confocal scanning microscope to observe growth and texture change of ice due to the dissociation of methane gas clathrate hydrate (CH(4) hydrate). The experiments were done under CH(4) gas atmospheric pressure and isothermal conditions between 170 and 268 K. Above 193 K, the dissociation of CH(4) hydrate resulted in many small ice particles that covered the hydrate surface. These ice particles had roughly the same shape and density between 193 and 210 K. In contrast, above 230 K the ice particles developed into a sheet of ice that covered the hydrate surface. Moreover, the measured release of CH(4) gas decreased when the sheet of ice formed at the surface of the hydrate. These findings can explain the anomalous preservation behavior of CH(4) hydrate; that is, the known increase of storage stability of CH(4) hydrate above 240 K is likely related to the formation of the ice that we observed in the experiments.     

三元脂肪酸/膨胀石墨复合相变材料的制备、包覆定形及热性能

采用癸酸、 月桂酸和棕榈酸的三元共晶混合物作为相变材料, 以膨胀石墨为基体, 通过膨胀石墨多孔结构的毛细吸附和复合涂饰剂的包覆定形, 将多元相变材料固定在膨胀石墨的孔道结构中, 制备出结构稳定、 密封性能优异、 热稳定性好和高导热的新型三元脂肪酸/膨胀石墨复合定形相变材料. 膨胀石墨具有膨胀疏松的多孔结构和良好的吸附性能; 其熔融潜热为95.6 J/g, 结晶焓为82.8 J/g, 说明其具有很好的相变蓄热特性和热循环稳定性; 材料的导热性能可增加至0.738 W/(m·K), 与脂肪酸相比得到大幅度提高.     

Activation energy of methyl radical decay in methane hydrate

The thermal stability of gamma-ray-induced methyl radicals in methane hydrate was studied using the ESR method at atmospheric pressure and 210-260 K. The methyl radical decay proceeded with the second-order reaction, and ethane molecules were generated from the dimerization process. The methyl radical decay proceeds by two different temperature-dependent processes, that is, the respective activation energies of these processes are 20.0 +/- 1.6 kJ/mol for the lower temperature region of 210-230 K and 54.8 +/- 5.7 kJ/mol for the higher temperature region of 235-260 K. The former agrees well with the enthalpy change of methane hydrate dissociation into ice and gaseous methane, while the latter agrees well with the enthalpy change into liquid water and gaseous methane. The present findings reveal that methane hydrates dissociate into liquid (supercooled) water and gaseous methane in the temperature range of 235-260 K.     

Hydrate kinetics study in the presence of nonaqueous liquid by nuclear magnetic resonance spectroscopy and imaging

The dynamics of methane hydrate growth and decomposition were studied by nuclear magnetic resonance (NMR) spectroscopy and imaging (MRI). Three well-known large molecule guest substances (LMGS) were used as structure H hydrate formers: 2,2-dimethylbutane (NH), methylcyclohexane (MCH), tert-butyl methyl ether (TBME). In addition, the impact of a non-hydrate former (n-heptane/nC7) was studied. The methane diffusion and hydrate growth were monitored by recording the 2H NMR spectra at 253 K and approximately 4.5 MPa for 20 h. The results revealed that methane diffuses faster in TBME and NH, slower in nC7, and slowest in MCH. The TBME system gives the fastest hydrate formation kinetics followed by NH, MCH, and nC7. The conversion of water into hydrate was also observed. The imaging study showed that TBME has a strong affinity toward ice, which is not the case for the NH and MCH systems. The degree of ice packing was also found to affect the LMGS distribution between ice particles. Highly packed ice increases the mass transfer resistance and hence limits the contact between LMGS and ice. It was also found that "temperature ramping" above the ice point improves the conversion significantly. Finally, hydrates were found to dissociate quickly within the first hour at atmospheric pressure and subsequently at a much slower rate. Methane dissolved in LMGS was also seen. The residual methane in hydrate phase and dissolved in LMGS phase explain the faster kinetics during hydrate re-formation.     

Gas hydrates of argon and methane synthesized at high pressures: composition, thermal expansion, and self-preservation

For the first time, the compositions of argon and methane high-pressure gas hydrates have been directly determined. The studied samples of the gas hydrates were prepared under high-pressure conditions and quenched at 77 K. The composition of the argon hydrate (structure H, stable at 460-770 MPa) was found to be Ar.(3.27 +/- 0.17)H(2)O. This result shows a good agreement with the refinement of the argon hydrate structure using neutron powder diffraction data and helps to rationalize the evolution of hydrate structures in the Ar-H(2)O system at high pressures. The quenched argon hydrate was found to dissociate in two steps. The first step (170-190 K) corresponds to a partial dissociation of the hydrate and the self-preservation of a residual part of the hydrate with an ice cover. Presumably, significant amounts of ice Ic form at this stage. The second step (210-230 K) corresponds to the dissociation of the residual part of the hydrate. The composition of the methane hydrate (cubic structure I, stable up to 620 MPa) was found to be CH(4).5.76H(2)O. Temperature dependence of the unit cell parameters for both hydrates has been also studied. Calculated from these results, the thermal expansivities for the structure H argon hydrate are alpha(a) = 76.6 K(-1) and alpha(c) = 77.4 K(-1) (in the 100-250 K temperature range) and for the cubic structure I methane hydrate are alpha(a) = 32.2 K(-1), alpha(a) = 53.0 K(-1), and alpha(a) = 73.5 K(-1) at 100, 150, and 200 K, respectively.     

石蜡/膨胀石墨/石墨烯复合相变储热材料的制备及性能

以石蜡(PA)作为相变储热材料、 膨胀石墨(EG)作为主导热材料和支撑材料, 石墨烯气凝胶(GA)作为导热增强材料和辅支撑材料制备了PA/EG/GA复合相变材料, 研究了GA添加量对复合相变材料相变温度、 相变潜热、 导热性能以及循环稳定性的影响. 结果表明, 所制备的80%PA-17%EG-3%GA复合相变材料导热性能良好, 循环稳定性出色. 与80%PA-20%EG复合材料相比, 该材料的相变温度、 相变潜热以及循环稳定性无明显变化, 但导热系数由4.089 W/(m·K)提升到了5.336 W/(m·K), 显示出良好的应用前景.     

Preparation of porous polyimide/in‐situ reduced graphene oxide composite films for electromagnetic interference shielding

Herein we report an easy and efficient approach to prepare lightweight porous polyimide (PI)/reduced graphene oxide (RGO) composite films. First, porous poly (amic acid) (PAA)/graphene oxide (GO) composite films were prepared via non‐solvent induced phase separation (NIPS) process. Afterwards PAA was converted into PI through thermal imidization and simultaneously GO dispersed in PAA matrix was in situ thermally reduced to RGO. The GO undergoing the same thermal treatment process as thermal imidization was characterized with thermogravimetric analysis, Raman spectra, X‐ray photoelectron spectroscopy and X‐ray diffraction to demonstrate that GO was in situ reduced during thermal imidization process. The resultant porous PI/RGO composite film (500‐µm thickness), which was prepared from pristine PAA/GO composite with 8 wt% GO, exhibited effective electrical conductivity of 0.015 S m^?1 and excellent specific shielding efficiency value of 693 dB cm² g^?1. In addition, the thermal stability of the porous PI/RGO composite films was also dramatically enhanced. Compared with that of porous PI film, the 5% weight loss temperature of the composite film mentioned above was improved from 525°C to 538°C. Moreover, tensile test showed that the composite film mentioned above possessed a tensile strength of 6.97 MPa and Young's modulus of 545 MPa, respectively. Copyright © 2016 John Wiley & Sons, Ltd.     

Thermodynamic properties of methane hydrate in quartz powder

Using the experimental method of precision adiabatic calorimetry, the thermodynamic (equilibrium) properties of methane hydrate in quartz sand with a grain size of 90-100 microm have been studied in the temperature range of 260-290 K and at pressures up to 10 MPa. The equilibrium curves for the water-methane hydrate-gas and ice-methane hydrate-gas transitions, hydration number, latent heat of hydrate decomposition along the equilibrium three-phase curves, and the specific heat capacity of the hydrate have been obtained. It has been experimentally shown that the equilibrium three-phase curves of the methane hydrate in porous media are shifted to the lower temperature and high pressure with respect to the equilibrium curves of the bulk hydrate. In these experiments, we have found that the specific heat capacity of the hydrate, within the accuracy of our measurements, coincides with the heat capacity of ice. The latent heat of the hydrate dissociation for the ice-hydrate-gas transition is equal to 143 +/- 10 J/g, whereas, for the transition from hydrate to water and gas, the latent heat is 415 +/- 15 J/g. The hydration number has been evaluated in the different hydrate conditions and has been found to be equal to n = 6.16 +/- 0.06. In addition, the influence of the water saturation of the porous media and its distribution over the porous space on the measured parameters has been experimentally studied.     

Polymorphic Mn/W/Na(K,Rb,Cs)/SiO<Subscript>2</Subscript> catalysts for oxidative coupling of methane

The phase composition of Mn/W/Na(K,Rb,Cs)/SiO₂ composites-catalysts was determined by X-ray diffraction. Polymorphic transformations involving tridymite (a SiO₂-based nonstoichiometric crystalline phase) were shown to be a reason for the high catalytic activity of Mn/W/Na(K,Rb,Cs)/SiO₂ in heterogeneous oxidative coupling of methane.     

Formation of methane hydrate from polydisperse ice powders

Neutron diffraction runs and gas-consumption experiments based on pressure-volume-temperature measurements are conducted to study the kinetics of methane hydrate formation from hydrogenated and deuterated ice powder samples in the temperature range of 245-270 K up to high degrees of transformation. An improved theory of the hydrate growth in a polydisperse ensemble of randomly packed ice spheres is developed to provide a quantitative interpretation of the data in terms of kinetic model parameters. This paper continues the research line of our earlier study which was limited to the monodisperse case and shorter reaction times (Staykova et al., 2003). As before, we distinguish the process of initial hydrate film spreading over the ice particle surface (stage I) and the subsequent hydrate shell growth (stage II) which includes two steps, i.e., an interfacial clathration reaction and the gas and water transport (diffusion) through the hydrate layer surrounding the shrinking ice cores. Although kinetics of hydrate formation at stage II is clearly dominated by the diffusion mechanism which becomes the limiting step at temperatures above 263 K, both steps are shown to be essential at lower temperatures. The permeation coefficient D is estimated as (1.46 +/- 0.44) x 10(-12) m2/h at 263 K with an activation energy Q(D) approximately 52.1 kJ/mol. This value is close to the energy of breaking hydrogen bonds in ice Ih and suggests that this process is the rate-limiting step in hydrate formation from ice in the slower diffusion-controlled part of the reaction.     

Spherical filler-promoting thermally conductive pathway in graphite-containing polymer composites for high heat radiation

Graphite is an efficient and affordable filler for polymer composites, allowing the control of thermal conductivity. In comparison to other thermally conductive fillers, graphite is lightweight and flexible but affords anisotropic thermal conductivity. Herein, the control of thermal conductivity of graphite-containing polymer composite sheet using spherical polymer particles as additional fillers is described. The thermal conductivity in the through-plane direction (λ_t) of the composite sheet is enhanced by varying the composition ratio of the two fillers (flaky graphite and spherical particles), and optimizing the forming temperature and pressure. Graphite-containing (25 wt%) polymer composite sheet formed by compression at 150 °C and 10 MPa exhibits λ _t value of 0.66 W/m K. Upon mixing of polystyrene microspheres, λ _t is successfully increased. The maximum value of thermal conductivity for a composite sheet with 35 wt% of graphite and 50 wt% of spherical particles is 7.51 W/m K, at 180 °C and 10 MPa. The graphite-containing polymer matrix forms a sequentially connected network-like structure in the composite sheet. Excess polymer microspheres lead to the formation of void structures inside the composite sheet, reducing the thermal conductivity. Thermo-camera observations proved that the composite sheets with higher λ _t value showed comparably high heat radiations. © 2020 Wiley Periodicals, Inc. J. Polym. Sci. 2020 , 58, 607–615