关于几个二阶显式差分格式的评述

本文着重讨论了MacCormack格式计算头部混合流场时出现的一些问题,当用不同的边界点算法与之配合时对计算稳定性的影响以及格式的稳定性在计算过程中的表现形式,比较了MacCormack格式与其他二阶格式的计算结果,在相同的边界算法的情况下各种格式的效能。本文也提出了一种新的物面边界的算法,其中用了特征概念,但又避免了传统特征线法的迭代过程。     

固体推进剂最大固体填充量的计算

本文提出了复合固体推进剂中最大固体颗粒填充量的概念,并通过多组分固体颗粒混合物料平均空隙率的计算来预估最大固体填充体积分数。通过有限级配的球形过氯酸铵进行了验证,表明了计算结果和实测他一致。     

一种简化算法的航空发动机全状态数学模型

通过一种简化的燃气空间模型对某型双转子混排加力涡扇发动机各部件进行了建模和仿真计算,并给出了详细的计算流程以及数学公式.在仿真过程中结合实际工程应用经验,对涡轮处采用一个局部的小循环迭代以减小涡轮部件内燃气空间计算的时间步长,提高计算精度.与其他数值方法相比,该简化算法避免了求解非线性方程组时大范围的循环迭代,提高了全状态数学模型的计算速度.通过仿真结果与发动机实际试车数据的对比,证明该方法有足够的计算精度.      

纪念诸德超主编

20 0 3年 1月 8日 ,《航空学报》编辑委员会主任委员、主编 ,中国著名计算力学和结构动力学专家 ,航空教育家诸德超教授 ,因病抢救无效 ,在北京逝世 ,享年 69岁。在他生命的最后几天 ,诸德超教授作为本刊主编 ,仍在病床上坚持批阅稿件。今天他的夫人林翠虹老师返回了他在 1月 7日上午一边输液 ,一边批阅的二十多篇《航空学报》稿件。看着这些批阅意见 ,我们的心情非常沉重 ,仿佛看到一位教育家 ,一缕烛火 ,一边在与病魔顽强搏斗 ,一边仍关心着他的事业 ,照亮着后人的前路。诸主编走了 ,但是他的智慧 ,严谨和奉献 ,仍将鼓舞着我们这些还在航…     

多级轴流压气机展向掺混的数值模拟

在紊流扩散模型的基础上, 通过引入二次流造成的气流过转折/欠转折影响, 完成了S₂流面的展向掺混数值模拟。分别用Galimore-Cumpsty掺混模型或Baldwin-Lomax紊流模型计及紊流扩散的作用;落后角的预估在二维Carter公式基础上引入了Roberts二次流修正模型。应用上述方法对两台3级轴流压气机进行了计算。计算的出口总温、总压等参数与实验较吻合, 与其他单独考虑紊流扩散和二次流效应的掺混计算结果比较表明了本文方法的有效性。      

对二维复杂湍流剪切流动数学模型的新探索

在复杂湍流剪切流中经常会遇到雷诺剪应力与法向速度梯度的反号问题,此时涡粘性假设不再适用。Hinze从湍流结构中存在大涡的概念出发提出了包含台劳级数二阶项的剪应力表达式,取得了正确的定性结果。作者进一步分析了Hinze公式,认为他未能取得合理的定量结果的主要原因,在于他的公式隐含着假设流体微团在大涡的输送作用下能对原有的平均动量有足够的记忆。本文从实验与计算的比较以及近代对湍流结构的知识出发,提出了新的假设。认为较小尺度的湍流结构在大涡的输送作用下能对所包含的剪应力有足够的记忆,在此假设下,推出了新的湍流剪应力表达式,计算结果与实验数据的初步比较说明了此假设的合理性。     

高超声速化学非平衡流动Navier-Stokes/DSMC耦合算法

基于相同的化学反应模型,在已有计算流体力学（CFD）和直接仿真蒙特卡罗（DSMC）方法及程序的基础上,采用Modular Particle-Continuum（MPC）耦合技术,建立了包含化学非平衡Navier-Stokes/DSMC耦合算法。算法结构中DSMC计算区域在CFD计算结果上根据当地克努森数自动选取。发展了适用于流场分区信息交换的亚松弛技术,抑制DSMC方法对CFD计算的影响。把DSMC方法和CFD的应用范围拓展到过渡流区,为复杂飞行器近连续过渡流区高超声速化学非平衡流动数值模拟研究提供了一种工程适用的预测分析手段。通过对二维圆柱高超声速化学非平衡绕流的算例与其他结果的比较研究,表明耦合算法不论在流场结构、流场非平衡现象,还是飞行器表面参数、整体气动力/热特性方面,都能够得到与全DSMC计算吻合的结果,证实了所建立的Navier-Stokes/DSMC耦合计算模型与方法的有效性和可靠性。仿真了某航天器解体碎片在过渡区的化学非平衡流动,得到碎片在过渡区的气动力/热特性,为碎片的陨落计算提供依据。     

