锥头圆柱体高速入水空泡数值模拟

采用流体体积法（VOF,Volume of Fluid）多相流模型,引入动网格技术,对锥头圆柱体垂直自由高速入水问题开展了数值模拟研究,分析了航行体入水后速度及入水空泡形态的发展规律,并与文献基于能量守恒定律的理论结果进行比较,两者吻合较好.研究结果表明,高速入水初期,结构受到极高的冲击载荷,导致速度迅速衰减;随着入水深度的增加,空泡不断拉长并向径向扩张,空泡内空化现象明显,空泡内混合相对空泡壁及空泡分离点附近压力分布影响显著.对于不同速度入水的研究结果表明,入水速度越高,航行体头部压力峰值越大,相同时刻入水空泡最大直径也越大.      

圆柱体垂直入水空泡形态试验

在不同攻角和不同速度条件下,对圆柱形运动体垂直入水空泡形态进行试验研究.对140°锥角头型圆柱体在零攻角和带攻角条件下进行垂直入水试验,分析了低速垂直入水所引起的入水空泡生成、发展和闭合流动过程以及空泡稳定性等;利用平头圆柱体进行不同速度条件下垂直入水试验,分析了入水速度对空泡闭合方式和形态的影响.研究结果表明:在试验速度范围内均能形成稳定的入水空泡,并经历生成、充分发展、空泡闭合以及最后完全脱落的过程,不同发展阶段具有不同的稳定性,其闭合方式和形态与入水速度密切相关.     

开放空腔壳体入水流场结构及流体动力特征研究

基于混合介质雷诺平均Navier-Stokes方程,对开放空腔壳体垂直入水运动过程开展了数值研究,得到了压力场、速度场分布,空泡波动、闭合特征,空腔气体涨缩规律,以及流体动力变化规律,并分析了空腔结构在入水运动过程中对流场结构和流体动力的影响。结果表明:液体随气体涨缩同步进出开放端;开放端局部形成波动的压力源和周期性的压力场、速度场分布;入水空泡呈现波动形态,其扩展程度与开放端液体流速相关;空泡内形成气体漩涡,随空腔涨缩往返进出空泡,对空泡闭合具有抑制作用;流体动力呈波动变化规律,频率与气体涨缩频率一致,幅值与气体涨缩程度成正比。开放空腔结构在入水过程中空腔内气体发生涨缩运动,对流场结构和流体动力产生周期性扰动作用,在一定程度上可以减缓冲击、维持空泡及运动的稳定性。      

空气域压力对高速射弹入水流场影响

为研究空气域压力对射弹入水流场的影响,采用流体体积法(VOF)多相流模型对锥头圆柱体高速入水过程进行数值模拟,得到射弹入水速度与入水深度以及空泡形态的分析结果,并将数值模拟结果与基于牛顿第二定律和文献的理论结果进行对比,验证数值方法的正确性.基于该方法针对不同空气域压力条件对入水过程进行分析,结果表明空气域压力对自由液面上方的喷溅形态、空泡表面闭合时间以及空泡内部空化效应影响较大.空气域压力越大,空泡敞开阶段的空泡口处半径越小,空泡发生表面闭合的时间越早,空泡内部空化效应越强.     

通气位置对潜射航行体流体动力特性影响分析

基于均质平衡流理论,通过求解混合介质的RANS方程、SST湍流输运方程和各相之间的质量输运方程,开展了通气位置对潜射航行体流体动力特性影响的三维数值模拟研究,对比分析了不同通气位置条件下空泡形态特性、表面压力分布以及阻力变化特性.结果表明:当两个通气口间距增加到一定值后,在两个通气口之间区域空泡发生断裂,在断裂区域迎流面和背流面表面压力值升高,且在空泡断裂闭合位置出现压力峰值;在相同通气量条件下,随着通气口间距的增大,压差阻力系数和粘性阻力系数均呈减小趋势变化.     

空泡摆动对超空泡航行体尾拍影响分析

针对超空泡航行体存在尾拍运动与空泡形态变化互相耦合作用的特点,建立了超空泡航行体尾拍耦合运动方程,并对超空泡摆动对超空泡航行体尾拍相互作用过程进行了数值研究,对比分析了不同速度和角速度的超空泡航行体运动过程及空泡摆动对尾拍的影响规律.研究结果表明:空泡摆动使航行体转角、整体运动周期、尾拍升力和阻力增大;在不同速度下,空泡摆动对航行体转角、整体运动周期及尾拍周期的影响随速度增大而逐渐减小;不同初始角速度下,空泡摆动对航行体转角、整体运动周期、航行体尾拍升力和阻力的影响随角速度增大而逐渐增大.     

航行器低速斜入水运动规律

