瞬时起动圆柱的初期流动

 本文将ADI方法和Poisson方程直接法结合起来求解涡量-流函数形式的N-S方程,研究瞬时起动圆柱的初期流动(Re=3000,5000,9500)。成功地算出实验观察到的α流动结构、β流动结构,以及这些结构在流动发展过程中交替出现的情况。本文方法效率高,计算结果与实验结果及他人的数值计算结果相符较好。可期望用该方法计算其它复杂流动。     

翼型动态失速的数值研究

用不可压缩流动的求解算法，结合WilcoxDC提出的k－ω模式和k－ωSST湍流模式，对翼型的动态失速进行了数值模拟。通过对典型的振荡翼型轻失速和深失速算例的计算结果分析可以看出：（1）绕动态失速翼型的流场结构十分复杂，轻失速和深失速在流动特性上有很大区别。计算结果显示：轻失速主要是由于后缘分离引起，分离涡的影响范围主要是在后缘附近。而深失速则首先形成很大的前缘分离涡，该分离涡在翼型表面上运动，并诱发出二次分离涡，引起翼型升、阻力系数的显著变化。（2）对于动态失速的翼型绕流，k－ωSST湍流模式是较为有效的，计算出的气动力系数迟带曲线变化趋势与实验结果符合得比较好。     

纳秒等离子体激励控制翼型流动分离机理研究

为研究纳秒介质阻挡放电(NSDBD)等离子体控制翼型流动分离的物理机理,采用已建立的NSDBD唯象学模型耦合非定常Navier-Stokes方程模拟纳秒等离子体对流场的作用。使用非定常雷诺平均NavierStokes方程(URANS)和大涡模拟(LES)两种求解方法,研究纳秒等离子体激励对NACA0015翼型流动分离控制。结果表明:NSDBD等离子体激励促使边界层提前转捩,转捩对控制流动分离起重要作用;NSDBD激励开始时在翼型前缘形成展向涡,展向涡促使分离剪切层失稳并最终进入尾迹,展向涡贴近壁面运动,将外区的高能气流带入近壁区,使上翼面流场结构发生变化,然后翼型前缘流动提前转捩促使流动经过一个小层流分离泡后发生湍流再附,最终在上翼面形成稳定的附着流动。     

NACA0012翼型低雷诺数绕流的实验研究

通过水槽氧气泡流动显示和PIV测速实验研究了NACA0012翼型在雷诺数为8200时的流动特性,重点炎注了翼型绕流结构随迎角的变化。研究发脱：分离点和分离翦切层形成旋涡的位置随迎角的增大而向上游移动,同时翼型上表面流动分离后形成的回流区尺寸随着翼利迎角的增加而增大。当流动再附于翼型上表面时,在再附点附近能够观测到展向涡的三维演化过程,并能观测到展向涡的局部配对现象。     

大迎角水平振动翼型绕流的数值模拟

从涡量流函数形式的N-S方程出发,数值计算了大迎角水平振动翼型的粘性绕流问题,得到近场涡结构及非定常演化过程。探讨了振动频率、振幅及迎角对流场涡结构及翼型升、阻特性的影响。结果表明,振动频率及振幅增加有助于提高翼型的平均升力;过大的迎角不利于升力提高。     

椭圆形控制环量翼型气动力的计算

本文用边界层理论与计及分离尾迹影响的位流理论相结合的方法计算控制环量翼型的绕流,用离散涡模型模拟分离尾迹的影响。由上、下翼面分离点处引入离散涡。只是翼型附近尾迹中的离散涡可以自由移动。离翼型较远的离散涡假设为按来流速度移动。忽略射流出口上游的边界层对射流的影响。这些简化使本文的方法十分有效。与实验比较表明,本文的方法能较准确地计算出翼型的分离点位置和分离点前的压强分布。计算出的升力系数与射流动量系数的关系与实验相符很好。     

翼型近尾迹流中大尺度涡结构的理论模型

