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

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

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全隐式求解器的良好鲁棒性和快速的收敛性。对于这种特殊部件细微尺度变化的流场模拟,风洞试验显得无能为力,只有通过高精度流体仿真软件来模拟优化设计的效果,为优化设计提供参考的依据。     

J2项摄动条件下线性化的编队卫星动力学方程

考虑J2项摄动力对参考卫星轨道的绝对影响和对伴随卫星轨道的相对影响,应用拉格朗日方法,在笛卡尔坐标系下,给出了一组考虑J2项摄动的线性化的编队卫星相对动力学方程。该方程考虑了J2项摄动力带来的相对运动平面内外的耦合和短周期项的影响,有效地提高了相对动力学模型的精度。仿真结果表明,该方程组可有效地描述J2项摄动条件下编队卫星的相对运动,其误差只有二体非线性模型的10%,利用其设计的水平圆编队半径的相对误差不超过1%,并且由于是线性化的微分方程描述,可方便地应用到编队卫星导航、制导与控制系统的设计中。     

应用于DSMC方法的直角网格技术

研究了应用于直接模拟蒙特卡洛DSMC(Direct Simulation Monte Carlo)方法的直角网格技术.通过对边界三角形表面微元进行统计,建立直角网格缺陷单元与表面微元之间的映射关系,并进一步将缺陷单元区分为一级缺陷单元以及二级缺陷单元,同时结合可变时间步长技术,从根本上解决了在DSMC方法中利用直角网格技术贴体光滑地匹配复杂几何型面的问题.给出了相应的自适应网格调整方法.基于以上策略,编制了通用DSMC计算程序RGD-Tool(Rarefied Gas Dynamic Tool).利用RGD-Tool程序对典型的圆球绕流以及无限薄斜挡板绕流两个算例进行了验证计算,结果表明所讨论直角网格技术是有效的.     

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

研究了非结构化直角坐标网格全贴体求解粘性不可压流体自然对流换热问题。该方法在于利用直角坐标网格但通过在边界附近保留不规则控制体,使得算法是完全贴体的。利用此数值方法对二维倾斜空腔内的自然对流及上固体边界移动时的对流换热问题进行计算,并与STAR-CD的计算结果进行比较。对比结果表明,该方法计算的结果与STAR-CD的结果很接近,可以满足一定精度要求。      

基于背景网格簇的动网格生成方法

提出了基于背景网格簇实现网格变形的新方法。背景网格簇由任一内边界节点和远场边界角点生成的背景网格构成,并随边界运动而变化。通过保持计算网格节点在每一个背景网格中所在单元的面积坐标不变(三维时保持体积坐标不变)而求出一组期望位置,加权平均该组期望位置,确定出该计算网格节点的新位置。其中权值与内边界节点和计算网格节点距离的倒数相关。新方法引入的背景网格簇可以有效地改进单个背景网格的不足之处。算例结果表明:基于背景网格簇的动网格生成方法实现简单,与弹簧近似法相比,新方法因不需要迭代求解方程组而非常高效,且拥有更强的网格变形能力;与Delaunay图映射法相比,该方法背景网格的单元个数极少,因此易于定位,且不会出现背景网格单元交叉的现象,网格变形能力更强,变形后的网格质量更好。     

利用二、三维嵌套技术数值模拟复杂边界下的流场

通过数值模拟复杂几何区域和边界条件下二、三维嵌套紊流模型,获得了计算区域内任意一点的流场信息。这一方法克服了常规的单纯模拟二维或三维紊流控制方程组而无法获得准确的自由表面和三维紊流场的缺点,通过两者有机的结合,并采用动边界处理等先进的计算技术,得出了精度更高、更为合理的紊流场。该方法曾被用于三峡工程大江截流完成之后导流明渠通航的流场模拟中,计算结果合理可信,部分成果的精度也已得到了实测资料的验证,成果为导流明渠通航路线的选择提供了依据     

基于混合笛卡儿网格方法的可压流动数值模拟

以可压缩黏性流动的数值模拟为研究背景,发展了一套自适应混合笛卡儿网格(AHCG)方法以及基于有限体积方法的雷诺数平均Navier-Stokes(RANS)的数值求解方法.为更好地模拟边界层的黏性流动在近壁面处采用贴体结构网格,剩余计算区域自动生成与之相重叠的笛卡儿网格,并同时发展了基于流场特征的笛卡儿网格自适应技术.结合ADT(alternating digital tree)算法显著减少了网格生成中“挖洞”和“贡献单元”搜索的消耗机时,50万左右的网格数目下,搜索耗时为0.062s,仅为普通遍历方法的1/1847.通过二维圆柱与两段翼型绕流的数值算例显示,定常AHCG方法能够准确地预测物面压力分布与升阻力系数并且具备处理复杂外形的能力;同时通过二维非定常圆柱绕流问题与NACA0015矩形机翼翼尖尾涡的捕捉算例显示,结合了动态自适应网格加密的非定常AHCG方法尤其适用于旋涡主导流动.      

基于混合笛卡尔网格的翼型绕流数值模拟

在对重叠网格技术研究的基础上,一套以笛卡尔网格为主体、贴体的结构网格为特殊固壁边界的混合笛卡尔网格(Hybrid Cartesian Grid,HCG)生成方法得到了发展,同时基于流场特征的笛卡尔网格自适应技术和有限体积方法的雷诺平均Navier-Stokes(RANS)的数值求解方法也得到了开发。通过对经典跨音速RAE2822翼型、双NACA0012翼型和高升力两段NLR-7301翼型进行绕流数值模拟,获得了较好的结果,验证了方法的可行性与精确性,为研究二维复杂外形流动问题提供了一种解决方法。     

三维非结构贴体直角网格的N-S方程数值模拟

基于八叉树数据结构,采用直接分层光顺的方法,解决了非结构直角网格在物面不正交贴体的难点,发展了任意复杂外形的三维非结构贴体直角网格,并且可以根据流场特征进行自适应处理。在此基础上采用最小二乘法对N-S方程进行了MUSCL有限体积法的非结构网格的粘性求解,提出了带约束边界条件的最小二乘法的处理办法。计算表明,三维非结构贴体直角网格的粘性求解,可以适应各种高超声速的复杂外形,计算结果与试验吻合较好。