海岸TIBL的数值模拟──E─ε闭合模式的应用

利用二维的非线性ＰＢＬ方程组，采用工程上常用的Ｅ－ε闭合方法模拟了海岸地区热力内边界层的形成和发展过程中的背景个均风、温场以及湍流场的特征，分析了ＴＩＢＬ廓线变化规律并获得与实际一致的结果     

宽顶堰紊流场的数值模拟

本文采用k-τ紊流模式求解完整的Reynolds方程,模拟了有坎宽顶堰自由堰流.在求解中,发展了一套用边界拟合坐标、控制体积法求解二维恒定有自由表面紊流的数值方法,提出了分区光滑变换、反馈迭代等新概念.结果表明,本文提出的计算堰流的数值方法,计算量较小,具有良好的稳定性和收敛性,能够满足工程实际需要.     

有曲率影响的湍流流场中几种k—ε模型的 …

采用圆柱坐标系下简化的高阶各向异性ｋ－ε模型和各向异性ｋ－ε模型分别对有曲率影响的弯管内流流动进行数值分析,并与实验数据作比较,结果表明,由于ＭＡＫＥ模型考虑了曲率对雷诺应力的影响,对雷诺应力的预估更与实验值符合,对近壁区主流速度、雷诺应力的计算,ＭＡＫＥ模型比ＳＫＥ模型、ＡＫＥ模型有进一步的改善,证明ＭＡＫＥ模型预测有曲率影响的湍流流动行之有效。     

用非线性κ—ε模型对狭缝冲击射流进行数值计算

采用Ｓｐｅｚｉａｌｅ的非线性κ－ε模型和标准的κ－ε的模型计算了狭缝冲击射流,给出了径向平均速度和脉动速度沿壁面法向的工同实验结果作了对比。结果表明,当其他条件都相同时,线性和非线性模型在冲击射流的计算中相差不大,今后的改进应当在近壁区域对紊流模型本身进行。     

波浪为形数值模拟

从Navier-Stokes方程出发，采用修正的k－ε模型来封闭雷诺方程，作为紊波波浪场的控制方程，采用修正的VOF方法跟踪自由表面，应用施主与受主单元模型求解流体体积函数F的输运方程，在差分离散方法的基础上，采用通度概念，对不规则边界进行处理，在近壁区范围内，采用壁函数技术解决壁面效应问题，通过数值计算成功地模拟出椭余波 通过潜堤时波面的变形过程。     

数值模拟气泡在粘性液体中上升运动的Euler—Lagrange方法

在考虑界面张力的条件下,用边界层方法通过把液体粘性对无粘势流的影响化为修改势流在自由面上边界条件的办法,导出了数值模拟二维气泡在粘性液体中上升和变形的边界积分方法,势流计算采用Ｌａｇｒａｎｇｅ坐标,粘性修正采用Ｅｕｌｅｒ坐标,边界层近似的误差为δ的三次方量阶。数值模拟结果显示,在大气泡情形,粘性的影响很小,随着气泡尺寸的减小,表面张力和粘性的影响都增强,但与表面张力相比,粘性的影响处于比较次要的地     

连铸热工过程的数值模拟与热／力模拟研究

采用数值模拟方法计算连铸坯成形过程的温度场，将数值模拟与热／力 模拟方法相结合，分析力学性能随不同热工过程变化的动态演变过程。得到 无表面裂纹的铸坯表面温度模型为：表面温度为缓慢下降过程，并在低塑性 温区的上限温度９００℃进行矫直。     

戈壁夜间大气边界层的数值模拟研究

本文采用一维非定常边界层模式研究了戈壁下垫面夜间边界层的结构和发展规律，得到了戈壁下垫面的夜间冷却条件下辐射逆温层的形成和发展以及下波风建立的典型时空图象。并与实测资料进行了比较。文中还对局地逆向流和地面降温率对辐射逆温层高度和下坡风的影响进行了试验。     

黑河地区边界层平均结构和通量的数值模拟

利用ＮＭＣ２．５°×２．５°同化资料，采用改进的高分辨Ｂｌａｃｋａｄａｒ边界层模式，对１９９１年７月７日０８ＢＬＴ（北京时）至次日０８ＢＬＴ黑河地区边界层平均结构和通量进行了２４ｈ数值模拟，得到一些较为合理的结果，给出了不同下垫面上晴天和夜间的边界层结构特征和通量演变，边界层结构与下垫面和天气条件关系甚密，并与已有观测结果基本一致．     

基于RNG k-ε湍流模型钝体绕流的数值模拟

采用RNG k-ε湍流模型对钝体绕流流场进行了数值模拟．计算采用非交错网格的有限体积法求解二维不可压N-S方程．计算结果表明，RNG k-ε湍流模型可以成功地模拟绕钝体的不稳定、非定常和剧烈分离流动。     

用改进的RNG模式数值模拟湍流冲击射流流动

