超声速进气道可压及不可压流动数值模拟

采用Chio-Merkle预处理矩阵对可压NS(Navier-Stokes)方程组进行时间预处理,分析了预处理方法的物理和数学背景.用有限体积方法,结合LU-SGS(Lower-Upper-Symmetric-Gauss-Seidel)隐式时间积分和AUSM+(P)(Precondition Advection Upstream Splitting Method)格式、中心差分2种空间离散格式,求解该预处理NS方程组.通过圆弧凸包流动和二维方腔顶盖驱动流动的数值试验表明,该方法克服了传统时间迭代方法模拟低速流动时的刚性问题,加速了收敛过程,可以同时有效地模拟可压及不可压流动.用此方法与块结构化网格技术相结合,进行了在不同飞行马赫数、攻角和出口反压条件下,混压式轴对称超声速进气道三维流场数值模拟.计算结果表明:该方法能准确地捕获复杂激波系,清楚地揭示进气道三维流动现象;进气道性能随工作条件的变化规律,与理论分析一致.      

复杂外形下的无网格算法研究

研究了二维、三维复杂外形下的无网格算法。在一种布置点云方法的基础上,发展一种曲面拟合的重构方式构造流场物理量,并应用于AUSM+_up格式计算欧拉方程的数值通量;计算采用了隐式时间推进,并引入当地时间步长和残值光顺等加速收敛措施。通过对某三段翼型低速流动、M6机翼跨音速流动、某全机跨音速流动进行了数值模拟,表明本文发展的无网格算法能有效地模拟复杂外形的无粘绕流。     

反射器夹层面板精密成形原理

对高精度、大面积不可展双曲反射面的制造问题,提出了基于"离散钉模、真空负压、蜂窝夹层、应力释放、回弹补偿"的高精度反射面板成形原理;建立了四边自由面内多点约束的薄板大挠度弯曲解法及夹层板成形过程仿真模型,用于确定成形工艺参数;利用动力显式和静力隐式算法相结合的方法实现成形过程仿真;通过开缝和开槽方法释放面内压应力,确定了面板开缝和开槽方案;通过逆向修正离散模包络面,使其与成形模拟得到的夹层板型面误差最小,得到模具修正因子,补偿回弹误差.在此基础上开发了夹层板精密成形CAE系统.该方法制造的面板型面精度RMS达到25 μm以下,已成功地应用于多项大型紧缩场和毫米波天线的制造.      

基于活动标架对Gosper曲线的研究

基于活动标架定义离散曲率^[1-3],并借此研究了Peano-Gosper分形曲线,得出一种新的生成Peano-Gosper空间填充曲线算法。通过建立曲线弯曲点位置号码与弯曲点曲率的映射,对于任意正整数n,能够输出对应弯曲点的曲率并画出相应弯曲点图象结构。相比于传统生成Peano-Gosper分形曲线的方法,本文算法减少迭代次数提高效率。     

燃烧过程辐射换热离散坐标模型的应用

对离散坐标(DO)辐射传热模型分别在长方体空间和旋流燃烧室内进行了验证与应用。长方体空间的计算结果表明:DO模型的纯辐射传热预报结果比六热流(SHF)模型准确,与计算精确的区域法的计算结果吻合较好,而SHF模型的预报结果却有较大的误差。对于旋流燃烧室内温度场的模拟预报也表明,DO模型的计算结果和实验值相吻合,明显高于用SHF模型计算所得的温度值。另外,DO模型在计算时间和边界处理等方面比区域法优越;在收敛速度和计算准确度等方面也比工程常用的六通量模型优越。      

高速旋转冲压增程弹用进气道复杂流场数值模拟

