长尾喷管两相流流场计算

为了给长尾喷管热防护结构的烧蚀和温度场计算提供准确的热边界条件，从二维轴对称N－S方程出发，采用颗粒轨道模型和有限体积的Jameson格式计算了长层喷管纯气相流场和两相流流场，分析和比较表明：在长尾喷管喉部附近以及扩散段的相同轴向位置，两相流情况下轴线上的燃气马赫数小于纯气相情况下的燃气温度；在进行长尾喷管热防护结构设计时应以两相流流场作为热边界条件。     

导弹水下点火的流体动力研究

对轴对称导弹水下点火这一非线性瞬态流动提出以下数学模型及求解方法。水流场假设为理想不可压势流,用边界积分方程法求解;喷管中的燃气流动按准一维非定常无粘完全气体流动处理,用逆步有限差分格式的特征线法求解;喷入水中的燃气流采用简化的等压泡模型,用Euler-Lagrange方法模拟气泡的演化。按照时间步进方法实现了水、气流动和物体运动的数值耦合求解。数值实验结果反映出导弹水下点火时水气流动的一些重要特性。     

绕近地面二维平板非定常不可压缩位势流的精确解

Yih给出了极端地效条件下二维水平平板垂直下落时有势绕流的首阶近似解;给出了同一问题在平板同时还有水平速度情况下的近似解。本文利用椭圆函数将平板和地面围成的区域保角变换至矩形内部,然后在变换平面上求解,得到了有地面效应时水平平板铅垂运动所诱导的势流精确解,以及平板同时还有水平速度情况下的“准定常”连续绕流近似解。     

再入飞行器高超声速绕流的一步和二步有限元计算

本文针对再入飞行器的空气动力问题，采用一步和二步有限元格式求解可缩守恒型Ｅｕｌｅｒ方程，并结合不同的人工粘性模型来捕捉激波，给出了高超声速二维钝头绕流的结果。     

非结构网格高超声速化学非平衡MHD流场数值模拟

 在非结构网格上对考虑化学非平衡效应的二维高超声速磁流体绕钝头体流动进行了数值模拟。控制方程由二维理想磁流体动力学（MHD）方程和组元连续方程两部分组成,化学动力学模型为5组元17反应模型。MHD方程空间离散采用AUSM格式,时间推进采用显式5步龙格-库塔格式,并通过弱耦合的方式与化学反应控制方程结合在一起。计算模型为二维钝头体,外加磁场为偶极子场,磁场源位于钝头体内部。在高超声速来流条件下,对有、无磁场干扰,是否考虑化学反应下的4种工况进行了数值计算,得到了满意的结果,并与有限的参考文献进行了对比。结果表明本文发展的方法能准确地模拟考虑化学非平衡效应的高超声速MHD流场。     

串装双燃速药柱发动机的内流场计算

对一维的加质燃烧室流场、轴对称的二维喷管内流场和高燃速药柱燃烧完成后的燃烧室二维内流场以及三维药柱几何做了一体化计算。应用轴对称的N－S方程、TVD有限体积法的显式MacCormack格式，对燃烧室后段和长尾喷管跨声速流场进行求解，得到了双燃速内弹道性能和轴对称二维流场中燃气参数的分布。结果具有较高的精度，并为进一步的燃烧室和喷管热结构分析提供数值依据。     

用N—S方程计算翼型非定常粘性大攻角绕流

本文给出了用ＮＳ方程对翼型非定常粘性大攻角绕流的数值模拟。控制方程为二维时均可压缩完全ＮＳ方程；湍流模型采用双层代数涡粘性模型。使用近似因式分解ＡＤＩ差分格式离散求解，网格是用保角变换方法生成的相对翼型固定的Ｃ型网格。本文给出两类典型非定常绕流数值模拟结果：翼型过失速常攻角周期流动和大攻角强迫俯仰谐振非定常绕流。并与国外的实验和计算结果进行了比较，表明了本方法准确、高效的特点。     

二维非定常欧拉方程的有限元解法

本文通过求解时间相关的欧拉方程，对振荡翼型的非定常二维无粘可压缩绕流作了数值模拟。求解方程时采用了有限体积积分的方法，并用ＶａｎＬｅｅｒ逆风通量分裂格式来计算控制体通量。隐式格式的采用，大大提高了计算效率；同时，通过插值使得空间离散达到了二阶精度。另外，本文采用基于非结构三角形网格的动态网格来处理刚性物体的移动，从而顺利地解决了复杂几何外形物体的非定常绕流的问题。文章最后给出了ＮＡＣＡ００１２翼型绕四分之一弦线作简谐俯仰振动的计算结果。一个周期内的瞬时压强分布以及升力和力矩系数都与实验结果吻合良好，     

高超音速化学非平衡钝劈绕流数值计算

 自1970年Davis提出粘性激波层方法以来,用数值方法求解高超音速轴对称钝体绕流问题国内外已做了大量工作,但至今仍未见到关于平面问题的计算结果。对于平面问题,虽然方程在形式上比轴对称简单,但由于二维效应,激波层较厚,用文献[2]的方法向下游区推进有困难。另外,在驻点线上采用极限关系式虽能克服方程的奇性,但驻点解对流向步长Δξ有依赖。     

弦向缝隙叶栅对边界层分离的控制

