板式表面张力贮箱导流板的理论计算与试验验证

板式表面张力贮箱的PMD通常包括导流板与出口附近的蓄液器,导流板是影响板式贮箱性能的关键因素之一,由于地面没有失重环境无法直接进行导流板的工作原理地面验证。为了确认导流板的可用性,通过理论分析建立了导流板液体传输能力的计算方法,通过中性浮力试验对导流板的液体收集能力进行了试验验证。经过多次试验的试验结果与理论计算结果的对比分析,得到了一种可应用于工程设计的计算方法,利用该方法编程形成了一套计算软件,可指导板式表面张力贮箱内导流板的设计。     

蓄液器在板式贮箱中的应用及性能分析

蓄液器是板式表面张力贮箱最重要的组件之一,它即可以和导流板搭配使用,也可以单独使用,来对推进剂进行有效管理。介绍了几种典型蓄液器的结构,在此基础上,介绍了蓄液器的工作原理,以及蓄液器在微重力下工作的基本方程,同时给出了蓄液器的临界加速度、最大蓄液量、临界蓄液量、推进剂残留量、有效蓄液量的求解方法,并采用这些方法对一典型的蓄液器的性能进行了分析。通过这些性能分析,可以为蓄液器的优化设计和微重力试验提供指导和依据。     

微重力下板式贮箱内液体晃动性能研究

板式贮箱内推进剂管理装置(PMD)以板式结构部件为主,其抑制液体晃动性能应能够满足卫星平台高稳定度和快速机动的需求。针对微重力下板式贮箱内液体晃动性能,建立基于有限体积(VOF)方法的气液两相流流动模型,通过数值模拟计算手段,比较分析研究了微重力环境下有防晃叶片和无防晃叶片的板式PMD晃动性能。计算分析结果表明,有防晃叶片的板式PMD抑制液体晃动性能明显优于无防晃叶片的板式PMD,使液面快速趋于稳定。同时,将计算分析结果与液体晃动微重力试验结果进行了分析比较,微重力试验结果与计算分析结果非常相近。研究结果表明,有防晃叶片的板式PMD具有更好的抑制液体晃动功能,能够显著地抑制贮箱内推进剂晃动,以满足卫星平台的高稳定度和快速机动的需求。     

微重力下球形贮箱内液体晃动有限元模拟

以卫星球形贮箱为研究对象,用任意拉格朗日-欧拉(ALE)法建立了数学模型,给出了计算的基本方程,并基于MSC-Dytran有限元程序,对常重力和微重力环境中贮箱内推进剂的晃动进行了数值模拟。结果表明,微重力环境中液体更易发生大幅晃动,甚至产生飞溅,工程中常用的等效力学模型已不再适用,而有限元数值模拟方法可较好地解决该问题,为卫星及大充液复杂结构航天器控制系统的设计和分析提供理论参考。     

带多充液圆柱贮箱航天器刚-液耦合动力学研究

文中以在低重环境下带多充液圆柱贮箱刚性航天器中刚-液耦合方程的建立和求解为主要研究目的。推导航天器中充液圆柱贮箱内任意点的牵连运动方程,根据壁面边界条件给出了贮箱内液体牵连晃动势的表达式;利用第二类边界条件下的傅立叶-贝塞尔级数展开法对低重力环境下的弯曲自由液面处的复杂动力学边界条件进行处理,建立以液体相对晃动势的模态坐标和晃动波高的模态坐标为状态向量的液体耦合晃动力学方程,通过积分分别得到了耦合晃动力和耦合晃动力矩的解析式;运用准坐标系下的拉格朗日方程建立以航天器主刚体姿态坐标和轨道坐标为状态向量的刚体耦合运动动力学方程,进一步联立上述耦合方程得到航天器整体系统的刚-液耦合动力学状态方程;最后,编制出适用于带多充液圆柱贮箱航天器内刚-液耦合动力学计算的模块化计算程序,通过计算实例验证所编程序的准确性的同时,研究了携带多充液箱航天器系统贮箱布局、外激励方式对航天器刚-液耦合系统动力学特性的影响。     

