非结构网格上欧拉方程的自适应逆风隐式有限元解法

本文发展了一种解二维欧拉方程的隐式逆风有限元格式。空间近似应用的是二阶精度逆风方法；时间近似上根据数值通量函数的线性化处理导出了有限元的隐格式。为减少数值振荡，引入了基于特征变量的限制因子。在经自适应处理的非结构网络上，本文对跨、超、高超声速绕流进行了数值模拟。计算实践表明发展的隐式有限元格式与显格式相比具有较高的计算效率和较强的适应性。     

飞机蒙皮机器人化制造曲面分片算法研究

为实现飞机曲面机器人加工的自动轨迹规划,研究了基于模型几何特征的飞机曲面分片算法。通过分析飞机构件加工的工艺参数信息,以表面加工允许的最大曲率或最大偏角为阈值,提出基于近似曲率的三角网格模型分片方法。以三角网格中任意面片为起点,沿网格邻接关系向外拓展联接,形成若干近似平面化的模型分片结构,当超过阈值时则停止拓展,从而得到一整片曲面。通过对飞机CAD模型进行分片,发现采用本文提出方法可以保证分片结果的均匀性,为后期加工轨迹自动规划规划提供保障。     

立式捏合机桨叶结构参数对混合釜流场影响的仿真分析

为揭示立式捏合机桨叶几何结构参数与混合釜流场之间的关系,以1 L两桨立式捏合机为研究对象,采用计算流体力学ANSYS Fluent 14.5软件对立式捏合机进行了三维模型仿真研究,分析了桨叶结构参数(桨桨间隙、桨叶螺旋角)的改变对混合物料挤压应力和剪切应力的影响。结果表明,增大桨桨间隙可减小桨叶对混合物料的挤压应力和剪切应力,空心桨螺旋角分别为40°、45°、50°,桨桨间隙由1.0 mm增加到3.0 mm,挤压应力、剪切应力分别平均减小了82.4%、57.23%;增大桨叶螺旋角可增大桨叶对混合物料的挤压应力和剪切应力,空心桨间隙分别为1.0、2.0 mm,桨叶螺旋角由35°增加到55°,挤压应力、剪切应力分别平均增加了92.8%、55.7%。     

无网格方法关于导数计算的程序验证

为了提高无网格方法对空间导数的数值计算能力,提出并论证了新的适用于直角坐标系、柱坐标系和球坐标系的空间导数无网格算法。基于数值解与理论解的对比分析,对数值误差和收敛性进行了后验评估。评估结果表明:该算法对函数和导数的估计均为2阶精度,支撑域尺度对数值误差的大小有一定影响,但不影响数值计算的精度等级。在离散尺度h为0.01的条件下,对所选函数及其导数的数值计算相对误差不大于0.65%。      

自适应直角坐标网格方法在对流换热数值模拟中的应用

复杂几何边界是对流换热计算面临的主要难点之一。自适应直角坐标网格方法可以比较方便的处理复杂几何边界,但其应用目前主要集中于可压无粘流的Euler方程、可压粘性流的Navier-Stokes方程的外部绕流计算。本文应用自适应直角坐标网格方法计算了两个不可压内流对流换热算例,并把计算结果同商用CFD软件STAR-CD的计算结果进行了对比。自适应直角坐标网格方法在本文的两个算例中得到了同STAR-CD一致的结果。这说明采用该方法数值模拟对流换热问题是适用的。      

发动机喷管内流场对流换热系数影响因素的数值分析

基于Roe-FDS格式,采用标准κ-e湍流模型对某液氧煤油发动机喷管内部燃气流对喷管内壁的对流传热换热系数进行了数值仿真与研究,分析了网格雷诺数以及恒温壁壁温对喷管内壁对流换热系数的影响.当近壁面网格第一层法向高度为10-3 mm数量级时,仿真结果与工程估算公式的计算结果相吻合.同时,在一定范围内,对流换热系数随来流速...     

固体姿轨控针栓推力器内流场数值仿真

为研究针栓推力器内.流场的特征和规律,采用二维轴对称模型,开展了针栓运动过程的非稳态及稳态流场仿真,分析了不同开度下的流场结构、针栓所受气动力及流场对针栓几何参数的敏感性.结果表明:由于针栓的干扰作用,流场结构呈现出复杂激波系、顶端流动分离等现象;气动力系数与针栓头部锥角大小呈正相关关系,与集气室压强无关;节流面形状对气流参数分布、激波强度及位置、气动力等表征的流动模态影响很大.     

变构型航天器参数化前处理方法

针对变构型航天器结构不固定的问题,基于VC 6.0软件对Trail及Session进程文件开发的方式,研究了利用Pro/E及Patran软件对航天器进行参数化建模及网格划分的具体操作方法,通过IGES中间文件的形式,将Pro/E及Patran相连接,利用VC开发工具建立了一套变构型航天器参数化前处理平台。采用两舱构型航天器对参数化前处理方法及平台进行了验证,结果表明:该参数化建模及网格划分方法切实可行,可为变构型航天器的静力学、动力学、热、结构等分析过程的前处理方法提供参考。     

壁面函数在激波诱导分离流动中的应用

以数值模拟激波-附面层干扰引起的流动分离问题为研究背景,发展了基于有限体积方法的雷诺平均Navier—Stokes（RANS）方程的流场数值模拟方法。利用壁面函数模型得到壁面剪切应力,通过修正壁面粘性通量,构造了一种新的湍流边界处理方法,并将其耦合到RANS方程和SSTk-ω湍流模型的数值求解中;同时,针对激波诱导引起的附面层流动分离问题,提出一种附面层网格加密技术,能够自适应加密分离区内附面层网格,使得在流动分离区域也能够使用壁面函数模型。数值算例表明,壁面函数模型能够降低数值模拟结果对网格的依赖性;同时也验证了壁面函数耦合附面层网格自适应方法,在处理激波诱导引起的附面层流动分离问题时的有效性和准确性。     

稳定性导数计算方法研究

稳定性导数的准确计算对于飞行器的操稳特性具有重要意义.应用计算流体力学(CFD)方法中的非结构化动网格技术建立能够模拟飞行器做周期性俯仰运动的强迫振荡法,以国际动导数标模Finner导弹为验证算例,获得不同马赫数下俯仰力矩系数的迟滞环曲线,进而计算Finner导弹在不同马赫数下的静稳定性导数和动稳定性导数.结果表明:本文计算得到的俯仰静稳定性导数与试验数据非常接近;在亚音速和超音速范围内,动稳定性导数计算结果与试验值很接近,但在跨音速范围内,本文计算结果与试验曲线的规律性一致,但误差较大.