有限平板绕流Monte—Carlo方法仿真

本文采用直接仿真Monte-Carlo方法求解有限长平板绕流等问题。本方法是通过计算机跟踪仿真分子的运动来实现数值模拟的。仿真分子间的碰撞计算由统计抽样确定。碰撞模型分别选用了硬球分子模型和负幂律分子模型。仿真分子与固壁作用采用由完全扩散反射和镜面反射按比例混合组成的模型。为了检验方法的可靠性,还计算了激波结构和Rayleigh问题等一维流动,二维计算采用同步并行程序。数值结果表明直接仿真Monte-Carlo方法能够较好地模拟稀薄气体力学中的一些问题。对于二维计算,所花费的机时和所需的内存均在国内机器所允许的范围内。     

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

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

轴流压气机子午流面计算中三维剪切作用的湍流模型

本文发展了一个轴流压气机子午流面二维计算中模拟三维剪切结构对子午流场作用的湍流模型。试验与这一模型进行了初步的比较。模型表明,在多级轴流压气机子午流场计算中,必须计入实际流场三维剪切结构对子午流场的湍流扩散作用。     

轮毂处理扩稳的机理探讨

在单级轴流压气机上,用周向槽式处理机匣在静子内径处进行轮毂处理,与转子尖部机匣处理作类比。用小型5孔探针测试技术详细测量了在最佳工作状态和近失速边界状态下静子叶片排下游气流的三元流场,特别是其端部流场,发现对于转子叶尖失速型的级,在其静子轮毂处理也可以提高级的稳定裕度。并提出其扩稳机理与失速转移的关系,提出了周向诱导速度影响二次流的物理模型,进一步剖析了周向槽机匣处理的扩稳机理。     

等速上仰翼型绕流结构的数值模拟

通过数值求解涡量-流函数形式的N-S方程研究等速上仰翼型的分离流动结构。用ADI方法解涡量方程;用Poisson方程直接法解流函数方程。在计算得到的流场中除可看到实验中观察到的多涡结构外,还发现一些实验中未观察到的现象。将分离流动结构的计算结果与翼型气动力变化相联系,并结合涡动力学理论,定性地解释了上仰翼型产生高升力的原因。     

转动坐标系中三维跨声速欧拉流的有限体积 TVD格式

在非惯性转动坐标系中,本文采用贴体网格、有限体积法离散和修正数值通量技术,将Harten的一维TVD格式推广到三维。由于转动使方程出现源项,文中通过对源项的巧妙处理,使修改后的格式能用于非齐次双曲守恒律方程组高分辨率的数值计算。为了加速解的收敛,提高显式时间推进的CFL数,本文采用隐式残值光顺技术。三维跨声速带非齐次源项欧拉方程的典型算例表明:捕捉的激波分辨率较高;激波前后没有发现大的数值波动和伪振荡现象;所得的跨声速流场解与实验较接近。     

冷热两股气流在圆管内掺混的速度场、温度场与辐射通量场的数值方法

用数值方法求解冷热两股轴对称平行射流在圆管内掺混的流场、温度场和辐射通量场。计算采用一种通用性的数值方法,其特点是强烈地以物理本质为基础,而不仅仅是单纯的数学推导。这样做,能方便地分析和解释计算结果,使之更真实地接近于实际。 通过计算,分析了当引射系数改变时,掺混的流场,温度场和辐射通量场的变化规律,其结果与实验数据比较一致。     

新型高超音速静压探针的设计和实验

参考了近40年来国外有关超音速和高超音速静压探针的研究成果,设计了多种并发展了二种高超音速静压探针。它们可精确地用于持续工作流场仅为30ms的流场静压测量。实测静压与自由流静压比p_(s,m)/p_(s,∞)=1.009;马赫数比M_(∞m)/M_(∞,)=0.993。实验在加拿大多伦多大学宇航研究院的高超音速炮风洞里完成。炮风洞喷管出口气流特性是:(空气)自由流总压比_(t,∞)=25.5MPa,总温T_(t,∞)=1000K;马赫数M_∞=8.30;自由流静压p_(s,∞)=2.13kPa;雷诺数R_e=3.2×10~7。以探针直径D为基准的R_(e,D)=34,000～52,800。这种探针己被用于“管内高超音速流场发展实验研究”中了。     

跨音速非定常流动的Euler方程隐式解

在随时间变化的贴本坐标系中给出求解非定常Euler方程的连续通量分裂法。在此基础上建立了可用于跨音速非定常流动的Euler方程隐式求解法。采用特征向量变换,可在保证原方程组离散化精度的条件下使计算大为简化。针对振动翼绕流特点建立了固连于物体的动坐标与固定坐标间的关系。数值计算在动坐标中进行,既简化网格生成又保证在物面上满足边晃条件。对NACA64A-10冀型绕1/4弦点做简谐俯仰振动的非定常气动力进行了计算,给出了与实验结果基本相符的计算结果。此外,还给出翼型做沉浮及同时进行沉浮与俯仰二自由度振动的非定常气动力的计算结果。     

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

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