首页 如何成为会员 意见反馈
主办 沈阳鼓风机研究所 /《风机技术》杂志社
   新闻  |   技术纵横  |  论坛  |  沈阳鼓风机研究所  |  风机协会  |  质检中心  |  风机标委会  |  风机技术杂志  |  企业商铺  |  供求  |  产品  |  书籍  |  招聘
当前位置:中国风机技术网 → 技术纵横 → 技术前沿

基于CFD的轴流通风机叶片的流场分析与改进设计

万 福/成都华川电装有限责任公司    

摘要 :通过计算流体力学( CFD )方法对轴流通风机叶片的流场进行了虚拟样机的数值模拟,不仅得到了流场的工作特性数据,而且提出了对叶片叶型的改进设计方案,并通过真实样机的试验验证了数值模拟分析的正确性和改进设计的可行性。最后,还对数值模拟与真实试验数据之间的差异原因进行了讨论。

关键词 :轴流式通风机;叶片; CFD ;流场分析;改进设计

中图分类号: TH432.1 文献标识码: B

文章编号 : 1006 - 8155 ( 2008 )02-0021-05

The Fl ow Field Analysis and Improved Design for Axial- fl ow Fan Leaf Based on CFD

Abstract : This paper numerically simulated the fl ow field for virtual protot yp e of the leaf of axial- fl ow fan by computational fl uid dynamics (CFD) method. The working characteristics data of fl ow field are not only obtained, but also the improved design plans for the leaf t yp es are presented and the correctness of the simulated analysis and the feasibility of the improved design are verified based on the real test for actual protot yp e. Finally, some possible reasons for the differences between simulated data and real test data are discussed.

Key words : axial- fl ow fan; leaf; CFD; fl ow field analysis; improved design

0  引言

  轴流通风机的传统设计方法主要有两种:一种是利用孤立翼型进行空气动力试验所得到的数据进行孤立翼型设计,称为孤立翼型设计方法;另一种是利用平面叶栅的理论和叶栅的吹风试验所得到的数据进行设计,称为叶栅设计方法 [1] 。试验测量方法所得到的试验结果真实可信,但往往受模型尺寸、流场扰动、人身安全和测量精度等的限制,有可能很难通过试验方法得到结果。此外试验还会遇到经费投入、人力和物力的巨大耗费及周期长等许多困难。计算流体力学( CFD )的计算方法是近年来发展起来的新型独立学科,它兼有理论性和实践性的双重特点,建立了许多理论和方法,为现代科学中许多复杂流动与传热问题提供了有效的计算技术 [2] 。轴流通风机叶片作为关键部件,其性能直接影响着风机的性能。轴流通风机设计的主要任务就是设计出能保证各项性能要求的高效率叶片。

  本文介绍的是采用现今先进的 CFD 方法,以我公司生产的一款汽车用冷凝器风扇的叶片为例,进行探索性的流场分析与改进设计研究。

1  叶片的 CFD 流场分析

1.1  对象描述

  该风扇总成的整体三维图如图 1 所示。叶片直径为 250mm ,材料为 PP ,其技术要求:在静压 p = - 50Pa (风机进口处的压力比周围空气低 50Pa ),转速为 2600r/min 的情况下,风扇总成在一个标准大气压、温度为 20 ℃、相对湿度为 50% 的空气,用标准电机在 12V 的电压下进行送风测试时,其送风量应≥ 900m 3 /h ,标准电机工作电流应≤ 7A 。


