NHT 第12篇(续).pptVIP

对 其系数可仿此写出。 （2）边界条件 在本算法中认为所有变量都给出为第一类边界 条件，亦即 值给定的边界条件，因而边界中的修 正值应为0, 无论 还是 其边界值均为零： 12-8-4 双块、多块修正 双块修正：计算区域中有两个物性截然不同的 区域时单块修正可能不适用，需要双块修正 多块修正：与多重网格类似。 12-9 带中心射流的通道内的紊流换热 － 湍流模型及壁面函数法 已知：如图10所示，带一般中心射流的流体进入一 平行板通道。流动为紊流， ＝10－6 , Pr=0.7。 求： 采用 标准 模型及 壁面函数法， 确定 突扩区域中的流场、温度场。 12-9-1 物理问题及其数学描写 T=400 V=100 Tw=0 1.0 4.0 T=100 V=10 图10 控制方程为： 其中 各个变量的扩散系数等如下： 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 边界条件为： （1）入口处：k 取来流平均动能的1％； 按以下公式计算 其中 （2）壁面处按壁面函数法； （3）出口处按局部单向化处理 （4）对称线上，u=0, 其余法向一阶导数为零。 12-9-2 数值方法及程序说明 （1） 源项的处理 采用局部单向化条件时，为防止迭代过程中出口法向流速出现负值对计算区域内的影响，采用提升方法避免之： outflow Vmin （2）出口速度的提升 （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 ******************************************* RELAX(1)=0.8 RELAX(2)=0.8 RELAX(5)=0.6 RELAX(6)=0.6 RELAX(13)=0.6 LSOLVE(1)=.TRUE. LSOLVE(5)=.TRUE. LSOLVE(6)=.TRUE. LPRINT(1)=.TRUE. LPRINT(2)=.TRUE. LPRINT(3)=.TRUE. LPRINT(4)=.TRUE. LPRINT(5)=.TRUE. LPRINT(6)=.TRUE. LPRINT(7)=.TRUE. LPRINT(11)=.TRUE. LAST=50 XL=1. YL=4. L1=7 M1=9 CALL UGRID RETURN ENTRY START DO 100 J=1,M1 DO 100 I=1,L1 U(I,J)=0. V(I,J)=10. V(1,J)=0. V(I,2)=10. IF(I.GT.4) V(I,2)=100. T(I,J)=100. T(1,J)=0. IF(I.GT.4) T(I,1)=400. AKE(I,J)=0.005*V(I,2)**2 ！进口动能的1％ DIS(I,J)=0.1*AKE(I,J)**2 10