低速圆管流动的粒子仿真

运用信息保存法对低速圆管的流动现象进行了模拟,并与实验结果进行了比较,粒子仿真结果与实验结果吻合较好,且优于NS方程的结果。研究表明,在对低速圆管的模拟过程中,运用IP法在获得较好的结果的同时,具有比DSMC方法更高的计算效率。IP算法是解决低速圆管流动问题的有效途径。     

预处理方法在全速度流场数值模拟中的应用

基于预处理方法,采用有限差分方法,LU-SGS隐式求解了二维可压Navier-Stokes方程,数值模拟了无粘突包、粘性空腔、化学非平衡流场,并与试验结果进行了对比。计算结果表明,所采用的数值模拟方法能够求解低速及亚、跨、超全速度流场和化学非平衡流场,所采用的预处理方法在亚声速范围内保持数值精度的前提下,加速了收敛。     

飞船再入舱三维化学非平衡流数值模拟

对飞船再入舱高超声速化学非平衡三维流场进行数值模拟。控制方程为含化学反应源项的全Ｎａｖｉｅｒ－ Ｓｔｏｋｅｓ 方程组, 化学模型为高温空气的七组元模型, 其组元成分为Ｎ２ 、Ｏ２、ＮＯ、Ｎ、Ｏ、ＮＯ＋ 、ｅ－ 。差分格式采用张涵信院士提出的ＮＮＤ激波捕捉格式并用全耦合方法处理流动方程及化学反应方程, 对化学反应源项采用全隐式处理, 运用时间预处理技术以加快收敛。同时, 也对完全气体模型进行计算, 以分析真实气体效应。与实验数据对比, 说明本文结果是可信的。     

涡流阀流量调节特性数值模拟

用数值模拟研究了涡流阀的流量调节特性。在来流总压恒定的条件下,用时间相关法求解了雷诺平均的Ｎａｖｉｅｒ—Ｓｔｏｋｅｓ方程。计算了无控制流和有控制流两种情况,得到了流动参数在流场中的分布和控制流量与喷管流量的关系。计算表明,涡流阀有良好的流量调节特性。     

二维水翼空化流动的数值模拟

基于RANS方程,结合不同湍流模型与空化模型对粘性流中二维NACA66012模型水翼附近流场进行数值计算与模拟.计算得到固定攻角、固定来流速度时升力系数和阻力系数随空化数的变化规律以及固定空化数时升力系数随攻角的变化规律,并对相应流场进行分析.在二维水翼空化流动的准稳态模拟中,重点对比分析了不同湍流模型、空化模型对空化流动模拟结果的影响,并将计算结果与实验数据进行比较.研究结果表明,在二维水翼附近空化流动的数值模拟中,采用k-ω湍流模型和均相流空化模型计算得到的结果与实验数据吻合较好.最后,采用分离涡模型对不稳定的空化流动进行了瞬态模拟,并得到空化流动的变化规律.     

等截面隔离段中激波串结构的数值模拟

采用二阶TVD格式的有限体积法偶合Baldwin Lomax代数湍流模型求解雷诺平均Navier stokes方程 ,数值模拟了二维矩形等截面隔离段中的流动现象 ,很好地模拟出由激波 /附面层干扰所形成的复杂的激波串流场结构。计算结果与国外有关文献的结果进行了比较。     

用动网格法计算理想平板的颤振导数

提出了一种数值模拟振动的理想平板绕流场 ,并提取其气动导数的方法。时间相关的不可压N S方程采用Projetion 2格式解耦 ,关于中间速度的动量方程的时间和空间离散采用二阶半隐格式 ,压力Possion方程迭代用多层网格法加速收敛。分别计算平板作竖向强迫振动或扭转强迫振动的气动力 ,用动网格法考虑平板和气流的耦合作用。由计算得到的气动力用最小二乘法确定 8个气动导数。计算结果和理想平板的Theodorsen理论值有较好的一致性     

飞船轴对称热化学非平衡流场数值求解

本文从轴对称热化学非平衡Ｎ－Ｓ方程出发,利用时间相关法,采用双温度、七组元反应气体模型,利用隐式ＮＮＤ有限差分格式和时间预处理技术数值求解了ＦＩＲＥＩ飞船热化学非平衡流场,得到了较为准确的结果,并分析了热力非平衡对流场的影响。     

谐波平衡法数值模拟周期性非定常流动

谐波平衡法是目前非定常流动数值研究中颇受关注的频域计算方法,该方法计算效率高,特别是对于周期性非定常流动问题具有明显的计算优势。采用谐波平衡法数值模拟高超声速导弹标模模型在不同参数下的强迫俯仰振动,重建气动力/力矩系数迟滞曲线,将计算结果与时域方法相比较,详细考察谐波平衡法的计算精度、效率和内存需求,分析不同参数对谐波平衡法计算的影响。计算结果表明,谐波平衡法适用于中小振幅下的非定常流动计算,且计算时间基本不受减缩频率的影响,对于长周期低频问题,较时域方法有明显计算优势。     

激波诱导异质气体界面失稳的数值模拟

