研一 流体力学课件 传热与流体流动的数值计算实例.pptVIP

1、CFD步骤 第九章 应用举例 第九章，看一些应用的例子。在一些参考文献中看到大量的例子。 要注意的是，仅应用三个通用的计算机程序就可适应广泛的应用。这三个程序是： 1) 二维椭圆型问题； 2) 二维抛物型问题； 3) 三维椭圆型问题； 每一个程序都设计成可用于处理直角坐标与圆柱坐标两种不同的坐标系。在应用中应当辅以有关物理过程的适当的数学模型（如紊流描述、化学反应），并引入有关问题的具体特性（如几何形状、流体物性及边界条件等） 实例1 热传导问题（变量T） 实例1 热传导问题（变量T） 实例1 热传导问题（变量T） 实例2 厚圆筒内的热传导（变量T） 实例2 厚圆筒内的热传导（变量T） 实例1 热传导问题（变量T） 实例1 热传导问题（变量T） §9.2 水平管内的复合对流 沿周向不均匀加热时的水平管内的充分发展层流流动与传热。 §9.3 围绕垂直圆管的熔化 §9.4 带内肋的圆管内的紊流与传热 带纵向直内肋的圆管是一种有效的强化传热的方法。 §9.5 一股吹斜流的紊流射流 由一圆孔喷出的紊流射流可当成一二维抛物流来分析。但若射流被来自于与轴相垂直的一股来流吹偏斜时，流动即表现出有趣的三维椭圆状态；这种流动状态示意图9.10。由烟囱排出的烟羽，垂直起降飞机底下排出的射流及某些气模冷却状态场可算作是这种吹斜流的例子。 § 9.6助推喷射器内的超混射流 图 9.13 超混射流的几何形状和流程 § 9.7 周期性的充分发展通道流 图 9.14 横放的平板组 § 9.8 蒸发器的热工水力分析 * * 传热与流体流动数值计算实例 课程性质：研究生课程 ? 授课时数：40 学时 授课教师：张建文 博士 物理图景的数学描述 （Mathematical Model） 离散化、求离散化系数 ( Discritization Coefficients ) 解方程，获得各场变量分布 ( Distribution Profile ) 求各传递系数 (Transfer Coefficient) 误解：CFD就是把计算问题交给计算机 MODE控制 2、坐标模式 1）Mode=1，x-y 坐标 2）Mode=2，x-r 坐标 3）Mode=3，θ-r 坐标 (1) (2) (3) MODE控制 2、坐标模式 统一模式 (4) r r uθ r θ 3 r 1 uz r z 2 1 1 ux y x 1 R(J) SX(J) ux Y X Mode SX（J）是X（I）的放大尺度，是J的函数；R（J）必须从圆心开始计算，x, y可任取 3、变量 F（I，J，NF） NF：不同的变量，1~10 NF= 1 2 3 4…….10 11 12 13 u v P’ 自给确定 P ρ Γ 4、控制体 XCVS(I)=X(I)—X(I-1) XCV(I) x方向C.V.大小 YCVR(J)=R(J) Δy, y向面积元 4、控制体 极坐标面积元 ΔA=R(J)ΔθΔr 图1 二维平平板热传导问题 Mode=1 ①方程 ②边界条件 上：T=3x+2 下：T=x 左：T=y 右：T=2y+1 各有5个控制体。 k=1，S=0 ③划分网格 L1=M1=7 计算出 XU(I) I=2，L1; YV(J) J=2，M1 均匀网格: CALL UGRID ④START ⑤GAMSOR 各有10个控制体？ k=10? S=10-5T? ⑥课堂讨论 T（I，J）=0.0（内部）实际设定内部（因为边界条件BOUND在START之后）本例在START中给定，BOUND为空，DENSE亦为空（但密度默认为1）。 k=1，S=0 不需更多处理。 仅需赋一次值即可。为节省机时或者也可在START中给出，空出GAMSOR 图2 厚圆筒内的二维热传导问题 Mode=2 坐标：柱坐标 r-x R(1)=1.0 ③边界条件处理 XL处边界： qL=h(Tf-TL)=5(100-TL) R(1)处边界：qR1=50 各有5个控制体。 ②划分网格 L1=M1=7 计算出 XU(I) I=2，L1; YV(J) J=2，M1 均匀网格: CALL UGRID Method=1 直接算 ② 收敛慢! ①TR1=TR2*+qR1/aW= TR2*+50/aW ⑤GAMSOR 各有1