基于马蹄涡模型的飞机尾流介电常数分布的快速建模方法

为便于研究飞机尾流电磁散射特性,提出了一种飞机尾流介电常数分布的快速建模方法。根据马蹄涡模型得到飞机尾流的诱导速度场,并利用Euler方程和等熵流方法求出了尾流的热力学参数分布,进而得到介电常数分布。同时,提供了一种将针对一种飞机计算出的结果外推到另一种飞机上的方法,避免了重复计算。该方法具有计算量小、涉及参数少的特点,是一种快速简便地模拟飞机尾流介电常数分布的方法。     

七叶大侧斜螺旋桨空泡特性及梢涡演化数值模拟研究

为研究七叶大侧斜螺旋桨尾流场及梢涡特性,本文基于DDES（延迟分离涡方法）建立了螺旋桨空化流场数值预报模型。为验证所建立数值模型的准确性,进行了多套不同尺度网格的不确定性分析,同时将非定常流动中螺旋桨（VP1304）螺旋桨空泡及梢涡特性计算结果与试验结果进行了对比。随后基于该模型对七叶大侧斜螺旋桨的尾流场及梢涡特性进行了数值分析。计算结果表明：结本文所建立数值模型精度较高,可以准确地捕捉到梢涡空泡及梢涡尾流场特性;同时DDES方法相比RANS方法在对复杂湍流流动的捕捉能力更强,更适用于螺旋桨梢涡捕捉;尾流场网格加密对尾流场模拟及准确捕捉梢涡十分重要,但对螺旋桨水动力性能影响不大;尾流区域的轴向速度场可以分为加速流动区和自由流动区,进速系数越小,自由流动区与加速流动区之间的界限向外扩张越明显;E1619桨在重载工况有明显的梢涡空泡产生,而轻载工况空泡面积较小且无梢涡空泡发生;梢涡在向下游发展过程中会发生相互融合,进速系数越大,融合发生的越晚,梢涡强度也越小;七叶大侧斜螺旋桨在尾流区域会产生一个分支涡,分支涡起始于吸力弯曲面下缘,且与梢涡呈一定的夹角,进速系数越大,夹角越小。     

有无翼尖涡扩散器的民航机翼数值模拟计算

采用ANSYS FLUENT软件对有无翼尖涡扩散器机翼翼尖涡的形成和消散进行基于Realizable k-ε涡粘模型的数值模拟计算,以探究翼尖涡扩散器对尾流的影响。通过对比分析两种机翼的静压系数、轴向涡量、速度矢量可知:加装翼尖涡扩散器不仅可以改变翼尖处静压系数的分布,使升力系数增大;还可以阻挡下翼面高压气流向上翼面流动,将翼尖涡分隔成涡量相反的四个涡,这四个涡在流向下游的过程中彼此消耗能量,最终减小了尾流的范围。     

旋翼地面效应的气动建模与特性

地面对旋翼气动特性影响明显，且导致旋翼流场更加复杂。为分析地面效应下的旋翼桨尖涡和流场变化特性，基于涡面和无滑移边界条件，求解第2类Fredholm方程获得地面涡面矢量分布，且将涡面矢量按涡扩散方程扩散到流体中，建立考虑黏性效应的地面气动模型，并耦合非定常面元/黏性涡粒子混合法以体现旋翼桨叶气动特性和旋翼尾迹的非定常效应，构建旋翼地面效应气动分析方法。通过计算Lynx尾桨地面效应下的性能和桨尖涡轨迹，并计算Maryland大学模型旋翼和NASA缩比旋翼地面效应下的垂向、径向速度分布，且与试验和CFD计算结果对比，验证了本文方法能较好捕捉地面效应下的旋翼尾迹变化特性和复杂速度场特性，且结果表明本文方法能较好模拟地面效应下旋翼桨尖涡的收缩、扩散、井喷、地面射流等物理现象。     

适用于整体求解吊舱推进器的定常面元法