针对介质跨越航行器控制困难的问题,提出一种空中控制水下非控的单一控制策略;为了分析航行器非控状态下斜入水运动的规律,构建了航行器低速入水动力学模型,并分别使用数值仿真方法和理论模型求解方法进行同一条件下的航行体入水运动仿真,通过对2种方法的仿真结果对比验证本文所构建航行体斜入水动力学模型的正确性。利用构建的入水动力学模型,分别对不同初始速度、角度、攻角条件下的入水过程进行了运动状态仿真并分析,得出了航行体在入水运动过程中的姿态位置变化规律。此入水规律将指导介质跨越航行器后续的水下航行、进而出水的一系列研究。      

几何和流动参数对上翘后体阻力系数的影响

利用N-S方程对9个上翘后体模型进行了气动力计算,主要研究后体几何参数和流动参数对上翘后体阻力系数的影响.研究结果表明, 后体的压差阻力系数分别随上翘角、收缩比的增加及迎角的减小而明显增加;后体摩擦阻力系数分别随后体的长细比的增加和雷诺数的减小而增加;后体越扁平,其压差阻力系数越大;在跨音速时,波阻系数也与上翘角有关,上翘角增加会导致波阻系数进一步加大.     

吹风比和曲率对旋转曲面气膜冷却效率影响

采用数值模拟方法,对旋转状态下曲率表面的气膜冷却进行研究.通过不同的曲率半径和吹风比,得到了冷却效率的分布情况,从而获得吹风比和曲率两因素对气膜冷却效率的影响规律.研究结果表明:低吹风比下,冷却效率较好;高吹风比下,气膜容易脱离壁面;凸表面的冷却效率随着曲率半径的增加而逐渐减小;而凹表面的冷却效率随着曲率半径的增加而逐渐增大;凸表面上曲率的影响作用随着旋转数的增加而逐渐弱化;而凹表面上曲率的作用随着旋转数的增大而逐渐增强.     

上升气泡与塑料平板在纯水中的碰撞黏附行为

为了揭示气泡尺寸和表面疏水性对塑料浮选过程的影响规律,采用高速摄影技术观测上升气泡与塑料平板在纯水介质中的碰撞黏附行为。根据气泡速率的变化情况,将碰撞黏附行为细分为碰撞、液膜排液、三相接触线扩散等3个阶段,基于图像处理技术定量分析了气泡碰撞形变、黏附时间、三相接触线扩散特性及其影响因素。结果表明:变形因子随气泡直径的增大由略大于1.00逐渐演化为略小于1.00,疏水性强的塑料平板使相同直径气泡的变形因子更大,碰撞时间随气泡直径的增大而增加,随疏水性的增加而减小。塑料平板疏水性越强,形成三相接触所用时间越短,当气泡直径为1.0 mm时,聚四氟乙烯（PTFE）和聚甲基丙烯酸甲酯（PMMA）平板的液膜排液时间均出现极小值,分别为4.8 ms和56 ms。三相接触线扩散时间随气泡直径的增大而增大,随塑料平板疏水性的增大而减小,在气泡直径大小相同时,三相接触线扩散直径随塑料平板疏水性的减小而减小。      

飞行器动态下俯过程中的负阻力现象

利用水洞中动态测力/流动显示一体化实验系统,研究飞行器简易模型等速俯仰过程动态气动特性.测力天平为三分量内式尾支形式;利用氢气泡方法观测机翼上表面流动,测力和流动显示结果同时记录.当下俯无量纲角速度k≥0.11时,观测到正攻角时的负阻力现象.分析测力和流动显示结果,表明负阻力是下俯过程中上翼面流动由分离向附着恢复过程中的迟滞效应引起的,负阻力的大小与下俯无量纲角速度和下俯起始攻角有直接关系.      

操纵面加装Gurney襟翼对无人机纵向气动特性的影响

通过风洞测力实验研究了矩形Gurney襟翼对多操纵面飞翼布局无人机纵向气动特性的影响.升降副翼和襟副翼加装Gurney襟翼都会增加全机升力系数以及低头力矩系数.通过比较发现,升降副翼加装Gurney襟翼的增升效果好于襟副翼,外侧操纵面控制效果好于内侧.升降副翼和襟副翼分别组合加装Gurney襟翼会进一步提高增升效果,并且升力系数随着Gurney襟翼高度的增加而增加.鸭翼上加装Gurney襟翼可以减小全机低头力矩系数,因此鸭翼和升降副翼、襟副翼组合加装Gurney襟翼在线性段提高全机升力系数的同时,可以保持低头力矩系数增量基本为0,该特性对改善飞行器性能尤为有利.     

翼伞系统纵向飞行性能分析

翼伞系统的飞行性能不仅取决于翼伞本身的气动特性,而且与安装角、伞绳长度、回收物阻力特征、翼载荷等系统参数密切相关。文章应用拉格朗日方程建立翼伞系统的纵向飞行力学模型,对翼伞系统进行飞行力学数值仿真,深入分析了系统参数以及开伞状态对翼伞系统纵向飞行性能的影响规律。结果表明：只有安装角在0°～20°时,翼伞系统才能达到稳定的滑翔状态,且安装角在4°～6°时对应两个稳定的滑翔状态,具体由开伞姿态和速度决定;伞绳特征长度的增加使系统的静稳定性增加;回收物的阻力特征增加6m2,翼伞系统的稳定滑翔角增加15°左右,而迎角减小不到1°;翼伞飞行速度随着翼载荷的增加而增加,其平方与回收物质量成正比。上述结论可为翼伞系统的工程实际应用提供参考。     