1.2  划分网格

  计算流体力学作为工程应用的有效工具,所面临的关键技术之一就是生成网格的质量的好坏,它直接影响到模拟结果的精度和所耗用的 CPU 时间。在计算敏感区域(壁面附近、尾流块、外形曲率大的表面)参数变化梯度大,如果网格太稀疏,则不能捕捉到流场的重要信息,造成误差大,甚至解不能收敛,故需取较密的一些网格;而在非计算敏感区域参数变化梯度较小,如果网格太稠密,则所耗用的 CPU 时间长,故应取较稀一些的网格。因此,应根据需要安排网格疏密。另外,曲线应尽量光滑,不能过分扭曲。在 CFD 的实际应用中,计算量一般非常大,合理划分网格可以大大节省机时,还可以避免自动网格划分中带来的过度疏密。高正交网格和高效展玄比是高质量网格的两个基本要素,它直接影响到计算的收敛时间和稳定性。

  显然该叶片是周期对称的,因而选用 1/5 的叶片进行计算。根据经验及反复调整,得到 36 万节点的叶片周围气体的网格如图 2 所示,其网格参数的评价结果(见表 1 )显示该生成的叶片网格是高效网格。为了避免因为网格密度不足而造成计算不可靠的影响,对改进前后的叶片的计算域均进行加大网格密度的验证,由原来 36 万节点的网格加密为 92 万节点的网格。对加密的网格进行计算得到的结果与原结果对比,误差不超过万分之一,说明原计算所用的网格密度足够,其计算结果是可靠的。

1.3  边界条件及求解参数设置

  在进行 CFD 流场分析时所采取的部分边界条件及参数设置如下:

  ( 1 )对流方式选用 k- w SST (剪切应力输运模型 [3] )模型;

  ( 2 )采用旋转坐标系的方法,设置计算域转速为技术要求的工作转速 2600r/min ;

  ( 3 )取参考压力为 1013 25Pa ,通过试算设定进口总压条件,调节进口总压的相对值使进口静压达到- 50Pa ,此时总压相对值约为- 29Pa ;

  ( 4 )由于在距叶片相对较近的一段距离内气体总压变化不大,且出口与大气相连,因而对风机出口设定出口平均静压相对值为 0Pa ;

  ( 5 )对于流线罩,考虑到实际状态为静止部件,其与叶片的距离较近,有相互作用,因而在边界条件中设为反向旋转的墙,即在旋转坐标系中流线罩是运动的。但由于流线罩是圆周对称的圆柱面,不会造成计算区域不稳定的问题;

  ( 6 )由于叶片顶部与流线罩相互作用强烈,生成叶片网格时将叶片与流线罩间单独分两个重合的面。利用两个重合的面建立一液体对液体的接触面边界条件。该边界条件充分利用了其可允许压力突变的特点,将叶片顶部与流线罩的复杂流动,通过两个重合面与叶片面分成流动相对平稳的两个区域。对于周期面自然选用周期边界条件,其余的面采用缺省的无滑移壁面边界条件。

  ( 7 )收敛控制与收敛准则:求解器的时间步长设为“自动调整时间步长”;最大迭代次数为 500 次;残差类型为均方根;前后计算的残差余量设为 0.00001 。

1.4  仿真结果与讨论

  为了得出风机的工作特性曲线,在只改变转速而不改变其余边界条件的情况下,得到如表 2 的仿真结果。

表 2 风扇总成的仿真特性表

转速/

(r/min)

进口静压/

Pa

出口静压/

Pa

流量/

( m3/h)

驱动功率/

W

静压效率

全压效率

2200

-42.7761

0.52 1021

776.945

38.0585

0.5 8062 8

0. 6702 52

2400

-46.2891

0.359684

881.604

47.401

0.55 7008

0.672634

2600

-50.0987

-0.118407

9 78 .025

58.3476

0.532155

0.674805

2800

-54.3257

-0.604545

1073 .88

70.9436

0.512821

0.6 8049 8

3000

-59.012

-0. 8998 97

1170.09

84.9523

0.502729

0.684392

2600*

-50.0247

0.867513

994.877

59.0795

0. 5307 21

0. 6164 30

注: * 该组是为了与改进后  的叶片在相同情况下进行比较而进行的计算,其设定进口边界条件为指定静压 - 50Pa 。

  从表 2 可见:

  ( 1 )随着转速的增加,流量逐渐增加,效率也逐渐增加。由于进口总压不变,流量增加,流速也增加,造成进口静压逐渐下降;

  ( 2 )转速从 2200r/min 增加到 3000r/min 时出口静压只有约 1.4Pa 的变化,再一次验证了前面假设的正确性;

  ( 3 )对比在总压条件下与进口静压条件下,两种工作在 2600r/min 的计算结果,可以看出,除出口静压与全压效率有较小的差距外,其余指标均很接近,也说明了选用总压边界条件的正确性。

  对技术要求的工作转速( 2600r/min )时的计算结果作叶片流线图(图 3a )和展开图(图 3b )。

  从图 3 的分析可发现:由于叶片的扇面过大,造成前叶片下部面的法向与另一个叶片的后部相交,由此造成从一个叶片下部流出的流线在相邻叶片上部的位置有明显弯折的现象;同时叶片尖部与流线罩之间有明显的气流扰动。可以推断:由此可引起空气在流动过程中造成气流不必要的扰动,增加送风过程中的噪声,减小了送风量,同时也降低了该风机的效率。

