一种基于求解椭圆型方程的结构动网格生成方法

基于椭圆型网格生成法,实现了一种简单高效的贴体结构动网格生成方法,可用于具有移动边界问题的非定常流动数值模拟。该方法提出,在网格变形过程中,Poisson方程需要的控制网格间距和正交性的源项可以通过提取已知的静态网格源项直接得到,并在整个动网格生成过程中保持不变。因此,在椭圆型网格生成中需要通过外迭代确定源项的过程可以得到省略,而且该方法不需要人工指定参数。这使得方法具有高效和易于嵌入到已有程序中的特点。数值模拟结果证明,采用这种方法获得的网格能够较好地保持静态网格原有的正交性和光滑性,在相同迭代步数约束下,网格求解效率低于传统弹簧模拟法,但鲁棒性优于弹簧模拟法。     

旋翼桨尖涡生成及演化机理的高精度数值研究

为了细致捕捉直升机旋翼桨尖涡的生成和演化过程,建立了一个基于高精度网格和高阶通量计算格式的旋翼桨尖涡计算流体力学(CFD)求解方法。在该方法中,流场求解选取旋转坐标系下的Navier-Stokes方程为控制方程;空间离散采用迎风Roe格式,并采用低耗散的5阶WENO(Weighted Essentially Non-Osciltatory)格式进行对流通量的计算;时间推进则采用双时间法,在伪时间步上使用隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)推进格式;应用嵌套网格方法实现桨叶网格和背景网格的数据交换。应用所建立的方法对悬停状态的旋翼桨尖涡流场进行了高精度模拟,在桨叶网格上精细地捕捉到了桨尖涡的具体生成过程,在背景网格上捕捉到了更多圈数的桨尖涡尾迹,并对桨尖涡的演化机理进行了深入研究。结果表明:建立的高精度数值方法能够有效地对旋翼桨尖涡的生成和演化过程进行细致模拟;悬停状态下旋翼桨尖气流在上下表面压力梯度的作用下经历了边界层增厚、逐渐卷起形成涡核以及最终脱离桨叶形成桨尖涡的过程。     

不对称电网故障下双馈风力发电机的比例-积分-谐振滑模控制器设计

针对不对称电网故障下,双馈风力发电机的控制策略进行了研究。讨论了电网故障时转子侧变流器(RSC)和网侧变流器(GSC)的控制目标,合理安排了电网严重故障时的控制优先级。基于比例-积分-谐振滑模控制原理设计了机侧和网侧变流器协同控制方案。通过MATLAB/Simulink仿真软件搭建了所设计控制方案的仿真模型。仿真结果表明,在电网正常运行的情况下,与传统的矢量控制策略相比,比例-积分-谐振滑模控制策略对输出功率、电流变化的响应更迅速,具有更好的动态性能和抗扰动能力;在电网不对称故障下,能够有效地抑制电磁转矩振荡和直流母线电压波动,提高了双馈式风力发电机的低电压穿越能力。     

单向阀三维动态流场稳定性仿真研究

对单向阀自激振动机理进行研究,采用线性分析方法提出了单向阀的临界稳定曲线。应用泵阀仿真软件PumpLinx的动网格技术,开展了单向阀三维动态流场仿真分析,数值研究了不同工作压力、不同质量流量条件下的单向阀动态稳定特性。结果表明:在某一工作压力条件下,小流量时单向阀工作不稳定,会形成周期性的阀门颤振,随着质量流量增大,单向阀动态稳定性增强;当质量流量足够大时,单向阀工作稳定,单向阀保持稳定开度工作;质量流量一定时,工作压力降低,单向阀动态稳定性增强。单向阀三维动态流场仿真结果与线性稳定性结论相一致,从三维动态流场角度验证了单向阀自激振动机理及线性稳定性分析方法的有效性。     

网格拓扑对DLR-F6构型数值模拟的影响

