液氢贮箱停放过程中的力热分析

为了寻找到运载火箭长时间停放过程中液氢贮箱的最经济液位,用于制定合理的发射流程以及紧急处置方法,采用计算流体力学(CFD)技术,对某型运载火箭停放期间液氢贮箱的力、热情况进行了仿真计算和分析。计算选用了VOF(Volume-of-fluid)两相流模型以及Lee相变模型,为了提高Lee模型在不同压力情况下对相变过程的模拟精度,采用安托因方程修正了该模型。修正后的模型首先由试验数据校验了其精确性,随后开展的液氢贮箱停放过程仿真结果表明：贮箱的竖直方向与径向均存在温度分层的现象,液相内会形成大的漩涡,该漩涡会使得冷热流体不断进行热交换,并导致贮箱内部的液氢出现气化。贮箱停放期间蒸发率最大值超过2 m 3/h,发生在停放4 h左右;而贮箱液位充填至37 m 3以上或17 m 3以下时蒸发率较低,最小值接近1 m 3/h。     

多臂空间机器人的视觉伺服与协调控制

针对多臂空间机器人自主目标抓捕任务,首先建立多臂空间机器人的运动模型和其与目标的相对运动模型,采用Kane方法建立多臂空间机器人的动力学模型;其次,研究基于视觉伺服的机械臂在线轨迹规划算法,并引入零反作用机动,消除机械臂运动对平台姿态的扰动;再次,在不使用零反作用机动功能时,分别使用基于角动量前馈补偿的协调控制算法和逆动力学方法设计了协调控制器,在机械臂运动时保持平台姿态和相对目标的位置。最后,开发了基于Matlab的仿真软件MASS（多臂空间机器人仿真）,仿真结果校验了上述方法的有效性。     

固体火箭发动机燃烧室三维流动数值计算

采用ｋε双方程湍流模型,以ＳＩＭＰＬＥ 计算程式求解翼柱型装药固体火箭发动机燃烧室内三维非定常不可压流ＮＳ 方程．在建立翼柱型装药简化模型的基础上采用边界标志法来表达燃面推移,采用从二维到三维的初场给定方法结合多重网格法求出了多个时间步的非稳态流场结构．计算结果表明燃烧室内旋涡运动呈现一定的空间与时间分布,周向加质导致燃气通道横断面上轴向速度分布的极不均衡．     

基于释热率分析的钝体贫油熄火过程分析

针对锥形钝体稳定的甲烷-空气预混湍流火焰复杂的熄火过程,采用大涡模拟(LES)与输运方程概率密度函数(TPDF)湍流燃烧模型相结合的模拟方法,研究远离熄火、近熄火及熄火点的火焰和释热率(HRR)数值变化情况,定量分析熄火判据。结果表明:冷态速度模拟结果和实验的相对均方根误差在10%以内,热态误差在20%以内;释热率是否出现在OH和CH₂O重叠的区域,是判断熄火的一个重要参数;远离熄火时,释热率高的区域主要沿内侧剪切层出现;近熄火工况下,释热率在流向轴上闭合,回流区下游也出现较大的释热率;熄火点工况下,较大释热率的区域在回流区下游和上游均有出现;模拟预测的熄火情况和实验PLIF结果一致;平均释热率可作为判断熄火的定量依据,即当钝体后方0.2d处内侧剪切层平均释热率与回流区平均释热率的比值小于4时,发生熄火。      

多芯片组件（MCM）用金锡焊带材料电阻率及热导率温度特性研究

对牌号为Au₈₀Sn₂₀的金锡焊带材料在208~423 K的电阻率及热导率与温度的函数关系进行了研究,并对其在多芯片组件（MCM）中的传热效果进行了评估。分别对材料在208~423 K中5个温度点的电阻率及4个温度点的热导率进行了测试,基于理论模型建立电阻率/热导率随温度变化的函数关系,最终采用模拟热扩散数值方法评估材料在高温下的传热能力。结果表明,采用修正函数模型后,金锡焊带材料在208~423 K下热导率与电阻率的关系符合测试结果,随芯片表面温度的边界条件从208 K升高至423 K,采用变温热导率模型得到的热流密度模拟计算结果相比理想化恒定热导率模型的差异性逐渐升高至5.5%。综上,金锡焊带材料热导率与电阻率的关系符合Wiedemann-Franz法则修正后的Smith-Palmer方程,在该材料传热设计时应考虑其热导率温度效应。     