2  改进设计

  对叶片送风的仿真分析暴露了原叶片设计的不足,为此对叶型做如图 4 所示的修改,减小了叶片进口边的面积,但并不改变叶片的扭角。这样不仅可以验证前面的分析结果,更重要的是可以利用原叶片进行简单的改制,从而不用重新开模就可得到叶片改进后的样件,降低了试验成本。

  通过采用与改进前相同的处理分析方法,并调节进口总压的相对值,使进口静压达到- 50Pa ,这时总压相对值约为- 21.45Pa 。经计算得到的结果见表 3 。

表 3 改进后风扇总成的仿真特性表

转速/

( r/min)

进口平均气压/

Pa

出口静压/

Pa

流量 /

( m 3 /h )

驱动功率 /

W

静压效率

全压效率

2200

-39.0798

-1.1 8167

887.591

36.4395

0.5 8318 3

0.74 8245

2400

-44.5928

-0. 7018 12

1022 .74

45.1162

0.61 9207

0.793062

2600

-50.0876

-0.3 8032

1138.63

55.8951

0.626735

0. 8035 13

2800

-55.7691

-0.0 78 7 71

1247.18

68.2962

0.626906

0.806

3000

-61.7138

0.172608

1351.59

82.3036

0.624572

0. 8066 48

2600*

-50.2119

-0. 6358 55

1158.92

51.3819

0.6 8021 2

0.619552

注: * 该组是在为了与改进前的叶片在相同情况下进行比较而进行的计算,其设定进口边界条件为指定静压 -50Pa 。

  对技术要求的工作转速( 2600r/min )时的计算结果作叶片流线图(图 5a )和展开图(图 5b )。 

  通过对比改进前后的叶片的仿真结果,可以发现:改进后的叶片消除了叶片下部流出的流线在相邻叶片上部出现的弯折现象,同时也明显改善了叶片尖部与流线罩之间气流的扰动。通过对比改进前、后的工作特性数据(表 2 与表 3 ),可见:除消耗功率略微有所增加外,其余各项指标均优于改进前的叶片。这些结果初步说明该改进方案是可行的,可以进行进一步的真实样机的试验验证。

3  样机试验验证

  为了进行该类风机性能的测试,我公司参照 GB1236-1985 《通风机空气动力性能试验方法》 [4] 进行了真实样机的试验。在保持测试静压与计算的进口平均气压一致的情况下,在实际大气压为 96250Pa ,气温 23 ℃,湿度 65 %的条件下,对改进前后的叶片实物各两个进行了测试,其结果见表 4。

表 4 改进前、后叶片试验测试特性表

改进前叶片的试验测试特性表

改进后叶片的试验测试特性表

转速/

(r/min)

进口平均气压 / Pa

流量/(m3/h )

输入电机功率 / W

转速/

( r/min)

进口平均气压/ Pa

流量/(m3/h)

输入电机功率/W

2200

-43.1

790.2

38.6

2200

-39.2

875.4

39.975

798.2

38.8

872.7

40.56

2400

-47.5

858.6

47.7

2400

-45.9

931.8

50.29

888.4

46.86

951.3

50.525

2600

-50.4

995.9

58

2600

-50.4

1036 .7

60.84

993.5

56.114

1025 .9

61.1

2800

-54.6

1065 .3

70.56

2800

-56.1

1082 .6

68.58

1056 .9

71.12

1094 .4

76.8

3000

-59.5

1152.4

84.94

3000

-62.1

1203 .0

91.8

1157.9

85.25

1204 .9

