二维平板水漂运动数值模拟

探索平板水漂运动的流体力学现象机理对飞行器着水问题的研究有重要的参考价值。基于有限体积法和k-ε RNG湍流模型求解非定常雷诺平均Navier-Stokes（URANS）方程,采用流体体积分数（VOF）模型与速度入口造波法构造数值波浪水池,并结合整体动网格方法,完成二维平板在静水面与波浪水面的水漂运动数值模拟。在与试验值和理论值对比的基础上,讨论初始姿态角、投掷角度与投掷速度对水漂运动的影响。进一步研究不同波浪参数和波浪位置对平板水漂的影响,并从能量角度展开分析。结果表明：20°姿态角平板能以最小的投掷速度实现水漂运动;运动中平板能量的相对损失率受初始投掷速度的影响较小,主要受投掷角和姿态角作用,随投掷角或姿态角的增加而增大。波浪情况下,在上行波位置触水的平板能够获得更大的接触面积,而更长的触水时间发生在波谷处;因此,这两个位置触水的平板在水漂运动过程中的相对能量损失大,其数值变化比波峰位置大5%左右;在上行波位置触水时,平板的速度衰减与相对能量损失随着波高的增加而增大,而波峰位置则有相反的变化趋势。     

着水初始条件对水陆两栖飞机着水性能的影响

针对着水初始条件对水陆两栖飞机静水面降落性能影响的问题,使用非定常雷诺平均Navier-Stokes方程和标准k-ω两方程模型求解非定常流场,耦合流体体积模型捕捉自由液面,结合改进的流场域网格划分策略,并借助任意拉格朗日欧拉方法处理机体、网格和水气交界面三者关系。基于上述方法,探讨初始姿态角、下降速度和前飞速度对着水性能的影响,其中驾驶舱与重心位置加速度、气动力与水动力、俯仰角与吃水深度等参数变化被考虑在内,逐步得到最优着水初始条件。结果表明：增大初始姿态角有利于减缓飞机各位置除初次触水过载峰值外过载峰值的变化。当飞机俯仰力矩出现极值时,驾驶舱垂向过载显著大于其他方向和位置过载。减小初始下降速度,各加速度值也随之减小,值得注意的是重心垂向加速度值与初始垂向下降速度平方值存在线性变化关系。着水前飞速度的选取应使得气动升力略大于重力,飞机才能获得更佳的着水性能。     

楔形体入波浪水面数值模拟

数值模拟分析楔形体入波浪水面流体力学现象和运动姿态变化历程，为水上飞机在波浪水面起降提供技术和理论支持。基于有限体积法求解非定常雷诺平均Navier-Stokes（URANS）方程和SST k-ω湍流模型，结合整体动网格法的速度入口造波方法、阻尼消波法和流体体积（VOF）液面捕捉技术构造数值波浪水池，耦合三/六自由度模型，实现楔形体入线性规则波和不规则波的数值模拟。在楔形体入静水面算例验证的基础上，详细研究二维楔形体在线性规则波波峰、波谷、平衡位置-上行速度和下行速度处入水过程的受力特性和运动姿态变化。结果表明：速度入口造波方法的波浪数值解与理论解析解较吻合，偏差为1%；在二维楔形体入规则波过程中，楔形体所受垂向力和垂向速度变化趋势相同，对横向位移影响较小，其数值变化少于0.01；在平衡位置入水过程中，模型的滚转角和横向速度变化明显，其数值变化约为波峰、波谷位置处的8倍和10倍。分析原因：平衡位置处波浪内流速度变化较快，模型所受波浪力冲量较大，且楔形体两侧斜边与波浪的相互作用力也参与到模型的运动中。数值模拟二维楔形体在线性不规则波5个随机时间点入水及三维楔形体在波峰、平衡位置-上行速度处入水过程，其变化规律及形成原因与入规则波相似。     

高速平板着水数值模拟

探索和揭示物体入水冲击的流体力学现象与机理对飞行器水上迫降问题的研究有重要的参考价值。对高速平板着水涉及到的复杂物理问题展开数值模拟,采用有限体积法求解非定常雷诺平均Navier-Stokes(URANS)方程和标准k-ε湍流模型,流体体积(VOF)模型捕捉水气交界面,整体动网格技术处理平板与水面的相对运动。在二维楔形体入水冲击的算例验证基础上,详细研究平板高速着水引起流体喷溅、射流、空气垫等现象和平板底面压力变化历程,结果表明:空气垫现象明显,俯仰角4°平板下表面出现规律的空气泡,10°时则不存在;平板下表面的水体沿壁面运动,当俯仰角为10°时,壁面水体的运动速度显著增加;在大俯仰角的情况下明显出现负压区。