含内热源的多孔方腔自然对流非正交MRT-LBM模拟

为了增强多孔方腔内流体流动与传热效果，采用非正交多松弛格子Boltzmann方法（MRT-LBM）对含有内热源的多孔方腔自然对流传热现象进行了数值模拟。研究了不同冷源布置方案（Scheme A~Scheme F）、内热源结构形式（Case 1、Case 2、Case 3）、内热源位置（a、b）、Darcy数、Rayleigh数等对多孔方腔内流体流动与传热的影响。计算结果表明:冷源布置方案对腔内流体流动与传热具有重要影响，当冷源左右对称布置时，腔内温度场及流场亦对称分布；在高Rayleigh数下采用Scheme A的双上部冷源布置方案能明显提高腔内的传热强度；内热源的形状对腔内对流传热影响很大，高Rayleigh数下，Case 3的布置方式更好。内热源的位置a和b对腔内的传热影响明显，提出了热壁面平均Nusselt数与位置a的拟合关系式，存在最佳的位置a（a=0.25），使得腔内的对流传热最强；热壁面平均Nusselt数亦随b值变化表现出特定的变化规律。随着b值的增加，热壁面平均Nusselt数呈现先增后减再增的趋势。      

二维有限长度柔性壁面上T-S波演化的数值研究

受自然界鸟类羽毛的柔性特征启发,利用数值模拟的方式开展了柔性壁面对亚声速边界层中T-S波演化的影响研究。刚性壁面上的数值结果与线性理论吻合得很好,验证了数值方法的可靠性。在此基础上,将部分刚性壁面替换为柔性壁面,结果表明,柔性壁面可以抑制T-S波在空间上的增长,从而推迟边界层流动转捩。壁面的变形不只跟随T-S波的波形,还因为柔性段与刚性段相接的前缘和后缘引起与扰动源频率相同的更大尺度的壁面波动,壁面的实际变形由这几种波叠加而成。开展的参数研究结果表明,增大表面的质量密度对于柔性壁面衰减扰动的效果几乎没有影响;增大表面张力和增加底部支撑的弹性系数可以增加壁面的刚性,减小壁面变形的幅度;增加阻尼可以抑制柔性段前后缘产生的大尺度壁面波动的传播,而对跟随T-S波的变形影响不大。总体上,柔性壁面的变形程度越大,其扰动的抑制效果越强。      

含内热源的多孔方腔流热耦合非正交 MRT-LBM数值模拟

针对含内热源的多孔方腔内自然对流现象问题,采用非正交多弛豫时间（MRT）格子Boltzmann方法进行了研究。分析了Rayleigh数（104 ≤ Ra ≤ 106）、内热源布局方式（水平、垂直及对角布局）、内热源几何尺寸大小（A=1/16,1/8,3/16,1/4）及两内热源间的间距（S=5/64,13/64,21/64）对流动传热的影响。结果表明:在Ra=104,105和S=5/64的情况下,任意内热源几何尺寸,内热源采用对角布局方式可获得更好的对流换热效果;在Ra=105,106和S=13/64,21/64的情况下,水平布局方式更优;在内热源采用水平布局,Ra=104的情况下,任意内热源几何尺寸,对流换热效果均随着内热源间距的增大而增强;而随着Ra增大,内热源几何尺寸减小,对流换热效果随着内热源间距的增大先增大后减小,而后随着内热源间距增大其对流换热效果减弱;对角布局也有相似规律,在其他条件一致的情况下,随着内热源几何尺寸的增加,其对流换热效果增强。      

二维微槽道电渗流的增强混合

微尺度流动中的混合一直是微流系统研制中的一个大问题,对壁面zeta势非均匀分布情况下的二维微槽道电渗流进行了数值分析,重点分析了这种情况下的微流动混合增强.基于电渗流的Helmholtz-Smoluchowski滑移速度模型,应用SIMPLE方法计算了4种滑移速度在整个壁面上的分布不均匀一致时的微槽道电渗流动的流场及被动标量在其中的混合.计算结果表明,滑移速度的非均匀一致导致微槽道电渗流中的流线发生剧烈的转折,从而使得微槽道流动的混合效率有很大的提高.流场中流线的转折越剧烈、转折的次数越多,微槽道流动混合效率的提高就越大.     

