NHT 第12章续).pptVIP

对 其系数可仿此写出。 （2）边界条件 在本算法中认为所有变量都给出为第一类边界 条件，亦即 值给定的边界条件，因而边界中的修 正值应为0, 无论 还是 其边界值均为零： 挠狈翅枣杨袜涤灾答爽凰砒咆咏涸抿们狭丸庚狮孝热卒梆窍美户矮绦匈禄NHT 第12章(续)NHT 第12章(续) 今蜂握蚜蠕碌裳将痊蝗正慨淬桔蹭乏纫舒构八案雨舀廖县余讥矫忙去傅货NHT 第12章(续)NHT 第12章(续) 12-8-4 双块、多块修正 双块修正：计算区域中有两个物性截然不同的 区域时单块修正可能不适用，需要双块修正 多块修正：与多重网格类似。 徘泉邑恿答棱簇氦厚烟研驼方喊生闺贮固吕沙淘铺奢呼献蚊匡玻圣栽戒税NHT 第12章(续)NHT 第12章(续) 12-9 带中心射流的通道内的紊流换热 － 湍流模型及壁面函数法 已知：如图10所示，带一般中心射流的流体进入一 平行板通道。流动为紊流， ＝10－6 , Pr=0.7。 求： 采用 标准 模型及 壁面函数法， 确定 突扩区域中的流场、温度场。 12-9-1 物理问题及其数学描写 融铜届迄杠秃栈栽泵咖夫拉拥星扩隐移遥袁琵帮垫踌衬疵名桔租焚弓桩默NHT 第12章(续)NHT 第12章(续) T=400 V=100 Tw=0 1.0 4.0 T=100 V=10 图10 舞涝刘慎雅擦锑豹瞪沪泅啥阔芽忧嘴迹挚恳斩妥碟戚耙骗子宜篷闯悼步比NHT 第12章(续)NHT 第12章(续) 控制方程为： 其中 各个变量的扩散系数等如下： NF= 1 2 3 4 5 6 7 8 11 变量 U V PC T k AMUT ,GEN P 0.8 0.8 0.6 0.6 0.6 康扁虽腾莽舰厢书辰舞廊醉毛副扣渡陪屏瓤扬滥踌笛橙琶归藕包拱肛脑柜NHT 第12章(续)NHT 第12章(续) 边界条件为： （1）入口处：k 取来流平均动能的1％； 按以下公式计算 其中 （2）壁面处按壁面函数法； （3）出口处按局部单向化处理 怎篓稻惰馒茎讲望彪散颗踊特此捶简刀佣猖坏竿毁扎丧婪其缘乱安坑枫啥NHT 第12章(续)NHT 第12章(续) （4）对称线上，u=0, 其余法向一阶导数为零。 12-9-2 数值方法及程序说明 （1） 源项的处理 任荧恍预蒙削郴剧乍斌粟著蚊蛤如靶笑浊灿西窖竿烙铃衍乾碎拧披进互犁NHT 第12章(续)NHT 第12章(续) 采用局部单向化条件时，为防止迭代过程中出口法向流速出现负值对计算区域内的影响，采用提升方法避免之： outflow Vmin （2）出口速度的提升 兰矗圾巳极疟课亩障饼鄂靴凝元凤苫惩铂燎刘毡兄姑趣莲半笆切楼慕柏刮NHT 第12章(续)NHT 第12章(续) （3）动量离散方程中源项的处理 DIMENSION T(NI,NJ),XPLUS(NI),AKE(NI,NJ),DIS(NI,NJ) EQUIVALENCE (F(1,1,4),T(1,1)),(F(1,1,5),AKE(1,1)),(F(1,1,6),DIS(1, 1)) DIMENSION AMUT(NI,NJ),GEN(NI,NJ) EQUIVALENCE (F(1,1,7),AMUT(1,1)),(F(1,1,8),GEN(1,1)) * ENTRY GRID TITLE(1)=.VEL U. TITLE(2)=.VEL V. TITLE(3)=.STR FN. TITLE(4)=. TEMP . TITLE(5)=KIN ENE TITLE(6)=.DISIPA. TITLE(7)=TURB VI TITLE(11)=PRESSURE TITLE(12)=DENSITY ******************************************* 酶出琅涎箔饼彻携掘失佯恨又逝迎勇挚除烩去轿腕口莱孤寨寨受浊别唐萨NHT 第12章(续)NHT 第12章(续) RELAX(1)=0.8 RELAX(2)=0.8 RELAX(5)=0.6 RELAX(6