模拟高度复杂流动的自适应网格算法

提出一种基于对等弧分布律进行扭转修正的自适应网格算法。该方法将三维网格系统分裂为３个单方向自适应网格子系统,在计算中与主控方程的求解格式相耦合,交替地使用。应用所构造的算法结合ＭａｃＣｏｒｍａｃｋ显式格式对喷流／绕流干扰流场进行了雷诺平均Ｎ－Ｓ方程模拟,给出了清晰的剪切层与激波波系结构。计算结果表明,该自适应网格算法可较大程度地提高ＭａｃＣｏｒｍａｃｋ格式的分辨率。     

关于建立高阶差分格式的问题

为了能在不太密的网格上捕捉到流场的细致结构。通常采用高阶精度的差分 式进行数值模拟。为能抑制计算中和在激波附近产生的虚假波动,本文从物理构思出发,提出了建立高阶格式的两个基本原则,作者称之为抑制波动的原则和稳定性原则。     

不同流速比超声速混合层气动光学效应研究

针对不同流速比混合层流场气动光学效应问题,首先采用大涡模拟数值方法进行了数值仿真,其次用光线追迹方法进行了气动光学效应仿真,最后对混合层的气动光学效应进行了评价。结果表明,混合层流场涡结构与光程差极小值存在一一对应关系,并且其流速比低的混合层气动光学效应更加严重。同时,研究了不同光学参数对斯特列尔比的影响,得到了相应的规律。     

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

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

一种高光谱图像的双压缩感知模型

高光谱图像因其海量数据性,给存储、传输及后续分析处理带来了挑战。压缩感知理论提供了一种全新的信号采集框架。针对高光谱数据的三维特性,提出一种双压缩感知的采样与重构模型。该模型在采样阶段兼顾高光谱数据的空间和谱间稀疏特性,构造了能同时实现空间和谱间压缩采样的感知矩阵;重构阶段不同于传统的压缩感知重构方法直接重构高光谱数据,而是将高光谱数据分离成端元和丰度分别进行重构,然后利用重构的端元和丰度信息合成高光谱数据。实验结果表明,所提双压缩感知在低采样率下重构精度较三维压缩采样提高了10 dB以上,更为显著的是运算速度提升了3个数量级,同时该方法还便于获得端元和丰度信息。     

波瓣数对波瓣S型混合二元喷管气动热力性能影响

依据某型涡扇发动机波瓣S型混合二元喷管,保持波瓣混合器长度、内扩张角、外扩张角以及宽高比不变,依次取波瓣混合器波瓣数为12,14,16,18,20,建立了一组具有不同波瓣数的波瓣S型混合二元喷管模型.采用经过验证的CFD方法,研究了波瓣数对波瓣S型混合二元喷管气动热力性能的影响规律.结果表明:在波瓣尾缘截面至第1个S弯截面区域,波瓣数对流体混合程度产生很大影响,并且热混合效率近乎随波瓣数增加而增加.在第1个S弯截面至波瓣S型混合二元喷管出口截面区域,波瓣数为16的波瓣S型混合二元喷管模型的总压恢复系数始终最低,其余模型的总压恢复系数以及热混合效率没有明显差别.在波瓣S型混合二元喷管出口截面上,波瓣数为16的波瓣S型混合二元喷管模型的热混合效率最高,达到0.850,然而其总压恢复系数相对于该截面上最高值下降了0.289%.此外,波瓣S型混合二元喷管的渐缩型流道能够提高流向涡强迫混合效果,但同时也加速流向涡的耗散速率.      

超声速尾迹-剪切流的混合增强

超声速条件下燃料和空气之间的高效混合是超然冲压发动机技术上的主要挑战。基于大涡模拟和流动稳定性分析,针对超声速尾迹-剪切流动开展了混合增强方法研究。尾迹的存在改变了混合层流动的速度剖面,对流动稳定性产生了重要影响,使混合层由三维最不稳定变为二维最不稳定,最不稳定扰动波频率和增长率增大。基于流动稳定性结果引入扰动的混合增强方式依然有效,根据稳定性结果设计了波纹隔板。数值结果表明：二维波纹壁引入的扰动未能增长,不具备混合强化效果,而三维波纹壁引入的扰动能够快速增长,具有混合强化效果,且波纹壁参数越接近最不稳定扰动波参数,混合强化效果越明显。     

二维微槽道电渗流的增强混合

微尺度流动中的混合一直是微流系统研制中的一个大问题,对壁面zeta势非均匀分布情况下的二维微槽道电渗流进行了数值分析,重点分析了这种情况下的微流动混合增强.基于电渗流的Helmholtz-Smoluchowski滑移速度模型,应用SIMPLE方法计算了4种滑移速度在整个壁面上的分布不均匀一致时的微槽道电渗流动的流场及被动标量在其中的混合.计算结果表明,滑移速度的非均匀一致导致微槽道电渗流中的流线发生剧烈的转折,从而使得微槽道流动的混合效率有很大的提高.流场中流线的转折越剧烈、转折的次数越多,微槽道流动混合效率的提高就越大.     

超声速复杂流动中湍流模型的性能评估

对8°,16°,20°,24°超声速二维压缩拐角进行了数值模拟,系统评估了SA(Spalart-Allmaras), k-ω, SST(Shear Stress Transport)3种工程常用湍流模型对激波/边界层干扰复杂流动的模拟适用能力.对于"弱"干扰、小分离,工程常用湍流模型能够较准确地预测其压力、摩擦力、热流分布,而对于"强"干扰、大分离则会造成较大的偏差.同时,随着分离区的增大,各湍流模型的适用性能差别更加明显,数值方法也有一定的影响.      

基于边界层转捩的高超声速进气道特性研究

为了探索边界层非强迫转捩对进气道性能的影响,采用数值计算的方法开展了边界层转捩对轴对称混压式高超声速进气道流场特性的研究。研究表明:随着进气道中心锥锥尖钝化半径增大,边界层转捩先推迟。当锥尖钝度大到一定程度时,边界层转捩位置前移。随着钝化半径进一步增大,边界层转捩再次推迟,转捩位置逐渐后移。来流湍流度越大,边界层越不稳定,边界层转捩越易发生。与湍流边界层相比,考虑边界层转捩时进气道的总压恢复系数及流量系数较高、热载荷及阻力系数较小,Ma=6.5时喉道处总压恢复系数最高上升17.3%,进气道阻力最大下降17.4%。边界层转捩对壁面热流密度分布影响较大,但对壁面压力分布影响较小。钝化影响进气道的自起动性能,随着钝化半径增大,自起动马赫数升高,而边界层转捩对进气道自起动性能影响较小。