一种混合粒度奇偶校验故障注入检测方法

为了实现高效的抗故障注入攻击,提出了一种混合粒度奇偶校验故障注入检测方法。传统奇偶校验检测方法为每n比特设置一个奇偶位,表示该n比特的奇偶性。随着n的减小,奇偶位个数增加,资源消耗增加,检测率提高。为了实现故障检测率和资源消耗的折中,对电路故障注入敏感部分或关键部分处理的数据采用细粒度奇偶校验（即n值较小）,对其他部分采用粗粒度奇偶校验。以RC5加密算法为例,阐述了混合粒度奇偶校验故障检测方法的原理和应用,并对不同粒度奇偶校验方法的故障检测率及资源使用进行了理论分析。实验结果表明,与整个RC5电路都采用字（n=32 bit）奇偶校验相比,混合粒度奇偶校验故障注入检测方法可以提高故障检测率29.44%,仅增加资源消耗2.48%。      

着陆器足垫冲击月壤动态行为离散元仿真分析

研究着陆器与月壤的相互作用对着陆稳定性、安全性和着陆机构设计具有重要的意义。建立了具有角度的着陆器足垫冲击月壤的动态作用力学模型,并通过离散单元法分析了足垫在0°、8°、15°和22°冲击角度和1 m/s、2 m/s和4 m/s冲击速度下月壤的动态变化情况,得出了在冲击速度4 m/s、冲击角度为0°时,冲击力最大为5 080.6 N;在竖向位移相同的情况下,冲击角度越大冲击力越小。另外,足垫冲击月壤角度越大对月壤的压实作用越小,其中22°时,前后孔隙率变化为26.5%。冲击角度一定时,速度越大对月壤的压实较滞后,而压实程度较大,冲击速度为4 m/s时,前后孔隙率变化为78.1%。研究结果对着陆机构设计和着陆控制具有参考意义。     

不同热解温度下酚醛树脂复合材料渗透率测试

为了获得不同热解温度下酚醛树脂复合材料渗透率,设计搭建了复合材料气体渗透过程实验测量装置,提出了一种测量复合材料渗透率的方法,基于达西定律获得了复合材料渗透率。以不同热解温度下酚醛树脂复合材料为研究对象,进行了气体渗透过程实验测量,获得了材料上下表面气体压力变化曲线,同时得到了渗透过材料的气体流量,进而通过达西定律得出材料的渗透率。结果表明,该实验装置能够实现复杂孔隙复合材料的渗透率的测量。整体上,热解温度越高,渗透率越大。热解温度为400℃,渗透率量级在10-13;600℃和800℃时,渗透率量级在10-11。材料渗透率K和热解温度T满足K=9.7×10-14T-4×10-11。该实验结果为该类材料的渗透和热质扩散性能分析提供了基础物性数据。      

基于红外热像仪对珠状凝结换热过程实验

通过非接触式测温装置红外热像仪,观测冷凝液滴的形成、生长以及液滴之间汇聚、滚落过程.实验发现在珠状冷凝过程中,其表面温度变化显著.随着冷凝过程的进行,小液滴逐渐汇聚成大液滴,其液滴单位面积换热量明显下降,液滴表面温度逐渐降低.微纳米级小液滴在珠状冷凝换热中具有至关重要的作用.液滴在洁净表面上的珠状凝结和液滴在汇合后自清理表面上的珠状凝结表现出不同的温度分布,即珠状凝结汇聚后的自清理表面一直被前驱膜所覆盖.      

高压高剪切率下密封环缝隙流体动态润滑特性

为减少高压高剪切率下热效应对密封环缝隙流体润滑性能的影响，以密封环缝隙流体剪切速度梯度、壁面剪切应力和油膜温升表征摩擦副润滑性能，采用RNG k-ε湍流模型，通过Workbench建立流-固-热多场模型，计算得到不同主轴转角下缝隙流动形态变化规律和温升前后油膜厚度变化量、不同L型槽内外径比及长径比下流速流型分布规律、壁面剪切应力、温度分布和热变形量。研究结果表明:密封环L型槽黏性底层随剪切速度增大而增厚，油膜厚度随运动周期增加而变薄；当L型槽内外径比小于1.07时，随比值减小密封环壁面剪切应力不断增大，最大变化率为7%；密封环L型槽长径比在0.19时平均壁面剪切应力达到最大，当长径比继续增大时，油膜区剪切速度梯度逐渐减小，油膜温升和热变形亦随之减小。研究结果可为马达优化以减少能量损失和改善密封环润滑条件提供理论指导和依据。      

