加权变系数的瞬态导热有限元法

研究了两点时间差分格式的瞬态导热温度场有限元法数值解的稳定性和振荡现象,推导出消除振荡的方法:将时间差分格式的加权系数σ与时间步长△t分别和温度振荡联系起来,得出变σ和变△t的函数关系式。      

应用分支和突变理论对飞机空间机动稳定性的研究

 本文应用分支和突变理论,研究了飞机空间机动全局稳定性,导出了求解滚转速率临界稳定值的解析公式,分析了非线性跳跃现象和周期性振荡现象,并找出了这两种现象产生的原因。     

出口封闭的冲压发动机进气道激波振荡现象

采用计算流体力学和风洞试验两种手段,对冲压发动机出口封闭的多种超声速进气道内通道气流振荡现象进行了研究。结果表明,在一定的超声速范围内,封闭的进气道内通道发生了较强气流振荡现象,振荡频率、幅值和内通道长度、飞行马赫数有关,该振荡现象是一种自激振荡现象。利用亥姆霍兹共振器频率公式可以对封闭的进气道内通道振荡频率进行预估,并对其归纳出一种修正方法,可以结合模型试验结果对进气道和燃烧室内流振荡频率和幅值进行分析,为飞行试验提供预估和参考。     

断续磨削时的脉动温度场解析

文章构造了周期变化的移动热源模型,并引进卷积概念建立了计算断续磨削时工件表层非稳态脉动温度场的理论公式。该公式可以包容连续磨削的情况,且可计算任意时刻的瞬态温度分布。文中利用推得的公式对断续磨削时的温度场及其降温机理作了深入的研究     

球面螺旋槽动静压气体轴承动态特性试验

根据轴径偏离稳态平衡位置进行变位运动时气膜力的变化规律,建立气膜刚度和阻尼系数的计算矩阵方程,使用滚动迭代算法计算气膜瞬态刚度和阻尼系数。分析升、降速阶段气体轴承在周期1、周期2和混沌运动过程中瞬态刚度和阻尼的振动特征,研究气膜涡动和振荡的力学机理,以改善气体轴承动态特性、减小气膜涡动和振荡,提高气体轴承运行的稳定性。结果表明:通过气膜刚度阻尼系数变化特征可以判断转子的运行状态。当气膜力波动,刚度阻尼系数出现周期性波动变化,引起气膜振荡;当气膜瞬时刚度系数波动幅值增大10倍以上,轴承转子接近临界失稳进入混沌运行状态;当气膜刚度阻尼系数发生畸变,轴承出现碰磨失效。     

飞机刹车副温度场的瞬态有限元模型

飞机制动过程中刹车装置产生的高温可使摩擦材料发生热降解,导致摩擦系数变化,出现热衰退现象,影响飞机的安全性能,因此飞机刹车副瞬态温度场分析日益受到关注。针对飞机刹车副瞬态温度场建立了具有对流和热辐射两种边界条件的有限元计算模型,分析了边界条件对温度场的影响,并用ANSYS对其进行仿真计算。仿真计算得到的温度与刹车副经验温度基本吻合,证明了所建立的有限元计算模型合理、可行,为飞机刹车装置温度场分析、计算和研究提供理论支持。     

发动机对进气温度畸／瞬变响应的理论研究

建立了模拟发动机工作稳定性的“激盘-滞后-容积”模型，并借助于平行压气机理论，利用数值积分法计算了发动机对进气温度畸/瞬变的响应。重点研究了稳态畸变对喘振工作边界及温度突升对发动机工作稳定性的影响，得到了瞬态稳定工作边界；还研究了稳态畸变与温度突升组合情况对发动机工作稳定性的影响。计算结果与国外资料比较，曲线的趋势一致，结论基本相符。     

带侵蚀效应的三维双燃速装药固体火箭发动机点火过程

针对三维双燃速装药固体火箭发动机,采用商用CFD计算软件的用户自定义方程对点火瞬态过程中的侵蚀燃速、燃面温度、质量源项、动量源项、能量源项及点火器入口等进行了二次开发,并进行了考虑侵蚀效应的点火瞬态过程流场仿真计算,与地面点火试验内弹道数据进行了对比分析。结果表明:点火初期发动机头部的气流流动十分复杂,且头部压强随时间变化十分剧烈;堵盖打开会造成头部升压梯度振荡,并产生负升压梯度;仿真与试验结果符合较好,说明采用一维非定常固相传热方程计算三维燃面温度可以满足当前工程要求,计算方法正确有效。      

喷管喉衬结构瞬态温度测量及导热反问题方法应用

