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

搭建薄液膜厚度测量与控制装置,采用微距参比电极后置法组建三电极体系,测量了厚度为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个网格点.      

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

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

武汉上空Na层季节和夜间变化的观测研究

根据武汉大学激光雷达2001年3月到2002年10月的观测数据,对武汉上空(30.5°N,114.4°E)的背景Na层结构进行了研究,并与其他观测站尤其是Urbana(40°N,88°W)的观测结果进行了比较。发现武汉上空11月的Na层柱密度最大,是4月的2.5倍,质心高度的极小值出现在5月和11月,极大值在4月、9月和12月;全年平均Na层均方根宽度为4.47km.Na层各参数呈现出周期性变化规律,其中Na层均方很宽度呈现出明显的半年周期性变化规律。夜间Na层柱密度呈现总体增加的趋势,质心高度在入夜后立即迅速上升,在2000LT以后上升速度趋缓,凌晨0500LT以后出现下降;Na层均方根宽度在夜间约有1.3km的起伏变化,它入夜后迅速减小,在午夜0035LT时出现一低谷,然后随时间逐渐变宽。     

运用爆破震动来识别地下裂隙和软弱夹层

从爆破震动信号主要来源于爆破区邻近范围内地质结构的自振的观点出发,提出了一种探测地下裂隙和软弱夹层的新技术:把震源布置在裂隙和软弱夹层附近进行爆破,由于介质中裂隙和软弱夹层的存在使得局部结构的刚度减小,直接影响到爆破介质振动的速度幅值和频率,因此可以通过地表接收器接收的振动信号(包括振动速度的幅值和频率)识别出裂隙和软弱夹层的大小和形状等特征.通过无网格的数值算例分别从水平裂隙和斜裂隙的情况下总结规律说明了这种方法的可行性.     

时域有限差分法连接边界电磁泄漏

用时域有限差分法计算目标的雷达散射截面时,一般用连接边界来引入平面入射波.理想情况下,当总场区没有散射目标时,该区域仅有入射波,散射场区电磁波为0.但在实际计算过程中,散射场区的电磁波一般不会严格等于0,这是因为在连接边界引入入射波时产生了电磁泄漏.一维情形下,用散射场区电场的平方和来衡量电磁泄漏程度.二维情形下,用等效原理将散射场区的电磁场进行远场外推,得到雷达散射截面,以此衡量电磁泄漏的大小.研究表明:时间步长、入射角度都能影响电磁泄漏大小.为使电磁泄漏较小,时间步长应接近于稳定性要求的最小步长,入射方向应避免垂直于计算区域边界.