多重网格法求解二元跨音速欧拉流

 近年来,采用时间相关欧拉方程或N-S方程求解跨音速流得到了迅速发展。为了获得二元或三元复杂形体正确的表面气动参数和足够高分辨率的流场物理特性,必须采用十分细密的网格,因而导致计算时间成倍增长。为此,除继续发展各种差分格式,探讨各种自适应网格技术外,运用多重网格技术加速收敛过程得到了国外普遍关注并取得了一定的成功。国内也已有采用多重网格法解跨音速欧拉流的文章,然而均局限于两重网格的计算。为探索多重网格法和焓阻尼技术对加速收敛的作用,我们以NACA0012翼型的跨音速流动为例,采用MacCormack两步显式格式进行了数值实验。     

基于蒙特卡洛方法的冰晶撞击特性计算

航空发动机摄入冰晶易发生内部结冰、造成发动机推力损失,严重时引发航空事故。冰晶结冰数值模拟是研究冰晶结冰的重要手段,而撞击收集系数计算是冰晶结冰数值模拟的关键环节。基于NNWICE平台开展冰晶结冰前期的阻力系数、运动轨迹和撞击收集系数的数值模拟与计算研究,实现了基于蒙特卡洛方法的冰晶撞击收集系数的准确计算;在此基础上发展了拉格朗日框架下冰晶粒子运动、传热传质耦合数值计算方法;通过完成二维NACA0012翼型的相关算例,对比分析了不同阻力模型和冰晶相变对撞击收集系数的影响。结果表明阻力模型对冰晶撞击收集系数影响较小,冰晶相变过程会对冰晶运动轨迹及撞击收集系数产生显著影响,为进一步发展冰晶结冰数值模拟计算方法奠定了基础。     

正反顺序加载下安装边螺栓预紧力分散度的数值模拟与试验验证

在航空发动机机匣安装边装配过程中,不同位置螺栓之间会产生弹性交互作用,不同加载方案会使装配后的螺栓预紧力存在一定的分散度。针对这一问题,建立了3层带止口机匣安装边3维接触有限元模型,研究不同加载方案下螺栓预紧力的施加方法。基于有限元模型数值模拟了3层机匣安装边分别在顺序加载方案、交叉加载方案和正反顺序2轮加载方案下螺栓预紧力的分散度。搭建了相应的试验装置,分别对3种加载方案的预紧力分散度进行试验验证。结果表明：采用正反顺序2轮加载后的螺栓预紧力分散度最小,在3层止口机匣安装边上采用正反顺序2轮加载方案后的螺栓预紧力分散度较小,在有限元仿真中螺栓预紧力分散度为2.43%,在试验中螺栓预紧力分散度为4.98%。采用此方案加载后,止口机匣安装边能达到更好的密封效果和较好的刚度对称性。     

不同壁温及差分格式下超燃冲压发动机的仿真

针对圆截面超燃冲压发动机相关实验进行数值仿真,采用了不同的壁面温度条件和差分格式,并探究其对仿真结果的影响规律。结果表明采用常用的绝热壁面边界条件会导致激波边界层交互位置处的回流区变大,进而导致流道收缩,壁面压力升高。采用高精度格式不仅不能更为精确地计算回流区的大小,反而会使回流区面积变得更大。通过改变壁面边界条件,采用300 K的等温壁面,可以使激波边界层交互位置处的回流区变得极小,使仿真结果与实验结果十分吻合。计算结果表明,针对当前构型与实验条件,盲目提高数值格式的精度并不一定会得到更好的结果,相反可能会使仿真与实验更加偏离,而适当修改壁面温度条件,即使数值格式精度较低,仍然可以得到很好的仿真结果。最后,针对与实验结果吻合的仿真结果,分析了圆截面超燃冲压发动机的流场特性,重点研究了凹腔剪切层及其质量交换特性。     

守恒型可压缩一维湍流方法及其在超声速标量混合层中的应用

