第五章 边界层理论单页

1. 引言像交织着爱、恨及相互需要的情侣关系一样，科学和技术之间的历史关系也是那样的令人惊讶和热烈．明显地，我们深深着迷于自亚里士多德、伽利略、牛顿、莱布尼兹一直到爱因斯坦在研究运动时通过思考获得的惊人科学奇迹，同时从车轮、天文望远镜、航天器到计算机的成功技术也是那样地引人入胜．超越世俗的孰优孰劣的疑问，科学和技术难道不是智慧和理性的两张面孔吗？对于引导我们到达宏观世界知识最外缘的数学模型，物理学家能对此满意吗？显然不，人们更需要去认识客观事物、去检验理论、去实验、去模拟和去探索，同时人们也需要去寻求、创造和理解．如今，运动力学（the science of motion-the mechanics-）是建立在数学建模、数值模拟和实验基础上的，后三者相互支撑并确保其平衡．但是，现今实验的成本、建模的难度、数值计算强度的增加都破坏了这个漂亮的结构，使之凌乱不堪．由物理学家构建的数学模型和数学之间的紧密联系要求数学模型的解能够使我们在绝大多数情况下忽略模型的分析而偏重于它的数值解，但有时这是非常困难的．很显然，力学家虽然没有耐心等待数学家对数学模型进行慢慢地分析,但是他们却一定准备沿着一条充满了严谨的启发推理的数学之径前行．自莱布尼兹和十分精确的几何领域中的分析出现后，许多数学工具得到使用．在数学模型的发展和寻求其解中，数学力量对物理学的进步贡献卓越．间或，一些令人惊讶的成果在被物理学家普遍称为“近似理论”中取得．因此，在众多的近似和分析理论中，发散级数被长期使用．数学家对这些级数的情有独钟绝非毫无根据．从明确定义的函数中计算出的级数的各项一定蕴含了扩展函数的信息．一般的，发散级数仅仅是渐进级数．和收敛级数不同之处是渐进级数是部分和，而这种部分和是扩展函数在当某个确定参数很小时的较好的表示．当参数为零时，函数被级数的首项精确表示．当参数非零但很小时，任意部分和是函数的近似．小参数一般用ε来表示．在物理学中,想将一个研究中的数学模型推导成其解是原模型近似解的简化模型，小参数是个决定因素，除了渐进级数的概念外，这还是个渐进展开的概念，也许在更一般的意义来讲，这也是一个近似的概念，而近似正是我们思考的核心．像“理论”这个词的含义有多种理解一样，“近似”一词也有多种解释．即使我们在数理物理学方面自我要求苛刻，但一些歧义依然存在．对照于由欧几里得规范的严格的推理要求的说法，“近似”一词有两种含义．当一种渐进近似获得之后，对于数学家而言，依据数学公式让ε取值足够小，近似的精确度就会非常明确．另一方面对于物理学家来说，近似就是寻找一个参数特定的值但精确度事前却未知．本书的主要目的就是通过提出一种连续的完备的扩张方法使得这两种定义一致．而这种SCEM 方法考虑的就是遵循一些严格的数学程序不得不解决的具体问题．而SCEM 方法讲述的正是贯穿本书的主题，被称之为奇异摄动问题．在这些问题中，当0→ε时，解不会一致趋于让0=ε时的相应的简化方程的解．值得注意的是非一致性是发生在维数小于原始区域的区域中．这也是这类问题被称为边界层的原因．当一个参数很小时，近似解的非一致性是个数学问题．但现在幸运的是，我们作为物理学家能分辨出哪些是已知的物理量，哪些是未知的．这些属于物理问题性质的基础知识能使我们更好地掌握数学模型．这是个具有特征尺度的无量纲化过程的例子， 它能使我们决定一些特定参数是否取得很小．事实上，正是通过被物理描述提供的多种选择去无量纲化才使得奇异摄动受到质疑．在机翼附近的流体只有远离机翼时才是真正意义上的无粘性．然而，对于稳定的不可压流体，作为控制方程的Navier-Stokes 方程，在无量纲的形式中，其物理参数只有雷诺数．现在，如果远离机翼，特征长度尺度就是雷诺数的倒数，和整体相比是非常小的．忽略了雷诺数倒数的项，我们得到的欧拉方程好像粘性消失了．并非流体的粘性系数取了另外的值，而是远离了机翼后速度的梯度足够小，粘性的影响可以被忽略了．这就意味着特征长度尺度的改变将使我们考虑气流壁附近的粘性影响是必不可少的．因此，基于后一种的长度尺度上的雷诺数将不再很大．靠近机翼，Navier-Stokes方程推导成边界层方程，即使这个模型比满足于边界条件的Navier-Stokes模型要简单．如何用在远离机翼时有效的欧拉方程的解和靠近机翼时有效的边界层方程的解去构造一种Navier-Stokes方程的一致有效的近似解呢？这是想解决这类特殊问题的关键之所在．即使研究高雷诺数流体外的其他问题，也是这种主要的思路．如何去寻找一些特征退化的问题和它们的有效区域，如何将它们联系起来，并最终构造一个初始问题的近似，这些都是统领全书的议题的关键点．诚然，应用的主要范围是流体动力学，但是2到6章的部分应用更是广泛且对物理学家来说非常有用，更一般意义上的应用是当我们面对大小参数模型时，我们能想起奇异摄动问题．第2章主要介绍这些问题．甚至最为简单的线性震动的例子都能说明方程的无量纲化过程是能让我们把握数学模型本质的首要关键．在这种框架下，物理学家这种去理解自己研究主题并建构模型的技能显然是解决问题的最强有力的工具．Friedrich's模型问题的简化使精确解很直接，并且对于奇异摄动问题来说，它也是刻画其解主要方法的教学案例．事实上，下一章主要讲两种方法，其中一种是众所周知的匹配渐进展开法，记作MMAE．另外一个方法鲜有人知，就是SCEM，我们也会明白它将是本书剩下的一个中心．第3章处理的是边界层的结构．一般地，物理的考察能给出寻找边界层位置的一些必要的线索．然而，一个简单的问题，一个精确解未知的二阶线性常微分方程，我们能将边界层的位置作为稳定问题来研究．这里给出了几个例子，是寻找一种解的近似和与之要求相应的边界层结构的．在所有的例子中，与初始值问题相比，我们最关注的是至少在局部存在定理无效的边值问题．在讲述正文时，附录能给予补充．在每章的结尾，复杂的问题都能让读者在相应的章节里找到结果．在本书的结尾给出了非常详尽的解．一些问题是真正的研究主题和从未发表的结论导出来的．本书是在法国出版的名为《渐进分析和库什极限》的英文版本．两个版本的大多数章节内容是一样的．其中第9章增补了针对于空气动力学流体的IBL方法的应用，第12章是全新的，用来处理通道流体．这些增补的内容进一步证明了SCEM的有效性．我们非常希望能通过这本书给读者提供一些必要的基本原理（包括数学的、实际的，等等）去理解和应用一些专用于边界层的标准渐进方法．在数理物理学中的许多问题中，这些方法是清楚理解解的结构的基础，而理解解的结构对于得到适当的数值解是非常关键的．此外，对于寻找含有边界层问题解的一致有效近似，我们希望SCEM能给予一个新的解释．在正则形式上，对于这种有效工具，SCEM与MMAE都提供了等效的补充观点．应用这种广义展开的补充，对现在还是不清楚的交互式边界层，SCEM能给我们带来合理的论证．如果推广的SCEM能应用于一些未知领域，那么我们将感到非常欣慰，因为我们工作的目的正是如此．例如，在流体动力学中，不稳定的或者三维的边界层，不稳定性和它们的控制都是将来研究的重要课题．2.奇异摄动问题的介绍物理学中使用的数学模型通常导致一个没有显示解的问题，当小参数出现或当计算区域是很大时，计算它们的数值解也变得相当困难．在这种情形下，往往通过令参数为零或者将研究限制在一个较小的区域，就能简化模型，这就是摄动问题．当小参数趋于零，记为0→ε，可能出现两种情况，一种是原问题的解当0→ε时并不在其定义域内一致地趋向退化问题（Reduced problem ）的解，另一种是在其定义域内一致地趋向退化问题的解．为了解决比较困难的数学问题，奇异摄动问题产生了．为了清楚地描述奇异摄动问题，下面我们来考虑一个积分微分算子（integro-differentialoperator ）εL 并求方程0)],([=Φεεεx L 的一个解),(εεx Φ，这里x 是区域D 中的变量，00εε≤<，0ε是一个固定的充分小的正常数．参数ε是一个无量纲，其蕴含着整个问题都是用无量纲变量来表示的．设0)]([00=Φx L 就是所谓的简化问题，先考虑简单问题，假定范数0Φ-Φε在研究区域D 中很小．用最大模范数，我们有),(0εδεK Max D<Φ-Φ 其中K 是一个不依赖于ε的正常数，)(εδ是一个正函数且.0)(lim 0=→εδε 如果这种性质满足，这个问题就称之为一个正则摄动问题（见问题2-4）．如果一些问题在整个区域D 中不满足上述性质，并且在比区域D 小的一个区域中一个奇点出现，那它就称为一个奇异摄动问题．在本章考虑的模型中，εΦ是已知的．这些教学问题都是用来描述主要的概念性难点以及解决这些难点的各种方法．2.1 正则和奇异问题2.1.1 线性震荡线性振子是正则摄动问题的典型例子．为了进一步地研究，我们来考虑下面的方程0222=++=y dxdy dx y d y L εε (2.1a) 且服从于初始条件.1,000====x x dx dyy (2.1b)函数),(εx y 要求0>x ，并且ε是一个充分小的正参数．所有量都是无量纲．当阻尼很小时，这个方程模拟了质点在阻尼弹性运动系统中的运动．这里“小”的含义在下面的分析中非常重要．下面一个小质量的物理问题是很有趣的．设),,,,(0I k m t y β*是质点离平衡位置的关于时间的位置函数．κ是弹性常数，β是阻尼系数．如果质点是从具有冲量0I 的平衡位置开始运动，那么(2.1a)能记为,022=*++**ky dtdy dt y d m β (2.2a) 且服从于初始条件.,0000I dt dy m y t t ===*=*(2.2b)设y 和x 都是无量纲变量Ly y *=，T t x =， 其中L 和T 分别是未被定义的长度和时间尺度．物体运动源自于冲量，因此令0I m L T =是合理的．有了这些新的变量，(2.2a)就可以写成无量纲形式,0022220=++y dx dy mLk I dx y d k mL I β (2.3a) 且服从于初始条件.1,000====x x dx dyy (2.3b)这样两个无量纲组产生了，并且都包含任意的长度L ，L 可以被m k I L 0=和mk I L 0β=两种方法定义．物理上，两个系数不是同一数量级的，当研究其中一个的小性时，另一个为)1(ο数量级，即可以用渐进分析．下面举两个例子：1.如果是弹性阻尼，第一组k mL I 220比第二组mLk I 0β要大的多，并且 m kI L 0=和k m T =， 因此小参数可以定义为.2m k βε=只要x 有界，下面我们就能看出相应的问题就是典型的正则摄动问题．这还是个小阻尼的例子．方程(2.3a) 化成.0222=++y dxdy dx y d ε根据Poincaré，当0→ε时，解的渐进行为可以按ε的幂展开为.)()()(),(2210 +++=x y x y x y x y εεε对于一个泰勒级数展开，式子中“ ”意味着省略项要比2ε小,并且近似性随着ε越来越小而越来越好．将展开式代入到初始方程并使得相同ε幂的系数相等，下面的方程来自于ε的零次幂系数和一次幂系数．)a 00202=+y dx y d ,1,00000====x x dx dy y)b dx dy y dx y d 012122-=+ .1,00101====x x dx dy y关于0y 的第一个问题是退化问题，能够得出无阻尼解.sin 0x y =关于1y 的第二问题能得到一个修正,sin 1x x y -=于是，一个近似解为.sin )1( +-=x x y ε可以看出，在有限时间区间τ<<x 0（其中τ不依赖于ε）中近似是一致有效的，修正很小．如果时间区间很大，近似无效．通过令1=ετ就能看出．由于时间区间太大，在展开式中出现奇点，所以这类问题被称之为“无穷时间”问题 ．这个专业术语来自于行星运动轨迹的研究．在小的时间尺度下由摄动方法获得的解是有效的，但是“无穷时间”项如果超过以一百年为阶的时间尺度，那它将没有什么实际意义．拿上述的近似解与精确解做比较，我们可以获得启发．由（2.6）给出的近似解正是精确解x e x y x221sin 1),(εεεε--=-泰勒级数展开的第一项．2.第二种情况．质量小，长度和时间尺度是mk I L 0β= 和.k T β=小参数ε是由2βεmk=定义的，(2.3a) 化成022=++y dx dy dx y d ε，.1,000====x x dx dyy （2.7）这个问题是个典型的奇异摄动问题，这奇异摄动问题恰好就是本书的主题．2.1.2 无穷时间问题我们考虑这个方程,0=+=y dxdy y L εε (2.8a) 服从初始条件 10==x y (2.8b)我们求它在0≥x 时的解．用如同（2.5）一样的展开，我们找到一个形如+++++=)()()()(),(2210x y x y x y x y x y n n εεεε的y 的近似值．将这个展开式代入（2.8a ）并使得相同ε幂的系数相等，得出下列相继的方程结果：1.00=dxdy 具有初始条件.100==x y 2.01y dxdy -=具有初始条件.101==x y 3.1--=n n y dx dy 具有初始条件.10==x n y 整理关于n y y y ,,10的这些解，得出.!)1(21),(22+-+++-=n x x x x y n n εεεε （2.9） 从精确解x e x y εε-=),( (2.10)我们可以看出困难之所在了．当x 变得很大时，对于一些项的考察，上面的展开将不再有效．显著的特点就是当ε很小和x 有界时，无穷级数能收敛到精确解，赋予ε一些值，则级数的部分和就是近似解．这时的展开式是一个收敛级数，而其部分和就是渐进展开的最简化形式．当x 大于原点的邻域时，为了转移奇性我们进行坐标变换.11+=x t图2.1 （2.8a ）的近似解由（2.9）给出的；y 的精确解由（2.10）给出的．令),(),(εεx y t Y ≡，我们能将（2.8a ）写成02=-=Y dtdY t Y L εε （2.11a ） 并服从于条件 .11==t Y （2.11b ）一个直接的展开.)()()(),(2210 +++=t Y t Y t Y t Y εεε导出近似解 .)21121()11(1),(22 ++-+-+=tt t t Y εεε （2.12） 相继的近似在原点附近有越来越多的奇性（图2.2）．通过展开下面的精确解可以看出．)]11(exp[),(--=t t Y εε （2.13）这个特点在相似的问题中也存在，但我们使用一个特殊的方法可以处理．下面考虑方程0)(=++=y dxdy y x y L εε .11==x y （2.14） 我们在区间10≤≤x 中求解．图 2.2 （2.11a ）的近似解由（2.12）给出的；Y 的精确解由（2.13）给出的．展开式 )()(),(10x y x y x y εε+=导出下列方程 1.000=+y dxdy x.110==x y 2.dxdy y y dx dy x 0011-=+ .011==x y 结果 +-+=)11(211),(2x x x x y εε （2.15） 能清楚地看出在原点附近，第二个近似比第一个近似的奇性多（图2.3）．精确解12),(22+++-=εεεεx xx y （2.16）在原点是有界的，对于ε大于0的任意值，12),0(+=εεy ，这是一个典型的“无穷时间”问题．图 2.3 （2.14）的近似解由（2.15）给出的；y 的精确解由（2.16）给出的．2.1.3 奇异问题奇异问题的原型是由Friedrichs[36]介绍的，用来证明由Prandtl[78]提出的边界层和粘性流体的一个匹配的．我们考虑下面一个方程 022=-+=a dxdy dx y d y L εε ,10<<a （2.17a ） 服从于边界条件,1,010====x x y y （2.17b ）我们在区间10≤≤x 中求解．这是一个比初始值问题复杂一些的边界值问题，像本章的其他研究问题一样，它的精确解是已知的．0=ε时，它的退化问题是0000=-=a dxdy y L ， 其解由A ax y +=0给出．这里的A 是一个常数，由两个边界条件决定．一般地，同时满足两个边界条件是不可能的．这个特点是某些奇异问题所特有的．当0=ε时，退化问题的阶要比初始问题的阶要低．如果边界条件0=x 是强行赋予的，解则变成.0ax y =由于a y =)1(0，所以这种近似不是一致有效的．类似地，如果强行赋给一个边界条件1=x ，则解变成,10a ax y -+= （2.18）它使得.1)0(0a y -=边界条件在原点不满足一定表明这是一个非一致收敛区域。

