叶轮机械流场计算中流线曲率不同的处理方法及其分析

本文介绍了在叶轮机械流场计算中,用不同的方法计算流线曲率(或径向速度变化)的方法。为了便于比较,对具有较大曲率的单级压气机的进口流场的变化趋势先做了理论上的分析,接着用五种不同的方法计算了该进口流场,并与试验得出的进口流场做了比较。     

超声速流场的流线-特征线坐标变换与应用

为了求解二维平面弯曲激波波后流场,讨论并发展了一系列基于流线-特征线坐标系变换的流场代数计算方法。该系列方法依据二维平面流场中气流角、静压沿特征线近似不变的特点,可快速求解平面弯曲激波波后流场。其中,等气流角近似法适合模拟流线轨迹,等静压近似法适合求解波后流场参数。在此基础上,又提出了一种改进的等气流角-等静压混合方法用于计算弯曲激波波后流场。和特征线法对比,等气流角-等静压混合法计算得到的流场基本特征与特征线法得到的结果相同,在来流马赫数Ma=6和Ma=4情况下误差分别仅为0.5%和0.15%,证实了该方法在求解二维平面弯曲激波流场中的适用性。     

跨音速压气机转子内三维流场数值模拟

采用多步Runge-Kutta与有限差分相结合的格式,离散叶轮机械内三维欧拉方程,应用代数法生成了压气机转子三维网格。用此方法模拟了NASA Rotor37跨音速压气机转子流场,计算得到的流场与试验结果进行了对比。结果表明,该方法收敛速度快,实用可行,计算简单省时。并且,该方法为研究叶轮机械内部流场计算提供了很好的基础。在此基础上进行更深层次的研究会变得很容易。     

一种适合于旋翼前飞非定常流场计算的新型运动嵌套网格方法

针对旋翼流场计算中关键的运动嵌套网格技术,提出了一个新的网格单元关系快速判断方法。该方法通过计算特殊网格点在背景网格上的位置来建立桨叶网格在背景网格中所占据范围的关系,并用于嵌套网格洞单元的识别,又进一步将该方法推广用于Inverse Map的建立,在此基础上,建立了一套用于前飞流场计算的运动嵌套方法。通过多种前飞旋翼测试算例,分别评估了背景网格、桨叶网格尺寸改变对计算时间的影响,同时验证了该方法的可靠性,结果表明构建网格嵌套关系所用时间随桨叶网格和背景网格尺寸的变化幅度较小,且具有良好的鲁棒性,能够满足旋翼前飞状态计算的需要。最后,在基于前飞非定常N-S/Euler流场求解器中,采用该方法模拟了UH-60A和7A旋翼悬停/前飞状态的气动特性,通过计算时间的对比体现出该方法的高效性,并且流场计算结果与试验数据对比表明该方法可以有效地用于旋翼前飞非定常流场的计算。     

一种高精度气动弹性计算方法

提出了一种实用高精度静气动弹性计算方法,该方法的气动力计算基于FLUENT软件,结构变形的求解和流场网格调整在ANSYS平台下完成,结构流场网格划分使用界面节点一致化方案。气动力计算和结构变形均使用专业软件完成,可以建立精细的计算模型并获得相应的精确计算结果。界面节点一致化方案降低了流固耦合难度,气动力和结构变形的传递简单而精确。最后以三维机翼模型为算例验证了该方法的可行性有效性。     

弹体穿越平面激波三维非定常流场的数值模拟

用数值方法求解三维非定常可压缩Navier-Stokes方程，计算了复杂弹体在相同马赫数、不同迎角下的流场和气动力特征，计算结果与实验结果吻合，证明本文采用的方法能够准确计算此类复杂流场；然后采用该方法模拟了弹体穿越正激波的非定常过程，得到了弹体周围流场及气动力特征随穿越过程的非定常变化，阻力特征穿越后急剧下降。     

旋转机翼悬停气动特性研究

