一种三维结构网格并行算法

基于三维结构网格,以MPICH2作为并行开发环境,研究了对等模式下流场的并行计算问题。在三维结构网格已经生成的前提下,针对计算区域的划分,网格输出格式和数据通信,提出了一种并行算法。算法的核心思想是将需要通讯的其他区域的单元按区号从1到n排列,并获得本区域通信单元的I-J-K编号与其他区域通讯单元的排列序号的映射关系。最后,在IBM BladeCenter HS22刀片机平台上进行了算例验证,结果表明算法具有很好的可靠性与通用性。     

神经网络在高空长航时无人机翼型多点优化中的应用

NSGAII算法在翼型多点设计中有着广泛的应用价值,然而其巨大的计算资源和计算成本限制了它的使用。为了解决这个问题,本文引入具有较强非线性映射能力的神经网络代理模型,采用实验设计结合BP法训练神经网络响应面来代替N-S方程求解翼型的性能。在实验点的数值模拟中,为了进一步节省计算资源,提高计算效率,采用网格的变形代替网格的重新划分,使得计算网格的更新速度提高了约50%。在翼型的参数化过程中,采用改进的PARSEC方法,用较少的参数实现了翼型的精确控制。为了增强神经网络的泛化能力,采用12-7-4-3-1的隐层结构。对NLF1015翼型的多点优化算例表明,此方法不仅显著降低了整个优化过程的计算量,而且对翼型的气动性能预测也具有较高的可信度,在高空长航时无人机的翼型设计中具有一定的潜力。     

不同尾翼鸭式布局远程弹在跨／超声速的气动特性实验研究

为了研究鸭式布局远程弹尾翼对气动特性的影响,设计了无尾翼,“T”型尾8翼,栅格尾翼三种尾翼布局,通过风洞测力实验研究不同布局在不同马赫数及迎角状态下对远程弹气动特性尤其是滚转特性的影响。实验结果显示：安装“T”型尾翼的模型和安装栅格尾翼的模型相比,在跨声速阶段,其升力特性优于栅格尾翼,也更利于滚转控制,但在超声速区域,栅格尾翼模型具有明显的升力特性优势,同时也容易进行滚转控制,而减小阻力是栅格翼将来需要解决的问题。     

二维翼型结冰数值计算

翼型结冰计算是随着结冰时间步长变化的准定常过程,因此在每个时间步长后都要更新各个模块的计算工作,如此反复直到达总结冰计算时间.在各模块计算中,分别采用了泊松方程法生成网格、求解雷诺时均方程计算流场、拉格朗日追踪法求解水滴轨迹、基于Messinger模型计算冰形、采用变距滤点法重构边界.引入动态链接库技术将各模块集成为结冰计算平台,并将三种冰形(楔形冰、槽状冰及混合冰)计算结果与国外软件及试验数据进行了对比.      

复合材料格栅结构的强度分析(英文)

复合材料格栅结构是由连续纤维缠绕的斜向及环向肋和蒙皮组成的结构。对ANSYS软件进行了二次开发,分别采用层合板和层合梁单元模拟复合材料格栅结构的蒙皮和肋。根据复合材料格栅结构的几何特征及其载荷分布特征,采用周期对称有限元模型,分别以最大应力准则、最大应变准则和蔡-吴准则和Chang刚度退化准则对轴压载荷作用下的复合材料格栅结构破坏过程进行了强度预测,计算结果表明采用蔡-吴准则预测并结合Chang刚度退化准则的计算结果与试验结果的一致性较好。     

环境减灾卫星国际服务模式探索与成效

