喷嘴结构对高压水射流影响及结构参数优化设计

喷嘴结构对高压水射流影响及结构参数优化设计收稿日期：2015-11-20基金项目：国家自然科学基金青年科学基金项目（E031303）第二炮兵工程大学科研基金青年项目（2015QNJJ034）作者简介：韩启龙（1979），男，甘肃宁县人，硕士，副教授，E-mail：longfeng.061106马洋（通讯作者），男，湖南澧县人，博士，讲师，E-mail：mldy0612韩启龙，马洋(第二炮兵工程大学 动力工程系， 陕西 西安 710025)摘要：喷嘴是产生高压水射流的关键部件，其结构形式对射流动力学性能有很大影响。以圆柱形喷嘴为对象，进行喷嘴结构对高压水射流的影响分析及结构参数优化设计。采用两相流计算流体力学模型进行喷嘴内外的射流流场分析。为节省计算资源，在优化设计时引入Kriging代理模型替代计算流体力学模型。分别采用改进的非劣分类遗传算法（Nondominated Sorting Genetic Algorithm, NSGA-II）和基于分解的多目标进化算法（MultiObjective Evolutionary Algorithms based on Decomposition, MOEA/D）进行单目标和多目标优化设计。研究结果表明：直线型喷嘴总体性能较优，凹型喷嘴次之，凸型喷嘴性能最差。以直线型喷嘴为设计对象，以射流初始段长度和流量为目标，得到了单目标和多目标优化设计结果，单目标优化时，两个指标较基准外形分别提高14.71%和27.56%。多目标优化时，优化得到的半锥角处于15.4, 89.8区间上。基于代理模型和进化算法的全局优化方法在进行喷嘴的优化设计时是有效的。关键词：高压水射流；喷嘴；全局优化；两相流；代理模型；MOEA/D中图分类号：V411 文献标志码：A文章编号：Influence of nozzle structure on high pressure water jet and optimization design of nozzle structure parameterHAN Qilong, MA Yang(Second Artillery Engineering University, Xian 710025, China）Abstract: Nozzle is the crucial component used to generate high pressure water jet, and its structure form has large influence on dynamic performance of high pressure water jet, so the influence of nozzle structure on high pressure water jet is analyzed, and the optimization design of nozzle structure parameter is implemented in this paper. Two phase flow computational fluid dynamics model is employed to analyze flow field. The Kriging surrogate model is used to replace the computational fluid dynamics model in the process of optimization design for reducing the computational resources. The Nondominated Sorting Genetic Algorithm (NSGA-II) and MultiObjective Evolutionary Algorithms based on Decomposition (MOEA/D) are respectively employed to carry out single and multi objective optimization design. The research results were summarized as follows. First, the general capability of line-form nozzle was the best, then the concavity-form nozzle, and the protruding-form nozzle has the worst capability. Second, the single and multi objective optimization design of line-form nozzle was implemented, in which core zone length and mass flux of water jet were taken as optimization objectives. Compared to the baseline, the two indexes increased by 14.71% and 27.56% respectively after the single objective optimization. The optimal semi-cone angle after multi objective optimization located on 15.4, 89.8. Third, the global optimization algorithm based on surrogate model and evolutionary algorithm was proved to be effective.Key words: high pressure water jet；nozzle；global optimization；two phase flow；surrogate model；MOEA/D由于具有清洗质量好、清洗速度快、绿色环保、安全性能高等优点1，高压水射流在固体发动机推进剂的清洗和切割中具有很好的应用前景2,3。喷嘴是产生高压水射流的关键部件，其结构形式对射流动力学性能和内部流场有很大影响。为了研究喷嘴结构的影响，广大学者采用不同的方法手段进行研究。刘庭成采用理论分析方法定性论述了喷嘴结构参数及其对清洗作业的影响4，以及喷嘴直径、胶管直径与清洗机压力、流量相匹配的重要性5。王洪伦6通过公式推导得到了影响喷嘴的冲击压力的几个因素。蒋大勇7利用自研的实验平台开展了高压水射流切割HTPB推进剂的实验研究，得到了适合切割推进剂的最佳喷嘴直径。为解决扇形喷嘴强度低、磨损快、能耗高的缺点，喻峰等8设计了一种新型多孔圆柱喷嘴，并通过数值模拟和试验的方法验证了该喷嘴在钢材除磷中良好的性能。在研究喷嘴构型对高压水射流流动影响时，理论分析的应用范围有限，试验研究尽管结果最为可靠，但成本太大。随着两相流动计算技术的发展，越来越多的学者选择数值仿真的方法进行研究1,9-12。这些研究主要关注某一类特定结构形式的喷嘴对射流性能的影响，但对该类型喷嘴的具体型面研究较少。还有一些学者采用正交试验11,13、遗传算法14对喷嘴的关键参数进行设计，在引入优化思想进行高压水射流喷嘴设计方面作出了很好的尝试。本文以常见的圆柱形喷嘴为对象，对比分析三种内部型面对射流流动特性的影响，并以表现最好的一种型面为基准，采用基于代理模型的全局优化算法对喷嘴的半锥角进行优化设计，得到了考虑单目标和多目标情况下的优化结果，这种优化设计方法能很好地兼顾设计效率和精度，可以为喷嘴的型面设计提供一定的参考。1 问题描述按照形状划分，喷嘴主要有圆柱形喷嘴、扇形喷嘴、异形喷嘴等，按照喷嘴的孔数划分有单孔喷嘴、多孔喷嘴。本文以常见单孔圆柱形喷嘴为研究对象，研究喷嘴型面及关键参数对高压水射流的影响。1.1 喷嘴形式圆柱形喷嘴的内部型面为旋成曲面，其母线的结构示意图如图1所示，喷嘴的结构参数主要有：入口直径和出口直径，收缩段半锥角，接入段长度、收缩段长度和平直延伸段长度。图1 圆柱形喷嘴结构示意图Fig.1 Schematic diagram of configuration of cylindric nozzle基于上述喷嘴结构，本文提出了另外两种内部型面的喷嘴，它们的、和与图1中的完全一样，不同之处在于其收缩段母线由直线变成了曲线。如图2所示为三种不同母线形式的圆柱形喷嘴示意图。图中的直线段代表的就是图1中的常见圆柱形喷嘴，凹型和凸型曲线代表的就是本文提出的两种喷嘴，图2中凹型和凸型曲线在两端分别与水平和竖直线相切。为表述方便，往下将这三种喷嘴分别简称为直线型喷嘴、凹型喷嘴和凸型喷嘴。图2 三种不同母线形式的圆柱形喷嘴示意图Fig.2 Schematic diagram of cylindric nozzle with three different kind of generatrix1.2 高压水射流流动特性分析按照流动状态可将射流分为层流射流和紊流射流。本文研究的高压水射流流速较大，流动应处于紊流状态，因而属于自由紊动射流15。水射流以初始速度射入周围静止的空气（如图3），由于存在较大的速度差，两种流体之间存在间断面，并且该间断面是不稳定的，一旦受到扰动就会产生波动，形成旋涡。这些旋涡会把周围原来静止的空气卷吸到射流中，随着紊动的发展，被卷吸并和射流一起运动的空气越来越多，射流边界也逐渐向两侧扩展。同时由于空气与射流的掺混，必然会产生对射流的阻力，使得射流边缘部分的流速降低，难以保持初始速度。水射流与空气的掺混自边缘逐渐向中心发展，在离开喷嘴出口一段距离之后，掺混现象发展到射流中心，从此往后，射流整个横断面上都发展成为紊流16。从上面的分析可知，高压水射流的流动结构十分复杂，为便于分析，可以按照不同的流动特性将水射流划分为几个区域（如图3）。从喷嘴出口开始向内外扩展的掺混区称为剪切层，其中心未受掺混影响，仍然保持初始速度的区域称为势流核。沿射流方向从喷嘴出口至势流核末端为初始段，初始段下游绝大部分为充分发展的紊动掺混区，称为主体段，在初始段和主体段之间存在很短的过渡段。图3 高压水射流流动特性示意图Fig.3 Schematic diagram of flow character of high pressure water jet2 研究方法2.1 流场分析方法采用Fluent 6.3.26软件进行水射流的流场分析。水射流与空气之间会发生剧烈的动量交换和紊动掺混，为了模拟气液两相混合运动，采用Mixture两相流混合模型。在考虑射流的湍流运动时，采用湍流模型。由于射流流场具有很好的轴对称特性，为了节省计算资源，采用如图4所示二维结构网格进行流场区域的剖分，其中(a)为全流场域网格，(b)为喷嘴附近网格。(a) 全流场域网格(a) Computing grid of the whole flow field(b) 喷嘴附近网格(b) Computing grid around nozzle图4 计算网格Fig.4 Computing grid计算区域的边界设置如图5。AB为中心轴线，设置为轴对称边界条件；喷嘴进口AH设置为压力入口条件；外边界IM、MK和KB设置为压力出口条件；HGFEJI设置为无滑移壁面边界条件。图5 计算边界示意图Fig.5 Schematic diagram of computing boundary采用本文的计算模型对文献12中半锥角为30度的圆柱形喷嘴的内外流场进行计算流体力学（Computational Fluid Dynamics, CFD）模拟，得到了射流沿轴线的速度分布。表1给出了其与文献结果的对比，可见结果吻合得很好，说明本文的流场仿真模型精度较好，计算结果可信。表1 圆柱形喷嘴射流轴线速度验证Table1 Validation of axis velocity of high pressure water jet with cylindric nozzle轴向位置(mm)405075125200文献结果(m/s)281.3277.2271.2268.5260.1本文结果(m/s)282.5279.3273.5270.2264.32.2 代理模型构建及验证在进行喷嘴几何参数的优化设计时，需要反复获取射流的流场特性，鉴于两相流场计算的耗时性，为了节省计算资源，引入Kriging代理模型来替代真实的CFD分析模型。本文构建代理模型的基本流程为：首先采用均匀采样的方式选取10k（k为设计变量的个数）个初始样本点，并采用CFD方法获取样本点上的流场特性，然后基于样本点上结果构建初始代理模型。一般来说此时代理模型的精度较差，并不能满足优化设计的要求。为了减少CFD分析的次数，提高设计效率，作者采用之前已经成功运用的改进EI加点方法17进行加点迭代，每一步迭代时，同时加入EI值最大的点和Kriging模型预测方差最大的点，不断修正代理模型的精度。采用逐一交叉验证方法18检验每次迭代后的代理模型的精度。经过8次左右的加点迭代，标准交叉验证残差全部在之间，且代理模型的最大误差不超过10%，这样的模型精度完全可以替代CFD分析进行优化设计研究。2.3 优化算法改进的非劣分类遗传算法（Nondominated Sorting Genetic Algorithm, NSGA-II）由于适应性强、计算效率高，得到了广泛的应用，也经常被用来作为很多优化算法比较的基准。本文采用NSGA-II进行单目标优化，可以获得较为合理的优化解。在进行多目标优化设计时，采用计算效率更高的基于分解的多目标进化算法（MultiObjective Evolutionary Algorithms based on Decomposition, MOEA/D）进行优化设计。该算法的基本思想19是将多目标优化问题采用一定的方法分解成若干单目标子问题，并对它们同时进行求解。在优化每一个子问题时只需要用到与其相邻的子问题的信息，这样MOEA/D每代的计算复杂度都比经典进化算法（如NSGA-II）的计算复杂度小，计算效率也更高19,20。3 喷嘴形式对射流特性的影响本节采用CFD方法研究图2所示的三种喷嘴形式对射流特性的影响，为下节锥角的优化设计奠定基础。假设喷嘴的基本几何参数相同：D=5mm、d=1mm、l3=3mm。以三种半锥角条件下的喷嘴为研究对象，具体的几何参数见表2。喷嘴入口的压力均为40MPa，外界环境为一个大气压。表2 不同半锥角条件下的喷嘴几何参数Table2 Configuration parameters of nozzle with different semi-cone angle(mm)(mm)Shape17186.34Shape2151011.31Shape323245图6给出了不同半锥角条件下，三种喷嘴形式产生的射流中心轴线速度及距喷嘴出口0.3m截面上（即图3中x=0.3m）速度的径向分布。可以看出射流中心轴线速度在喷嘴内靠近出口附近急剧增大，喷嘴出口处的最大速度与理论计算速度12相符。由于粘性的作用，随着射流与周围空气的相互掺混，轴线速度逐渐减小，且在离开喷嘴一定距离后存在一个小的突降，这个距离即为初始段长度。三种情况下轴向距离为0.3m处射流速度的径向分布较为类似，各种情况下最大的区别在于截面上中心轴线处速度的不同。Shape1和Shape2的射流特性较为相似，它们的半锥角较小，在小锥角情况下，直线型喷嘴的射流初始段最长、中心轴线速度和0.3m处射流速度最大，因而其产生的射流性能最优；凸型喷嘴射流初始段最短、中心轴线速度和0.3m处射流速度最小，射流性能最差；凹型喷嘴产生的射流性能介于二者之间，且与直线型喷嘴较为接近。当半锥角为45度时，凸型喷嘴射流的性能仍然最差，其中心轴线速度较小锥角时甚至更小；凹型喷嘴的产生的射流性能最佳，其中心轴线速度较小锥角时直线型喷嘴产生的射流中心速度甚至更大；直线型喷嘴的射流性能略差于凹型喷嘴的射流性能。从图中还可以发现，随着锥角的增大，喷嘴内部流动损失增大，直线型和凸型喷嘴的射流性能略微变差，但在大锥角情况下，凹型喷嘴使得射流向中心轴线集中的效果更加明显，因而其产生的射流性能反而变好。总之在小锥角条件下适宜采用直线型喷嘴，大锥角条件下适宜采用凹型喷嘴。直线型喷嘴总体性能较优且易于加工，凹型喷嘴次之，凸型喷嘴性能最差。(a) 外形1射流中心轴线速度(a) Central axis velocity of shape1(b) x=0.3m截面上外形1射流速度沿径向分布(b) Velocity distribution along radius direction of shape1 when x=0.3m(c) 外形2射流中心轴线速度(c) Central axis velocity of shape2(d) x=0.3m截面上外形2射流速度沿径向分布(d) Velocity distribution along radius direction of shape2 when x=0.3m(e) 外形3射流中心轴线速度(e) Central axis velocity of shape3(f) x=0.3m截面上外形3射流速度沿径向分布(f) Velocity distribution along radius direction of shape3 when x=0.3m图6 三种喷嘴形式下射流中心轴线速度及x=0.3m处速度的径向分布Fig.6 Central axis velocity and velocity distribution along radius direction when x=0.3m of high pressure water jet with three different nozzle4 圆柱形喷嘴锥角优化设计4.1 优化问题以上节分析的直线型喷嘴为研究对象。半锥角为唯一的设计变量，其变化范围为。半锥角取最小值时，对应于图1中时的情况；当半锥角取最大值时，对应于的情况。初始段长度L对高压水射流清洗和切割作业都十分关键，另外在进行固体发动机推进剂的清洗和切割时，要求射流流量Q尽可能小，以便于废液的处理，因而L和Q作为本文设计的目标。由于初始段长度定义为喷嘴出口至势流核末端距离，因而L可以从射流中心轴向流速分布中获取，实际操作时，在轴线上射流流速与出口中心流速之差小于1m/s的区域都为初始段。另外为了便于比较，取表2中的Shape1（直线型喷嘴）为基准外形。4.2 单目标优化分别以L和Q为目标进行单目标优化设计，得到的设计结果如表3所示。可见当半锥角为15.44度时，初始段长度最大，较基准外形时增大14.71%，但此时射流流量也略微增大；当半锥角为89.88度时，射流流量最小，较基准外形时减小27.56%，但此时初始段长度却剧减76.22%。可见两个优化目标存在一定的冲突，很难找到合适的设计变量使得两个目标同时达到最优，因而需要进行多目标优化。表3 单目标优化设计结果Table3 Single objective optimization results以L最大的设计结果以Q最小的设计结果优化值与基准外形差别优化值与基准外形差别L(mm)35.3814.71%7.336-76.22%Q(kg/s)0.2140.66%0.154-27.56%15.4489.884.3 多目标优化为了兼顾射流的整体性能，采用MOEA/D算法进行两目标优化设计，得到的设计结果如图7所示。图中黑色圆点为优化前沿，是优化结果在目标空间的映射。前沿上两个端点对应的解即为单目标优化结果。设计者可以根据具体的设计需求选择合适的优化解。另外从设计结果来看，前沿对应的设计结果（半锥角）处于区间上，且半锥角越小对应于初始段长度L越大、射流流量Q也越大，半锥角越大则L越小、但Q越小。图7 多目标优化结果Fig.7 Multi-objective optimization results5 结论本文首先分析了高压水射流的流场结构特性，为喷嘴型面的影响分析和设计奠定基础。然后通过对典型半锥角情况下直线型、凹型和凸型喷嘴内外流场的一体化求解，分析三种构型喷嘴对射流性能的影响。最后以直线型喷嘴为对象，以半锥角为目标进行了单目标和多目标优化设计。通过研究得到以下结论：(1) 小锥角情况下，直线型喷嘴产生的射流性能最佳，大锥角情况下，凹型喷嘴最合适。直线型喷嘴总体性能较优且易于加工，凹型喷嘴次之，凸型喷嘴性能最差。(2) 以直线型喷嘴为设计对象，得到了以射流初始段长度和流量为目标的单目标优化设计结果，两个单项目标分别提高14.71%和27.56%。得到了多目标优化设计前沿，优化得到的半锥角处于区间上。(3) 本文用于流场分析的两相流CFD模型具有较好的精度，可以满足优化设计的需求。基于代理模型和进化算法的全局优化方法在进行喷嘴的优化设计时是有效的。参考文献1陈亮. 高压水射流扇形喷嘴内外流场仿真分析 D. 兰州: 兰州理工大学, 2007.CHEN Liang. The simulation of high-pressure water jet nozzles of internal and external flow field D. Lanzhou: Lanzhou University of Technology, 2007.2杨林, 唐川林, 张凤华. 高压水射流技术的发展及应用 J. 洗净技术, 2004, 2(1): 9-14.YANG Lin, TANG Chuanlin, ZHANG Fenghua. Development and application of high pressure water jet J. Cleaning Technology, 2004, 2(1): 9-14. (in Chinese)3承文, 柳玉杰. 水射流技术在火箭发动机研制生产中的应用 J. 化学清洗, 1999, 15(5): 39-43.CHENG Wen, LIU Yujie. Application of water jet technology on research and production of rocket engine J. Chemical Cleaning, 1999, 15(5): 39-43. (in Chinese)4刘庭成, 范晓红, 刘焱. 高压水射流清洗机喷嘴的结构与参数 J. 清洗世界, 2010, 26(9): 32-37.LIU Tingcheng, FAN Xiaohong, LIU Yan. Structure and parameters of high pressure water-jet cleaner nozzles J. Cleaning World, 2010, 26(9): 32-37. (in Chinese)5刘庭成. 高压水射流流动过程中能量损失的分析 J. 清洗世界, 2004, 20(5): 9-11.LIU Tingcheng. Analysis on the energy loss in the flowing process of high pressure water jet J. Cleaning World, 2004, 20(5): 9-11. (in Chinese)6王洪伦, 龚烈航, 姚笛. 高压水切割喷嘴的研究 J. 机床与液压, 2005, 4: 42-43.WANG Honglun, GONG Liehang, YAO Di. Investigation of cutting nozzle in high pressure water jet J. Machine Tools and Hydraulic, 2005, 4: 42-43. (in Chinese)7蒋大勇. 高压水射流切割HTPB推进剂的实验分析 J. 解放军理工大学学报, 2013, 14(4): 430-435.JIANG Dayong. Chipping performance of HTPB propellant by high-pressure water jet J. Journal of PLA University of Science and Technology, 2013, 14(4): 430-435. (in Chinese)8喻峰, 祝爽, 李璐, 等. 新型高压水除磷喷嘴流场仿真与试验研究 J. 机械设计与制造, 2014, 1: 91-96.YU Feng, ZHU Shuang, LI Lu, et al. The flow-field simulation and experimental research of the new-type high-pressure water jet nozzle in descaling J. Machinery Design and Manufacture, 2014, 1: 91-96. (in Chinese)9易松林, 汪志明, 孙怡红, 等. 径向钻井高压水射流喷嘴内外流场分析 J. 石油机械, 2013, 41(3): 15-20.YI Songlin, WANG Mingzhi, SUN Yihong, et al. Analysis of the internal and external flow fields of high pressure water jet nozzle in radial drilling J. China Petroleum Machinary, 2013, 41(3): 15-20. (in Chinese)10王东, 吴雨川, 罗维平, 等. 高压水射流切割流场的数值模拟研究 J. 武汉科技学院学报, 2005, 18(3): 15-18.WANG Dong, WU Yuchuan, LUO Weiping, et al. Study of flow field on cutting with high pressure water jet J. Journal OF Wuhan University of Science and Engineering, 2005, 18(3): 15-18. (in Chinese)11于洪, 陆庭侃. 高压水射流切割喷嘴的结构设计和参数优化数值模拟研究 J. 机床与液压, 2009, 37(11): 90-92.YU Hong, LU Tingkan. Numerical simulation study on structure design and parameters optimization of high pressure water jet cutting nozzle J. Machine Tools and Hydraulic, 2009, 37(11): 90-92. (in Chinese)12周文会. 高压水射流喷嘴内外部流场的数值模拟研究 D. 兰州: 兰州理工大学, 2008.ZHOU Wenhui. Numerical simulation research on the high pressure water jet nozzle in inside and outside efflux flow field D. Lanzhou: Lanzhou University of Technology, 2008.13邵忠良. 纯高压水射流除锈喷嘴的数值模拟及优化试验研究 J. 机床与液压, 2015, 43(13): 78-81.SHAO Zhongliang. Numerical simulation and optimization experiment research of high-pressure water jet descaling nozzle J. Machine Tools and Hydraulic, 2015, 43(13): 78-81. (in Chinese)14朱左明, 高鑫, 王煊军, 等. 固体火箭发动机水射流清理工艺的多目标优化 J. 含能材料, 2014, 22(5): 624-629.ZHU Zuoming, GAO Xin, WANG Xuanjun, et al. Multi-objective optimization for waterjet cleaning process of solid rocket motor J. Chinese Journal of Energetic Materials, 2014, 22(5): 624-629. (in Chinese)15刘沛清. 自由紊动射流理论 M. 北京: 北京航空航天大学出版社, 2008.LIU Peiqing. Free turbulent jet theory M. Beijing: Beijing University of Aeronautics and Astronautics Press, 2008. (in Chinese)16董志勇.射流力学 M. 北京: 科学出版社, 2005.DONG Zhiyong. Jet dynamics M. Beijing: Science Press, 2005. (in Chinese)17马洋, 杨涛, 张青斌. 高超声速滑翔式升力体外形设计与优化 J. 国防科技大学学报, 2014, 36(2): 34-40.MA Yang, YANG Tao, ZHANG Qingbin. Configuration optimization design of hypersonic gliding lifting body J. Journal of National University of Defense Technology, 2014, 36(2): 34-40. (in Chinese)18Donald R J, Matthias S, William J W. Efficient global optimization of expensive black-box functions J. Journal of Global Optimization, 1998, 13: 455492.19ZHANG Q, LI H. MOEA/D: a multiobjective evolutionary algorithm based on decomposition J. IEEE Transactions on Evolutionary Computation, 2007, 11(6): 712-31.20马洋. 基于代理模型和MOEA/D的飞行器气动外形优化设计研究 D. 长沙: 国防科学技术大学, 2015.MA Yang. Research on the Aerodynamic Shape Optimization of Vehicle Based on Surrogate Model and MOEA/D D. Changsha: National University of Defence Technology, 2015.收稿日期：2015-11-20基金项目：国家自然科学基金青年科学基金项目（E031303）第二炮兵工程大学科研基金青年项目（2015QNJJ034）作者简介：韩启龙（1979），男，甘肃宁县人，硕士，副教授，E-mail：longfeng.061106马洋（通讯作者），男，湖南澧县人，博士，讲师，E-mail：mldy0612