重叠网格方法研究多子弹抛撒气动问题

子母弹抛撒气动问题是一类典型的有相对运动的多体气动问题,也是重叠网格方法最有价值的应用领域之一。采用MI-Grid重叠网格软件系统,耦合流动力学方程和动力学/运动学方程,研究子弹分离过程的气动性能和运动轨迹,并结合正交试验设计,讨论和评估了出舱条件及气动干扰对子弹分离运动的影响,指出子弹在抛撒流场中受激波干扰对姿态运动影响很大,单子弹抛撒时母弹激波扫射使子弹俯仰特性剧烈波动,多子弹抛撒时弹舱流动受子弹挤压反射激波使子弹姿态抬头。     

带制退器的膛口流场数值模拟

带制退器的膛口流场是一种具有高度瞬变性、强烈方向性的复杂波系结构流场.采用分块网格划分的整体运动处理方法和基于动态层变方法的结构动网格技术处理弹丸高速运动造成的网格变化,模拟某型带炮口制退器的火炮发射过程中的非定常燃气射流流场,根据数值结果,详细讨论了初始流场、火药燃气流场的形成与发展以及与弹丸的耦合和相互作用过程.数值计算结果可以为炮口制退器的设计与实验提供一定理论指导,为深入研究奠定良好基础.     

基于两种模型的燃气轮机级内流场数值模拟的分析比较

借助Fluent软件,采用混合平面模型和滑移网格模型分别对某船用增压锅炉燃气轮机级内流场进行了数值模拟(CFD)分析,并对比了两种模型条件下燃气轮机级出口参数的变化规律.对比分析结果表明:非定常干扰因素主要来自动静叶片排间的相对运动和每个叶片排流场周向非均匀性引起的非定常流动,动静叶栅之间的非定常干扰因素是导致级出口参数不同于定常情况下的根本原因,从而进一步揭示了燃气轮机级内流动具有固有非定常性的物理本质.     

运载火箭芯级与助推级分离的数值模拟

采用动态重叠网格方法和双欧法,耦合求解非定常流动控制NS方程和6DOF运动方程,数值模拟运载火箭助推级和芯级的自由分离过程和附加外力作用下的强迫分离过程,并对芯级、助推级分离的主要影响因素进行分析。数值结果与试验结果吻合,表明动态重叠网格方法能正确预测火箭芯级与助推级间复杂的流动现象。根据数值模拟结果,给出助推级在不同分离工况下的分离参数及分离规律。     

PVM平台下三维超声速钝头体底部流场的并行数值模拟及其拓扑分析

用拼接法划分网格，并运用基于消息传递的并行软件平台ＰＶＭ分配计算任务，实现了较大规模的三维超声速钝头体底部流场的并行数值模拟。结果表明，底部压强特性与实验符合较好，流场的拓扑结构合理。     

CFD在高速滑行表面性能预报中的应用

针对滑行表面数值预报结果中船底出现异常水气分布的问题,基于三种常用CFD软件FLUENT、CFX及STAR-CCM+研究了VOF界面插值格式、网格类型及数目等因素对滑行表面性能数值预报的影响,并提出相应改善数值计算模拟的具体措施。研究表明:基于相同的结构化网格及计算设置时,采用几何重构界面插值格式预报滑行块的阻力及底部水气分布上最为准确,CICSAM格式次之,修正HRIC格式最弱。通过选择合理的网格类型及网格数目也可以改善滑行船体性能预报精度。CFX采用四面体网格,STAR-CCM+采用多面体网格可以较为准确预报滑行艇体的流场和阻力。建议在滑行艇的阻力及流场预报中采用几何重构的VOF界面插值方式,选择合理的网格类型并进行适当的网格加密可以提高水气界面捕捉的精度,减小阻力预报的误差。     

基于迎风型有限体积法的非结构化网格流场格心数据高质量体绘制方法

在计算流体力学领域,对复杂流场进行数值模拟时常用非结构化网格划分空间,并广泛采用迎风型有限体积法求解流场.该方法将数值解存储在非结构化网格单元(控制体)中心,产生格心数据,并假设物理量在相邻单元邻接面上形成"间断".然而,在可视化此类数据时,目前的重构方法破坏了间断约束,导致"流间断"信息丢失,严重制约了绘制质量.为此...     

潜艇上浮出水运动及粘性流场直接数值模拟方法

为得到潜艇上浮出水运动及流场特性,以带流水孔柱形潜体为研究对象,基于RANS方程及SSTk-ω湍流模型,并应用重叠网格技术实现了不同工况下潜体上浮的直接数值模拟,获取了模型上浮过程中横倾、纵倾及上浮速度等参数的时历曲线及流场细节信息。研究结果表明:该方法较合理地反映了潜体的上浮运动规律,能够直观展现潜体出水及倾覆过程,可为潜艇上浮出水稳性的预报及分析提供可行途径。     

三维高超音速复杂流场的分区数值计算

基于拼接网格系统 ,推导了一种适用于高超音速复杂流场分区数值计算的区域边界格式 ,通过区域边界格式对区域边界上的通量进行守恒性处理 ;在此基础上 ,发展了一个高超音速复杂流场分区算法。对高超音速粘性绕流流场进行了分区数值模拟 ,取得了较好的结果。     

船用液力偶合器内部流场的CFD数值模拟

针对目前国内对液力偶合器流场三维仿真计算均采用的只取泵轮和涡轮等过度简化几何模型导致不能真实模拟内部流场,以及使部件强度分析无法开展的问题,依据"泵轮+涡轮+转动壳体"的原几何标准,采用CFD技术,并通过几何建模、网格划分和数值计算等步骤对某液力偶合器在多种滑差工况的内部流场进行了数值计算。数值计算结果表明:计算值与台架试验值的误差均满足工程要求,由此证明了该方法的正确性。     

