螺旋桨飞机滑流非定常数值模拟研究

采用非定常数值模拟方法,基于滑移网格技术,分析研究螺旋桨滑流对全机气动特性的干扰影响。用多块结构网格数值模拟软件,从搭接边界建立、搭接网格生成、流场插值传递3个方面对螺旋桨滑移搭接网格的数值模拟方法进行全面详细的介绍,并且阐述基于该技术的非定常数值模拟方法。在此基础上,开展螺旋桨飞机模型验证计算,分析研究螺旋桨滑流对全机流场和气动性能的影响,通过与风洞实验数据的对比分析,验证滑移搭接网格技术的可行性与准确性。研究结果表明,基于滑移搭接网格技术的螺旋桨滑流数值模拟,可以真实反映螺旋桨滑流流场特性,为螺旋桨类飞机的设计评估提供一种有效的数值模拟手段,在工程方面具有重要的实用价值。     

氧化层非均匀增长对热障涂层应力分布的影响

利用扫描电镜对不同氧化程度的热障涂层(TBCs)进行观察,分析发现热生长氧化层(TGO)各部位并非均匀增长。假设热生长氧化层的形状为正弦曲线,用统计方法处理扫描电镜图片可以得到热生长氧化层的基本几何参数。由此建立单波形二维轴对称有限元模型,计算并分析了非均匀增长模型的应力状况。结果表明:相同温度下,非均匀增长模型的轴向、剪切和Mises应力分量的变化趋势与均匀增长模型的对应趋势一致;但是非均匀增长模型的各个应力分量最大值均小于均匀增长模型的应力分量最大值。      

旋转弹体马格努斯效应数值模拟方法研究

旋转的弹体通常会存在马格努斯效应,严重影响飞行稳定性和弹道轨迹。提出基于旋转壁面法的多对称面弹体旋转效应数值模拟方法,无需借助动网格和多参考系模型,仅需要将旋转速度附加壁面,即可使用定常方法模拟马格努斯效应;同时,使用SOCBT标准弹体模型对该方法进行验证分析。结果表明:计算值与试验值吻合较好;由于旋转引起的附面层堆积以及大攻角分离效应、激波的干扰,计算时附面层的网格节点要足够密,并且应该针对不同的攻角选择精度与效率兼顾的计算模型。     

多电极偶接对金属大气腐蚀影响的试验与仿真

搭建薄液膜厚度测量与控制装置,采用微距参比电极后置法组建三电极体系,测量了厚度为100 μm的3.5wt% NaCl液膜下2A12铝合金、TC18钛合金和30CrMnSiNi2A高强钢的极化曲线和耦合体系的电偶电流;建立了基于薄壳电流分布的Comsol腐蚀仿真模型,得到了100 μm液膜下电极电位分布和电流密度分布图,通过对电极表面进行局部电流面积分计算得到了其电偶电流,对不同电极进行了面积的参数化扫描,讨论分析了不同面积比对于电偶电流的影响。结果表明：利用仿真模型得到的电偶电流值与试验值吻合较好,在2A12/30CrMnSiNi2A/TC18多电极体系中2A12充当阳极,TC18/30CrMnSiNi2A充当阴极,且阴阳极极性不随面积比的变化而发生转变,电偶电流主要集中在30CrMnSiNi2A与2A12之间,电偶电流与TC18面积呈线性关系,其对数与2A12和30CrMnSiNi2A面积的对数之间呈线性关系。利用建立的仿真模型对模拟搭接件大气腐蚀进行了仿真研究,仿真发现垫圈、2A12与30CrMnSiNi2A偶接接缝处及2A12与TC18偶接铆钉处电流密度分布更为集中,易发生腐蚀,与暴晒试验结果相符,说明该仿真模型对于多电极偶接的腐蚀分布具有很好的预测作用。     

细长体飞行器自由边界热模态试验与仿真

随着耐热承载一体化材料在新型高超声速飞行器上的应用,承力结构的工作温度不断提高,各类热模态特性逐渐得到关注。针对非平面形状的细长体飞行器自由边界条件下的热模态特性开展了研究。通过研究模拟气动加热条件的圆筒形加热笼、模拟自由边界的耐高温柔性支撑、非接触激光多普勒测振、耐高温激振杆激励等试验方法,获得了细长体结构自由边界条件下随温度变化的前3阶模态变化情况。结果表明：对此类薄壁长圆筒类结构,温度升高对模态频率影响可以超过6 Hz。开展有限元仿真,并与试验取得的热模态结果的变化规律进行对比。结果表明：建立考虑温度对结构弹性模量、热应力影响的壳单元模型,能够较好地预测出前3阶模态频率在全部受热时间范围内的最大下降量,可为高超声速飞行器控制系统设计时的拉偏范围提供参考。     

