应用滑移界面方法求解动／静相干 叶排内非定常流动

基于滑移界面技术，发展了求解二维动／静相干叶排内非定常粘性流动的计算方法和软件，该方法克服了混合平面方法所引入的方法误差．压力与速度耦合采用非迭代PISO算法；计算区域的离散采用非交错斜交网格．计算了二维不可压缩Couette流动及无限长相对运动动／静叶排内的非定常流动，验证了用该方法和软件计算动／静相干叶排内非定常流动的有效性．     

鱼类游动的N—S方程数值模拟

通过求解二维守恒型非定常流动N-S方程，对鱼类运动时鱼体周围二维非定常流场进行数值模拟，数值计算空间采用有限体积离散，时间推进格式采用Jameson的双时间法。描述了各状态鱼体模型、计算网格的生成方法和数值计算方法。结合计算结果对鱼类游动时反卡门涡街的形成和鱼体推力的产生机理进行了分析，并对不同斯特哈尔数对推力的影响进行了研究。结果表明：非定常N-S方程数值方法能够正确模拟鱼体运动时的流场和反卡门涡街的形成；产生明显推力的St数的范围是0．23～0．53，而且在0．35处推力达到最大。     

求解非定常方柱绕流的可变周期时间谱方法

针对周期性非定常流动特点，采用时间谱方法求解了无限大区域内低雷诺数下二维不可压缩层流绕方柱流动问题．由于方柱绕流的非定常运动周期事先未知，在计算过程中加入了基于残差梯度的可变周期计算方法，使变量和周期可以同时迭代计算．研究了时间谱方法的计算精度和求解速度，分析了不同求解参数对最终结果的影响．计算结果表明，时间谱方法可以准确地模拟方柱后漩涡的非定常运动现象，而计算所需时间只是传统的基于有限差分的双时间推进方法的10％～25％，流体非定常运动的Strouhal数和方柱的阻力特性证实了这种方法的正确性．可变周期计算方法对求解参数的变化不敏感，可以准确地预测运动周期．     

应用滑移界面方法求解动/静相干叶排内非定常流动

基于滑移界面技术,发展了求解二维动／静相干叶排内非定常粘性流动的计算方法和软件,该方法克服了混合平面方法所引入的方法误差．压力与速度耦合采用非迭代ＰＩＳＯ算法;计算区域的离散采用非交错斜交网格．计算了二维不可压缩Ｃｏｕｅｔｅ流动及无限长相对运动动／静叶排内的非定常流动,验证了用该方法和软件计算动／静相干叶排内非定常流动的有效性     

非等距网格高精度差分方法用于气动声学问题计算

研究了非等距网格下高阶精度有限差分方法用于气动声学问题的可行性．通过Taylor级数展开法构造了不等距网格下的七点六阶精度空间离散格式，分析了格式适用的波数范围，时间积分采用显式四阶精度推进格式，可直接计算非定常欧拉方程用于气动声学问题．算例采用随机变化的网格间距，对一维单波方程和球形波方程进行了计算，验证了非等距网格格式模拟波动问题的能力．对二维情况，计算了亚音速均匀流中初始声、涡和熵脉冲波问题，得到了很好的结果．     

收一扩喷管内二维紊流的数值研究

本文应用隐式有限差分理论对喷管内二维紊流流动进行了数值求解．采用了贴体曲线坐标系统，代数紊流方程模化紊流，整个流场用时间推进迭代方法获得．研究表明，控制方程右端粘性项矢是主变量欠的一个线性函数．文中所提供的计算方法与计算结果将为三维非定常计算打下基础．     

预处理算法在非定常流动计算中的应用

介绍了一种可应用于全速度范围的非定常N-S方程数值求解方法,并全面考核该方法的适用性。该方法通过引进伪时间导数项,对伪时间导数进行预处理以解决低速时方程系数矩阵刚性过大的问题;同时用双时间步推进算法保证时间精度,将该方法推广应用到非定常问题的求解。文中给出了若干算例,覆盖了从极低马赫数到超音速、无粘/层流/湍流、二维/三维、定常/非定常情况。算例结果表明,该方法在较宽的使用范围内均能得到理想的计算结果。     

二维非定常对流扩散方程的对角占优隐格式及多重网格算法