表面张力贮箱的微重力实验验证--静平衡与重定位

新研制的表面张力贮箱保证在微重力或无重力环境下、在任何时候都能向系统提供不含气的推进剂，因而需要确定液体的空间分布以及在给定微重力加速度下的重定位过程，以便为贮箱的设计提供依据。为此，按照相似理论用缩比模型进行了一系列微重力落塔实验，然后根据测定的重定位时间计算出原型表面张力贮箱的重定位时间。     

微重力环境下推进剂贮箱中三维气液平衡界面的数值模拟

以部分管理表面张力贮箱的管理舱为研究对象,利用三维气液平衡界面计算程序Surface Evolver,在无重力和微重力且几何边界条件比较复杂的环境下对管理舱内的气液平衡界面进行数值模拟;计算结果与已经应用卫星的理论计算完全吻合。     

萤火一号火星探测器防晃型贮箱技术

针对萤火一号(YH-1)火星探测器贮箱需具有防晃和蓄液排气等工作特点,对防晃型贮箱设计进行了研究.根据表面张力和微重力流体力学原理,设计了实现该功能的系列推进剂管理装置(PMD).对层流式和叶片式两种PMD进行了微重力流体仿真和试验.经计算机仿真和落塔试验验证,优选出最佳PMD方案,设计的防晃贮箱满足探测器要求.     

液氧贮箱增压过程研究

采用数值模拟方法对液氧贮箱增压过程进行研究.贮箱内流场采用流体体积函数（VOF）多相流模型考虑,选择标准双方程k-ε湍流模型分析湍流效应,气液两项之间的热量、质量转移通过自定义程序（UDF）求解.获得了贮箱压力、排液流量、气垫温度、液氧温度对贮箱内流场温度分布的影响.计算结果表明,在稳定增压过程中,贮箱液面无扰动,贮箱内温度分层分布;各参数变化时,对贮箱内温度分布的影响主要是温度梯度的变化,并且各工况下液面附近和扩散器附近温度梯度基本相同.     

微重力环境中重力梯度加速度引起的流体晃动运动

本文研究了关于旋转轴在贮箱的非对称轴上且远离贮箱的几何中心情况下，流体在微重力环境中由重力梯度加速度诱发的晃动特性，建立了问题的数学模型并对模型进行了数值模拟。以高级Ｘ射线天文物理实验卫星（简称ＡＸＡＦ─Ｓ）作为研究对象，获得了由旋转运动引起的重力梯度加速度的数学表达式。晃动问题的数值计算以与卫星固连的非惯性坐标系为基础，目的是寻求一种较为易处理且适合于流体力学方程的边界和初始条件。通过数值计算获得了流体作用于卫星贮箱上的力和力矩。     

第二代表面张力贮箱的研究与应用进展

为了对国外新型卫星推进剂贮箱的研究动态进行跟踪，对第二代表面张力贮箱的研究与应用进展做了阐述。首先对比了不同类型推进剂贮箱的特点，在此基础上，从试验研究、数值模拟和实际应用三个角度对第：二代表面张力贮箱——以板式结构为主的表面张力贮箱的发展做了综述。研究表明，国际上已经对第二代表面张力贮箱的性能做了大量的数值模拟和搭载试验，如今第二代表面张力贮箱已经成为当今国际大型卫星推进剂贮箱的主流。这为新型卫星推进剂贮箱的研究指明方向并提供有价值的参考。     

斜切反喷管性能分析

固体火箭发动机前端斜切反喷管，其结构简单、作用时间短、气动型面具有尖点，并在超音速区有台阶，喷管内存在一系列激波，并伴有流动分离现象。本文从雷诺平均的非定常Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程出发，结合采用Ｂｏｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型，利用时间相关法及ＭａｃＯｃｏｒｍａｒｋ两步显格式求解，模拟了斜切反喷管流场。计算得到的壁面压强分布与风洞吹风实验测得的压强分布相一致。该方法可应用于斜切反喷管的     