鸭式旋转机翼(CRW)是一种先进的高速直升机方案,旋翼同时也是机翼,是其最关键气动部件之一。采用结构分区拼接网格技术进行空间离散,分区建立参考系,通过求解多重参考系下的N-S方程来计算旋翼流场,首先以传统的Caradonna-Tung实验旋翼的亚、跨声速悬停流场分析为例验证该方法的可靠性,进而采用该方法对旋转机翼悬停流场进行了数值计算,旋翼拉力计算值和地面实验值吻合较好,结果分析表明旋转机翼的悬停流场有着不同于传统旋翼的流场特性。     

单级跨声速压气机进气畸变数值模拟

基于彻体力模型的基本思想,发展了一套能够反映进气畸变对风扇/压气机气动性能影响的三维数值计算程序CSAC.首先详细研究了彻体力模型的具体建模方法,然后发展了相应的三维数值计算程序.利用该程序针对某单级跨声速压气机在均匀进气下的流场进行了计算验证,计算结果与雷诺平均Navier-Stokes(RANS)计算吻合得很好,流场主要参数分布与相关实验数据也有很好的一致性.利用该程序分别计算分析了进口存在稳态周向总压畸变及周向总温畸变的全环三维畸变流场,结果表明:压气机进口畸变流场存在大尺度、强三维特征,将导致压气机总压比及稳定裕度下降.所有计算结果均显示该模型能够在降低对计算资源和工程经验依赖的同时,有效分析进气畸变对压气机特性的影响.      

矢通量分裂隐式有限体积法解超音速喷流强干扰流场

发展了用于轴对称流场,任意非正交曲线网格下的矢通量分裂隐式有限体积法;统一了二维和轴对称流场计算的表达式。该方法对粘性通量项进行了分解和归类,考虑了包括交叉导数项在内的所有项对隐式增量的贡献。对带中心喷流的流动计算表明该方法收敛速度和计算精度均优于ＭａｃＣｏｒｍａｃｋ显式格式。横向喷流强干扰流场的计算初步揭示了姿控发动机喷流产生间接推力的机理,由此,设计中可降低对主发动机额定推力的要求     

模型超音速燃烧室流场和性能的数值模拟

研究了简化模型超音速燃烧室流场和性能的数值模拟方法，该燃烧室在支板后缘设置一缝隙式喷嘴，平行于超音速空气流喷射氢气。用椭圆型偏微分方程数学模型，ＭａｃＣｏｒｍａｃｋ差分格式成功地计算了两种不同后绿尺寸的流场及其性能，为了对照比较，还计算了另外两种流场。计算结果表明，支板后的流场存在回流区，因喷氢的压力高于超音速空气流的压力和壁面附面层的影响，燃烧室内将出现斜冲波和膨胀波，压力沿横向变化明显，由此证明流场与“边界层流动”有性质上区别。提出了全流场按性质不同分段进行数值模拟的方法，利用计算得到的流场节点状态参数（如温度、压力等），积分获得燃烧室的性能参数，和文献的实验数据对比，计算结果合理、可靠。     

高超音速平板前缘干扰数值解

本文采用了文献〔1〕中提出的一个数值求解Navier-Stokes方程的空间二阶精度单步差分格式,数值计算了超音速和高超音速平板前缘干扰问题。通过简单的模型方程对格式进行了描述和分析。用这一方法计算了M_∞=20,Tw/T_0=0.06,Re_(∞L)=4×10~3和M_∞=3,Tw=T_0,Re_(∞L)=10~3的平板前缘干扰问题。计算结果与实验值比较符合的很好,通过这些计算结果可使我们更清楚地了解平板前缘干扰的流动特性。     

无波动,无自由参数的耗散差分格式（NND）在喷流计算中的应用

本文用文献[1]建立的无波动、无自由参数的耗散差分格式(NND格式),对无粘喷流流动进行了数值模拟,并就一典型的自由射流流动,与国外高精度、高分辨率的PLM差分格式(Piecewise-Linear Method)的计算结果和实验结果作了比较,得到了很一致的结果,且所用的计算时间少于PLM格式。     

变结构系统的递阶切换模式

 本文研究关于线性多变量变结构系统的递阶切换模式及其切换函数c与变结构控制u（x,s）的设计问题。文中提出了两种递阶一切换模式，并按模式的到达条件及对滑动模态的极点配置的原则给出了各种情况下切换函数及变结构控制的设计方法。此种模式的最大优点是系统的控制容许失误。     

