三维曲线坐标系N—S方程通用形式和数字解

运用张量分析理论，比较简捷地导出了三维任意曲线坐系上椭圆型Ｎ－Ｓ方程的通用形式。在对这些方程进行数值离散时，采用了有限控制体积的非交错风格设置，该算法双ＳＩＭＰＬＥ为基础，为了抑制压力振荡场的了出现，在压力修正方程中对压力梯度的离散附加一项１－δｘ的差分。作为算例，本文给出了三维９０°强曲率弯曲管道中气流的不可压缩紊流状态的数值计算结果，并与相关的实验数据进行了对比，结果是满意的。     

二维突扩管流的有限分析解

运用有限分析法求解Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，获得了二维突扩管流的数值解。对于雷诺数为２５，５０，１００，１５０，２００的来流，预报的流动图谱具有相似性，回流区域及再附距离的计算也获得了比较满意的结果     

任意迎角二维平板的超声速和高超声速俯仰振动导数

应用频率的摄动法,本文得到了任意迎角二维平板机翼的超声速和高超声速俯仰振动导数。该导数是由包括直到缩减频率的五阶项所组成的。数值结果表明高阶项的贡献在较低马赫数或大迎角情况下是很重要的。作为平均迎角趋于零的特殊情况,本文的数值结果与超声速线化的位流理论是一致的,在高超声速情况与活塞理论一致。     

用Euler方程计算某些特殊的二解分离流动

采用非定常Ｅｕｌｅｒ方程、总体时间步长，隐式时间推进，计算了几种特殊类型本的二维有分离绕流。一类是光滑的圆柱体，在跨音速时计算出因激波诱导分离而周期性地脱落在尾迹上的非对称涡排。另一类是带尖角的物体，如三角形柱体，计算出从尖角处脱落在尾迹上的分离旋涡，形成卡门涡街，对这两类物体，计算的平均阻力系数，分离涡脱落的Ｓｔｒｏｕｈａｌ数以及流态都与实验接近。结果表明，对于这类绕流问题，用非定常Ｅｕｌｅｒ方     

双重互易边界元法在求解N—S方程中的应用

提出了用双互易边界元法求解Ｎ－Ｓ方程的新思路。用Ｐｅａｃｅｍａｎ－Ｒａｃｈｆｏｒｄ算子分理解将时间相依的Ｎ－Ｓ方程分解为线性和非线性子问题，线性问题用共轭梯度法解除压强－速度的耦合；非线性问题进行局部线化。对所得的Ｐａｓｓｉｏｎ方程及相关类型方程采用双重互易的边界元法求解，消除了传统Ｎ－Ｓ方程边界元解法的区域积分问题。     

非结构网格的高超声速粘性MHD流场数值模拟

应用数值方法在非结构网格上对磁场干扰下的二维高超声速钝头体粘性绕流进行了数值模拟。控制方程为N-S方程耦合Maxwell方程的粘性MHD方程组,空间离散采用有限体积方法,对流项用AUSM格式计算,粘性项用中心格式求解,时间推进用显式5步Runge-Kutta格式,引用双曲型散度清除技术加强▽.B=0的条件。计算模型为二维钝头体,在高超声速来流条件下,分别对有、无均匀磁场干扰下的流动进行了数值模拟。计算结果表明,在均匀磁场干扰下,激波脱体距离显著增加,物体表面压力急剧下降。对比表明文中发展的计算方法可以准确地进行二维粘性MHD流场的数值模拟。     

螺旋桨滑流对飞机绕流影响的试验研究

在国产双涡桨支线飞机发展过程中，曾专门研究过螺旋桨滑流对飞机表面压力分布的影响。本文以此为基础，从滑流对飞机绕流影响出发，研究了螺旋桨滑流引起飞机气动特性的变化，重点讨论了气动特性变化的机理。     

快响应PSP技术用于高超声速圆柱绕流的非定常压力测量

