可压缩性粘流绕翼型的摩阻计算

本文研究了一种可压缩性粘流绕翼型的摩擦阻力计算方法。这种方法包括用修正了的Thwaites~[1]法计算头部层流边界层以及按Green~[4]法计算湍流边界层。按照上述方法计算了在不同雷诺数和不同马赫数下两个翼型剖面的摩阻系数。经与文献[5]比较,吻合较好。     

有弓形冲波影响情况下的层流附面层方程的数值解

本文给出了一种轴对称层流附面层方程的数值解法。利用Levy-Lees变换,将物理平面上的附面层方程变换到新的坐标平面内,再利用三点隐式有限差分法来解这组方程。在考虑熵吞影响的情况下,本文发展了一种与国内外惯用的流线吞嚥法不同的不需迭代即可给出附面层外缘参数的方法,从而使计算时间缩短了三分之二。本方法给出的结果同其他作者的数值解和实验结果作了比较,结果是令人满意的。本文还给出了驻点热流在等熵和变熵情况下随雷诺数和马赫数变化的计算结果,并同实验结果作了比较。     

双基推进剂热分解反应动力学参数的计算

本文根据在CDR-1型热分析仪上测得的一条非等温DSC曲线,确定了双基推进剂热分解反应的最可几机理函数,提出了计算三个动力学参数(E,A和n)的数值方法.该法用线性最小二乘法、迭代法以及二分法与最小二乘法相结合的方法,以积分方程、微分方程和放热速率方程拟合实验数据.在逻辑选择建立了微分形式或积分形式的机理函数的最可几一般表达式后,用放热速率方程计算相应的动力学参数E,A和n的值.     

用N-S方程有限体积法计算弹体绕流/底喷流

用单元中心有限体积法在物理空间曲线网格离散积分形式N-S方程，用坐标变换计算粘性通量项中的各偏导数值，用Runge-Kutta三步显式格式进行时间推进求解。采用了人工粘性、当地时间步长及隐式残值光顺。计算了旋成体大迎角涡流和带底部喷流的绕流场。计算中发现，各向异性隐式残值光顺较之各向同性隐式残值光顺能更快地加速迭代过程的收敛。计算的物面压力分布与实验基本符合，计算的旋成体大迎角涡流图谱合理。计算结果说明底部喷流对旋成体绕流及空气动力载荷有显著影响。     

双时间法求解大迎角翼型绕流非定常N-S方程

本文采用双时间推进法、顶点格式有限体积法求解非定常ＮＳ方程、应用多层网格技术计算了ＮＡＣＡ００１２翼型大迎角振荡绕流,并对伪时间推进收敛精度,物理时间步长,多层网格迭代方式等对计算结果与计算ＣＰＵ时间的影响进行了讨论,另外还讨论了采用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ湍流模型时用固定转捩或全湍流模式对计算结果的影响,并与实验结果进行了比较。     

一种二元高超声速进气道起动特性的尺度效应研究

对一种二元混压式进气道三维流场进行了数值和实验考察,研究了不同尺度进气道模型自起动性能的变化。结果表明,在相同来流单位雷诺数条件下,随着模型尺度的减小,进气道自起动马赫数有所提高,起动性能有所降低。同时对不同尺度模型进行雷诺数匹配,发现在相同雷诺数下,不同尺度模型的起动性能相近,表明雷诺数是影响不同缩尺模型起动性能不同的主要原因。在可获得的实验结果范围内,数值模拟所得到的自起动结果基本与之相符。此外,对实验中发现在低雷诺数下进气道反而呈现出自起动特征的异常现象进行了初步分析,通过数值模拟比较指出了低雷诺数下来流偏向层流流态,可能会导致进气道呈现一种“起动”状态。     

直升机桨叶扬起下坠碰撞动响应计算

