流体力学——平板边界层编程
《流体力学》课件 4.3 普朗特边界层方程
2u 0
y 02
?
?
1
St
0
t 0
1
u0
0
x 0
0
0
y 0
假设：1. 设在我们所研究的问题中，
Eu
?
p y
0 0
1 Re
2
x 0
0 2
1
2 0
y 02
当地导数与局部导数相当（或更小）；
2. 压力梯度力作为被动力，与方程中的
惯性力或粘性力中的大者相当；
u 0 x0
0
u 0 y 0
Eu
p0 x0
1 2u0 Re x02
2u 0 y 0 2
St
0
t 0
u0
0
x0
0
0
y 0
Eu
p0 y 0
1 2 0
Re x02
2 0
y 0 2
2. 量阶的概念及一般运算法则
指某个物理量在整个区域内相对于标准小参数而言的平均水平。
取 为估阶的标准。
**U 2 uU ud y
0
uU ud y
0
U
2
0
d
y
0
1
u U
d
y
0
u 2
d
y
* , 1 u d y
0 U
动量通量损失：U 2 * u2 d y
0
理想流体通过流管ⅠⅡ动量通量为— U 2 *
粘性流体通过流管ⅠⅡ 动量通量为—
u 2
3. 在边界层内惯性力与粘性力同阶。
1
1
u 0
x
0
0
y 0
＝0
1
St
平板边界层实验指导1213
平板边界层实验指导一．实验目的1）测量平板边界层流速剖面，加深对边界层概念的认识；了解层流和湍流边界层的差异。
2）掌握热线风速仪和皮托管测速技术。
二．实验原理U 大Re 数绕平板流动，在平板边界附近有一个薄层，流速从平板处的零值，经过该层迅速增大到接近来流速度U ，此薄层被称为边界层。
通常定义0.99u =处到平板的距离为边界层厚度。
在平板的前段，边界层内流动呈层流状态，即层流边界层。
建立直角坐标系如图1，原点在平板前端，x 轴沿来流方向，轴垂直平板。
定义局地雷诺数y Re x Uxν=，ν为流体的运动学粘性系数。
从平板前端向后，在某个x 位置以后，Re x 足够大，边界层内流动变得不稳定；继续向后，当Re x 超过临界值Re xc 后，边界层内流动发展为湍流。
Re xc 被称为转捩雷诺数，其大小受多种因素影响，包括来流湍流度、平板粗糙度和其他扰动等。
对光滑平板边界层的观测研究表明，在低湍流度风洞中（湍流度低于1%），Re xc 可达；对于较大的来流湍流度，Re 610xc 也可以低至几千甚至几百。
在层流边界层中，粘性力与惯性力同量级。
除平板前端外（Re 100x <），层流边界层流速剖面满足Blasius 解，即()u Uf η′=，f满足200,0,1f ff f f f ηη′′′′′+=⎧⎪′===⎨⎪′=∞=⎩--------------------(1)其中η=该速度剖面如图2所示。
相应地，层流边界层厚度c δ≈从固壁向外，湍流边界层可分为粘性底层、过渡区和湍流核心区。
在粘性底层内，分子粘性应力远大于湍应力，流速呈线性分布。
在湍流核心区，情况正好相反，分子粘性可略，流速呈对数分布。
设u u u +∗=，yu y ν∗+=，其中u为脉动平均流速，u ∗=为摩擦风速，wτ为壁面上的切应力，ρ为流体密度。
在粘性底层u y +=+，-------（2-1） 在湍流核心区1ln u y κ++=C +，-------（2-2）常数和由实验确定。
平板层流边界层内速度分布规律
在流体力学中，平板层流边界层是一个非常重要的概念，它描述了流体在平板表面附近的流动情况。
在平板层流边界层内，流体的速度分布呈现一种特定的规律，这种规律可以用数学公式来描述。
根据实验和理论分析，我们发现平板层流边界层内的速度分布呈现线性分布的特点。
也就是说，在边界层内，流体的速度随着离开平板表面的距离的增加而线性增加。
这种线性分布规律可以用公式表示为：u = u0 + βx，其中u是x位置处的速度，u0是平板表面处的速度，β是速度梯度，x是距离平板表面的距离。
这个公式非常简单，但它却准确地描述了平板层流边界层内速度分布的基本规律。
这个规律是通过大量的实验和理论分析得出的，具有很高的可信度。
通过这个公式，我们可以了解到流体的速度是如何随着离开平板表面的距离而变化的，这对于理解流体动力学的基本规律和解决实际工程问题具有重要的意义。
此外，平板层流边界层的形成还受到多种因素的影响，如流体本身的性质、平板表面的粗糙度以及流体的流动条件等。
不同的流体和流动条件下，平板层流边界层的形成机制和速度分布规律可能会有所不同。
因此，在实际应用中，我们需要根据具体情况对平板层流边界层的速度分布进行测量和计算，以便更好地理解和控制流体流动。
总之，平板层流边界层内速度分布的线性规律是一个非常重要的流体动力学概念，它对于理解流体动力学的基本规律和解决实际工程问题具有重要的意义。
通过深入研究和探索这个规律，我们可以更好地掌握流体动力学的本质，为未来的科学研究和技术创新提供更加
坚实的基础。
平板层流边界层编程说明
g z n f 0 1 3
2. 求解关于无量纲流函数f的常微分方程组，得到无量纲速度；
1 df1 f f f 0 f2 2 d f1 0 0 转换为常微分方程组 f 0 0 df 2 f3 f 2 0 0 f 0 0 d 1 3 f 0 3 df 3 1 f f 1 f 0 1 3 3 d 2 u ? 四阶 Runge- Kutta法求解 U v ? U
五、程序框图
设定η合理的分析范围和步长
采用四阶Runge-Kutta法求解置换函数 g，得到待定系数β
根据β，采用四阶Runge-Kutta法求解 无量纲流函数f
确定无量纲速度u、v
画出u、v的分布曲线
六、参考示例
u U f v 1 f f U 2 U x
四、程序的原理说明
设定η合理的分析范围和步长，采用改进的四阶Runge-Kutta 法进行求解： 1. 求解关于置换函数g的常微分方程组，确定待定系数β；
dg1 1 dz g 2 g 2 gg 0 g1 0 0 转换为常微分方程组 g 0 0 dg 2 g3 g 2 0 0 g 0 0 dz g 3 0 1 dg 1 3 g1 g 3 g 0 1 2 dz
平板层流用四阶Runge-Kutta法分析二维平板层流边界层流动。 求解：流场内无量纲流函数及其各阶导数，画出速度分布曲 线。
二、问题的描述
U∞
η
f 1
f 0 0,
《水力学》课件——第九章 边界层理论基础
位移厚度 1
因为有了边界层，使通
y
过断面的流量比理想流体
流动时减少了
(U ux ) d y
0
δ
0.99U ux
把这些流量折合成理想
流体流动通过一个厚度 1
δ
的流量，这个厚度就叫做
1
位移厚度。
根据定义
u
1 = (1
0
x )d y U
y
0.99U
边界层使来流的流线
向外排挤了位移厚度的
δ
ux
距离，所以位移厚度也
u x (U
0
根据定义
u
2=
x (1 0U
ux) d y u x)d y U
显然, 2< 1
§9—4 平板边界层动量积分方程
对平板绕流的如图区域应用动量方程，进口断面选在平板前缘 处，出口断面离前缘距离为x，出口断面厚度为当地边界层厚度 δ(x)，进口断面厚度取为出口断面的δ(x)-δ1(x)，这样通过进 口断面和出口断面的流量是相等的，必有一条流线可以连接两 个断面的厚度，用它作为区域的上边界。
一侧摩擦力
Cf =
摩阻系数
1
D
= 1.328 el
U 2 (bl)
R 1/2
2
二.平板紊流边界层
平板紊流边界层兼有 壁面紊流和自由紊流的
① 粘性底层 0 < y+ < 5 ② 过渡区 5 < y+ < 70
性质，在边界层的外 区，流动特性与圆管紊 流有所不同。
③ 紊流区
+>
<
④ 不稳定区
y 0.4
由于平板首部转捩点前必有一段层流边界层，所以不存在全 程为紊流的边界层，只能是混合边界层。按全程为紊流边界层 的摩擦阻力计算应作修正。
流体力学的边界层
流体力学的边界层
哎呀呀，一听到“流体力学的边界层”这个词，是不是感觉脑袋都大啦？其实我一开始也是这样的，觉得这东西可难可难啦！
你想想啊，咱们平常看到的水流、气流，好像就是那么随意地流来流去，可谁能想到这里面还有这么个神秘的“边界层”呢？
就比如说，咱们在河里玩水，水从咱们脚边流过，感觉凉凉的、滑滑的。
可你有没有想过，靠近咱们脚的那一层水，和远处的水流动的方式不太一样呢？这靠近咱们脚的这一层，就是边界层啦！
有一次，我和小伙伴们做了一个有趣的实验。
我们弄了一个小小的水槽，里面装满了水，然后让水慢慢流出来。
我们发现，在水槽的边上，水好像流得特别慢，就像是被什么东西拖住了一样。
这难道不神奇吗？
再打个比方，就像咱们跑步的时候，身边的风呼呼地吹。
靠近咱们身体的那一层风，跑起来就没那么顺畅，这也像是一种边界层的现象呢！
那这个边界层到底有啥用呢？这可太重要啦！比如说飞机的翅膀，要是不考虑边界层的影响，飞机可能就飞不起来啦！还有汽车的外形设计，也得考虑边界层，要不然风阻大得吓人，得多费油啊！
我就好奇地问老师：“老师，这边界层咋这么神奇啊？”老师笑着说：“孩子，这世界上神奇的东西多着呢，边界层只是其中一小部分。
”我又接着问：“那怎么才能更好地研究它，让它为我们服务呢？”老师摸摸我的头说：“只要你们好好学习，以后就能明白啦！”
我心里就想，哼，我一定要把这神秘的边界层搞清楚！
你说，这流体力学的边界层是不是特别有趣？它就像是一个隐藏在我们身边的小秘密，等着我们去发现，去探索！我觉得啊，只要我们用心去观察，去学习，就能揭开它神秘的面纱，让它为我们的生活带来更多的便利和惊喜！。
流体力学中的边界层理论
流体力学中的边界层理论流体力学是研究流体运动和相互作用的学科。
在流体力学中，边界层理论是一个重要的概念，它描述了流体靠近固体壁面时的流动特性。
本文将介绍流体力学中的边界层理论，从基本原理到应用实例，全面探讨这一理论的重要性和实际价值。
一、边界层现象的定义和意义在流体力学中，边界层是指流体流动中靠近固体表面的一层，其流动特性与远离边界的无限远处的流体不同。
边界层现象的产生和发展对于很多实际问题都具有重要意义。
例如，当空气流过汽车的外表面时，边界层的存在会对气流的分离和阻力产生影响。
准确理解和掌握边界层理论，对于优化设计和改善物体运动性能具有重要作用。
二、边界层理论的基本原理1. 平衡条件边界层理论的基本假设是边界层内的流动是定常流动和局部平衡的。
在这一假设下，可以利用物理量的守恒方程和牛顿运动定律来进行分析和计算。
2. 边界层方程边界层方程是描述边界层内流体运动的关键方程组。
它包括连续性方程、动量方程和能量方程。
这些方程考虑了流体内部各个物理量的平衡和变化，并通过求解边界层方程组可以得到流体在边界层内的运动状态。
3. 粘性效应粘性是边界层理论考虑的一个重要因素。
由于流体的粘性特性，边界层会出现剪切应力和速度剖面变化。
这些粘性效应对于固体表面的摩擦力和阻力产生重要影响，因此必须在边界层理论中加以考虑。
三、边界层理论的应用实例1. 空气动力学在航空航天工程中，边界层理论被广泛应用于翼型设计和气动力分析。
通过准确计算边界层内的流动特性，可以优化飞行器的升力和阻力性能，提高飞行效率。
2. 水力学在水力学领域，边界层理论被用于河流和水泥工程的设计和分析。
通过控制边界层内的水流运动，可以减小底摩擦阻力，提高水流的输送能力。
3. 汽车工程在汽车设计中，边界层理论被用于研究车体表面的空气流动。
通过优化车体形状和减小边界层厚度，可以降低空气阻力，提高汽车的燃油经济性。
四、结语流体力学中的边界层理论是研究流体流动与固体界面相互作用的重要理论框架。
平板边界层实验报告
平板边界层实验报告引言平板边界层实验是一种常见的流体力学实验方法，用于研究在流体与固体界面发生的各种现象和特性。
通过实验可以获取边界层厚度、速度剖面、摩擦系数等参数，对于理解流体边界层的特性具有重要意义。
本实验报告将详细介绍平板边界层实验的原理、实验装置、实验过程和实验结果，并对实验结果进行分析和讨论。
实验原理在实验中，我们使用平板边界层实验装置对流体的边界层进行研究。
其原理基于以下几点：1.边界层理论：边界层是指流体流动过程中处于流体与固体物体之间的一层流动区域，其特点是速度梯度较大、流动剪切应力较高。
边界层的特性对于流体的运动、传热和传质等过程具有重要影响。
2.平板边界层：平板边界层是指位于平板表面附近的边界层，它是边界层研究中最常见的情况之一。
通过对平板边界层的研究，可以深入理解边界层的结构、特性及其对流体流动的影响。
3.流动速度剖面：边界层中流体的速度随距离平板表面的距离而变化，一般呈现一定的速度剖面形态。
通过测量流体速度剖面，可以确定边界层的厚度和速度分布特性。
实验装置实验装置由以下几个主要部分组成：1.平板：平板用于产生平板边界层。
通常采用光滑的表面，材质多为金属或塑料。
2.流体：实验中常使用空气或水作为流体介质。
流体通过输送装置注入到实验装置中。
3.流量计：流量计用于精确测量流体的流量，以保证实验条件的准确性。
4.速度测量装置：速度测量装置用于测量流体在平板边界层中的速度。
常见的测量方法包括热线法、激光多普勒测速法等。
5.数据记录系统：数据记录系统用于记录实验过程中获得的各项数据，包括流体流量、速度剖面等。
实验步骤本实验的具体步骤如下：1.准备工作：清洁实验装置，确保平板表面光滑且无杂质。
2.实验装置搭建：按照实验要求搭建实验装置，包括安装平板、连接流体输送装置和速度测量装置。
3.流体注入：启动流体输送装置，将流体注入实验装置中，并调节流量控制阀以控制流体的流量。
4.测速：使用速度测量装置对流体在平板边界层中的速度进行测量。
平板湍流边界层课件
图 11 6 紊流边界层自由边界示意图
第20页/共49页
湍流平板边界层的湍动特性
为了深入理解边界层中的紊流结构，常对紊流中两 个相邻测点同时进行脉动流速的量测，以分析紊流 的空间特性。 空间相关函数（space correlation function）:
层中间歇系数 的分布规律如图
11-5所示并可用下式表示：
1 2
1
erf
5
y
0.78
（11-11）
第19页/共49页
图 11 5 紊流平板边界层间歇系数[3]
湍流平板边界层的湍动特性
边界层内紊流与边界层外势流的交界面有时称为边界 层的自由（freeboundary）。图11-6为自由边界的示意 图。自由边界随时间而变动，具有随机的性质。光滑
中无量纲量采用 u' v' 表示单位质量切应力的无量纲量。在紧靠壁面处未
U2
能量测到有关数据。
在边界层的外边界，即紊流边界 层与上部势流的交界面处紊流具 有间歇性质。克莱巴诺夫[3]测得 的资料显示，在 y 0.8 处，平
板紊流边界层即具有明显的间歇 性质，而当 y 1.2 时则流速基本 上不再呈现脉 动。平板紊流边界
一般用y 表示x2 ，认为：
粘性底层： 0 y 5 10 过渡区：5 10 y 30 70
对数区： 30 70 y , y 0.2
以上三个区域统称内区。 尾流区（外区）： 0.2 y 1.0
第8页/共49页
湍流平板边界层的流速分布与分区结构
第9页/共49页
流体力学第六章 边界层理论
流体力学第六章
流体力学第六章
Q
v
uv
u dy
udy U
y x 0 0 x
x 0
而
0
uK1
v y
dy
0
uK1
u x
dy
1 K
2
0
x
uK2dy
1 K
2
x
0
uK2dy
U K2
于是第二个积分
vuKudy
v
0
y K10 y
uK1
dyK1(x10u(dyU uK2)U dyK1UK2)
流体力学第六章
u
u x
v
u y
p x
2u y 2
已知普朗特方程组
p y
0
u x
v y
0
0
uk 1
udy x
0
ukv
udy y
p x
0
uk dy
0
uk
2u y2 dy
积分一
积分二
积分三
其中 (x)
(6 2 1)
流体力学第六章
b(x) a(x)
ddxx(x)dx
x 0
0
uk1
u y
2
dy
uk2dy Uk1
udy
k 1 x 0
k 1 x 0
p x
0
uk
dy
k
0
uk1
u y
2 dy
(6-2-3)
流体力学第六章
uk2dyUk1 udy
k1 x 0
k1x0
px0ukdyk0uk1uy2dy
(6-2-3)
上式为哥路别夫积分方程。
流体力学教案第8章边界层理论
第八章 边界层理论§8-1 边界层的基本概念实际流体和理想流体的本质区别就是前者具有粘性。
对层流而言，单位面积摩擦力的大小yud d μτ=，可以看出，对于确定的流体的等温流场，摩擦力的大小与速度梯度有关，其比例函数即动力粘度。
速度梯度yud d 大，粘性力也大，此时的流场称为粘性流场。
若速度梯度yud d 很小，则粘性力可以忽略，称为非粘性流场。
对于非粘性流场，则可按理想流体来处理。
则N-S 方程可由欧拉方程代替，从而使问题大为简化。
Vlv l lV v A y u V l tVl t u mρρμρρ======2223d d d d 粘性力惯性力当空气、蒸汽，水等小粘度的流体与其它物体作高速相对运动时，一般雷诺数很大。
由vVl==粘性力惯性力Re ，则在这些流动中，惯性力>>粘性力，所以可略去粘性力。
但在紧靠物体壁面存在一流体薄层，粘性力却与惯性力为同一数量级。
所以，在这一薄层中，两者均不能略去。
这一薄层就叫边界层，或叫速度边界层，由普朗特在1904年发现。
a ．流体流过固体壁面，紧贴壁面处速度从零迅速增至主流速度，这一流体薄层，就叫边界层或速度边界层。
b ．整个流场分为两部分 层外，0=∂∂yu，粘性忽略，无旋流动。
层内，粘性流，主要速度降在此，有旋流动。
c ．由边界层外边界上∞=V u %99，来定义δ，δ为边界层厚度。
d ．按流动状态，边界层又分为层流边界层和紊流边界层。
由于在边界层内，流体在物体表面法线方向(即yu∂∂)速度梯度很大，所以，边界层内的流体具有相当大的旋涡强度；而在层外，由于速度梯度很小。
所以，即使对于粘度很大的流体，粘性力也很小，故可忽略不计，所以可认为，图8-2空气沿平板边界层速度分布外部区域边界层边界层外的流动是无旋的势流。
边界层的基本特征有： (1)1<<Lδ⇒薄层性质，其中L 为物体的长度；沿流方向↑↑→δx 。
(2) 层内yu∂∂很大， 边界层内存在层流和紊流两种流态。
二维不可压缩定常流动(平板)边界层方程
二维不可压缩定常流动(平板)边界层方程下载温馨提示：该文档是我店铺精心编制而成，希望大家下载以后，能够帮助大家解决实际的问题。
文档下载后可定制随意修改，请根据实际需要进行相应的调整和使用，谢谢!并且，本店铺为大家提供各种各样类型的实用资料，如教育随笔、日记赏析、句子摘抄、古诗大全、经典美文、话题作文、工作总结、词语解析、文案摘录、其他资料等等，如想了解不同资料格式和写法，敬请关注!Download tips: This document is carefully compiled by the editor. I hope that after you download them, they can help you solve practical problems. The document can be customized and modified after downloading, please adjust and use it according to actual needs, thank you!In addition, our shop provides you with various types of practical materials, such as educational essays, diary appreciation, sentence excerpts, ancient poems, classic articles, topic composition, work summary, word parsing, copy excerpts, other materials and so on, want to know different data formats and writing methods, please pay attention!二维不可压缩定常流动(平板)边界层方程引言在流体力学领域，研究流体在固体表面附近的行为至关重要。
平板边界层实验
基本实验一（物理概念类）：平板边界层实验
通过本实验认识固壁上存在流动边界层，理解边界层内的速度分布特点。
当实际流体以较大的速度U 绕物体表面流动时，直接与固壁接触的流体速度为零，然后以一定的速度梯度增加，逐步过渡到层外势流速度U , 这一流体剪切层叫边界层。
气体绕平板作定常流动时，边界层沿流动方向在平板上的变化如下图所示。
边界层沿平板渐增厚，开始是层流，经过一段距离后层流转变成湍流。
当地雷诺数是决定流态和流态转换的一个无量纲参数 Re x Ux ν
= （1）
上式中x 为从平板前缘点算起的位置坐标，ν为流体的运动粘度。
从层流转变为湍流的临界
雷诺数一般是在5×105 ~ 3×106 范围内。
1. 边界层厚度
(1) 名义厚度δ定义为在边界层的外边界流速达到外部势流速度U 的99%时的厚度。
对平板层流边界层和湍流边界层（5Re 10x ≤），其名义厚度分别为 1/54/55.0
0.382l t x x U
U ννδδ⎛⎫== ⎪⎝⎭ （2） (2) 位移厚度δ*设平板边界层内速度分布为()u y ，位移厚度定义为
0(1)d u y U δδ*=-⎰ （3）
(3) 动量损失厚度θ
0(1)d u u y U U δθ=-⎰ （4）。
流体力学第十章边界层理论
梯度较大，一定距离后尾迹逐渐扩散，速度梯度减小，最终消失在主流
区中。 （三）渐扩管中的流动
图10.5所示为渐扩管中的流动。由于 流道截面逐渐增大，主流区中压力不断增
尾迹区
高，流体便需要消耗动能来补充压力能，
但是在边界层中由于粘性摩擦力的影响而 损失的动能较主流区大，因此其动能不足
图10.4 绕过流线型机翼的流动
动是层流，称为层流边界层，受粘性力的控制。
当流体沿平板继续流动，边界层逐渐增厚，
扰动便会发展起来，边界层中的流动变成紊 W 层流区 过渡区
紊流区
流，此时边界层厚度 增加很快，称为紊流
w
w
边界层。边界层由层流向紊流转变时，不是
突然发生的，中间有一过渡区，称作变流区。
δ
在与板面直接接触的地方，还有一层极薄的 δ
层流底层（对光滑板尤其明显）。边界层由
层流向紊流的转变，取决于雷诺数Re 的大小。
x W
层流底层
对绕流流场， Re与主流流速 W 、流体运动粘 度 和自板端向后流过的距离 x 有关，即
图10.7 流体绕过流线型锐端平板
Re Wx Wx
（10.1）
第10页
退出 返回
第十章 边界层理论
第一节 边界层特性
在绕流流场中，边界层的流动同样也有由层流转入紊流的现象。如 图10.7所示为处在均速主流流场中的流线型锐端平板。刚接触板端时，流 速 W 是均匀的。进入平板后，由于粘性作用，在壁面处便出现一层极薄的 边界层。
第9页
退出 返回
第十章 边界层理论
第一节 边界层特性
因为边界层厚度 极小，扰动在其中不易发展，所以此时边界层中的流
第十章 边界层理论
第一节 边界层特性
平板边界层实验报告
流体力学实验平板边界层实验报告班级姓名实验日期指导教师北京航空航天大学流体力学研究所流体力学实验平板边界层实验报告一、实验目的测定平板边界层内的流速分布，并比较层流边界层及紊流边界层的速度分布的差别。
二、实验设备本实验使用的是一个二维开路闭口低速风洞，在该风洞实验段中装有两块平板，以分别测量层流及紊流边界层的速度分布。
为测量速度分布，在平板板面上安装有总压排管及静压管。
这些测压管分别用橡皮管连接到多管压力计上，通过测量多管压力计液柱高度推算出速度来，具体原理见后。
为测出实验段风速，在实验段侧壁上装有风速管，风速管的总压孔及静压孔也分别用橡皮管连接于多管压力计上，装备情况见图1。
图1三、实验原理当气流流过平板时由于粘性作用使紧贴平板表面处的流速为零，离开板面速度就逐渐增大，最后达到相当于无粘时的气流速度。
对平板来说，就等于来流速度了。
由于空气粘性很小，只要来流速度不是很小时，流速变化大的区域只局限在靠近板面很薄的一层气流内，这一薄层气流通常叫作边界层。
人为地规定，自板面起，沿着它的法线方向，至达到99%无粘时的速度处的距离，称为边界层厚度δ。
不可压流场中，每一点处的总压P 0，等于该点处的静压和动压122ρv 之和。
p p v 0212=+ρ 则 v p p =-20()ρ（1）因此只需测出边界层内各点处的静压p ，总压p 0，就可计算出各点的速度来。
但考虑到垂直平板方向的静压梯度等于零（即∂∂p y /=0），我们只需在平板表面开一静压孔，所测的静压就等于该点所在的平板法线方向上各点的静压。
要测边界层内的速度分布就只要测出沿平板法线上各点的总压即可。
p i 0──为各测点的总压。
p i ──为各测点的静压。
v i ──为各测点的速度。
γ ──为多管压力计所使用的液体重度（公斤/米3）。
∆h i ──为各测点总压管与静压管的液柱高度差。
ρ ──为空气的密度，实验时可依据当时室温及大气压强由表查出。
平板壁面上流体边界层求解的数学模型
青海大学综合型设计实验文献综述平板壁面上流体边界层求解的数学模型姓名：王晓明学号： 0821105128院系：化工学院化学工程系专业：化学工程与工艺班级： 081指导教师：李晓昆老师摘要：我们生活中常见的一些流体，如水、空气，它们的粘性力相对于其他形式的力，如内力、重力和压力，都是相对较小的，看起来似乎是可以被忽略的。
然而在1744年的时候，d' Alembert发现，如果我们考虑流体时把粘性完全忽略，那么理论上得到的结果和实验中的结果会相差非常大。
在1904年的第三届世界数学家大会上，Prantle.L提出了边界层的概念，也为我们解决边界层的问题提出了一种方法。
关键词：边界层形成重要性理论依据边界层发展边界层（boundary layer）是高雷诺数绕流中紧贴物面的粘性力不可忽略的流动薄层，又称流动边界层、附面层。
这个概念由近代流体力学的奠基人，德国人Ludw ig Prandtl于（普朗特）1904年首先提出。
从那时起，边界层研究就成为流体力学中的一个重要课题和领域。
在边界层内，紧贴物面的流体由于分子引力的作用，完全粘附于物面上，与物体的相对速度为零。
一、边界层的形成流体在大雷诺数下作绕流流动时，在离固体壁面较远处，粘性力比惯性力小得多，可以忽略；但在固体壁面附近的薄层中,粘性力的影响则不能忽略,沿壁面法线方向存在相当大的速度梯度，这一薄层叫做边界层。
流体的雷诺数越大，边界层越薄。
从边界层内的流动过渡到外部流动是渐变的，所以边界层的厚度δ通常定义为从物面到约等于99％的外部流动速度处的垂直距离，它随着离物体前缘的距离增加而增大。
根据雷诺数的大小，边界层内的流动有层流与湍流两种形态。
一般上游为层流边界层，下游从某处以后转变为湍流，且边界层急剧增厚。
层流和湍流之间有一过渡区。
当所绕流的物体被加热（或冷却）或高速气流掠过物体时，在邻近物面的薄层区域有很大的温度梯度，这一薄层称为热边界层。
二、边界层理论物体在雷诺数很大的流体中以较高的速度相对运动时，沿物体表面的法线方向，得到如图（1）所示的速度分布曲线。
工程流体力学 第六版 第7章 边界层理论
1
2
1+ ? 0
1
? ~ 即：y ,
2 y 2 y y y x
y2 x2 y x
y
x y
x
x
x 2 x 2x
y x2 y2
2
12
1
1ε
1ε
1 1
1
1 12
1
2
简化N-S方程：
x
x x
y
x y
1
p x
v(
2 x x 2
2 x y 2
)
1
11
ε
1
1
1 （2 12
1
7.1.1 边界层概念 7.1.2 边界层内的流态
7.1 边界层概念
边界层：（1904年，第三届国际数学家学会，普朗特第一次提出）
实际流体绕过物体流动时，由于流体粘性的影响在物 体表面附近形成沿面的法线方向速度变化很快的薄层。
常见绕流现象
飞机/汽车阻力、 炮弹/球体飞行、 建筑、叶片绕流...
y 无黏性区
Fsx
p x
(
δ 0
ρυxdy )dx
动量：e
x
(
0
x dy )dx
e 边界层外边界上的速度
平板: υ∞ 曲面:υe(x)
流出动量 -流入动量 =
x
( δ 0
ρυx2dy )dx
υe
x
(
δ 0
ρυxdy )dx
➢ x方向的表面力：
AB面： p
y A
p 1 p dx
dl 2 x
θ
C
d
BD面： τwdx
即：(
p y
0)
，
流体力学——平板边界层编程
对于本次编程编程作业，小组运用matlab 和c++两种程序对平板边界层问题和绕过楔形体边界层流动问题进行分析研究。
以下是运用matlab 解决问题的过程。
一、 平板边界层问题该问题可以归结为在已知边界层条件下解一个高阶微分方程,即解0''5.0'''=+ff f 。
Matlab 提供了解微分方程的方法，运用换元法将高阶微分方程降阶，然后运用“ode45”函数进行求解。
函数其难点在于如何将边界条件中1',→∞→f η运用好，由四阶龙格-库塔方法知其核心是换元试算匹配，故在运用函数时通过二分法实现1',→∞→f η是可行的。
程序如下：第一问m 函数function dy = rigid(t,y) dy = zeros(3,1); dy(1) = y(2); dy(2) = y(3);dy(3) = -0.5*y(1)*y(3);%第一问main 程序[T,Y] = ode45('rigid',[0 5],[0 0 0]);plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+') %二分法试算f ’’的初始值以满足f ’趋向无穷时的边界条件，图像上可以清晰看出f ’无穷时的结果 >> [T,Y] = ode45('rigid',[0 5],[0 0 1]); plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+') >> [T,Y] = ode45('rigid',[0 5],[0 0 0.5]); plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+') >> [T,Y] = ode45('rigid',[0 5],[0 0 0.25]); plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+')>> [T,Y] = ode45('rigid',[0 5],[0 0 0.375]); plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+') >> grid on>> [T,Y] = ode45('rigid',[0 5],[0 0 0.3125]); plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+') >> grid on>> [T,Y] = ode45('rigid',[0 5],[0 0 0.34375]); plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+') grid on>> [T,Y] = ode45('rigid',[0 5],[0 0 0.328125]); plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+') grid on>> [T,Y] = ode45('rigid',[0 10],[0 0 0.328125]);%当f ’’为0.328125时，逼近结果已经很好，在0到5的变化范围内已经非常接近精确解 plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+') grid on>> [T,Y] = ode45('rigid',[0 5],[0 0 0.335975]); plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+') grid on选取f’’=0.335975时的数据展示，T代表η的变化，Y的第一二三列代表ff,f,'''为形象展示所得结果，如图-1所示：二、楔形绕流边界层问题解决绕过楔形的边界层流动问题与平板问题的不同之处在于微分方程变为0α，其中α可取1，并不失其普遍性，由于β+ffβ+ff'''2=)'1(''-小于-0.199时会发生分离，该微分方程不再适用，故β取值大于-0.199。
第七章 边界层及其基本计算
人们常采用近似解法，其中应用的较为广泛的是边界层动量积 分方程解法。
海南大学机电学院工程流体Biblioteka 学7.3 边界层的动量积分方程
二、边界层动量积分方程的推导
ve
➢ 粘性不可压流体绕物体定常二维流动；
➢ 取微元ABCDA，对其x方向的动量变化及 受力情况分析。
海南大学 第七章 粘性流体绕物体的流动
主编：孙文策 教师：马庆芬
第六章 粘性流体管内流动
1 边界层概念 2 层流边界层微分方程 3 边界层动量积分方程 4 平板层流边界层的计算 5 圆管内流动的边界层 6 边界层分离与卡门涡街
7 绕流阻力与阻力系数
海南大学机电学院
工程流体力学
7.1 边界层概念
一、边界层定义
工程流体力学
第六章 粘性流体管内流动
1 边界层概念 2 层流边界层微分方程 3 边界层动量积分方程 4 平板层流边界层的计算 5 圆管内流动的边界层 6 边界层分离与卡门涡街
7 绕流阻力与阻力系数
海南大学机电学院
工程流体力学
7.7 绕流阻力与阻力系数
一、绕流阻力的组成
绕流阻力
摩擦阻力
流体作用在物体表面上的切应力引起。
7.1 边界层概念
二、边界层的基本特征
与物体的特征长度相比，边界层的厚度很小, x
边界层内沿厚度方向，存在很大的速度梯度。 边界层厚度沿流体流动方向是增加的。 由于边界层很薄，可以近似认为边界层中各截面上的压强等于
同一截面上边界层外边界上的压强值。
在边界层内，黏性力与惯性力同一数量级。 边界层内的流态，也有层流和紊流两种流态。
7 绕流阻力与阻力系数
平板层流边界层的近似计算
平板层流边界层的近似计算§8-4平板层流边界层的近似计算作为应用边界层的积分关系式来决实际问题的例子，下面我们来研究不可压粘性流体定常流流经平板的问题。
如图所示：设x轴沿着平板，y轴为平板法线方向。
坐标原点在平板前缘点上，来流的沿x轴，板长为l。
假定来流流经平板时，平板上下两层形成层流边界层，如图所示。
现在要求的是边界的厚度的变化规律和摩擦阻力F D。
由于顺来流方向放置的平板很薄，可以认为不引起流动的改变。
所以，在边界层外边界上，，由势流的伯努利方程：两边对x求导，则：即：p=常数，即边界层外边界上压力为常数。
而边界层内，。
所以整个边界层内向点压力相同。
即整个流场压力处处相等。
代入上式则变成：(1)(1)式中有三个未知数u，，δ，所以再补充两个方程。
①当x固定时，假设边界层内速度u的分布为：(2)可以看出层内随y↑—>u↑，这和实际情况是符合的。
边界条件：1) 壁面外，y=0，u=0；2) 边界层外边界处，y=δ，u= V∞；3) 边界层外边界处，y=δ，；4) 边界层外边界处，由于u=V∞，由层流边界层微分方程(即普朗特边界层方程),在边界层的外边界上：5) 在平板壁面处，y=0，u=υ=0，又由上式(普朗特边界层方程)，得：；把边界条件代入(2)式，得：再把上面的五个系数代入(2)式，得第一个补充关系式，即层流边界层中的速度分布规律为：再对上式求导，并利用牛顿内摩擦定律，得：(3) 再将上式代入(1)式求积分，则得到：(4)(5) 将(3)，(4)，(5)代入(1)式，得：，积分得：确定积分常数C，x=0， =0，C=0，于是得：，它的精确解为，并且的表达式为的三次方时，得出的解比四次方精确。
其系数为4.64。
因此，不能认为选择速度分布时，多项式数越多越好。
由上式可看出：x—>δ；V∞—>δ↓。
将δ表达式，代入(c)式，得切向应力：从上式可以看出：沿平板长度方向(x方向)，越来越小，这是因随x，速度边界层越来越厚，边界层内速度变化渐趋缓和之故。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
对于本次编程编程作业，小组运用matlab 和c++两种程序对平板边界层问题和绕过楔形体边界层流动问题进行分析研究。
以下是运用matlab 解决问题的过程。
一、 平板边界层问题该问题可以归结为在已知边界层条件下解一个高阶微分方程,即解0''5.0'''=+ff f 。
Matlab 提供了解微分方程的方法，运用换元法将高阶微分方程降阶，然后运用“ode45”函数进行求解。
函数其难点在于如何将边界条件中1',→∞→f η运用好，由四阶龙格-库塔方法知其核心是换元试算匹配，故在运用函数时通过二分法实现1',→∞→f η是可行的。
程序如下：第一问m 函数function dy = rigid(t,y) dy = zeros(3,1); dy(1) = y(2); dy(2) = y(3);dy(3) = -0.5*y(1)*y(3);%第一问main 程序[T,Y] = ode45('rigid',[0 5],[0 0 0]);plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+') %二分法试算f ’’的初始值以满足f ’趋向无穷时的边界条件，图像上可以清晰看出f ’无穷时的结果 >> [T,Y] = ode45('rigid',[0 5],[0 0 1]); plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+') >> [T,Y] = ode45('rigid',[0 5],[0 0 0.5]); plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+') >> [T,Y] = ode45('rigid',[0 5],[0 0 0.25]); plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+')>> [T,Y] = ode45('rigid',[0 5],[0 0 0.375]); plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+') >> grid on>> [T,Y] = ode45('rigid',[0 5],[0 0 0.3125]); plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+') >> grid on>> [T,Y] = ode45('rigid',[0 5],[0 0 0.34375]); plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+') grid on>> [T,Y] = ode45('rigid',[0 5],[0 0 0.328125]); plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+') grid on>> [T,Y] = ode45('rigid',[0 10],[0 0 0.328125]);%当f ’’为0.328125时，逼近结果已经很好，在0到5的变化范围内已经非常接近精确解 plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+') grid on>> [T,Y] = ode45('rigid',[0 5],[0 0 0.335975]); plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+') grid on选取f’’=0.335975时的数据展示，T代表η的变化，Y的第一二三列代表ff,f,'''为形象展示所得结果，如图-1所示：二、楔形绕流边界层问题解决绕过楔形的边界层流动问题与平板问题的不同之处在于微分方程变为0α，其中α可取1，并不失其普遍性，由于β+ffβ+ff'''2=)'1(''-小于-0.199时会发生分离，该微分方程不再适用，故β取值大于-0.199。
其解法与平板绕流类似，在平板绕流基础上增加β的变化即可。
程序如下：第二问m函数function dz = rigid1(t,z)dz = zeros(3,1);t=-0.199; %其中t表示β的变化，可由-0.199至无穷大dz(1) = z(2);dz(2) = z(3);dz(3) = -z(1)*z(3)-t*(1-z(2)^2);第二问main函数[T,Y] = ode45('rigid1',[0 5],[0 0 0]);plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+') %二分法试算f’’的初始值以满足f’趋向无穷时的边界条件，图像上可以清晰看出f’无穷时的结果>> [T,Y] = ode45('rigid1',[0 5],[0 0 1]);plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+')>> [T,Y] = ode45('rigid1',[0 5],[0 0 0.5]);plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+')>> [T,Y] = ode45('rigid1',[0 5],[0 0 0.25]);plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+')>> [T,Y] = ode45('rigid1',[0 5],[0 0 0.15]);plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+')>> grid on>> [T,Y] = ode45('rigid1',[0 5],[0 0 0.05]);plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+')>> grid on>> [T,Y] = ode45('rigid1',[0 5],[0 0 0.025]);plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+')grid on%当f’’为0.025时，逼近结果已经很好，在0到5的变化范围内已经非常接近精确解plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+')grid on>> [T,Y] = ode45('rigid1',[0 5],[0 0 0]);plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+')grid on选取β=-0.199,f’’=0.025时的数据进行展示，程序代码中T代表η的变化，y(i)的代表''f的值。
第二个问题的数据如下表：f,f,'当β变化时，只需要更改主程序中t的值即可，再运用main程序解微分方程即可得到相应的解。
三、优化楔形绕流问题的算法初值的选取对于减少计算量有很大帮助，首先考虑运用拟合得到经验性公式，尝试如下：在β变化时按此式（其中t代表β）f''=-0.2877t^2+1.1892t+0.3915变化f’’相应的初值利于快速得到结果。
该式子是由二阶拟合多组数据得到的结果。
具体程序为x = [-0.1988 -0.19 -0.18 -0.16 -0.14 -0.10 0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.80 1.00 1.20 1.60 2.00];y = [0 0.0860 0.1285 0.1905 0.2395 0.3191 0.4696 0.5870 0.6869 0.77680.8542 0.9277 0.9960 1.120 1.233 1.336 1.521 1.687];A = polyfit(x,y,2)z = polyval(A,x);plot(x,y,'k+',x,z,'r')%f''=0.2133*t^3-0.8160*t^2+1.4229*t+0.4189 三阶公式%f''=0.2826*t^5-1.4303*t^4+2.6049*t^3-2.1971*t^2+1.5037*t+0.4723 五阶公式图像如图-3因此，在β发生变化时只需要在主程序中相应的变化t值即可。
例如取β=1.2时，由f''=-0.2877t^2+1.1892t+0.3915式子可得初值为1.4030,更改m函数中t值为1.2，再调用main函数，并赋予f’’新的初值，如下：[T,Y] = ode45('rigid1',[0 5],[0 0 1.3199]);plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+')grid on应用过程中发现f’并不能较快的趋近于1，即使运用四次拟合公式（f''=0.2826*t^5-1.4303*t^4+2.6049*t^3-2.1971*t^2+1.5037*t+0.4723）也不能较快的得到收敛的结果，其主要原因是精度达不到要求和拟合本身的误差过大。