1Transport Phenomena, Xu Jian, 2009第五章边界层理论边界层概念 边界层方程 边界层分离2Transport Phenomena, Xu Jian, 20095.1 边界层概念在上述层流动量传递的若干实例的分析中，(1)形状简单；(2)引入了假设：管道无限长、忽略进口段影响。

实际问题要复杂得多。

边界层理论，粘滞力对动量传递影响的一般理论，是粘性流体力学的基础，也与热量传递过程和质量传递过程有着密切的关系。

3Transport Phenomena, Xu Jian, 20095.1 边界层概念Prandtl(1904)提出边界层概念，把统一的流场，划分成两个区域，边界层和外流区；其流体流动（沿流动方向和沿与流动方向垂直的方向）有不同的特点。

边界层：流体速度分布明显受到固体壁面影响的区域。

边界层的形成：¾壁面处流体的“不滑脱”no-slip ¾流体的“内摩擦”作用 边界层厚度δ¾U ＝0Æ0.99 U 04Transport Phenomena, Xu Jian, 20095Transport Phenomena, Xu Jian, 20095.1 边界层概念流过一物体壁面的流体分成两部分¾边界层，粘性流体，不能忽略粘滞力¾外流区，理想流体，可以忽略粘滞力6Transport Phenomena, Xu Jian, 2009边界层理论的要点边界层厚度δ的变化¾前缘处，δ＝0¾x ↑, δ↑；沿壁面的法向将有更多的流体被阻滞¾δ<<x边界层内，δ<<x （距离很小）；0Æ0.99 U 0（速度变化大）¾速度梯度很大，剪切力很大¾流体速度减慢Æ惯性力<<层外，惯性力与粘性力数量级相当7Transport Phenomena, Xu Jian, 2009边界层流动的转变x<x c (临界距离)层流边界层 过渡区 湍流边界层转变判据：¾临界值：5×105；¾特征长度：距前缘的距离；¾特征速度：来流速度0Re xU ρμ=8Transport Phenomena, Xu Jian, 2009圆管进口段效应靠近管壁部分：边界层，速度减慢；厚度不断增大，进口段长度之后，汇交在管中心处；充分发展段的流动状态取决于交汇处边界层的流动状态；进口段的中心部分：无粘性流动区，速度均匀，区域不断缩小，在边界层汇交时消失；沿程速度不断增大Î压降增大（附加压降）；9Transport Phenomena, Xu Jian, 20095.2 边界层方程普兰德边界层方程：量级比较 边界层积分动量方程：动量衡算沿平壁层流边界层的计算：动量积分方程的应用10Transport Phenomena, Xu Jian, 20095.2.1 普兰德边界层方程2222222211x x x x xy y y y y x y u u u u P u u x y x x y u u u u P u u x y y x y μρρμρρ⎛⎞∂∂∂∂∂+=−++⎜⎟∂∂∂∂∂⎝⎠⎛⎞∂∂∂∂∂+=−++⎜⎟⎜⎟∂∂∂∂∂⎝⎠讨论不可压缩流体在平板壁面上的稳态二维层流2222221x x x x Du u u u PDt x x y z υρ⎛⎞∂∂∂∂=−+++⎜⎟∂∂∂∂⎝⎠2222221y y y yDu u u uPDtyx y z υρ⎛⎞∂∂∂∂=−+++⎜⎟⎜⎟∂∂∂∂⎝⎠不可压缩流体的Navier-Stocks 方程不可压缩流体在边界层中作稳态二维流动，方程简化为：y0x u u x y∂∂+=∂∂连续性方程：11Transport Phenomena, Xu Jian, 20095.2.1 普兰德边界层方程普兰德首先发现可以通过比较数量级简化方程：¾Re 较大时，边界层的厚度δ<<x¾边界层内的惯性力和粘性力数量级相当 标准数量级：¾x 为距离的标准数量级，记为x=O(1)¾u 0为速度的标准数量级，记为u 0=O(1)¾边界层厚度δ的数量级记为δ= O(δ)，远远小于O(1) 其他物理量的数量级：¾u x 与u 0是一个数量级，记为u x =O(1)¾y 与u 0是一个数量级，记为u x =O(1)12Transport Phenomena, Xu Jian, 2009其他物理量的数量级(1)(1)(1)x x u u O O x x O ∂Δ≈==∂Δ()222(1)(1)(1)(1)x x u u O O x O O x ∂Δ≈==∂Δyx u u x y ∂∂+=∂∂(1)x u O x∂=∂+(1)y u O y∂=∂()y u O δ=(1)1()()x x u u O O y y O δδ∂Δ≈==∂Δ()22222(1)1()()x x u u O O y O y δδ∂Δ≈==∂Δ22221x x x x xy u u u u P u u x y x x y μρρ⎛⎞∂∂∂∂∂+=−++⎜⎟∂∂∂∂∂⎝⎠数量级(1)(1)×1()()δδ×(1)21()δ13Transport Phenomena, Xu Jian, 2009其他物理量的数量级22221x x x x xy u u u u P u u x y x x y μρρ⎛⎞∂∂∂∂∂+=−++⎜⎟∂∂∂∂∂⎝⎠(1)(1)×1()()δδ×(1)21()δInertial Force=Viscous Force:2()O μδρ=1(1)PO xρ∂≤∂22221y yy yx y u u u u P u u x y y xy μρρ⎛⎞∂∂∂∂∂+=−++⎜⎟⎜⎟∂∂∂∂∂⎝⎠(1)()δ×()(1)δ×()δ1()δ2()δ()δ≤()δ(1)14Transport Phenomena, Xu Jian, 2009普兰德边界层方程2210x x xxy yx u u u dP u u x ydx y u u x yμρρ∂∂∂+=−+∂∂∂∂∂+=∂∂000x y x y u u y u u ====∞=时，时，普兰德边界层方程B.C.通过数量级比较得到的简化方程：应用条件：不可压缩流体在边界层中作稳态二维流动，而且Re 比较大15Transport Phenomena, Xu Jian, 20095.2.2 边界层积分动量方程卡门避开使用N-S 方程，直接对边界层进行衡算x 方向质量衡算：¾左侧进入：¾右侧流出：¾上部外流区进入yxz dxdy 1个单位距离δlyu 0, ρμlx u dy ρ∫()00ll x xu dy u dy dxxρρ∂+∂∫∫()lx u dy dxxρ∂∂∫()()2220000000u (-u )ll l l x x x xlx x u dy u dy dx u dy u dy dxx xdx u u dyx ρρρρρ∂∂+−−∂∂∂=∂∫∫∫∫∫x 向净动量变化率：不可压缩流体沿平板壁面的稳态二维流动16Transport Phenomena, Xu Jian, 2009边界层积分动量方程作用于控制体的x 向外力¾壁面剪切力：¾作用在左右侧面的压力差：1s dx τ−⋅⋅1Pdx l x∂−⋅⋅∂00(u )l x x s Pu u dy l x xρτ∂∂−=+∂∂∫0[,]x y l u u δ∈=00(u )x x sP u u dy x xδρδτ∂∂−=+∂∂∫只考虑x 方向的流动00(u )x x s d dPu u dy dx dxδρδτ−=+∫边界层内外压力近似相等00(u )x x sd u u dy dx δρτ−=∫卡门边界层积分动量方程17Transport Phenomena, Xu Jian, 2009边界层积分动量方程可以求出边界层厚度、流体阻力、曳力系数等；方程有u x ，τw ，δ三个变量，需要补充u x =f 1(y)，τw =f 2(δ)的关系；需要预先假定一个速度分布方程才能求解，故只能算是一种近似的方法。

