含内热源的多孔方腔自然对流非正交MRT-LBM模拟

为了增强多孔方腔内流体流动与传热效果,采用非正交多松弛格子Boltzmann方法（MRT-LBM）对含有内热源的多孔方腔自然对流传热现象进行了数值模拟。研究了不同冷源布置方案（Scheme A~Scheme F）、内热源结构形式（Case 1、Case 2、Case 3）、内热源位置（a、b）、Darcy数、Rayleigh数等对多孔方腔内流体流动与传热的影响。计算结果表明:冷源布置方案对腔内流体流动与传热具有重要影响,当冷源左右对称布置时,腔内温度场及流场亦对称分布;在高Rayleigh数下采用Scheme A的双上部冷源布置方案能明显提高腔内的传热强度;内热源的形状对腔内对流传热影响很大,高Rayleigh数下,Case 3的布置方式更好。内热源的位置a和b对腔内的传热影响明显,提出了热壁面平均Nusselt数与位置a的拟合关系式,存在最佳的位置a（a=0.25）,使得腔内的对流传热最强;热壁面平均Nusselt数亦随b值变化表现出特定的变化规律。随着b值的增加,热壁面平均Nusselt数呈现先增后减再增的趋势。      

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

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

轴向通流旋转盘腔流动换热稳态实验与数值模拟

通过开展稳态实验及数值模拟探究了轴向通流旋转盘腔的流动结构与换热特性。通过改变流量系数、旋转雷诺数等参数,探究不同工况下旋转盘两侧及盘罩内侧壁面温度和努塞尔数的径向分布规律。结果表明:在轮缘加热的状态下,旋转盘两侧温度径向分布均呈凹函数形态,且旋转盘迎风面换热强度普遍高于背风面;后轴颈盘罩向两端旋转盘导热,其壁面温度径向分布呈"中间高、两侧低"的状态;随着轴向流量系数的增大,盘腔内部气体对流加剧,径向臂及涡对结构更加明显,旋转盘及轴颈表面换热效果增强;旋转盘腔内的流动换热特性受强迫对流和类Rayleigh-Benard对流2种机理的共同影响。      

PANS方法在双圆柱绕流数值模拟中的性能分析

PANS方法以传统的雷诺平均Navier-Stokes(RANS)方程为模板,引入控制参数对原方程进行修改,是一类能够有效模拟分离湍流流动的混合方法.采用基于Menter 剪切应力输运(SST)湍流模型的PANS方法,选取双圆柱绕流作为研究对象, 分析双圆柱流动结构,验证方法的准确性与可靠性,考察PANS方法中模型参数选取对于计算结果的影响,评估不同类型PANS方法的性能.研究表明:模化湍动能比例这一模型参数对计算结果有显著影响,且计算域内统一的模型参数取值难以处理复杂流动问题.在此基础上使方法中模化湍动能比例取值随流场物理信息和网格尺度变化,所得结果与SST DES计算结果及实验结果吻合较好,表明该方法适用于复杂湍流数值模拟.      

基于CE/SE法的翼型流场激波捕捉计算

在二维四边形结构网格下,针对时空守恒元和解元(CE/SE, space-time Conservation Element and Solution Element)方法捕捉激波需要双倍空间网格点的问题,重新对守恒元和解元进行了定义,计算点同时包含单元中心和网格节点,并推导得到了新的CE/SE方法的计算公式.以此为基础,结合非敏感克朗数(CNIS, Courant Number Insensitive Scheme)计算格式和当地时间步长方法,对激波翼型流场进行数值模拟,并与原方法及AGARD报告进行对比.结果表明,对原CE/SE方法的改进是有效的,可明显提高激波分辨能力,并且激波前后无明显的数值震荡发生,新方法适合应用于翼型流场中捕捉激波.     

Von Karman模型三维大气紊流仿真理论与方法

基于三维紊流场相关函数矩阵,用蒙特卡罗法仿真生成Von Karman模型三维空间大气紊流场,实现了对复杂的Von Karman模型的精确推导和应用,仿真流场的相关特性在理论上可以任意逼近Von Karman理论模型.通过提高白噪声序列的白化程度及优化蒙特卡罗算法的方法提高了仿真精度和效率.数值仿真结果证明:所生成的空间大气紊流场符合模型特性,具有较好的统计特性和仿真精度,满足实时飞行仿真对紊流场数据的要求.