基于SA和SSTγ-■模型对边界层转捩预测的比较

利用作者共同研发的In-house代码TRANS3D平台,在SA和SST两种常用湍流模型框架下构建了γ-■转捩模型,并以二维低速S809层流翼型和三维小后掠跨音速F5层流机翼为对象,比较了两种不同湍流模型构架下γ-■转捩模型预测的气动力特性和流场分布。结果表明:两种转捩模型均能预测航空转捩计算中常见的分离流转捩和自然转捩类型,明显改善了中低雷诺数流场下的预测精度,但由于选取基准湍流模型的不同,基于SA和SST的γ-■转捩模型在流动细节上依然存在着一些差异。     

非均匀来流对新型高超弯曲激波二维进气道的影响

针对一种采用新型压升规律的曲面压缩面所设计的高超弯曲激波二维进气道,利用数值模拟方法研究了6种不同长度平板发展而来的非均匀流对其性能参数的影响,并和两种常规二维进气道吞非均匀来流的能力进行了比较,着重对比了设计Ma数下三种不同压缩形式进气道喉道截面的流场畸变程度。数值模拟结果表明,采用新型压升规律的高超弯曲激波二维进气道性能受来流气流不均匀性的影响较小,且喉道截面的流场畸变指数小,对非均匀来流具有一定的校正作用。     

基于人工粘性的间断Galerkin有限元方法

发展了一种基于人工粘性的间断Galerkin有限元方法,作出改进以增强适应性,并向非均匀网格推广.选取了典型算例对方法进行验证.一维激波管算例表明,改进的方法在保持计算结果高分辨率的同时,能够更好地抑制非物理振荡.分析得知,当人工粘性方法用于二阶DG格式时,所得计算结果的数值耗散较大,而当格式精度大于二阶时,采用人工粘性方法所得的结果的分辨率较高.通过计算圆柱高超声速粘性绕流,将三阶DG格式与三阶MUSCL格式和五阶WENO格式的结果进行对比,结果表明,该人工粘性方法对于高超声速流动计算也具有一定的优势.     

基于本征正交分解和代理模型的流场预测方法

 为了实现流场的快速求解,基于本征正交分解(POD)与代理模型提出了一种全新的流场预测方法。其基本原理为:首先采用本征正交分解将一定数量的样本流场分解为同等数量的基模态流场;然后用少数包含了绝大部分样本流场特征的基模态流场拟合所有的样本流场;最后用代理模型建立起决定样本流场的输入参数与拟合系数之间的近似函数关系。针对几何外形不同的二维翼型定常流场预测结果表明:在亚声速情况下,预测误差收敛的模态数量不超过20个,继续增加模态的使用数量不能明显提高预测精度;在跨声速情况下,预测误差收敛的模态数量为26个,当使用的模态数量达到前10个时,继续增加模态的使用数量能提高绝大部分流场区域的预测精度,但同时会在激波附近引入"噪声"激波特征而降低该区域的局部预测精度。在这两种情况下,预测流场所需时间均不到高精度计算流体力学(CFD)方法的1/200。     

An empirical model to predict the 1-AU arrival of interplanetary shocks

We extend the empirical coronal mass ejection (CME) arrival model of Gopalswamy et al. [Gopalswamy, N. et al. Predicting the 1-AU arrival times of coronal mass ejections, J. Geophys. Res. 106, 29207, 2001] to predict the 1-AU arrival of interplanetary (IP) shocks. A set of 29 IP shocks and the associated magnetic clouds observed by the Wind spacecraft are used for this study. The primary input to this empirical shock arrival model is the initial speed of white-light CMEs obtained using coronagraphs. We use the gas dynamic piston–shock relationship to derive the ESA model which provides a simple means of obtaining the 1-AU speed and arrival times of interplanetary shocks using CME speeds.     

高背压进气道中内外流耦合作用的大涡模拟