一维湍流（ODT）方法是一种能在一维计算域上遵循湍流基本物理规律的湍流建模方法。通过结合确定性和随机性求解方法,能够在一维计算域上准确捕捉到湍流统计规律,且降维建模可显著减小计算量。ODT方法主要被广泛用于不可压湍流和湍流燃烧研究,若要将其拓展用于模拟高速可压缩湍流,需对建模方法进行深度改进。相比于不可压ODT方法,本文基于欧拉参考框架,针对可压缩湍流的特性,将因变量由原始变量改为有利于减小可压缩湍流模拟误差的守恒通量,并加入了组分求解模块。对确定性和随机性求解模块均进行了相应的深度改进,开发出具有标量混合模拟功能的守恒型可压缩ODT方法。在确定性模块中改为求解以守恒通量为变量的一维截断控制方程,在随机性模块中构造一维涡时,将三联映射的作用对象也相应地由原始变量改为守恒通量,并选用了可保证变密度情况下动量守恒的双核变换。通过模拟空间发展超声速平面湍流混合层并将自相似阶段结果与实验结果比对,验证该方法对可压缩剪切湍流场中标量混合的捕捉精度。守恒型可压缩ODT方法模拟得到的速度场和组分场的平均剖面和脉动强度分布与实验结果准确吻合,精度明显优于传统的耦合梯度扩散亚格子模型的大涡模拟方法（LES-GRAD.DIFF.）以及耦合线性涡（LEM）亚格子模型的大涡模拟方法（LES-LEM）,且该方法的降维处理使其在降低计算成本方面具有显著优势。     

临近空间环境下封闭方腔内耦合换热特性

以临近空间浮空器载荷舱为应用背景,对复杂热边界条件下含热源的三维封闭方腔内自然对流、表面辐射和导热的耦合问题进行了数值模拟。综合考虑对流换热、长波辐射、太阳辐射等因素的影响,建立了临近空间热环境模型。通过Fluent软件用户自定义函数（UDF）引入外部非定常的辐射-对流耦合热边界条件,对腔内换热特性的昼夜变化进行研究,并分析了腔壁厚度、发射率和导热系数对其的影响。数值结果表明,腔内平均温度昼夜变化很小,约为12.9 K,但温度场分布随太阳方位变化而变化;腔内对流换热较弱,同一时刻最大温差约为71.3 K;腔壁热阻和发射率增加会削弱自然对流的强度。      

非均质壁面对液滴俘获能力的数值模拟研究

为研究不同非均质壁面对液滴的俘获能力,采用界面追踪法（FTM）结合广义滑移边界建立了接触角模型,对液滴在非均匀润湿的非均质壁面上的运动过程进行了数值模拟研究。液滴在倾斜壁面上受到重力作用由均匀润湿部分下滑至非均匀润湿部分,通过改变Bo数、Oh数、非均匀润湿程度研究了液滴在非均匀润湿区域的运动规律。研究表明:Bo数越大,液滴运动受壁面阻力影响越小,液滴下滑的速度越快,液滴越难以被俘获;Oh数越大,液滴运动受壁面阻力影响越小,液滴越难以被俘获;非均匀润湿程度越大,非均质壁面对液滴的阻力越大,液滴越易被俘获。      

低速大转矩电机冷却系统数值计算

在低速大转矩电机系列设计中,考虑到机座的通用性和外形尺寸的规范性,采用了轴向水路的水冷机座结构。通过对电机冷却系统的流体场计算,得到了电机的散热情况,耦合温度场计算,得到了电机冷却水和机座的温升,并通过对机座的模态分析得到了电机的固有振动频率。     

三结太阳电池高能电子损伤模拟及仿真分析

利用wx AMPS软件构建了高效Ga In P/Ga As/Ge太阳电池的中电池模型,并对电池抗辐照性能进行模拟研究。模拟发现,当辐照缺陷密度较小时,缺陷对中电池的电性能影响较小;当缺陷密度较大时,电性能的下降与电子注量值的对数成正比。计算电池的I-V和量子效率谱(QE曲线)可知,电池电性能的下降直接对应于量子效率的下降、饱和暗电流的增强以及并联电阻的衰降。模拟结果与试验结果的对比显示,在各子电池均匀损伤的假定下,1 Me V电子辐照的缺陷引入率约为0.81。