基于能量观点的混合层流优化设计

为了合理地在混合层流设计中减小阻力,降低能量消耗,利用吸气控制功率消耗与阻力、吸气速度的关系式,建立了考虑以吸气功率最小为优化目标的优化设计方法。该优化设计方法采用了自由变形（FFD）参数化方法,紧支型的径向基函数（RBF）动网格技术,改进的微分进化（DE）算法,以及耦合基于eN转捩预测的RANS流场高精度求解器。针对25°后掠角的跨声速无限展长后掠翼,进行了以阻力最小为优化目标的均匀吸气和以功率消耗最小为优化目标的分布式吸气的混合层流优化设计。优化结果表明,基于能量观点的优化结果在雷诺数10×106下可以达到均匀吸气的阻力收益,相比初始构型,阻力降低了29.1%,上下翼面转捩位置分别推迟了18%和15%弦长,功耗降低了1.7%;而在雷诺数20×106状态下,相比初始构型,阻力减小了41.3%,比均匀吸气阻力优化结果提高了4.5%,上下翼面转捩位置分别推迟了52%和14%弦长,功耗降低了8.14%。优化结果表明,建立的基于能量观点的混合层流优化方法是可行的。      

二维翼型绕流的CE/SE方法

采用时空守恒元和解元CE/SE(space-time Conservation Element and Solution Element method)法,求解二维Euler方程,开展了翼型绕流的无粘数值模拟研究.用非敏感克朗数计算格式消除克朗数过小引起的数值耗散对解的污染,结合当地时间步长法,解决网格不均匀引起的当地克朗数变化跨度大的问题.对NACA0012翼型的无激波流场进行了二维数值模拟,并与AGARD算例做了对比.结果表明:CE/SE方法的计算结果与AGARD结果吻合得很好,为该数值计算方法对翼型绕流数值模拟的进一步应用奠定了基础.      

多级降压调节阀内部流动特性数值模拟

为研究多级降压调节阀内部流动过程和各降压级的压降影响因素及规律,基于Mixture两相流模型和Schnerr-Sauer空化模型的数值计算方法,建立了多级降压调节阀内部流场的计算模型。通过数值模拟计算,获得了调节阀内部的流动特性,分析和探讨不同压差和不同阀门开度下各降压级压降特性。数值计算结果表明:流体流经阀门节流元件时,由于多股流体的强烈掺混和剪切使得阀门内存在大量的漩涡;各降压级的压降存在一定的差异,第3级套筒处的空化程度最大,阀座处更易发生空化;阀门在同一开度下,工质沿各降压级的压降在不同压差下的变化趋势相同;第3降压级的压降随阀门开度的减小而增大,在25%开度下阀门第3降压级的平均压降最大,达到了8 MPa,相同流量下该降压级小孔中工质的速度也增加。研究结果为降压阀的合理设计提供了理论依据和参考价值。      

单组元推力器毛细管两相流影响分析

单组元推力器在工作过程中产生的高温向毛细管处传递,使毛细管内产生气液两相流动。随着推进剂流量的减小,气相组分将不断增加,并会产生明显的气阻。文章就气阻对推力器的影响进行了专题研究,采用FLUENT软件,基于体积分数模型和多相流控制方程,对84种条件下毛细管内两相流的分布情况和两相流对推力器毛细管压差的影响情况进行了数值模拟,并通过两种推力器热试车开展了毛细管内两相流对推力器燃压影响的专项试验。数值模拟结果和试验结果表明,在毛细管内径足够小时,毛细管内的两相流引起的气阻均会导致毛细管内部压降的突然变化,从而引起推力器燃压阶跃,燃压逐渐降低,并进一步影响推力,使推力逐渐下降。     

蜻蜓爬升过程飞行特征实验研究