超声速混合层气动声场的数值模拟

超声速混合层的噪声是气动声学研究的经典问题,制约着航空发动机性能的提升。本文用数值方法对超声速混合层声辐射特性进行研究,具体方法:用抛物化稳定性方程(PSE)计算混合层近场大尺度扰动演化,结合一种基于渐进展开和稳定性分析的积分预测远场声压,得到快慢两种模态马赫波辐射的强度和方向。分析了上下侧来流温度对气动声场的影响,结果表明,上下侧温度比小于1时,慢模态占优;而上下侧温度比大于1时,快模态占优。     

广义Burgers方程声爆传播模型高阶格式离散

声爆预测技术是研制新一代环保型超声速民机的关键技术之一。针对当前基于广义Burgers方程的声爆远场传播模型的数值求解格式精度较低的情况,开展了高阶格式离散方法研究。通过分析该模型中各效应项的物理性质,应用合适的高阶精度格式对各项分别进行离散求解,并采用经典几何声学射线法计算声爆传播路径,实现了对地面声爆波形的精准预测。通过对美国F-5E机型的声爆飞行试验和第二届国际声爆预测研讨会的典型算例进行数值模拟,验证了本文高阶离散方法的可靠性。进一步的结果分析表明,采用高阶离散格式的预测方法相比传统二阶精度预测方法具有明显的低耗散特性,在同等网格数目上可获得更高分辨率的计算结果,且其网格收敛性远优于传统二阶精度格式,同时在计算效率上也有一定优势;在Burgers方程的各效应项中,非线性项的影响更为明显,因而采用高阶离散格式时其优势更为突出;热黏性吸收项对数值结果影响较小,采用高阶离散格式对计算精度的提升效果也不明显,实际计算中仍可以采用传统离散格式,甚至忽略该项的贡献。     

马赫数对侧压式高超音速进气道及等直隔离段三维内流场的影响的数值分析

用 N-S方程和 RNG k-ε紊流模型计算了 RBCC用侧压式高超音速进气道三维内流场 ,重点分析了马赫数对流场的影响。从出口截面上气流参数的均匀程度来看 ,Ma3 时比 Ma6时的更均匀 ;从压缩气流的主要部位来看 ,Ma6情况下在支板前肩点之后到喉部截面之前的这个等宽度收敛通道内 ,而 Ma3 时主要是在前肩点之前的收敛通道内 ;从对气流的压缩程度来看 ,Ma6比 Ma3 更大 ;从唇口激波的强度来看 ,Ma6比 Ma3 更强 ;从进气道的流量捕获率来看 ,Ma6比 Ma3 更高 ;高超音速进气道内靠近侧壁以及侧壁和顶板的角区附近是最有可能出现亚音速的区域。     

用非结构直角网格和欧拉方程计算飞机绕流

 采用八叉树结构 ,生成复杂外形绕流计算的非结构直角网格。物面附近用投影方法 ,使网格贴体。并将Jameson的有限体积法推广用于这种网格的欧拉方程计算。对歼击机模型的绕流计算表明 ,网格生成的机时花费很少 ,总体质量好 ,因而欧拉方程解算的收敛质量也好。     

基于欧拉方程的跨音速外形优化设计方法

将流场分析和优化方法耦合起来，形成一种基于飞机轴向截面积分布的设计方法，用以减小跨音速和超音速飞机的波阻，在本文方法中，流场用欧拉方程进行解算，以阻力参数构成目标函数，以多个截面面积（或当量旋成体半径）为约束，选择合适的优化变量，采用传统的梯度法进行优化设计，文中利用非精确搜索措施，避免了用传统一维搜索所需的大量机时，本算法所需机时少，优化效果明显，可用于飞机的初步设计。