多体分离非定常气动特性数值模拟

发展建立了用于研究多体分离的非定常流场数值模拟方法,该方法采用网格弹性变形的结构网格动网格技术,采用有限体积法求解ALE(Arbitrary Lagrangian-Eulerian)描述形式的三维非定常Euler方程,与刚体运动方程耦合,数值模拟了子弹穿越母弹头部激波的子母弹分离过程,并与采用攻角修正的准定常计算结果进行了对比分析,表明多数情况下两种数值方法得到的结果比较接近。     

薄片分离过程流场特性的数值仿真研究

利用弹簧近似法和局部网格重构相结合的非结构动网格方法,耦合求解Navier-Stokes方程和六自由度刚体动力学方程,数值模拟了三维薄片非定常分离过程。数值仿真结果有效反映出薄片分离过程中运动速度、轨迹、受力变化规律及流场特性的发展与变化。根据数值仿真结果对薄片分离机理、多体气动干扰作用进行了分析。分离过程中单个薄片与主体段周围形成复杂流动,并伴随着尾迹涡的运动、发展等非定常过程。气动载荷的非线性变化是薄片分离的主要原因。     

基于预估校正和嵌套网格的虚拟飞行数值模拟

 针对导弹虚拟飞行数值模拟问题,发展了空气动力学/飞行力学数值计算方法和软件。控制方程为非定常雷诺时均Navier-Stoker(RANS)方程和刚体六自由度运动方程;流场求解器为有限体积法结构网格求解器,时间推进采用双时间步法,湍流模型为Spalart-Allmaras一方程模型;采用Adams预估校正法实现飞行力学方程与流场控制方程的耦合计算;使用嵌套网格方法模拟多体运动。首先模拟了美国国家航空航天局(NASA)窄条翼导弹模型纵向虚拟飞行,研究耦合方式和时间步长的影响。仿真结果表明,双时间步三阶Adams耦合方法,同等精度下可以显著增大时间步长,缩短仿真时间。最后,采用该方法模拟了导弹自由摇滚特性和纵向虚拟飞行,模拟结果与试验值吻合较好。     

多体全自由飞过程的数值模拟

以三维非结构粘性直角网格重构为构架,数值求解了固定空间网格瞬间运动边界的N-S方程,逼近了流场的非定常过程,结合六自由度的弹道方程,对单体或多体运动过程进行了全自由飞数值模拟,数值模拟结果和级间分离的风洞自由飞试验进行了对比,初步的计算结果说明,全自由飞数值模拟在一定程度上可以描述风洞试验的飞行轨迹,并且适应能力更强,整个全自由飞数值模拟系统可以为更复杂的单体机动飞行或多体运动过程提供快速实用的数值模拟方法.     

基于数值虚拟飞行技术的飞行器动态特性分析

基于结合结构重叠网格、闭环PID控制器、舵偏控制律、刚体六自由度运动和非定常Navier-Stokes方程求解等模块的数值虚拟飞行技术,对“起源号”返回舱、基本带翼导弹外形的多自由度非定常运动、受控特性及控制参数的整定开展了模拟。分析了不同自由度（DOF）下飞行器的运动特性,飞行器受扰动后的稳定性及控制参数的整定。计算结果表明：利用数值虚拟飞行技术可有效地开展复杂流动下飞行物体非线性运动问题的研究,对研究飞行器在非线性流动下的动态特性、受控特性、流动机理研究以及控制律的设计检验具有较高的工程价值和应用前景。     

用非结构动网格方法模拟有相对运动的多体绕流

发展了用于模拟有相对运动的多体非定常绕流，并计算由流场决定的外弹道的数值方法。该方法采用非结构动网格，显式NND有限体 格式，求解非定常Euler方程，并耦合了刚体动力学方程，作为例证，计算了俯仰振动的NACA0012翼型绕流，计算结果与实验及文献结果吻合，最后，整个策略用于计算二维机翼－外挂物分离流场，并确定了外挂物下落轨迹。     

