《河北科技大学学报》模版

第36卷第期 河北科技大学学报 Vol.36,No.2015年月 JournalofHebeiUniversityofScienceandTechnology 2015文章编号:1008-1542(2015)00-0001-00 doi:10.7535/hbkd.2015yx铁基换热面CaCO3析晶污垢结合力

分子动力学模拟陈帅帅L程延海L徐金龙2，彭玉兴L韩东太1(1.中国矿业大学机电工程学院，江苏徐州 221116;2.湛江出入境检验检疫局技术中心，广东湛江524022)摘要：为了研究腐蚀产物Fe2O3在污垢与换热表面相互作用时的作用机理，采用分子动力学(MD)方法模拟了Fe和Fe2O3吸附方解石(1-10)晶面的体系，并且对2个体系进行吸附过程的平衡处理，通过计算2个体系吸附后的各种能量推导出吸附结构中结合能的大小，将Fe和Fe2O3吸附方解石(1-10)的2个体系的结合能进行了比较。计算结果表明：Fe2O3吸附方解石(1-10)晶面的结合能要比Fe吸附方解石(1-10)晶面的结合能高，同时，径向分布函数的计算也进一步证明了该结果。该理论计算得出，Fe的腐蚀产物Fe2O3吸附的污垢比Fe吸附的污垢更难以清理。关键词：方解石；腐蚀产物；Fe2O3；分子动力学模拟；吸附；结合能中图分类号：TK172 文献标志码：AMoleculardynamicssimulationofCaCO3crystallizationfoulingadheredoniron-basedmaterialsurfaceCHENShuaishuai1,CHENGYanhai1,XUJinlong2,PENGMixing1,HANDongtai1(1.SchoolofMechatronicEngineering,ChinaUniversityofMiningandTechnology,Xuzhou,Jiangsu221116,China;2.ScienceandTechnologyCenter,ZhanjiangEntry-ExitInspectionandQuarantineBureau,Zhanjiang,Guangdong524022,China)Abstract： InordertostudytheformationmechanismoffoulingwiththeactionofFe2O3corrosionproduct,theabsorptionsystemofFeandFe2O3oncalcite(1-10)crystalfacewassimulatedandtreatedtobebalancedintermsofmoleculardynamicsmethod.Theadhesionenergyoftwoabsorptionsystemswasdeducedandcomparedbycalculatingallenergyoftwoabsorptionsystems.TheresultsshowthattheadhesionenergyofFeoncalcite(1-10)crystalfaceishigherthanthatofFe2O3oncalcite(1-10)crystalface.Meanwhile,theresultsarefurtherprovedbythecalculationofradialdistributionfunction.ThetheoreticalcalculationindicatesitishardtocleanupthecrystallizationfoulingadheredonFesurfacethanthatonFe2O3surface.KeyW0rds:calcite;corrosionproduct;Fe2O3;moleculardynamicssimulation;adsorption;adhesionenergy收稿日期：2015-00-00；修回日期：2015-00-00；责任编辑：基金项目：国家自然科学基金（51006117）；中国博士后特别资助项目（2012T50497）作者简介：陈帅帅（1988一），男，山东泰安人，硕士研究生，主要从事换热面污垢形成机理及防垢方面的研究。通讯作者：程延海副教授。E-mail：chyh1007@陈帅帅，程延海，徐金龙，等.铁基换热面CaCO3析晶污垢结合力分子动力学模拟[J].河北科技大学学报，2015,36(0)：000-000.CHENShuaishuai，CHENGYanhai,XUJinlong，etal.MoleculardynamicssimulationofCaCC^crystallizationfoulingadheredoniron-basedmaterialsurface[J].JournalofHebeiUniversityofScienceandTechnology2014,36(0)：000-000.污垢广泛存在于工业生产的许多领域中，据统计，90%以上的换热设备都存在着不同程度的污垢问题［1］,污垢的存在降低了传热效率，同时，增大了流动阻力，造成了巨大的经济损央1-5］。同时，在安全方面，换热面上污垢的积聚也经常引起局部过热或超温而导致机械性能下降，甚至引发事故。污垢的积聚还常引起换热面的局部腐蚀乃至穿孔，严重威胁了换热设备的安全运行。在热交换器污垢的对策方面，污垢与基体间的结合力一直是一个难以预计的问题。工业生产中的一些工况下污垢难以清理，而一些工况下也出现污垢自然脱落的现象。目前研究结合力常用的方法有电极法和划痕试验法［6-8］，电极法常用在计算气体间的结合力，划痕试验法常用在计算镀层与基体间的结合力，但这2种实验方法在计算基体与污垢的结合力方面应用效果很差，不能得到精确的结合力数据中任］。而在理论方法方面，计算结合力的有量子力学方法和分子动力学方法，其中分子动力学理论在计算基体与污垢结合力方面比较精确，因此得到普遍的应用［1-7］。HDICKE等［1］采用分子动力学（MD）方法模拟了聚合物添加剂抑制水垢的过程，ZHANG等采用分子动力学方法模拟了低聚物抑制剂抑制方解石成形的过程2］,TANG等团利用MD方法模拟了硫酸盐和添加剂影响碳酸钙结晶的过程，张曙光等利用分子动力学方法模拟了水溶性聚合物与方解石晶体相互作用【7】。COETZEE等【16研究表明，Zn2+比Ni2+，Ag+，Mg2+，Cu2+，Fe2+等更能促进文石生长，溶液中Zn2+和Ni2+自由离子含量的不同导致了碳酸钙晶相的不同。CHENG等研究了Zn2+在CaCO3表面吸附的微观结垢过程【17】。KATZ【18】认为，Fe3+比Fe2+更能抑制方解石的生长，并给出2种离子在pH值下的抑制作用。而HERZOG等［鸟认为，Fe2+比Fe3+对方解石生长的抑制作用更强。TAKASAKI等【20研究了Fe3+抑制方解石生长的过程以及温度、pH值、初始过饱和度等对其的影响。RAKITIN等［21］使用电化学方法研究了阴离子阻碍污垢层形成的作用是按照Cl-，NO3-，Br-，I-的顺序递减。可见，在离子影响碳酸钙晶型方面，存在着不同的研究结果。笔者采用分子动力学方法研究了含有Fe3+的Fe2O3对Fe基体CaCO3析晶污垢的影响，通过计算Fe2O3和Fe吸附污垢的结合能大小并进行比较，同时利用径向分布函数［11］进行了计算验证。1模型构建和模拟方法1.1模拟力场分子动力学（MD）模拟部分使用的是materialstudio软件的Discover模块，模拟采用第一性原理原子间相互作用势函数，模拟力场采用的是Compass力场［12】Compass力场可以模拟金属离子、金属氧化物和金属单体，还有一些小分子和高分子。与其他的力场（CVFF，PCFF）相比，该力场测量的准确性比较大，而且对于2种晶体的混合体系，也可以采用比较合理的模型进行描述。1.2模型构建方解石属于R3-C空间群，晶胞参数是a=b=0.4988nm,c=1.7061nm,a=&=900，Y=1200，由文献［1］得，（1-10）晶面是方解石主要的生长晶面，模拟过程中建立超级晶胞的体积为L4964nmx5.1183nm乂L08572nm,Ca2+和C°；各3层，其中包含630个原子（Ca：126,C：126,0:378）。建好模型后采用Discover中SmartMinimizer的优化方式进行能量最小优化，优化后的方解石（1-10）晶面的一个超级晶胞如图1所示。

