- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
2003’全国流体力学青年研讨会论文集
壁面展向周期振动对雷诺应力输运的影响‘
黄伟希许春晓崔桂香张兆顺
清华大学工程力学系,北京100084
摘要:本文对壁面在展向周期运动的槽道湍流进行了直接数值模拟,通过对施加控制胡始阶段
和统计稳定阶段的雷诺应力输运方程的分析,研究了壁面展向周期振动的减阻机理。发现了压
力变形项在减阻中的关键作用,揭示出该项的持续减小是湍流受到抑制的主要原因。
关键词:展向壁面振动.壁湍流,减阻书1L堙
0引言
以减阻为目的的湍流的主动控制因其广阔的应用前景而成为湍流研究的热点。壁面在展
向的周期振动可使壁面摩擦阻力有效降低,而且作为一种主动控制方法,它不需要信息的反
馈,近年来获得了人们的重视。
自从上世纪90年代以来,人们利用计算和实验研究了展向壁面振动对壁湍流的影响…,
都得到了持续的减阻效果。此外还得到了一些共同的结论.如近壁区平均流向速度梯度减小,
缓冲区增厚,对数区上移,雷诺应力和湍动能减小和峰值外移等.这些结论反映了壁面展向
周期振动后的湍流边界层在统计性质上的变化,但是它们却不能揭示减阻的根本原因。
研究表明.近壁区流向涡结构在壁湍流的产生和维持中起到了重要作用,它所诱发的高
速流体的下扫运动是产生壁面摩擦阻力的重要原因,所以通过干扰近壁流向涡结构可以达到
碱阻的目的,壁面展向周期振动正是这样一种控制手段。目前对壁面展向周期振动减阻机理
的解释主要有两种模型:一种是根据壁面展向周期振动对流向涡的扭曲而提出的展向涡模型,
展向涡使得近壁处平均流向速度梯度减小,并由于流向涡受扭曲而带来流向涡脉动量的减小,
使得触发事件的强度碱小,从而带来壁面阻力的减小”1;另一种是近壁流向涡和高低速条带相
互作用的模型,根据条件平均得知壁面展向周期振动将使低速条带插入流向涡下方,并使高
速条带向流向涡上方运动,从而在雷诺应力的象限分析中一三象限的比重大大增加,这就造
成了雷诺应力的减小”。…。壁面展向振动的槽道湍流也可以看成是非定常的三维边界层,Moin
等”1利用直接数值模拟研究了施加展向压力梯度而形成的三维边界层的性质,并对雷诺应力输
运方程进行了平衡分析,提出压力变形项在减阻中起关键的作用。
本文利用直接数值模拟对壁面展向周期振动的槽道湍流进行研究,通过对振动初始阶段
和统计定常阶段雷诺应力输运方程的分析,研究振动对湍流的影响,从而进一步揭示湍流受
到抑制和壁面切应力减小的内在机理。
1数值方法
本文对壁面在展向作用期振动的槽道湍流进行了直接数值模拟,为了进行对比,上壁面
在展向进行周期运动,下壁面保持静止。计算域和坐标系如图l所示,图中x,Y,z分别表示流
向、横向和展向坐标,“,v,W为相应方向的速度分量。控制方程为不可压缩牛顿流体的N-S
方程:
翌+v.vv=一v口+胛2V (1)
V.V=0 (2)
0072032)
自然科学基金资助项目(批准号:l
黄伟希、许舂A-g-:壁面展向周期振动对雷诺应力输运的彰响
式中V为速度向量,P
3自E,JJ,V为运动粘性系数。在流向和展向采用周期边界条件.在壁
面上采用无滑移条件。上壁面展向周期振动通过边界条件直接引入:
(3)
和展向采用Fourier变换,在横向采用Chebyshev变换。在时间方向采用三阶精度时间分裂格
式。采用谱截断方法(也称3/2规则)消除混淆误差。具体方法参见【10】。
图1计算域及坐标系
幽2减阻率髓时间的演化
本文在计算时保持流向流量不变,雷诺数Re.=U。日/v=2666,其中u。为截面平均速度。
计算域在并,Y,z方向的尺寸为4FⅣ×2H×2FH。三个方向的网格数分别取64×65×64。为验
证程序的正确
文档评论(0)