混合尺寸目标电磁计算的时域伪谱/有限体积混合方法

针对单一算法对混合尺寸目标进行时域电磁分析的困难,提出一种时域伪谱(PSTD)同时域有限体积(FVTD)混合方法。FVTD可方便地分析复杂的几何结构和材料,但是难以计算电大尺寸的目标,PSTD则特别适合计算电大尺寸的规则结构,但在模拟复杂的几何结构尤其是带有曲边结构以及电大、电小共存结构时存在困难。混合方法克服了单独算法的缺点,融合各自的优势,提高了算法的求解能力和应用范围。为了减小两种算法连接边界带来的反射,采用了FVTD计算面均值的二次函数重构方法,给出了交叠网格和非交叠网格两种混合方案。数值试验表明,混合方法有较高的精度,具有时域分析混合尺寸目标电磁问题的能力。     

燃烧室-喷管流场数值解

本文用时间相关法计算了固体火箭发动机燃烧室─喷管亚跨声速流场数值解,控制方程用ＭａｃＣｏｒｍａｃｋ二步显格式：边界点参数用物理边界条件和参考平面上的特征方程计算。计算表明,达到收敛的数值积分步数比纯喷管的跨声速计算要多得多。虽然喷管壁上和轴线上的马赫数分布与纯喷管计算类似,但喷管中的等马赫线分布与纯喷管计算的结果［３］相差较远。     

椭圆型方程曲线网格的快速生成方法

本文在用椭圆型方程的网格生成方法生成曲线网格时,通过比较四种常用的计算格式,得到了收敛速度较快的生成方法。在计算流体力学领域,对于边界形状随时间变化的流动问题的数值计算,及自适应网格的生成等同题,生成网格的收敛速度直接影响流场的计算速度。在上述情况下使用本文的结果,能有效的减少机时,提高计算效率。     

水下应用栅格翼动态展开参数预示方法

基于水下应用栅格翼动态展开过程动力学模型,根据展开特征角度定常水动力方法,拟合获得展开全程受到的流体力矩,并引入考虑相对速度影响的修正因子,形成水下应用栅格翼动态展开过程参数工程预示方法,以典型展开时序点航行体运动参数为设计输入,对栅格翼展开过程运动参数进行预示。通过与水下应用栅格翼非定常流场仿真计算数据以及水下航行体弹射试验数据对比,验证了上述预示方法的正确性及工程适用性,为水下应用栅格翼方案设计的优化及展开不同步性分析提供设计参考。     

用LIIF法测量DF/HF化学激光器喷管混合性能的数值分析

激光诱导碘荧光(LIIF)是一种很好的流场测量技术,为了验证该方法测量化学激光器流场混合性能的精度,模拟得到连续波DF HF化学激光器的HYLTE喷管冷流场,并计算了用氦碘混合气体作为氧化剂流、用氦气作为燃料流的HYLTE喷管的流场。通过比较前一流场中的氟原子和后一流场中的碘分子的摩尔百分比分布情况得知,当作为氧化剂流的氦碘混合物中碘的质量百分比在54%左右时,用LIIF法测量DF HF化学激光器喷管的混合情况具有较高的精度。     

航天器编队飞行构形保持与重构的继电型控制

研究小偏心参考轨道航天器编队飞行构形保持与重构控制的工程实现技术。首先给出相对运动状态转移矩阵,并推导出常推力控制情况下相对运动状态递推的解析表达式;进而给出脉冲推力、连续变推力和继电型推力三种发动机的推力模型和示意图。分别将等速度冲量的三种推力模型代入相对运动状态转移矩阵中,比较相同条件下相对运动控制作用效果的差异,理论推导结果表明:在一阶意义上,三种推力模型对相对运动控制作用等效,因而航天器编队构形保持与重构控制可以基于继电型推力模型来实现。     

高超音速流场的隐式TVD格式并行计算

本文采用区域分割技术和拼接网格的并行策略，发展了一个适合于分布式存贮多机系统的ＴＶＤ隐式有限体积并行算法；并在ＰＶＭ并行环境下，对三维高超音速绕流流场实现了多机并行计算，通过负载平衡等方法得到了较高的加速比（在二处理机系统上加速比为１∶８４，在四处理机系统上为３∶４４）。     

平面运动激波绕障碍物流动的非结构Euler解

用三角形非结构网格及有限体积法求解Euler方程 ,对运动激波跨越不同截面形状的二维障碍物绕流进行了数值模拟。计算结果给出了不同时刻各流场的复杂结构 ,该结构与实验结果一致。计算表明用非结构算法模拟复杂外形绕流具有简单、准确的特点 ,并能避免尖角处非物理结果的产生     

视频传感器覆盖范围精细化计算方法

视频传感器覆盖范围计算是视频传感器网络初始配置、优化部署等的核心基础,其精度严重影响覆盖分析的结果,然而现有很多方法并未考虑计算方法执行效率与精度的影响,而方法的精度与效率在大范围监控视频覆盖分析应用中显得尤为重要。因此,提出一种兼顾速度与精度的视频传感器覆盖范围计算方法,将地理空间进行格网化;计算每个格网角点被视频传感器覆盖的状况,用0或1表示;接着,以格网四个角点的覆盖情况描述格网的覆盖情况,将格网分为16种覆盖状态;逐步细化整个格网覆盖状态非(0000)和(1111)的格网,直至达到预设精度;获得视频传感器的覆盖范围。实验结果表明,该方法能够更为精确地获得部署在地理场景中视频传感器的覆盖范围,同时兼顾了计算方法的速度和精度。