首页/文章/ 详情

【经典教材翻译】18-机翼的激波膨胀技术

1月前浏览2006

二维超音速流动的一般解可以认为是均匀流动、激波和膨胀波的组合。前几节中讲解了解决正常激波、斜激波和膨胀波的方法,可以用来解决超音速机翼流动。这称为激波-膨胀技术Shock-Expansion Technique。它在有激波的地方使用激波关系,在有膨胀的地方使用普朗特-迈耶膨胀关系,并且假定两者之间的流动是均匀的。


平板机翼


     

   
 

图43:零攻角的平板机翼

考虑将一个薄平板放置在超音速流中,如图43所示。对于零攻角,平板上任何地方都没有流动转向。因此,上下表面的压力很大均匀。阻力和升力都为零。(注意:这种方法忽略了摩擦损失)。

如果平板以等于α的攻角放置,那么板的上下面之间就会产生压力差。(图44)

图44:攻角的平板机翼

流动将上表面的前缘视为凸角。因此产生普朗特-迈耶膨胀。在尾缘,流动通过激波被压缩回水平流方向。下表面的前缘是激波,因为它形成了一个凹角。流动通过膨胀扇从尾缘流出。仔细观察尾缘后的流动显示有两个气流 - 一个是通过上表面的膨胀和激波处理的,另一个是通过压力面上类似特性处理的气体。两个激波的强度不一样。因此,气流的密度和温度不同。然而,在尾缘,压力和流动角度被等化。这在尾缘引起一个滑流。从所示的压力分布,可以计算升力和阻力为

其中c是弦长。请注意,这种阻力不是由粘性引起的,如不可压缩或亚音速流动中那样。这是由波(激波和膨胀)引起的,这些波是超音速流动特有的。这是超音速波阻力的一个例子。同样,粘性、表面摩擦效应被忽略了。升力和阻力可以表示为系数CL和CD。这些是通过用升力或阻力作用区域的来无量纲化相应的力得到的,其中A是升力或阻力作用的面积。这也可以写成。因此,

其中Cplower和Cpupper是下表面和上表面的压力系数。


钻石形机翼


     

   
 

图45:钻石形机翼周围的流动

考虑一个典型的超音速流动机翼,即图45所示的钻石形机翼。零攻角下的流动产生如图所示的特征。在前缘,每一侧都有一个激波。然后在最大厚度处有膨胀波。流动通过另一个激波系统离开尾缘。

流动关于机翼中心线对称,升力为零。但是存在阻力,由下式给出,

可以概括这个机翼结果并为阻力和升力开发一个公式。如果机翼的半楔角为δw,θ是机翼任何一侧的取向,那么现在可以对每一侧的压力求和,以确定升力和阻力系数如下

用每一侧的Cp表示,

同样对于阻力,

用每一侧的Cp表示,


激波和膨胀波之间的相互作用


     

   

在钻石形机翼的情况下,远场流动中可能会发生激波和膨胀之间的相互作用。一般来说,这些对表面压力的影响不大,不需要对整个流场进行完整分析就可以得到结果。图46显示了流动相互作用。

图46:膨胀波和激波之间的相互作用


薄机翼理论


     

   

激波-膨胀技术是准确的,但它需要单独求和表面组分,因此是一个复杂的数值过程。如果考虑机翼很薄,并且只考虑弱斜波激波,可以找到封闭形式的解。这是超音速薄机翼理论Supersonic Thin Aerofoil Theory。假设流动在通过机翼时只从自由流方向轻微偏转。因此,由于激波属于弱激波类别,那么流动中任何地方的压力变化由下式给出,

假设压力P与P∞不远, 并且机翼上的局部马赫数与M∞不远, 使上述方程简化为

将所有压力参照P∞, 并且将流动方向参照自由流,

这给出,

因此,我们有一个简单的表达式来计算放置在流动中的任何表面上的Cp。有趣的特点是Cp取决于局部流动倾斜θ,而不管什么特性导致流动转向。


攻角下平板机翼周围的流动


     

   

考虑上面处理过的平板机翼。流动在两个表面上以α角倾斜。因此,

升力和阻力系数由下式给出,

并注意对于小α,cos(α) = 1,则


钻石形机翼


     

   

对于钻石形机翼,对于激波后的流动,

对于膨胀波后的流动,

使用这个方程时,对于压缩使用正号,对于膨胀使用负号。阻力系数由下式给出,

这可以写成


任意机翼


     

   
 

图47:任意机翼周围的流动

放置在超音速流动中的一般机翼_图47_,可以认为是α攻角的平板、由于弯曲引起的流动偏转αc(x),和由于厚度引起的流动偏转h(x)的线性叠加。对于这种情况,

其中是截面的平均属性。


二阶理论


     

   

上面显示的近似理论是一阶的,因为它只保留了Cp展开中涉及θ的第一个重要项。作为更准确的替代方案,布塞曼提供了一个二阶理论,其中也包括θ2项。根据这个理论,

这也写成

对于这个方程,如果流动正在经历压缩,则使用正号,如果流动正在膨胀,则使用负号。

