超声速条件下内埋式武器分离特性的数值分析

采用非结构动网格技术,耦合求解欧拉方程和六自由度(6-DOF)运动方程,对内埋式导弹在亚声速和超声速条件下的分离轨迹和姿态变化进行了比较,分析了不同马赫数下气动力对导弹分离特性的影响,重点讨论了超声速条件下导弹分离过程中相对于载机的位置和导弹的姿态变化情况,分析了超声速飞行环境下导弹安全分离的条件.     

三维复杂几何外形的表面三角形网格自动生成

按点－线－面的层次数据结构,选取几种有代表性的曲线曲面,实现了三维复杂外形的表面描述。采用曲面到平面的保形变换,在平面内由阵面推进法自动生成三角形网格,实现了三维复杂外形的表面三角形网格自动生成。     

粘性流体换热问题用直角坐标网格的贴体解法

采用直角坐标网格但通过在边界附近保留不规则控制体的非结构化直角坐标网格方法求解粘性不可压流体流动和换热问题.在边界附近切割单元的求解采用常规的方法处理,而不采用插值方法处理,并在同位网格上对控制方程组进行离散.因而该算法是贴体的,并可适用于解决复杂边界流动换热问题.通过对同轴环行空间和二维倾斜空腔内的对流换热问题的计算对该数值方法在传热问题中的应用进行验证.对比结果表明,计算结果与精确解和STAR-CD的结果在一定Re数和网格数时是很接近的,可以满足一定的精度要求,说明该数值计算方法可用于传热问题.      

基于非结构重叠网格的二维外挂物投放模拟

本文发展了基于非结构重叠网格的动态网格技术,并提出了新的相邻单元搜索算法。非结构重叠网格集非结构网格和重叠网格的优点于一身,适合处理复杂外形和多体相对运动问题。新的相邻单元搜索算法在碰到物面时能继续搜索下去,搜索效率大大提高。非定常Euler方程的求解,采用有限体积法,时间上采用双时间步长推进。为验证算法的有效性,先对绕振荡NACA0012翼形的流动进行了模拟;然后通过耦合求解刚体三自由度运动方程,模拟了二维外挂物投放过程,结果表明非结构重叠网格在处理多体相对问题时具有优越性。     

非结构网格上N—S方程解法研究及在多段翼型绕流中的应用

在非结构网格上求解N－S方程时，靠近物面的网格伸展比相当大，用有限体积法进行显式求解时不易收敛并且雷诺数有一定的限制。本文采用的人工耗散项中含有网络的方向尺度，因而在不同区域选取不同：在粘性区内该项比物理粘性项小得多，在粘性区外与Jameson的人工耗散项相同。生成非结构网格采用的是层推进技术与传统的阵面推进法相结合的方法，这样能保证在粘中性区内网格有一定的方向规则，而在粘性区外网格完全非结构化。     

二维有限元网格全自动生成

利用ＰＥＤＳ（ＰａｒａｍｅｔｒｉｃＥｎｇｉｎｅｅｒｉｎｇＤｅｓｉｇｎＳｙｓｅｍ）造型系统的几何，拓扑信息实现了二维有限元网格自动生成，在实现过程中，改进了Ｗａｓｔｏｎ三角化算法，把此方法应用到有限元剖分中，还利用了迭代的Ｌａｐｌａｃｉａｎ算法对网格进行了均匀化处理，另外，还提供了可靠，快速对网格编辑，局部细化的方法，最终使有限元网格生成完全自动化，而且速度快，质量好，密度易于控制。     

CFX在军机微小部件优化设计中的应用

针对某型飞机在跨声速飞行中出现垂尾方向舵的铆钉被吸出的问题,本文采用CFX 5.7对全机加细微尺度的扰流条构型进行了一定状态的精细数值模拟。详细分析扰流条改进前后的流场差异,对原来的扰流条改进设计,通过流场的变化判断改进设计是否达到预期目的。计算前处理的网格采用ICEMCFD 5.1网格生成器生成HEX六面体结构网格,网格规模为450万。由于跨声速流场是风洞试验和数值计算模拟的难点,为了准确地反映全机复杂流动的物理现象,本文数值模拟选取M a=0.95,4°迎角,SST湍流模型,二阶差分格式,采用小的物理时间步长在600个时间步以内即可获得跨声速流场的收敛解。通过对流场计算结果的详细分析,证实了方案优化设计的合理性和有效性。本文的计算是CFX 5.7对复杂构型细微尺度变化流场模拟的一个典型算例,充分反映了ICEM 5.1中HEX六面体结构网格对流场细微变化的捕捉能力和CFX 5.7全隐式求解器的良好鲁棒性和快速的收敛性。对于这种特殊部件细微尺度变化的流场模拟,风洞试验显得无能为力,只有通过高精度流体仿真软件来模拟优化设计的效果,为优化设计提供参考的依据。     

基于三维非结构网格的连续伴随优化方法

基于非结构网格,采用Euler方程、连续伴随方法以及自由形面变形技术(Free-form deformation,FFD)结合序列二次规划(Sequential quadratic programming,SQP)优化算法,开发出了三维飞行器的气动外形优化代码。应用Fortran 90语言,动态存储数组和指针,节省了内存。采用Roe形式的二阶迎风格式和LU-SGS隐式时间离散格式求解连续伴随方程,同时采用一种新的参数化坐标求解方法,使得FFD技术更加鲁棒。和有限差分方法求得的梯度进行比较,验证了连续伴随方法所求梯度的正确性,实现了ONERA M6机翼的减阻优化,证明本文开发的优化代码能够实现三维气动外形的减阻优化。     

并行化非结构重叠网格隐式装配技术

重叠网格技术广泛应用于复杂外形及多体分离问题的数值模拟,但是超大规模重叠网格装配的鲁棒性和效率依然是实际应用中的瓶颈。基于格心型的有限体积格式,发展了一种并行化的非结构重叠网格隐式装配技术。该技术采用节点的壁面距离作为插值边界的判断准则,并通过物理边界推进确定出活跃区域。此外,基于网格分区策略实现了隐式装配技术的并行化,为了能够适应超大规模的计算网格,在每一个网格子区内单独建立用于查询的交替数字树（ADT）数据结构,整体网格点在各个子区内进行查询,并通过信息通讯接口（MPI）传递查询结果并进行逻辑判断。该并行化的非结构重叠网格隐式装配技术具有较高的自动化程度,通过多体装配算例以及机翼外挂物分离算例,证明了该方法的实用性。     

非结构网格下可压缩黏性流体的隐式算法

采用牛顿线化和Roe近似Jacobians矩阵隐式算法,实现基于非结构网格的三维定常黏性流动问题的求解.分别对平板、跨声涡轮叶片Mark-Ⅱ和CW-22这三种具有代表性的流动进行数值模拟,并详细分析了流动细节,模拟结果与理论以及实验值吻合良好,准确预测了流场内激波、二次流等重要的流动现象.验证了基于非结构网格的可压缩流体隐式计算方法对于黏性流动问题的计算是正确的.