提出了一种数值求解二维非定常变系数对流扩散方程的对角占优、空间为二阶精度的隐格式。利用Fourier分析方法证明了该格式是无条件稳定的．由于格式具有对角占优性，因此适用于大梯度（高雷诺数）问题的数值求解．另一方面，为了克服传统迭代法在求解隐格式时收敛速度慢的缺陷，采用了多重网格加速技术，大大加快了迭代收敛速度，提高了求解效率．数值实验结果证明了该格式的精确性、稳定性和对高网格雷诺数问题的强适应性．     

Navier-Stokes方程的最佳非线性谱Galerkin算法

提出了求解二维非定常ＮａｖｉｅｒＳｔｏｋｅｓ方程的最佳非线性谱Ｇａｌｅｒｋｉｎ算法,分析了近似解的一致收敛速度．和标准的谱Ｇａｌｅｒｋｉｎ算法与非线性谱Ｇａｌｅｒｋｉｎ算法相比,该算法具有节省计算量的优点．     

贯流风机气动噪声数值预估

通过精细求解二维非定常Reynolds平均的Navier-Stokes方程,数值模拟了贯流风机内部的复杂流场。随后从流场的数值结果中提取出叶片、涡墙和后墙的脉动压力作为声源,进行声场计算。以声学中的Ffowcs Williams-Hawk-ings(FW-H)方程作为出发方程,数值求解贯流风机的噪声场。计算结果表明在贯流风机中,后墙的压力脉动与涡墙的压力脉动是主要的噪声源。该文的数值预估不仅在贯流风机的总体气动性能上与实验测试结果吻合,同时气动噪声场的预估结果也与实验测试结果吻合良好。     

基于伴随方法的动态非定常气动外形优化设计

为实现动态非定常条件下的气动外形优化设计,基于非定常形式的欧拉方程及伴随方程,建立了一套完整的针对非定常动态问题的气动外形优化设计方法,并采用所建立的方法,对俯仰振荡的NACA0012翼型在跨声速(Ma=0.8)及限制条件下进行了翼型时均减阻外形优化. 对优化设计的结果进行了分析. 研究结果表明:该动态非定常气动外形优化设计方法具有较高的可靠性和实用性.      

基于神经网络的跨音速非定常气动力的辨识

利用递归神经网络(RNN)模型具有时间记忆性,且会考虑之前的输入输出对当前输出影响的特点,以递归神经网络方法建立了NACA0012翼型在跨音速阶段的非定常气动力模型;利用CFD计算NACA0012翼型绕其刚心作变频俯仰运动的跨音速气动力系数为训练数据,建立跨音速非定常气动力模型。以建立的跨音速非定常气动力模型预测NACA0012翼型作俯仰简谐振动的气动力系数,并与CFD计算的气动力系数进行对比。结果表明,该模型具备优良的逼近非线性非定常气动力的能力;针对跨音速二维翼型,该模型相比CFD可以更快速地构建,并能迅速且较为准确地预测不同频率下作简谐振动时的气动力。     

翼型颤振分析的时域方法

使用全耦合法预测了翼型NLR7301的颤振边界.在流体域的动网格上,通过求解非定常的雷诺平均N-S方程,计算出翼型表面的非定常气动力.在耦合了翼型的结构动力学方程后,得到翼型表面的变形位移,通过变形网格技术对流体域的网格进行更新.在此基础上利用时间推进的方法模拟了翼型的自激振动过程.     

风力机翼型非定常气动力分析(英文)

基于状态空间描述的Beddoes-Leishman动态失速模型,对风力机翼型进行非定常气动力分析.考虑到风力机翼型工作时的实际情况,在模型中忽略了可压缩性效应和起始于翼型前缘的流动分离.模型考虑了气流的近尾流效应和在失速区域的后缘分离效应.模型用4个气动状态来描述非定常气动力系数动力学,其中两个用于描述近尾流效应中的时间迟滞,另两个用于描述后缘分离效应.采用模型对做简谐变桨运动的FFA-W3-241风力机专用翼型进行了非定常气动力分析.数值计算结果与实验值吻合良好,说明模型能较好地描述风力机翼型的非定常气动特性.将动态失速数值计算模块与已有的风力机气动与结构分析软件集成,利用软件对一台1.5MW变速变桨距风力机的发电工况进行了仿真.仿真结果表明,翼型的非定常气动特性对动态载荷计算结果影响很大,因此在风力机的设计过程中应该予以充分考虑.     