采用块结构网格与二阶精度流场分区求解技术,对高速旋转、含侧向支柱冲压增程弹丸进气道内外复杂流场进行了数值模拟,得到了高速旋转工况下对应于不同来流马赫数和攻角,以及临界工况时超音速进气道内外流场复杂的波系结构。在旋转工况下,进气道流场结构与不旋转时相似,但旋转速度以及攻角的存在对冲压发动机进气道的总体性能产生了负面影响,进气道总压恢复系数和动能系数均有所降低,而流场畸变指数则显著增大;冲压发动机进气道在较低马赫数工作时的综合性能优于在较高马赫数工作时的综合性能。     

应用CCSDS空间数据链路安全协议的卫星遥控安全方案探讨

针对国内卫星遥控体制国际化的需求,结合已有的空间数据系统咨询委员会(CCSDS)遥控协议体系架构,设计了应用CCSDS空间数据链路安全协议(SDLS)的遥控安全方案。该方案选用认证和加密认证2种业务,并针对不同的遥控安全业务,分别设计了定制遥控帧格式、安全关联、算法管理及密钥管理业务。以目标机为核心搭建验证系统,通过2个典型应用场景验证遥控安全方案的可行性,结果表明:该方案可提高遥控安全系统的灵活性、可扩展性,并且能够有效保护密钥强度,延长密钥生命周期。遥控安全方案的设计及其在遥控单元中的实现,可为我国卫星遥控体制与国际接轨提供参考。     

二维复值Ginzburg-Landau方程的一个高阶紧致ADI差分格式

对二维复值金兹堡朗道(Ginzburg-Landau,GL)方程提出一个基于时间分裂的高阶紧致交替方向隐式有限差分格式。本文通过时间分裂法将GL方程分裂成一个非线性子问题及两个线性子问题,对非线性子问题以及其中一个线性子问题均通过精确积分进行计算,并对另一线性子问题构造紧致交替方向隐式差分格式进行数值计算。实际计算中,在每一时间步,利用追赶法求解一族常系数三对角线性代数方程组,从而使得算法既具有较高精度又拥有较快的计算速度。数值实验表明该算法在时间和空间方向分别具有二阶和四阶精度,并模拟了方程的一些动力学行为。     

基于PIV速度场测量的压强梯度计算

介绍了基于PIV速度场测量结果计算压强梯度的基本原理和3种不同的计算方法,并通过对流高斯涡模型系统研究了离散格式选取、PIV参数设置和流场特征等对压强梯度计算的影响.结果表明:PIV速度测量的精度是压强梯度计算的关键性影响因素,1%噪声对结果的影响比其他过程高约2个数量级;时间分辨率和空间分辨率的设置并非越高越好,而是有个合适的中间值;拉格朗日方法在大多数情况下优于欧拉方法,但当流动结构较复杂时性能较差,应根据流动类型和流场特征合理选取压强梯度计算方法.      

辐射换热对平板气膜冷却性能影响

为了分析辐射换热对平板气膜冷却性能的影响,分别对不考虑辐射换热(耦合计算)与考虑辐射换热两种条件下的平板隔热屏进行了三维流热耦合数值模拟研究,并与绝热壁(不考虑隔热屏固壁导热)的计算结果进行了对比。辐射换热采用离散坐标法。得到了3种不同主/次流总温比(Rt)条件下气膜冷却效果、流场、单位面积冷流体热负荷以及流量系数的对比结果与变化规律。结果表明:气膜冷却效果沿流向逐渐降低,且随总温比的增加而降低;在展向上气膜冷却效果呈现中间高两边低的分布规律,最大与最小相对值在绝热壁条件下达到最大,为94.5%,辐射换热次之,为13.5%,不考虑辐射换热时最小,为9.8%。与不考虑辐射换热结果相比,辐射换热对气膜冷却效果的影响在流动方向逐渐变大,气膜孔附近展向温度梯度增大,冷流体热负荷增加了92.8%,二者计算的流量系数相差不大,为1.1%,与绝热壁计算结果相比,流量系数减小了13.1%。随着总温比增加流量系数减小,最大值与最小值相差不大,为1.1%。