第五章边界层理论及层流边界层中的传递现象5.1 边界层理论的要点5.1.1 问题的提出前述，Re∝惯性力/粘性力当Re<1时，惯性力<<粘性力，可用“爬流”模型，略去惯性力项，N-S方程==>爬流方程（stokes近似），解决一些实际问题（沉降、润滑、渗流等），获得比较满意的结果。

但工程流动问题，绝大多数的Re很大。

这时，是否可以完全略去粘性力，使Navier-Stokes方程==>Euler方程（理想流体）。

但是，这样的结果与实际情况相差很大。

突出的一例即“达朗倍尔佯谬（paradox）——在流体中作等速运动的物体不受阻力”。

究竟应当怎样才能正确地处理大Re数的流动呢？这个矛盾一直到1904年，德国流体学家普兰德（Prandtl）提出了著名的边界层理论（大Re数的流动中，大部分区域的惯性力>>粘性力，但在紧靠固——流边界的极薄流层中，惯性力≈粘性力），才令人满意地解决了大Re数的流动的阻力问题。

后人把Prandtl 提出的流动边界层概念，推广到流动系统的传热边界层和传质边界层，从而确定传热、传质的速率以及了解有关的影响因素。

还有人研究了边界层中的化学反应，解决了一些实际问题。

因此，边界层理论被认为是近代流体力学的奠基石。

