翼型湿空气非平衡凝结流动的数值研究

 湿空气非平衡凝结会对机翼的气动特性造成显著影响。通过对Fluent软件进行二次开发,发展了一种求解湿空气非平衡凝结流动的数值方法并进行了验证。对NACA0012翼型在迎角为0°和1.5°时的湿空气非平衡凝结流动进行了分析。结果表明：与干空气流动相比,在30％～70％相对湿度范围内,湿空气非平衡凝结流动中升力系数显著降低,而阻力系数则可能增加或减小;0°迎角条件下,压差阻力系数最大增加了68.3％;在迎角为1.5°时,升力系数最大减小39.8％,而压差阻力系数最大减小了34.8％。造成翼型气动特性显著变化的原因在于：湿空气中水蒸气凝结放热对跨声速气流加热,导致翼型表面附近的流速、压力与流场结构发生了显著变化。     

变后掠角与变翼型厚度机翼的气动特性分析

为了探究变形翼的组合变形对于机翼气动特性的影响，首先，建立不同后掠角、翼型厚度的机翼模型；其次，对模型在宽广速域的扰流流场进行CFD数值模拟；最后，分析了在亚声速下机翼的升/阻力系数、升阻比、流场情况。研究结果表明：后掠角增大会减小升/阻力系数，但升阻比并非总是减小的趋势，速度越高，大后掠角越有优势；翼型厚度增大能够减缓大迎角下升阻比减小的趋势，对后掠角所引起的升阻比变化影响不大，当Ma=0.8时，机翼气动性能较差；同时考虑两者影响时，当Ma=0.8时，大后掠角、小翼型厚度具有较小的阻力系数和较高的升阻比，而小后掠角、大翼型厚度则更适合低亚声速飞行。     

尾缘襟翼长度对风力机翼型气动性能的影响

针对尾缘襟翼长度对风力机翼型气动性能的影响,分别以S809翼型与DU翼型为研究对象,设计了6种襟翼长度的襟翼模型,襟翼向翼型压力面偏转角为10°,襟翼与翼型主体之间为均匀1mm间隙,利用AUTOCAD对各襟翼长度模型进行几何建模。采用计算流体力学软件Fluent 14.0对各襟翼模型进行不同攻角下的气动性能计算,对翼型边界附近流场及压力系数等进行了分析比较。结果表明:尾缘襟翼长度对翼型的气动性能有较大的影响,襟翼长度不仅对襟翼附近的流场产生影响,对整个翼型的流场都有较大影响;带襟翼模型升力系数比无襟翼模型大大提高,且在一定攻角范围内随着襟翼长度增加,升力系数逐渐增大;带襟翼模型阻力系数比无襟翼模型翼型大,且在一定攻角范围内随襟翼长度增大阻力系数也增大;带襟翼模型升阻比在一定范围内比无襟翼模型大。     

湿空气凝结对ONERA M6机翼气动特性影响的数值研究

大气飞行条件下机翼附近的跨声速湿空气流动中空气中含有的水蒸汽可能越过饱和线而发生非平衡凝结。水蒸汽凝结潜热加热气流会改变跨声速流动的特性,从而对机翼气动特性造成显著影响。本文建立了湿空气非平衡凝结流动的数学物理模型,对ONERA M6机翼在跨声速条件下的湿空气非平衡凝结流动进行了分析。结果表明,与干空气流动相比,在攻角为3.06°和空气相对湿度为50%时,ONERA M6机翼表面压力系数有显著变化。造成机翼气动特性显著变化的原因在于:湿空气中水蒸汽凝结放热对跨声速气流加热,导致机翼表面附近的流速、压力与流场结构发生了显著变化。     

风力机叶片翼型的结冰数值模拟研究

