分子动力学模拟压水反应堆中联氨对水的影响

本文采用分子动力学方法模拟在常温常压下(1 atm,298 K)和在压水堆环境下(155 atm,626 K),水分子数为256,联氨(N₂H₄)分子数为0,25,50,75等不同数目时,水和联氨粒子系统的动力性质和微观结构.同时探讨了联氨分子的引入对水中溶解氧的影响.从模拟结果可知,在常温常压下,当联氨的分子数为0,25,50,75时,粒子系统的均方位移会随联氨分子数的增加而增加;联氨分子数为0与为25,50,75比较时会少一个数量级;压水堆环境下,联氨分子数 关键词： 分子动力学 压水堆 联氨     

丙烷-氢气-空气预混层流燃烧特性研究

利用高速纹影摄像法和球型发展火焰研究了常温常压下丙烷-氢气-空气预混层流燃烧特性,获得了不同氢气体积分数和当量比下混合气的层流燃烧速率、Markstein数、Zeldovich数和Lewis数.结果表明:随着氢气比例的提高,层流燃烧速率增加,火焰厚度降低;当氢气体积分数小于60%时,随着当量比的增加,Markstein数降低,当氢气体积分数大于60%时,随着当量比的增加,Markstein数增加.当量比小于1.2时,随着氢气比例增加,Markstein数降低.当量比大于1.2时,随着氢气比例增加,Markstein数增加.随着氢气比例的增加,Zeldovich数降低,全局Lewis数降低.     

多壁碳纳米管储氢的物理吸附特性

采用巨正则蒙特卡罗方法 ,模拟常温、1 0MPa下氢在扶手椅型多壁壁碳纳米管中的物理吸附过程 .氢分子之间、氢分子与碳原子之间的相互作用采用Lennard Jones势能模型 .研究了双壁碳纳米管外 (内 )径固定而内 (外 )径改变时的物理吸附储氢情况 ,发现氢分子主要储存在双壁碳纳米管的管壁附近 ,当双壁碳纳米管的内外管壁间距由 0 .34nm增大到 0 .6 1或 0 .88nm时可有效增加物理吸附储氢量 ,并给出了相应的理论解释 .在此基础上 ,计算了管壁间距为 0 .34、0 .6 1和 0 .88nm时的三壁碳纳米管的物理吸附储氢量 ,并与相同条件下单壁和双壁碳纳米管的物理吸附储氢量作了比较 ,发现多壁碳纳米管的物理吸附储氢量随碳管层数的增加而减小 .     

离子化程度对硼取代富勒烯C19Bn+（n=0-3）体系储氢性能的影响

本文采用密度泛函理论方法,研究了氢气分子在硼取代富勒烯C19Bn+（n=0-3）表面的非解离吸附态、化学吸附态、及吸附态之间转换过渡态的结构和能量,得到了吸附态之间转换反应的势能面. 研究表明通过调整体系的离子化电荷可以调整吸附位点的电荷布居数,随着硼取代富勒烯离子化体系C19Bn+的吸氢活性位点电荷布居数增大,氢分子被极化的程度以及氢分子在材料表面的吸附能均逐渐增大,同时改变电荷布居数也对氢气分子在表面吸附态间转变势垒高度产生了量级为几个Kcal/mol的影响. 希望我们的研究能够对富勒烯储氢的进一步发展提供参考.     

交流非平衡等离子体裂解甲烷制C2烃

用介质势垒放电研究了甲烷在无氮条件下的直接转化。天然气（甲烷含量94%）在常温常压下通过介质势垒放电反应器形成非平衡态等离子本。反应器两电极间施圆满的交流电压为15-30kV，甲烷分子在此环境中裂解为CH3、CH2、CH和H，并偶联生成C2烃和氢气。实验表明，甲烷转化率和C2烃收率随放电电压和停留时间增加而增加，虽然C2烃总的选择性变化不大，但乙炔、乙烯各自的选择性随电压升高而增加。     

钴修饰单壁碳纳米管的贮氢性能

本文利用密度泛函理论中的广义梯度近似系统地研究了Co原子修饰的(5,5)单壁碳纳米管(SWCNT)的贮氢性能.结果表明:Co原子以桥位的形式吸附于(5,5)SWCNT表面能量最低,在Co原子周围可以吸附3个完整的氢分子,Co原子的修饰提高了(5,5)SWCNT对氢分子的吸附能力.态密度图显示(5,5)SWCNT具有金属性,磁矩为0,Co-SWCNT和Co-SWCNT·n H_2(n=1-3)的磁矩为1μB,氢分子的吸附未改变体系的磁性.Co-SWCNT·3H_2的氢分子平均吸附能为0.51 e V,适中的吸附能有利于常温常压下实现可逆贮氢.     