向量计算及其在刀具上的应用

董仁杨同志是一九六○年大学毕业,分配担任工具设计员。他以毛主席的光辉哲学思想为指针,用矛盾对立统一观点,分析了切削运动的各个因素,抓住了刀刃和切削速度这一对基本矛盾,创新地提出了刀具几何角度的向量定义,以及进行计算的一套比较系统的新概念、新方法,用以解决刀具设计、制造、检验、安装中各种平面和空间的计算问题,尤其是解决运动状态下复杂的刀具几何角度的计算,能使复杂问题简化,解     

适合机动飞行的旋翼气动模型研究

为改善对旋翼空气动力他轴响应的预测,在旋翼气动模型中引入时变的尾迹弯曲参数K_R,以体现机动旋翼的尾迹畸变效应.与已有方法相比,该模型的计算结果更为接近实验数据.还建立了动态尾迹畸变模型,并分析了计入尾迹收缩效应的必要性.结果表明,为了正确计算机动飞行中的旋翼气动力,特别是操纵的他轴响应,应当在模型中采用随状态变化的尾迹弯曲参数K_R.      

车载大视场视觉惯性组合导航算法研究

车载视觉导航在大角度旋转和高速运动等场景下表现不佳,由于视觉和惯性的误差特性互补,所以视觉惯性组合导航算法有良好的应用前景。大视场相机获取外界信息丰富,但大视场相机采样频率低影响视觉数据采集,而且基于优化的组合导航算法计算速度慢,故提出了一种适用于大视场相机的车载视觉惯性组合导航算法。该算法侧重于数据预处理：设计了基于伽马函数的对比度受限自适应直方图均衡化来降低图像噪声,通过FAST角点与光流法的组合完成特征提取与跟踪,进一步加快了计算速度。引入了惯性数据等效旋转矢量二子样算法来减少数据量和计算量,通过预积分算法为误差优化减少计算量,通过视觉误差、惯性误差以及边缘化误差联合优化来提高精度。与其他组合导航算法相比,该算法精度较高,计算速度最快,与VINS-Mono相比节省时间19.7%。     

高超声速流场与结构温度场一体化计算方法

在对国内外流场与结构温度场一体化计算方法的不足进行细致讨论的基础上,提出了一种高超声速流场与结构温度场一体化计算方法.采用统一的积分方程组作为气动加热和结构传热物理过程的控制方程,对整个物理场进行统一的迎风格式有限体积方法离散,给出了流场与结构交界面上温度、温度梯度及导热系数等参数的计算方法.在时间推进方面,定常状态采用多步龙格库塔迭代格式,非定常状态则采用双时间步长方法.采用发展的一体化算法对二维圆管模型的气动加热和结构传热问题进行了数值模拟.结果表明:2s时驻点物面温度为390K,与其他文献的误差在3.1K范围内;稳定时驻点物面温度为647K.      

利用有限域上二元多项式构造带仲裁的认证码

在传统的无条件安全认证码中,假定发方和收方是相互信任的,他们不会互相欺骗。但有些情况也可能不是这样,发方在发出一个消息后,不承认是他发的,收方捏造了一个消息,声称来自发方,在这种情况下就需要增添仲裁方。带仲裁的认证码能解决通信系统中发方与收方互不信任的问题。利用有限域上二元多项式构造了一个带仲裁的认证码,并计算了所构造码的参数。当编码规则按等概率分布选取时,计算出了敌方与收方成功模仿攻击和成功替换攻击的最大概率．以及发方成功模仿攻击的最大概率。     

气体运动论BGK格式的翼型绕流数值模拟

以气体运动论BGK(Bhatnagaar-Gross-Krook)格式为基础,通过改进非均匀网格下的通量计算方法,把气体运动论BGK方法应用到二维贴体网格上模拟了可压缩的无黏和黏性流动。在重构阶段,以交界面周围4个网格单元为一组,先把每个网格单元的宏观量(密度、动量和总能量)转换在相对于交界面的坐标系下,然后使用van Leer限制器插值求出交界面的宏观量,在求解阶段再使用BGK理论求出流过交界面的通量。通过上述方法使贴体网格和均匀网格,在计算格式上一致,与其他求解二维贴体网格的气体运动论BGK方法相比,该方法相对简单,计算量减小。结合贴体网格通量计算和气体运动论BGK格式方法理论,分别对亚、跨、超声速下的NACA0012翼型绕流进行了数值模拟和对比研究,所得可压缩无黏流动的数值结果与北大西洋公约组织航天研究与发展咨询组(AGARD)的结果吻合良好。本文还模拟了低雷诺数状态下的可压缩黏性流动,所得结果与文献中的数值计算结果吻合良好,证明了该方法在数值模拟及激波捕捉方面具有鲁棒性,为在复杂流动中的应用提供了一种简单的途径。     