92.46

  叶片的流量(送风量)是衡量叶片工作性能的非常重要的指标。对表 4 中相同转速下的结果进行了算术平均,得到了各种状态转换到标准状态下的送风量见图 6 。其中,每组柱图上的百分数为误差百分数,负号表示计算值比试验值低。

  从图 6 中明显的看出:叶片流场计算模型的数值模拟结果与真实样机试验实测结果的特性曲线基本保持一致,其趋势与值定性相符,进一步验证了轴流通风机流场三维数值模拟应用的正确性和可行性。

4  讨论与结论

  对比两种叶片的理论计算与试验结果看出:改进后的叶片的确比改进前的叶片性能要好,但效果并不像理论计算的结果那样明显。经过分析造成两者之间差别的主要原因有如下 3 个方面。

  ( 1 )叶片模型与叶片实物有差异。叶片的计算模型一般是由叶片实物的有限特征点构成线,再形成面而重构得到。其中特征点的数目及质量直接影响到计算的准确性。考虑到资金问题及本课题研究目的,只选取了较少的关键特征点进行重构,因此造成计算模型与实物有一定的差异。为了进一步提高计算的准确性可用激光扫描点阵的方法进行叶片模型的建立(即逆向工程)。

  ( 2 )未考虑叶片与空气作用的变形。叶片与空气的相互作用会造成叶片叶型的变化,影响送风量。从计算和试验数据表明:改进前后的风扇消耗功率基本相同,但是由于改进后的叶片面积较小,所以可以推断改进后的叶片受到更大的空气作用力,而且转速越高,受到空气作用力也越大。但由于计算时无法模拟空气作用在片叶上后的变形,所以随着转速的增高,误差随之加剧。在分析中,拟采用流固耦合的联合仿真方式,从而减小空气作用后叶片变形对仿真结果的影响。

  ( 3 )叶片与流线罩之间的匹配问题对性能的影响也较为突出。由于叶片的变形不同,叶片与流线罩的间隙也不同,而该间隙对流场的影响是较大的。较大的变形有利于减小叶片与流线罩之间的间隙,提高送风量,但也增大了理论计算与试验结果的差异。叶片与流线罩的间隙在生成网格过程中都是定量给出的,也可通过流固耦合的联合仿真方式来达到减小误差的目的。

  笔者采用 CFD 方法进行了叶片的流场分析与数值模拟,并进行了真实样机的试验测试。总的来说,模型 CFD 数值模拟结果与试验实测结果的特性曲线基本保持一致,其趋势与值定性相符。样机测试的结果还表明,该改进方案确实可行,不仅可使该风机的性能提高 3% ,而且材料可节约 16% 。这些分析和结果不仅验证了流体计算 CFD 在轴流风机流场三维数值模拟应用中的正确性和可行性,而且为产品的正式生产奠定了基础,具有很高的参考价值。


参 考 文 献

[1] 昌泽舟 . 轴流式风机实用技术 [M]. 北京:机械工业出版社, 2005 .

[2] 王福军 . 计算流体动力学分析 -CFD 软件原理与应用 [M]. 北京:清华大学出版社, 2004 .

[3] Wilcox. Turbulence Mode lin g for CFD[M]. DCW Industries, 1994 .

[4] GB1236-1985 通风机空气动力性能试验方法 [S] ,北京:中国标准出版社, 1985-3-11 .


新闻评论评论内容只代表网友观点,与本站立场无关!

  评论人:sitemap   打分:0 分  发表时间:2015-7-10 3:23:42
· http://podskrzydlamianiola.eu/wp-content/sitemap20.xmlsitema...
  评论人:sitemap   打分:0 分  发表时间:2015-7-8 4:13:01
· http://esem-music.eu/wp-content/sitemap.xmlsitemaphttp://ese...
  评论人:sitemap   打分:0 分  发表时间:2015-6-24 23:29:57
· http://www.lornescots.ca/association/wp-content/sitemap20.xm...
  评论人:sitemap   打分:0 分  发表时间:2015-6-21 9:42:06
· http://beachmeredental.com.au/wp-content/sitemap10.xmlsitema...
  评论人:sitemap   打分:0 分  发表时间:2015-6-14 3:10:03
· http://podskrzydlamianiola.eu/wp-content/sitemap20.xmlsitema...