二维弹性圆柱涡致振动的尾涡模态

对均匀来流情况下Re数为200的二维单自由度弹性支撑圆柱涡致振动问题进行了数值模拟。采用时空有限元方法求解二维不可压N-S方程,利用新型显示积分方法求解圆柱振动方程,并用弹簧动网格技术实现网格更新。针对低折合阻尼比涡致振动,分析了不同圆柱自振频率下升力系数、阻力系数以及圆柱位移等参数的变化趋势,成功地捕捉了"锁定"、"拍"和"相位转换"等现象。讨论了圆柱脱落尾涡模态随时间的演变,观测到新的2P+2S尾涡模态,当圆柱振动幅值最大时,尾涡在下游不同位置可能出现2S、2P+2S和2P模态。     

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

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

旋转FGMs层合圆柱壳行波模态频率分析

为推进功能梯度（FGMs）材料在发动机、导弹和火箭等领域的应用，旨在研究旋转FGMs层合圆柱壳的行波模态频率。采用Voigt模型和Sigmoid体积分数描述FGMs层合圆柱壳的材料属性，考虑科里奥利力、离心惯性力、环向初应力以及热内力推导了FGMs层合圆柱壳的能量表达式。采用切比雪夫正交多项式构造位移容许函数，建立了任意边界旋转FGMs层合圆柱壳的模态频率方程，并探讨了组分含量、夹层厚度、温度梯度和弹簧刚度系数等对FGMs层合圆柱壳模态频率的影响。结果表明:夹层厚度相比Sigmoid体积分数对模态频率的变化更为敏感；高旋转短薄壁圆柱壳相比长薄壁圆柱壳对边界条件和失稳现象的影响更为敏感；轴向弹簧相比其他弹簧对模态频率的影响更大。      

竖直平板自然对流强化换热的实验

研究各种强化换热结构对竖直平板自然对流换热的强化作用.通过比较圆柱扰流、圆锥扰流、三棱柱扰流、三棱锥扰流、四棱柱扰流、四棱锥扰流和直角三角翼纵向涡发生器LVG(Longitudinal vortex gener-ator)等对自然对流换热的强化效果,发现:单种换热结构中直角三角翼纵向涡发生器的强化换热效果最好,组合结构中直角三角翼叉排与圆柱扰流的混合排列换热效果最好,其能将竖直平板自然对流换热效果提高24.%.      

周期性变形圆柱绕流特征的数值分析

应用有限体积法求解不可压缩二维N-S方程，对在均匀来流中圆柱-椭圆柱周期性变形运动进行了数值模拟。对于两种典型状态，即雷诺数为40和100，本文描述了圆柱变形运动所引发的非定常流动随时间变化的过程，揭示了变形过程中柱体附近流场的变化主要取决于由壁面压力梯度和壁面的切向运动加速度所产生的壁面涡量流，分析了不同变形频率和不同变形量对升力和阻力、旋涡形成及其运动的影响，并发现了频率锁定现象。     

斜劈式翼型扰流柱冷却通道流动与换热数值研究

基于NACA翼型和斜劈式纵向涡发生器结构设计了一种新型斜劈式翼型扰流柱及其在冷却通道内的排布形式,采用数值模拟方法对含有该种翼型扰流柱的涡轮叶片冷却通道流动换热特性进行了研究,得到了流场和温度场的详细特征及含有斜劈式翼型扰流柱冷却通道的换热能力。分析结果表明:斜劈式翼型扰流柱有良好的强化换热效果。该翼型扰流柱的斜劈式结构可以在其尾缘后方形成持久而稳定的二次流纵向涡,有助于加强流场的扰动作用,从而增强通道换热能力。      

谱单元法及其在多圆柱绕流分析中的应用

系统阐明谱单元方法,基于谱单元方法对低雷诺数Re=200时不同间距下的顺排两圆柱和Re=150正方形排列的四圆柱绕流及其阻力系数、升力系数等进行数值模拟。研究比较不同间距比L/D(两圆柱圆心距离与圆柱直径之比)对两圆柱和四圆柱绕流的影响,计算分析了涡量图分布、平均阻力系数和斯托罗哈数随间距比的变化。研究表明,间距比对顺排两圆柱和正方形四圆柱绕流影响显著;顺排两圆柱绕流存在临界间距比,在Re=200时临界间距比约为3.6。正方形排列四圆柱存在三种流态。当流场从一种流动形态变成另外一种流动形态时,力学参数发生显著变化,在某些间距比区间内出现骤升或骤降现象。     

圆柱跨声速绕流中的激波/湍流相互作用大涡模拟研究

采用大涡模拟方法研究了圆柱跨声速绕流中的激波/湍流相互作用问题,来流马赫数M∞取为0.75,基于圆柱直径D的雷诺数为2×105。计算结果表明,圆柱分离点处出现一道斜激波,并且以与涡脱泻Strouhal数一致的特征频率向上游传播。激波运动导致流场中出现反对称的流动模态,剪切层中压力信号的功率谱曲线中存在0.4、-1和-5次方的斜率关系,剪切层中的剪切应力角约为0°,脉动速度以流向脉动速度为主,并且沿剪切层的大尺度结构组织性减小。     

