低速大迎角三角翼旋涡流动的有粘怀无粘数值模拟

通过人工压缩性方法求解了三维定常不可压Ｅｕｌｅｒ方程和雷诺平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，应用Ｂｅａｍ－Ｗａｒｍｉｎｇ近似因子分解格式及其对角化形式，采用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型。计算模型为７０°三角翼，分析了有粘与无粘流动旋涡特性的影响，尤其是对Ｅｕｌｅｒ方程模拟大迎角分离涡、涡破裂的能力及其局限性提供了数值验证，计算结果与实验结果相一致。     

旋成体大迎角分离流的Euler方程数值模拟

用三维Euler方程计算旋成体的大迎角绕流,我们发现:亚声速自由流时,计算得到的是无分离流动,与实际的分离流不相符合,因而计算的法向力大大小于有分离的实验值;而当自由流为超声速,且横向自由流M数较高时,由于背风表面出现激波,无粘Euler方程的数值计算能自动捕捉到背风面上的分离涡,计算法向力与实验值接近。为了对亚声速自由流用Euler方程数值模拟大迎角旋成体背风面上分离涡,根据问题的物理图像,在假设的分离线上,提出强加Kutta条件的方法,成功地计算出有分离流动,使计算的法向力与实验值接近。     

螺旋桨流场Euler/NS数值解对比分析

采用多块点对点网格生成技术,生成螺旋桨流场空间计算网格。针对不同前进比,分别求解Euler和NS方程对螺旋桨流场进行数值模拟,并与实验结果进行了对比分析。研究发现无粘和有粘的计算结果均接近于实验结果,因此基于Euler方程的流场数值模拟更适用于工程设计。     

悬停旋翼无粘流场数值模拟中的多重网格方法

发展了一种基于不同空间离散格式的多重网格算法,并应用于悬停旋翼无粘绕流的Euler方程数值模拟。由于悬停旋翼流场中存在不可压区域,同时旋翼尾涡系统的发展需要较长的时间,使得旋翼流场的收敛速度远低于固定翼流场,因此研究旋翼流场的多重网格算法具有重要意义。空间离散采用了Roe s FDS格式和Jameson中心有限体积格式,时间推进应用了五步Runge-Kutta方法。采用多重网格V循环方式,对一跨声速悬停旋翼无粘流场进行了数值计算。计算结果表明:多重网格算法可以显著加速悬停旋翼无粘流场的数值计算收敛速度;无论在激波分辨率还是在计算精度上,Roe s FDS格式都优于JST格式。     

曲边单元间断有限元方法求解二维Euler方程

考虑到间断有限元方法对边界的敏感性,采用基于八节点曲边四边形单元的间断有限元方法求解了Euler方程的圆柱绕流问题。详细推导了八节点四边形单元的变换关系,给出了Jacobi矩阵行列式的具体表达式。对比直边单元和曲边单元的计算结果,采用曲边单元后,计算结果符合Euler方程的无粘假设,得出了八节点四边形单元间断有限元方法求解Euler方程是合适的结论。     

在无粘跨音速流动中影响机翼颤振的参数研究

随着基于气动Euler方程的人工无粘计算流体动力学谐波平衡代码的应用，一些结构参数和自由流跨音速马赫数具有二自由度典型跨音速翼型的颤振，所用谐波平衡气动模型和时间有关项采用线性离散，这样处理以后，和以前计算模型的效率相比，可以获得更高的计算效率。     

弯头钝锥的超音速无粘绕流数值计算方法

本文用差分法给出了弯头钝锥的超音速无粘绕流值计算方法。为解决物体体轴弯曲给推进计算带来的困难,本文通过Euler方程自变量的转换,将方程的求解面变为可以同推进方向倾斜的平面,即对Euler方程进行“斜推”或“平推”求解。由数值算例表明,本方法能得到较好的结果。     

跨声速翼型绕流的Euler／边界层方程干扰数值解

本文利用Euler方程和可压缩湍流边界层积分方程研究绕跨声速翼型的有粘与无粘强干扰流动。应用有限差分法在贴体的网格上求解时间相关的Euler方程,以剪功积分方法求解翼面贴附和分离湍流边界层流动,并引入一个松弛方程描述剪应力对上游湍流历程的延迟响应。有粘/无粘干扰采用表面源模型。计算结果表明,对翼面存在强干扰流动情况,获得了与实验值基本吻合的结果。     

应用重迭网格技术求解复杂组合体无粘流场

本文数值求解了捆绑火箭这一复杂组合体的无粘干扰流场。该计算利用重迭网格(Overlapping Grid)技术求解了分离系数的Euler方程,获得了复杂组合体相互干扰的重要流场特性。通过与风洞试验的比较,计算结果与实验有较好的一致性。     

不同马赫数的无粘和粘性流动高阶精度隐式计算方法

 应用隐式时间推进法对不同马赫数的无粘和粘性流动进行数值分析,给出了基于预处理方法的高阶精度隐式求解方法。方程离散采用改进的高阶精度对流迎风分裂格式。该方法通过求解曲线坐标系可压缩Euler和 Navier-Stokes方程,对低速到超音速范围内的无粘和粘性流动的典型问题进行了数值分析。计算结果与文献的数值结果或实验数据对比分析表明,该方法对低速到超音速范围内的无粘和粘性流动问题进行数值分析是可行而有效的。     

基于控制理论的Euler方程翼型减阻优化设计