为了方便开展吊舱推进器设计和混合式CRP(Contra-Rotating propeller)的面元法分析,提出了一种将吊舱和螺旋桨作为整体进行定常面元法求解的可行办法,对相关的面元法基本积分公式进行了分析,说明了定常面元法整体求解的结果与实际结果之间的关系。从吊舱桨尾流区的实际流场状况和计算效果出发,分析了整体求解吊舱推进器的面元法尾涡模型,并对此种尾涡模型的修正予以详细说明。在进行相关的理论分析后,编制了整体求解吊舱推进器的定常面元法程序,对吊舱推进器进行计算分析,并与试验结果和RANS(Reynolds-averaged Navier-Stokes)方法计算结果进行了对比,最后对本文所提出的定常整体面元法求解方法进行了误差分析。计算结果表明,整体面元法求解吊舱推进器的误差与试验值相比在5%以内;在计算耗时方面,将吊舱和螺旋桨分开的迭代求解面元法是整体求解面元法的1.7倍。     

反弹自消散系统尾涡控制机制研究

为了研究加快尾涡消散的机理,以获得更合理的尾流间隔,进行了壁面反弹实验的数值模拟计算。实验采取在模拟机翼后缘安置反弹面的方法,诱发反弹二次涡与主涡相交出现不稳定性,以加快尾流衰减。数值计算运用N-S方程,对雷诺应力项采用Realizableκ-ε涡粘模型进行计算,分析经二次涡干涉后尾涡轴向涡量衰减、涡核下沉运动等参数的变化,以探究反弹二次涡对飞机主涡消散的演变影响。实验结果表明,反弹涡对主涡衰减有明显的促进作用。该实验为人工干预缩减机场尾流间隔的研究奠定了理论基础。     

三维钝体结构风压分布与绕流特性的大涡模拟研究

针对风场作用下三维钝体结构的绕流特性与空间风压分布特性,采用大涡模型对其进行了数值模拟研究,以平均风速和脉动风速作为入口边界条件,在稳态分析的基础上对9种不同来流方向的工况进行了瞬态动力分析,再现了钝体绕流的冲撞、分离、重附着等现象,揭示了钝体绕流主、次涡循环区的形成、发展机理,提出了空间风压场分布的一般规律。研究表明:钝体绕流对流场风速剖面影响的强弱程度依次为尾流区、屋顶区、来流区,且对中心区域流场风速剖面的影响要大于两侧流场,并由中心区域向两侧逐渐减小;来流区钝体绕流特性与空间风压分布特性较屋顶区和尾流区稳定,尾流区涡旋的非对称生成、脱落对钝体左右两侧流场的冲撞、分离、重附着等具有重要的影响;来流方向与钝体呈45°或135°夹角时,对尾流风场的影响范围最小,但随着来流方向逐渐垂直于钝体,尾流影响区域也逐渐增大。     

舰船空气尾流场对直升机着舰的影响研究

舰船空气尾流场是直升机舰上起降时的主要环境条件,对直升机的操稳性能及飞行安全有很大的影响。通过某型舰模风洞试验的PIV结果,介绍了舰船空气尾流场特性,如:等速度场和截面流线图、舰船机库脱体涡等,分析了下冲气流、涡流区及开/关机库大门等对直升机着舰的影响,对保障直升机飞行安全有重要价值。     

水平轴风力机偏航与俯仰尾流特性

风力机的尾流干扰效应是影响风场下游风力机性能和风场出力的主要因素，通过调整风轮姿态可有效改善尾流干扰效应，但同时也引入了更为复杂的非定常气动特性和尾流场分布。为了研究风轮姿态控制的尾流特性，基于升力线和自由涡尾迹模型建立了水平轴风力机气动和尾流仿真模型，以NREL 5MW风力机为研究对象，在额定风速条件下，研究了叶轮偏航和俯仰工况下的尾流结构特性。在此基础上，对存在风切变和塔架影响时的风力机偏航和俯仰工况产生的尾流干扰进行了分析。研究结果表明：风轮的偏航和俯仰会引起尾流涡结构失稳提前、远尾流区湍流度增大，同时会导致尾流中心发生偏转和尾流作用区域减小，并加快速度亏损恢复，降低尾流效应对下游的干扰。风切变会产生类似风力机俯仰的效果，使得尾流抬升，作用区域减小。塔架对近尾流区的干扰较为显著，而俯仰控制会降低塔架对尾流的干扰。     

风-浪-舰耦合对舰船尾流场及直升机气动的影响