环境与灾害监测预报小卫星星座是我国对地观测系统的重要组成部分。2008年9月6日,我国成功发射环境与灾害监测预报小卫星星座A、B两颗光学小卫星（简称“环境减灾-1A、1B卫星”）。环境减灾-1A、1B卫星成功发射以来,民政部国家减灾中心（民政部卫星减灾应用中心）积极组织开展卫星在轨测试、运行管理、减灾应用、用户服务与宣传推广等工作。目前,卫星运行状态良好,各项指标正常。一年来,在成功应对国内数十场新发灾害监测与评估,积极开展国内应用服务的同时,也积极探索国际服务模式,拓展为国外救灾工作提供灾害监测产品服务,并先后成功应对澳大利亚森林大火、加拿大火灾、中缅边境火灾、海地特大地震等灾害应急监测与评估工作,在国际上获得好评,展示了我国空间技术在减灾救灾工作中的良好应用,树立了我国和平利用空间技术的大国形象。     

某中空长航时无人机机载浮标分离安全性研究

研究机载浮标投放过程中气动干扰特性及浮标运动规律,对机载浮标安全投放方案的设计具有重要意义。采用计算流体力学方法对机载浮标投放过程进行数值模拟。通过动态嵌套网格技术模拟浮标运动,结合刚体六自由度运动方程求解三维非定常 N-S ( Navier-Stokes )方程,模拟基于察/打一体无人机平台的浮标自由投放过程,获得详细的包括浮标运动姿态、运动轨迹和受力情况等在内的数据信息;通过算例,验证本文数值模拟程序的准确性。结果表明:浮标投放过程中,轻质浮标口盖自由分离会对螺旋桨后推式无人机的安全性造成重大影响;考虑浮标口盖固连在浮标降落伞上的状态下,初速度为 0m / s 时投放浮标整体(包含浮标、降落伞、口盖)离螺旋桨最近距离为 0.151m ,安全余量不够;初速度为 1.8m / s 时投放浮标整体离螺旋桨最近距离为 0.671m ,满足投放安全性要求。     

一种适用于三维混合网格的GMRES加速收敛新方法

为提高流场计算收敛效率,发展了一套适用于三维混合网格Naiver-Stokes方程求解的并行广义最小残差（GMRES）隐式时间推进方法。该方法由科学计算可移植扩展工具包（PETSc）中的Krylov子空间求解器实现,线性方程系统中的系数矩阵直接以显式给出以提高算法的稳定性。为进一步提高GMRES方法的收敛速度,对非结构网格的序号进行了重排序,使得系数矩阵的非零元素尽量向主对角线靠近。利用所发展的GMRES方法,完成了对ONERA-M6机翼、AIAA阻力预测会议通用研究模型（CRM）等算例的计算,计算结果与试验结果吻合良好。通过与其他隐式推进方法进行比较,对算法的收敛特性进行了研究。结果表明,所发展的GMRES方法计算更加稳定,残差下降速度相对LU-SGS（Lower-Upper Symmetric Gauss-Seidel）方法更快,尤其是气动力系数向着收敛解逼近的速度更加明显,提高了计算效率。     

基于径向基函数与混合背景网格的动态网格变形方法

发展了一种基于径向基函数(RBF)插值和混合背景网格映射的动态网格变形方法。混合背景网格由运动物面附近的各向异性网格和偏离运动物面的各向同性网格组成。物面变形或运动后,包含在各向异性背景网格中的所有物面网格点的位移由指定的物面运动规律得到,而各向同性背景网格点的位移通过径向基函数插值获得。然后变形后的计算网格利用计算网格与混合背景网格之间固定的映射关系代数插值得到。各向同性网格距离物面一定的距离,消除了物面保形对径向基函数插值的要求,从而可以减少径向基函数基点的数量,以提高插值的效率和可靠性。最后,通过二维和三维算例对该方法的变形网格质量、参数影响、网格拓扑依赖性、复杂外形变形能力及网格变形效率进行了测试。测试结果表明,该方法能够高效生成大变形下高质量的变形网格。     

多体分离数值模拟及方法验证

利用自主研发的混合网格Euler／RANs方程解算软件WoF90以及六自由度动力学与运动学方程解算程序,采用脚本语言技术实现网格的自动生成,建立了一套可用于多体分离问题的准定常Euler流动非结构网格方法。通过AEDC算例对非结构准定常Euler方程方法解决多体分离问题进行了计算研究,验证了方法的可行性。