5.1.2 流动边界层（速度边界层）以平板流动为例，x方向一维稳态流动，在垂直壁面的y方向上，流动可划分为性质不同的两个区域：（1）y<δ（边界层）：受壁面影响，法向速度变化急剧，du/dy很大，粘性力大（与惯性同阶），不能忽略。

（2）y>δ(层外主流层)：壁面影响很弱，法向速度基本不变，du/dy≈0。

所以可忽略粘性力（即忽略法向动量传递）。

按理想流体处理，Euler方程适用。

这两个区域在边界层的外缘衔接起来，由于层内的流动趋近于外流是渐进的，不是突变的，因此，通常约定：在流动边界层的外缘处（即y＝δ处），u x＝0.99u∞，δ——流动边界层厚度，δ＝δ(x)。

对于实际流体的流动，无论流动形态是层流还是紊流，真正能求解的问题很少。

这主要是由于流体流动的控制方程本身是非线性的偏微分方程，处理非线性偏微分方程的问题是当今科学界的一大难题，至今还没有找到一套完整的求解方案。

但在实际工程中的大多数问题，是流体在固体容器或管道限定的区域内的流动，这种流动除靠近固体表面的一薄层流体速度变化较大之外，其余的大部分区域内速度的梯度很小。

对于具有这样特点的流动，控制方程可以简化。

