变截面三维机翼的气动性能及噪声研究.pdf
《变截面三维机翼的气动性能及噪声研究.pdf》由会员分享,可在线阅读,更多相关《变截面三维机翼的气动性能及噪声研究.pdf(6页珍藏版)》请在咨信网上搜索。
1、第 61 卷 第 8 期Vol.61 No.82023 年 8 月August 2023农业装备与车辆工程AGRICULTURAL EQUIPMENT&VEHICLE ENGINEERING0 引言机翼绕流的现象属于三维的湍流流动,并且是绕流现象中十分常见的问题。对于机翼绕流的研究不仅可用于航空航天业,而且可为风力机以及水利行业提供借鉴。研究机翼绕流对于提高翼型元件的气动性能、降低气动噪声等方面具有十分重要的意义。随着近些年来翼型元件的绕流速度越来越快、计算机技术水平的不断提高,其相关的数值研究也越来越被重视。对于翼型绕流的研究不仅局限于翼型的气动性能,也更多地关注到翼型的气动噪声方面。计算气
2、动噪声的方法目前比较通用的是使用基于 Lighthill 声类比模型1的 FfowcsWilliams-Hawkings 方程2(FW-H 方程)。大涡模拟作为 CFD 中兼具精度性和经济性的数值模拟方案,在近些年来得到了广泛的应用和发展。2015 年,樊艳红3使用基于 WALL 亚格子模型的大涡模拟方法对风力机翼型绕流进行数值模拟,并把模拟结果和实验结果进行对比,发现模拟结果能很好地吻合实验结果,说明使用 WALL 亚格子模型的大涡模拟能准确模拟翼型绕流。本文使用 ANSYSFluent,采用基于 WALL 亚格子模型4的大涡模拟(LES)方法对三维变截面机翼的气动特性机翼气动噪声进行数值模
3、拟研究,旨在为三维变截面机翼绕流方向上的气动性能研究提供一些参考,使对于翼型气动特性和气动噪声的产生机理及现象做进一步的了解和认识。doi:10.3969/j.issn.1673-3142.2023.08.023变截面三维机翼的气动性能及噪声研究曹开强(200093上海市上海理工大学机械工程学院)摘要 为了深入研究接近于真实翼型的变截面三维翼型的气动性能以及噪声表现,构建了一个基于 NACA0018 的变截面三维机翼模型,使用 ANSYSFluent 软件对该模型进行气动性能以及气动噪声的研究。数值仿真模拟采用大涡模拟(LES)方法,使用 WALL 亚格子模型,研究了该三维翼型在攻角为 15、
4、来流速度为 40m/s 情况下的绕流特性。使用基于 Lighthill 声学类比模型的 FW-H 方程进行气动噪声的计算。结果显示,在机翼上表面某点之后开始不断形成分离泡产生流动不稳定性,造成压力和速度波动并产生波动的气动噪声,气动噪声在前缘驻点位置以及后缘位置达到最大值。机翼上表面开始产生分离泡之后,流动出现了在翼型展向位置上的横跨流动,这将加剧流动的不稳定性。关键词 变截面;三维机翼;气动性能;气动噪声;大涡模拟;FW-H 中图分类号 V211.5 文献标志码 A 文章编号 1673-3142(2023)08-0114-06引用格式:曹开强.变截面三维机翼的气动性能及噪声研究 J.农业装备
5、与车辆工程,2023,61(8):114-119.Research on aerodynamic performance and noise of variable-section three-dimensional airfoilCAOKaiqiang(SchoolofMechanicalEngineering,UniversityofShanghaiforScienceandTechnology,Shanghai200093,China)AbstractInordertohaveadeeperunderstandingoftheaerodynamicperformanceandnoisep
6、erformanceofthevariable-section3Dairfoilthatwasclosetotherealairfoil,a3DwingmodelwithvariablesectionbasedonNACA0018wasconstructed,andANSYSFluentsoftwarewasusedtostudytheaerodynamicperformanceandaerodynamicnoiseofthemodel.Thenumericalsimulationadoptedthelargeeddysimulation(LES)methodbasedontheWALLmod
7、eltostudytheflowcharacteristicsofthethree-dimensionalairfoilwhentheangleofattackwas15andtheinflowvelocitywas40m/s.AerodynamicnoisewascalculatedusingtheFW-HequationbasedontheLighthillacousticanalogymodel.Theresultsshowedthat,afteracertainpointontheuppersurfaceofthewing,thecontinuousformationofseparat
8、ionbubblesproducedflowinstability,whichcausedpressureandvelocityfluctuationsandgeneratedfluctuatingaerodynamicnoise.Theaerodynamicnoisereachedthemaximumattheleadingedgestagnationpointandthetrailingedge.Afterseparationbubblesbegantoformontheuppersurfaceoftheairfoil,theflowappearedspanwiseacrosstheair
9、foil,whichwouldexacerbateflowinstability.Key wordsvariablesection;3Dwing;aerodynamicperformance;aerodynamicnoise;LES;FW-H收稿日期:2022-06-17115第 61 卷第 8 期1 数值模拟方法1.1 湍流数值模拟方法LES 是目前适用非常广泛的一种流场数值模拟方法,能根据网格的精密度而自由地调整数值模拟的精度。其相比于直接数值模拟(DNS)需要的网格分辨率和计算量都更少,能节约大量的计算资源,同时又能够获得比雷诺时均模拟(RANS)方法更多的湍流信息,例如大尺度涡流的速度
10、和压强脉动等。因此 LES 是结合目前计算机能力下非常具有应用价值的数值模拟方法,也是CFD的一个热门领域。LES 的基本思想是根据网格尺度把湍流分为大尺度湍流和小尺度湍流5,对大尺度的湍流直接使用 N-S 方程进行求解计算。而小尺度的湍流因为其具有各向同性,并且小涡对大涡的影响在运动方程中体现为类似于雷诺应力的应力项,称之为亚格子雷诺应力,所以在 LES 中利用了亚格子模型模拟小涡对大涡的影响,而不直接求解小尺度涡,以此节省大量计算,使求解变得可能。假设在流场中某瞬态变量为 (x,t),对该变量采用滤波操作,如式(1)、式(2)6:(,)(,)(,)dx tG r txr tr=-(1)(,
11、)dG r xr1=#(2)经过这个滤波操作后,瞬态变量 (x,t)被分解为 2 部分:(x,t)=(x,t)+(x,t),其中(x,t)代表大尺度的变量,(x,t)代表小尺度的变量。本文使用 ANSYS 软件,采用三维、不可压缩的非定常模拟方法进行仿真计算,控制方程的离散方法是有限体积法7。使用大涡模拟方法求解流场信息,LES 亚格子模型采用 WALL 模型4,WALL的优点在于黏性系数不随计算网格而改变,也不受周围湍流情况的影响,目前已经得到越来越多的应用。速度和压力耦合方案选用 SIMPLEC8,空间离散格式为 2 阶中心差分格式,时间项离散格式为 2 阶隐式计算格式。1.2 FW-H
12、方程气动噪声的产生和传播可直接通过求解 N-S 方程得到,但是这就要求计算所使用的网格必须具有很高的精度,这意味着计算量大、耗时长。为了弥补这些缺点,在 Fluent 里采用基于 Lighthill 声学近似模型的 FW-H 方程,该方法可将声音的产生和传播过程分开计算,相比直接求解 N-S 方程减少了工作量,对计算机的要求也降低,使求解变得可能。FW-H 方程是 FfowcsWilliiams 和 Hawkings 在1969 年引入广义函数根据 N-S 方程推导出来的,在静止流体中运动物体发声的控制方程,表示为(,)()()()ctxp x ttufxP nu ufx xH fT1iinn
13、niijjinnijij2222202222222222 22t otodtod-=+-+-+leojj99CC&$0.0(3)式中:ctx1i222222222-波动算子;(,)p x til观测点在 t 时刻的声压值;Tij Lighthill 张量,Tij=-Pij+uiuj-c2ij;H(f)Heaviside 函数。在计算时,需要用到的流场变量通常包括速度分量、压强分布等信息,这些信息可以使用 LES或者 DES 等方法获得。在本文中使用 LES 方法求解流场信息,结合FW-H方程求解流动的声场信息,设置机翼表面为声源面。2 模型建立和网格划分2.1 模型建立本文是基于 NACA00
14、18 对称翼型建立的三维变截面机翼模型,其中大弦长 C1=100mm,小弦长C2=95mm,展向长度 20mm,建立模型如图 1 所示。流体域建立时,将机翼的左右两侧面与流场模型的耦合边界完全贴合,即流体域的宽度等于翼型展向长度为 20mm。以机翼前缘点作为参考,流体域入口到前缘参考点的距离 L1=200mm,流体域进口设置为速度入口条件。流体域出口到前缘参考点的距离 L2=400mm,流体域上边界到前缘参考点的距离 H1和流体域下边界到前缘参考点的距离 H2均为 130mm,可认为上下壁面对翼型附近的流动不产生影响。机翼前缘参考点距流场出口的距离要大于距流场入口的距离,以消除可能存在的回流对
15、于流场计算的影响。设置机翼的迎角为 15。流体域模型的建立如图 2 所示。图 1 变截面三维机翼Fig.1 Variable section 3D wing图 2 三维流体域示意图Fig.2 Schematic diagram of 3D fluid domainL1H2H1L2YX曹开强:变截面三维机翼的气动性能及噪声研究116农业装备与车辆工程 2023 年2.2 网格划分网格划分必须根据所使用的湍流数值模拟方法进行,因为不同的湍流数值模拟方法对网格和边界层的要求不同。LES 对网格要求足够精细,网格分辨率要能满足计算所需捕捉的最小湍流的尺度。在 LES 中,大于网格尺度的湍流将被直接求解
16、,而小于网格尺度的湍流将使用亚格子模型表征其对大涡的影响,所以网格尺度体现了LES的计算精度。不同的数值模拟方法对于边界层的要求也是不一样的,总体而言,边界层的处理方法包括壁面函数法和近壁面模型,本文 LES 要求使用近壁面模型法。不同的边界层处理方法要求不同的 y+值标准,y+作为边界层网格划分的衡量依据,其值直接和边界层的第 1 层网格高度相关,是一个无量纲数,表达式为yu y*o=+(4)式中:u*近壁面摩擦速度;y第 1 层网格距离壁面的距离;o 流体的运动粘度。壁面函数法和近壁面模型对 y+的要求是不同的,近壁面模型方法要求在近壁面内要分布足够细密的网格,通常要求 y+在 1 左右。
17、根据式(2),可大致推算出本文模型所需要的第 1 层网格高度为 0.013mm,边界层设置 15 层,增长率为 1.2,能满足本文 LES 的 y+要求。其它网格尺度为 1mm,网格总数为 460 万,网格质量全部在 0.5 以上且大部分都在 0.7 及以上,能满足本文模型对于网格的要求。网格划分详情如图3所示。3 Fluent 仿真计算使用 Fluent2021R2 大涡模拟对模型进行瞬态仿真,亚网格尺度模型使用 WALE 模型,声学模型使用 FW-H 模型。Fluent 求解器类型设置为采用压力修正算法的压力基求解器,压力速度耦合算法使用 SIMPLEC,空间离散格式为 2 阶中心差分格式
18、,时间项离散格式采用 2 阶隐式算法,设置流场速度入口流速为 40m/s。设置时间步长为 0.00001s,计算时间步数为 10000 步。机翼模型中尺度最小的网格为边界层加密网格,其流向长度为 0.5mm,而流体流速为 40m/s,流体流过一个网格单元需要 0.0000125s,因此本文选取的 0.00001s 能满足计算要求。机翼模型的大弦长为0.1m,流体完整流过机翼需要0.0025s,即 250 个时间步,因此本文计算所采用的时间步数均满足计算要求。求解具体设置如表 1 所示。表 1 求解参数及设置Tab.1 Solving parameters and settings数值求解方法L
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 截面 三维 机翼 气动 性能 噪声 研究
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。