细长锥型射弹不可压缩超空泡流的数值计算方法

假定理想不可压缩流体呈无旋定常流动、起空泡尾部Riabushinsky为闭合方式,基于水动力学势漉理论及细长体理论.建立了描述水下细长锥型射弹超空泡流动的积分微分方程.谊方法发展了求解谊类方程的数值离散方法,提出了多种超空泡外形初值解.分析了它们对超空泡形态计算结果的影响,优化了计算过程,简化了初始遮代条件.计算得到的...     

超空泡航行体锥形空化器优化设计

分析比较超空泡外形计算的不同经验公式,结果表明Guzevsky模型具有较好的精度与适用范围。建立了以减小阻力、增大空化数、增加舱内容积为优化目标的超空泡航行体锥形空化器设计模型,以权重因子法和遗传算法求解了优化问题,验证性计算结果与已有文献吻合。研究了空化数、航行体长径比等因素对超空泡航行体空化器设计的影响规律,结果表明空化数是影响航行阻力的首要因素,空化数决定了超空泡长径比并使之与航行体外形匹配,通过优化设计能够在不增加阻力的情况下有效提升航行体舱内容积。     

势流理论的一种数值方法在空泡形态中的应用

基于势流理论的假设,采取常规边界元法对具有自由面的定常平板和轴对称空泡流动数值方法进行了研究,运用Riabourichinsky镜像空化模型,建立了一种空泡自由面求解的优化迭代方法,以自由流线长度为参数求解模拟了定常无重工况下流经二维平板和轴对称圆盘的空化流动.数值计算结果与解析解的相对误差不超过3%,证实了轴对称体空泡外形的近似回转椭球体假说在空化器附近是不成立的.计算结果与同类文献数据的对比显示了边界元法进行空化流场计算具有着显著的优越性,所采用的迭代方法精度高、收敛速度快,所得结论对深入了解空化特性具有参考价值.     

复杂外形导弹出水过程空化流数值计算研究

基于Zwart-Gerber-Belamri空化模型,采用动网格和复合网格技术,提出了一种复杂外形潜射导弹出水过程中空化流数值计算方法。对头肩部、舵面及突起物的空化生成、演化及影响因素进行了分析。结果表明,出水过程初期空化数降低使附体空泡迅速扩大,随后空泡脱落与振荡的影响将更为显著;肩部在15m/s出水时产生空化并随速度上升而迅速增长, - 攻角时肩空化非对称性显现;舵面空化受出水速度影响较小, 攻角以下时空化面积小于全舵面积的10%,但攻角超过 后空化面积迅速增长至50%以上;减小突起物尺寸有利于避免空化产生,突起物空化受攻角影响较小,但在高速下可诱导临近弹体产生空化。     

高速水下射弹空泡形态及其CFD分析

射弹在水下空泡形态高速运动,实验研究存在很多限制。本文应用两种方法对高速水下射弹的空泡形态及阻力性质进行了研究,一是应用CFD仿真模拟软件,二是应用基于空泡截面独立扩张原理所建立的公式。研究表明：仿真计算结果与公式计算结果吻合很好;随着空化数的减小,空泡的相对长度、相对直径和长细比都在增大;在小空化数下,空泡前部形态基本不变,阻力系数主要为压差阻力系数,其值基本不变。研究结果为射弹的外形设计提供参考。     

潜射导弹出水过程空化流数值计算

运用Zwart-Gerber-Belamri空化模型,采用动网格和复合网格技术,提出复杂外形潜射导弹出水过程中空化流数值计算方法。对头肩部、舵面及突起物的空化生成、演化及影响因素进行分析。仿真结果表明:出水过程初期空化数降低使附体空泡迅速扩大,随后对空泡脱落与振荡的影响将更为显著;肩部以15 m/s的速度出水时产生空化并随速度上升而迅速增长,5°~10°攻角时肩空化非对称性显现;舵面空化受出水速度影响较小,攻角小于等于5°时空化面积小于全舵面积的10%,但攻角大于5°后空化面积迅速增长至50%以上;减小突起物尺寸有利于避免空化产生,突起物空化受攻角影响较小,但在高速下可诱导临近弹体产生空化。     

KP458螺旋桨水动力及空泡特性研究

以螺旋桨(KP458)为研究对象,研究其空泡流场特性。首先通过梢涡区域划分及网格加密,对螺旋桨(KP458)在无空化均匀流条件下进行敞水性能数值模拟,建立螺旋桨流场数值预报模型,并对4套网格进行不确定度分析。将计算结果与实验数据进行对比,然后应用Schnerr-Sauer空化模型数值模拟。通过将试验结果与螺旋桨(VP1304)水动力及空泡特性仿真结果进行对比,无空化流场与空化流场中的推力、扭矩系数误差均低于8%左右,以此验证数值计算模型的可行性。研究表明：采用的数值模型具有良好的精度,螺旋网格加密能够捕捉螺旋桨产生的梢涡;KP458螺旋桨在进速系数较低时有较为明显的梢涡空泡产生,随着进速系数增大空泡面积逐渐减小且梢涡空泡消失,螺旋桨的梢涡强度逐渐减小且梢涡融合现象发生越早。     

轴对称物体空泡流动的数值计算

基于势流理论和边界元方法,对水下运动的回转体的空泡流进行了数值模拟.将半经验空泡封闭模型由平面流动推广至轴对称流动.对头部为圆锥的圆柱体空泡流动进行了数值计算.分析了空化数、半锥角、空泡长度间的关系.计算与实验结果符合较好.     

无尾舵布局超空泡航行体弹道优化设计