压敏涂料(PSP)技术是飞行器风洞实验大面积定量测压和流动显示的重要工具。为了发展非定常流动压力测量和脉冲风洞中全局压力测量能力,航天空气动力技术研究院开发了快速响应压敏涂料技术,与中科院化学所共同合作开发的压敏涂料采用PtTFPP作为发光基团,稳定性较强,持续光照下发光强度衰减为1.5%/h。采用自主研制的静态标定设备在标定腔内测量不同温度和压力条件下20mm×20mm涂料样片的表面发光强度来确定涂料的压力响应特性和Stern-Volmer公式,并设计制作了2套快响应时间动态标定设备,测得涂料典型响应时间在0.2ms。在中国航天空气动力技术研究院FD-03高超声速风洞中对平板圆柱装置进行了Ma=5的 PSP 试验,利用高速相机采集图像数据后,经过批量数据处理,采集频率为250Hz,得到了间隔时间4ms 的连续压力场数据。结合纹影图像和油流图对得到的PSP结果进行了分析,利用同时采集的测压孔数据对PSP 结果进行了比较。试验结果表明,快响应PSP技术可以更详细地显示流场结构,并能同时得到很好的定量压力数据。     

三圆柱绕流的时均压力分布和气动力

亚临界状态下绕等边三角形排列的三圆柱的流动,其中心距的变化范围为1.5～4.0d(d为圆柱直径)。测量了不同间距下各圆柱的压力分布并推算了气动力。结果表明,在小间距时,三圆柱间的干扰是严重的,圆柱表面的压力分布有相当大的改变,局部地区有很强的负压;除阻力之外,还产生了相当大的横向力。     

一种快速求解Orr—Sommerfeld方程中性稳定性曲线的方法

提出一种快速准确求解Ｏｒｒ－Ｓｏｍｍｅｒｆｅｌｄ（Ｏ－Ｓ）方程的中性曲线和临界雷诺数ＲｅＣｒ的方法。在求解过程中,先分析α０,ｃｒ０的范围,利用网格自适应自动搜索α０,ｃｒ０,采用三次样条函数计算并迭代到一定精度从而获得雷诺数Ｒｅ与波数α,ｃｒ的不足。算例表明,本方法计算速度快,结果准确,花较少的机时就能计算出与谱方法精度相当的临界ＲｅＣｒ。     

尖前缘任意翼型的超音速和高超音速俯仰安定性导数

本文推广了文[1]的欧拉方程摄动法,讨论了尖前缘任意翼型俯仰安定性导数的计算方法,并应用该方法计算了双圆弧翼型的俯仰安定性导数,其数值结果与文[2]的超音速线化位流方程有限差分法的数值结果比较,更接近于实验值。在激波附体的条件下,本文的方法可用于计算大迎角尖前缘任意翼型的超音速和高超音速俯仰安定性导数。     

混合网格在求解三维N—S方程中的应用

用三棱柱／四面体组成的混合网格求解三维Navier-Stokes方程，在物面附近采用三陵柱网格，而其他地方则采用四面体网格，与单一类型网格相比，混合网格易于调节粘性层内网点数目，并能减少存储量，提高解的效率，Navier-Stokes方程的求解采用有限体积法，对人工耗散项进行了改进，使得求解的结果更接近真实的物理解，本文对上述方法开展了数值实验研究。结果表明，用三棱柱／四面体组成的混合网格求解N－S方程是非常有效的。     

后掠机翼欧拉方程数值模拟

本文对后掠机翼用保角转绘加剪切变换生成C-H型网格,并在这类网格拓扑上,用Jameson的三维欧拉方程有限体积法,四步Runge-Kutta时间推进格式,研制出可供分析三维后掠机翼亚、跨、超声速绕流的计算程序。本方法的特色是改进了机翼表面网格点的分布,使机翼后缘与网格线一致。在数值计算中,通过附面层位移厚度计算,引入了附面层粘性修正,改善了计算结果。     

超声速流中激波/湍流附面层干扰数值模拟