根据流动稳定性理论相干结构的相关理论 ,提出了一种二维波模型来描述翼型近尾迹流动中的大尺度涡结构 ,通过采用对近尾迹流动区域分区处理的方法计算了近尾迹区域的大尺度涡结构。计算所得的流线图 ,涡量分布以及雷诺应力分布具有与实验结果相似的分布规律。说明该二维模型能够较好的描述近尾迹流动中这种大尺度涡结构的运动学和动力学特性      

低雷诺数翼型层流分离现象大涡模拟方法

为准确预测翼型低雷诺数条件下出现的层流分离流动现象,发展了基于结构化拼接网格的大涡模拟方法。控制方程为Favre滤波Naver-Stokes方程,并选用多种亚格子模型。空间离散采用AUSM格式以及高阶WENO格式,时间推进采用显式方法和隐式方法。以SD7003翼型在雷诺数60 000及4°迎角下的层流分离流动为研究对象,对比分析了数值格式、亚格子模型、网格尺度对流场预测结果的影响。数值结果表明:划分的计算网格能够有效解析小尺度流动结构,基于隐式亚格子模型、采用AUSM格式和双时间步推进的大涡模拟方法能够准确预测流动分离、转捩、再附等复杂流动现象,计算得到的平均压强系数与雷诺应力分布与文献数据吻合较好。与转捩模式计算结果对比进一步表明:发展的大涡模拟方法能够准确预测翼型低雷诺数层流分离流动非定常演化过程,为下一步的研究工作建立了有效的数值模拟手段。     

使用GAO-YONG湍流方程组计算翼型分离流

采用SIMPLE方法求解GAO-YONG不可压湍流方程组, 对不同来流迎角下的NACA0012翼型绕流结构进行了数值模拟, 给出了翼型绕流分离流结构随迎角的变化特征和翼型在分离绕流中的气动力参数.与实验数据以及大涡模拟结果的比较表明, GAO-YONG不可压湍流方程组能够对翼型绕流的分离点、分离涡形态、表面压力分布、升阻特性做出较好的预测, 能够模拟翼型大攻角分离流动, 计算结果优于FLU-ENT软件中的k-ε、k-ω、k-ωsst模型的计算结果.      

基于RANS-LES混合方法的翼型大迎角非定常分离流动研究

使用雷诺平均Navier-Stokes方程-大涡模拟（RANS-LES）混合方法中的延迟分离涡模拟（DDES）方法,模拟了NACA 0015翼型在大迎角下的静态绕流和强迫振荡运动并与实验值进行了比较。在大迎角静态翼型大分离流动模拟中,DDES方法捕获了非定常RANS计算未能获得的机翼背风面的涡脱落现象。在所采用的RANS和DDES模型中,基于剪切应力输运（SST）湍流模型的SST-DDES混合方法给出的时均压力系数分布与实验吻合得最好。在大迎角强迫振荡翼型绕流模拟中,DDES方法得到的非定常气动载荷与实验值吻合得很好,正确地反映了最大迎角处阻力和俯仰力矩的阶跃性突变,而非定常RANS计算则给出了完全错误的趋势。     

跨音速非定常流动的Euler方程隐式解

在随时间变化的贴本坐标系中给出求解非定常Euler方程的连续通量分裂法。在此基础上建立了可用于跨音速非定常流动的Euler方程隐式求解法。采用特征向量变换,可在保证原方程组离散化精度的条件下使计算大为简化。针对振动翼绕流特点建立了固连于物体的动坐标与固定坐标间的关系。数值计算在动坐标中进行,既简化网格生成又保证在物面上满足边晃条件。对NACA64A-10冀型绕1/4弦点做简谐俯仰振动的非定常气动力进行了计算,给出了与实验结果基本相符的计算结果。此外,还给出翼型做沉浮及同时进行沉浮与俯仰二自由度振动的非定常气动力的计算结果。     

基于非结构网格的气动弹性数值方法研究

基于非结构网格和网格变形技术,研究了耦合求解流体和弹性结构体互相作用问题的数值方法.气动弹性模型采用了耦合求解可压缩雷诺平均的Navier-Stokes方程和线性结构动力学方程的方法.为了能处理在复杂几何域内的问题,采用了基于混合单元的非结构网格有限体积方法.开发的方法在二维翼型简谐振动进行了验证并应用于AGARD 445.6翼型颤振边界的数值模拟.      

