周帅帅1a,1b,1d，李 静1a，杨 浩1b，乔聪震1d，卢春喜2，张经纬1a,1c
（1.河南大学 a.纳米杂化材料应用技术国家地方联合工程研究中心;b.化学化工学院；c.纳米功能材料及其应用河南省协同创新中心； d.河南省催化反应工程中心，河南 开封 475000；2.中国石油大学(北京) 化学工程与环境学院,北京 昌平 102249）
DOI:10.13732/j.issn.1008-5548.2022.02.011
收稿日期： 2021-10-06，修回日期：2021-10-19，在线出版时间：2022-02-17 。
基金项目：NSFC-联合基金项目，编号：U2004176；河南省自然科学基金项目，编号：202300410063。
第一作者简介：周帅帅(1989—)，男，讲师，博士，研究方向为过程强化。E-mail：10150173@vip.henu.edu.cn。
通信作者简介：
卢春喜(1963—)，男，教授，博士，博士生导师，研究方向为过程强化。E-mail：lcx725@sina.com；
张经纬(1975—)，男，教授，博士，博士生导师，研究方向为纳米材料工程。E-mail：jwzhang@henu.edu.cn。
摘要： 针对液相法石墨烯生产关键设备搅拌釜反应器,采用实验和模拟的方法对其内部多相复杂流场特性和规律进行研究。采用流体容积模型两相模型实现搅拌釜内多相流场的预测，并实验验证模型的适用性；采用计算流体动力学和离散元法建模耦合方法，计算片状颗粒在搅拌釜内的运动特性。结果表明：剪切速率在桨叶边缘及搅拌槽边壁附近区域较高，而颗粒主要集中于搅拌槽底部靠近边壁区域，并在底部形成颗粒堆积现象；单颗粒运动轨迹表明，非球形颗粒主体流动沿着流线方向，受底部颗粒堆积的影响，颗粒在运动过程中会在堆积区域停留较长时间；单颗粒受力结果分析发现，颗粒在运动过程中受搅拌桨和边壁的碰撞作用，会受到较大的脉冲力，而该合力持续时间较短，是颗粒运动轨迹改变的主要作用力。
关键词： 液相法；石墨烯；搅拌釜反应器；非球形颗粒；计算流体动力学；离散元建模
Abstract:The methods of experiment and simulation were used to reveal characteristcis of multi-phase complex flow field inside the stirred tank reactor,which was the key unit in graphere production by using liquied-phase method.Volume of fluent model of two-phase was used to predict the multi-phase flow field in the stirred tank,and suitability the model was verified by experimental results.Computational fluid dynamics-discrete element modeling coupling method was used to calculate the movement characteristics of the non-spherical particles in the stirred tank.The results show that the shear rate is much higher in the region near the edge of blades and the reactor wall,the particles are mainly concentrated in the region near the reactor wall at the bottom of the reactor,and the phenomenon of particle accumulation is formed at the bottom.The single particle trajectory shows that the main flow of non-spherical particles is along the streamline,and is affected by the accumulation of particles at the bottom.The particles will stay in the accumulation area for a long time during the movement.The forces exerted on the single particles show that the particles will be subjected to a larger pulse force due to the collision of the blades and the reactor wall during the movement,and the resultant force lasts for a short time,which is the main force for the change of the particle trajectory.
Keywords:liquid phase method;graphene；stirred tank reactor；non-spherical particle;computational fluid dynamics;discrete element modeling
参考文献（References）：
[1]陈玉华.新型清洁能源技术[M].北京：知识产权出版社，2019.
[2]朱宏文,段正康,张蕾,等.氧化石墨烯的制备及结构研究进展[J].材料科学与工艺,2017,25(6):82-88.
[3]李渝,胡志波,文兴青,等.石墨烯材料制备技术与创新及发展趋势[J].中国非金属矿工业导刊,2021(2):19-23.
[4]HUANG J,ZHANG J.Scalable production of few layered neat graphene by soft ball-microsphere rolling transfer[J].Carbon,2019,154:402-409.
[5]帅骁睿,张鹏程,张正卿,等.氧化还原法石墨烯绿色制备技术研究进展[J].化工管理,2020(29):31-34.
[6]SUTTER P.Epitaxial graphene:how silicon leaves the scene[J].Nature Materials，2009，8:171-172.
[7]涂溶,胡志颖,章嵩,等.一种Si/SiC/石墨烯材料的制备方法:110323126B[P].2021-09-14.
[8]白瑞,李振宇,慕苗,等.石墨烯@Fe3O4纳米复合材料的制备及表征[J].化工科技,2021,29(4):16-18.
[9]张岳,慕晓刚，吕晓猛，等.石墨烯及化学改性石墨烯的胶状悬浮物合成研究进展[J].化学推进剂与高分子材料,2017(1):40-44.
[10]ARAO Y,MIZUNO Y,ARAKI K,et al.Mass production of high aspect ratio few layer graphene by high speed laminar flow[J].Carbon,2016,102:330-338.
[11]RACCICHINI R,VARZI A,PASSERINI S,et al.The role of graphene for electrochemical energy storage[J].Nature Materials,2015,14:271-279.
[12]CAI X,LUO Y,LIU B,et al.Preparation of 2D material dispersions and their applications[J].Chemical Society Reviews,2018,47:6224-6266.
[13]YI M,SHEN Z G.A review on mechanical exfoliation for the scalable production of graphene[J].Journal of Materials Chemistry：A,2015,3(22):11700-11715.
[14]MO J Y,GAO Z M,BAO Y Y,et al.Suspending a solid sphere in laminar inertial liquid flow-experiments and simulations[J].AIChE Journal,2015,61(4):1455-1469.
[15]莫君媛.固液搅拌槽内颗粒运动行为及悬浮机理的研究[D].北京：北京化工大学,2015.
[16]侯家鑫.立式换热管内椭球形颗粒运动和传热特性研究[D].湘潭：湘潭大学,2020.
[17]SHAO T,HU Y,WANG W,et al.Simulation of solid suspension in a stirred tank using CFD-DEM coupled appraoch[J].Chinese Journal of Chemical Engineering,2013,21 (10):1069-1081.
[18]李文华,白亮,虞上长.电除尘器流场、颗粒场及电场的数值模拟研究[J].发电技术,2021,42(3):336-342.
[19]侯淑晴,王智勇,王伟.基于Hertz-Mindlin接触模型的摩擦界面颗粒剪切膨胀特性研究[J].合肥工业大学学报(自然科学版),2018,41(12):1601-1605.
[20]文愿运,刘马林,刘荣正,等.颗粒离散单元法数值模拟与典型实验对比研究[J].中国粉体技术,2015,21(3):1-5.
[21]KATAGIRI J,HAYASHI N,KOYANAKA S.Elementary-volume-scale simulations of inertial flow in sphere pack:improvement of Di Felice drag model in high porosity[J].Materials Transactions,2020,61 (5):1026-1031.
[22]MEI R.An approximate expression for the shear lift force on a spherical particle at finite Reynolds numbers[J].International Journal of Multiphase Flow,1992,18 (1):145-147.
[23]OESTERLEE B,BUI D T.Experiments on the lift of a spinning sphere in a range of intermediate Reynolds numbers[J].Experiments in Fluids,1998,25(1):16-22.
[24]张永震.搅拌釜式生物反应器的计算流体力学模拟[D].天津：天津大学,2005.
[25]舒雷.基于CFD的搅拌反应器流场模拟和设计优化研究[D].重庆：西南大学,2020.