南京地区斜向加热电离层效应数值模拟

与垂直加热相比,斜向加热电离层具有灵活性高、影响范围广和便于实际操作等优势.在非偏区考虑欧姆吸收,基于电子能量方程和连续性方程构建电波斜向加热低电离层的物理模型,并利用IRI-2007和NRLMSISE-00经验模型提供的背景参数对南京地区斜向加热低电离层进行数值模拟,对比不同加热条件下电子温度和电子密度的扰动情况.研究结果表明:电波加热效果随入射仰角和有效辐射功率的增大而增大;电子温度和电子密度增幅随电波频率增大而减小;X波模比O波模造成的电子温度扰动幅度和电子密度扰动幅度更大,同时X波模比O波模能更快地使电子温度和电子密度达到稳定状态;一定范围内较小仰角、较低频率、较大有效辐射功率的电波能使电子密度更快达到稳定,后两者还能加快电子温度达到稳定的过程;电子温度达到稳定所需时间随入射电波仰角呈单峰变化,仰角为62°时达到最大.      

Command generator tracker based direct model reference adaptive tracking guidance for Mars atmospheric entry

In order to accurately deliver an entry vehicle through the Martian atmosphere to the prescribed parachute deployment point, active Mars entry guidance is essential. This paper addresses the issue of Mars atmospheric entry guidance using the command generator tracker (CGT) based direct model reference adaptive control to reduce the adverse effect of the bounded uncertainties on atmospheric density and aerodynamic coefficients. Firstly, the nominal drag acceleration profile meeting a variety of constraints is planned off-line in the longitudinal plane as the reference model to track. Then, the CGT based direct model reference adaptive controller and the feed-forward compensator are designed to robustly track the aforementioned reference drag acceleration profile and to effectively reduce the downrange error. Afterwards, the heading alignment logic is adopted in the lateral plane to reduce the crossrange error. Finally, the validity of the guidance algorithm proposed in this paper is confirmed by Monte Carlo simulation analysis.     

联翼布局俯仰力矩非线性变化特性的数值模拟

联翼布局飞机具有优良的升阻特性,是下一代亚声速飞机优先选择的气动布局型式之一,但在某些情况下其俯仰力矩随迎角的增大会表现出较强的非线性变化特点.针对该问题,在Ma=0.75条件下,采用数值模拟方法对某亚声速联翼布局气动性能及其绕流流场进行研究,通过对各部件气动特性分析,结合不同前翼绕流流动状态下前/后翼绕流场特点及截面气动力分布特点,揭示了前翼对后翼绕流流场干扰是引起其俯仰力矩非线性变化的主要原因.计算结果表明:在一定迎角下,该联翼布局飞机前翼绕流发生分离,从而影响后翼绕流流场,引起后翼气动效率下降,导致全机俯仰力矩随飞机迎角发生非线性上扬,对该机飞行性能的提高带来严重影响.     

火星大气进入段抗饱和固定时间阻力加速度跟踪制导律设计

研究了具有固定时间收敛特性的火星探测器大气进入段的标称轨迹跟踪制导问题。首先，针对横向运动，给出与速度成线性关系的航向误差漏斗走廊形式，完成了倾侧角的反号逻辑设计。与横程漏斗走廊反号逻辑相比，该逻辑计算量小，更适用于宇航计算机。与航向误差宽度走廊反号逻辑相比，该逻辑在高速状态下能够避免倾侧角的频繁切换，可提高任务成功的概率。其次，针对纵向运动，通过RBF神经网络补偿了倾侧角饱和问题，利用积分滑模设计了阻力加速度固定时间饱和跟踪制导律，其不仅可有效消除滑模控制的抖振问题，且将跟踪误差以两种不同形式引入制导律，能够加速收敛，能够保证跟踪误差在固定时间内快速收敛至0。最后，通过数值仿真验证了所设计的横向倾侧角切换逻辑和纵向制导律对标称轨迹的快速、精确跟踪能力。     

开裂式阻力方向舵流固耦合机理分析

为了探究阻力方向舵开裂状态下的流场形态和流固耦合运动机理,采用计算流体力学(CFD)方法开展了不同开裂角下的二维阻力方向舵的流场计算。基于动力学模态分解(DMD)方法对各流场进行模态分解,分析了各模态的流动特征及频率变化。结果表明,在20°开裂角的范围内,机翼绕流的流场结构以开裂区内的驻涡及后缘脱落涡为主,流场各阶模态频率随来流速度的增大而增大,随开裂角的增大而减小。同时,对不同开裂角的二维翼型开展了流固耦合计算。结果表明, 随着折减速度的增加,系统的流固耦合运动形式由涡致振动发展为流动失稳,系统的失稳边界随着开裂角的增大而提高。