扩散抛物化(DP)NS方程组的意义及其在计算流体力学中的应用

本文首先讨论扩散抛物化（DP）NS方程组的早期研究工作：它的提出、数学性质、意义和在CFD中的应用，然后讨论扩散抛物化理论的一些新发展。这些新发展是对NS方程组数值计算进行物理分析的基础上得到的，其中包括NS方程组差分计算时，粘性剪切流对网格间距和格式精度的要求；粘性项只保留剪切项的广义扩散抛物化（GDP）NS方程组，它的性质和应用。由于高Re数流动在NS方程组的差分计算中，网格Re数彼此相差悬殊的特点，产生了计算离散单元守恒方程组的新的算法思路，即离散流体力学（DFD）算法。在DFD算法中需要同时计算三种不同的守恒方程组（Euler,DPNS和NS方程组）。本文讨论了DFD算法中需要同时计算三种不同的守恒方程组（Euler,DPNS和NS方程组）。本文讨论了DFD格式的构造、它的优点和应用。并以超声速绕前后台阶流动为算例，来说明GDPNS方程组的用处和DFD算法的优点。DPNS方程组、GDPNS方程组、DFD算法是高智提出的，对这些问题，他的合作者从理论、算法、数值验证和某些应用又作了系统的研究。     

基于注意力模型的视觉/惯性组合里程计算法研究

由于外界环境的干扰和传感器精度的限制,视觉/惯性组合里程计的输入数据存在一定的噪声,这会增加里程计的解算误差,而且误差会随着时间积累。针对以上问题,设计了一种基于注意力模型的视觉/惯性组合里程计算法。该算法使用卷积神经网络和长短时记忆网络分别构建了视觉特征提取器与惯导信息特征提取器,同时引入了两种注意力模型:加权组合网络以及开关组合网络,对视觉特征信息和惯导特征信息的融合噪声进行降噪处理。通过在组合里程计算法中添加闭环校正环节,有效地抑制了里程计误差随时间的积累。对比实验结果表明,设计的组合里程计算法与其他算法相比,无论在性能上还是在精度上都有明显的提升。     

机动飞行中旋翼尾迹畸变效应实验

为准确计算机动飞行时的旋翼他轴响应,气动模型应计入尾迹畸变效应。体现此效应的关键在于推出时变的尾迹弯曲参数KR。本文通过烟流实验,得出随旋翼前进比、拉力系数和机动角速率变化的KR表达式,并分析指出,国外文献采用固定KR值不足以体现机动效应。为了给理论预测提供参考,本文还测量了机动旋翼的空气动力时间历程,表明旋翼轴接受准阶跃输入后,同轴及他轴的空气动力响应有过冲,且俯仰速率愈大,过冲幅值愈大;本实验中他轴响应为主响应的15%。     

飞行仿真中的机动生成器

讨论了作为数字飞行仿真计算中的关键环节的轨迹跟踪和指不乔量两种机动生成器的原理与应用，并根据常赠模型他机动生成器设计参数的选取问题，机动生成器由于采用统一的比例－积分－微分（ＰＩＤ）反馈形式，简化并改进了以往数字飞行仿真中有关环节的处理，在使用中具有良好的效果。     

一种新的三维四步编织复合材料几何模型及其在宏观弹性性能预测中的应用

 为了准确预测三维四步(1×1)编织复合材料的宏观弹性常数,提出了一种考虑纤维束三维几何结构和纤维束与基体相互作用、便于有限元网格划分的单胞几何模型,该模型精确地反映了纤维体积比。以该几何模型为基础,用六面体单元分别对纤维束和基体进行了网格划分,满足了纤维束与基体界面间的位移连续条件。采用有限元模型计算了编织复合材料的宏观弹性常数,与其他计算结果和实验结果的比较表明,该几何模型和有限元模型比较合理。     

基于等效靶方法的飞机单击中易损性通用计算模型

飞机单击中易损性/生存力经常表示为飞机在给定一次威胁(例如碎片)随机打击条件下的杀伤概率及易损面积结果.将弹药威力评估领域中的"等效靶"方法引入到飞机易损性评估中,提出了一种飞机在弹丸或导弹碎片单次打击下的杀伤概率及易损面积的通用计算模型.为了更加反映真实情况,该模型考虑了一枚碎片(弹丸)连续穿透部件后运动状态变化的影响.应用表明,所提出的通用模型解决了易损性评估中Pk/h(部件给定打击下的杀伤概率)的计算问题,与其他通常采用的单击中模型相比更容易实现计算机化.     

具有离体涡的边条布局的低速气动力数值计算

使用等效集中涡和输入涡来模拟边缘离体涡的涡格法,在本文中应用到边条机翼布局的气动力计算中。在进行边条机翼的“混合流型”的气动力计算时,对边条采用了产生离体涡的计算模式,在主翼上则采用了附着流动的计算模式。使用这种方法,对两组边条布局进行了计算,计算结果与实验数据进行了比较。