给出了无尾舵布局超空泡航行体在纵垂面内的受力分析及运动方程,建立了变深度机动最优轨迹问题的数学模型。针对空泡扩张与收缩带来的延时效应,提出改进的延时型高斯伪谱法,利用所有高斯节点的运动状态插值获得整个运动过程的状态量,由此得到各个时刻的实时超空泡航行微分方程。通过计算各节点处运动方程约束,对其变深度机动轨迹进行数值求解,得到相应的优化弹道结果及控制策略。对给定布局方案的超空泡航行体下潜机动最优弹道问题进行求解,计算结果表明,采用空化器偏转控制策略,无尾舵布局超空泡航行体可满足运动稳定性与机动性要求,所提出的延时型高斯伪谱法能够处理超空泡航行弹道优化中的延时效应,具有较好的快速收敛性。     

高速射弹超空泡的形态特性

为了掌握高速射弹超空泡的形态特性,利用某口径火炮进行了射击试验,射弹最大初速达900m/s。获得了射弹超空泡形态的高速影像;分析了射弹超空泡形态演变过程及超空泡形态参量(长度、直径、长细比等)随空化数变化的规律。研究结果表明:空泡直径变化分为快速扩张阶段、缓变阶段和快速收缩阶段,三个阶段耗时比约为7∶12∶6;空泡的扩张过程与收缩过程耗时比约为2∶1;空泡无量纲长度、无量纲直径和长细比均随着空化数的增大而减小。     

高超声速乘波飞行器三维流场的并行数值模拟

针对高超声速乘波飞行器三维绕流流场,在基于LINUX+MPI系统的分布式并行计算平台上,并行求解了三维雷诺平均的N-S方程。并行数值方法采用的是有限体积方法(FVM)、OC-TVD差分格式、B-L代数湍流模型及流场分区的并行方法。计算结果表明,所采用的并行数值模拟方法能够求解包含强激波的流场,激波穿越区域边界时无断层、错位等通量不守恒的现象。并行计算效率高,8个处理机计算时的并行加速比达到了6 8。     

流体力学中基本方程的特征方程及特征曲面

本文根据偏微分方程组的特征理论,讨论可压缩流体力学基本方程的特征方程和特征曲面,同时结合可压缩流体流动的实用例子,阐述确立其特征方程和特征曲面的具体方法。     

预处理方法求解低马赫数流场的边界条件

通过求解预处理后的二维可压缩Euler方程 ,数值试验了常见的三种进口边界条件的提法及其对稳定性和收敛性的影响。数值试验结果表明 ,预处理后 ,在低马赫数时基于黎曼不变量的边界条件为反射边界条件 ,而另两种边界条件为无反射边界条件     

三维超声速隔离段湍流内流场旋涡结构的数值模拟

从基于雷诺平均的N S方程出发 ,采用有限体积方法离散控制方程 ,数值模拟了三维超声速隔离段湍流内流场。计算中采用了二阶OC TVD差分格式、LU隐式算法和Baldwin Lomax代数湍流模型。数值结果与实验做了对比 ,并结合隔离段中的激波串结构分析了其横截面上旋涡结构的发展过程及不同外形条件下旋涡的不同结构。计算结果表明采用本文发展的方法模拟隔离段湍流流场是可行的 ,截面为正方形与长方形的隔离段内的涡旋结构截然不同     

预处理方法在全速度流场数值模拟中的应用

基于预处理方法,采用有限差分方法,LU-SGS隐式求解了二维可压Navier-Stokes方程,数值模拟了无粘突包、粘性空腔、化学非平衡流场,并与试验结果进行了对比。计算结果表明,所采用的数值模拟方法能够求解低速及亚、跨、超全速度流场和化学非平衡流场,所采用的预处理方法在亚声速范围内保持数值精度的前提下,加速了收敛。     

等截面隔离段中激波串结构的数值模拟

采用二阶TVD格式的有限体积法偶合Baldwin Lomax代数湍流模型求解雷诺平均Navier stokes方程 ,数值模拟了二维矩形等截面隔离段中的流动现象 ,很好地模拟出由激波 /附面层干扰所形成的复杂的激波串流场结构。计算结果与国外有关文献的结果进行了比较。     

含新型间断探测器的混合WCNS格式在间断无粘可压流的应用

为了让高精度数值格式在含间断和小尺度涡等复杂结构的超声速无粘可压缩流动情况下,仍能鲁棒地捕捉激波并快速得到流场高保真的模拟结果,研究了以子模板导数组合为基础的光滑度量算法,构造了精度与鲁棒性兼顾的新型间断探测器,使间断识别对小尺度涡也具有高分辨率;研究了混合加权紧致非线性格式(weighted compact nonlinear scheme, WCNS)方法,对流场中的光滑与间断区域分别使用线性与非线性加权格式求解,从而克服单一非线性格式在光滑区分辨率难以达到设计精度的问题。数值实验表明,使用新型间断探测器的混合WCNS格式对一维、二维Euler方程模拟结果良好,并且相比于在全流场使用局部特征分解的原始WCNS方法有计算效率的提高。     

交叉激波湍流边界层干扰问题数值模拟

基于BANS方程,采用有限体积法离散计算区域,对流项使用Harten TVD格式进行流场求解,湍流部分应用Menter SST两方程湍流模型,建立了隐式LU-SGS方法的三维CFD计算程序,数值模拟了来流马赫数4.96条件下的18°×18°对称双楔产生的交叉激波/湍流边界层干扰问题,计算了不同剖面沿展向压强分布以及流场的详细结构,并将计算结果与相关文献的实验与计算结果进行了比较.