Fig.1OptimizedcrystalofcalciteFig.1OptimizedcrystalofcalciteFe属于IM-3M空间群，晶胞参数是a=b=c=0.28644nm，a=&=y=90o，取其主要的生长面（001）晶面，建立超级晶胞的体积大小基本与方解石（1-10）相同。Fe2O3属于R-3C空间群，晶胞参数是a=b=0.5035nm，c=1.372nm，a=&=90o，y=120o，（001）晶面是Fe2O3的主要生长面⑷，建立超级晶胞的体积大小也与方解石（1-10）基本相同，尽管Fe2O3属于离子型晶体，但是在分布好力场之后，为了能够对体系进行正确的计算，优化前需要删除Fe-O键。建立好模型之后采用Discover模块中的SmartMinimizer的优化方式分别对其进行能量最小优化，优化后的最终模型如图2所示。图2Fe和Fe2O3最终优化构型Fig.2FinallyoptimizedconfigurationofFeandFC2O31.3模拟方法将优化好的Fe和Fe2O3的超级晶胞分别置于方解石（1-10）晶面上，组成一个混合吸附体系模型，然后对体系进行平衡。为了避免模拟过程中周期边界的影响，在2个模型的第2层上方增加5nm真空层。在对体系进行平衡过程中，可以忽略方解石晶面的晶格振动，因此冻结方解石晶面的全部原子9-10］，只允许Fe和Fe2O3进行弛豫到最优化状态，由于考虑到研究过程中压力因素的非关键性，因此采用正则系统（NVT），模拟在Andersen恒温热浴下进行，温度设定为298K，库仑力和范德华力这2种非键作用采用0.95nm的球形截断，时间步长设定为1fs，总的模拟时间为100ps，总模拟步数为10万步。图3为100ps体系平衡后Fe吸附方解石（1-10）和Fe2O3吸附方解石（1-10）的稳定构型。