采用修正的B／L湍流模型以及多块结构化网格求解了二维N-S方程。分别对超声速流和高超声速流中的激波／湍流附面层干扰进行了数值研究。本文首先研究了进口马赫数为2．96的超声速流。计算结果准确预测了入射斜激波在平直壁面引起湍流附面层分离的流动特征：分离点的反射激波、分离包引起的膨胀扇以及再附点的反射激波。计算的壁面压力分布与实验值吻合较好，计算的分离区长度与实验值比较有一定误差。本文还对进口马赫数为9．22的高超声速流中压缩角引起的激波／湍流附面层干扰进行了数值研究。计算结果与实验结果吻合较好。     

三维曲线坐标系N－S方程通用形式和数字解

运用张量分析理论，比较简捷地导出了三维任意曲线坐标系上椭圆型Ｎ－Ｓ方程的通用形式。在对这些方程进行数值离散时，采用了有限控制体积的非交错网格设置，该算法以ＳＩＭＰＬＥ为基础，为了抑制压力振荡场的出现，在压力修正方程中对压力梯度的离散附加一项１—δｘ的差分。作为算例，本文给出了三维９０°强曲率弯曲管道中气流的不可压缩紊流状态的数值计算结果，并与相关的实验数据进行了对比，结果是满意的。     

热完全气体的热力学特性及其N-S方程的求解

首先用五次多项式拟合给出了温度在50～3000K范围内的热完全空气的焓值与温度之间的函数关系式，导出了其它热力参数e，cp，cv和γ的表达式。接着提出了热完全空气总温、总压的计算方法，并将其计算结果与量热完全空气的结果进行了比较。最后，将该热完全空气模型用于N-S方程求解，对NAPA软件进行了改进，并用该软件时高马赫数钝体绕流流场和乘波体流场进行了计算，分析了气体真实效应时流场结构及参数的影响。结果表明，本文提出的总温、总压计算方法及时NAPA软件的改进是成功的，可较准确地模拟高超声速流动的主要特征。     

串列叶栅粘性紊流流场有限分析数值研究

本文将近年发展起来的有限分析方法(Finite Analytic Method)用于曲线坐标系上紊流N-S方程的数值计算。计算了来流雷诺数为2.0×10~5,气流攻角分别为0°,10°,-10°三种情况的二维粘性紊流串列叶栅流场。文中用k—ε紊流模型模化紊流,以壁面函数方法处理近壁区流动参数。数值计算结果与实验结果相比较,吻合程度令人满意。     

辐射及化学非平衡流耦合场计算方法

把一种新的三维热辐射计算方法应用到高超声速非平衡流场的数值模拟过程中:在非结构网格上对高超声速化学非平衡流场的数值模拟过程中,耦合了辐射输运方程来求解辐射热流。控制方程为含有化学反应源项和辐射源项的三维N-S方程,数值离散格式采用了Jam eson有限体积法,辐射输运方程采用有限体积法求解。化学反应模型为7组元7反应模型,并采用了化学反应特征时间步长。本文分别对RAM C-Ⅱ和MU SES-C模型进行了计算,计算结果与参考文献的结果吻合。     

增强能量增量法求解无轴承开关磁阻电机电感

有别于磁链法和需要求解二阶偏导数的能量增量法,本文采用一种称作增强能量增量法计算非线性系统的磁链、静态电感(割线电感)和动态电感(差分电感).该法首先给定各回路的工作点,其次在工作点附近给予一较小的扰动电流,则系统的磁场也随之扰动.结合磁性材料的B-H曲线,磁场能量增量可以通过对磁场强度H、磁感应强度B增量的积分获得.由此,在不需要微分的条件下,各回路的磁链、静态电感和动态电感可以分别求得.本文以12/8 7.5 kW无轴承开关磁阻电机为例,在认为三相绕组完全分时导通时,可忽略相间互感;以线性铁芯材料为例,说明了主、悬浮绕组间互感在合适的连接下也可以忽略.因此,仅需求解一相主绕组、悬浮绕组的静态、动态自感即可.求解这些电感随位置角和工作电流的变化,取得了与实验一致的结果.