针对铰接式旋翼桨叶扬起下坠碰撞问题,采用中等变形梁模型,同时引入有限转角模拟桨叶绕铰的刚性运动。运用Hamilton原理建立桨叶扬起下坠碰撞动力学方程,采用隐式Newmark法计算桨叶在有初始挥舞角下的下坠碰撞响应。计算结果与试验数据的对比表明:如果采用定步长的方法计算响应,计算时间较长,计算结果会产生明显的延迟。为了解决上述问题,在碰撞前后设置了一碰撞缓冲区,采用了变步长的方法强制计算点在跳过此缓冲区时必须进入该区。此法在计算桨叶扬起下坠碰撞动响应时,能较好的解决计算延迟问题,计算时间大大减少,计算结果与试验数据吻合较好。     

基于γ-Re_(θt)转捩模型的低雷诺数翼型数值分析

通过综合运用Michel转捩判据和γ-Re_(θt)转捩模型的方法实现对低雷诺数翼型的气动分析。首先通过Michel转捩判据和层流分离位置分别估算转捩动量厚度雷诺数,然后分别用Langtry和Tomac提出的经验关系式计算转捩动量厚度雷诺数和转捩区长度,通过UDF(User-Defined Functions)将不同组的3个经验关联值写入Fluent软件对E387翼型(雷诺数为3×10~5)进行数值分析,并与试验值进行了对比。结果表明:基于Michel转捩判据比基于层流分离更能准确预测转捩位置,基于层流分离预测的转捩位置过于靠前,但是2种方法的气动力计算结果与试验值都比较吻合,Langtry和Tomac的经验关系式计算结果之间差别不大,Tomac方法能够计算出分离泡,且与Selig油流试验吻合较好。     

基于变步长CFD/RBD方法的旋转弹轨迹仿真

采用变步长方法实现非定常Navier-Stokes方程和刚体六自由度运动方程的耦合求解,并对美国ARL旋转弹试验模型的飞行轨迹进行了仿真。旋转弹自由飞行过程中受到空气阻力作用,转速下降,气动力变化趋缓。因此轨迹仿真时先选取较小的时间步长,随着转速的下降逐渐增大步长,这样可以减少计算时间。刚体六自由度运动方程的数值求解采用变步长Adams预估校正法。将定步长和变步长两种方法应用于旋转弹实验模型轨迹仿真,并将计算得到的旋转弹空间位置、姿态角与试验结果进行了对比,结果显示:方法能较精确地仿真旋转弹飞行轨迹。在相同初始步长条件下,变步长方法较定步长方法计算时间减半;在相同计算时间条件下,变步长方法较定步长方法精度更高。     

翼型低雷诺数层流分离泡数值研究

翼型在不可压和低雷诺数条件下工作时的性能越来越引起了人们的研究兴趣。本文采用Rogers发展的双时间步Roe格式,求解拟压缩性修正不可压N-S方程。数值模拟了低雷诺数条件下(Re=6.0×104,1.0×105,2.0×106),E387翼型在不同攻角条件下(α=0°、4°、7°)翼上表面后缘部分的层流分离现象,对应于长分离泡。时均化的结果同McGhee[12]的试验结果进行比较,验证了计算结果的正确性。研究了低雷诺数层流分离的非定常和时均化特性,并对漩涡脱落中主涡、二次涡以及出现的漩涡对并的周期性过程和对气动力脉动造成的影响进行了较为细致的研究。非定常计算结果表明,低雷诺数条件下的层流分离现象,是周期性的旋涡脱落过程。所谓长层流分离泡是其时均化积分的结果。     

螺旋形强制谐波振荡对涡破碎位置的影响

本文用非线性波理论分析并计算了螺旋形振荡对涡破碎位置的影响,并在旋涡发生器中进行了施加不同频率的振荡对涡破碎位置影响的实验。实验与计算结果的趋势基本一致。     

机翼大迎角涡破裂时气动特性计算方法的研究