基于闭环控制的带翼导弹虚拟飞行数值模拟

通过耦合求解非定常Navier-Stokes方程、刚体六自由度（6DOF）运动方程及舵偏控制律,发展了一套虚拟飞行数值模拟平台,对基本带翼导弹外形进行了虚拟飞行数值模拟研究。流场求解采用基于结构网格的有限体积法,时间推进采用双时间步法,使用刚性重叠动网格技术模拟舵与弹身的相对运动,采用四阶Runger-Kutta方法求解运动方程组,舵偏控制律采用基于比例-积分-微分（PID）的反馈控制。模拟结果表明：发展的虚拟飞行数值模拟平台能够处理包含复杂运动边界的非定常运动问题,可以为非定常气动力建模提供数据、预示飞行器稳定性及检验飞行控制律,该模拟平台具备一定的工程应用价值。     

考虑喷流效应的机载导弹发射及气动干扰数值模拟

研究比较了机载导弹发射时是否考虑导弹尾喷流情况下导弹运动特性及机-弹气动干扰影响.通过动态结构嵌套网格技术实现了导弹运动的模拟,结合刚体六自由度运动方程求解三维非定常Euler方程得到流场信息.导弹推力由发动机燃烧室喷流计算得到;当不考虑喷流效应时,推力通过直接在导弹尾部给定非定常作用力实现.应用此方法模拟了基于类"全...     

基于CFD/RBD方法的旋转弹轨迹仿真

旋转弹在飞行过程中,由于表面气动载荷分布不对称,可能会出现耦合共振、灾难偏航等现象,使导弹飞行失常偏离目标。针对此类问题,采用非定常N-S方程和刚体动力学方程耦合求解的方法,仿真旋转弹飞行轨迹。非定常气动力计算使用非结构混合网格N-S方程求解程序HUNS3D,刚体六自由度运动方程采用改进的四阶Adams预估校正法求解。将这种耦合求解方法应用于美国ARL旋转弹试验模型的轨迹仿真,并将计算得到的旋转弹姿态角、空间位置和气动力与相关文献进行对比分析,结果表明:本文方法能较精确地仿真旋转弹的飞行轨迹。     

安装角度对内埋式导弹分离特性的影响

基于弹簧变形理论和局部重构的非结构动网格技术,采用有限体积法求解ALE描述下的三维非定常Euler方程,并耦合刚体六自由度运动方程,模拟了超声速条件下三种不同安装状态导弹从弹舱中弹射分离的过程,分析了导弹质心相对位移以及姿态变化情况,发现增大初始安装角度能够改善导弹在超声速环境下的分离特性。     

基于非结构动网格的非定常流数值模拟方法

给出了基于非结构动网格的非定常流动问题的数值模拟方法,重点说明两种动态非结构网格方法的原理以及在非定常流模拟中的算法.一种是基于非结构网格的重叠网格方法,该方法结合了非结构网格方法和重叠网格方法的优点,以物面距为准则参数对网格节点进行分类,进而建立重叠网格的网格问插值关系,通过采用动态非结构重叠网格算法,可以处理多体间...     

A CFD-based numerical virtual flight simulator and its application in control law design of a maneuverable missile model

A CFD-based Numerical Virtual Flight (NVF) simulator is presented, which integrates an unsteady flow solver on moving hybrid grids, a Rigid-Body Dynamics (RBD) solver and a module of the Flight Control System (FCS). A technique of dynamic hybrid grids is developed to control the active control surfaces with body morphing, with a technique of parallel unstructured dynamic overlapping grids generating proper moving grids over the deflecting control surfaces (e.g. the afterbody rudders of a missile). For the flow/kinematic coupled problems, the 6 Degree-Of-Freedom (DOF) equations are solved by an explicit or implicit method coupled with the URANS CFD solver. The module of the control law is explicitly coupled into the NVF simulator and then improved by the simulation of the pitching maneuver process of a maneuverable missile model. A nonlinear dynamic inversion method is then implemented to design the control law for the pitching process of the maneuverable missile model. Simulations and analysis of the pitching maneuver process are carried out by the NVF simulator to improve the flight control law. Higher control response performance is obtained by adjusting the gain factors and adding an integrator into the control loop.     