锂原子修饰碳原子链团簇的结构和储氢性能研究北大核心CSCD

利用密度泛函理论研究锂原子修饰线型碳原子链团簇Li2Cm（m＝2－8）的结构及其储氢性能．结果表明, Li原子可键合于碳链团簇的两端,Li原子本身不发生团聚,氢在Li2 Cm （ m＝2－8）中能以分子形式吸附,每一个Li原子最多可吸附5个氢分子,氢分子的平均吸附能为0．460～2．276 kcal·mol^－1．其中Li原子修饰C2团簇的质量储氢分数最大,为34．72 wt％,表明了它在常温常压条件下作为储氢材料的可行性．     

平面星形Li_6Si_6团簇的结构及其储氢性能研究

采用密度泛函理论(DFT)方法研究平面星形Li_6Si_6团簇的结构及其储氢性能.结果表明,氢分子能在平面星形Li_6Si_6团簇表面发生吸附,每个Li原子周围均可有效吸附三个氢分子,结构的稳定性及合适的吸氢条件表明平面星形Li_6Si_6团簇在常温常压条件下可以作为储氢媒介.     

H2的振动及转动受激喇曼散射研究

本文研究了不同泵浦光偏振态及氢气气压时,氢的振动及转动受激喇曼散射的竞争效应.在5atm的氢喇曼池中,以波长为532nm的椭偏激光为泵浦源,当椭圆的压缩系数为tg26°～tg32°时,获得了60多条Q(1)和S(1)混合谱线.并在低压氢及线偏或准线偏光泵浦条件下,观察到迄今尚未见报道的受激R(1)喇曼谱线.     

锂原子修饰碳原子链团簇的结构和储氢性能研究

利用密度泛函理论研究锂原子修饰线型碳原子链团簇Li2Cm(m=2—8)的结构及其储氢性能. 结果表明, Li原子可键合于碳链团簇的两端,Li原子本身不发生团聚,氢在Li2Cm (m=2—8)中能以分子形式吸附,每一个Li原子最多可吸附5个氢分子,氢分子的平均吸附能为0.460 ~ 2.276 kcal.mol-1. 其中Li原子修饰C2团簇的质量储氢分数最大,为34.72 wt%,表明了它在常温常压条件下作为储氢材料的可行性.     

离子浓度及表面结构对岩石孔隙内水流动特性的影响

酸性环境引发的岩石孔隙表面溶解增加了孔隙内水溶液的盐离子浓度,破坏了孔隙的表面结构.本文采用分子动力学模拟的方法研究了纳米级岩石孔隙内水溶液的流动特性,分析了盐离子浓度和孔隙表面结构对水流速度分布的影响及原因.研究结果表明:纳米级岩石孔隙内的水溶液流动符合泊肃叶流动特性,流速呈"抛物线"分布;随盐离子浓度增加,水溶液内部氢键网络变得更为致密,水黏度随其呈线性增长;水溶液中离子浓度越大,孔隙表面对水流动的阻力越大,最大流速越小,速度分布的"抛物线"曲率半径越大;岩石孔隙表面结构的破坏改变了流动表面的粗糙程度,增加了孔隙表面对H2O分子的吸引力.随表面结构破坏程度的增大,水溶液在近壁区域的密度增大,流速降低;当表面破坏程度达到50%时,水溶液在近壁区域出现了明显的负边界滑移现象.     

Molecular dynamics simulation studies of cryoprotective agent solutions: the relation between melting temperature and the ratio of hydrogen bonding acceptor to donor number

Molecular dynamics simulations have been carried out for glycerol–water–sodium chloride ternary solution due to its important role in cryopreservation engineering. The radial distribution functions for atom pairs potentially related to C–H ··· O and O–H ··· O hydrogen bonds were calculated. The radial distribution functions for the H (connected to C)–O atom pair do not exhibit peaks between 2 and 3 Å, whereas the radial distribution functions for the C–O atom pair exhibit distinct peaks between 3 and 4 Å. The reason for this is because most C–H ··· O geometries are bent and deviate from linearity. The ratios of acceptor to donor numbers for water and glycerol molecules decrease as the solute concentration increases. A characteristic concentration has been found that divides solutions with different mechanisms. Below the characteristic concentration, the melting temperature is linearly related to the ratio of acceptor to donor number for water molecules, whereas above the characteristic concentration, the melting temperature is linearly related to the ratio of acceptor to donor number for glycerol molecules. Further studies indicate that the relations are independent of hydrogen bonding criteria and temperature.     

