基于四叉树切削网格的N-S方程求解方法

以四叉树非结构化网格为基础,提出了背景区域采用正方形四叉树网格、边界区域采用切削网格的一种可以表达复杂几何形状的网格生成方法,该网格具有生成过程简单,正交性好等优点.在这种网格的基础上,采用非结构网格有限体积法进行离散得到了多种形状切削网格并存时Navier-Stokes（N-S）方程的求解算法,并以顶盖驱动斜方腔流和方腔内热圆柱自然对流为例,应用上述算法实现了网格生成和流动数值模拟,与基准解进行了比较,一致性较好.计算结果表明这种切削网格方法及其N-S方程求解方法具有可靠性和应用前景.     

一种三维切削网格上的N-S方程离散方法

针对一种计算流体力学的三维切削网格自动生成方法Kitta Cube,提出一种新的离散方法(6+N模型)来描述和处理Kitta Cube中的任意三维切削网格.该方法能避免对各种类型的切削网格进行穷举分析,实现采用统一模式在三维切削网格上离散和求解N-S方程.用该方法对三维同心球和偏心球之间环形通道内的自然对流换热进行模拟,数值结果表明:通过Kitta Cube能较精确地表达曲面,面积误差在1.0%以内;利用6+N模型能在该网格上实现复杂区域内的流动和传热模拟,其数值精度与贴体网格相当,误差在0.5%以内;降低偏心球环形通道内部球心位置,能增加自然对流的传热性能.     

一种求解N-S方程的自适应直角网格方法

提出了一种用直角网格表达背景、切削网格表达边界的非结构化自适应直角网格方法.该方法采用四叉树保存网格数据,将切削简化成6种类型,用速度的旋度和散度作为自适应加密标准,从而可实现任意二维区域网格的自动生成和自适应加密.通过将极小网格边界化处理,利用SIM-PLE算法处理速度和压力的耦合,实现了该网格上N-S方程的离散和求解.算例表明,该方法网格生成简单,可以用于任意形状上的流动和传热模拟,相比非自适应方法,用一半的网格数目即可达到相同的计算精度.     

不同大气稳定度下绿化街谷流动的大涡模拟

耦合大涡模拟与植被阻力、热源项模型, 提出了城市灌木绿化街谷内风场和污染物浓度分布的热、动力数学模型. 首先利用均匀植被层边界层湍流算例考核了模型的合理性以及程序的稳定性, 表明本文模型能较好模拟植被层内部速度分布, 并能对温度变化做出正确响应. 对形状因子为0.5 的理想绿化街谷, 研究了街谷大气不同稳定度对街谷内风场和污染物浓度分布的影响. 结果表明, 相比于裸露街谷, 植被层的引入减弱了街谷内环流风场强度, 同时减小了街谷顶部污染物置换速率(pollutant exchange rate, PER), 导致绿化街谷地面、背风面和迎风面附近污染物浓度增加. 对于街谷大气稳定度对绿化街谷内风场和污染物浓度分布的影响, 街谷大气强不稳定情形下, 街谷顶部、地面、背风面和迎风面附近, 风速比街谷大气中性层结时明显降低, 污染物浓度相应增加; 街谷大气弱不稳定情形下, 街谷顶部和迎风面附近, 风速有轻微减弱, 同时, 街谷地面和背风面附近, 风速有加强趋势, 导致绿化街谷内污染物浓度在迎风面和背风面比中性层结时减小, 而在两排植被层之间有所加强. 关于街谷大气稳定度对街谷顶部PER 的影响, 街谷大气弱不稳定情形街谷顶部PER 要比街谷大气中性层结时小, 而街谷大气强不稳定情形街谷顶部PER 最小.