飞行器地面效应定常/非定常流场数值模拟

采用重叠网格技术,通过求解Navier-Stokes方程,对飞行器起飞状态下增升装置的地面效应进行了数值模拟.选取GAW-1为基本翼型进行增升装置设计,进一步生成重叠网格,分别采用“虚拟运动地面”模式边界条件进行了计及地面效应定常流场数值模拟,以及飞行器运动模式进行了计及地面效应的非定常流场数值模拟.对比分析了不同数值模拟方法之间对气动特性以及流场的影响,为飞行器地面效应数值模拟提供了参考.     

基于BPOD的气动弹性降阶及其在主动控制中的应用

流体动力系统是一个高阶的非线性复杂系统,这严重限制了优化和控制在该系统中的应用。模型降阶技术是解决这一问题的有效工具。将特征正交分解方法（Proper Orthogonal Decomposition,POD）和平衡截断方法结合起来,形成平衡特征正交分解方法（Balanced Proper Orthogonal Decomposition,BPOD）,并应用到气动伺服弹性模型降阶中。该方法将由POD快照得到系统可控和可观Gramian矩阵的近似表达,通过该近似得到平衡降阶模型。以一个二维翼段的气动弹性系统为例,首先将流体控制方程线化;然后利用BPOD方法得到非定常气动力的降阶模型以及降阶的气动弹性系统;最后,对降阶系统设计主动控制律,通过控制面偏转来抑制翼型颤振。数值仿真结果表明BPOD降阶模型可以精确地模拟高阶非线性的流体动力系统并且可以有效地应用于气动弹性主动控制中。     

小幅振动旋成体跨音速绕流的数值分析

假设绕小幅振动旋成体所导致的物理量的变化为平均定常扰动量的高阶小量,那么将物理量(速度、密度及压力)摄动展开后代入Euler方程,便导出了前两阶摄动量的控制方程和相应的边界条件,用TVD格式进行数值求解,得到了平均定常跨音速流(零阶)和基频(一阶)非定常分量的数值解。     

扑翼飞行二维情形数值分析

基于计算流体力学(CFD)方法对微扑翼飞行的非定常空气动力学问题进行研究,已被大多数学者所认可.在对昆虫飞行运动仿真模拟的基础上,建立了扑翼运动二维简化模型及运动方程,并利用FLUENT软件及其嵌入用户自定义函数(UDF)和动网格技术,分析了二维翼型在不同控制飞行参数(包括扑翼频率、最大转动幅值、平动和转动相位差等)下的平均升力和阻力系数.数值模拟所得结果为微扑翼飞行器的设计提供了参考依据.     

非均匀来流下飞翼布局气动特性模拟方法研究

飞翼布局的非均匀来流会使得飞翼的流场不均匀,从而影响其气动性能。通过引入"网格速度"来计入非均匀来流的影响,求解非定常Navier-Stokes方程,实现了非均匀来流下飞翼布局气动特性的数值模拟。首先采用该方法对NACA0006翼型非均匀来流下的气动力响应进行了计算,计算结果与文献计算结果吻合良好。进一步对某飞翼布局飞机在非均匀来流下的气动力响应过程进行了数值模拟,得到了不同振荡速度比率下飞翼布局气动力系数随相位角变化曲线,结果表明,振荡速度比率越大,升力响应的幅值越大,非均匀来流对飞行器气动特性影响越明显。     

二维非定常有涡流动的数值模拟方法

为有效地模拟二维有分离现象的粘性流动,发展了求解涡量流函数方程的数值方法。对涡量输运方程的时间项运用四阶Runge-Kutta法离散,从而将方程分解为4个计算步分别求解。在每一计算步对方程进行Steger-Warming近似因式分解,从而使涡量的计算在空间的两个方向上分别进行。对流项采用Chakravaythy-Oscher总变差减少(TVD)格式离散。流函数的Poisson方程运用Tschebyscheff SLOR方法交替方向迭代求解。运用该方法对突然起动圆柱和高雷诺数时弯曲薄翼的非定常有涡流动进行了数值模拟,并将计算结果与其它方法的计算结果和试验结果进行了对比,结果表明本方法具有精度高、数值稳定性好和计算效率高的优点。