型腔充填过程中气液两相流计算及三维自由表面的确定

针对液体充填型腔过程的特点,在固定计算域中数值模拟型腔充填过程中气液不相混两相流流场。采用了体积守恒形式的连续性方程,以防止由于气液两相密度差大而引起的数值求解中的不收敛现象。将PHOENICS软件的流场计算与守恒标量求解大变形自由表面的方法相结合,在固定计算域中完成了数值模拟三维液态充填型腔的过程。针对两种情况下的计算结果表明,本文提出的方法可用于高速充填过程和复杂型腔充填过程的数值模拟。     

飞机外形网格拓扑和多块网格拼接技术

多块网格方法成功的关键在于网格拓扑的选择和块与决之间的拼接方法。针对真实飞机外形多块网格生成,本文分析研究了翼类和体类部件的各种网格拓扑;论述了翼身组合体的多块网格非结构连续拼接技术;指出连续拼接方法的缺点,以及相应的解决方法;同时,还讨论了分块策略和确保块与块之间网格线完全连续拼接的方法;通过实例进一步阐述了本文工作的另外一些独到之处。     

关于两方程湍流模型的考虑

通过对两方程湍流模型的分析,提出了在主流区或远壁区用两方程模型,近壁区用代数湍流模型的思路,即所谓的“双层”模型概念。通过对平板湍流附面层的计算表明,这种分层计算措施获得的速度分布、壁面剪切力较单纯应用低雷诺数两方程湍流模型时的结果准确。对非均匀来流条件下90°弯管内高湍流度流场的计算表明,双层模型可获得更加准确的速度分布。从而表明本文提出的分层计算措施是有效的,可满足工程计算要求。      

两端固支输流管的稳定性和混沌运动分析

本文详细讨论了两端固支输流管在其支承简谐运动激励下的稳定性及混沌运动,推导了管道运动方程,运用相平面图,Lyapunov指数,Poincare映射以及运动时间历程等分析手段重点考察了系统参数u(流体流速）和n(激振频率）变化对管道运动的影响,结果表明,在所研究的系统中存在发生混沌运动的参数区域,而且发现有3条通向混沌运动的途径,即：（1）周期倍化;（2）间歇阵发混沌;（3）概周期运动。     

挤压油膜阻尼器-滑动轴承-转子系统非线性动力特性的数值分析及实验研究

 采用非线性模型,研究挤压油膜阻尼器 -滑动轴承 -转子系统的运动稳定性和分岔特性。利用有限元离散变分不等原理和基于线性规划余原理的迭代方法,求解自由边值条件下的油膜力。采用模态综合技术中的缩减系统线性自由度的方法,并提出基于 Newmark法的预估 -校正 -局部迭代的方法,求解系统响应。采用 Poincaré映射法,将周期解的确定转化为 Poincaré截面上不动点的确定,然后分析不动点的稳定性。最后通过实验转子,从实验与数值计算进行对比分析,验证了挤压油膜阻尼器的减振特性及其提高系统运动稳定性的特性,同时也说明了上述方法可以有效地分析这类系统的局部非线性动力特性。     

块结构的MUSCL算法求解三维进气道流场

本文将ＭＵＳＣＬ格式扩展为三维算法,并用于湍流流场的数值模拟。通过求解超声双斜楔压缩面流动,得到了含有激波、嵌入激波、剪切流和湍流边界层等流动特征的复杂流动结构,数值解与实验数据吻合得很好。文中将此算法与块结构网格技术相结合,使之能求解具有三维复杂几何形状的流场。并应用该算法模拟了Ｆ－１６战斗机进气道在５种典型飞行状态下的三维湍流流场,计算结果清楚地揭示了这些情况下的三维流动现象,得出了进气道的性     

一种同位网格上的迎风有限元方法

通过Ｇｒｅｅｎ定理将对流项变量从微分算子中分离出来,从插值函数入手引入迎风格式,是对强对流问题有限元计算中对流项变量的一种新的处理方法。按这种方法采用局部斜迎风格式及速－压同位网格公式,构成了一种对高Ｒｅｙｎｏｌｄｓ数流体流动数值模拟比较有效的有限元方法。     

用双流体—轨道模型模拟四角喷燃模型炉内三维湍流两相流动和煤粉燃烧

目前炉内两相流动和煤粉燃烧数值模拟中多半用颗粒随机轨道模型和单流体无滑移模型，这些模型都难以完整地给出三维空间内颗粒速度，浓度，湍流度分布的信息。主采用双流体－轨道模型（颗粒相连续介质－轨道模型）对一个四角喷燃模型炉内三维湍流两相流动及煤粉燃烧进行了模拟。此模型基于欧拉气相方程组、欧拉颗粒连续方程组和动量方程组以及拉氏颗粒能量方程和质量变化的方程，并使用k-ε-kp两相湍流模型，EBU－Arrhenius湍流燃烧模型，离散坐标辐射传热模型，煤粉颗粒的水分蒸发，热解挥发模型和焦炭燃烧的扩散－动力模型等。热态模拟中，为了减小为信散造成的影响，采用了扭转坐标法（将坐标扭转一定的角度使之与煤粉射流方程一致）。为了检验数值模拟，采用三维相位移普勒测速仪（PDPA）对于冷态模型炉内湍流两相流场进行了测量，得到了两相速度，湍流脉动及颗粒浓度的分布。分别对冷态模炉内两相流动和热态模型炉内三维两相流动和煤粉进行了模拟，冷态两相流动的计算与实验结果的对比表明预报的两相流场是合理的，热态模拟的结果给两相速度，气相温度、组分浓度及壁面热，显示出靠近出口处气相速度和温度分布不对称，造成一个局部高温区。     

俯仰阻尼导数分量的CFD数值模拟

提出一种直接求解直接阻尼导数的方法,该方法不仅适用于轴对称外形,也适用于非轴对称外形.数值模拟飞行器的非定常强迫沉浮运动和强迫角振动,并在Etkin的非定常气动力模型基础上,辨识得到飞行器的洗流时差导数和俯仰阻尼导数,研究飞行器俯仰通道各阻尼导数的数值计算方法.对弹道外形和基本带翼导弹标模外形及Hyflex升力体外形进行研究分析表明,将阻尼导数分量相加得到的俯仰阻尼导数与直接求解强迫角振动得到的俯仰阻尼导数与试验结果吻合很好,各阻尼导数分量随质心位置的变化趋势也与理论预测相符.对于带翼飞行器,超声速条件下,洗流时差导数在俯仰阻尼导数中占主导作用.      

基于CFD的斜盘/滑靴副油膜特性分析

结合三维Navier-Stokes方程和任意拉格朗日-欧拉(ALE)描述方法,首先,提出了一种基于计算流体力学(CFD)的滑靴副油膜特性分析方法,该方法能综合考虑滑靴副结构参数、柱塞泵工况参数对油膜特性的影响。然后,针对某种滑靴副结构,仿真得到了24种不同工况下的油膜厚度,分析了工况(温度、转速和出口压力)与滑靴副油膜的定量关系。本文还提出了一个与油膜特性相关的液动力参数的描述公式,并基于某固定结构尺寸的滑靴副CFD仿真结果,研究证明了对于固定结构的滑靴副结构,该参数仅与油液黏度(温度)相关,与其他工况参数无关。根据液动力参数公式,可以方便地给出滑靴副油膜厚度的解析方法。最后,将基于解析方法与基于CFD仿真方法得到的油膜厚度结果进行了对比,证明了解析方法的准确性。