固体火箭发动机燃烧室三维流动数值计算

采用ｋε双方程湍流模型，以ＳＩＭＰＬＥ 计算程式求解翼柱型装药固体火箭发动机燃烧室内三维非定常不可压流ＮＳ 方程．在建立翼柱型装药简化模型的基础上采用边界标志法来表达燃面推移，采用从二维到三维的初场给定方法结合多重网格法求出了多个时间步的非稳态流场结构．计算结果表明燃烧室内旋涡运动呈现一定的空间与时间分布，周向加质导致燃气通道横断面上轴向速度分布的极不均衡．     

三维抛物化N-S方程全隐格式矢通量分裂求解

采用全隐格式的矢通量分裂算法,在通用曲线坐标系内求解三维无因次抛物化强守恒型N-S(Navier-Stokes)方程.所建立的计算程序,可计算较为复杂边界的三维流场,如求解亚、跨、超音速共存的流场,分区模拟内外流共存流动及其混合剪切尾流.给出了一轴对称亚跨超音速喷管内流场的计算结果,包括计算网格,流矢分布,压力、温度和马赫数等值线分布.也给出了一喷管内外流及其混合尾流的计算结果,包括计算网格,流矢分布,无因次压力、密度和温度等值线分布.     

三维四向编织复合材料弹性性能的有限元预报

在四步法方型三维编织复合材料细观结构单胞模型的基础上,假设纱线具有六边形横截面,应用ANSYS软件,建立了实体有限元模型,该模型表面与预制件表面平行,更加符合三维编织复合材料的实际结构,且有利于力学性能的分析,同时,该模型中还考虑了织物的编织角与内部编织角的几何关系.基于该模型,计算了材料的弹性常数,分析了编织角和纤维体积含量对弹性常数的影响规律.结果表明,数值计算与实验结果吻合较好.此外,还确定了材料内部的应力场分布,为材料的强度预测奠定了基础.      

基于SA湍流模型的DES方法应用

使用基于SA(Spalart-Allmaras)湍流模型的分离涡模拟(DES,Detached-Eddy Simulation)方法,对圆柱绕流和超声速圆柱底部流动等高雷诺数大分离流动进行了数值模拟.无粘通量的计算采用二阶迎风格式,时间推进使用二阶隐式双时间步法.计算表明相对于雷诺平均N-S方程(RANS,Reynolds-Averaged Navier-Stokes)模拟,时空均为二阶精度的DES模拟能够分辨更为精细的分离区三维旋涡结构,对分离区主要湍流结构的捕捉和平均压力分布的预测取得了满意的结果.     

无振荡、无自由参数格式在三维太阳风流动数值模拟中的应用

从快速、有效地进行空间天气数值预报的需要出发,针对1998年5月份行星际太阳风暴事件,采用全三维流体力学(HD)模型进行数值试验,考察了无振荡、无自由参数(NND)格式在三维太阳风流动数值模拟中的应用.计算中内边界的密度分布由观测K日冕数据来确定,由此得到的源表面密度分布具有和源表面电流片相似的结构.数值试验表明: 虽然该格式不需要加入人工粘性,但是具有很好的计算稳定性.在扰动计算时激波间断耗散小,在三维任何方向上间断所占的网格数比较少,没有数值色散现象.      

固体发动机实际成型药柱燃面退移快速算法

针对固体发动机电子计算机断层扫描(CT)图像数据具有伪影噪声的问题，和实际成型药柱燃面粗糙度大、退移计算难度大的问题，提出一种CT数据的快速最小距离函数(CTFMDF)法。固体火箭发动机CT图像中存在伪影噪声，采用非均值滤波(NLM)算法对CT图像进行去噪处理，采用Scharr算子对去伪影后的图像进行Canny边缘检测，提取装药燃面。最大类间方差(OTSU)算法将装药分离，并建立三维装药体数据模型，对燃面数据建立多个并行K-d树，快速检索出装药到燃面的最小距离。实验结果证明：对于不同的装药结构，可以完成任意燃去厚度时的燃面位置，且CT-FMDF法运行时间更短。基于CT实测数据，对于带有初始燃面缺陷的装药，可计算出燃烧时缺陷对燃面的影响。     