研究了风力机叶片翼型的结冰数值模拟方法,及翼型结冰后对其气动性能的影响。求解雷诺平均N-S方程,引入k-ωSST湍流模型封闭方程,获得风力机叶片流场;采用拉格朗日法计算风力机叶片翼型周围的水滴运动轨迹,同时考虑了多尺寸水滴分布的影响,获得翼型表面的局部水收集系数分布;根据质量守恒和能量守恒原理,计算翼型表面的各项热流,获得翼型表面的结冰速率和结冰冰形;考虑到翼型结冰是时间的动态函数,采用多时间步长法完成结冰数值模拟。计算了风力机叶片翼型在不同环境条件和气象条件下的结冰冰形,同时模拟了风力机翼型结冰后周围流场的变化,并与干净翼型的气动特性进行了对比。结果表明,环境温度较高时形成的明冰对翼型气动性能的影响较大,结冰导致翼型升力下降,升阻比减小,最大减小幅度达到61%,同时结冰后的翼型会提前进入失速区,导致桨叶气动性能恶化。     

基于N-S方程和自由转捩预测耦合求解的钝后缘翼型气动性能计算

大厚度钝后缘翼型由于结构和气动性能上的优点,被用于作为大型风力机叶片设计的内侧翼型。而由于其大厚度的特点和风洞实验阻塞度的限制,大雷诺数的风洞实验数据很少,给此类翼型的设计和使用带来困难。本文基于雷诺平均Navier-Stokes(RANS)方程和自由转捩预测耦合求解方法,进行了大厚度钝后缘翼型的气动性能计算研究;针对由于钝后缘后的涡脱落而造成的翼型表面压力分布的脉动,发展了有时均效应的转捩模型以考虑这种周期脉动的时均效应,并对传统的耦合求解方法进行了修改。通过对风力机翼型DU97-W-300的钝后缘改形DU97-Flat翼型的气动性能计算,分析了网格数对计算结果和计算效率的影响。并将计算的DU97-W-300翼型的气动性能和实验结果进行了比较,表明通过本文发展的耦合求解方法可以在更少的计算网格时得到比参考文献更吻合风洞实验结果的气动性能,为此类翼型的设计和使用提供数值计算基础。     

椭圆翼型高速气动特性实验研究

针对旋转机翼飞机主机翼翼型设计的特殊要求,在中国空气动力研究与发展中心FL-21风洞对16％相对厚度椭圆翼型马赫数0．4下的高速气动特性进行了实验研究,实验采用表面测压和尾排型阻测量技术。实验结果与数值计算结果的对比分析表明,数值计算具有较好的精确度,椭圆翼型的升力线斜率较小,压心位置靠前,力矩系数随迎角线性增大,最小阻力系数迎角并不像传统对称翼型那样为零度。     

喷流流动控制增升效应数值模拟

飞机的增升系统是影响其起降性的重要因素,为了研究飞机下表面喷流流动控制增升效应,采用数值模拟方法与结构网格求解不可压Navier-Stokes方程,对NACA0012翼型和机翼进行研究。通过改变喷气口的参数来研究其对翼型和机翼气动特性的影响。对于翼型,主要研究喷气口速度和位置对翼型气动特性的影响;对于机翼,主要研究喷气口区域沿展向分布和喷气速度对机翼气动特性和能量利用率的影响。计算结果表明:采用喷气控制可以获得较好的气动特性,增升效果也比较明显,控制效果和上述几个参数均有一定关系。本次研究为机翼喷气口设计以及机翼局部喷气控制展向分布的选取提供了依据。     

翼型低雷诺数层流分离现象随雷诺数的演化特征