为了获得实时可靠的固体火箭发动机喉衬结构瞬态温度,采用内埋快响应烧蚀热电偶对发动机喉衬结构瞬态温度进行测量,并应用导热反问题方法(Beck序列函数法)计算其喉部内壁面温度和热流密度。采用的试验方案有效避免了传感器探头与喉衬内壁面烧蚀速率不同而干扰流场或中途烧毁等现象发生。试验结果显示越靠近内壁面的测量点温度峰值越高,测点温度瞬态变化与发动机工作状态同步性越好。计算结果显示,计算的平均误差为5.6%,最大误差为7.9%,其内壁温度最高出现在发动机工作结束时,热流密度在0.1s内迅速上升,而后又迅速下降,最大值约为30MW/m2。在发动机停止工作2s后,喉衬开始向外界传热。     

不规则地形条件下双向DMFT电波传播特性算法研究

 在传统中心差分离散混合傅里叶变换(DMFT)抛物方程模型基础上,针对传统算法计算动态阻抗边界条件时的数值振荡及忽略后向传播与散射造成计算误差增大的问题,提出基于前-后向混合差分方法的递归双向DMFT(TW-DMFT)模型,以一阶前-后向差分方程代替传统的中心差分方程拟合阻抗边界条件,以递归双向传播模型改进传统单向预测模型.仿真分析了双向DMFT模型计算不同边界及媒介条件下单刃峰和实际地形的电波传播分布特性,与传统算法及几何绕射理论(GTD)进行比较,新算法提高稳定性及精度的同时弥补了传统算法的不足,更具有普适性.     

ENO-AUSMPW格式对激波涡干扰的数值模拟

高阶ENO格式在复杂流场计算中具有重要的地位,它理论上可以达到任何阶的精度。本文采用二维高阶ENO插值方法,构造了高阶AUSMPW格式,并对包含复杂激波结构的激波涡干扰流场进行了数值模拟。计算结果表明了高阶ENO插值下的AUSMPW格式具有较高的激波分辨率和较低的数值耗散。     

结冰风洞扩散段损失分析研究

针对在建的高亚声速结冰风洞,用数值方法研究了风洞的部段损失。利用CFD技术,数值模拟了不同设计状态下风洞收缩段、试验段、扩散段的损失,并与试验结果进行对比,分析了扩散角对扩压损失的影响。对扩散段进行了优化设计改造,成功达到了风洞性能指标。研究结果表明：利用数值模拟方法,可以方便地进行各种设计方案的损失比较,既经济便捷又可靠,在风洞设计中可大力推广。     

更准确、更精确、更高效——高超声速流动数值模拟研究进展

从准度、精度和效率3方面回顾了近几十年来高超声速流动数值模拟研究的进展。在物理模型方面,介绍了高超声速数值模拟中高温气体效应、稀薄气体效应以及湍流效应的建模与模拟,基于雷诺平均Navier-Stokes(RANS)方程重点对现阶段较为关注的高超声速边界层转捩的模式理论研究进行了介绍。在空间离散算法方面,主要介绍了高超声速数值模拟中常用的二阶精度迎风格式以及高阶精度格式的发展及其应用。在时间推进方面,主要回顾了隐式时间推进方法的发展及其应用。在误差和不确定度估计方面,主要介绍了其概念、来源以及常用的分析方法,同时给出了迭代误差估计、Richardson外插法以及敏感性导数方法等初步研究结果。最后,讨论了高超声速流动数值模拟中下一步需关注的问题。     

基于BVD原理的高保真空间重构方法

简要综述了一类基于单元边界变差最小化(Boundary Variation Diminishing,BVD)原理,设计双曲守恒律高保真数值格式的空间重构方法。BVD原理要求尽量减少通过重构得到的网格边界两侧物理量之间的差,从而能够有效地控制黎曼求解器中的数值黏性。BVD方法针对数值解的空间分布特征,选择多个函数作为空间重构的候补函数,并根据BVD判定准则从候补函数中选取最合适的函数进行空间重构。BVD判据不需要根据求解对象进行经验参数(阈值)的调整。选用适当的候补函数和BVD准则,可以完全避免现有算法中为抑制数值振荡而必须采用的非线性限制。BVD格式能在抑制数值振荡的同时,有效地控制数值耗散,可以对光滑解与间断解都获得高保真的计算结果。本文概述了BVD方法的基本思想、设计相关格式的基本思路,以及一些具有很强实用价值的BVD格式。并通过单相和两相可压缩流动的一些典型算例验证BVD格式的特点和优势。     

Optimal guidance against active defense ballistic missiles via differential game strategies