二维机翼二自由度跨音速颤振分析

 <正> On the assumption that airfoils take small amplitude ha-rmonic oscillating and middle reduced frequency, the two-dimensional, inviscid, small-disturbance nonlinear unsteady transonic flow equation can be split into two parts. One is the nonlinear steady small-disturbance equation, the other is the time-linearized unsteady small-disturbance equation. The steady equation is solved by Carlson code. The unsteady equation is solved by integral method. The standard v-g method is used to solve the flutter eigenvalue equations. A transonic flutter analysis is performed for a NACA 64A006 airfoil with pitching and plunging degrees of freedom.The Mach numbers considered are 0.7, 0.8, and 0.85. The aerodynamic coefficients are obtained by the above method. For each Mach number, the flutter speed and the corresponding flutter reduced frequencies are achieved by varying the airfoil-airmass ratio, plunge-pitch frequency ratio. The results are compared with other computational method and good agreement has been observe     

翼型舵面偏转非定常流动数值模拟

针对NACA0012翼型舵面偏转问题,数值模拟了不同参数对翼型气动特性的影响。基于非结构动网格技术,采用ALE有限体积描述下的二维可压缩非定常N-S方程,计算通量采用Vanleer格式、时空二阶格式,利用Venkatakrishnan限制器抑制数值振荡。非定常计算结果表明,NACA0012翼型绕1/4弦点作周期性俯仰振动的升力系数和俯仰力矩系数结果与实验数据吻合良好,验证了数值方法的准确性;在翼型舵面表面有分离区产生,升力系数和俯仰力矩系数形成滞回环,在亚声速情况下,滞回环幅值较小,进入超声速阶段以后,幅值增大,随着翼型间缝隙宽度逐渐增加,翼型升力系数和俯仰力矩系数与无缝翼型相比逐渐降低。     

基于联合射流的翼型动态失速流动控制研究

动态失速控制对于提高翼型气动特性具有重要意义。采用联合射流方法对翼型俯仰动态失速控制进行数值模拟,完成两方面的研究:一是射流关闭时射流通道对动态失速特性的影响,二是射流开启时不同射流动量系数对动态失速控制的影响和分析。结果表明:射流关闭时,射流通道的存在对翼型上仰过程中附着流阶段的气动特性产生不利影响,使得升力系数明显下降,但是对翼型下俯过程中失速分离流阶段的气动特性影响不明显;射流开启后,动态失速特性得到极大改善,迟滞环面积显著减小,升力增加,阻力减小,且阻力和力矩的峰值显著减小,原基准翼型力矩曲线的负阻尼区域消失。     

二维非定常欧拉方程的有限元解法

本文通过求解时间相关的欧拉方程，对振荡翼型的非定常二维无粘可压缩绕流作了数值模拟。求解方程时采用了有限体积积分的方法，并用ＶａｎＬｅｅｒ逆风通量分裂格式来计算控制体通量。隐式格式的采用，大大提高了计算效率；同时，通过插值使得空间离散达到了二阶精度。另外，本文采用基于非结构三角形网格的动态网格来处理刚性物体的移动，从而顺利地解决了复杂几何外形物体的非定常绕流的问题。文章最后给出了ＮＡＣＡ００１２翼型绕四分之一弦线作简谐俯仰振动的计算结果。一个周期内的瞬时压强分布以及升力和力矩系数都与实验结果吻合良好，     

翼型前缘变形对动态失速效应影响的数值计算