A molecular level study of liquid water with a flexible water model

The structural properties of water at different temperatures and pressures have been investigated by using a flexible water model and the inherent structure mechanism. The presence of 60$^\circ$ peak in the O--O--O angle distribution function and the behaviours of the hydrogen bonds in the first shell indicate that some water molecules in the second shell move toward the central molecules through the bending (not breaking) of hydrogen bonds and even become first-shell molecules of the central molecule on the basis of the O--O cutoff distance but not first-shell molecules by means of the hydrogen-bond criterion. The inherent-structure analysis of the O--O radial distribution functions at different pressures shows that the first peak is almost independent of the pressure; the position of second peak moves from 0.45 to 0.32nm as the pressure increases from $1\times 10^5$Pa to $1\times 10^9$Pa. This particularly evident pressure effect, i.e. the constant nearest-neighbours and the transformation of outer-neighbours on the basis of O--O distance, together with the results at different temperatures, gives a positive evidence for the two-state outer-neighbour mixture model: liquid water is a mixture of Ice-Ih-type-bonding and Ice-II-type bonding structures.     

Investigation of the local structure variance of water molecules in laser-induced thermal desorption process

This paper presents the use of molecular dynamics simulation in the study of laser-induced thermal desorption (LITD) of water molecules adjacent to a laser-heated Au substrate. The local structure of the water molecules is investigated by considering the densities of the oxygen and hydrogen atoms, the average number of neighbors, n_NN, and the average number of H-bonds, n_HB. At an equilibrium temperature of 300 K, the simulation results show that three adsorption water layers are formed in the immediate vicinity of the Au surface, and that each four-fold hollow site on the uppermost Au(0 0 1) surface is occupied by a single water molecule. Following laser-induced heating of the Au substrate with a sub-picosecond laser pulse of 350 fs, the substrate temperature increases to 1000 K. This causes a gradual heating of the adjacent water film, which is accompanied by a decrease in the values of n_NN and n_HB. Hence, it can be concluded that an increase in the water film temperature destroys the hydrogen-bonding network throughout the water film. Although the maximum local temperature of the water film occurs in the region immediately adjacent to the Au substrate, it is determined that the attractive energy between the Au atoms and the water molecules in this region causes the water molecules to aggregate together to form three-dimensional water clusters. Furthermore, this energy prevents the hydrogen bonds in this region from breaking apart as violently as those within the phase explosion region. Finally, it is observed that the phase explosion phenomenon occurs in the region of the water film where the values of n_NN and n_HB are at a minimum.     

界面氢键对受限水结构和动态特性的影响

应用反应力场分子动力学方法, 模拟了水限制在全羟基化二氧化硅晶体表面间的弛豫过程, 研究了基底表面与水形成的界面氢键, 及其对受限水结构和动态特性行为的影响. 当基底表面硅醇固定时, 靠近基底表面水分子中的氧原子与基底表面的氢原子形成强氢键, 这使得靠近表面水分子中的氧原子比对应的氢原子更靠近基底表面, 从而水分子的偶极矩远离表面. 当基底表面硅醇可动时, 靠近基底表面水分子与基底表面原子形成两种强氢键, 一种是水分子中的氧原子与表面的氢原子形成的强氢键, 数量较少, 另一种是水分子中的氢原子与表面的氧原子形成的强氢键, 数量较多, 这使得靠近表面水分子中的氢原子比对应的氧原子更靠近表面, 从而水分子的偶极矩指向表面. 在相同几何间距下, 当基底表面硅醇可动时, 表面的活动性使得几何限制作用减弱, 导致了受限水分层现象没有固定表面限制下的明显. 此外, 固定表面比可动表面与水形成的界面氢键作用较强, 数量较多, 导致了可动表面限制下水的运动更为剧烈.     

The effects of water dilution on hydrogen,syngas, and ethylene flames at elevated pressure

