1.一种离心泵叶轮的优化设计方法，其特征在于：按照如下步骤进行：

1)首先根据设计要求，依据设计经验来具体确定一组需要进行优化的参数：叶轮进口宽度b1、叶轮出口宽度b2、叶片进口角β1、叶片出口角β2、叶片数z，这一组需要进行优化的参数并在下面给出了约束范围，其它叶轮参数则作为设计常量；

所述需要进行优化的参数的约束范围为：

30°＜β1＜40°

15°＜β2＜40°

上式中

Km1-叶轮进口轴面系数，

ns-泵的比转速，

Fr-流量，

h-扬程，

n-叶轮转速，

g-重力加速度；

2)建模和网格划分：按照离心泵叶轮的参数，使用三维造型软件在确定的坐标系下对离心泵叶轮的计算域进行造型；利用ANSYS ICEM软件对三维模型进行网格划分；

3)CFX仿真验证：将网格导入CFX软件中，在给定的边界条件下进行离心泵叶轮在全流量工况下的瞬态数值模拟计算，由软件自动生成结果；再对数值模拟结果进行后处理以获得离心泵的流量-扬程、流量-效率性能曲线；根据性能曲线对比在所绘制的网格下扬程和效率是否达到设计要求，若不满足设计要求，则重复步骤2)对三维模型重新进行网格划分，直至得到满足设计要求的网格；

4)涡核结构提取：以步骤3)得到的模拟结果为基础，基于Q方法在计算域中进行涡核结构提取，得到涡核强度随流量变化的数据，其中，Q值表示流体微元旋转与变形的强度大小关系；

5)熵产分析：以步骤3)得到的网格为基础，根据公式(1.1)、(1.2)、(1.3)，由CFD计算结果的速度分布计算离心泵内部能量损失的分布情况；然后由公式(1.4)对式(1.2)、(1.3)在叶轮区域进行体积积分可得叶轮总熵产率；

上式中

Spro，t-单位体积叶轮总熵产生率，

-单位体积叶轮时均熵产生率，

Spro，d′-单位体积叶轮脉动熵产生率，

T-流体温度，

ρ-流体密度，

ε-湍流耗散率，

u，v，w-流体分别在X、Y、Z方向上的速度分量， 表示时均速度值，u′，v′，w′表示脉动速度值；

6)高效点匹配：将步骤4)、步骤5)中得到的数据进行无量纲化处理，导入到数据处理软件中进行绘图，得到叶轮总熵产率和涡核强度随流量变化的曲线；其中涡核强度最小时对应的无量纲流量为Fr1，叶轮总熵产率最小时对应的无量纲流量为Fr2；

7)采用模拟退火优化算法，对叶轮的设计参数进行优化，使得涡核强度最小时对应的无量纲流量和叶轮总熵产率最小时对应的无量纲流量均为F′r，且F′r＝Fr1＝Fr2＝1，此时得到叶轮的最优设计参数。

2.根据权利要求1所述的一种离心泵叶轮的优化设计方法，其特征在于：所述步骤4)中，Q方法提取涡结构的步骤为：根据公式(1.5)、(1.6)、(1.7)确定Q值：变形率张量-

涡量-

Q值-

上式简化到笛卡尔坐标下Q值如下：

其中，Q方法定义Q＞0的区域为涡管所在区域，上式中变形率张量和涡量则是由数值模拟结果在后处理软件CFD-POST中通过编辑CEL表达式计算得到；提取Fr＝Frd的微元等值面确定旋涡形态，其中，Frd为设计流量；确定Q值之后，在CFD-POST中，作出该Q值下的涡量等值面，采用旋涡强度系数进行着色，再由软件定量计算出涡核强度，并记录数据。

3.根据权利要求1所述的一种离心泵叶轮的优化设计方法，其特征在于：所述步骤5)中，对叶轮总熵产率进行无量纲化处理，是采用每一个叶轮熵产率的数值除以叶轮熵产率的最小值；对涡核强度进行无量纲化处理，是采用每一个涡核强度的数值除以涡核强度的最小值；对流量进行无量纲化处理，是采用每一个流量的数值除以设计流量Frd。

4.根据权利要求1所述的离心泵叶轮的优化设计方法，其特征在于：所述步骤3)中，边界条件包括进口压力、出口流量、壁面和动静交界面设置以及叶轮转速。