翼型或机翼的动态失速效应所引起的低头力矩和正气动阻尼限制了飞行器气动性能的提高,甚至可能诱导发生不稳定运动。应用于小尺寸机翼的前缘动态变形（DDLE）技术,通过实时改变前缘形状,能够改善翼型前缘区域的速度梯度,进而抑制动态失速效应。采用转捩剪切应力输运（SST）黏性模型结合分区混合动态网格技术,研究了这种前缘变形对机翼俯仰运动所引起的非定常流动的影响,得到通过小幅度前缘变形抑制和延迟动态失速的方法,从而提高翼型的气动性能。翼型NAC A0012的数值模拟结果与动态失速风洞试验结果比较表明：所使用的数值计算方法能够较为准确地模拟翼型在动态失速过程中升力系数与俯仰力矩系数的变化情况,可用于研究前缘变形对翼型俯仰运动所引起的非定常流动的影响。前缘动态变形翼型俯仰运动过程的非定常流场的数值模拟表明：在大迎角下不同幅度的前缘下垂运动能够抑制流动分离的发生,从而抑制动态失速,但在大迎角下小幅度高频率的前缘下垂变形能更高效地抑制动态失速;前缘变形幅度以及变形沿中弧线的分布对升力系数和俯仰力矩系数的影响并不明显。     

基于天线安装的菱形翼无人机翼型优化设计

以菱形翼布局高空长航时(HALE)传感器无人机研究为背景,分析了翼型内部安装平面相控阵雷达天线的基本原理。通过推导雷达探测性能的估算方程,建立了在翼型中安装平面相控阵雷达天线的数学模型。采用数值模拟的方法,对典型雷诺数下的RAE2822翼型的气动特性进行了研究并与试验数据进行了对比,验证了使用有限体积法k-ωSST(Shear Stress Transport)湍流模型求解雷诺平均Navier-Stokes方程在这一状态下的可靠性和适用性;同时对不同厚度的NACA 64A系列翼型的流场结构和流动机理进行了分析,在此基础上提出了一种满足平面相控阵雷达天线安装情况下翼型的优化设计思想。优化结果验证了这一设计思想的可行性。     

Stall flutter prediction based on multi-layer GRU neural network

The modeling of dynamic stall aerodynamics is essential to stall flutter, due to the flow separation in a large-amplitude pitching oscillation process. A newly neural network based Reduced Order Model (ROM) framework for predicting the aerodynamic forces of an airfoil undergoing large-amplitude pitching oscillation at various velocities is presented in this work. First, the dynamic stall aerodynamics is calculated by solving RANS equations and the transitional SST-γ model. Afterwards, the stall flutter bifurcation behavior is calculated by the above CFD solver coupled with structural dynamic equation. The critical flutter speed and limit-cycle oscillation amplitudes are consistent with those obtained by experiments. A newly multi-layer Gated Recurrent Unit (GRU) neural network based ROM is constructed to accelerate the calculation of aerodynamic forces. The training and validation process are carried out upon the unsteady aerodynamic data obtained by the proposed CFD method. The well-trained ROM is then coupled with the structure equation at a specific velocity, the Limit-Cycle Oscillation (LCO) of stall flutter under this flow condition is predicted precisely and more quickly. In order to predict both the critical flutter velocity and LCO amplitudes after bifurcation at different velocities, a new ROM with GRU neural network considering the variation of flow velocities is developed. The stall flutter results predicted by ROM agree well with the CFD ones at different velocities. Finally, a brief sensitivity analysis of two structural parameters of ROM is carried out. It infers the potential of the presented modeling method to depict the nonlinearity of dynamic stall and stall flutter phenomenon.     

扑动翼型的低雷诺数气动特性分析

通过求解引入拟压缩项的不可压Navier-Stokes方程,数值模拟了绕扑动翼型的低雷诺数非定常流动。针对厚度在4%-12%之间的NACA对称翼型,分析了翼型厚度等参数对扑动翼型气动特性的影响。在低雷诺数条件下,对于纯俯仰运动,随着翼型厚度的减小,平均阻力系数也变小。而对于纯沉浮运动,发现翼型厚度对气动特性的影响和俯仰运动有很大的差别,平均阻力系数随着翼型厚度的减小而变大。通过对沉浮运动一个周期流线图的分析,认为这是翼型前缘涡的影响造成的。由于前缘涡的影响,翼型厚度增加,平均压差阻力系数变小,甚至会出现负值。雷诺数的影响研究表明,随着雷诺数的增加,扑动翼型的阻力系数减小的趋势越缓慢。