为研究波浪和舰船运动对高海况下舰载直升机起降的影响，建立了一套耦合风-浪-舰和直升机的多相流数值模拟方法。使用动态流体作用模拟舰船运动，采用动量源方法模拟直升机旋翼，对3～6级海况下ONRT舰船摇荡运动时的尾流场非定常气动特性进行了研究，并计算了UH-60A直升机旋翼气动载荷的时间历程。研究表明：6级海况下流场速度分量和旋翼气动载荷的波动达到5级海况的2倍，而3～5级海况之间的差别较小；悬停跟进时直升机不断出没于舰船运动导致的尾流中速度和涡量不同的区域，使得旋翼气动载荷有极强的非定常特性；在6级海况下，旋翼拉力与波浪-舰船耦合运动同周期的波动达到20%,并伴随着穿越尾流剪切层时产生的5～10 Hz的小扰动。     

基于黏性涡模型的旋翼流场数值方法

 建立了一种适用于旋翼非定常流场特性分析的黏性涡数值方法。在该方法中:流场中的大尺度涡被离散为若干微小的涡元,通过求解涡量-速度形式的Navier-Stokes方程模拟涡元的输运等过程;黏性扩散效应采用高精度的粒子强度交换法进行计算,而桨叶附着涡以及新生涡环量采用了Weissinger-L升力面理论进行求解;为显著提高计算效率,在诱导速度及其梯度的计算中还引入了快速多极子算法(FMM)。应用上述方法,对悬停和前飞状态下的多个旋翼流场算例进行了计算,通过对比旋翼尾迹涡量特征和诱导速度分布等,验证了该方法的有效性。此外,还将本方法与旋翼计算流体力学(CFD)方法及传统的自由尾迹方法进行了比较,结果表明黏性涡方法在兼顾效率的同时,还能够更好地捕捉旋翼尾迹运动。     

虚拟体法数值模拟双圆柱强迫振荡系统

用虚拟体方法研究了均匀流及尾流中圆柱横向振荡的流动特征.单圆柱结果与相关试验吻合,得到了"锁定"区域,发现涡量跳跃现象.对小间距双圆柱的数值模拟发现,在"锁定"区域附近随激励频率的增加流场特征会出现若干显著的变化,包括涡量的跳跃,涡模态的改变,正负涡量带的形成和消失.在大间距情况下发现下游圆柱的振荡对上游尾涡脱落没有明显的影响,下游圆柱尾流场没有发现"锁定"现象.     

绕水翼超空化阶段汽液混合体动态特性的实验观测

为深入了解超空化流动机理,采用高速录像和数字粒子图像测速系统（DPIV）对绕hydronautics水翼的超空化流场进行了观测。结果表明：在两相共存超空化阶段,水汽混合相和汽相的分布决定了超空化区域速度场和涡量场的分布特性;水汽混合区的速度相对较低,而汽相区则与主流区有着相近的速度分布,而且空化区与主流区交界面处有着较大的速度梯度;在水汽混合区与主流区的交界处形成了上下涡带,其涡量分布表现出了明显的时间依赖性。     

垂直交叉双圆柱绕流强迫振动数值模拟

采用虚拟体法数值模拟了平面单圆柱绕流强迫振动,以及空间垂直交叉双圆柱绕流下游圆柱强迫振动流场。验证了单圆柱强迫振动中的锁定状态以及相位突变现象,从而证实了该数值方法模拟振动流场的可靠性。研究了在雷诺数Re=150、间距为5倍圆柱直径、下游圆柱按正弦曲线振动时,对尾流场的影响。下游圆柱两端尾流场在振动的作用下,涡的横向间距增大,而中心尾流场由于受上游圆柱尾流的影响而保持原先的状态。下游圆柱在锁定区的振动使尾流场变得稳定,表现在流向二次涡结构的减少,以及尾流场的速度场频谱趋向单一化。     

旋翼桨-涡干扰气动特性计算及参数影响研究