本文计算了小展弦比机翼纵向和横侧气动特性,并计及涡破裂的影响。本方法以Purvis升力面理论和Polhamus吸力比拟为基础。并参照Lan方法进行大迎角特性修正。算例表明,本方法适用于小展弦比机翼气动特性的计算。     

近壁湍流相干结构的数值模拟

提出了近壁湍流对称单个相干结构的理论模型,采用紧致有限差分和Fourier谱展开相结合的数值方法研究它的演化问题,得到雷诺应力分布及高剪切层的形成过程与实验结果一致。     

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

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

圆柱突然起动的快速涡解法

 文中将Chorin的涡方法、格子涡(VIC)方法和快速Fourier变换(FFT)相结合,发展了一套快速有效的求解高Re数下非定常N-S方程的算法(快速涡方法)。通过对二维圆柱突然起动问题的数值计算,得到流动分离及其演化过程(包括二次涡及二次涡对的生成);此外还得到回流区对称轴线上速度分布。这些结果与实验数据及其它数值解相比,是令人满意的。该程序对其它高Re数非定常流动问题也是适用的。     

喷管扩张段烧蚀的“台阶效应”计算

当喉衬和扩张段内衬采用具有不同烧蚀率的材料时,常出现“台阶效应”.本文以气动热化学烧蚀机理为基楚,提出应用适当的对流换热和加剧的化学反应计算对台阶效应区进行修正的方法,使台阶效应区的烧蚀率预示更切合实际工况.     

可压混合层涡并放热效应的涡动力分析

 基于双元单步不可逆燃烧动力学简化模型 ,采用紧致差分算法直接求解二维可压 Navier-Stokes方程 ,研究了时间发展非预混扩散火焰平面反应混合层流动。拟序结构由数值线性稳定性分析方法设定的粘性可压最不稳定小扰动初场促生。随放热强度增大 ,得到了基频涡卷撕裂现象。涡动力分析表明 ,放热主要通过膨胀和斜压这两种强非线性过程抑制大尺度剪切涡的演化     

某涡轮级三维流场非定常数值研究

基于耦合Lantry-Menter转捩模型的SST(Shear stress transport)湍流模型,对某涡轮级在3个雷诺数工况下的全三维流场分别进行了定常、非定常数值模拟,目的在于对比分析两种计算方法所得结果之间的差异。计算结果表明,进口雷诺数为4×104时和5×105时,定常与非定常计算的涡轮级效率的差异较大,最大相差0.6个百分点,而当雷诺数为8×104时,两者差异很小;非定常计算能够较好的模拟叶尖泄漏涡的发展过程、低雷诺数下上游尾迹与下游转子叶片吸力面分离边界层干涉诱导卷起涡的形成、输运及其所导致的叶片吸力面表面压力大幅波动等非定常流动现象。最后,给出了转子出口截面上周向与径向的非定常压力脉动的分布,为基于Lighthill声类比方法的计算气动声学研究奠定了基础。     

求解非定常N-S方程的一种自适应时间步长的罚函数有限元方法

 本文提出了一种求解速度-压力变量形式的非定常不可压N-S方程的数值方法。其中采用罚函数的Galerkin有限元进行空间离散,改进了梯形法的时间积分格式并结合自适应时间步长技术使计算效率最优。文中对圆柱突然启动到Re=100和200时的流场进行了数值模拟,所得到的Karman涡街的诱发过程是满意的。     

雷诺数效应对小流量多级轴流压气机的性能影响

通过试验和数值模拟,对一台设计压比为6,效率0.84的小流量轴流压气机进行了雷诺数效应的分析和研究。介绍了这台压气机在设计转速下的性能曲线、不同雷诺数下折合转速n=0.85的试验结果,以及二维平面叶栅和三维数值模拟分析结果。试验和数值模拟的分析表明,随着雷诺数降低,会增加附面层厚度,改变压气机的三维流场结构,降低压气机的性能。根据折合转速n=0.85的试验和数值模拟分析结果,对这台压气机在折合转速n=1.0性能进行了雷诺数修正。