预处理法求解定常/非定常混合网格的全速流场

运用时间导数预处理法在混合网格上求解三维定常/非定常N-S方程。预处理后的N-S方程,用二阶迎风格式的有限体积法离散,隐式、LU-SGS迭代求解,双时间步推进求解非定常流。本文以一维预处理欧拉方程为基础进行分析,提出了一种特征变量的边界条件以适应预处理后的特征系统。从极低速到跨声速范围的众多算例表明,本文的边界条件处理稳定有效,预处理法在各种速度下都具备同样快速的收敛性,在较宽的使用范围内均能得到理想的计算结果。     

红外热像仪作用距离现场检测装置的研究

针对目前红外热像仪现场性能检测手段的不足,设计一套便携式红外热像仪最小可分辨温差(MRTD)检测装置,首先现场测量出红外热像仪的MRTD,再由大气透过率模型计算得到目标的表观温差,最后通过Johnson准则分析计算得到红外热像仪探测或识别目标的能力.试验结果表明该方法准确可行,可用于各类红外热像仪作用距离的检测.     

基于STL文件的笛卡尔网格生成方法研究

采用从STL文件中直接获取模型几何数据的方法,利用ADT数据结构快速剔除相同顶点,提出一种新的“推进查找”方法替代以往网格单元是否在物面内部的判定方法,避免了对物面的自封闭要求,无需进行几何重构,直接完成网格生成,更进一步提高了自适应笛卡尔网格生成速度和自动化程度,为实现CAD与CFD的之间的数据对接提供了一种新的方式。     

电弧喷射推力器电磁场模拟

为了揭示电弧喷射推力器内部电磁场特征及其与高温电离气体间的相互作用机理,建立了比较完善的电磁场模型,详细讨论了电磁场数值模拟方法及其计算效果.模型同时考虑了电场力、洛仑兹力和电子压强梯度对传导电流的贡献,由电子动力学方程和粒子微观碰撞理论推导出电流方程,由麦克斯韦方程组推导出电磁场控制方程,采用包括强隐过程(SIP)迭代法在内的十种迭代方法求解离散的电磁场方程.研究表明,SIP迭代法是适用于电弧喷射推力器电磁场模拟的快速有效的方法,模型能够准确地反映推力器内部电磁场与高温电离气体间的相互作用机理.     

低功率氩电弧喷射推力器中的辐射

为了研究低功率氩电弧喷射推力器中的辐射损失,采用辐射模型和非辐射模型对其工作过程进行了对比数值模拟分析。采用二阶精度无波动、无自由参数的耗散差分格式(NND格式)数值求解耦合电磁源项和辐射源项的N-S方程组,并采用隐式残差光滑法加速收敛;采用有限控制容积积分方法离散求解椭圆型偏微分电磁场方程,并采用逐点超松弛迭代方法加快收敛速度。数值模拟结果对比给出了辐射模型和非辐射模型流动分布情况,并比较了两种模型推力器推力、比冲和推进效率。研究结果表明,与非辐射模型相比,辐射模型比冲低0.137%,推进效率低1.03%,辐射对低功率电弧喷射推力器性能影响比较小。      

低功率电弧喷射推力器实验研究

为深入了解低功率电弧喷射推力器的启动特性和工作性能，对不同推进剂和不同结构尺寸的推力器进行了广泛的多启动实验研究。研究中建成了推进剂供给系统和真空模拟系统，研制了电源处理单元和小推力测量装置，开发了实验数据采集系统，完成了实验推力器结构设计和材料选择，并分别采用氩、氮、氨为推进剂成功地进行了点火实验。研究表明，推进剂性质对推力器启动特性和工作特性具有重大影响，推力器优化设计应该针对具体的推进剂进行并考虑适当宽广的流量范围。     

雷诺应力模型在三维湍流流场计算中的应用

从雷诺应力模型出发,通过求解雷诺平均N-S方程组获得三维湍流流场的数值解。计算中比较了多种湍流模式,并进行了相应的流场计算。本文完成了两个典型算例。从算例1的计算结果与实验值比较中发现:采用雷诺应力模型(RSM)计算的三维流场比采用k-ε模型更贴近实验值;算例2采用了RSM模型及三维非结构网格,对一典型内流问题进行了三维流场计算。算例的数值实践表明:采用雷诺应力模型可以有效的计算各向异性的湍流流场;另外,发展非结构网格有助于模拟壁面附近的流动,并节省计算机内存。      

化学反应动力学对超声速燃烧的影响

非平衡化学反应是超声速燃烧中非常重要的过程,其反应动力学模拟直接影响其内部温度分布。针对某支板氢燃料超声速燃烧室,分别采用了单步化学反应动力学、7组分8步反应和9组分19步反应动力学模型,对比了燃烧反应速率、内部温度场和速度场的分布。结果表明,采用详细化学反应动力学,燃烧热释放缓慢,在流场中下游自由基参与反应,其温度场与试验值吻合更好。     

特高压输电塔线体系气弹模型设计与风洞试验

针对特高压输电线路塔线体系风致响应较大的特点,以某输电塔线工程为原型,按照空气动力学弹性相似准则,设计制作了输电线路塔线气动弹性模型,进行了多个风向角、多个风速下单塔与塔线体系在紊流场中的风洞试验.试验结果表明,对于跨度不大的塔线体系,较低风速下(风速约<20-25 m/s)塔线中塔的加速度要高于单塔,当风速继续增加时,单塔的加速度将超过塔线体系.塔线体系总的风振响应需要考虑塔线耦合作用的影响.     

再入飞行中DSMC与Navier-Stokes两种模型的计算与分析

首次提出了再入飞行过程中"小Knudsen数特征区(即[Kn1,Kn2]域)"的概念,并用直接模拟蒙特卡罗(DSMC)算法与求解Navier-Stokes方程两种模型详细研究了在这个区域中的流动.数值计算表明:采用王保国教授课题组自己编制的DSMC源程序和求解多组分、考虑热力学非平衡以及化学反应非平衡、弱守恒型Navi...