压气机叶栅叶顶间隙流的动力学模态分解

为了研究低速孤立压气机叶栅叶顶间隙流的非定常运动型态,采用大涡模拟(LES)技术对流场进行数值模拟,并运用动力学模态分解(DMD)技术对x/c=1.0677弦长处S3截面叶顶二次流速度场进行模态分析。结果表明:动力学模态分解(DMD)能够得到速度场的定常模态和不同频率下的主要振荡特征。同时揭示了频率为f1=824.9Hz的一阶振荡主要表现为泄漏涡的低频周期性生成与退化;而频率为f2=9807.2Hz的二阶振荡是由叶顶移动端壁形成剪切层内的不稳定造成的,主要表现为泄漏涡的振荡与分裂。     

统计能量分析（SEA）在复杂结构分析中的低频限研究

借助于AutoSEA软件和实验手段，对类似于飞机机身舱段的圆柱壳结构的内部声压级进行了计算和实验测量，对结果作比较，并通过建立不同SEA模型，探讨复杂结构的SEA分析中子系统模态数和下限分析频率的关系，借用功率流分析，发现复杂结构的SEA低频分析，受给目标子系统主要传递能量的子系统的支配，这对SEA下限分析频率的确定和低频SEA模型的划分具有一定的意义。     

转子叶片气弹稳定性与强迫响应分析

针对由叶轮机械转子叶片流致振动引起的叶片高循环疲劳失效的突出问题,采用数值模拟方法对叶片流致振动特性进 行预测分析。常见的流致振动问题包括气弹稳定性和强迫响应,分别以压气机和高压涡轮转子叶片为例给出了气弹稳定性和强 迫振动响应的分析方法及其流程。对于压气机叶片,采用能量法和特征值法进行了气弹稳定性计算。结果表明：气弹稳定性结果 相互吻合,在高阶模态下的气动阻尼比在低阶模态下的更小,根据振型特征确定危险点位于叶尖近前缘和近尾缘部位,在该部位 可能导致叶片掉角;对于高压涡轮叶片,分析了其强迫振动响应特征。结果表明：叶片非定常气动激励主要来源于上游导叶,在设 计状态下其频率距离转子叶片第5阶模态频率较近,并对危险模态阶次的共振应力进行分析,在0.14%阻尼下振动应力最大可达 134 MPa。     

典型及衍生激波针构型的减阻降热流动特性

为探索激波针对超声速钝头飞行器进行减阻降热时的更优衍生构型,采用数值模拟方法对3种典型单扰流物构型、6种双扰流物构型、两类多扰流物构型和钝锥型激波针的流动特性进行了研究,认为加装激波针后的几何本质相当于“镂空式”的锥型钝头体。模拟结果显示：激波针头部扰流物相对直径较大时,减阻率随激波针相对长度的变化曲线没有明显的峰值点,而是存在一个变动幅度很小的峰值段,且相对直径在0.3～0.4左右时减阻效果最佳;典型激波针的最大减阻率约为50%,采用双扰流物构型时略有提升;中部增加多个扰流物时减阻率随扰流物数量增多而增大,最大减阻率超过60%,但气动加热问题较严重。相比而言,钝锥型激波针减阻降温的综合性能最好,最大减阻率可达60%左右,降温率约为7%。     

基于POD和DMD方法的跨声速抖振模态分析

跨声速抖振现象是由于非定常跨声速流动中激波的自激振荡而引起的结构强迫振荡,这种现象在跨声速飞行器中普遍存在,对飞机的结构强度和疲劳寿命有不利影响。基于模态分解的分析方法是进一步发展抖振控制手段的有效工具。本文通过两类典型模态分析方法（本征正交分解（POD）和动态模态分解（DMD））对OAT15A翼型的跨声速抖振现象进行分析,通过对模态频率、翼面压力分布、流场重构误差等方面的研究,将两种模态分解方法进行对比。发现基于频率特征的DMD方法能够准确捕捉抖振的临界稳定特征和抖振主频的典型模态,同时能够更准确反映流场变量在激波间断附近随时间的变化过程;而POD方法尽管在流场重构时具有较小的总体误差,但对激波附近压强随时间的变化历程拟合较差。     

斜、直圆柱绕流的CFD模拟

采用CFD软件CFX5.5对直圆柱和斜圆柱绕流进行数值模拟,计算了直、斜置圆柱的阻力系数、升力系数、周向平均压力系数及其脉动压力系数的均方根值,并与试验结果进行了比较.计算表明,斜圆柱绕流与直圆柱绕流有很大的差异:斜圆柱绕流的平均阻力系数减小,平均升力系数不等于零,升力系数的频域除包含有传统的卡门涡脱频率外,还包含有许多低频成份.     

应用等离子体激励器抑制串列双圆柱噪声的数值仿真研究