基于危害度航空发动机可靠性评估模型及应用

针对航空发动机可靠性指标评估的需要,提出了基于危害度故障统计分析模型.该模型考虑各故障样本对发动机危害性的影响,将出现的故障按照对发动机安全、性能、任务及维修等指标的影响进行等级分类,从故障的失效机理出发建立相应的分布模型,采用分布计算和二次分布等算法进行整体可靠性指标评估.应用该模型对英产和国产的2种同类型航空发动机进行了故障统计分析,计算了2种发动机的可靠性指标,并进行了对比分析.最后以累积工作时间和日历工作时间进行了故障的时间变化分析,提出了故障样本的时效性.      

面向地面运动目标的无人机协同航迹规划

由于现有无人机协同航迹规划主要局限于固定目标,提出了一种面向地面运动目标的无人机协同航迹规划方法.给出面向地面运动目标的无人机协同航迹规划总体方案;采用按时间推进的自规划与同步规划交替进行的协同规划机制,满足协同信息交互存在时间间隔的约束;在规划中将无人机运动与地面目标运动结合并构造适合于地面运动目标的启发函数,利用启发式搜索算法快速生成航迹段;采用分布计算协同变量的方法降低了协同计算量.仿真结果表明该方法可以使无人机通过信息交互和自主分布计算快速生成面向地面运动目标的三维协同航迹,且满足规划的实时性要求.     

Hydrothermal instability of thermocapillary convection in large-Prandtl-number liquid bridges under microgravity

Linear stability analysis was performed to study the mechanism of transition of thermocapillary convection in liquid bridges with liquid volume ratios ranging from 0.4 to 1.2, aspect ratio of 0.75 and Prandtl number of 100. 2-D governing equations were solved to obtain the steady axi-symmetric basic flow and temperature distributions. 3-D perturbation equations were discretized at the collocation grid points using the Chebyshev-collocation method. Eigenvalues and eigenfunctions were obtained by using the Q–R method.The predicted critical Marangoni numbers and critical frequencies were compared with data from space experiments. The disturbance of the temperature distribution on the free surface causes the onset of oscillatory convection. It is shown that the origin of instability is related to the hydrothermal origin for convections in large-Prandtl-number liquid bridges.     

集气腔总压对电弧喷射推力器工作过程的影响

为了研究集气腔总压对电弧喷射推力器工作过程的影响,在分析其工作机理的基础上采用化学非平衡流动和稳态电磁流体电磁场模型对不同集气腔总压下推力器工作过程进行了数值模拟.流体力学方程组和电磁场方程考虑了多种流动机理及电磁场与高温电离气体的相互作用,化学动力学模型考虑了各种碰撞反应.采用二阶精度NND格式求解流体力学方程组,采用有限速率化学反应模型计算组分生成率,采用交替方向隐式(ADI)超松弛迭代法求解电磁场离散方程.给出了不同集气腔总压下推力器内部参数分布及其宏观性能.研究表明,集气腔总压对推力器工作过程具有多方面的影响,在保证电弧稳定的情况下,适当提高集气腔总压可同时提高比冲和推进效率.     

COSMIC PLASMA DYNAMO

A new dynamo model based on the polarization of plasma is presented in this paper. From the Maxwell equations in a moving medium, a magnetization vector can be caused with Rongon current. The steady solar magnetic field is solved from the equations. On the assumption that the meridianal flow is ignored, the distribution of magnetic field is put out. In the model, there is no additional parameter considered. The intensity of magnetic field inside the sun ranges from 1-6 T. The surface magnetic field around the pole is in the order of 1 ×10-3 T, at low latitude the calculated surface magnetic field has the order of 1×10-2 T. The maximum magnetic field is around 30° in latitude.     