轨控喷流干扰气动-运动一体化非定常数值模拟

针对运动条件对反作用控制喷流控制气动干扰的影响,基于雷诺平均N-S方程,建立非定常气动--运动一体化方法,通过计算流体力学(CFD)/刚体动力学(RBD)耦合求解,对一种采用轨控喷流直接力控制的锥柱裙外形进行了数值模拟,研究了干扰流场结构、力和力矩与飞行器状态随时间变化特性,分析了运动非定常计算和定常计算的对比结果以及外流参数(马赫数,攻角等)变化对气动运动特性的影响。研究结果表明：轨控喷流使飞行器运动时存在强烈的非定常特性,定常计算的法向力放大系数和干扰俯仰力矩大于运动非定常计算结果。     

动网格在固体火箭发动机非稳态工作过程中的应用

利用Fluent流场计算软件、动网格技术、UDF文件,用DEFINE_GRID_MOTION定义燃面边界的移动,用DEFINE_PROFILE定义边界类型,考虑侵蚀燃烧、压强变化率对推进剂燃速的影响,对轴对称变截面固体火箭发动机的非稳态工作过程内流场进行了瞬态分析.得到了变截面轴对称固体发动机稳态工作过程中装药燃面推移图像,并得到了发动机内弹道参数分布云图及其随时间的变化规律.     

基于预处理的非稳态流场计算方法

采用预处理方法可提高Navier-Stokes方程或Eu ler方程计算低马赫数流场的准确度和收敛性能。采用双时间步和预处理相结合的方法计算非稳态流场,保证了采用预处理方法后计算结果的时间精确性。在每一物理时间步,控制方程虚拟时间导数项采用LU隐式时间格式,对流通量采用AUSM 和MUSCL格式进行求解。通过对前台阶流和圆柱绕流的计算,验证了该方法的广泛适用性和计算结果的可靠性。     

推力矢量发动机燃气舵舵间干扰的数值分析

通过求解N-S方程,对某空空导弹燃气舵在不同状况下的三维湍流干扰流场进行数值仿真,得到了燃气舵在不同舵偏角下绕流流场的复杂波系结构,计算了测量舵在无干扰舵和有干扰舵情况下升力和阻力等参量随舵偏角的变化曲线。所得舵间干扰相对误差不超过2%,与实际要求相吻合。结果表明,所采用的流动模型和数值模拟方法对燃气舵的气动性能进行预示是有效的。     

马赫数振荡状态下带抽吸槽进气道非定常数值模拟

研究了在来流马赫数振荡状态下带抽吸槽的二维混压式超音速进气道的气动特性,通过给定非定常边界条件,对飞行马赫数为2．2、振荡马赫数幅值为0．154的进气道非定常流场进行了数值模拟。与定常条件下数值模拟结果进行对比,结果显示非定常流场与定常流场有较大差异。在振荡状态下,进气道的性能发生周期性变化,存在一个椭圆形或类似椭圆形的迟滞回路。     

应用动网格技术计算着陆过程中缓冲发动机喷流流场

利用Fluent软件对着陆过程中缓冲发动机的喷流流场进行数值计算,发动机与地面相对位置的不断变化使计算域形状不断改变.为了保证计算的正确性,必须采用动网格技术.介绍了实现移动网格的基本方法,并结合使用网格弹性平滑移动和网格重构方法实现了着陆过程中缓冲发动机喷流流场的非定常数值计算;得到了随发动机着陆产生的喷流流场的压强...     

燃气源控制气流非稳态耦合传热问题研究(Ⅰ)--计算模型和方法

针对某燃气源控制气流与导热套环之间的非稳态耦合传热问题,建立了数学模型,并采用整体求解的方法进行求解,即把流体区域的流动和对流换热与固体区域的导热过程组合起来,作为一个统一换热过程来处理。为确保计算结果在物理上的真实性,采用“假密度”法求解以温度T为求解变量的能量方程。对于非稳态控制方程离散后的离散方程采用了PISO方法进行求解。