华中科技大学热科学与工程实验室 EP》 HUST Lab of Thermal Science Engineering 1时间与空间的离散化 当进行数值求解时,首先要 W 做的事情是在所研究的时间和 空间区域内把时间和空间分割 成为有限大小的小区域。图92 表示了长柱体矩形截面上区域 K1时刻 离散化的情况。 对于给定的空间区域, K时刻 在x方向上的步长为△x,在 y方向上的步长为4y,用它 K+1时刻 们作为空间尺度可以将矩形 图9.2计算区域的离散化 区域划分成纵横交错的网格 系统,计算区域就被这些网格线分隔成一系列的小的区域,称 为控制面积,对于三维情况则为控制体积或控制容积,因而常 在一般意义上称之为控制体;控制体的中心点称为节点。 20033-2
º]^TU'(i äåÅX7æâãz bz2ç è wxßéé4 /ê³ëìíÆîÛ8 Åï4 YÇNÝ 7 }8²³ 7 }8²³ s´ ð ñqµíÆ ¶2èòóôõö3 ÷øÏùg×ö3Ð2úèd÷ûéü ýþÛ_YǼÅï ýþì_`ýþ_£W 7d¥¦8üÌ ýþì±ýþìÀ ü 4
华中科技大学热科学与工程实验室 EP》 HUST Lab of Thermal Science Engineering 控制体的形状是随着坐标系的不同而改变的,这里的控制 体是一个个的矩形面积。网格的步长在每一个方向上可以均匀 划分,也可以不均匀的划分;所得到网格,相应地被称为均匀 网格或者非均匀网格。 选用不同的步长和不同的划分方法,可以将同一区域划分 出不同大小、不同数目的控制区域,以及不同数目的节点数。 获得每个节点上的温度值,就是导热数值计算的目的。显 然,随着步长的不断减小,节点数目的不断增加,由节点温度 表示的离散的温度场就会更加接近连续的温度场,但计算工作 量也会随之增加。 在时间方向上离散化的步长常用Ar来表示,4r的选取 也是可大可小的,也可以随时间的进程而变化。显然,无限小 的时间步长Δτ亦会使得离散温度变化接近连续的温度改变, 但随之而来的是相应的计算工作量将会增加。 20033-2
ýþìÆX¾÷tu£ghýþ ìXd··íÆÛ_4ö3²³7 d· }8qª ¶2¹qtª ¶2±æÉoö3¸Õùü ª ö3`a+ª ö34 stu²³ztu¶2 %qµud¶2 OtußétuTýþ*tuT T4 ÙÉ · 8UÏXTU4. 5¾²³t é Tt /¯Ï\¯ð ¹¾Ì4 7 }8²³Ws/ Þ ¹Xqßqé¹q¾]!£4.5xé ²³É\ ¾Ì£X¸Õðµ4
华中科技大学热科学与工程实验室 EP》 HUST Lab of Thermal Science Engineering powe 2节点方程的建立 建立节点差分方程可以采用不同的方法,主要分为两大类: 第一类包括泰勒( Taylor)级数展开法和多项式拟合法,它 偏重于从数学的角度进行推导,其优点是便于对离散方程进行 数学特性分析,但缺点是变步长网格的离散方程形式复杂、导 出过程的物理概念不清晰、不能保证差分方程具有守恒特性。 第二类包括控制体热平衡法和控制容积积分法,其优点是 推导过程的物理概念清晰、离散方程系数具有一定物理意义、 保证差分方程具有守恒特性,但缺点是不便于对离散方程进行 数学特性分析。 下面我们采用控制体热平衡法来建立节点方程。 20033-2
#$ 12 !qrstu %i2 ßÊ d T!>%z"#"Q%´ $%ÇnTv&]^Á»' X(ÇY !]^ Tv)Ñ2* X²³ö3 !Æ"+, O)!t-.tØ/012 !1w23)Ñ4 È ýþì©4%zýþ__2%»' X Á)!-. !÷T1wdÝ¥¦ /012 !1w23)Ñ* Xt(ÇY !]^ Tv)Ñ24 ½Ûrsýþì©4%#$ !4
华中科技大学热科学与工程实验室 EP》 HUST Lab of Thermal Science Engineering power ①内节点方程 控制体热平衡法建立节点方程的过程是将能量守恒方程应用 于控制体,建立该节点与周围节点之间的能量平衡关系式,再利 用傅里叶导热定律,最后获得控制体节点温度与周围节点温度之 间的关系式。 考察图9-2中的节点P及其控制体,由能量平衡关系应有 y+g+ds+Φx+①=△E 式中,Φy、①p、Φ、和Φ分别为邻近节点W、E、S和N通过传 导方式传给节点P的热流量;①为单位时间控制体内热源的发热 量;ΔE为控制体单位时间内热能的增加量。由导热傅里叶定 律,在线性温度分布的假设下,时刻K周围节点传给节点P的热 流量分别为: N w=t(w-Tp)Ay. 1 E =*(TE -TP)Ay.1 w x P K 》 -(7-TA)x104y2r1 20033-2
ýþì©4%#$ !)!XµØ23 !Õs Çýþì#$5 67 ÌØ©48÷"fÔ s9h:ÝÓ;6ÙÉýþì 67 Ì 8÷"4 <= *»ýþìØ©48÷Õw "> >!> z>"2° ? zV)B "BN @±># @AýþìB? ± ýþì@AØ49h:Ý Ó7ÐÑ2ÚCD½E67 BN @2° Ê
华中科技大学热科学与工程实验室 EP》 HUST Lab of Thermal Science Engineering 控制体的发热流量y=9Ay1,其中q为内热源强度,即单位 时间单位体积的内热源发热量。 控制体单位时间的内能增加量为 K+1 C 1或△E=p KPt Ax△y·l 前者为时间上的向前差分,而后者为时间上向后差分。以上 关系式中温度T的上标为所在时刻,下标为所在空间位置。 假设4x=4jy,经整理可以得出二维非稳态导热问题的内节点 的两种差分格式的差分方程,即 a)显式差分格式 △ K T+T+T+T)+(1-4 △ pc 定义网格傅里叶数F0、=x2,其物理意义是表征控制体的导 热性能与热储蓄性能之间的对比关系,反映控制体温度随时间 变化的动态特性。显式差分格式简化为 20033-2
ýþì?@ » BFá@A @Aì_B?4 ýþì@AØ ¿a 8}¿12£6a 8}61248 8÷"8 æ7E½ æ7AG4 CD HIqÉOÈ+,- j123"12 !á ݦö39h:T »¥¦XJýþì ÑØ KLÑØÌYM8÷ÍÎýþì¾ A-)Ñ4."123"