This work investigates experimentally and numerically the kinetic effects of water vapor addition on the burning rates of H₂, H₂/CO mixtures, and C₂H₄ from 1 atm to 10 atm at flame temperatures between 1600 K and 1800 K. Burning rates were measured using outwardly propagating spherical flames in a nearly constant pressure chamber. Results show good agreement with newly updated kinetic models for H₂ flames. However, there is considerable disagreement between simulations and measurements for H₂/CO and C₂H₄ flames at high pressure and high water vapor dilution. Both experiments and simulations show that water vapor addition causes a monotonic decrease in mass burning rate and the inhibitory effect increases with pressure. For hydrogen flames, water vapor addition reduces the critical pressure above which a negative pressure dependence of the burning rate is observed. However, for C₂H₄ flames, the burning rate always increases with pressure. The results also show that water vapor addition has the same effect as a pressure increase for H₂ and H₂/CO flames, shifting the reaction zone into a narrower window at higher temperatures. For all fuels, water vapor addition increases OH formation via H₂O + O while reducing the overall active radical pool for hydrogen flames. For C₂H₄, the additional HO₂ production pathway through HCO results in a dramatic difference in pressure dependence of the burning rate from that observed for hydrogen. The present work provides important additions to the experimental database for syngas and C₀–C₂ high pressure kinetic model validations.     

Molecular dynamics study on water self-diffusion in aqueous mixtures of methanol,ethylene glycol and glycerol: investigations from the point of view of hydrogen bonding

Molecular dynamics simulations have been performed to investigate the aqueous binary mixtures of alcohols, including methanol, ethylene glycol (EG) and glycerol of molalities ranging from 1 to 5 m at the temperatures of 273, 288 and 298 K, respectively. The primary purpose of this paper is to investigate the mechanism of water self-diffusion in water-alcohol mixtures from the point of view of hydrogen bonding. The effects of temperature and concentration on water self-diffusion coefficient are evaluated quantitatively in this work. Temperature and concentration to some extent affect the hydrogen bonding statistics and dynamics of the binary mixtures. It is shown that the self-diffusion coefficient of water molecules decreases as the concentration increases or the temperature decreases. Moreover, calculations of mean square displacements of water molecules initially with different number n of H-bonds indicate that the water self-diffusion coefficient decreases as n increases. We also studied the aggregation of alcohol molecules by the hydrophobic alkyl groups. The largest cluster size of the alkyl groups clearly increases as the concentration increases, implying the emergence of a closely connected network of water and alcohols. The clusters of water and alcohol that interacted could block the movement of water molecules in binary mixtures. These findings provide insight into the mechanisms of water self-diffusion in aqueous binary mixtures of methanol, EG and glycerol.     

绝缘油在水和酸共同作用下的分子动力学模拟

为探讨水分和温度对变压器油中酸扩散行为的影响,在分子模拟软件内建立了水分子、矿物油分子、甲酸分子的复合模型,水的体积分数分别为1%,3%,5%,甲酸的质量分数为3%.通过分子动力学方法模拟计算了甲酸分子在不同含水量油模型中的扩散系数、径向分布函数以及复合模型的总动能;此外,选取了水体积分数为1%的复合模型为研究对象,以温度为变量进行动力学模拟.结果表明:随着水含量的增多,甲酸的扩散系数并没有逐渐增大,甲酸分子与水分子之间的氢键数目呈逐渐减少的趋势.而随着温度的升高,油模型具有的动能越大,酸分子的扩散逐渐增强.油模型中,含水量越多,水和酸在热运动下碰撞油分子的能量越大;温度越高,相互之间的碰撞越激烈.模拟结果说明水分增加或温度增高,均对油分子的结构稳定性造成一定程度的破坏.     

关于压水堆产武器级钚的模拟计算

防止核扩散是国际社会共同努力的目标,其中武器级核材料的防扩散是重中之重.钚是反应堆的副产品,如果不计较经济效益,利用铀为核燃料的反应堆都可以生产武器级钚.本文基于日本Takahama-3压水堆建立了五个模型,并进行中子和燃耗计算,得到两种燃料棒产武器级钚的条件、燃料棒轴向的燃耗分布、组件内燃料棒燃耗的变化区间和全堆芯燃料棒径向燃耗分布.基于上述模型和计算数据给出压水堆堆芯内含有武器级钚的准确位置和UO_2燃料棒中武器级钚的产量.这种低燃耗的乏燃料给国际核不扩散带来了巨大风险,国际社会应该加强对此类乏燃料的监管.