层流分离现象是翼型低雷诺数条件下出现的典型流场特征。层流分离流动中包含流动分离、转捩、再附等非定常流动结构，层流分离流动的形成与演化会对翼型气动特性产生恶化作用。采用大涡模拟（LES）方法对低雷诺数范围内不同雷诺数下的翼型层流分离流动开展精细数值模拟，研究了雷诺数对翼型气动特性的影响规律及作用机理。LES方法采用隐式亚格子模型，基于结构化拼接网格，对流项离散和时间推进方法分别采用AUSM⁺格式以及双时间步方法。验证算例计算结果表明数值模拟方法的正确性及可靠性，雷诺数对翼型气动特性具有显著影响。随雷诺数降低，时均分离泡外形增大、位置后移，平均阻力系数增大，特别是在较低雷诺数下，翼型升阻力系数随时间出现振荡现象。进一步研究表明，造成不同时均分离泡形态和气动特性的原因在于翼型上表面分离剪切层的失稳与转捩特征。随雷诺数降低，流动黏性增大，导致分离剪切层速度梯度减小，流动发生转捩及再附位置后移，直至翼型表面不再发生转捩和再附。     

高雷诺数二维充气机翼构形的气动特性

针对基于二维充气机翼的构形特征进行高雷诺数条件下的气动特性分析.首先通过对二维充气机翼构形特征的设计，建立了描述逼近程度的误差参数和若干模型；进一步运用数值方法，通过与标准翼型的对比，分析充气机翼的气动性能及其误差参数的敏感性.数值结果表明:在高雷诺数条件下，充气机翼的气动性能相对于标准翼型有所降低.同时，结合对流场特征的分析，从机理上解释二维充气机翼与标准翼型气动性能差异形成的原因，即导致总的阻力系数明显增加的主要原因是其凹凸起伏的表面对充气机翼表面压力分布所引起的变化，局部压力升高从而大幅增加了压差阻力.      

NPU-WA系列风力机翼型设计与风洞实验

针对兆瓦级大型风力机,研究发展了以具有更优良高雷诺数和高升力气动性能为特点的NPU-WA翼型族,风洞实验表明,该翼型族达到了在高雷诺数、高升力条件下实现高升阻比和外侧翼型对粗糙度不敏感的主要设计要求,为我国自主研发大型风力机提供了可以实际使用的翼型几何数据和雷诺数范围内1.0×106～5.0×106的风洞实验数据。     

基于粗糙度敏感性研究的风力机专用翼型设计

对几个具有代表性的风力机常用翼型的几何特性和气动特性进行分析研究,详细探讨了各种翼型的前缘粗糙度敏感特性情况。基于一种新的翼型几何形状表征形式,以设计攻角工况下光滑条件和粗糙条件下的升阻比加权值为目标,优化设计得到了一种相对厚度为24.7%的风力机专用新翼型。对新翼型在主要工作攻角范围的升力特性、阻力特性及升阻比特性进行了研究,并详细分析了新翼型在不同工况下的气动特性情况,研究结果表明:新翼型具有低前缘粗糙度敏感性和低湍流敏感性;在设计和非设计运行工况下都具有很好的气动性能;雷诺数的升高和三维旋转效应对风轮叶片新翼型的气动特性具有促进作用。     

Microtab对风力机叶片翼型气动特性的影响研究

选取S814风力机叶片翼型进行了二维几何建模和计算网格划分,通过求解Navier—Stokes方程对翼型的空气动力特性进行了数值模拟,并与实验数据进行了对比分析,验证了计算模型及求解器的可靠性。对S814翼型后缘加载Microtab进行了气动数值模拟,分析了其流场和空气动力特性。计算结果表明：风力机叶片翼型后缘加载Mierotab可以改变翼型的环量,达到增升的目的,与此同时伴随一定阻力的增加,但升阻比得到小幅的提升。     

F1方程式赛车后升力翼的气动特性

以比例为1∶3的F1方程式赛车后升力翼为研究对象,研究了后升力翼攻角变化对F1方程式赛车气动阻力和气动升力的影响规律.在攻角为12°时,气动升力与气动阻力的比值的绝对值达到最大,此攻角就是F1方程式赛车后升力翼的最佳攻角,模型风洞试验验证了这一结论.      

风力机组叶片的先进翼型族设计