首先，由于远离固体壁面的大部分流动区域流体的速度梯度很小，可略去速度的变化，这部分流体之间将不考虑粘性力的存在，视为理想流体，用欧拉方程或伯努利方程就可求解。

而靠近固体壁面的一个薄层——称为流动边界层，在它内部由于速度梯度较大，不能略去粘性力的作用，但可以利用边界层很薄的特点，在边界层内把控制方程简化后再去求解。

这样对整个区域求解的问题就转化为求解主流区内理想流体的流动问题和靠近壁面的边界层内的流动问题。

第一节边界层理论的基本概念一、边界层的定义流体流经固体表面时，靠近表面总会形成一个薄层，在此薄层中紧贴表面的流体流速为零，但在垂直固体表面的方程（法向）上速度增加得很快，即具有很大的速度梯度，甚至对粘度很小的流体，也不能忽略它表现出来的粘性力。

（因此，流体在绕流过固体壁面流动时，紧靠固体壁面形成速度梯度较大的流体薄层称为边界层。

）而在此边界层外，流体的速度梯度很小，甚至对粘度很大的流体，其粘性力的影响也可忽略，流体的流速与绕流固体表面前的流速v0一样。

可以把这部分在边界层外流动的流体运动视为理想流体运动，不考虑粘性力的影响。