基于TRIP3.0软件平台和多块结构网格技术,开展了控制方程和网格拓扑等因素对DLR-F6构型数值模拟的影响研究。数值计算采用与试验相同的参数,采用了RANS/TLNS方程,生成了O型和H型2种拓扑的粗、中、细网格进行模拟。分别从气动特性、压力系数分布曲线和表面流态3个方面对结果进行分析。通过与试验数据的对比表明,在模拟小分离流动时,采用H型拓扑网格和RANS方程模型,获得的计算结果会更为准确。     

可压缩流动激波装配在格心型有限体积法中的应用

发展了一种基于格心型有限体积方法(FVM)的激波装配算法。通过定义网格节点属性可以灵活调用激波装配和激波捕捉计算方法。在使用激波装配方法时,激波节点运动速度和下游运动速度通过Rankine-Hugoniot(R-H)关系式获得,同时采用非结构动网格技术描述激波的运动以及调整其他网格节点的位置。流过激波面元的通量为上游单元的基本通量,物理概念更加清晰,通量计算也更为准确。在计算过程中,网格节点属性可以发生变化,以此实现对带有拓扑变化流场的描述。数值试验表明:本文提出的计算方法不但具有较高的计算精度,同时能有效地避免由于捕捉激波而出现的数值问题。     

基于多核计算机集群系统的网格分区策略

网格分区技术是提高并行计算效率的重要手段。基于图的分区技术(如Metis等)已经可以完全做到负载平衡,但是随着计算机集群系统的发展,节点间的通讯成为并行计算中的重要时间消耗部分。为此,提出一种适用于非结构网格的新分区策略(MC Partition method),减少了节点间的通讯规模,结果表明,这种分区方法可以较大的提高并行计算效率。     

亚跨超CFD软件平台研制

亚跨ＣＦＤ软件平台研制的宗旨是综合国内外计算流体动力学的最新研究成果，用计算机软件技术的新发展，建立一个面向对象的计算流体动力学的研究开发环境，加快计算流体动力学最新成果向工程应用的转换速度。本文介绍初级版本的研究进展和研究成果，主要匝数值网格生成技术，流场解算器，面向对象的界面设计等三部分组成。     

Determination of optimal samples for robot calibration based on error similarity

Industrial robots are used for automatic drilling and riveting.The absolute position accuracy of an industrial robot is one of the key performance indexes in aircraft assembly,and can be improved through error compensation to meet aircraft assembly requirements.The achievable accuracy and the difficulty of accuracy compensation implementation are closely related to the choice of sampling points.Therefore,based on the error similarity error compensation method,a method for choosing sampling points on a uniform grid is proposed.A simulation is conducted to analyze the influence of the sample point locations on error compensation.In addition,the grid steps of the sampling points are optimized using a statistical analysis method.The method is used to generate grids and optimize the grid steps of a Kuka KR-210 robot.The experimental results show that the method for planning sampling data can be used to effectively optimize the sampling grid.After error compensation,the position accuracy of the robot meets the position accuracy requirements.     

Numerical study on flow fields and aerodynamics of tilt rotor aircraft in conversion mode based on embedded grid and actuator model

A method combining rotor actuator disk model and embedded grid technique is presented in this paper, aimed at predicting the flow fields and aerodynamic characteristics of tilt rotor aircraft in conversion mode more efficiently and effectively. In this method, rotor's influence is considered in terms of the momentum it impacts to the fluid around it; transformation matrixes among different coordinate systems are deduced to extend actuator method's utility to conversion mode flow fields' calculation. Meanwhile, an embedded grid system is designed, in which grids generated around fuselage and actuator disk are regarded as background grid and minor grid respectively, and a new method is presented for ‘donor searching' and ‘hole cutting' during grid assembling. Based on the above methods, flow fields of tilt rotor aircraft in conversion mode are simulated, with threedimensional Navier–Stokes equations discretized by a second-order upwind finite-volume scheme and an implicit lower–upper symmetric Gauss–Seidel(LU-SGS) time-stepping scheme. Numerical results demonstrate that the proposed CFD method is very effective in simulating the conversion mode flow fields of tilt rotor aircraft.