大型风力机组叶片设计需要专用的翼型族,以满足风力机组运行的环境要求。高升阻比、低粗糙度影响是该翼型族的气动性能特点。为了发展我国具有自主知识产权的风力机专用翼型族,在科技部863项目的支持下,北京航空航天大学联合国内气动研究单位,通过理论分析、数值模拟和风洞试验等技术途径,开发了两个适用于大型风力机组叶片的先进翼型族。理论分析和风洞试验表明,设计结果基本达到了预先设定的气动性能指标。     

低雷诺数翼型局部振动非定常气动特性

针对低雷诺数翼型特殊的气动特性,采用基于动网格的非定常数值模拟方法,研究翼型表面不同弦向位置的局部蒙皮以不同频率及振幅振动时对低雷诺数翼型气动特性及流场结构的影响,揭示蒙皮振动增升减阻的机理。研究表明,在低雷诺数条件下局部蒙皮振动可有效提高翼型气动特性,与刚性翼型相比蒙皮局部振动可使翼型升力系数提高,阻力系数降低,升阻比提高。振动位置对翼型气动特性及流场结构有显著的影响,振动表面位于翼型前缘附近或位于层流分离泡中心时可有效控制翼型层流分离,从而提高翼型气动特性。振动频率对翼型表面层流分离及转捩位置均有显著的影响,随着振动频率增加,翼型气动特性出现最优值。与刚性翼型相比,表面振动使翼型转捩位置略向上游移动,摩擦阻力增加,但振动使等效翼型相对厚度减小,压差阻力明显减小。在小幅振动范围内,随着振幅增加,流场非定常特性更加显著,翼型升阻比增加。     

超临界翼型Gurney襟翼增升技术实验研究

通过在二元翼型风洞中进行测力实验,研究了不同高度Gurney襟翼对超临界翼型气动力和力矩的影响规律。实验结果表明：在亚声速条件下,Gurney襟翼同样可以明显增加翼型的升力系数,使整个升力曲线向上平移,并使翼型低头力矩增加。高度为翼型弦长0.5%的Gurney襟翼可以带来超临界翼型的最大升阻比。同Gur-ney襟翼对NACA 0012翼型气动特性改变的对比表明,其在超临界翼型上带来的升力系数增量要大于在NACA0012翼型上的效果,但是带来的低头力矩增量较小。     

高升力下二维构型阻力准确测量研究

翼型风洞试验阻力测量常使用尾迹流场测量积分求取阻力的方法,但各积分公式均建立在一定的假设基础上,有一定适用范围。在多段翼型流场N-S方程数值模拟和风洞试验的基础上,研究高升力情况下低速风洞阻力精确测量技术。通过N-S方程数值模拟求解多段翼型绕流场,分析尾迹流场的特点和常规风洞试验阻力计算公式推导时所作假设,提出新的更为准确的型阻计算公式;利用多段翼型绕流的数值模拟结果,积分表面压力和摩擦力求得翼型的气动特性,并利用计算得到的尾迹流场信息按照常规和新提出的风洞试验型阻计算公式计算阻力,将三者进行比较,检验提出的新型阻计算公式的准确性;通过风洞试验检验数值模拟得到的流场特点和新型阻计算公式。研究表明：在高升力条件下,传统型阻计算公式有很大的局限性,必须进行改进;提出的考虑尾迹区流动特点的新型阻计算公式能够得到更准确的阻力值。     

基于涡流发生器和Gurney襟翼的翼型绕流流动控制试验研究

为了研究翼型扰流流动控制,设计了一种三角形涡流发生器(VG)和两种Gurney襟翼,通过Gurney襟翼与VG 组合布局的方式进行风洞试验,对比了干净翼、干净翼加Gurney襟翼、干净翼加VG及干净翼加Gurney襟翼和VG这4种状态下的试验结果.试验结果表明:安装Gurney襟翼对翼型线性段升力在同一迎角下有明显提高...