旋翼/机身干扰非定常流场数值模拟(英文)

基于非结构动态嵌套网格方法对直升机旋翼/机身干扰非定常流场进行了三维Euler方程数值模拟。嵌套网格系统由包含旋翼的运动区网格和包含机身的静止区网格组成,两区网格均来自于同一套初始的整体非结构网格,人工边界以及初始的插值贡献单元对应关系在网格分区时同时确定。在非定常计算中两区网格始终存在一定的重叠区域,无需重新建立人工边界。高效的贡献单元搜索算法可以根据初始的插值贡献单元对应关系很快更新每一个真实时间步上的贡献单元对应关系。为了准确计入旋翼桨叶的挥舞和变距变形,在运动区网格中引入线性弹簧方法。采用非定常Euler方程和双时间推进方法模拟了美国佐治亚理工学院旋翼/机身干扰实验外形非定常流场,计算结果与实验值以及文献的计算值相吻合。     

共轴式双旋翼悬停流场和气动力的CFD计算

考虑到共轴式双旋翼流场特征高度复杂、桨叶承受非定常气动载荷的特点,为更好地预测共轴式双旋翼的气动特性,把非结构嵌套网格方法和网格的自适应技术相结合,发展了一套适合于共轴式双旋翼流场数值计算的求解器.在该求解器中,采用非结构嵌套网格方案来模拟桨叶之间存在的相对运动,自适应网格技术用来捕捉尾迹对流场和气动特性的影响,求解惯性坐标系下的非定常N-S主控方程来模拟流场的非定常特性.应用该求解器,首先计算了有试验结果可供对比的一副试验旋翼的诱导速度场分布,在此基础上,计算了共轴式双旋翼的桨尖涡轨迹和拉力分布特性,并与单旋翼的计算结果进行了对比和分析,得到了一些有意义的结论.     

A novel high-order scheme for numerical simulation of wake flow over helicopter rotors in hover

Accurate prediction of tip vortices is crucial for predicting the hovering performance of a helicopter rotor. A new high-order scheme(we call it WENO-K) proposed by our research group is employed to minimize numerical dissipation and extended to numerical simulation of unsteady compressible viscous flows dominated by tip vortices over hovering rotors. WENO-K is referred to as an adaptively optimized WENO scheme with Gauss-Kriging reconstruction, and its advantage is to reduce dissipation in smoo...     

Application of immersed boundary method in turbomachines

Simulating unsteady turbulent flow in turbomachines is still challenging due to the complexity of blade geometry and relative motion between rotor and stator. This study presents an Immersed Boundary Method (IBM) for high-Reynolds turbomachinery internal flows, and shows the advantage of the automatic grid generation techniques and flexible moving boundary treatments. The wall functions are used in the present method to alleviate the wall resolution restriction of turbulence simulation. The Two-Dimensional (2-D) IBM solver, which was previously developed and tested for a low-speed compressor, is further validated for a well-documented Low-Pressure Turbine (LPT) cascade. Both the blade loading and the total pressure losses in the wake are well captured by the present 2-D solver. The complex Three-Dimensional (3-D) effects in turbomachines motivate the further development of an extended 3-D IBM solver by using a curvilinear-coordinate system that facilitates the hub and casing boundary treatment. The good performance of the 3-D solver is demonstrated through comparison with CFX solver solutions for the rotor configuration of Advanced Noise Control Fan (ANCF). Further effects of the grid resolution on capturing the blade wake are discussed. The results indicate that the present 3-D solver is capable of reproducing the evolution of the blade wake with suitable computational grid.     

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

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