发明名称 一种基于CFD的可调导叶对泵能量性能预测的方法
摘要 本发明涉及一种基于CFD(Computational fluid dynamics)的可调导叶对泵能量性能预测的方法,属于水利工程和机械工程技术领域。一种基于CFD的可调导叶对泵能量性能预测的方法,其特征是:具体包括以下步骤:步骤1、基于CFD数值计算构建不同可调导叶安放角时泵的能量性能数据库;步骤2、基于Visual Fortran程序绘制不同可调导叶时泵的综合特性曲线并输出CFD数值计算的可调导叶安放角时任意数量的泵能量性能数据;步骤3、构建可调导叶对泵能量性能预测的数学模型;通过本发明,提供了一种简便可靠的基于CFD的可调导叶对泵能量性能预测的方法。
申请公布号 CN104636542A 申请公布日期 2015.05.20
申请号 CN201410852277.7 申请日期 2014.12.31
申请人 扬州大学 发明人 杨帆;刘超;周济人;王芃也;邵正杰
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 扬州苏中专利事务所(普通合伙) 32222 代理人 许必元
主权项 一种基于CFD的可调导叶对泵能量性能预测的方法,其特征是,基于泵的CFD数值计算构建泵的能量性能数据库;采用泵综合特性绘制Visual Fortran程序对泵能量性能数据进行处理绘制泵的综合特性曲线,并输出CFD数值计算的可调导叶安放角时任意数量的泵性能数据;以泵的效率和泵的扬程为因变量,以可调导叶的安放角和泵的流量为自变量采用多元非线性回归分析方法分别构建通过2个自变量预测因变量的方法;具体包括以下步骤:1)基于CFD数值计算构建不同可调导叶安放角时泵的能量性能数据库;采用CFD数值计算程序,对配不同可调导叶安放角的泵进行8~10个工况点的能量性能的三维定常数值计算,基于CFD数值计算结果,分别求出泵的扬程、效率;泵的扬程H的计算式:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>H</mi><mo>=</mo><msub><mrow><mo>[</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>outlet</mi></msub></munderover><mfrac><msub><mrow><mo>(</mo><mfrac><msub><mi>P</mi><mi>soutlet</mi></msub><mi>&rho;g</mi></mfrac><mo>)</mo></mrow><mi>i</mi></msub><msub><mi>N</mi><mi>outlet</mi></msub></mfrac><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>outlet</mi></msub></munderover><mfrac><msub><mrow><mo>(</mo><mfrac><msup><msub><mi>v</mi><mi>outlet</mi></msub><mn>2</mn></msup><mrow><mn>2</mn><mi>g</mi></mrow></mfrac><mo>)</mo></mrow><mi>i</mi></msub><msub><mi>N</mi><mi>outlet</mi></msub></mfrac><mo>]</mo></mrow><mi>out</mi></msub><mo>-</mo><msub><mrow><mo>[</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>in</mi></msub></munderover><mfrac><msub><mrow><mo>(</mo><mfrac><msub><mi>P</mi><mi>sin</mi></msub><mi>&rho;g</mi></mfrac><mo>)</mo></mrow><mi>i</mi></msub><msub><mi>N</mi><mi>in</mi></msub></mfrac><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>in</mi></msub></munderover><mfrac><msub><mrow><mo>(</mo><mfrac><msubsup><mi>v</mi><mi>in</mi><mn>2</mn></msubsup><mrow><mn>2</mn><mi>g</mi></mrow></mfrac><mo>)</mo></mrow><mi>i</mi></msub><msub><mi>N</mi><mi>in</mi></msub></mfrac><mo>]</mo></mrow><mi>in</mi></msub><mo>+</mo><msub><mi>Z</mi><mi>out</mi></msub><mo>-</mo><msub><mi>Z</mi><mi>in</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000649751840000011.GIF" wi="1716" he="341" /></maths>式(1)中:P<sub>s</sub>为断面各节点的静压值;Z为断面几何中心的位能;v为断面各网格节点的绝对速度;N为断面网格节点数总和;下表outlet表示出口断面;下表in表示进口断面;泵的效率η的计算式:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>&eta;</mi><mo>=</mo><mfrac><mrow><mn>30</mn><mi>&rho;g</mi><msub><mi>QH</mi><mi>ZZ</mi></msub></mrow><mrow><mi>&pi;n</mi><msub><mi>T</mi><mi>p</mi></msub></mrow></mfrac><mo>;</mo><msub><mi>T</mi><mi>P</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mi>&Delta;</mi><msub><mi>A</mi><mi>i</mi></msub><mo>{</mo><mo>[</mo><mover><mi>r</mi><mo>&RightArrow;</mo></mover><mo>&times;</mo><mrow><mo>(</mo><mover><mi>&tau;</mi><mo>&RightArrow;</mo></mover><mo>&CenterDot;</mo><mover><mi>n</mi><mo>&RightArrow;</mo></mover><mo>)</mo></mrow><mo>]</mo><mo>&CenterDot;</mo><msub><mi>e</mi><mi>z</mi></msub><mo>}</mo><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mi>P&Delta;</mi><msub><mi>A</mi><mi>i</mi></msub><mo>[</mo><mrow><mo>(</mo><mover><mi>r</mi><mo>&RightArrow;</mo></mover><mo>&times;</mo><mover><mi>n</mi><mo>&RightArrow;</mo></mover><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mover><mi>e</mi><mo>&RightArrow;</mo></mover><mi>z</mi></msub><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000649751840000012.GIF" wi="1771" he="161" /></maths>式(2)中:<img file="FDA0000649751840000013.GIF" wi="57" he="73" />为转轴方向的单位向量;△A<sub>i</sub>为压力面或吸力面上第i单元的面积;<img file="FDA0000649751840000014.GIF" wi="48" he="62" />为△A<sub>i</sub>上的单位向量;<img file="FDA0000649751840000015.GIF" wi="40" he="63" />为向径;<img file="FDA0000649751840000016.GIF" wi="44" he="63" />为不含静压力P的应力张量;Q为泵装置流量;n为叶轮转速。2)基于Visual Fortran程序绘制不同可调导叶时泵的综合特性曲线,并输出CFD数值计算的可调导叶安放角时任意数量的泵能量性能数据;为绘制不同可调导叶安放角时泵的能量性能综合特性曲线,规定:可调导叶出口水流方向与转轮旋转方向相同时,可调导叶安放角为正,反之为负;采用Visual Fortran程序对步骤1建立的泵能量性能数据库进行自动处理绘制出泵的能量性能综合特性曲线,即泵的扬程、流量和效率3者的复杂关系曲线,并采用内插值的方法给出CFD数值计算的不同可调导叶安放角时任意数量的泵能量性能数据,为最终提出可调导叶对泵能量性能预测的数学模型奠定足够的泵能量性能样本数据;3)构建可调导叶对泵能量性能预测的数学模型;采用多元非线性回归分析方法(Multiple Non‑liner Regression)构造可调导叶对泵能量性能预测的数学模型,该数学模型一种快捷方法;采用Gauss‑Newton算法对非线性回归方程系数进行求解;将可调导叶安放角θ记为x<sub>1</sub>,流量Q记为x<sub>2</sub>,扬程H记为y<sub>1</sub>,可得非线性回归方程1;将进口导叶片安放角θ记为x<sub>1</sub>,流量H记为x<sub>2</sub>,泵效率η记为y<sub>2</sub>,可得非线性回归方程2,回归预测模型的准确度可通过数学模型的判定系数R<sub>2</sub>进行评判(判定系数R<sub>2</sub>需大于0.98),若判定系数R<sub>2</sub>不满足要求,则需回到步骤2,重新输出更多的泵性能数据样本;建立的2元3次非线性回归预测数学模型的通用表达式如下:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>y</mi><mo>=</mo><mi>a</mi><mo>+</mo><msub><mi>bx</mi><mn>1</mn></msub><mo>+</mo><msub><mi>cx</mi><mn>2</mn></msub><mo>+</mo><msubsup><mi>dx</mi><mn>1</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>ex</mi><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><msub><mi>fx</mi><mn>1</mn></msub><msub><mi>x</mi><mn>2</mn></msub><mo>+</mo><msubsup><mi>gx</mi><mn>1</mn><mn>3</mn></msubsup><mo>+</mo><msubsup><mi>hx</mi><mn>2</mn><mn>3</mn></msubsup><mo>+</mo><msub><mi>ix</mi><mn>1</mn></msub><msubsup><mi>x</mi><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>jx</mi><mn>1</mn><mn>2</mn></msubsup><msub><mi>x</mi><mn>2</mn></msub></mrow>]]></math><img file="FDA0000649751840000021.GIF" wi="1657" he="95" /></maths>式中:a<sub>1</sub>、b<sub>1</sub>、c<sub>1</sub>、d<sub>1</sub>、e<sub>1</sub>、f<sub>1</sub>、g<sub>1</sub>、h<sub>1</sub>、i<sub>1</sub>、j<sub>1</sub>为回归方程的常数及系数。
地址 225009 江苏省扬州市大学南路88号