采用大涡模拟研究了出口堵塞比为50.8%的轴对称进气道流动,重点考察了内外流耦合作用下流动的非定常特性。采用国家数值风洞（NNW）工程仿真软件进行数值模拟,得到的壁面平均压力、瞬时压力分布与试验数据符合良好。分析表明：为匹配出口背压,进气道在喉道区域形成激波串结构,使内流道流场分为上游超声速区、中部激波串区以及下游亚声速区;在激波串区,剧烈的逆压梯度产生了分离激波、激波串、分离区及分离剪切层等复杂结构;伴随着激波串运动和边界层大尺度分离,进气道壁面压力出现宽频脉动特征。脉动压力的时空分布表明：内流道脉动压力以扰动波的形式传播,为此建立的声反馈模型能较好地预测亚声速区的主导频率。相关性分析表明：激波串运动受上下游流动耦合作用,其中,频率为St=0.7的运动主要受上游流动影响,频率为St=0.9的运动主要受下游压力扰动波影响。     

一种确定模拟的行星际MHD激波局地参数的新方法

简要阐述了分析模拟的行星际磁流体力学（MHD）激波的局部性质时，采用无厚度局部平面激波这一假设的合理性，说明了在激波未扰动区域（激波上游），物理量在几个小时内的变化很小这一事实，利用平面激波的分析方法，提出了分析模拟的行星际MHD激波的新方法，包括激波位置的确定，上下游状态参数的选择，激波局部参数的计算以及激波的分类，最后应用这种方法对一个二维的MHD模拟结果进行了分析。结果证实了过去文献关于磁流体力学混合激波空间连接和时间演化的链式规则，而且说明位于太阳赤道附近的慢激波和中间激波最终会发展为快激波。     

冲击与气膜的组合形式对冷却效果的影响

通过数值模拟,研究了涡轮叶片弦中区所采用的新型双层腔冷却结构的冷却特性,系统分析了冲击与气膜的组合冷却流动换热的机理,讨论了冲击孔与气膜孔的组合形式对组合冷却效果的影响.计算参数范围是:吹风比M=0.6~2.0,冷气进口雷诺数Re=2000~5000.计算结果表明:①气膜孔与冲击孔的位置及其排列方式对双层腔结构的冷却效果的影响是非常明显的,且存在一个最佳的组合冷却形式;②在狭小的封闭空间内,冲击靶面的努塞尔数分布呈明显的双峰结构,冲击滞止点处于两个峰值之间的峰谷.      

用GAO-YONG湍流模式对翼型跨音粘流的数值模拟

应用GAO-YONG可压缩湍流模式,数值模拟了NACA0012,RAE2822翼型的定常跨音速粘流算例.对流项采用三阶ROE格式,扩散项采用二阶中心格式,用多步Runge-Kutta显式时间推进法求解空间离散后的控制方程.计算结果很好地预测了翼型表面的压力系数的分布、激波的位置、马赫数等值线的分布等情况,并且对翼型表面激波与边界层相互干扰以及层流向湍流的转捩问题进行了分析计算.计算结果与实验值符合很好,表明GAO-YONG可压缩湍流模式应用合适的计算方法能够高精度模拟翼型跨音粘性流动问题,并且基于GAO-YONG可压缩湍流模式各向异性湍流粘性的机理,提供了一种预测转捩起始位置的判别方法.      

航空铝合金铣削加工中切削力的数值模拟研究

为了弥补当前斜角切削数值模拟多采用直线刀刃的不足,结合立铣加工的实际情况,提出了适合立铣加工的螺旋齿单刃斜角切削有限元模型,进而对航空铝合金7050-T7451进行了铣削加工切削力的数值模拟研究,得到了切削力值。通过铣削力实验测得了同样切削条件下的铣削力值,数值模拟结果与实验值比较吻合,从而证明所建立的有限元模型是正确的,可用于预报铣削力值。铣削加工切削力的数值模拟研究为航空铝合金切削加工的工艺参数优化、刀具的合理选择及其优化设计奠定了基础,同时也为进一步有效控制整体结构件的加工变形提供了新的研究手段。