2结果与讨2.12种结,3分别吸附方解石(1-12结果与讨2.12种结,3分别吸附方解石(1-1n!|gurationsofFeandFC2O3abso判断结构体系有没有达到平衡是以温度和能量平衡为依据的［刀。以Fe2O3吸附方解石(1-10)分子动力学(MD)模拟为例，图4a)和图4b)分别给出了体系平衡阶段最后10ps的温度平衡曲线和能量平衡曲线，从图4a)中温度变化曲线可以看出，温度曲线在298K上下波动，而且波动的范围非常小，说明体系温度图4 Fe2图4 Fe2O3吸附方解石体系温度和能量波动曲线Fig.4TemperatureandenergyfluctuatingcurveofFe2O3absorbingcalcite2.2Fe和Fe2O3吸附方解石晶面的结合能当Fe和Fe2O3吸附方解石（1-10）的2个体系达到平衡稳定状态时，分别计算它们相互作用后最终体系的总能量，然后排除Fe或者是Fe2O3的影响计算方解石（1-10）晶面的单点能，排除方解石（1-10）晶面的影响计算Fe和Fe2O3各自的单点能。由固体理论可得，吸附能为负表示能够吸附，负值越多说明越容易吸附，吸附能为正表示很难吸附，定义体系的结合能气为吸附能的负值［刀。2个体系的吸附能计算公式⑵为（1）△E=气-（气+Ea）（1）式中：△E为吸附能；E为相互作用后体系的总能量；Es为相互作用后方解石晶面的单点能；Ea为相互作用后Fe或者Fe2O3的单点能。表1给出了2个体系相互作用后的总能量和单点能。根据表1的计算结果可以看出，Fe吸附方解石（1-10）的结合能为76kJ/moloFe2O3吸附方解石（1-10）的结合能为32622kJ/mol，其值远大于Fe吸附方解石（1-10）的结合能，从而说明Fe2O3吸附污垢后，污垢更加难以清理。表1Fe和Fe2O3吸附方解石相互作用后的能量值Tab.1InteractionenergyofFeandFC2O3absorbingcalcitekJ/mol构型△EEtEsEaEbFe吸附方解石-7663218502364219963154290076Fe2O3吸附方解石-3262224152425786141869250326222.3Fe和Fe2O3吸附方解石的径向分布函数为了进一步的分析Fe2O3对Fe吸附污垢的影响，计算了2个平衡体系的径向分布函数g（r）口】，图5给出了Fe吸附方解石和Fe2O3吸附方解石2个平衡体系的径向分布函数，径向分布函数g（r）指的是距离一个原子为r时找到另一个原子的概率。如图5可以看出，Fe2O3吸附方解石的径向分布函数图中特征值的波峰要比Fe吸附方解石的径向分布函数图中特征值的波峰要高，说明Fe2O3吸附方解石体系发现方解石的概率要比Fe吸附方解石体系要大，因此Fe2O3吸附方解石体系的结合能要比Fe吸附方解石体系要高，这与前面结合能的分析结果相一致。2.5「0.60.6图5Fe和Fe2O3吸附方解石晶面的径向分布函数图TOC\o"1-5"\h\z1~1 II IIIIIIL IIIIIIILII IILIILIII ILII1I 1II1I ILILI ILILLI1I ILILIILIILILIILII 1ILI ILIII I1II 1I1IILIIIILI I1I,-Fig 5LRadial distribution ◎ncQon&f席and Fe2O3?bsorbing calcite crystalI~i iri czi kO O5 I~i kO H l£：iOkOOklS O G CZ' 必CD GL L CM「7 E E 寸 寸gggCO I I 00 006 6间距r3结语采用分子动力学(MD)模拟方法模拟了Fe和Fe2O3吸附方解石(1-10)晶面的体系，通过计算2个体系相互作用后各种能量，进而推导出2个体系吸附后的结合能，结果显示含Fe3+的Fe2O3吸附方解石(1-10)晶面的结合能比Fe吸附方解石(1-10)晶面的结合能大，利用径向分布函数也验证了这一结果。证明了Fe基材料吸附污垢比较难以处理的原因是腐蚀产物Fe2O3吸附污垢的结合力比较大。本研究的结果为以后碳钢表面的除垢和防垢提供了理论依据。参考文献/References:HADICKEE,RIEGERJ,RAUIU,etal.Moleculardynamicssimulationsoftheincrustationinhibitionbypolymericadditives[J].PhysicalChemistryChemicalPhysics,1999,1(17):3891-3898.ZHANGQiuyu,RENHua,GUILLOTB,etal.Molecularsimulationofoligomerinhibitorsforcalcitescale[J].JournalParticuology,2012,10(3):266-275.TANGYongming,ZHANGFan,CAOZiyi,etal.CrystallizationofCaCQinthepresenceofsulfateandadditives:Experimentalandmoleculardynamicssimulationstudies[J].JournalofColloidandInterfaceScience,2012,377(1):430-437.陈立桥，刘伟平，陈家林，等.磁铁矿与赤铁矿型纳米晶的形貌和尺寸控制生长及其性能调变J].中国科学：化学，2011,41(7)：1220-1221.CHENLiqiao,LIUWeiping,CHENJialin,etal.Controllinggrowthofnanocrystalmorphologyandsizeandmodulatingperformanceonmagnetiteandhematite[J].ChinaScience：Chemistry,2011,41(7):1220-1221.贾长治，殷军辉，薛文星.MDADAMS虚拟样机从入门到精通[M].北京：机械工业出版社,2010.张延庆.板壳有限元计算几何理论[D].徐州:中国矿业大学,1990.ZHANGYanqing.ComputationalGeometricTheoryofFiniteElementofPlateandShell[D].Xuzhou:ChinaUniversityofMiningandTechnology1990.王文生.非线性规划在可烯物配置中的应用[C]//中国运筹学会第五届大会论文集.西安：西安电子科技大学出版社，1996：468-471.张曙光，石文艳，雷武，等.水溶性聚合物与方解石晶体相互作用的MD模拟[J].物理化学学报，2005,21(11)：1198-1204.ZHANGShuguang,SHIWenyan,LEIWu,etal.Moleculardynamicssimulationofinteractionbetweencalcitecrystalandwater-solublepolymers[J].ActaPhysico-chimicaSinica,2005,21(11):1198-1204.王世燕，张军，卢贵武，等.聚合物阻垢机理的分子动力学模拟研究[J].中国石油大学学报(自然科学版)，2007,31(5)：144-147.WANGShiyan,ZHANGJun,LUGuiwu,etal.Moleculardynamicssimulationonmechanismofscaleinhibitionofpolymer[J].JournalofChinaUniversityofPetroleum(EditionofNationalScience),2007,31(5):144-147.GUILLOTB,SATORN.Carbondioxideinsilicatemelts:Amoleculardynamicssimulationstudy[J].GeochimicaEtCosmochimicaActa,2011,75(7):1829-1857.LINQ,PLAGENTZV,SCHAFERD,etal.Remediationoftrichloroethyleneandmonochlorobenzene-contaminatedaquifersusingtheORC-GAC-Fe-CaCO3system:Volatilization,precipitation,andporositylosses[J].Pedosphere,2007,17(1):109-116.杨虹，卢贵武，张雪粉，等.油、水、活性剂混合体系在方解石表面吸附的微观机理[J].中国石油大学学报(自然科学版)，2009,33(1)：131-135.YANGHong,LUGuiwu，ZHANGXuefen，etal.Microscopicmechanismofadsorptionofoil-water-surfactantmixtureoncalcitesurface[J].JournalofChinaUniversityofPetroleum(EditionofNationalScience),2009,33(1):131-135.FAKeqing,NGUYENAV,MILLERJD.InteractionofcalciumdioleatecollectorcolloidswithcalciteandfluoritesurfacesasrevealedbyAFMforcemeasurementsandmoleculardynamicssimulation[J].InternationalJournalofMineralProcessing,2006,81(3):166-177.SPAGNOLID,KERISITS,PARKERSC.Atomisticsimulationofthefreeenergiesofdissolutionofionsfromflatandsteppedcalcitesurfaces[J].JournalofCrystalGrowth,2006,294(1):103-110.任世伟，杜雄杰，侯春旺，等.硅团簇拉伸/压缩性能的分子动力学模拟[J].河北科技大学学报，2010，31(5)：399-402.RENShiwei,DUXiongjie,HOUChunwang，etal.Moleculardynamicssimulationoftension/compressionpropertiesofSicluster[J].JournalofHebeiUniversityofScienceandTechnology,2010,31(5):399-402.张娟，王春芳，宋秭龙，等.竹炭对亚甲基蓝的吸附平衡和动力学研究[J].河北科技大学学报，2012,33(5)：416-421.ZHANGJuan，WANGChunfangSONGZilong,etal.Adsorptionequilibriumandkine