《一种基于CFD的离心泵多工况水力优化方法》是江苏大学于2010年10月27日申请的发明专利，该专利的申请号为2010105205616，公布号为CN101956711A，公布日为2011年1月26日，发明人是刘厚林、王凯、袁寿其、谈明高、王勇，该专利属于离心泵水力优化领域。
《一种基于CFD的离心泵多工况水力优化方法》包含以下步骤：第一，基于外特性实验建立离心泵多工况水力优化的初始样本库；第二，建立基于CFD的离心泵多工况水力优化模型；第三，基于单点水力设计和均匀设计法，采用CFD技术生成离心泵多工况水力优化的样本库，并构建响应面近似模型；第四，采用全局优化算法对近似模型进行全局优化；第五，根据优化结果，采用PCAD对离心泵进行水力设计，从而得到离心泵的多工况水力模型。
2014年11月6日，《一种基于CFD的离心泵多工况水力优化方法》获得第十六届中国专利优秀奖。
（概述图为《一种基于CFD的离心泵多工况水力优化方法》摘要附图 ）
截至2010年10月，传统的离心泵设计方法是单工况（设计工况）设计方法，无法兼顾离心泵多个工况之间的性能，其多个工况下的性能主要是通过反复的“设计—计算（或试验）—修改”实现的。在设计开发过程中没有充分利用CFD技术进行离心泵的辅助分析和优化，经验仍占主导地位，CFD仅是用来检验设计的结果，而不是用来驱动产品设计。因此，该方法不但设计效率低，而且高效区范围窄，主要原因是其主要几何参数不是最优解。
若能将离心泵传统的“设计—计算（或试验）—修改”过程通过一集成平台自动探索最优方案，则离心泵将会突破传统设计的瓶颈，实现离心泵多工况水力优化设计。随着计算机技术和优化算法的快速发展，上述假设已变成现实。iSIGHT优化设计平台可以实现各种CFD软件和自编程序的集成，并提供先进的优化设计方法，实现了设计过程的自动化和数据的可视化。因此，利用iSIGHT平台集成CFD软件和自编程序，实现离心泵多工况水力优化设计，最终达到离心泵“设计—计算—修改”过程的自动化，具有重要的现实意义。
经检索，截至2010年10月尚未见关于基于CFD的离心泵多工况水力优化方法的文献和申报专利，仅有一些学者做离心泵非设计工况下CFD数值模拟等研究工作。
《一种基于CFD的离心泵多工况水力优化方法》的目的是提供一种基于CFD的离心泵多工况水力优化方法，通过采用外特性实验、CFD技术和全局优化算法来解离心泵多工况水力优化问题，从而求解出一组关键几何参数的最优解集。
采用外特性实验建立初始样本库，并采用均匀设计方法和CFD技术补充样本库，从而构建起响应面近似模型；并在此基础上，采用iSIGHT软件集成Pro/E、Gambit、Fluent以及自编的离心泵多工况性能计算程序pmcc.exe，选用自适应模拟退火算法（或基于初始种群渐进漂移的自适应遗传算法）来求解离心泵多工况水力优化问题。
其水力设计方法步骤如下：
（1）基于外特性实验建立离心泵多工况水力优化的初始样本库；
搭建离心泵外特性试验台；离心泵扬程H由离心泵进、出口的压力表测量得到；采用电测法测量离心泵的功率P；离心泵流量Q由离心泵出口管路系统上的电磁流量计读出；离心泵的效率η由公式η=ρgQH/P；
通过测量得到的不同比转数、不同工况下离心泵进、出口之间的压力差，建立关键几何参数（D2、β2、b2、z、D1、β1、b1、Dj、φ、Ft、D3、φ0、b3、Ds、Dd）与不同工况下离心泵扬程Hi之间的关系；通过测量得到的不同比转数、不同工况下离心泵功率Pi，建立关键几何参数与功率Pi之间的关系；通过测量出口管路的流量Qi，建立关键几何参数与不同工况下离心泵效率ηi之间的关系；
基于上述外特性实验，建立离心泵多工况水力优化的初始样本库。
（2）建立基于CFD的离心泵多工况水力优化模型；
求设计变量x=，使
，
。式中：i为各工况点，i≥3；Δz为泵进、出口总水头之差；M为叶轮扭矩；w为叶轮的角速度；D2为叶轮出口直径；β2为叶片出口安放角；b2为叶轮出口宽度；z为叶片数；D1为叶片进口直径；β1为叶片进口安放角；b1为叶片进口宽度；Dj为叶轮进口直径；φ为叶片包角；Ft为蜗壳喉部面积；D3为蜗壳基圆直径；φ0为蜗壳隔舌角；b3为蜗壳进口宽度；Ds为泵进口直径；Dd为泵出口直径。
（3）基于单点水力设计和均匀设计法，采用CFD技术生成离心泵多工况水力优化的样本库，并构建响应面近似模型；
基于单点设计方法对离心泵进行初始设计；选择均匀设计方法作为实验设计方法，以减少CFD试验次数；选择关键几何参数D2、β2、b2、z、D1、β1、b1、Dj、φ、Ft、D3、φ0、b3、Ds、Dd作为设计参数，根据均匀设计方法确定了50组离心泵CFD试验方案；
根据上述50组CFD实验方案，分别采用Pro/E对其进行三维造型，并输出关键几何参数的关系表达式input.txt作为下一步优化的输入文件；分别采用Gambit对其进行网格化分；分别采用Fluent数值计算，得到其样本点的性能，并保存到样本库中；构建样本点的几何参数与离心泵性能之间的近似模型。
（4）采用全局优化算法对近似模型进行多工况多目标全局优化；
采用VisualC++2005将其编写成离心泵多工况性能计算程序pmcc.exe，该程序根据不同工况下离心泵CFD数值计算得到的扬程H和扭矩N,计算出其功率P以及效率η，并将扬程H、功率P以及效率η保存到output.txt中；
采用iSIGHT集成Pro/E、Gambit、Fluent以及自编的离心泵多工况性能计算程序，并以离心泵单点设计得到的关键几何参数值为设计变量、多个工况下的扬程和功率为约束条件、多个工况下效率最大为目标函数，同时使用自适应遗传算法、或基于初始种群渐进漂移的自适应遗传算法，对近似模型进行多目标全局优化；
如果不满足收敛准则，则把CFD计算得到的性能值保存到样本库，重新构建近似模型，并改变设计变量值进行全局优化，直至满足收敛准则为止；收敛后，将关键几何参数的最优解及离心泵的性能值分别保存在input.txt和output.txt中。
（5）根据优化结果，采用泵水力设计软件PCAD2010对离心泵进行水力设计，从而得到离心泵多工况水力优化模型。
1、能够根据离心泵多个工况下的扬程和功率要求，以并以各工况点的效率最大为目标函数，对其进行多工况水力优化设计；
2、还可以根据多个工况点下的扬程要求，并以各工况点的功率最小、效率最大为目标函数，对离心泵进行节能改造；
3、还可以为主动控制小流量工况下不稳定流动和抑制大流量工况下汽蚀发生提供借鉴。
图1为离心泵多工况水力优化设计方法的流程图。
图2为《一种基于CFD的离心泵多工况水力优化方法》实施例的优化系统图。
|
|
1.《一种基于CFD的离心泵多工况水力优化方法》其特征在于，根据给定多个工况点下的扬程和功率值，基于外特性实验和CFD数值计算得到的样本库构建的近似模型，采用合适的全局优化算法对离心泵多工况水力优化问题进行求解，从而计算出离心泵关键几何参数的最优值，其中，工况点i≥3；具体步骤如下：
（A）基于外特性实验建立离心泵多工况水力优化的初始样本库；搭建离心泵外特性试验台；离心泵扬程H由离心泵进、出口的压力表测量得到；采用电测法测量离心泵的功率P；离心泵流量Q由离心泵出口管路系统上的电磁流量计读出；离心泵的效率η由公式η=ρgQH/P；通过测量得到的不同比转数、不同工况下离心泵进、出口之间的压力差，建立关键几何参数（D2、β2、b2、z、D1、β1、b1、Dj、φ、Ft、D3、φ0、b3、Ds、Dd）与不同工况下离心泵扬程Hi之间的关系；通过测量得到的不同比转数、不同工况下离心泵功率Pi，建立关键几何参数与功率Pi之间的关系；通过测量出口管路的流量Qi，建立关键几何参数与不同工况下离心泵效率ηi之间的关系；基于上述外特性实验，建立离心泵多工况水力优化的初始样本库；
（B）建立基于CFD的离心泵多工况水力优化模型；求设计变量x=，使
，。式中：i为各工况点，i≥3；Δz为泵进、出口总水头之差；M为叶轮扭矩；w为叶轮的角速度；D2为叶轮出口直径；β2为叶片出口安放角；b2为叶轮出口宽度；z为叶片数；D1为叶片进口直径；β1为叶片进口安放角；b1为叶片进口宽度；Dj为叶轮进口直径；φ为叶片包角；Ft为蜗壳喉部面积；D3为蜗壳基圆直径；φ0为蜗壳隔舌角；b3为蜗壳进口宽度；Ds为泵进口直径；Dd为泵出口直径；
（C）基于单点水力设计和均匀设计法，采用CFD技术生成离心泵多工况水力优化的样本库，并构建响应面近似模型；基于单点设计方法对离心泵进行初始设计；选择均匀设计方法作为实验设计方法，以减少CFD试验次数；选择关键几何参数D2、β2、b2、z、D1、β1、b1、Dj、φ、Ft、D3、φ0、b3、Ds、Dd作为设计参数，根据均匀设计方法确定了50组离心泵CFD试验方案；根据上述50组CFD实验方案，分别采用Pro/E对其进行三维造型，并输出关键几何参数的关系表达式input.txt作为下一步优化的输入文件；分别采用Gambit对其进行网格化分；分别采用Fluent数值计算，得到其样本点的性能，并保存到样本库中；构建样本点的几何参数与离心泵性能之间的近似模型；
（D）采用全局优化算法对近似模型进行多工况多目标全局优化；采用VisualC++2005将其编写成离心泵多工况性能计算程序pmcc.exe，该程序根据不同工况下离心泵CFD数值计算得到的扬程H和扭矩N,计算出其功率P以及效率η，并将扬程H、功率P以及效率η保存到output.txt中；采用iSIGHT集成Pro/E、Gambit、Fluent以及自编的离心泵多工况性能计算程序，并以离心泵单点设计得到的关键几何参数值为设计变量、多个工况下的扬程和功率为约束条件、多个工况下效率最大为目标函数，同时使用自适应遗传算法、或基于初始种群渐进漂移的自适应遗传算法，对近似模型进行多目标全局优化；如果不满足收敛准则，则把CFD计算得到的性能值保存到样本库，重新回到步骤（C）构建近似模型，并改变设计变量值进行全局优化，直至满足收敛准则为止；收敛后，将关键几何参数的最优解及离心泵的性能值分别保存在input.txt和output.txt中；
（E）根据优化结果，采用泵水力设计软件对离心泵进行水力设计，从而得到离心泵多工况水力优化模型。
2.根据权利要求1所述的一种基于CFD的离心泵多工况水力优化方法，其特征在于：所述步骤（B），离心泵多工况水力优化模型是在离心泵单点设计的基础上，使用CFD技术对离心泵进行多工况全流场数值计算得到的离心泵多工况水力优化模型。
该实施例基于CFD的离心泵多工况水力优化方法的流程如图1所示，包括如下步骤：
（1）基于外特性实验建立离心泵多工况水力优化的初始样本库。
搭建离心泵外特性试验台；离心泵扬程H由离心泵进、出口的压力表测量得到；采用电测法测量离心泵的功率P；离心泵流量Q由离心泵出口管路系统上的电磁流量计读出；离心泵的效率η由公式η=ρgQH/P。
通过测量得到的不同比转数、不同工况下离心泵进、出口之间的压力差，建立关键几何参数（D2、β2、b2、z、D1、β1、b1、Dj、φ、Ft、D3、φ0、b3、Ds、Dd）与不同工况下离心泵扬程Hi之间的关系；通过测量得到的不同比转数、不同工况下离心泵功率Pi，建立关键几何参数与功率Pi之间的关系；通过测量出口管路的流量Qi，建立关键几何参数与不同工况下离心泵效率ηi之间的关系。
基于上述外特性实验，建立离心泵多工况水力优化的初始样本库。
（2）建立基于CFD的离心泵多工况水力优化模型。
求设计变量x=，使
，。式中：i为各工况点，i≥3；Δz为泵进、出口总水头之差；M为叶轮扭矩；w为叶轮的角速度；D2为叶轮出口直径；β2为叶片出口安放角；b2为叶轮出口宽度；z为叶片数；D1为叶片进口直径；β1为叶片进口安放角；b1为叶片进口宽度；Dj为叶轮进口直径；φ为叶片包角；Ft为蜗壳喉部面积；D3为蜗壳基圆直径；φ0为蜗壳隔舌角；b3为蜗壳进口宽度；Ds为泵进口直径；Dd为泵出口直径。
（3）基于单点水力设计和均匀设计法，采用CFD技术生成离心泵多工况水力优化的样本库，并构建响应面近似模型。
基于单点设计方法对离心泵进行初始设计；选择均匀设计方法作为实验设计方法，以减少CFD试验次数；选择关键几何参数D2、β2、b2、z、D1、β1、b1、Dj、φ、Ft、D3、φ0、b3、Ds、Dd作为设计参数，根据均匀设计方法确定了50组离心泵CFD试验方案。
根据上述50组CFD实验方案，分别采用Pro/E对其进行三维造型，并输出关键几何参数的关系表达式input.txt作为下一步优化的输入文件；分别采用Gambit对其进行网格化分；分别采用Fluent数值计算，得到其样本点的性能，并保存到样本库中；构建样本点的几何参数与离心泵性能之间的近似模型。
（4）采用VisualC++2005将其编写成离心泵多工况性能计算程序pmcc.exe。
该程序能够根据不同工况下离心泵CFD数值计算得到的扬程H和扭矩N,计算出其功率P以及效率η，并将扬程H、功率P以及效率η保存到output.txt中。
（5）采用全局优化算法对近似模型进行多工况多目标全局优化，如图2所示。
采用iSIGHT集成Pro/E、Gambit、Fluent以及自编的离心泵多工况性能计算程序，并以离心泵单点设计得到的关键几何参数值为设计变量、多个工况下的扬程和功率为约束条件、多个工况下效率最大为目标函数，同时使用自适应遗传算法、或基于初始种群渐进漂移的自适应遗传算法，对近似模型进行多目标全局优化；
①将设计变量写入Pro/E的参数化模型表达式的模板文件input.txt中，通过Pro/E的批处理文件调用trail.txt完成几何模型的生成，并输出pump.stp文件；
Pro/E的批处理文件：
deltrail.txt.*
"C:proeWildfire2binproe.exe"pro_wait-g:no_graphicsinput.txt
②根据Pro/E输出的pump.stp文件，通过Gambit的批处理文件读入命令流文件mesh.jou自动进行网格划分，同时输出为网格文件pump.msh；
Gambit的批处理文件:
C:Gambitntbinntx86gambit.exe-inpmesh.jou
③通过Fluent批处理文件完成读入各个工况下的命令流文件solve.jou自动设置计算模型和边界条件等的加载，并将各工况下的进、出口总压及叶轮扭矩输出至pump-i.txt文件中；
Fluent批处理文件:
"C:Fluentntbinntx86fluent.exe"3d-isolve.jou
④根据不同工况下离心泵CFD数值计算得到的扬程H和扭矩N，采用自编的离心泵多工况性能计算程序pmcc.exe，计算离心泵的功率P以及效率η，并将扬程H、功率P以及效率η保存到output.txt中；
⑤根据自适应模拟退火算法、或基于初始种群渐进漂移的自适应遗传算法，并以单点设计的参数值作为初始值、多个工况下的扬程和功率要求作为约束条件、多个工况下的效率值最大作为目标函数对其进行多目标全局优化，从而迭代出关键几何参数的最优解；
⑥收敛后，将关键几何参数的最优解及离心泵的性能值分别保存在input.txt和output.txt中。
（6）根据优化结果，采用具有自主知识产权的泵水力设计软件PCAD2010对离心泵进行水力设计，从而得到离心泵多工况水力优化模型。
2014年11月6日，《一种基于CFD的离心泵多工况水力优化方法》获得第十六届中国专利优秀奖。
温馨提示：