边界层内、外区域间没有明显的分界面，而把边界层边缘上的流体流速v x视为v x=0.99v0，因此从固体表面至v x=0.99v0处的垂直距离视为边界层的厚度δ。

二、边界层的形成与特点边界层内的流动可以是层流，也可以是带有层流底层的紊流，还可以是层流、紊流混合的过渡流。

评判边界层层流或紊流的参数为雷诺数Re=vxρ/η，式中v为边界层外边界上流体流速，x为距边界层起点的距离（即流体进入平板的长度）。

第八章 边界层理论§8-1 边界层的基本概念实际流体和理想流体的本质区别就是前者具有粘性。

对层流而言，单位面积摩擦力的大小yu d d μτ=，可以看出，对于确定的流体的等温流场，摩擦力的大小与速度梯度有关，其比例函数即动力粘度。

速度梯度yu d d 大，粘性力也大，此时的流场称为粘性流场。

若速度梯度yu d d 很小，则粘性力可以忽略，称为非粘性流场。

对于非粘性流场，则可按理想流体来处理。

则N-S 方程可由欧拉方程代替，从而使问题大为简化。

Vlv llV v A yu V l t V lt u m ρρμρρ======2223d d d d 粘性力惯性力当空气、蒸汽，水等小粘度的流体与其它物体作高速相对运动时，一般雷诺数很大。