昆虫真实飞行过程中的飞行特征是仿生流体力学机理研究的基础和关键。本文针对蜻蜓(黄蜻)大、小爬升角2种飞行状态的运动规律和动力学特性展开研究,根据蜻蜓的趋光特性诱导其进行爬升飞行。采用2台光轴相互垂直的高速摄像机进行拍摄,通过特征点匹配和三维重构方法准确地捕捉了2种爬升飞行过程中蜻蜓身体和翅膀的运动参数,并进行动力学特性对比分析。实验结果表明:蜻蜓在进行大爬升角爬升时需要的上升力大于向前的推力,最大拍动幅度比小爬升角爬升时约大40%;扑翼频率比小爬升角爬升时约大3.3 Hz,前后翅相位差相对于小爬升角爬升时减小20°以上;另外,大爬升角爬升过程中前翅前倾角度更大,这样能够使蜻蜓身体保持更大的俯仰角,翅膀能够获得更大的上升力。      

返回式卫星烧蚀热防护机理与数值模拟

详细讨论了低温炭化材料的质量损失和吸热机理。利用分层模型,给出具有不同物理和化学特性的各层的能量守恒关系。用积分法给出各层热防护的数值模拟,数值模拟计算结果与一维热传导解析解结果一致,亦与地面烧蚀实验結果一致。     

Lunar gravity field modeling critical analysis and challenges

This paper summarizes and provides a critical analysis of the historical developments of lunar gravitational models from the earliest use of ground based tracking systems of the Lunar Orbiter to the Lunar Prospector mission. This encompasses a comprehensive and critical analysis of the various methods used in the estimation of the gravity coefficients and the processing of large batches of diverse measurements and data types. It has been shown that weakness exists in the current models of the lunar gravity field, which is primarily due to the lack of far side lunar tracking data information, which makes the lunar potential modeling difficult but expected to be overcome as data from SELENE satellite-to-satellite tracking becomes available. Comparisons of various lunar models reveal an agreement in the low order coefficients of the spherical harmonics. However, substantial differences in the models exist in the higher-order harmonics. A numerical comparison has been presented showing the performance of all the contemporary lunar gravitational models used within the astrodynamics community and available in public domain. Improvements to the current models are part of a continuing process and the recent model improvements and future possibilities in lunar gravity modeling are discussed. A brief review of the recent missions has been presented. It is hoped that this critical review will benefit the researchers by presenting the historical as well as state of the art in this field.     

侧向微重力环境下板式表面张力贮箱内推进剂晃动行为分析

侧向微重力是航天器在轨飞行时在东西位置保持和南北位置保持状态时所处的加速度环境.由于航天器贮箱内推进剂在侧向加速度环境下的重定位过程易产生晃动，因此侧向加速度环境对贮箱内推进剂管理装置(PMD)的管理能力的要求更高.为确保板式贮箱对推进剂的在轨管理能力，需要开展一系列的数值仿真与试验对PMD的管理能力进行验证.本文以板式表面张力贮箱为研究对象，采用VOF两相流模型对侧向微重力环境下重定位过程中贮箱内推进剂的运动进行数值模拟，并对推进剂运动中的晃动行为进行分析.最终，通过数值模拟结果验证了板式表面贮箱的管理能力，为微重力落塔试验与空间站试验提供参考依据.     

模拟弹塑性碰撞的修正SPH方法

从基本的无网格光滑粒子法SPH(Smoothed Particle Hydrodynamics)近似出发,修正了模拟固体力学中大变形弹塑性碰撞的SPH方法.在边界处采用修正的边界条件,弹塑性分析过程中采用增量理论计算应力,迭代过程中用守恒光滑法进行滤波修正消除拉力不稳定.对SPH方法进行了程序实现,给出了杆弹塑性碰撞的算例.计算分析表明,SPH方法节点的影响域较大、精度较相同节点间距有限元法的结果有一定差距,但是通过增加粒子数量可以提高SPH的精度,保持了其简单性和计算大变形的特性.      

卫星规避方案量化分析方法研究

空间碎片数量日益增多,为保护在轨卫星的运行安全,需要针对碎片碰撞制定规避方案,用来规避碰撞风险.本文在调研国外卫星规避机动流程的基础上,分别研究了卫星规避相关的轨道机动方式、轨道预报模型以及交会风险计算模型,建立了卫星规避方案量化分析方法,为规避方案的选择提供参考.