在Yakhot和Orszag提出的重正化群（RNG）k-ε湍流模式的基础上，对模式系数Cц和近壁处理方法加以改进，并将其推广应用于湍流冲击射流的数值模拟。将计算得到的平均速度分布和湍动能分布与标准k-ε模式的计算结果以及实验数据进行了比较。结果表明，用改进的RNGk-ε模式得到的数值结果与标准k-ε模式的结果相比具有显著的提高，尤其对滞止区附近的湍动能分布预测结果与实验结果具有较好的一致性，这是因为PNGk-ε模式与改进的近壁处理方法相结合有效地降低了涡粘性的预测值，进而引起湍动能预测值的下降，且计算中采用了Cц和Cε1的改进表达式，适当地考虑了流动中各向异性湍流和不平衡应变率影响的结果。该模式在工程湍流计算中具有较好的应用前景。     

基于标准k-ε模型的平衡大气边界层研究

平衡态大气边界层的精确模拟是计算风工程研究中的一项重要且容易被忽视的工作。文章基于标准k-ε模型的控制方程,对大气边界层的平衡条件作出适当假设,总结了与剪应力相关的入流边界条件的通用形式。为实现大气边界层的自保持,在湍动能和耗散率输运方程中添加源项使入流边界条件与湍流模型相协调,并对顶部和底部边界条件施加剪应力修正;然后,分别以二维缩尺和足尺模型为例进行了平衡大气边界层的模拟。研究发现,在湍动能和耗散率输运方程中添加源项对改善大气边界层的自保持产生了不可忽略的影响,在顶部和底部施加剪应力能较好地改善顶部和底部大气边界层的自保持效果,通过同时添加源项和施加剪应力能很好地实现大气边界层的自保持。     

山谷环流的二维数值模拟

建立一个采用地形跟随坐标系,能适用于复杂地形的边界层数值模拟的二维静力平衡大气模式．模式只需要单点探空资料及高层地转风资料作为初始条件．利用该模式对祁连山引起的山谷风环流进行模拟研究,试验结果表明在早晨有较强的下坡风,厚度为１０００ｍ;而在下午表现为谷风,厚度为８００ｍ．     

搅拌槽内三维流动场的RNGκ—ε数值模拟

RNGκ-ε模型在耗散率方程中通过系C1^*引入描述流场畸变效应的附加源项后，在一定程度上会改善对旋转流，浮力流等较复杂湍流的预报能力。本文应用该模型对六直叶涡轮搅拌桨的三维流动场进行了数值模拟，并将计算结果与实验数据进行了比较。计算结果表明：RNGκ-ε模型对桨叶附近速度场的预报较κ-ε模型有一定程度改善，但对湍流动能的预报却要比κ-ε模型差。若要进一步改装对桨叶附近流动场的预报，必须放弃基于各向同性假设的湍流模型，转而采用能够反映各向异性的模型或采用先进的计算方法。     

水平井数值模拟过程中的井筒动态特性分析

为进一步明确水平井数值模拟过程中井筒内因素对生产动态特性的影响 ,利用水平井数学模型和多年来进行数值模拟研究的经验 ,对影响水平井动态特性的井筒模型网格划分、井筒热损失、井筒压降等 3个因素进行了研究。结果表明 ,水平井井筒模型的不同网格划分决定了影响油藏和井筒的压力等参数分布的端部效应 ;井筒内的热损失对水平井井长优选、油井产量等均有影响 ,其中干度的敏感性分析表明 :干度提高 ,产液量和油汽比均有所提高 ,适当提高注汽速度将减少热损失 ;给出了对井筒内造成的压力损失的基本条件 ,可以采取不同井筒半径、射孔方案、完井技术等手段对压力损失进行控制     

竖壁自然对流的数值模拟

为了研究竖壁自然对流换热的特点,将边界层的控制方程转化为相似方程.在此基础上,进一步采用引入中间变量的方法,将相似方程转化为五个常微分方程.运用MATLAB对相似解进行数值模拟,获得了层流边界层内速度和温度的分布特点,同时,得到了与其相关的其它物理量的变化规律.     

溢流反弧段紊流边界层变化规律的数值模拟

在溢流反弧段流线同心圆假定的基础上,通过理论推导,建立了表征反弧段动水压强、壁面势流流速、壁面切应力的相应表达式;将描述边界层发展变化规律的连续方程、动量积分方程、动能积分方程联立,并运用列主元高斯消去法将其转化为一阶微分方程组,采用双边法求其数值解,从而模拟出反弧段水深、边界层厚度以及流速分布指数的沿程变化规律。     

筛板洗选效果的数值模拟研究

选用雷诺平均Navier-Stokes方程,采用标准κ-ε模型,定义专门的边界模型条件,运用流体分析软件对圆弧波浪式筛板和直条式筛板的洗选过程分别进行数值模拟,得出了原煤颗粒入流速度为2m/s时,两种筛板在洗选过程中的速度分布图、密度分布图,根据模拟结果分析,对传统的直条式筛板进行优化设计,设计出圆弧波浪式筛板并进行工业性试验,与传统的直条式筛板的透筛物和分选效果进行比较,比较结果表明圆弧波浪式筛板在筛选效果上比直条式筛板更具有优越性,同时也验证了湍流模型选取的合理性。