注意,系数C1和C2仅是马赫数和γ的函数。

因此,有三种方法来计算转向超音速流动中的压力。其中,激波-膨胀技术最准确。其余两种仅适用于小的流动角度变化。布塞曼的方法Busemann's method更准确。


通过取消波减少阻力


     

   

在超音速流动中,波是阻力的主要来源。通过从系统中“移除”波,可以减少阻力。以激波冲击实体墙为例,这会产生反射激波。后者发生是为了使流动与墙面平行。如果墙本身在点(0)处以θ角旋转(图48),那么流动就会沿着墙面流动,就不需要反射激波了。这种现象也可以解释为在(0)处产生一个膨胀波,它抵消了反射激波。现在系统没有波,因此没有波阻力。

图48:波的取消

建立在波取消理念上的一个巧妙装置是布塞曼双翼飞机(图49)。几何形状和来流马赫数被如此安排,以产生一个完全对称的激波系统,在出口没有任何波。从理论上讲,这给出了零波阻力。如果布塞曼飞机在非设计条件下运行,如图50所示,出口流动不是无波的。有结果的波阻力。

图49:设计条件下的布塞曼双翼飞机

图50:非设计条件下的布塞曼双翼飞机


来源:CFD饭圈
理论科普
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2024-09-08
最近编辑:1月前
CFD饭圈
硕士 分享CFD文章,感谢关注
获赞 22粉丝 23文章 378课程 0
点赞
收藏
作者推荐

一个用了ANSYS十年的过来人,领悟出的25条仿真建议

使用了ANSYS软件十年了仿真项目完成了一个又一个从一个小白到现在的“老油条”作为过来人总结了以下25条仿真技巧有些东西点到为止需要亲身体验过才懂不知道下面有几条内容能引起您的共鸣?1.理解边界条件:准确的边界条件至关重要。花时间理解你的系统的物理行为,以定义适当的边界。2.网格质量很重要:在网格化上投入努力。高质量的网格化可以带来更准确的结果。ANSYS提供了多种网格化选项;探索并明智选择。3.收敛性检查:在仿真期间定期监控收敛性。如果出现收敛问题,调整设置或网格。4.利用远程求解:对于大型仿真,利用远程求解或高性能计算(HPC)资源以获得更快的结果。 5.参数分析:使用ANSYS Workbench的参数分析功能探索不同的设计方案,以了解设计变化的影响。6.材料属性精度:确保输入准确的材料属性。微小的偏差可能会显著影响仿真结果。7.使用自适应网格:ANSYS的自适应网格化功能可以在需要的地方动态细化网格,提高准确性而不增加过多的计算成本。8.预处理清理:在运行仿真之前,彻底检查模型是否有几何错误或重复,以避免分析过程中出现问题。9.探索多物理场仿真:ANSYS为多物理场仿真提供了强大的能力。尝试将不同的物理场耦合进行综合分析。10.后处理掌握:了解等值线图、路径图和结果探针等后处理工具有助于从仿真结果中提取有价值的见解。 11.优化工具:ANSYS提供了优化工具。熟悉这些工具,以便有效地根据各种标准优化设计。12.考虑时间依赖性分析:在具有瞬态行为的场景中,时间依赖性分析可以更准确地理解系统动态。13.文档和笔记:维护仿真设置和所做的假设的详细文档,以供将来参考或协作。14.利用远程可视化:对于复杂模型,使用远程可视化工具与结果互动,而不会使本地机器过载。15.热分析考虑:在热仿真中,准确考虑材料属性、边界条件和对流效应。 16.合理更换湍流模型:对于流体仿真,研究流动模式、湍流和边界层效应,以优化准确性。不同的湍流模型,仿真出来的结果相差很大,真的,需要下功夫了解各个湍流模型的适用范围,这,就叫做“专业”。17.备份和版本控制:定期备份仿真文件,并考虑使用版本控制以防止数据丢失并有效管理更改。18.模型简化策略:在复杂模型中,考虑简化几何形状或使用对称性以减少计算负载,同时不牺牲准确性。19.耦合场分析:探索ANSYS的耦合场分析能力,其中多个物理现象相互作用,以实现更真实的仿真。20.从教程和论坛中学习:ANSYS提供了广泛的教程和用户论坛。利用这些资源来提高你的技能。 21.考虑误差估计:使用ANSYS误差估计功能来衡量你的仿真结果的准确性。22.与实验数据验证:尽可能地,将仿真结果与实验数据进行验证,以确保准确性并建立对模型的信心。23.协作并寻求反馈:与同行或领域专家进行交流,以获取关于仿真设置和结果的反馈,以实现持续改进。24.保持更新:跟上ANSYS软件更新和新功能,以利用进步进行更好的仿真。25.持续学习:仿真是一个不断发展的领域。保持好奇并不断学习,以跟上新方法和最佳实践。 需要亲身体验过才懂,不知道上面有几条内容能引起您的共鸣? 来源:CFD饭圈

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习 福利任务 兑换礼品
下载APP
联系我们
帮助与反馈