The optimal guidance problem for an interceptor against a ballistic missile with active defense is investigated in this paper. A class of optimal guidance schemes are proposed based on linear quadratic differential game method and numerical solution of Riccati differential equation. By choosing proper parameters, the proposed guidance schemes are able to drive the interceptor to the target and away from the defender simultaneously. Additionally, fuel cost, control saturation, chattering phenomenon and parameters selection were taken into account. Satisfaction of the proposed guidance schemes of the saddle point condition is proven theoretically. Finally, nonlinear numerical examples are included to demonstrate the effectiveness and performance of the developed guidance approaches. Comparison of control performance between different guidance schemes are presented and analysis.     

非定常数值模拟方法的发展及其在动态绕流中的应用

基于混合通量分裂的思想,通过应用Gauss-Seidel迭代求解差分方程,构造了一种时空二阶精度、无条件稳定的隐式迭代NND算法,并讨论了时间精度与稳定性,亚迭代收敛判则,几何守恒律的应用以及动壁边界条件等相关问题.通过引入一种简便易行的加权函数来综合刚性动网格生成技术和超限插值生成动网格这两种方法的优点,发展了一种计算量小、比较实用的加权动网格生成技术.作为应用实例,本文给出了多个动态物体绕流的数值模拟算例,计算结果表明了本文数值方法的成功.     

旋转坐标系下分区计算的LU隐式方法

将隐式时间迭代方法应用于并行计算是CFD研究中的热点问题,研究发现LU-SGS格式在旋转坐标系下对分区边界比较敏感,对旋翼悬停状态进行分区并行计算时碰到了计算发散的问题。针对此问题,对基于LU分解的3种隐式时间迭代格式（LU-SGS、DP-LUR和HLU-SGS）进行了对比研究,设计了静止流场下的旋转网格算例对计算方法进行测试。结果表明,LU-SGS格式在网格边界采用简化处理方法,当边界的逆变速度增大时会引起数值误差放大,误差的不断积累导致计算发散。DP-LUR和HLU-SGS格式通过在边界单元采用雅克比迭代算法,能有效消除分区边界影响,使计算格式在大CFL数条件下保持稳定,其中HLU-SGS继承了LU-SGS迭代效率高的特点。在此基础上,采用并行化的LU隐式方法对Caradonna-Tung旋翼进行了并行计算,针对两套稀密程度不同的背景网格,在产生复杂分区边界的计算网格条件下,均获得了气动力和旋翼尾迹都充分收敛的流场。通过计算,分析了背景网格密度对预测气动力性能和捕捉尾迹流场的影响,计算结果与试验值吻合良好,验证了当前隐式计算方法适用于旋转坐标系并行计算,适合于推广至其他大规模并行分区的流场计算。     

SA和SST湍流模型对高超声速边界层强制转捩的适应性

凸起物是高超声速流动中常用的一种人工转捩装置。采用高阶精度算法模拟了高超声速进气道压缩面上的强制转捩流动,转捩装置为一排高度为1mm的钻石型凸起物和斜坡型凸起物。考察了Spalart-Allmaras(SA)模型和剪切应力输运(SST)湍流模型对该问题的适应性。在考察过程中通过丰富的算例分析了网格规模、可压缩修正和空间离散格式等对计算结果的影响。在层流区,计算能与试验取得非常一致的结果。但在湍流区,计算得到的热流通常高出试验数据。经分析发现其原因是强制转捩的湍流边界层与自然转捩的湍流边界层在涡结构上存在较大差别,使得湍流模型的效果较差。针对高超强制转捩湍流涡结构丰富的特点,对SST湍流模型进行了修改。计算结果表明,该修改方法对提高热流精度具有一定效果。     

LU-AUSMPW混合格式在可压缩无粘和粘性流动数值分析中的应用

本文首先将 AUSMPW格式与三阶 MUSCL格式融合, 给出了其在任意曲线坐标下的三维形式, 并将其与 LU-SGS格式结合, 应用于可压缩 Euler和 Navier-Stokes方程的求解。其次, 构造了一种新的通量限制器。最后, 为了验证 LU -AUSMPW混合格式的性能, 将平面叶栅跨音速无粘流动以及喷管超音速粘性流动作为算例。本文计算结果与文献计算结果和实验数据相符很好, 表明采用 LU-AUSMPW混合格式数值模拟可压缩流场, 具有较高的计算精度、较快的收敛速度和良好的稳定性。      

非结构网格上欧拉方程的自适应逆风隐式有限元解法

本文发展了一种解二维欧拉方程的隐式逆风有限元格式。空间近似应用的是二阶精度逆风方法；时间近似上根据数值通量函数的线性化处理导出了有限元的隐格式。为减少数值振荡，引入了基于特征变量的限制因子。在经自适应处理的非结构网络上，本文对跨、超、高超声速绕流进行了数值模拟。计算实践表明发展的隐式有限元格式与显格式相比具有较高的计算效率和较强的适应性。