由vVl ==粘性力惯性力Re ，则在这些流动中，惯性力>>粘性力，所以可略去粘性力。

但在紧靠物体壁面存在一流体薄层，粘性力却与惯性力为同一数量级。

所以，在这一薄层中，两者均不能略去。

这一薄层就叫边界层，或叫速度边界层，由普朗特在1904年发现。

a ．流体流过固体壁面，紧贴壁面处速度从零迅速增至主流速度，这一流体薄层，就叫边界层或速度边界层。

b ．整个流场分为两部分 层外，=∂∂yu ，粘性忽略，无旋流动。

层内，粘性流，主要速度降在此，有旋流动。

c ．由边界层外边界上∞=V u %99，来定义δ，δ为边界层厚度。

d ．按流动状态，边界层又分为层流边界层和紊流边界层。

由于在边界层内，流体在物体表面法线方向(即yu ∂∂)速度梯度很大，所以，边界层内的流体具有相当大的旋涡强度；而在层外，由于速度梯度很小。

所以，即使对于粘度很大的流体，粘性力也很小，故可忽略不计，所以可认为，图8-2空气沿平板边界层速度分布外部区域边界层边界层外的流动是无旋的势流。

边界层的基本特征有： (1)1<<Lδ⇒薄层性质，其中L 为物体的长度；沿流方向↑↑→δx 。