再人体复杂外形无粘流场数值计算研究

本文在文[1,2,3]的基础上研究了非定常Euler方程的推进迭代方法,并将无波动、无自由参数的耗散差分格式(下称NND格式)发展为隐式、迎风格式,用以计算有差曲控制翼的再入复杂飞行器的育攻角、有侧滑角的超声速无粘绕流流场。其结果是准确的。本文的方法具有以下特点: 1.可以在同一计算程序中实现亚跨超流场的计算:在超声速区域,采用推进技术;在亚跨声速区,则自动选择推进迭代求解技术。 2.NND格式不仅具有公式简单,无自由参数的特点,而且过激波时满足熵条件,不会产生非物理解,也不发生振荡。 3.隐式差分格式的无条件稳定和Gauss-Seidel迭代的快速收敛,计算省时。 本文的方法曾用于烧蚀头部凹陷外形和航天飞机简化外形的计算,能自动选择迭代和推进区域。可以预言,如计算机条件许可,本方法可以用于大攻角的无粘流场计算。     

跨音速大迎角Euler方程数值分析

应用有限体积法离散三维Euler方程,选用三步显式格式进行时间推进求解;并通过当地时间步长、残值光顺、焓修正等方法加速收敛。对一维波动方程的von Neumann稳定性分析说明这种三步格式的最大Courant数是2。对三角翼跨音速大迎角流动的数值分析表明文中的三步格式及焓修正方法优于人们经常采用的四步、五步格式以及焓阻尼技术。     

流-势函数法的人工密度格式在叶型设计中的应用

本文采用流-势函数法进行跨音速叶型设计,用人工密度格式取代类型相关格式。与相关格式相比,人工密度格式具有全场算法统一、收敛过程稳定等优点。本文提供的设计算例表明,该方法能有效地进行叶型的初始工程设计。     

影响液体火箭发动机涡轮效率的某些因素分析

分析了液体火箭发动机中所用涡轮的各种损失因素,并建立了相应的数学模型,给出了计算方法。对于喷嘴中的二相流损失,采用龙格-库塔法求解一维定常全耦合的二相流控制方程。用迭代的方法进行冲击式涡轮及有小反力度的涡轮的设计计算,以使设计出的涡轮参数与效率相统一。通过对计算路径的控制,使得本方法和程序能同时用于冲击式涡轮及有小反力度的涡轮的设计计算;在计算结果中,还着重对二相流损失对涡轮效率的影响进行了分析。     

用CSCM格式求解激波反射问题

本文对Euler方程采用CSCM格式,数值求解了激波反射问题,与目前流行的交替方向隐式(ADI)格式以及总变差递减(TVD)格式的计算结果进行了比较,并作了简要分析讨论。     

紧致格式数值模拟超音速粘性绕流问题

通过数值求解可压Navier-Stokes方程的方法计算了球锥粘性绕流问题。方程中粘性项按通常办法逼近,为改善收敛速度利用了文[1,2]中之算子附加修正方法。方程中之流向导数利用了四阶精度的紧致差分。这是由于该方向上流动参量光滑和边界条件易于处理的缘故。在法向方向上利用了中心差分。对差分方程之隐式部分利用了近似因式分解法。由于选用了文[3]中之特殊Jacobiall系数矩阵分裂法而使计算工作量大为减少。用这种方法计算了超音速粘性绕流问题。计算结果与实验结果相符很好。     

参数不确定空间机械臂系统的增广自适应控制

讨论了载体位置与姿态均不受控制的漂浮基两杆空间机械臂系统的控制问题,为此对系统的运动学、动力学作了分析。结果表明 :结合系统动量守恒及动量矩守恒关系得到的系统广义Jacobi关系、以及系统的动力学方程将为系统惯性参数的非线性函数。证明了借助于增广变量法可以将增广广义Jacobi矩阵以及系统动力学方程表示为一组适当选择的惯性参数的线性函数。在此基础上,给出了系统参数未知时,空间机械臂末端抓手跟踪惯性空间期望轨迹的增广自适应控制方案。通过仿真运算,证实了方法的有效性。