三维可压缩 Navier-Stokes 方程的间断Galerkin 有限元方法研究

拓展了二维间断 Galerkin(DG)有限元方法研究,将该数值方法用于三维可压缩欧拉方程和 Navier-Stokes方程的求解。基于六面体网格单元,采用插值方法将物面的四边形面网格单元构造为弯曲面网格单元,更好地表述了真实物面特征;物面边界相邻体网格单元相应构造为高阶体网格单元,其余体网格单元采用八节点六面体单元,以较小的计算代价使网格满足 DG 方法计算需求。通过对三维带 bump 管道内流、圆球绕流以及旋转流线体绕流进行的数值求解,验证了边界弯曲方法的可行性及 DG 方法的高精度特性。此外,由于采用了隐式计算方法,仅需较少的时间步就能迭代收敛。     

萤火一号火星探测器科学应用数据处理技术

对萤火一号(YH-1)火星探测器科学应用数据处理技术进行了研究.给出了用磁场数据确定磁堆积区边界(MPB)特性的最小方差分析(MVA)、离子数据各质量通道权重修正和利用星星掩星数据确定火星电离层电子密度等数据处理技术,有助于进一步了解火星空间环境的本质.     

抑制恒压变量柱塞泵压力低频振荡方法的研究

为解决恒压变量柱塞泵出口压力出现低频振动的问题,阐述了恒压变量机构的原理,从理论上分析了影响出口压力的因素,并指出产生振荡的原因,提出在泵装配过程中,采用工艺关键要点控制的方法可抑制压力低频振荡。试验表明,该方法有效地减小了泵出口压力的低频振荡,使产品装调合格率大大提高。     

间隙非线性气动弹性系统颤振及控制问题研究进展

含有间隙结构的气动弹性系统非线性颤振问题是飞行器气动弹性力学工程领域的研究热点和难点。根据目前现代飞行器结构轻量化设计及更大机动性能的发展趋势,非线性颤振问题日益突出,直接关系到飞行器的安全与性能。因此综述了近几十年来带间隙非线性的非线性气动弹性力学模型、非线性系统辨识及非线性动力学与控制等问题的研究进展。在已有相关研究成果的基础上提出了今后值得进一步解决和关注的研究问题。     

高效率的特征型紧致WENO混合格式

特征型紧致加权基本无振荡（WENO）混合格式HCW-R结合了迎风紧致格式CS5-P和WENO格式,具有十分优异的分辨率特性。但在求解多维方程组时,HCW-R格式需要求解块状三对角方程组,因而计算代价十分高昂。采用迎风紧致格式CS5-F代替CS5-P,构造了一个新的特征型紧致WENO混合格式HCW-E。由于HCW-E的特殊形式,其可沿迎风方向、由边界处向内推进求解,避免了处理三对角或块状三对角方程组,从而其计算代价与显式格式无异。虽然就分辨率而言,HCW-E稍逊于HCW-R,但前者的计算效率要显著高于后者。因此,当花费相同的计算代价,HCW-E格式可以获得更好的数值结果。一系列求解Euler方程组的数值试验验证了HCW-E的高分辨率特性和相比HCW-R更高的计算效率。HCW-E格式的效率优势在求解高维问题时更为明显。      

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

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

短螺旋型燃烧室旋流结构及气动边界发展规律数值模拟研究

短螺旋型燃烧室的头部轴线沿发动机周向与发动机轴线偏转一定角度,可以在保证燃烧效率不变的条件下有效缩短燃烧室轴向长度。短螺旋型燃烧室流场的最大特点是旋流流动单侧受限。为了研究头部安装角α变化对燃烧室旋流流动特性的影响,基于数值方法对短螺旋型燃烧室进行了计算分析。结果表明:随着α变化,旋流器下游旋涡依次出现对称环状、马蹄状、环状结构;随着α增加,气动边界逐渐出现并抑制旋流的周向扩张,导致流场出现不同的旋涡结构;不同α下的切向角动量随轴向距离增加而衰减,但α为35°和45°时,气动边界在非受限侧出现并对旋流产生约束,角动量衰减变慢;当α为0°,15°,35°,45°时,燃油液滴依次集中分布在旋流器下游两旋涡边缘、侧壁面和头部端壁、非受限旋涡边缘。本文研究了不同α下的旋涡结构及气动边界沿轴向的演化过程,为短螺旋型燃烧室进一步的设计与优化提供基础。     

过渡型凹腔超声速燃烧流动振荡抑制方法

针对长深比为10.0的过渡型凹腔在隔离段入口马赫数为3.0条件下存在的冷流自激振荡现象,提出了一种凹腔内增加肋条抑制振荡的方案。通过试验和数值计算,对该方案抑制振荡的效果进行了检验,并分析了肋条增加前后燃烧室流场结构和燃烧性能的差别。研究发现：通过在凹腔内增加肋条能够消除过渡型凹腔冷流工况下存在的175.8 Hz的自激振荡,燃烧流场也更加稳定;增加肋条后凹腔的稳焰能力有所降低,部分在凹腔未完全燃烧的煤油进入扩张段后继续发生反应,从而使燃烧区向下游延伸、增大,发动机的燃烧效率和净推力分别降低5.4%和8.9%,但推力更加平稳;燃烧室一维平均热流密度峰值由2.9 MW/m²降低至1.8 MW/m²,燃烧室的热环境大幅改善。     

新电位移边界条件下的压电介质断裂问题

应用边界元法,并结合一种特殊的复变函数基本解,研究了压电介质的断裂问题。由于该基本解完全遵循本构方程没有对电位移边界作假设,所以结论同非穿透性裂纹假设的结果有较大差异。通过计算表明,外加力载荷不但会引起裂纹尖端应力奇异性,也会引起电场奇异性,同样,外加电压也会引起裂纹尖端应力奇异性。     

蜂窝状薄壁管拉伸时的内陷分析

蜂窝材料因具备广阔的结构设计空间,已广泛应用于航空航天等工业领域。探索蜂窝元胞的组织结构、组分性质与材料整体性能之间的关系有助于新型功能性材料的设计开发。本文从非平面Vertex模型的势能形式出发,结合蜂窝薄壁圆管拉伸时的轴对称特征,通过变分法得到了圆管拉伸时母线满足的控制方程,揭示了边界效应是蜂窝状薄壁圆管受拉时产生内陷的原因,并结合母线控制方程的若干特解,考察了非平面Vertex模型势能形式中包含的材料性质参数与圆管拉伸后内陷程度及范围的关系。研究还指出,蜂窝材料的材料性质包含在非平面Vertex模型的势能形式中,为蜂窝状材料整体性能的研究提供了一种新的思路。