前混合水射流中丸粒速度的数值模拟与PIV试验研究.pdf

研究论文RESEARCH引文格式：朱泽华，刘兰荣，马泳涛，等。前混合水射流中丸粒速度的数值模拟与PIV试验研究.航空制造技术,2 0 2 3,6 6(13):8 0-8 6.ZHU Zehua,LIU Lanrong,MA Yongtao,et al.Numerical simulation and PIV experimental study of shot velocity in pre-mixedwaterjetJj.Aeronautical Manufacturing Technology,2023,66(13):80-86.前混合水射流中丸粒速度的数值模拟与PIV试验研究朱泽华,刘兰荣”,马泳涛”，卢春生”,张彬?（1.郑州大学河南省资源与材料工业技术研究院，郑州450 0 0 1；2.郑州大学机械与动力工程学院，郑州450 0 0 1;3.School of Civil and Mechanical Engineering,Curtin University,Perth WA 6845,Australia)【摘要前混合水射流中丸粒速度是决定混合射流性能的关键因素。采用CFD-DEM双向耦合方法构建了混合射流的三相仿真模型，气液相采用VOF模型并考虑离散相体积分数的影响，丸粒固态相采用拉格朗日坐标下的离散元方法进行建模，通过专用的耦合接口实现双向耦合，得到丸粒在出口确定位置范围内速度仿真值。通过专用的PIV试验台测出相应位置丸粒的速度值。通过对比，仿真结果与测试结果的差值在6%以内，证明了仿真模型的正确性。研究结果表明，喷嘴入口压力和喷嘴出口直径越大，丸粒速度越大；混合比（丸粒体积分数）、喷嘴长度和喷嘴锥角越大，丸粒速度越小，对后续制定改性工艺方案提供了理论基础。关键词：前混合水射流；改性；CFD-DEM；PI V；丸粒速度Numerical Simulation and PIV Experimental Study of Shot Velocity in Pre-Mixed WaterjetZHU Zehua,LIU Lanrong,MA Yongtao,LU Chunsheng,ZHANG Bin?(1.Henan Provincial Industrial Technology Research Institute of Resources and Materials,Zhengzhou University,Zhengzhou 450001,China;2.School of Mechanical and Power Engineering,Zhengzhou University,Zhengzhou 450001,China;3.School of Civil and Mechanical Engineering,Curtin University,Perth WA 6845,Australia)ABSTRACT The shot velocity is the key factor to determine the performance of pre-mixed waterjet.The three-phasesimulation model of a mixed jet was constructed by the CFD-DEM bidirectional coupling method.The VOF model wasused for a gas-liquid phase and the volume fraction influence of a discrete phase was considered.The shot solid phasewas modeled by the discrete element method under Lagrangian coordinates.The two-way coupling was realized througha special coupling interface,and the simulation value of shot velocity in a fixed position range of the exit was obtained.By using the particle image velocity(PIV)test system,the shot velocity at the corresponding position was measured.Bycomparison,the difference between simulation and testing results was less than 6%,implying correctness of the simulationmodel.The simulation results show that the larger the nozzle inlet pressure and nozzle outlet diameter,the greater the shotvelocity;The larger the mixing ratio(volume fraction of shots),nozzle aspect ratio and nozzle cone angle,the smaller theshot velocity.This provides a theoretical basis for the subsequent formulation of a modification process.Keywords:Pre-mixed waterjet;Modification;CFD-DEM;PIV;Shot velocityD0I:10.16080/j.issn1671-833x.2023.13.080*基金项目：国家自然科学基金（U1804254）；河南省高等学校重点科研项目（2 1A460024）。80航空制造技术2023年第6 6 卷第13期RESEARCH研究论文前混合水射流喷丸表层改性技术是一种新兴的金属工件表面处理技术，其改性机理是将预混合的丸粒通1计算模型过高压水进行加速，高速喷射到金属工件表面，从而使工件产生塑性变形,形成残余压应力层-2 。在前混合水射流中,射流束能量场的大小及分布直接影响改性效果，而射流束能量场与丸粒的速度密切相关3-4,因此研究前混合水射流束中丸粒的速度分布，对揭示前混合水射流表层改性机理，提高表层改性能力及制定改性工艺参数具有重要意义。Basha等5 采用Fluent软件中的VOF模型和DPM模型对后混合磨料水射流中喷嘴内部气液固三相的流动特性进行了数值模拟,基于VOF模型模拟了水-气两相流,基于DPM 模型模拟了颗粒轨迹。Long等6 采用了相同的耦合策略并研究了颗粒尺寸、密度、形状因子对后混合射流喷嘴出口颗粒速度的影响。Liu等也采用VOF+DPM的策略,并求解了超高速磨料水射流外流场的速度分布。Qiang等8 采用Fluent软件中欧拉模型和DPM模型模拟了后混合射流系统中切割头内部的射流流动和颗粒混合过程，分析了颗粒入口角度、人口位置、聚焦管收敛角对流体的加速过程的影响。在PIV测速方面，Thongkaew等9 结合PIV和LIF（激光诱导荧光）技术，对后混合磨料射流中的粒子轴向速度进行试验研究。Zelenak等10 也采用相同技术对后混合磨料水射流的颗粒速度进行试验研究,并通过PTV算法分析单个粒子的运动。章文峰等I基于PIV技术和图像处理技术对前混合水射流中颗粒速度进行试验研究，同时得到了磨料的速度和位置信息，发现在14mm喷距范围内，磨料速度先增大后减小，PIV试验采用锥直形喷嘴和球形陶瓷丸粒,喷嘴出口直径为6 mm，丸粒密度为2.7 g/cm。综上所述，上述的研究多为单一的通过仿真模拟或者PIV测速试验研究液固混合射流中丸粒的速度，研究对象多为后混合水射流。后混合水射流是指在水射流形成之后，利用混合腔内负压引入磨料，最终形成后混合水射流，因此，后混合水射流中磨料与水混合不均匀，能量交换效率低。数值模拟时对于颗粒的建模多采用DPM模型,并且只采用了单相耦合,未能考虑水丸之间的相互作用。本文通过耦合计算流体力学软件Fluent和离散元软件EDEM,对气液固三相射流束进行了数值模拟研究。结合PIV测试技术和图像处理技术对前混合水射流束的铸钢丸粒速度进行了试验测试，验证了仿真模型的可信度。基于CFD-DEM耦合方法探究喷嘴入口压力、混合比（丸粒体积分数）喷嘴结构参数对丸粒速度的影响，为制定工艺参数和优化喷嘴结构提供技术基础。1.1CFD模型的构建1.1.1数学模型高压射流通常处于流状态，所以系统除了求解连续性方程、N-S方程等基本控制方程外，还需要遵循额外的湍流输运方程12-13。采用雷诺平均法模拟湍流流动时，会使雷诺时均N-S方程（式(1）产生新的雷诺应力未知项，即puu，为了封闭方程组，本文采用RealizableK-e瑞流模型求解该项。U;-(u,)(u,uj)tX式中,p为流体密度；为动力黏性系数;p为压力；u,和u,为流体速度；t为时间；x、z 为位置；i和j的取值范围为1,2,3；S,为动量守恒方程的广义源项。此外，前混合水射流包含气液固3相，属于多相流动问题。在该多相流系统中，可以将气液两相看成是连续相,将丸粒看成是离散相5-6.14。离散相的求解过程在1.2 节中介绍。关于连续相的求解，对气液两相利用连续介质假设理论，在Euler坐标系下，采用有限体积法对其进行模拟。基于该理论，本文采用Fluent中的VOF模型对前混合水射流的连续相进行模拟，需要注意的是，在VOF模型中,连续相中的气液两相共享一组动量方程，但是考虑到离散相体积分数的影响，需要对VOF的连续性方程和N-S方程进行修改，修改后的连续性方程和N-S 方程15 分别为fu0t(u)f式中，&为网格孔隙率；F为表面张力；u为网格内流体速度；Pr为网格内流体密度；亏为黏性应力张量；g为重力加速度；Fr为流体颗粒相互作用项的反作用力。1.1.2物理模型及边界条件选取水射流技术常用的锥直形喷嘴，计算区域包括喷嘴内部和喷嘴外部直径2 0 mm、长度6 0 mm的圆柱区域。结构化网格模型及边界条件如图1所示，考虑到喷嘴壁面，在靠近喷嘴壁面处生成边界层网格,并在靠近气液两相分界处加密网格,用于改善VOF两相界面追踪效果。将Inlet设置为压力人口边界条件，Wall和Walll设置为无滑移速度壁面边界条件，Wall2和Outlet设置为压力出口边界条件，压力值设为0.10 MPa（大气压）。1.2DEM模型的构建在前混合水射流多相流系统中,对于离散相的求2023年第6 6 卷第13期航空制造技术8 1卫u.ujS,(1)X;(2)(uu)f（u,u,)F,r8pl Jor(3)前混合水射研究论文RESEARCH解，可以在拉格朗日坐标系下基于离散元（DEM）方法进行,该方法对所有颗粒进行追踪16 ,考虑颗粒与颗粒、颗粒与边界的相互作用。基于该方法,本文采用EDEM对前混合水射流的离散相进行模拟，并采用Hertz-Mindlin无滑动接触模型求解颗粒单元之间的接触力。在EDEM设置过程中，瑞利时间步取17%，总仿真时间为0.0 1s。试验中采用铸钢丸和硬质合金喷嘴，材料参数如表1所示。1.3相间作用力由前面的分析可知,对于前混合水射流多相流系统,通过CFD-DEM耦合方法在欧拉坐标系中描述连续相的运动，在拉格朗日坐标系下描述离散相的运动，通过相间作用力来实现连续相和离散相的双向耦合。在前混合水射流系统中,所考虑的颗粒-流体相互作用力如式(4)所示 15。FrFF,FumFFsurFMepfvm式中，Fr是颗粒k受到的所有流体-颗粒之间的相互作用力；F。是电力，通过自由流曳力模型求解17 ；Fs，是压力梯度力；Fvm是虚拟质量力；Fg是Basset力；Fsar是Saffiman 升力;FMag 是Magnus升力。2PIV测试试验2.1PIV试验装置介绍前混合水射流PIV测速专用试验台如图2 所示，包括两部分，一是前混合水射流喷射系统，用于提供前混合水射流；二是PIV测速系统（美国TSI公司）。PIVWalll测速的工作过程为：（1）对相机进行标定，计算出校准因子，得到图像坐标与真实坐标之间的转换关系；（2）喷射射流，相机与脉冲激光同步，在两个连续激光脉冲期间捕获射流图像，完成图像采集；（3）通过计算机系统对拍摄的一系列图片进行分析处理,进而得到速度矢量及速度相关量。相机的分辨率为2 36 0 17 7 6 像素，帧率为16 顿/s，最小跨帧时间为19 5ns；激光器的脉冲频率为15Hz,脉冲能量为 2 0 0 mJ。2.2PIV图像处理PIV试验的拍摄范围为从喷嘴出口到距离喷嘴出口约6 5mm的范围，通过PIV测速试验采集的射流原始图像的放大图如图3（a）和（b）所示。原始图像可能存在少量由于激光亮度过大导致的过度曝光点或丸粒边界与射流之间的噪音，本文通过Matlab图像处理(4)技术对原始图像进行预处理，从而削弱上述因素。最后通过基于PTV18算法的Matlab程序对预处理后的A顿、激光器相机数据控制及Outlet处理系统同步器图2 PIV测速试验设备Fig.2PIV velocity measurement equipment流喷射系统标定系统mm(a)A帧Inlet-Fig.1 Schematic diagram of grid and boundary conditionsWall2Wall图1网格和边界条件示意图mm表1材料参数设置Table1Parametersettingofmaterials密度/材料(gcm)丸粒7.5喷嘴7.882航空制造技术2023年第6 6 卷第13期(b)B剪切模量1 mm泊松比形状直径/mmGPa0.3800.2580（c）速度矢量图球形0.25图3PIV原始图像及速度矢量图Fig.3 PIV original images and velocity vector diagramRESEARCH研究论文B帧图像进行分析处理，得到丸粒的速度及位置信息，速度矢量图如图3（c）所示。经过上述处理后，得到一对图像中丸粒的速度相关信息。需要注意的是，丸粒的识别主要是通过给定的丸粒与背景之间的灰度值值和来判定的,但是在同一幅图像中，每个丸粒与背景的灰度值差值有所不同,所以可能会误识别丸粒，从而出现错误的速度值，因此需要对异常值进行剔除，本文采用基于拉依达准则的Matlab程序对速度异常值进行剔除19 。2.3试验内容通过对比三因素（喷嘴入口压力、混合比、喷嘴长度）三水平下的仿真和PIV测试结果，从多角度验证仿真模型的准确性。考虑到试验安全性及测试条件,喷嘴人口压力、混合比、喷嘴长度均限制在一定范围内。评价指标包括两个，一是外流场中丸粒的速度变化趋势，用于定性地说明二者的符合程度；二是距离喷嘴出口10mm30mm、50 m m截面处丸粒的平均速度，分析仿真值与试验值的误差大小，用于定量说明仿真模型的准确性。试验时在磨料罐中添加8 0 0 mL丸粒,每次拍摄100组图片。PIV验证试验的试验参数如表2 所示。2.4试验结果分析不同喷嘴入口压力、混合比和喷嘴长度下，丸粒速度在外流场中的分布情况分别如图46 所示。可以看出，在任一喷嘴入口压力、混合比和喷嘴长度下，丸粒速度变化趋势的仿真结果和试验结果基本相同,在6 0 mm的靶距范围内，两种方法得到的丸粒速度基本保持不变。需要注意的是,本文的结论和文献11 结论不同，这可能是由于喷嘴出口直径和丸粒密度不同引起的。进一步地，通过式（5）对丸粒仿真速度的误差率进行分析，分别计算不同参数下距离喷嘴出口10 mm、表2 PIV验证试验的参数设置Table 2 Parameter setting in PIV verification experiments试验喷嘴入口混合比/%序号压力/MPa14263844546474849430mm、50 m m 截面处丸粒的仿真速度误差率。通过分析可知，在4MPa、6 MPa、8 MPa 3组压力下丸粒仿真速度的误差率绝对值最大分别为3.53%、2.10%、2.9 7%；在2.41%、4.47%、6.10%3组混合比下，丸粒仿真速度的误差率绝对值最大分别为4.6 3%、4.8 4%、4.2 4%；在2 6mm、31m m、41m m 3组喷嘴长度下,丸粒仿真速度的误差率绝对值最大分别为4.8 7%、5.2 0%4.7 2%。3组结果均表明在距离喷嘴出口距离为10 mm、30 m m、50 m m处，丸粒的仿真结果与试验结果基本相同。Vave.sVave.WVave.e式中，w为丸粒仿真速度的误差率；Vave.s为丸粒仿真平均速度;vave。为丸粒试验平均速度。在计算平均速度时，8075(-s.u)/率*47065605550454009590(-S.)/*858075706560550喷嘴长度/喷嘴出口mm直径/mm4.47464.47464.47462.41414.47416.10414.47264.47314.4741(5)仿真值试验值102030405060距离喷嘴出口距离/mm(a)4 MPa仿真值试验值102030405060距离喷嘴出口距离/mm(b)6 MPa105100(-s.U)/率*1.61.61.61.61.61.61.21.21.2959085807570650图4不同喷嘴入口压力下试验与仿真的结果对比Fig.4 Comparison of experimental and simulation results underdifferent nozzle inlet pressures2023年第6 6 卷第13期航空制造技术8 3仿真值试验值上102030405060距离喷嘴出口距离/mm(c)8 MPa研究论文RESEARCH9085(r-s.w)/取8075706560555008580(-s.w)/率7570656055504508075(r-s.W)/*706560555045400图5不同混合比下试验与仿真的结果对比Fig.5Comparison of experimental and simulation results underdifferent mixing ratios由于在某一时刻落在指定截面上的丸粒数量过少甚至为0,所以将该截面前后1mm范围区域的丸粒均视为该截面上的丸粒。综上所述，通过对不同喷嘴人口压力、混合比和喷嘴长度下的仿真速度及试验速度大小、误差进行对比分析，均验证了仿真模型的正确性。3前前混合水射流中丸粒速度影响因素研究丸粒速度是影响前混合水射流改性效果的重要参数,在进行影响因素分析时,以经过PIV试验验证的仿真模型为基础,采用单因素分析法研究了喷嘴入口压力、混合比、喷嘴结构参数对前混合水射流中丸粒速度的影响规律,其中各参数的取值范围以涵盖实际前混合84航空制造技术2023年第6 6 卷第13期8580(r-s.W)/取*75706560仿真值55试验值5045102030405060距离喷嘴出口距离/mm(a)2.41%仿真值试验值102030405060距离喷嘴出口距离/mm(b)4.47%仿真值试验值102030距离喷嘴出口距离/mm(c)6.10%仿真值试验值1101020304050距离喷嘴出口距离/mm(a)26 mm8580(1-S.U)/氧7570656055504508580(1s.u)/率球75706560555045405060160仿真值试验值102030405060距离喷嘴出口距离/mm(b)31 mm仿真值试验值0102030405060距离喷嘴出口距离/mm(c)41 mm图6 不同喷嘴长度下仿真与试验的结果对比Fig.6 Comparison of experimental and simulation results underdifferent nozzle lengths水射流喷丸表层改性技术为基础。3.1喷嘴入口压力由于前混合水射流能量利用率高，工作压力通常在40MPa以内，喷嘴入口压力对丸粒最大速度的影响规律如图7 所示，丸粒最大速度随着喷嘴人口压力的增大而增大。由液体孔口出流理论可知，喷嘴人口压力增大使水的速度增大，高速水对丸粒起到加速作用，进而使丸粒速度增大。压力是影响丸粒速度的重要因素，通过调节压力来调节丸粒速度是较为简单的方式,因此对所测数据进行拟合，得到了在上述约束条件下丸粒最大速度mx.p与喷嘴入口压力p的关系式为mx.,=40.56/P(6)4结论研究论文RESEARCH3.2混合比确定，喷嘴出口直径对丸粒最大速度的影响规律如图在改性工艺中,混合比取值通常在10%以内，混合10所示，丸粒最大速度随着喷嘴出口直径的增大而增比对丸粒最大速度的影响规律如图8 所示，丸粒最大大。喷嘴出口直径从1mm增大到2.5mm，丸粒最大速速度随着混合比的增大而减小。混合比从1%增大到度从2 19 m/s增到2 40.33m/s，增加了9.7 4%。喷嘴出10%，丸粒最大速度从2 2 6.9 6 m/s降到2 0 1.46 m/s，减小口直径越大，在相同混合比的情况下，丸粒与壁面碰撞了11.2 4%。混合比越大，单位时间内高速水需要加速的概率就越小，速度损失就越小。的丸粒数量就越多，水的负载变大，这会使水的速度降3.5喷嘴锥角低，进而导致丸粒速度降低。在制定改性工艺时，若选在改性工艺中，为提高丸粒速度的同时减小喷嘴磨取较小的混合比，可得到较高的丸粒速度，但是单位时损，锥角不宜过大或过小，通常在30 左右，喷嘴锥角对间内打击至工件表面的丸粒数量较少，进而影响改性丸粒最大速度的影响规律如图11所示，丸粒最大速度效果。随着锥角的增大而减小。锥角从10 增大到6 3,丸粒3.3喷嘴长度最大速度从2 2 8.8 5m/s降到2 18.38 m/s，减小了4.58%。喷嘴长度对丸粒最大速度的影响规律如图9 所示，丸粒在喷嘴中主要是沿着喷嘴向前流动，锥角越大，丸丸粒最大速度随着圆柱段长度的增大而减小。圆柱段粒与壁面碰撞的概率就越大，速度损失就越大。长度从18 mm增大到6 6 mm，丸粒最大速度从2 2 9.8 3m/s降到2 11.7 5m/s，减小了7.8 7%。圆柱段长度越长，水和喷嘴壁面的摩擦距离越长,并且丸粒与壁面的碰撞概率增大，这都会导致丸粒速度减小。3.4喷嘴出口直径喷嘴出口直径通常根据改性零部件的具体尺寸来270255仿真值.-拟合曲线(r-s.)/取鲁2402252101951801651501351201012162025喷嘴人口压力/MPa图7 喷嘴入口压力对丸粒最大速度的影响规律Fig.7Effect of nozzle inlet pressure on the maximum velocity of shot230225(r-s.W)/率来22021521020520051 1.62.5图8 混合比对丸粒最大速度的影响规律Fig.8 Effect of mixing ratio on the maximum velocity of shot（1）基于CFD-DEM耦合方法模拟了前混合水射流喷嘴内外流场的三相流动，在考虑颗粒碰撞作用的前提下，实现了颗粒与流体之间相互交换动量、能量的双向耦合。234(-s.w)/单普*230226222218214210。324046.3混合比/%182127图9 喷嘴长度对丸粒最大速度的影响规律Fig.9Effect of nozzle lengths on the maximum velocity of shot245(r-s.U)/率鲁240235230225220215101图10 喷嘴出口直径对丸粒最大速度的影响规律Fig.10 Effect of nozzle outlet diameter on the maximum velocity of shot2023年第6 6 卷第13期航空制造技术8 540圆柱段长度/mm1.21.6出口直径/mm662.02.5研究论文RESEARCH230(r-s.U)/取半普*2282262242222202181016图11喷嘴锥角对丸粒最大速度的影响规律Fig.l1Effect of nozzle cone angle on the maximum velocity of shot（2）PI V 测速试验对仿真模型的验证结果表明，通过改变压力、混合比和长径比，在6 0 mm靶距范围内，仿真结果与试验结果吻合度较好，验证了该仿真模型的正确性。（3）前混合水射流中丸粒最大速度随着喷嘴入口压力、喷嘴出口直径的增大而增大，随着混合比、喷嘴长度、喷嘴锥角的增大而减小。改性效果不仅和丸粒速度相关,还和丸粒性质、分布相关,应综合考虑多方因素，制定合理的工艺路线，得到最优的改性效果。参考文献1马泳涛，翟林邦，朱泽华，等.水射流表层改性技术研究进展.航空制造技术，2 0 2 1,6 4(7):2 8-35.MA Yongtao,ZHAI Linbang,ZHU Zehua,et al.Research review onsurface modification with water jet peening technologyJJ.AeronauticalManufacturing Technology,2021,64(7):28-35.2 SADASIVAM B,HIZAL A,AROLA D.Abrasive waterjetpeening with elastic prestress:A parametric evaluationJJ.InternationalJournal of Machine Tools and Manufacture,2009,49(2):134-141.3ANASIEWICZ K,KUCZMASZEWSKI J.Effect of abrasivewaterjet peening surface treatment of steel plates on the strength of single-lapadhesive jointsJj.Advances in Science and Technology Research Journal,2017,11(3):270276.4 WANG F C,ZHOU D P,XU Q W,et al.Mathematicalmodel of rock stress under abrasive slurry jet impact based on contactmechanicsJ.International Journal of Rock Mechanics and MiningSciences,2018,107:1-8.5 BASHA A T,ANNONI M,MONNO M,et al.Investigationof the hydrodynamic characteristics of abrasive water jet cutting headJ.International Journal of Machining and Machinability of Materials,2013,14(1):105.6 LONG X P,RUAN X F,LIU Q,ET AL.Numerical investigationon the internal flow and the particle movement in the abrasive waterjetnozzleJJ.Powder Technology,2017,314:635-640.7 LIU H,WANG J,KELSON N,et al.A study of abrasive waterjetcharacteristics by CFD simulationJJ.Journal of Materials ProcessingTechnology,2004,153-154:488493.86航空制造技术2023年第6 6 卷第13期8 QIANG Z R,WU M P,MIAO X J,ET AL.CFD research onparticle movement and nozzle wear in the abrasive water jet cutting head.The International Journal of Advanced Manufacturing Technology,2018,95(9):40914100.9 THONGKAEW K,WANG J.An experimental study of theparticle velocities in abrasive waterjetsJ.International Journal of AbrasiveTechnology,2018,8(2):147-156.10 ZELENAK M,FOLDYNA J,LINDE M,et al.Measurementand analysis of abrasive particles velocities in AWSJJ.ProcediaEngineering,2016,149:77-86.2540锥角/（）6311章文峰，卢义玉，汤积仁,等.基于PIV技术的磨料水射流中固体磨料粒子速度分布实验研究1.振动与冲击，2 0 16,35(8):159-165,186.ZHANG Wenfeng,LU Yiyu,TANG Jiren,et al.Experimental studyof velocity distribution of solid abrasive particles in abrasive water jetsbased on PIV technologyJJ.Journal of Vibration and Shock,2016,35(8):159165,186.12 CAI C Z,HUANG Z W,LI G S,et al.Particle velocitydistributions of abrasive liquid nitrogen jet and parametric sensitivityanalysisJJ.Journal of Natural Gas Science and Engineering,2015,27:1657-1666.13王福军.计算流体动力学分析：CFD软件原理与应用M.北京：清华大学出版社，2 0 0 4.WANG Fujun.Computational fluid dynamics analysis:Principleand application of CFD softwareM.Beijing:Tsinghua University Press,2004.14】岳湘安.液-固两相流基础M.北京：石油工业出版社，1996.YUE Xiangan.Fundamentals of liquid-solid two-phase flowM.Beijing:Petroleum Industry Press,1996.15 WU L,GONG M,WANG J T.Development of a DEM-VOFmodel for the turbulent free-surface flows with particles and its applicationto stirred mixing systemJJ.Industrial&Engineering Chemistry Research,2018,57(5):1714-1725.16胡国明，颗粒系统的离散元素法分析仿真：离散元素法的工业应用与EDEM软件简介M.武汉：武汉理工大学出版社，2 0 10.HU Guoming.Analysis and simulation of granular system bydiscrete element method using EDEMMJ.Wuhan:Wuhan University ofTechnology Press,2010.17 NOROUZI H R,ZARGHAMI R,SOTUDEH-GHAREBAGHR,et al.Coupled CFD-DEM ModelingM.Iran:Wiley,2016.18 SHOKRI R,GHAEMI S,NOBES D S,et al.Investigation ofparticle-laden turbulent pipe flow at high-Reynolds-number using particleimage/tracking velocimetry(PIV/PTV)JJ.International Journal of Multi-phase Flow,2017,89:136-149.19王岩，隋思涟.试验设计与MATLAB数据分析M.北京：清华大学出版社，2 0 12.WANG Yan,SUI Silian.Design and MATLAB data analysisM.Beijing:Tsinghua University Press,2012.通讯作者：刘兰荣,实验师，硕士,研究方向为抗疲劳制造技术和表层改性技术。(责编阳光）