收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

边界元法在气热耦合计算及冷却结构优化中的应用研究

王振峰  
【摘要】: 航空发动机技术已成为衡量国家科技工业水平和综合国力的重要标志,为此各国竞相开展了航空发动机的高技术研究。随着航空工业的发展,对发动机提出了更高的要求,要求它有更高的推重比、更大的单位体积输出功率、更低的燃油消耗率及更好的可靠性能等。提高涡轮前燃气温度是提高发动机性能的有效途径,然而,涡轮进口温度的提高速度远远高于叶片材料耐温性能的发展速度,由此为保证涡轮叶片在高温环境下安全可靠地工作,必须采取有效的冷却技术和热防护措施。涡轮叶片温度场的准确预测可以有效地指导设计人员进行涡轮叶片冷却结构设计,以提高冷却效率、延长叶片工作寿命。本文的主要工作是在哈尔滨工业大学推进理论与技术研究所原有的流体计算程序基础上,搭建具有自主知识产权、适用于气冷涡轮叶栅气热耦合换热计算的数值仿真平台;并进一步深入分析湍流模型、冷却通道入口湍流特性条件以及换热系数实验准则式对气热耦合计算结果的影响。最后应用基于边界元方法(BEM)和有限差分方法(FDM)的气热耦合程序建立了气冷涡轮叶片冷却结构全耦合优化设计体系。 本文首先对边界元法进行了深入的理论研究,详细推导了二维、三维边界元法求解势函数问题的控制方程,并建立边界元法计算固体区域热传导问题的数学模型,编写了二维、三维边界元程序。利用解析解和数值解对边界元程序进行准确性校验,计算结果表明,边界元法仅仅计算模型的表面网格单元,使得数值模拟的前处理和求解时间缩短。同时,由于其具备解析与离散相结合的特性,因此计算结果准确、计算精度较高。对三维边界元程序进行网格适应性研究,研究表明边界元方法对网格质量的依赖性很小。 将边界元计算方法用于燃气涡轮叶栅气热耦合换热计算中,搭建起气热耦合数值仿真平台。并应用此气热耦合数值仿真平台数值模拟了NASA-MarkII涡轮叶栅具有典型流动特点的两个实验工况。分析数值模拟结果得出,三维气热耦合程序能够准确模拟不同流动特性的涡轮叶栅内部流动;叶栅内部温度场计算结果说明,在涡轮叶片表面的流体边界层内部温度梯度较大,传热过程剧烈,在流体边界层内部流体由顺压力梯度到逆压力梯度的过渡区域,由于叶栅内复杂的流动特点,采用B-L代数模型的三维气热耦合程序在预测边界层局部区域流体的传热特性时存在一定偏差。 采用商业软件Fluent和CFX对相同的实验叶栅进行气热耦合计算,着重分析气热耦合计算的网格适应性、冷却通道入口不同湍流特性条件对耦合换热结果的影响以及不同湍流模型对具有流动分离特性的涡轮叶栅的流动、传热及热应力特性的预测能力。计算结果表明,涡轮叶栅流动特性计算结果对边界层网格依赖性很小,对网格质量要求不高;而边界层不同的网格划分对叶栅内传热特性的预测略有差异,其主要影响因素为边界层网格厚度和流向网格加密程度。对于涡轮叶片冷却通道、冷却孔等内部流动,湍流强度完全依赖于上游流动的历史,冷却通道的入口湍流特性对通道内流动发展影响较大,从而冷却通道进口湍流特性条件能否正确给定对涡轮叶栅耦合换热计算结果的准确度具有较大影响。CFX提供的考虑转捩流动特性的k- -SST- -湍流模型能够准确分辨层流及转捩状态,对边界层内复杂流动和传热过程模拟较为准确,但其不足之处为,k- -SST- -湍流模型过高地估计了转捩区域的湍动能而造成转捩区的过大估计,从而使得对流换热系数计算值过大,温度偏高。而三维气热耦合程序采用的修正B-L代数模型能够更为准确模拟不存在分离或分离较小的涡轮叶栅热环境,计算结果显示,B-L代数模型除在叶片边界层转捩区域存在一定误差外,其余位置对涡轮叶栅传热特性的预测要强于CFX提供的k- -SST- -湍流模型计算结果。同时,从涡轮叶片热应力特性分析结果得出,在涡轮叶栅气热耦合计算中温度场预测结果的准确与否对涡轮叶片热应力结果预测、以至涡轮叶片寿命预估和冷却结构设计工作至关重要。 将换热系数实验准则式引入气冷涡轮叶栅气热耦合计算中,在计算过程中考虑了基于气冷涡轮叶片冷却通道不同几何模型和实际工况的修正因子,详细对比分析涡轮叶片内冷通道的不同换热系数实验准则式对耦合计算结果的影响;分析气冷涡轮叶片内冷通道采用、不采用换热系数实验准则式对气热耦合计算结果的影响以及对不同湍流模型预测叶栅流动、传热过程准确程度的影响。计算结果表明,气冷涡轮叶栅气热耦合计算中,涡轮叶片冷却通道、冷却孔等位置采用修正的对流换热系数实验准则式可以得到准确的叶栅流场及叶片温度场分布。换热系数实验准则式的采用绕开了湍流模型对冷却通道流体区域的求解过程,使得湍流模型能够专一求解叶栅主流燃气区域流动和传热过程,从而能够更为准确分析不同湍流模型对叶栅主流区域流动、传热特性的计算能力。同时避免了湍流模型计算多区域流场时求解误差的叠加,并且节省计算时间。 将三维气热耦合程序作为气冷涡轮叶栅数值模拟的主程序,建立了气冷涡轮叶片冷却结构全耦合优化设计体系。对NASA-MarkII涡轮叶栅10个径向冷却通道进行优化,在优化过程中将冷却通道孔径、冷却通道空间位置和流经冷却通道的冷气流量作为设计变量进行参数控制,并且对冷却通道孔径和位置进行必要的约束。优化目标函数包括涡轮叶片温度最大值、温度平均值以及冷却孔冷气流量。优化中通过加权求和的方法将多目标函数转化为单目标函数,寻求涡轮叶片冷却结构的最优解。得益于气热耦合程序中求解固体区域热传导问题的边界元程序模块,使得涡轮叶片冷却结构全耦合优化计算时间缩短,工作量减少。在优化过程中,边界元方法省去了固体区域网格重复生成、重复计算的过程,避免了固体区域网格与流体区域网格的插值误差,从而提高优化效率和问题求解的计算精度。优化结果表明,冷却通道位置以及冷却通道孔径的优化,使得叶片区域温度最大值降低;优化后冷气流量的减少,使得叶片区域平均温度略有提高,但提高幅度很小。另外,优化后冷却通道流体与涡轮叶片的换热受入口段的影响有所降低,从而使得涡轮叶片径向温度梯度减小。


