山西能源学院教案 授课班级能动1701-1704 授课时间 计2学时 课题(章节 及内容) 5.3边界层型对流传热问题的数学描写 掌握用数量级分析来简化层流边界层对流热微分方程的方法: 教学目的 掌握沿平板受迫流动的层流边界层积分方程组建立及其分析求 和要求 解: 理解流动边界层内的流态。 重点 数量级分析法的基本思想: 难点 边界层微分方程的导出。 教学进程 (含课堂 教学内容:流动边界层和热边界层的含义,边界层能量方程的推 教学内容、 导,流动边界层内的流态,二维稳态边界层型对流传热问题的数 教学方法、 学描述。 辅助手段等)》 教学方法:讲授与练习、启发讨论、诱导式、归纳总结法。 作业布置 复习题1,2 主要 《传热学》第四版,杨世铭,陶文铨, 参考资料 高等教育出版社,2006年8月 课后自我 总结分析
山西能源学院教案 授课班级 能动 1701-1704 授课时间 计 2 学时 课题(章节 及内容) 5.3 边界层型对流传热问题的数学描写 教学目的 和要求 掌握用数量级分析来简化层流边界层对流热微分方程的方法; 掌握沿平板受迫流动的层流边界层积分方程组建立及其分析求 解; 理解流动边界层内的流态。 重 点 难 点 数量级分析法的基本思想; 边界层微分方程的导出。 教学进程 (含课堂 教学内容、 教学方法、 辅助手段等) 教学内容:流动边界层和热边界层的含义,边界层能量方程的推 导,流动边界层内的流态,二维稳态边界层型对流传热问题的数 学描述。 教学方法:讲授与练习、启发讨论、诱导式、归纳总结法。 作业布置 复习题 1,2 主 要 参考资料 《传热学》第四版,杨世铭,陶文铨, 高等教育出版社,2006 年 8 月 课后自我 总结分析
山西能源学院教案 5-3边界层型对流传热问题的数学描写 一、边界层的概念 边界层的概念是1914年普朗特提出的。他认为对流换热热阻的大小主要取决于 靠近壁面附近流体的状况,因为粘滞性起作用的区域仅仅局限在靠近壁面的薄层 内。在此薄层以内,ù、t变化最为剧烈。为计算的简化,因此引入了速度边界 层和温度边界层的概念。 1.速度边界层 主流区 1)定义: 流体流过固体壁面时,由于壁面层流 边界层区 体分子的不滑移特性,在流体黏性力的作 用下,近壁流体流速在垂直于壁面的方向上会从壁面处的零速度逐步变化到来流 速度,如图5一10所示。流体流速变化的剧烈程度,即该方向上的速度梯度,与 流体的黏性力和速度的大小密切相关。普朗特通过观察发现,对于低黏度的流体, 如水和空气等,在以较大的流速流过固体壁面时,在壁面上流体速度发生显著变 化的流体层是非常薄的。因而他把在垂直于壁面的方向上流体流速发生显著变化 的流体薄层定义为速度边界层,而把边界层外流体速度变化比较小的流体流场视 为主流区(势流流动区域)。这样,引入边界层的概念之后,流体流过固体壁面的 流场就人为地分成两个不同的区域,其一是边界层流动区,这里流体的黏性力与 流体的惯性力共同作用,引起流体速度发生显著变化:其二是主流区(势流区), 这里流体黏性力的作用非常微弱,可视为无黏性的理想流体流动,也就是势流流 动。 2)边界层的厚度 我们说边界层是壁面上方流速发生显著变化的薄层,但其边缘所在的位置人 们却是模糊的。在实际分析边界层问题时通常规定,当速度变化达到主流速度的 99%处的空间位置y为速度边界层的外边缘,那么从这一点到壁面的距离就是边 界层的厚度δ 3)引入的意义:
山西能源学院教案 5-3 边界层型对流传热问题的数学描写 一、边界层的概念 边界层的概念是 1914 年普朗特提出的。他认为对流换热热阻的大小主要取决于 靠近壁面附近流体的状况,因为粘滞性起作用的区域仅仅局限在靠近壁面的薄层 内。在此薄层以内, u、t 变化最为剧烈。 为计算的简化,因此引入了速度边界 层和温度边界层的概念。 1. 速度边界层 1)定义: 流体流过固体壁面时,由于壁面层流 体分子的不滑移特性,在流体黏性力的作 用下,近壁流体流速在垂直于壁面的方向上会从壁面处的零速度逐步变化到来流 速度,如图 5-10 所示。流体流速变化的剧烈程度,即该方向上的速度梯度,与 流体的黏性力和速度的大小密切相关。普朗特通过观察发现,对于低黏度的流体, 如水和空气等,在以较大的流速流过固体壁面时,在壁面上流体速度发生显著变 化的流体层是非常薄的。因而他把在垂直于壁面的方向上流体流速发生显著变化 的流体薄层定义为速度边界层,而把边界层外流体速度变化比较小的流体流场视 为主流区(势流流动区域)。这样,引入边界层的概念之后,流体流过固体壁面的 流场就人为地分成两个不同的区域,其一是边界层流动区,这里流体的黏性力与 流体的惯性力共同作用,引起流体速度发生显著变化;其二是主流区(势流区), 这里流体黏性力的作用非常微弱,可视为无黏性的理想流体流动,也就是势流流 动。 2)边界层的厚度 我们说边界层是壁面上方流速发生显著变化的薄层,但其边缘所在的位置人 们却是模糊的。在实际分析边界层问题时通常规定,当速度变化达到主流速度的 99%处的空间位置 y 为速度边界层的外边缘,那么从这一点到壁面的距离就是边 界层的厚度 。 3)引入的意义:
a)缩小计算区域。以温度为20℃的空气沿平板的流动为例,在不同来流速 度”下,相对于平板长度1,δ是一个比1小一个数量级以上的小量。在这样小 的薄层内,流体的速度要从0s变化到接近于主流流速,所以流体在垂直于主 流方向上的速度变化是十分剧烈的。由于边界层内都很大,动量微分方程 oy'oy 式中的粘性力和惯性力及能量微分方程中的导热和对流项都需考虑,有了边界 层,可以将对流换热的研究集中在边界层内。 b)边界层内的流动和换热可以利用边界层的特点进行简化,即在边界层中 可以运用数量级的分析方法对纳维一斯托克斯方程进行简化。 例如:运用数量级分析的方法,边界层内粘性流体的稳态动量方程可简化为: +=+ p dx ay2 与二维稳态的纳维一斯托克斯方程相比,上述运动微分方程的特点是: (①)在u方程中略去了主流方向的二阶导数项:()略去了关于速度9的动量方程。 (认为边界层中乎=0。 y 4)临界雷诺数: 流体的流动可区别为层流和湍流两种。 边界层也呈现出层流和湍流。 流体以“,的流速沿平板流动。在 2/ 平扳的起始段,δ很薄。随着x的 层流 过渡 增加,由于壁面粘滞力的影响逐渐 向内部传递,δ逐渐增厚,同时黏 性力对流场的控制作用也逐步减弱,从而使边界层内的流动变得紊乱。在某 一距离x。以前一直保持层流特征。此时流体作有秩序的分层运动,各层之间 互不干扰。这时的边界层称层流边界层。随着边界层厚度的增加,本来很平顺的 层流流动状态就会变成紊乱无序的紊流流动状态。此时流体质点在沿x方向流动 的同时,又作着紊乱的不规则脉动,紊流流动状态的边界层则为紊流边界层,我 们把边界层从层流过渡到紊流的x值称为临界值,记为x©,其所对应的雷诺数称 为临界雷诺数,即Rec
a)缩小计算区域。以温度为 20℃的空气沿平板的流动为例,在不同来流速 度 u 下,相对于平板长度 l, 是一个比 l 小一个数量级以上的小量。在这样小 的薄层内,流体的速度要从 0 m/s 变化到接近于主流流速,所以流体在垂直于主 流方向上的速度变化是十分剧烈的。由于边界层内 y t yu , 都很大,动量微分方程 式中的粘性力和惯性力及能量微分方程中的导热和对流项都需考虑,有了边界 层,可以将对流换热的研究集中在边界层内。 b)边界层内的流动和换热可以利用边界层的特点进行简化,即在边界层中 可以运用数量级的分析方法对纳维-斯托克斯方程进行简化。 例如:运用数量级分析的方法,边界层内粘性流体的稳态动量方程可简化为: 与二维稳态的纳维—斯托克斯方程相比,上述运动微分方程的特点是: (i)在 u 方程中略去了主流方向的二阶导数项;(ii)略去了关于速度 的动量方程。 (iii)认为边界层中 0 y p 。 4)临界雷诺数: 流体的流动可区别为层流和湍流两种。 边界层也呈现出层流和湍流。 流体以 u 的流速沿平板流动。在 平扳的起始段, 很薄。随着 x 的 增加,由于壁面粘滞力的影响逐渐 向内部传递, 逐渐增厚,同时黏 性力对流 场的控制作用也逐步减弱,从而使边界层内 的流动变得紊乱。在某 一距离 xc 以前一直保 持层流特征。此时流体作有秩序的分层运动, 各层之间 互不干扰。这时的边界层称层流边界层。随着边界层厚度的增加,本来很平顺的 层流流动状态就会变成紊乱无序的紊流流动状态。此时流体质点在沿 x 方向流动 的同时,又作着紊乱的不规则脉动,紊流流动状态的边界层则为紊流边界层,我 们把边界层从层流过渡到紊流的 x 值称为临界值,记为 xc,其所对应的雷诺数称 为临界雷诺数,即 Rec
惯性力_PuaXe 粘性力 =4 实验研究的数据表明,流体平行流过平板的临界雷诺数大约是(一般对平板: R.=500000)。当然,这一数据与来流速度的紊流程度和平板前沿的几何形状 密切相关,因而临界雷诺数的数值会在一定范围内变动。 我们不难发现,要使边界层的厚度远小于流动方向上的尺度(即),也就是 所说的边界层是一个薄层,这就要求雷诺数必须足够的大。因此,对于流体流过 平板,满足边界层假设的条件就是雷诺数足够大。由此也就知道,当速度很小、 黏性很大时或在平板的前沿,边界层是难以满足薄层性条件。 4)边界层理论的四个基本要点 a.当粘性流体沿固体表面流动时,流场划为主流区(势流区)和边界层区。 在边界层区内,速度在垂直于壁面方向剧烈变化。而主流区速度低度几乎为零。 速度边界层成立的条件是Re>1。 b.主流区的流动视为理想流体的流动,用描述理想流体的方程求解。边界层 区应考虑粘性的影响,用粘性流体的边界层微分方程求解,其特点是主流方向流 速的二阶导数项忽略不计。 ℃.边界层的流动状态分为层流和紊流。紊流边界层分为层流底层、缓冲层和 紊流核心区。 d.边界层的厚度是与壁面尺寸相比很小的量,且远不只小于一个数量级。即: 6<1. 2.热(温度)边界层 1)定义 当流体流过平板而平板的温度tw与来流流体的温度to不相等时(即存在温度 差),对于上述的低黏性流体,如果流体的热扩散系数也很小,在壁面上方也能 形成温度发生显著变化的薄层,常称为热边界层。 2)热边界层厚度 仿照速度边界层的约定规则,以过余温度为来流过余温度的99%处定义为热
实验研究的数据表明,流体平行流过平板的临界雷诺数大约是(一般对平板: Rec=500000)。当然,这一数据与来流速度的紊流程度和平板前沿的几何形状 密切相关,因而临界雷诺数的数值会在一定范围内变动。 我们不难发现,要使边界层的厚度远小于流动方向上的尺度(即 l),也就是 所说的边界层是一个薄层,这就要求雷诺数必须足够的大。因此,对于流体流过 平板,满足边界层假设的条件就是雷诺数足够大。由此也就知道,当速度很小、 黏性很大时或在平板的前沿,边界层是难以满足薄层性条件。 4)边界层理论的四个基本要点 a. 当粘性流体沿固体表面流动时,流场划为主流区(势流区)和边界层区。 在边界层区内,速度在垂直于壁面方向剧烈变化。而主流区速度低度几乎为零。 速度边界层成立的条件是 Re>>1。 b. 主流区的流动视为理想流体的流动,用描述理想流体的方程求解。边界层 区应考虑粘性的影响,用粘性流体的边界层微分方程求解,其特点是主流方向流 速的二阶导数项忽略不计。 c. 边界层的流动状态分为层流和紊流。紊流边界层分为层流底层、缓冲层和 紊流核心区。 d. 边界层的厚度是与壁面尺寸相比很小的量,且远不只小于一个数量级。即: l 。 2. 热(温度)边界层 1)定义 当流体流过平板而平板的温度 tw 与来流流体的温度 t∞不相等时(即存在温度 差),对于上述的低黏性流体,如果流体的热扩散系数也很小,在壁面上方也能 形成温度发生显著变化的薄层,常称为热边界层。 2)热边界层厚度 仿照速度边界层的约定规则,以过余温度为来流过余温度的 99%处定义为热 c c c u x u x Re 粘性力 惯性力
边界层的外边缘。(当壁面与流体之间的温差达到壁面与来流流体之间的温差的 0.99倍时,即,此位置就是边界层的外边缘),而该点到壁面之间的距离则是 热边界层的厚度,记为t。如果整个平板都保持温度tw,那么,x=0时δtx)=O, 且随着x值的增大逐步增厚。在同一位置上热边界层厚度与速度边界层厚度的相 对大小与流体的普朗特数P有关,也就是与流体的热扩散特性和动量扩散特性 的相对大小有关。理论上B=型 St~fRex,Pr) 由此式可以看出,热边界层 是否满足薄层性的条件,除了 Free stream -6,x Rex足够大之外还取决于普朗特 Thermal boundary 数的大小,当普朗特数非常小 layer 时,热边界层相对于速度边界层 Tw 就很厚,反之则很薄。 二、边界层微分方程组 利用上述的边界层的概念,我们可以对流体流过平板的对流换热微分方程组进 行相应的简化。按照普朗特的边界层为一个薄层的假设,以及满足这一假设下的 流场特征,我们运用数量级分析的方法可以对上述微分方程组在边界层中得以简 化。 数量级分析:比较方程中各量或各项的量级的相对大小:保留量级较大的量 或项:舍去那些量级小的项,方程大大简化 由边界层假设,如果设定X的数量级为1,(x与1相当),那么Y的数量级 t0(1): 定义为△(一个小量);在设定主流方向上的速度u的数量级为1,温度: 的情况下,对动量微分方程进行数量级分析。 对于稳态过程而言,如不考虑体积力,则动量方程为: 8'u8'u dy pdx 会 o2v 8'v ++2
边界层的外边缘。(当壁面与流体之间的温差达到壁面与来流流体之间的温差的 0.99 倍时,即 ,此位置就是边界层的外边缘),而该点到壁面之间的距离则是 热边界层的厚度,记为δt 。如果整个平板都保持温度 tw,那么,x=0 时δt(x)=0, 且随着 x 值的增大逐步增厚。在同一位置上热边界层厚度与速度边界层厚度的相 对大小与流体的普朗特数 Pr 有关,也就是与流体的热扩散特性和动量扩散特性 的相对大小有关。理论上 p r c a P δt~f(Re×, Pr) 由此式可以看出,热边界层 是否满足薄层性的条件,除了 Re×足够大之外还取决于普朗特 数的大小,当普朗特数非常小 时,热边界层相对于速度边界层 就很厚,反之则很薄。 二、边界层微分方程组 利用上述的边界层的概念,我们可以对流体流过平板的对流换热微分方程组进 行相应的简化。按照普朗特的边界层为一个薄层的假设,以及满足这一假设下的 流场特征,我们运用数量级分析的方法可以对上述微分方程组在边界层中得以简 化。 数量级分析:比较方程中各量或各项的量级的相对大小;保留量级较大的量 或项;舍去那些量级小的项,方程大大简化 由边界层假设,如果设定 X 的数量级为 1,(x 与 l 相当),那么 Y 的数量级 定义为Δ(一个小量);在设定主流方向上的速度 u 的数量级为 1, 温度: 的情况下,对动量微分方程进行数量级分析。 对于稳态过程而言,如不考虑体积力,则动量方程为: ( ) 1 2 2 2 2 y u x u x p y u v x u u ( ) 1 2 2 2 2 y v x v y p y v v x v u t ~ 0(1); Tw
分析:因为6<1,则11,6~△;来流速度u~1 (1)比较u与a: 若架的数量级为}即1时,则多的数量级为公相应地整~尔,而染~1 ax 1 y 2 (2)分析流体垂直于板面方向的速度ⅴ相当于什么数量级? y u沿边界层厚度由0到w。: u~40~0(1) 由连续性方程: dy ouu0(1) ay Ox 即u必与m属于同一数量级,~1可以得出v的数量级为△。在边界层这 Bx 样一个极薄的薄层中,流体的垂直速度从平板表面处的零,增加到边界层外缘处 的v,而其变化率的数量级却仍与相同,这就表明速度v的数值一定大大 地低于u的数值,即v~△。这样v在x方向上的变化率也必远远小于u在x方 向上的变化率,既然u~1,则~△。同理,既然u~1,~△。同(1)分 0x2 析,,也必更加大大地大于y,既然v~△,则~1。 oy? 024 综上所述,动量方程的数量级可分别表示成: auau 1 op 82u 82u u- -+v +y( (a) O dy p ax dr? 11△1 1 1 △ A -+ +XG (b) poy 2 1△△△ 1 上两式右边括号内的两项相比,后一项要比前一项大得多,所以u和均 @x2 0x2 可略去。在边界层内,粘性力与其它力属于同一数量级,从上两式左边可以看出, 惯性力分别相当于1和△数量级,按照方程两边数量级应该一致的原则,因此粘
分析:因为 l ,则l ~1, ~;来流速度 u ~1 (1)比较 xu 与 yu : 若 xu 的数量级为 1 1 即 1 时,则 yu 的数量级为 1 ,相应地 2 2yu ~ 21 ,而 2 2xu ~1 (2)分析流体垂直于板面方向的速度 v 相当于什么数量级? u 沿边界层厚度由 0 到 u: 由连续性方程: u ~ u ~ 0(1) ~ ~ 0(1) l u x u y v 即 xu 必与 y v 属于同一数量级, y v ~1 可以得出 v 的数量级为Δ。在边界层这 样一个极薄的薄层中,流体的垂直速度从平板表面处的零,增加到边界层外缘处 的 v ,而其变化率 y v 的数量级却仍与 xu 相同,这就表明速度 v 的数值一定大大 地低于 u 的数值,即 v~Δ。这样 v 在 x 方向上的变化率也必远远小于 u 在 x 方 向上的变化率,既然 xu ~1,则 x v ~Δ。同理,既然 2 2xu ~1, 2 2x v ~Δ。同(1)分 析, 2 2y v 也必更加大大地大于 2 2x v ,既然 2 2x v ~Δ,则 2 2y v ~ 1 。 综上所述,动量方程的数量级可分别表示成: ( ) 1 2 2 2 2 y u x u x p y u v x u u (a) 1 1 1 Δ 1 2 1 1 2 1 ( ) 1 2 2 2 2 y v x v y p y v v x v u (b) 1 1 Δ 2 1 2 上两式右边括号内的两项相比,后一项要比前一项大得多,所以 2 2xu 和 2 2x v 均 可略去。在边界层内,粘性力与其它力属于同一数量级,从上两式左边可以看出, 惯性力分别相当于 1 和Δ数量级,按照方程两边数量级应该一致的原则,因此粘
性力也必相应地相当于1和△数量级,即运动粘度y~△?。 (3)考察压力梯度的数量级: 在a式中惯性力和粘性力都相当于1的数量级,假定p的数量级是1,那么 迎的数量级至多也是1,否则方程式无法成立。同理,即的数量级至多也是△, Bx (注意:乎与不同,紧贴固体表面的流体流速u等于零,但此处的压力并不 等于零。边界层断面上的压力变化并非从零变到边界层外的压力P,因而~上, y△ 但即的数量级却至多只是△。)这一结果告我们在边界层中压力不随Y的变化而 dy 变化,仅仅是X的函数。 (4)比较a式和b式的数量级后,发现b式的各项要比式a的各项小一个数量 级,因而b时可以舍去,于是边界层的动量微分方程就由两个变为一个,即 u0+0-12+y9 dy p ax dy2 而边界层内的压力梯度仅沿x方向变化,而边界层内法向的压力梯度极小。即 边界层内任一截面压力与y无关而等于主流压力,则乎,虫 dx dx 当流体纵掠平板时,边界层外主流中的流速没有变化,按照铂努力方程,主 流中的压力也不变,史虫0 可以得出: Ouou 8'u -+v9 =Y Ox Oy 02 同样对边界层中二维稳态能量方程进行数量级比较,有 821821 器+器=品(假》品(月 数量级 1片6日 (得加(合a 1 1 2 上述结果表明:(1)要使等号前后的项有相同的数量级,热扩散率α必须具 有82的数量级。实际上,除液态金属外的流体都满足这一分析。(2)等号 后方括号内的两个中,驴心爱,圆面可以能主流方向的二阶号数项器 略去。于是得到二维,稳态、无内热源的边界层能量方程为
性力也必相应地相当于 1 和Δ数量级,即运动粘度 ~Δ 2。 (3)考察压力梯度的数量级: 在 a 式中惯性力和粘性力都相当于 1 的数量级,假定 的数量级是 1,那么 xp 的数量级至多也是 1,否则方程式无法成立。同理, y p 的数量级至多也是Δ, (注意: y p 与 yu 不同,紧贴固体表面的流体流速 u 等于零,但此处的压力并不 等于零。边界层断面上的压力变化并非从零变到边界层外的压力 p,因而 yu ~ 1 , 但 y p 的数量级却至多只是Δ。)这一结果告我们在边界层中压力不随 Y 的变化而 变化,仅仅是 X 的函数。 (4)比较 a 式和 b 式的数量级后,发现 b 式的各项要比式 a 的各项小一个数量 级,因而 b 时可以舍去,于是边界层的动量微分方程就由两个变为一个,即 2 2 1 yu xp yu v xu u 而边界层内的压力梯度仅沿 x 方向变化,而边界层内法向的压力梯度极小。即 边界层内任一截面压力与 y 无关而等于主流压力,则 dxdp xp 2 2 1 y u dx dp y u v x u u 当流体纵掠平板时,边界层外主流中的流速没有变化,按照铂努力方程,主 流中的压力也不变, 0 dxdp xp 可以得出: 2 2yu yu v xu u 同样对边界层中二维稳态能量方程进行数量级比较,有 ( ) 2 2 2 2 y t x t a y t v x t u
8't u-+y-=a- dx dy dy2 至此,二维、稳态、无内热源的边界层换热微分方程组为: au+-0 dx ay 0+-1电+v p dx 3 0t01 021 -+v 微分方程组经过在边界层中简化后,由于动量方程和能量方程分别略去了主 流方向上的动量扩散项和热量扩散项,从而构成上游影响下游而下游不影响上 游的物理特征。这就使得动量方程和能量方程变成了抛物型的非线性微分方程: 且由于动量方程由两个变成为一个,而且项可在边界层的外边缘上利用伯努利 方程变成的形式。于是方程组在给定的边值条件下可以进行分析求解。布劳修 斯在引入相似变量下把动量方程变成常微分方程后由进行了求解,而普尔豪森也 对能量方程进行了相同的处理。获得相应的速度分布和温度分布,进而求得壁面 的摩擦系数和对流换热系数。 普朗特从他多年从事水力学试验中所观察到的事实出发,创造性的用数量级对 比法简化了原始的微分方程组,开拓了对流换热理论界的道路,成为流体力学和 传热学发展史的里程碑。在他发表边界层微分方程组时,许多著名数学家都不以 为然,认为他的这种舍弃微分方程中若干项以至整个方程的做法是荒谬的。(1904 年发表)直到1908年,他的学生布劳修斯用边界层方程获得了外掠平板的解, 并且其正确性被尔后的实验证实,普朗特的边界层理论才站住脚。所以实验研究 仍是技术发展必不可少的。 三、速度边界层厚度与热边界层厚的关系 1.比较边界层的动量方程和能量方程
2 2y t a y t v x t u 至此,二维、稳态、无内热源的边界层换热微分方程组为: 0 y v xu 2 2 1 y u dx dp y u v x u u 2 2 y t a y t v x t u 微分方程组经过在边界层中简化后,由于动量方程和能量方程分别略去了主 流方向上的动量扩散项 和热量扩散项 ,从而构成上游影响下游而下游不影响上 游的物理特征。这就使得动量方程和能量方程变成了抛物型的非线性微分方程; 且由于动量方程由两个变成为一个,而且 项可在边界层的外边缘上利用伯努利 方程变成 的形式。于是方程组在给定的边值条件下可以进行分析求解。布劳修 斯在引入相似变量下把动量方程变成常微分方程后由进行了求解,而普尔豪森也 对能量方程进行了相同的处理。获得相应的速度分布和温度分布,进而求得壁面 的摩擦系数和对流换热系数。 普朗特从他多年从事水力学试验中所观察到的事实出发,创造性的用数量级对 比法简化了原始的微分方程组,开拓了对流换热理论界的道路,成为流体力学和 传热学发展史的里程碑。在他发表边界层微分方程组时,许多著名数学家都不以 为然,认为他的这种舍弃微分方程中若干项以至整个方程的做法是荒谬的。(1904 年发表)直到 1908 年,他的学生布劳修斯用边界层方程获得了外掠平板的解, 并且其正确性被尔后的实验证实,普朗特的边界层理论才站住脚。所以实验研究 仍是技术发展必不可少的。 三、速度边界层厚度与热边界层厚的关系 1. 比较边界层的动量方程和能量方程
现在,我们来分析热边层厚度6、与流动边界层厚度6之间的关系。为 此,考虑一个掠过平板的强制对流换热问题。在这类强制对流中,重力场可 忽略不计,且压力梯度为零,于是式(5-15)简化为 u ax au-v3y2 3u+v7y (5-18) 将此式与边界层能量微分方程式(5-16)相比较,发现它们在形式上是完全 类似的。只要y=a,且4与t具有相同的边界条件,例如y=0时t=tw, =uw(uw=0并不影响讨论),及y一∞时t*tm,*u四,则式(5-16)与式 (518)有相同形式的无量纲形式的解,即《二“”与,二的分布完全相 4oo-uw too-tw 同。换句话说,当a=1时,如果热边界层的厚度的定义与流动边界层厚 度的定义相同(例如均取来流值的99%的位置作为边界层的外边界),则有 8,=8。可见比值la可以表征热边界层与流动边界层的相对厚度。我们 称v/a=cplA为背朗特数,记为P,,它反映了流体中动量扩散与热扩散能 力的对比。除液态金属的P,数为0.01的数量级外,常用流体的P数在 0.6-4000之间,例如各种气体的P一数大致在0.6一0.7之间。流体的运 动粘性反映了流体中由于分子运动而扩散动量的能力。 在忽略动量方程压力项后,当P=1时,动量方程与能量方程完全相同。即速度 分布的解与温度分布完全相同,此时速度边界层厚度等于温度边界层厚度。 当P>1时,Pr=u/a,o>,粘性扩散>热量扩散,速度边界层厚度>温度边界层 厚度。 当Pr1 层流:温度呈抛物线分布 湍流:温度呈幂函数分布 湍流边界层贴壁处的温度梯度明显大于层流
层流:温度呈抛物线分布 湍流:温度呈幂函数分布 湍流边界层贴壁处的温度梯度明显大于层流 动粘性反映了流体中由于分子运动而扩散动量的能力。 在忽略动量方程压力项后,当 Pr=1 时,动量方程与能量方程完全相同。即速度 分布的解与温度分布完全相同,此时速度边界层厚度等于温度边界层厚度。 当 Pr>1 时,Pr=υ/a,υ>a,粘性扩散 >热量扩散,速度边界层厚度>温度边界层 厚度。 当 Pr<1 时,Pr=υ/a,υ<a,粘性扩散 <热量扩散,速度边界层厚度<温度边 界层厚度。 这也可以从公式得出
0 0 层流 亲流 ().>(). 故:湍流换热比层流换热强! δ与的关系:分别反映流体分子和流体微团的动量和热量扩散的深度 2.求解示例一流体外掠等温平壁的对流换热: 例如:对于主流场均速”。、均温:,,并给定恒定壁温的情况下的流体纵掠平 板换热,即界条僬为=0,1=1 y=6u=4o, t=t 求解上述方程组(层流边界层对流换热微分方程组),可得局部表面传热系数h,的 表达式 注意:层流 =0=8 笑-03=8) Nu,=0.332 Re2.Pr1/3 3.特征方程及努塞尔特数 式中: hx 努塞尔Nusselt)数 λ Re ux 雷诺Reynolds)数 注意:特征尺度 V 普朗特数
y u u t t y u v t tw , 0 0, 0, 3 1 2 1 0.332 a u x x hx 注意:层流 3 1 2 1 0.332 a h x u x x 1 2 1 3 0.332 Re Pr Nu x x 与 t 的关系:分别反映流体分子和流体微团的动量 和热量扩散的深度 2.求解示例——流体外掠等温平壁的对流换热: 例如:对于主流场均速u、均温 t ,并给定恒定壁温的情况下的流体纵掠平 板换热,即边界条件为 求解上述方程组(层流边界层对流换热微分方程组),可得局部表面传热系数 x h 的 表达式 3.特征方程及努塞尔特数 w t w L 故:湍流换热比层流换热强! y T y T , , 式中: hx Nu x x 努塞尔(Nusselt)数 u x x Re 雷诺(Reynolds)数 a Pr 普朗特数