在Jameson的控制论气动优化设计思想下，应用Euler方程研究了翼型的反设计和减阻问题。设计过程中的梯度是通过求解一个伴随偏微分方程而得到的，每个设计迭代只需求解一次Euler方程和一次伴随方,一与设计变量数无关。在具体数值执行中：①利用分部积分将目标函数变分的计算转化成了类似于计算Euler方程残值的形式，节省了机器时间；②根据控制理论的要求，将减阻问题转化成了相应地修改伴随方程的物面边界条件和目标函数的变分表达，使问题得到了简化；③利用特征不变量分析法，处理了伴承方程的物面和远场边界条件。设计算例证明了本文方法可靠性好、收敛快、大大节约设计时间。     

基于遗传算法多段翼型外形优化设计

为提高多段翼型增升效能,开展襟翼外形和缝道参数同时优化设计研究。优化算法采用遗传算法,以求解RANS方程为气动特性分析方法,通过椭圆方程控制生成多段翼型外形,同时优化缝宽、搭接量、襟翼偏角等位置参数和襟翼外形控制参数,实现多段翼型优化设计。设计实践表明,与只优化位置参数相比,同时优化襟翼外形和位置参数得到的多段翼型有更大的升力系数,方法是可行的,具有一定的工程应用前景。     

基于控制理论的压气机叶型数值优化方法

将基于控制理论的气动优化方法应用于轴流压气机叶型设计.以Euler方程作为流动控制方程,具体推导得出了其相应的伴随方程,分析了边界条件,并给出求解方法.以给定压力分布作为目标函数,将参数化叶型作为设计变量,在求得目标函数对设计变量的梯度信息后,结合BFGS优化算法得到优化方向,更新设计变量完成叶型的优化设计.通过三个算例验证了该叶型优化设计方法的有效性.      

气动声学计算中一种声扰动方程的改进

研究了气动声学计算中的一种改进声扰动方程(IAPE)。采用三种不同的声传播模型,数值分析了声在均匀流以及剪切流中的传播问题。与线化欧拉方程(LEE)相比,除了个别位置预测的压力峰值较小以外,声扰动方程(APE)几乎能够预测到类似的物理特性;另外,APE方程具有无需求解密度方程,计算成本相对较少等优点。为弥补APE方程预测压力峰值较小的缺点,通过分析APE方程,这类方程中需要加入剪切流与声波相互作用的源项。为此,对APE方程进行改进,加入了一种声源项,计算结果表明,改进的APE方程(IAPE)能够和线化欧拉方程的计算结果保持一致。     

用椭圆型、抛物型两种方法计算平面叶栅湍流尾迹

本文用连续方程、稳定不可压缩椭圆型和抛物型时均N-S方程及K-ε双方程湍流模型计算了平面叶栅湍流尾迹。计算表明,上述两种方法都能得到较好的数值结果,但抛物方法可节省大量计算时间和内存空间。     

有二次涡的激波边界层干扰流动的DPNS和NS方程计算

本文针对具有二次涡复杂分离再附现象的激波边界层干扰流动,数值地考察了扩散抛物化Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ（ＤＰＮＳ）方程组的适用情况。壁面摩阻和压力、主涡和二次涡的涡高和涡长、分离再附位置以及流线图等特性的计算表明：ＤＰＮＳ方程组的数值结果均与ＮＳ方程组的数值结果很好相符。     

边界层转换阶段直接数值模拟的影响矩阵法

提出一种求解时间相关、三维不可压完全Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的影响矩阵法。考虑流动区域为上半空间的平极边界层流动,在水平方向,利用周期性条件,采用Ｆｏｕｒｉｅｒ离散,在半无穷方向,定义了一种特殊映像并利用Ｃｈｅｂｙｓｈｅｖ多项式展开。由此将Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程转化为以水平波数（ｋ１,ｋ２）为参数的一组法向Ｈｅｌｍｈｏｌｔｚ方程组,同时可更有效地求解压力。本文讨论了连续问题的影响矩阵法与原问题的等价性。影响矩阵法的离散形式,以及影响矩阵法的算法实现等。最后给出了一维Ｈｅｌｍｈｏｌｔｙ问题的数值实验结果。     

钝体俯仰阻尼导数数值计算

分别采用：（1）Euler方程的谐振摄动法；（2）薄层近似的非定常N－S方程模拟有粘扰动流场两种方法计算了钝体外形的高超声速俯仰阻尼导数。通过计算表明：对于复杂外形，方法（1）适用于零攻角附近的流动计算，而方法（2）考虑了流场的粘性作用，能够应用于一般流态（3＜M＜20，0＜a＜20°）的计算。此外，本文发展了文献［6］的动导数转换公式，使之适用于二维转轴情况，并在方法（2）中通过引入“亚迭代”过程提高了非定常流场的计算精度。     

半波GLM方法封闭的湍流方程对转捩过程的数值模拟

ＴＲＡＮＳＩＴＩＯＮＳＩＭＵＬＡＴＩＯＮＵＳＩＮＧＧＬＭ－ＥＭＢＡＳＥＤＣＬＯＳＥＤＥＱＵＡＴＩＯＮＳＯＦＴＵＲＢＵＬＥＮＣＥＧａｏＧｅ（４ｔｈＤｅｐｔ．ＢｅｉｊｉｎｇＵｎｉｖｅｒｓｉｔｙｏｆＡｅｒｏｎａｕｔｉｃｓａｎｄＡｓｔｒｏｎａｕｔｉｃｓ，Ｂｅ...     

解决欧拉方程奇异性的方法探讨

提出了一种解决欧拉方程奇异性的工程算法,从飞行器运动方程入手,推导出具体计算公式,并进行了数值计算验证。该方法无论是理论分析还是工程计算都证明是可行的,没有发现原理缺陷和方法误差,不会产生奇异性,在飞行器仿真、模拟和姿态控制中具有普遍的应用和参考价值。