使用Gao-Yong湍流方程计算翼体角偶流动

Gao-Yong湍流方程基于侧偏统计平均方法保留了湍流脉动量的一阶统计信息,并引入加权漂移速度对称性及正交各向异性,导出漂移流连续方程、动量方程和机械能方程,最后依据湍流物理的唯象论,使用动量传输链概念模化封闭了整个方程组.大量算例验证了此方程对广泛范围的复杂湍流问题的适定性.为了进一步的工程应用,将其加入开源计算流体力学软件OpenFOAM(Open Field Operation and Manipulation)中,并就翼体角偶流动进行了数值模拟.翼体角偶流动具有三维分离、马蹄涡等十分复杂的流动特征,通过与实验测量结果的对比研究表明,Gao-Yong湍流方程不论在定性还是定量上,都成功捕捉到这类复杂流动现象的主要特征,得到了令人满意的结果.     

电离层Es离子分布数值建模

从离子连续性方程和动量方程出发,比较全面地考虑了太阳光电离、星光电离、地冕电离、星际背景电离和流星离子流等电离源,综合分析Es层的主要动力学过程和光化学过程,以风剪切理论为基础,研究中性成分碰撞、电场作用、金属离子作用和分子离子作用等机制对离子分布的影响,建立了一维时变电离层Es物理模型.对离子产生率、垂直方向的离子速度在高度和时间上的变化进行仿真和计算,得到了在电场及风剪切作用下的离子密度剖面24h内的变化规律.根据建立的模型,以昆明相干散射雷达观测的径向风结果作为模型输入参考,仿真并得到了电离层Es电子密度的时空分布,据此反演出电离层Es临界频率f₀E_s.与昆明站电离层测高仪同时间观测的情况进行对比,结果比较一致,验证了建立的Es层物理模型正确性.      

基于N-S方程串并行计算的机翼优化设计

对典型的大展弦比和中等展弦比机翼在0.7~0.9马赫数区间进行了串行和并行计算的数值校验,证实了以N-S (Navier-Stokes)方程为主控方程的串行和并行流场求解器的正确性,并讨论了并行效率和加速比.结合Powell算法,讨论了在确定的机翼平面形状和翼型的条件下,以升阻比最大为目标的三维机翼截面翼型最大厚度与扭角的优化设计.算例结果表明,厚度的非线性分布和负的扭角会改善机翼流场的流动状态,使机翼的升阻比得到提高,优化设计方法是可行的.      

基于热-流-固耦合方法火焰筒壁温三维数值模拟

准确地确定火焰筒壁温分布是对其进行热结构分析的前提和基础.基于ANSYS/CFX程序平台,采用基于燃烧模型的热-流-固耦合有限元方法,对某环形回流燃烧室的整个流场、温度场进行数值模拟,得出火焰筒壁温三维温度场分布.计算中考虑了燃烧热源和热辐射,燃烧和流场变化对换热的影响,以及燃气、煤油和固体物性参数随温度的变化.数值模拟结果表明,将流场和固壁进行耦合求解,能更精确反映流场和温度场的整个形态,可以模拟出较为合理的流场和温度场分布以及相应的流动换热特性.所做工作对火焰筒结构设计具有重要的工程应用价值.      

使用GAO-YONG湍流方程组 对顶盖驱动方腔流的计算

运用SIMPLER方法、QUICK格式及交错网格技术求解GAO-YONG湍流方程组,对二维顶盖驱动方腔流进行数值模拟,得到方腔内的涡量分布和方腔中心水平剖线上的速度变化曲线,并与DNS(Directly Numerical Simulation)结果进行定性对比,有较好的一致性.结果表明,GAO-YONG湍流方程组中与湍流多尺度现象相对应的机械能方程对此湍流方程组能够合理地模拟出顶盖驱动方腔流中复杂的湍流粘性分布以及湍流能量逆转现象起到了关键作用.研究进一步证明了GAO-YONG湍流方程组能够得到复杂流动的真实粘性场.