基于二维N-S方程,利用有限差分数值离散方法,对激波诱导异质气体界面失稳的现象进行了数值模拟,与文献中实验结果和计算结果进行了定性比较,并进一步分析了整个流动的非定常动态变化特性和非线性特征.研究表明,本文数值模拟的非定常流场图谱与文献中的实验结果和数值结果吻合较好;数值结果捕捉到了六氟化硫界面的演变过程及流场中复杂的波系结构.     

三维超声速隔离段湍流内流场旋涡结构的数值模拟

从基于雷诺平均的N S方程出发 ,采用有限体积方法离散控制方程 ,数值模拟了三维超声速隔离段湍流内流场。计算中采用了二阶OC TVD差分格式、LU隐式算法和Baldwin Lomax代数湍流模型。数值结果与实验做了对比 ,并结合隔离段中的激波串结构分析了其横截面上旋涡结构的发展过程及不同外形条件下旋涡的不同结构。计算结果表明采用本文发展的方法模拟隔离段湍流流场是可行的 ,截面为正方形与长方形的隔离段内的涡旋结构截然不同     

高超声速变形飞行器翼面变形模式分析

为提高高超声速翼身组合式飞行器的射程,研究了采用不同翼面变形模式时,飞行器在马赫数3~8内的气动特性和翼面效率。针对典型的轴对称翼身组合式外形,采用Navier-Stokes方程进行数值模拟,对伸缩、变后掠和二维折叠三种变形模式下的外形在超声速及高超声速来流条件下进行模拟,并对升阻比、翼面单位面积升阻比和操稳特性进行分析。结果表明:在超声速及高超声速范围内,变后掠变形模式在宽速域内升阻比提高明显,同时具备优良的翼面效率及操稳特性,其在马赫数3~8范围内具有最优的综合性能。研究成果能对高超声速翼身组合式变形飞行器布局设计提供参考,具有一定的指导意义。     

用高阶高精度 WENO 格式求解二维激波附面层干扰流场

采用时间相关法求解二维 Navier-Stokes方程 ,数值模拟二维平板层流附面层与激波干扰流场 ,给出了物面压力分布和应力分布。计算中 ,对流项空间导数的差分离散采用高阶高精度 WENO格式 ,时间方向采用具有 TVD性质的 Runge-Kutta方法 ,粘性项采用二阶中心差分。所得压力分布和应力分布与国外实验结果吻合较好 ,计算实践表明高阶 WENO格式具有优异的性能 ,应用前景广阔     

球锥体高超声速绕流的俯仰阻尼导数的数值计算

采用Navier Stokes方程描述球锥体外形非定常振动流场 ,在Etkin理论下给出俯仰阻尼导数的计算公式。对定常流场的计算采用ADI形式的NND格式 ,对非定常流场的计算采用四步Runge Kutta方法 ,并引入变系数残值光顺技术加速收敛。将俯仰阻尼导数计算结果与实验及工程计算进行了比较 ,并数值研究了其随振动中心的变化规律。     

高超音速粘性流场数值计算中差分格式的对比分析

本文通过三维高超音速粘性流场的求解对三种常用的差分格式： Ｈａｒｔｅｎ 的 Ｔ Ｖ Ｄ 格式、修正后的 Ｏｓｈｅｒ Ｃｈａｋｒａｖａｒｔｈｙ 的 Ｔ Ｖ Ｄ 格式以及 Ｊａｎｅｓｏｎ 的当地极值递减 （ Ｌ Ｅ Ｄ） 格式对激波的分辨率、边界层计算的准确性和格式的计算效率等三个方面进行了较细致的对比分析， 为三维高超音速复杂流场的数值模拟打下基础。     

一种求解飞行器廓线间断处超音速复杂流场的方法

用有限体积法求解三维完全气体的Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,并将通量守恒的分区计算方法和Ｂ－Ｌ湍流模型修正结合在一起,捕捉到了间断处的分离涡、激波、膨胀波等流场信息     

Numerical simulation of base flow with hot base bleed for two jet models

In order to improve the benefits of base bleed in base flow field, the base flow with hot base bleed for two jet models is studied. Two-dimensional axisymmetric Navier–Stokes equations are computed by using a finite volume scheme. The base flow of a cylinder afterbody with base bleed is simulated. The simulation results are validated with the experimental data, and the experimental results are well reproduced. On this basis, the base flow fields with base bleed for a circular jet model and an annulus jet model are investigated by selecting the injection temperature from 830 K to 2200 K. The results show that the base pressure of the annular jet model is higher than that of the circular jet model with the changes of the injection parameter and the injection temperature. For the circular jet model, the hot gases are concentrated in the vicinity of the base. For the annular jet model, the bleed gases flow into the shear layer directly so that the hot gases are concentrated in the shear layer. The latter temperature distribution is better for the increase of base pressure.     

预处理方法求解低马赫数流场的边界条件

通过求解预处理后的二维可压缩Euler方程 ,数值试验了常见的三种进口边界条件的提法及其对稳定性和收敛性的影响。数值试验结果表明 ,预处理后 ,在低马赫数时基于黎曼不变量的边界条件为反射边界条件 ,而另两种边界条件为无反射边界条件