湍流数据同化技术及应用

近年来数据同化（DA）被引入湍流动力学研究中,通过融合实验测量和数值计算,提高了实验测量的精度和广度,改善了数值模拟的预测性能。实验观测、预测模型和同化算法是数据同化的三要素,湍流研究中的实验观测包括热线风速仪、粒子图像测速法（PIV）、压力传感器等局部测量数据,预测模型主要指流动控制方程及湍流封闭方程,而同化算法包括贝叶斯推断、集合卡尔曼滤波（EnKF）、伴随等。稳态数据同化一般结合雷诺平均Navier-Stokes （RANS）模型方程,从重新标定模型常数、修正涡黏模型方程形式误差、修正雷诺应力项等方面着手。非稳态的数据同化一般包括四维变分（4DVar）等时间连续的数据同化方式以及顺序数据同化。4DVar通过时间正向和逆向积分迭代,存储量和计算量都非常大。顺序数据同化不需要时间逆向积分,可以在若干时刻对实验观测进行间断地植入,正向求解整个系统。另外,随着人工智能的飞速发展,湍流数据同化研究也向智能化迈进。对于纯数据驱动的湍流机器学习,其缺乏物理本质的约束,而基于物理信息的机器学习在物理本质上与数据同化是一致的。     

前缘钝化和化学非平衡效应对斜爆震发动机进气道特性的影响

为了探究前缘钝化、化学非平衡效应对斜爆震发动机进气道的工作性能及出口温度边界层分布的影响,采用热完全气体、化学非平衡气体两种模型对不同顶板、唇口前缘钝化半径下斜爆震发动机进气道进行数值模拟,结果表明：相比基准进气道,钝化后进气道上壁面温度边界层厚度较大,下壁面温度边界层厚度较小;在化学非平衡气体模型下,顶板前缘钝化半径(R1)>2mm时进气道的顶板附近和分离区内离解反应较为明显,唇口前缘钝化半径(R2)>2mm时进气道的唇罩、唇口板附近离解反应较为明显;当钝化半径≥4mm时,两种气体模型下进气道出口总压恢复系数和静温的相对变化量绝对值大于0.5%,有必要考虑化学非平衡气体效应对进气道出口性能的影响。     

侧板边界层对正激波/边界层干扰特性影响的数值模拟研究

为了研究外压式进气道处于临界工况时结尾正激波同时与压缩面和侧板上边界层的相互干扰,专门设计了三组具有不同前伸长度侧板的简化构型作为研究对象,利用数值模拟手段评估了侧板边界层厚度对正激波/边界层干扰特性的影响。仿真结果表明,当没有侧板边界层参与干扰时流动呈现较好的准二维特性,但当侧板边界层参与干扰后将形成较强的角区干扰结构,该角区三维干扰结构与对称面上的准二维干扰结构存在此消彼长的关系。此外,波系的空间形态也将由“准二维λ波”结构变为“双λ波+角区压缩波”结构,波系形态的改变则进一步导致壁面回流区分布以及摩擦力线拓扑结构变得更加复杂。     

层流平板摩擦阻力的数值计算

为提高层流摩擦阻力数值计算的精度,采用数值求解N-S方程的方法,对层流平板摩擦阻力进行了计算.计算中采用4种先进上风格式,同时配合以14种不同的壁面法向网格分布和5种常用限制器,并将计算结果与布拉修斯解进行对比.结果表明,van Leer的FVS(Flux Vector Splitting)格式数值耗散大,对网格分布和限制器敏感,不适合用于摩擦阻力计算.限制器和计算格式的粘性分辨率越高,摩擦阻力的计算精度越高.在计算量允许的情况下,可以减小壁面法向第1层网格高度,以提高计算精度,但应同时兼顾网格长宽比.对于摩擦阻力计算,边界层内应至少布置4个网格点.      

基于水文剖面的三维地质建模方法

鉴于剖面是水文地质中常见的、有价值的数据形式,提出了基于水文剖面的三维地质建模方法,即"剖面信息提取",并给出了建模步骤.分析了基于剖面建模应该注意的问题和需要的预处理.详细介绍了剖面概化、二维剖面映射到三维空间、建立背景网格以及校正错误剖面信息的方法.设计了一种"顶面—尖灭线—透镜体"的三步建模法构建透镜体模型,该方法使复杂的多值建模问题转化为单值建模问题.描述了单值层面和透镜体的建模算法.最后给出应用实例,证明了基于剖面建模方法的有效性.