本文提出了一种在Ｓ１流面上进行势流与粘性流迭代求解弦向缝隙叶栅流场的计算方法，将Ｓ１流面上势流方程组的解与经Ｉ１１ｉｎｇｗｏｒｔｈ－Ｓｔｅｗａｒｔｓｏｎ变换后的边界层方程结合起来联立求解弦向缝隙叶栅流场。计算值同实验结果相当吻合，在正大攻角时，向缝隙叶栅能效地控制边界层分离。     

叶轮机械新流型探索雏议

从“流型”的角度, 对比飞行器空气动力学和叶轮机械气体动力学的发展历程, 认为叶轮机械内的流动仍在追求定常附体流型, 由于其流动及结构特性, 使得“有害”分离问题难以避免。本文提出应探索将外流领域内的定常脱体涡流型或下一代非定常脱体涡流型引入到叶轮机械的可能性与技术途径本文并综述了上述探索性研究的有关背景和现状, 以及所设想的基本途径和影响因素。      

平板扰流片振动诱导涡运动研究

 本文研究了平板上扰流片静止及振动时所诱导的分离涡的形成过程。当拢流片静止时。在拢流片后仅形成不规则的大尺度涡,一旦执流片振动。就诱导出有固定频率的分离涡。拢流片振动频率,对分离涡的尺度、强弱、成对、合并过程有重要影响。文中还对的流片振幅大小的影响进行了讨论。     

超音速绕凹角湍流分离流动的计算

本文用统一的Levy-Lees变换以及正算法与逆算法相结合,求解了超音速绕凹角湍流分离流动。 对附着流区用边界层正算法,压强分布用流过尖劈统一的高超音速与超音速公式,湍流模型取代数涡粘性模型;对凹角分离区用边界层逆算法,给定位移厚度δ~*分布,湍流模型取代数松弛模型;边界层计算采用Cebeci-Keller Box方法;计算成功地算得分离流场,较好地预估了分离点与重附点位置以及壁面压强分布与表面摩擦应力分布。     

跨音速叶栅粘流计算的多级广义有限元法

将 Taylor-Galerkin广义有限元法和多级有限元的思想结合起来 ,并在人工粘性的处理上作了改进 ,构成了收敛速度和稳定性均较好的多级广义有限元法。利用该方法 ,并基于 N-S方程研究了二维跨音速叶栅流动。计算结果与实验结果符合较好。通过计算表明 ,该方法计算稳定、收敛较快 ,是叶轮机械内部跨音速流动强有力的计算手段     

凹半球降落伞模型在加速和定常流中的流场观测

 在北京航空学院水槽中,用氢泡法对绕凹半球模型流动进行了研究。模型前后流速分布是用氢泡时间线测量的,它与在起动开始时由无旋流理论估算值接近,但当起动旋涡变大以后,在尾流中流速差别变大。在模型后部诱导出较大的反流速席。所以由无旋估计的虚质量只在开始时是可用的,并且在起动涡变大时,虚质量应更大。 在定常流中,在模型前有时形成一个大旋涡,这个旋涡是非定常的并且可能引起侧向力和伞的非定常运动。模型上适当开孔可减小甚至消除此旋涡。     

文氏效应在叶片根部马蹄涡控制中的应用

为研究文氏管在叶片根部马蹄涡控制作用的应用,针对简化了的叶片-平板结构在叶片上游平板上的某位置开狭缝,在平板的另一侧构造一文氏管(文氏管的喉道位于狭缝处),引射叶片根部的流体,进而控制马蹄涡.数值模拟了不同狭缝位置的湍流流动;结果表明了改变狭缝相对于叶片的位置能很好地控制马蹄涡,甚至能完全消除叶片前马蹄涡.提出的方法可以利用外势流的能量而无需额外消耗能量,为一种被动的控制手段,有很大的实际应用前景.利用文氏效应的控制方法简单实用,因此工程中采用该法是不错的一种选择.      

LTRAN2的一种改进方案

本文提出了LTRAN2的一种改进方案BTRAN2用来分析翼型的各种频率的非定常运动。用E-O调转换格式和ADI方法求解了完全的二维非定常跨音速小扰动位流方程,并用单调转换的AF2格式计算了定常跨音速小扰动方程,以此作为非定常计算的初场。本文给出了带后缘正弦振荡襟翼的NACA 64 A006翼型的绕流和做正弦俯仰振荡的NACA 64 A010翼型的绕流计算结果,它们与Euler方程解或实验数据很吻合。     

利用机上传感器测量燃油流量的方法研究

主要介绍了在小型商用飞机无法加装流量传感器的情况下,利用机上传感器实现流量测量的方法,通过该方法完成了该飞机在各种状态下燃油流量的测量。     

非定常N-S流中的瞬时临界点

本文用微分方程的相似理论把非定常轨线方程变换成自洽系统,考察变换后临界点邻域在非定常流场中的映象,进而对非定常流中“瞬时临界点”溉念的可用性及其限度作了讨论。     

小振辐振荡圆柱诱导的二次旋涡流动

本文给出充满粘性不可压缩流体时同心圆柱间内柱作小振幅振荡所诱导的二次旋涡流动的解析解,对典型的柱径比给出容器外边界对二次旋涡流动流线、涡量分布和动能分布的影响。给出了等涡量线和等动能线分布的鞍点。采用数值求解非定常Stokes流动考察了二次旋涡流动的发展过程。二次定常旋涡流动的显示实验表明,解析解、数值计算和实验显示所示流谱基本一致。