知网文化
【相似文献】
中国期刊全文数据库 前20条
1 李瑞阳,卓宁;L齿型螺旋肋片管表面放热系数沿其周向的变化规律的研究[J];上海理工大学学报;1985年04期
2 陈燕燕;罗二仓;戴巍;;热声回热器流动、传热新模型及其性能的实验表征[J];工程热物理学报;2008年03期
3 侯忠霖;姚山;王廷利;张兴国;金俊泽;;一种铝合金水冷界面换热系数反求方法的研究[J];材料热处理学报;2008年01期
4 朱彦伟;马争胜;高扬;;影响叶片表面换热主要因素的分析[J];硅谷;2008年20期
5 赵振明;吴宏伟;丁水汀;徐国强;;旋转状态下气膜冷却换热系数的实验[J];推进技术;2008年06期
6 李茂德;程惠尔;;过氧化氢溶液表面沸腾换热研究[J];航空动力学报;2009年03期
7 李广超;朱惠人;樊慧明;;双向扩张型孔射流角度对气膜冷却特性影响的实验[J];航空动力学报;2009年05期
8 朱会悦;;水浴除尘器设计[J];中国新技术新产品;2009年14期
9 管小燕;任家隆;;液滴低温气体雾化射流冷却换热理论分析与实验研究[J];现代制造工程;2009年11期
10 文耀普;郭舜;;重力热管换热的试验研究[J];流体机械;1983年09期
11 孙曾润;;大型两相热虹吸管传热性能的实验[J];南京工业大学学报(自然科学版);1983年02期
12 周子成,安晨晖;氟里昂卧式壳管式冷凝器的模拟与实验[J];流体机械;1987年12期
13 倪志军,刘松龄;涡轮叶片型面外换热系数的数值计算[J];航空发动机;1997年02期
14 王邦文,李谋谓,赵永忠,朱起健,王峰丽;钢板强适应横向冷却曲线[J];北京科技大学学报;2003年01期
15 吴玉庭,朱宏晔,任建勋,梁新刚;聚光条件下太阳电池的热电特性分析[J];太阳能学报;2004年03期
16 马兰,朱惠人,屈展,许都纯;涡轮动叶表面气膜冷却换热实验研究[J];机械设计与制造;2005年08期
17 王景智,吴坚;大功率中速柴油机活塞镶圈断裂的有限元分析[J];内燃机工程;2005年04期
18 郭涛;朱惠人;许都纯;白江涛;;热色液晶瞬态测量全表面换热系数的技术研究[J];测控技术;2006年09期
19 张魏;;多孔板换热的数值模拟研究[J];西安航空技术高等专科学校学报;2006年05期
20 夏利江;魏琪;冯志彬;;风冷垂直管吸收器外部冷却过程数值研究[J];山东建筑大学学报;2007年03期
中国重要会议论文全文数据库 前10条
1 何茂刚;王小飞;张颖;;制冷用水平管降膜蒸发器的研究进展及新技术[A];第五届全国制冷空调新技术研讨会论文集[C];2008年
2 蒋毅;张小松;殷勇高;;无机填料的直接蒸发冷却实验研究[A];制冷空调新技术进展——第四届全国制冷空调新技术研讨会论文集[C];2006年
3 成剑;金听祥;郑祖义;;平行流换热器换热性能影响因素的分析[A];第九届全国空调器、电冰箱(柜)及压缩机学术交流会论文集[C];2008年
4 钱强;朱跃钊;廖传华;张维薇;;海水淡化中水平管降膜蒸发的研究[A];2008中国水处理技术研讨会暨第28届年会论文集[C];2008年
5 何国庚;王忠衡;;冰浆流体流动换热的实验研究[A];第十二届全国冷(热)水机组与热泵技术研讨会论文集[C];2005年
6 陈莹;高飞;高冈大造;徐林虓;李维仲;左建国;;翅片管式换热器空气侧性能的数值模拟[A];中国制冷学会2009年学术年会论文集[C];2009年
7 伍兵;陈永;赵克文;杨素波;严学模;李茂林;张大德;王地君;黎建全;张均祥;;攀钢弧形板坯连铸机二冷状况的研究[A];1999中国钢铁年会论文集(上)[C];1999年
8 韩志涛;姚杨;;空气源热泵冷热水机组空气侧换热器换热系数的计算方法[A];全国暖通空调制冷2002年学术年会论文集[C];2002年
9 徐峰;郭烈锦;白博峰;;超临界压力下管内层流换热的数值模拟[A];第六届全国低温与制冷工程大会会议论文集[C];2003年
10 王希靖;包孔;韩晓辉;郭瑞杰;;2A12与5A02搅拌摩擦焊温度场模拟分析比较[A];第六届21世纪中国焊接技术研讨会暨第四届计算机在焊接中的应用交流会论文专刊[C];2004年
中国博士学位论文全文数据库 前10条
1 张超;燃气透平导叶气热耦合实验与数值研究[D];中国科学院研究生院(工程热物理研究所);2012年
2 卜晓兵;Al-Cu与Al-Si-Cu合金凝固微观组织的仿真模拟[D];湖南大学;2013年
3 万飞;立式连续热处理炉带钢加热与冷却数学模型理论研究与应用[D];重庆大学;2011年
4 赵蔚琳;高温热管翅性能及其强化传热过程的研究[D];南京工业大学;2004年
5 陶毓伽;大功率固体激光器冷却研究[D];中国科学院研究生院(工程热物理研究所);2010年
6 王振峰;边界元法在气热耦合计算及冷却结构优化中的应用研究[D];哈尔滨工业大学;2010年
7 王振峰;边界元法在气热耦合计算及冷却结构优化中的应用研究[D];哈尔滨工业大学;2010年
8 田晓虎;平行流冷凝器空气侧换热性能和压降的CFD仿真及试验研究[D];重庆大学;2007年
9 温志伟;基于数值分析的大型同步电机内温度场的研究[D];中国科学院研究生院(电工研究所);2006年
10 全晓军;硅基梯形微通道内水蒸汽凝结换热特性及等温同向喷射流研究[D];上海交通大学;2008年
中国硕士学位论文全文数据库 前10条
1 邢磊;基于瞬态法的界面接触换热系数实验研究[D];大连理工大学;2010年
2 廖铮玮;高强度钢板热成形换热系数估算及实验研究[D];大连理工大学;2013年
3 刘志新;18CrNiMo7-6齿轮热处理模拟及换热系数的测定[D];大连交通大学;2013年
4 许瑾;42CrMo钢热处理过程数值模拟及换热系数的测定[D];大连交通大学;2013年
5 郑泽花;钢淬火换热系数计算及温度与微观组织模拟[D];大连理工大学;2010年
6 周芝龙;铸件—模具界面换热机理的研究[D];重庆大学;2013年
7 常国光;GCr15喷雾淬火的数值模拟及换热系数的计算[D];大连交通大学;2013年
8 李彩文;6xxx铝合金在线淬火换热系数及数值模拟研究[D];中南大学;2011年
9 刘斌;喷淋和浸没对高温钢板冷却特性影响的实验研究[D];重庆大学;2012年
10 胡理中;复杂截面铝合金型材在线淬火过程界面换热研究[D];湖南大学;2013年
中国重要报纸全文数据库 前10条
1 ;铜材让空调更完美[N];中国有色金属报;2005年
2 由航天航空学院郑泉水教授、黄克智院士完成 由材料系黄勇教授等完成 由计算机系张尧学教授等完成 由精仪系殷纯永教授与天津大学叶声华院士等合作完成 由航天航空学院过增元院士等与山东大学合作完成 由精仪系、精华同方光盘股份有限公司副研究员裴京等完成 由土木系聂国教授等与中南大学合作完成;我校2004年度国家科学技术奖部分获奖项目简介[N];新清华;2005年
3 张欣欣 冯妍卉 姜泽毅;干法熄焦节能降耗技术研究新进展[N];中国冶金报;2007年
4 文涛;海尔产品占了一半[N];中国质量报;2000年
5 段金平;让浅层地温能造福于民[N];中国矿业报;2007年
6 魏双林;华北干燥领跑国内干燥机业[N];中国工业报;2005年
7 春华;海尔空调为什么[N];中国信息报;2000年
8 王一娟 ;我国中药产业快速增长[N];经济参考报;2005年
9 李玉峰   ;中药现代化天津屡夺全国冠军[N];天津日报;2005年
10 小于;海尔空调喜获国家节能认证[N];中国机电日报;2000年
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978