建立了一个适用于旋翼桨-涡干扰(BVI)气动特性分析的数值方法。在该方法中,控制方程采用惯性坐标系下的非定常Euler方程,以适合于分析BVI流场的特点;为便于前飞流场分区求解和信息传递,使用了嵌套网格方法;对于流场中涡线的模拟,建立了一种适用于有限体积格式的涡引入方法——广义网格速度法,以简化计算。应用上述方法对旋翼BVI流场进行了计算,并与可得到的试验数据进行对比,验证了方法的有效性。通过对比桨叶弦向位置压强的变化得出,在旋翼BVI过程中,气动载荷主要来自桨叶10%弦长内的前缘部分的压强突变。文中还进一步分析了涡强、干扰距离和干扰角度对BVI气动特性的影响。本文的BVI数值计算表明:当涡接近桨叶前缘时,升力达到最大;而涡位于后缘位置时,诱导速度改变桨叶环量分布,破坏了后缘的库塔条件,但随着涡的远离,桨叶表面环量开始重构。     

烧蚀环境下的圆柱绕流计算模型

基于多学科理论建立了高温高压气流环境下圆柱烧蚀、剥蚀的数理模型。利用离散涡方法计算流场与圆柱表面压力分布;采用三方程烧蚀模型计算热化学控制机制下的烧蚀速率,并将其与扩散控制机制下的结果相比较以确定烧蚀量;引入颗粒轨道模型求解剥蚀颗粒的运动。基于该模型对高温高压燃烧室内圆柱形烧蚀试件的绕流场、烧蚀量以及剥蚀颗粒的运动轨迹进行了编程计算分析。研究表明:低雷诺数条件下烧蚀环境对绕流流场的影响不大,而较高雷诺数条件下烧蚀绕流流场与无烧蚀绕流流场的涡街分布存在显著差别;在低雷诺数条件下流场对烧蚀速率的影响甚微,较高雷诺数条件下烧蚀速率略大于无流场烧蚀情况;涡量的存在决定着剥蚀颗粒的运动轨迹与分布,进而影响尾流热场与烧蚀速率。该模型可为热防护材料的烧蚀实验提供参考,并为烧蚀、剥蚀过程的非线性分析提供数据支持。     

某机载吊舱的气动特性分析

为满足航空遥感试验研究需求,设计并分析了一种飞机机腹吊舱。文中采用CFD软件对飞机机腹加装的吊舱进行气动外形特性分析,首先利用Icem软件对机载吊舱进行结构性网格划分,通过Flunet软件模拟飞机在多种高空环境下的飞行特性,仿真其受压力变化分布,并计算飞机重心的变化,提出吊舱设计应考虑压力调节需求和机腹结构加强需求,并评估飞机重心变化对飞行安全的影响。经飞行试验证明,该吊舱设计合理,可为同类产品设计提供参考。     

基于数值模拟的飞机近场尾涡特征参数计算

尾流间隔是飞机起降和飞行过程的安全保障。近场尾涡特征参数是尾涡远场演变和消散研究的基础，是缩减尾流间隔以提升机场及空域容量的理论依据。采用雷诺平均法数值模拟A320飞机尾涡的近场演化过程，对近场尾涡特征参数的计算方法进行了研究。结果表明：与Hallock-Burnham模型相比，Lamb-Oseen模型能更准确地描述近场尾涡速度分布；采用平均椭圆法计算的近场尾涡涡核半径比采用传统方法计算的精度更高；数值模拟中，A320机型更适合采用3～12 m的平均环量值来描述近场尾涡强度。     

大型水洞中螺旋桨尾流场PIV测试研究

以标准桨DTMB-P4119为试验对象,开展了螺旋桨尾流场测试分析。试验捕捉到了清晰的螺旋桨梢涡和尾涡的流动结构以及梢涡、尾涡片在尾流中的演化。分析比较了三种载况下的尾流特征,利用捕捉到的梢涡涡心坐标计算了部分涡线的螺距角。文章将测试结果与LDV测试数据进行了比较分析,结果表明,PIV与LDV对微观流动结构的测量取得了比较一致的结果。     

双缝激波诱导推力矢量喷管的非定常流动特性研究

采用分离涡模拟方法对双缝激波诱导推力矢量喷管的三维流场进行数值模拟,分析原始流场的流线、压力系数分布。运用傅里叶模态分解技术对喷管对称面压力系数进行模态分解和时间演化。结果表明:双缝激波诱导推力矢量喷管的流场中包含两个模态,一阶模态主要表现为下游二次流口激波处和尾流上剪切层处的压力振荡,二阶模态主要表现为剪切层处涡的脱落。本文对研究矢量喷管非定常流场特定频率下对应的潜在动力学信息具有重要意义。