针对等离子体激励下的串列双圆柱绕流噪声抑制问题，通过将等离子体体积力模型、脱落涡模拟、声比拟理论等技术相结合的数值模拟方法，研究不同来流速度下等离子体激励器安装位置对双圆柱分离流形态控制与远场噪声抑制效果的影响。结果表明，当所施加的等离子体激励位于圆柱流动分离点附近时，控制措施可有效减小分离涡尺度和湍流强度，并显著降低远场监测点的总声压级。随着来流速度增大，等离子体激励器的降噪效果增强，同时最优安装位置前移。当来流速度达到55m/s时获得最优降噪效果，其远场监测点声压级频谱峰值和总声压级分别降低11.5dB和8.3dB。而随着来流速度的进一步增大，等离子体激励器的降噪效果逐渐减弱。所得结果对于等离子体流动控制抑制串列圆柱噪声的实际应用有一定指导意义。     

基于 DDES 算法的有扰流片腔体气动噪声分析北大核心CSCD

为了进一步了解腔体流场的物理特性以及锯齿形扰流片这种被动控制方法对腔体噪声振荡的抑制效果,对马赫数0.85下、长深比为5∶1的矩形腔体进行了气动声学分析。分析过程采用了基于Spalart-Allmaras单方程模型的延迟分离涡模拟(DDES)的CFD方法。计算结果以声压级(SPL)形式与英国QinetiQ基地进行的实验研究结果进行比较。计算准确地预测出了第二阶和第三阶主导模态的声调噪声幅值,其误差均不超过2dB。为了验证腔体前缘安装扰流片这种广泛应用的腔体噪声被动控制方法,计算重点分析了带有锯齿形扰流片的腔体结构噪声等级。计算分析结果与空腔结果对比分析后认为前缘锯齿形扰流片可以降低第二阶主导模态的峰值噪声等级10dB以上,其它频率下的腔体噪声等级可降低接近5dB,降噪效果较好。     

局部贴敷黏弹性阻尼层圆柱壳振动频率与阻尼有限元分析

在引入黏弹性材料复模量模型和考虑圆柱壳弹性边界的基础上,研发了有限元程序来解决考虑频率依赖性的复合圆柱壳振动频率与阻尼的非线性计算问题。创建了一种4节点24自由度复合壳单元来模拟局部贴敷黏弹性阻尼层圆柱壳的力学行为,推导了单元的刚度及质量矩阵。提出由6个弹簧组组成的周向变刚度弹性约束模型来模拟圆柱壳底部的弹性边界条件。确定了复合圆柱壳的动力学有限元方程,并描述了用特征向量增值法求解其振动频率与阻尼的迭代计算过程。对贴敷ZN-1自由阻尼层圆柱壳进行了实例计算,结果表明:该算法的计算结果与实验测试值误差最大为3.69%,另外黏弹性材料的频率依赖性对复合结构固有频率影响小于0.01%,而对模态损耗因子影响最大为109.47%。     

螺旋桨飞机近地飞行时的纵向稳定性

根据某型飞机气动试验数据,分析了地面效应和螺旋桨滑流增量对飞机静稳定性的影响;根据小扰动线性化方程,计算了不同飞行条件下的飞机模态特性,分析了地面效应和螺旋桨滑流对飞机动稳定性的影响。结果表明：地面效应将使飞机纵向运动频率增加,而阻尼比降低;地效区内滑流增量使短周期震荡频率和阻尼比均增加,也使长周期阻尼比增加,而对长周期模态频率的影响不大;无地效时,滑流将使长短周期两种模态的频率均降低,而阻尼比增加。     

基于DDES算法的有扰流片腔体气动噪声分析

为了进一步了解腔体流场的物理特性以及锯齿形扰流片这种被动控制方法对腔体噪声振荡的抑制效果,对马赫数0.85下、长深比为5∶1的矩形腔体进行了气动声学分析。分析过程采用了基于Spalart-Allmaras单方程模型的延迟分离涡模拟(DDES)的CFD方法。计算结果以声压级(SPL)形式与英国QinetiQ基地进行的实验研究结果进行比较。计算准确地预测出了第二阶和第三阶主导模态的声调噪声幅值,其误差均不超过2dB。为了验证腔体前缘安装扰流片这种广泛应用的腔体噪声被动控制方法,计算重点分析了带有锯齿形扰流片的腔体结构噪声等级。计算分析结果与空腔结果对比分析后认为前缘锯齿形扰流片可以降低第二阶主导模态的峰值噪声等级10dB以上,其它频率下的腔体噪声等级可降低接近5dB,降噪效果较好。     

基于混合RANS/LES方法与FW-H方程的气动声学计算研究

在气动噪声的数值计算过程中,非定常流动的求解精度对声学计算结果有着重要的影响.以机体噪声计算的标准算例双圆柱绕流气动噪声问题为研究对象,采用基于非线性k-ε湍流模型的限制数值尺度(LNS)方法对双圆柱绕流进行了数值模拟,将计算得到的气动特性和流动特征与相应的试验结果进行了对比分析.为了求解远场观测点处的气动噪声,在精确求解双圆柱绕流流动的基础上结合基于FW-H(Ffowcs Williams-Hawkings)方程的声类比方法进行数值计算,并通过圆柱体的展向相关性对计算结果进行了修正,将得到的最终结果与相应的声学试验结果进行了对比,两者吻合良好,表明该数值方法是准确、可靠的.