旋成体射弹倾斜入水运动仿真

为了研究运动参数和弹头外形对弹体斜入水过程的影响规律,采用气液两相流体积分数和水汽空化模型,通过嵌套网格实现刚体三自由度运动学和动力学耦合,模拟了弹体以80～100 m/s速度倾斜入水开空泡阶段的运动过程。经文献实验验证,入水弹体速度与位移的误差为0～6%和-8%～0,转动角度误差为-6%～0。通过对入水速度和入水角度的多工况模拟研究,发现入水速度增大,弹体轴向冲击载荷增大,最大载荷与速度的平方呈线性关系,弹体速度非线性衰减率大;入水角增大,弹体转动角速率减小,运动稳定性强,速度衰减率不受入水角影响。与圆锥头部弹体相比,采用头部阶梯状修型后的弹体的平均速度衰减率、转动角速率和最大轴向冲击载荷分别降低到66.7%、40%和77.2%,显著提高了运动稳定性。     

水下航行体惯性类水动力系数的计算方法

为了系统地计算水下航行体全套惯性类水动力系数,提高计算效率和计算精度,通过基于无粘模型的计算方法对水下航行体的运动进行了预报;通过UDF及动网格技术,对匀速运动和匀加速运动的SUBOFF模型进行了分析.设计了单方向速度线性变化的匀加速直线运动和匀加速回转运动,并通过换算和数值拟合处理得到了潜艇所受的惯性力和惯性类水动力系数.该系数与试验误差保证在6％之内,验证了所提方法的可行性及准确性.     

球锥体高超声速绕流的俯仰阻尼导数的数值计算

采用Navier Stokes方程描述球锥体外形非定常振动流场 ,在Etkin理论下给出俯仰阻尼导数的计算公式。对定常流场的计算采用ADI形式的NND格式 ,对非定常流场的计算采用四步Runge Kutta方法 ,并引入变系数残值光顺技术加速收敛。将俯仰阻尼导数计算结果与实验及工程计算进行了比较 ,并数值研究了其随振动中心的变化规律。     

非定常Euler方程数值计算中高精度格式应用

采用五阶精度加权紧致非线性格式(WCNS)和非定常“双时间步”方法求解非定常Euler方程,模拟NACA0012翼型强迫俯仰振动流场,研究了高精度格式应用到非定常计算时“双时间步”方法物理时间步长、子迭代收敛判据、子迭代步数以及物理时间导数离散方法对计算精度和计算效率的影响.     

双时间隐式方法求解非定常绕流的相关问题

空间采用三阶精度的WeightedNND(WNND)格式,时间方向则采用具有二阶精度的双时间步隐式方法求解二维非定常Euler方程,数值模拟了NACA0012翼型俯仰强迫振动流场,讨论了内迭代步数、内迭代收敛判据以及物理时间步长对计算结果、计算效率与计算CPU时间的影响,并将计算结果与实验值进行了比较。     

固定翼无人机曲线路径跟踪的积分向量场方法

常规的向量场方法在处理无人机曲线路径跟踪问题时很容易受非定常风扰的影响而使得跟踪误差增加,因此很多方法采用用无人机的惯性坐标系(地速和方位角)替代机体坐标系(空速和偏航角)的方式来提高抗风性能。但是,这种方式只能处理大小和方向均恒定的风扰,这在实际飞行中是过于理想的假设。为了克服这些不足,提出了一种采用侧偏距的积分来主动抵消非定常风扰的积分向量场方法用于固定翼无人机曲线路径跟踪控制。根据期望路径的曲率及路径角,结合无人机自身的状态信息设计了曲线路径跟踪策略,并且使用李雅普诺夫理论证明了提出的方法能够确保闭环系统的全局渐进稳定。最后,使用高性能半实物仿真系统验证了提出方法的抗风跟踪性能。     

火箭发动机燃烧室内气体非定常流动数值模拟

应用新的ＰＩＳＯ算法对液体火箭发动机内非定常流动过程进行了数值模拟计算。算法采用一步隐式预测、两步显式校正完成每一时间层计算,而不是通常的多次迭代计算,因而大大缩短计算时间。     

吸气式高超声速飞行器助推分离过程数值仿真

针对吸气式高超声速冲压发动机验证性试验特殊的飞行环境和助推分离条件，以某轴对称吸气式高超声速飞行器级间分离问题为具体研究对象，采用非结构网格局部网格重构技术和非定常问题非定常六自由度问题仿真方法，对该复杂构型飞行器助推分离过程进行数值计算。研究得到弱干扰冷态分离状态下飞行器及助推器的运动参数和气动力参数在分离过程中的发展规律。对0.3 s内助推器的位移轨迹进行分析，判断分离方案的可行性，并给出最佳的分离工况条件。     

非定常超音速流动数值计算方法比较

为适用于非定常流场计算 ,在保证时、空二阶精度的前提下 ,通过引入变系数残值光顺技术 ,提高了四步Runge Kutta方法的计算效率 ;并针对球锥简谐振动过程 ,数值比较Runge Kutta多步法和“双时间法”的计算精度。     

非牛顿流体不定常环管流研究

本文采用积分变换的方法,找到了一类非牛顿流体在环形管道中不定常流动的解析解,并进行了数值计算,分析了非牛顿性系数和其他各参数对二阶流体不定常流动的影响,指出当二阶流体非牛顿系数相同时,环管流与一般管流比达到稳定的特征时间较短,并且相应的速度分布、平均速度分布数值均较小。在外半径相同时,环管流内壁的剪应力较之一般管流,其大小随内径而改变。环管流外壁剪应力总相应地小于内壁剪应力。