借宝地写几个小短文，介绍CFD的一些实际的入门知识。主要是因为这里支持Latex，写起来比较方便。
CFD，计算流体力学，是一个挺难的学科，涉及流体力学、数值分析和计算机算法，还有计算机图形学的一些知识。尤其是有关偏微分方程数值分析的东西，不是那么容易入门。大多数图书，片中数学原理而不重实际动手，因为作者都把读者当做已经掌握基础知识的科班学生了。所以数学基础不那么好的读者往往看得很吃力，看了还不知道怎么实现。本人当年虽说是学航天工程的，但是那时本科教育已经退步，基础的流体力学课被砍得只剩下一维气体动力学了，因此自学CFD的时候也是头晕眼花。不知道怎么实现，也很难找到教学代码——那时候网络还不发达，只在教研室的故纸堆里搜罗到一些完全没有注释，编程风格也不好的冗长代码，硬着头皮分析。后来网上淘到一些代码研读，结合书籍论文才慢慢入门。可以说中间没有老师教，后来赌博士为了混学分上过CFD专门课程，不过那时候我已经都掌握课堂上那些了。
回想自己入门艰辛，不免有一个想法——写点通俗易懂的CFD入门短文给师弟师妹们。本人不打算搞得很系统，而是希望能结合实际，阐明一些最基本的概念和手段，其中一些复杂的道理只是点到为止。目前也没有具体的计划，想到哪里写到哪里，因此可能会很零散。但是我争取让初学CFD 的人能够了解一些基本的东西，看过之后，会知道一个CFD代码怎么炼成的（这“炼”字好像很流行啊）。欢迎大家提出意见，这样我尽可能的可以追加一些修改和解释。
言归正传，第一部分，我打算介绍一个最基本的算例，一维激波管问题。说白了就是一根两端封闭的管子，中间有个隔板，隔板左边和右边的气体状态（密度、速度、压力）不一样，突然把隔板抽去，管子内面的气体怎么运动。这是个一维问题，被称作黎曼间断问题，好像是黎曼最初研究双曲微分方程的时候提出的一个问题，用一维无粘可压缩Euler方程就可以描述了。
这里
这个方程就是描述的气体密度、动量和能量随时间的变化（）
与它们各自的流量（密度流量，动量流量，能量流量
）随空间变化（）的关系。
在CFD中通常把这个方程写成矢量形式
这里
进一步可以写成散度形式
一定要熟悉这种矢量形式
以上是控制方程，下面说说求解思路。可压缩流动计算中，有限体积（FVM）是最广泛使用的算法，其他算法多多少少都和FVM有些联系或者共通的思路。了解的FVM，学习其他高级点的算法（比如目前比较热门的间断有限元、谱FVM、谱FDM）就好说点了。
针对一个微元控制体，把Euler方程在空间积分
用微积分知识可以得到
也就是说控制体内气体状态平均值的变化是控制体界面上流通量的结果。
因此我们要计算的演化，关键问题是计算控制体界面上的。
FVM就是以这个积分关系式出发，把整个流场划分为许多小控制体，每个控制体和周围相邻的某个控制体共享一个界面，通过计算每个界面上的通量来得到相邻控制体之间的影响，一旦每个控制体的变化得到，整个流场的变化也就知道了。
所以，再强调一次，关键问题是计算控制体界面上的。
初学者会说，这个不难，把界面上的插值得到，然后就可以计算。有道理！
咱们画个图，有三个小控制体i-1到i+1，中间的“|”表示界面，控制体i
右边的界面用表示，左边的就是。
| i-1 | i | i+1 |
好下个问题：每个小控制体长度都是如何插值计算界面上的
？
最自然的想法就是：取两边的平均值呗，
但是很不幸，这是不行的。
那么换个方法？直接平均得到？
还是很不行，这样也不行。
我靠，这是为什么？这明明是符合微积分里面的知识啊？
这个道理有点复杂，说开了去可以讲一本书，可以说从50年代到70年代，CFD科学家就在琢磨这个问题。这里，初学者只需要记住这个结论：对于流动问题，不可以这样简单取平均值来插值或者差分。如果你非要想知道这个究竟，我现在也不想给你讲清楚，因为我眼下的目的是让你快速上手，而且该不刨根问底的时候就不要刨根问底，这也是初学阶段一种重要的学习方法。
好了，既然目的只是为了求，我在这里，只告诉你一种计算方法，也是非常重要、非常流行的一种方法。简单的说，就是假设流动状态在界面
是不连续的，先计算出界面两边的值，和，再由
它们用某种方法计算出。上述方法是非常重要的，是由一个苏联人Godunov在50年代首创的，后来被发展成为通用Godunov方法，著名的ENO/WENO就是其中的一种。
好了，现在的问题是：
1 怎么确定和
2 怎么计算
对于第一个问题，Godunov在他的论文中，是假设每个控制体中是均
匀分布的，因此
第二个问题，Godunov采用了精确黎曼解来计算。什么是“精确黎曼解”，就是计算这个激波管问题的精确解。既然有精确解，那还费功夫搞这些FVM算法干什么？因为只有这种简单一维问题有精确解，稍微复杂一点就不行了。精确解也比较麻烦，要分析5种情况，用牛顿法迭代求解（牛顿法是什么？看数值计算的书去，哦，算了，现在暂时可以不必看）。这是最初Godunov的方法，后来在这个思想的基础上，各种变体都出来
了。也不过是在这两个问题上做文章，怎么确定，怎么计算。
Godunov假设的是每个小控制体内是均匀分布，也就是所谓分段常数(piecewise constant)，所以后来有分段线性(picewise linear)或者分段二次分布(picewise parabolic)，到后来ENO/WENO出来，那这个假设的多项式次数就继续往上走了。都是用多项式近似的，这是数值计算中的一个强大工具，你可以在很多地方看到这种近似。
Godunov用的四精确黎曼解，太复杂太慢，也不必要，所以后来就有各种近似黎曼解，最有名的是Roe求解器、HLL求解器和Osher求解器，都是对精确黎曼解做的简化。
这个多项式的阶数是和计算精度密切相关的，阶数越高，误差就越小。不过一般来说，分段线性就能得到不错的结果了，所以工程中都是用这个，Fluent、Fastran以及NASA的CFL3D、OverFlow都是用这个。而黎曼求解器对精度的影响不是那么大，但是对整个算法的物理适用性有影响，也就是说某种近似黎曼求解器可能对某些流动问题不合适，比如单纯的Roe对于钝头体的脱体激波会算得乱七八糟，后来加了熵修正才算搞定。
上次（https://www.doczj.com/doc/777255862.html,/node/399）说到了求解可压缩流动的一个重要算法，通用Godunov方法。其两个主要步骤就是
1 怎么确定和
2 怎么计算
这里我们给出第一点一个具体的实现方法，就是基于原始变量的MUSCL 格式（以下简称MUSCL）。它是一种很简单的格式，而且具有足够的精
度，NASA著名的CFL3D软件就是使用了这个格式，大家可以去它的主页（https://www.doczj.com/doc/777255862.html,/Cfl3dv6/cfl3dv6.html）上看手册，里面空间离散那一章清楚的写着。
MUSCL假设控制体内原始变量（就是）的分布是一次或者二次多项式，
如果得到了这个多项式，就可以求出控制体左右两个界面的一侧的值
和。
我们以压力p为例来说明怎么构造这个多项式。这里我只针对二次多项式来讲解，你看完之后肯定能自己推导出一次多项式的结果（如果你搞不定，那我对你的智商表示怀疑）。
OK，开始
假设，这个假设不影响最终结论，因为你总可以把一个区间线性的变换到长度为1的区间。
假设压力p在控制体i内部的分布是一个二次多项式
，控制体i的中心处于处，左右两个界面就是
和。
这里先强调一个问题，在FVM中，每个控制体内的求解出来的变量实
际上是这个控制体内的平均值。
所以，
。
我们知道，和，等距网格情况下和处的导数可以近似表示为
那么
由上述三个有关a，b和c的方程，我们可以得到
这样就可以得到f(x)的表达式了，由此可以算出和
通常MUSCL格式写成如下形式
对应我们的推导结果（二次多项式假设）。
但是这不是最终形式。如果直接用这个公式，就会导致流场在激波（间断）附近的振荡。因为直接用二次多项式去逼近一个间断，会导致这样的效果。所以科学家们提出要对间断附近的斜率有所限制，因此引入了一个非常重要的修改——斜率限制器。加入斜率限制器后，上述公式就有点变化。
这里是Van Albada限制器
是一个小数（），以防止分母为0。
密度和速度通过同样的方法来搞定。
密度、速度和压力被称作原始变量，所以上述方法是基于原始变量的MUSCL。此外还有基于特征变量的MUSCL，要复杂一点，但是被认为适合更高精度的格式。然而一般计算中，基于原始变量的MUSCL由于具有足够的精度、简单的形式和较低的代价而被广泛使用。
OK，搞定了。下面进入第二点，怎么求。关于这一点，我不打算做
详细介绍了，直接使用现有的近似黎曼解就可以了，都是通过和
计算得到。比如Roe因为形式简单，而非常流行。在CFL3D软件主页（https://www.doczj.com/doc/777255862.html,/Cfl3dv6/cfl3dv6.html）上看手册，附录C的C.1.3。
想了一下，还是把Roe求解器稍微说说吧，力求比较完整。但是不要指望我把Roe求解器解释清楚，因为这个不是很容易三言两语说清的。
Roe求解器的数学形式是这样的
显然这个公式的第一项是一个中心差分形式，先前说过简单的中心差分不可行，原因是耗散不足导致振荡，说得通俗点就像一个弹簧，如果缺乏耗
散（阻尼）它就会一直振荡。“耗散”这个术语在激波捕捉格式中是最常见的。第二项的作用就是提供足够的耗散了。
这里和已经用MUSCL求得了，的定义在第一讲中已经介绍了。只有是还没说过的。
这个矩阵可以写成特征矩阵和特征向量矩阵的形式
而
，，和的具体表达式在许多书上都有，而且这里的矩阵表达有问题，所以就不写了。
是由、和代入计算得到。而、和采用所谓Roe平均值
这才是Roe求解器关键的地方！
总结一下，就是用Roe平均计算界面上的气体状态，然后计算得到
，这样就可以得到了。如果有时间，我后面会找一个代码逐句分析一下。
总之，计算还是很不直接的。构造近似黎曼解是挺有学问的，需要对气体动力学的物理和数学方面有较深的理解。通常，如果不是做基础研究，你只需要知道它们的特点，会用它们就可以了，而不必深究它们怎么推导出来的。
1．流体的连续介质模型：研究流体的宏观运动，在远远大于分子运动尺度的范围里考察流体运动，而不考虑个别分子的行为，因此我们可以把流体视为连续介质。 它有如下性质： （1）流体是连续分布的物质，它可以无限分割为具有均布质量的宏观微元体。 （2）不发生化学反应和离解等非平衡热力学过程的运动流体中，微元体内流体状态服 从热力学关系 （3）除了特殊面外，流体的力学和热力学状态参数在时空中是连续分布的，并且通常 认为是无限可微的 2．应力：有限体的微元面积上单位面积的表面力称为表面力的局部强度，又称为应力，定义如下：=n T A F A δδδlim 0→ 3．流体的界面性质：微元界面两侧的流体的速度和温度相等，应力向量的大小相等.方向相反或应力分量相等。 4．流体具有易流行和压缩性。 5．应力张量具有对称性。 6．欧拉描述法：在任意指定的时间逐点描绘当地的运动特征量（如速度、加速度）及其它的物理量的分布（如压力、密度等）。 7．拉格朗日描述法：从某个时刻开始跟踪质点的位置、速度、加速度和物理参数的变化，这种方法是离散质点的运动描述法称为拉格朗日描述法。 8.流线：速度场的向量线，该曲线上的任意一点的切向量与当地的的速度向量重合。 迹线：流体质点点的运动迹象。 差别：迹线是同一质点在不同时刻的位移曲线。 流线是同一时刻、不同质点连接起来的速度场向量线。 流线微分方程：ω dz v dy u dx == 迹线微分方程：t x U i i ??= 9．质点加速度：质点速度向量随时间的变化率。 U U t U a )(??+??= 质点加速度=速度的局部导数+速度的迁移导数。 物理量的质点导数=物理量的局部导数+物理量的对流导数。
基于CFX的离心泵 内部流场数值模拟基于CFX的离心泵内部流场数值模拟 随着计算流体力学和计算机技术的快速发展，泵内部的流动特征成为热点研究方向，目前应用 CFX 软件的科研人员还较少，所以将CFX 使用的基本过程加以整理供初学者参考。如有不对之处敬请指教。 、CFX数值计算的完整流程 、基于ICEM CFD勺离心泵网格划分 2.1导入几何模型 2.2修整模型 2.3创建实体 2.4仓U建PRAT 2.5设置全局参数 2.6划分网格 2.7检查网格质量并光顺网格2.8导出网格—选择求解器2.9导出网格 、CFX-Pre设置过程 3.1基本步骤 3.2新建文件
3.3导入网格 3.4定义模拟类型3.5创建计算域3.6指定边界条件3.7建立交界面
3.8定义求解控制 3.9定义输出控制 3.10写求解器输入文件 3.11定义运行 3.12计算过程 四、CFX-Post 后处理 4.1计算泵的扬程和效率 4.2云图 4.3矢量图 4.4流线图 2.1导入几何模型 在ICEMCFD软件界面内，单击File宀Imort Geometry^STEP/IGES（—般将离心泵装配文件保存成STEP格式）， 将离心泵造型导入I C E M如图3所示。 图3导入几何模型界面
2.2 修整模型 单击Geometry^Repair Geometry 宀Build Topology，设置Tolerenee，然后单击Apply，如图 4 所示。拓扑 分析后生成的曲线颜色指示邻近表面的关系：gree n =自由边，yellow =单边，red =双边，blue =多边，线条 颜色显示的开/关Model tree T Geometry T Curves T Color by cou nt,Red curves 表示面之间的间隙在容差之 内，这是需要的物理模型, N41 f !孕ECHH 匚丁E> !1 Z-和-1 ：z? ...... ....................... 兰直卤* 百曲gw 卜宀-im * Q涕曲空JIT^J 厂社tt-sfri- Piwpe^ifl-5 CorFklr air^ i Cphcri s Quip^jr 匸* JO 匸叭和皈X XWM X ■an. y% wn- Yellow edges 通常是一些需要修补的几何。 亠 图4修整模型界面 2-3 创建实体单击Geometry^Creade Body,详细过程如图5所示。
《流体力学》典型例题（9大类） 例1~例3——牛顿内摩擦定律（牛顿剪切公式）应用 例4~例5——流体静力学基本方程式的应用——用流体静力学基本方程和等压面计算某点的压强或两点之间的压差。 例6~例8——液体的相对平衡——流体平衡微分方程中的质量力同时考虑重力和惯性力（补充内容） （1）等加速直线运动容器中液体的相对平衡（与坐标系选取有关） （2）等角速度旋转容器中液体的平衡（与坐标系选取有关） 例9——求流线、迹线方程；速度的随体导数（欧拉法中的加速度）；涡量计算及流动有旋、无旋判断 例10~16——速度势函数、流函数、速度场之间的互求 例17——计算流体微团的线变形率、角变形率及旋转角速度 例18~20——动量定理应用（课件中求弯管受力的例子） 例21~22——总流伯努利方程的应用 例23——综合：总流伯努利方程、真空度概念、平均流速概念、流态判断、管路系统沿程与局部损失计算 例题1：如图所示，质量为m ＝5 kg 、底面积为S ＝40 cm ×60 cm 的矩形平板，以U ＝1 m/s 的速度沿着与水平面成倾角θ＝ 30 的斜面作等速下滑运动。已知平板与斜面之间的油层厚度δ ＝1 mm ，假设由平板所带动的油层的运动速度呈线性分布。 求油的动力粘性系数。 U G=mg δ θ 解：由牛顿内摩擦定律，平板所受的剪切应力du U dy τμ μδ =＝ 又因等速运动，惯性力为零。根据牛顿第二定律： 0m ==∑F a ，即： gsin 0m S θτ-?= ()32 4gsin 59.8sin 301100.1021N s m 1406010m U S θδμ--?????==≈????? 粘性是流体在运动状态下，具有的抵抗产生剪切变形速率能力的量度；粘性是流体的一种固有物理属性；流体的粘性具 有传递运动和阻滞运动的双重性。 例题2：如图所示，转轴的直径d ＝0.36 m ，轴承的长度l ＝1 m ，轴与轴承的缝隙宽度δ＝0.23 mm ，缝隙中充满动力粘性系数0.73Pa s μ=?的油，若轴的转速200rpm n =。求克服油的粘性阻力所消耗的功率。 δ d l n 解：由牛顿内摩擦定律，轴与轴承之间的剪切应力 ()60d d n d u y πτμ μδ =＝ 粘性阻力（摩擦力）：F S dl ττπ=?=
计算流体力学教案 Teaching plan of computational fluid mechanics
计算流体力学教案 前言：本文档根据题材书写内容要求展开，具有实践指导意义，适用于组织或个人。便于学习和使用，本文档下载后内容可按需编辑修改及打印。 一、流体地基本特征 1.物质地三态 在地球上，物质存在地主要形式有：固体、液体和气体。 流体和固体地区别:从力学分析地意义上看，在于它们对外力抵抗地能力不同。 固体：既能承受压力，也能承受拉力与抵抗拉伸变形。 流体：只能承受压力，一般不能承受拉力与抵抗拉伸变形。 液体和气体地区别：气体易于压缩；而液体难于压缩； 液体有一定地体积，存在一个自由液面；气体能充满任意形状地容器，无一定地体积，不存在自由液面。 液体和气体地共同点：两者均具有易流动性，即在任何 微小切应力作用下都会发生变形或流动，故二者统称为流体。 2.流体地连续介质模型
微观：流体是由大量做无规则运动地分子组成地，分子之间存在空隙，但在标准状况下，1cm3液体中含有3.3×1022个左右地分子，相邻分子间地距离约为3.1×10-8cm。1cm3气体中含有2.7×1019个左右地分子，相邻分子间地距离约为3.2×10-7cm。 宏观：考虑宏观特性，在流动空间和时间上所采用地一切特征尺度和特征时间都比分子距离和分子碰撞时间大得多。 （1）概念 连续介质（continuum/continuous medium）：质点连续充满所占空间地流体或固体。 连续介质模型（continuum continuous medium model）：把流体视为没有间隙地充满它所占据地整个空间地一种连续介质，且其所有地物理量都是空间坐标和时间地连续函数地一种假设模型：u =u（t,x,y,z）。 （2）优点 排除了分子运动地复杂性。物理量作为时空连续函数，则可以利用连续函数这一数学工具来研究问题。 3.流体地分类
第二章计算流体力学的基本知识 流体流动现象大量存在于自然界及多种工程领域中，所有这些工程都受质量守恒、动量守恒和能量守恒等基本物理定律的支配。这章将首先介绍流体动力学的发展和流体力学中几个重要守恒定律及其数学表达式，最后介绍几种常用的商业软件。 2.1计算流体力学简介 2.1.1计算流体力学的发展 流体力学的基本方程组非常复杂，在考虑粘性作用时更是如此，如果不靠计算机，就只能对比较简单的情形或简化后的欧拉方程或N-S方程进行计算。20 世纪30～40 年代，对于复杂而又特别重要的流体力学问题，曾组织过人力用几个月甚至几年的时间做数值计算，比如圆锥做超声速飞行时周围的无粘流场就从1943 年一直算到1947 年。 数学的发展，计算机的不断进步，以及流体力学各种计算方法的发明，使许多原来无法用理论分析求解的复杂流体力学问题有了求得数值解的可能性，这又促进了流体力学计算方法的发展，并形成了"计算流体力学" 。 从20 世纪60 年代起，在飞行器和其他涉及流体运动的课题中，经常采用电子计算机做数值模拟，这可以和物理实验相辅相成。数值模拟和实验模拟相互配合，使科学技术的研究和工程设计的速度加快，并节省开支。数值计算方法最近发展很快，其重要性与日俱增。 自然界存在着大量复杂的流动现象，随着人类认识的深入，人们开始利用流动规律来改造自然界。最典型的例子是人类利用空气对运动中的机翼产生升力的机理发明了飞机。航空技术的发展强烈推动了流体力学的迅速发展。 流体运动的规律由一组控制方程描述。计算机没有发明前，流体力学家们在对方程经过大量简化后能够得到一些线形问题解读解。但实际的流动问题大都是复杂的强非线形问题，无法求得精确的解读解。计算机的出现以及计算技术的迅速发展使人们直接求解控制方程组的梦想逐步得到实现，从而催生了计算流体力
基于CFX的离心泵内部流场数值模拟 基于CFX的离心泵内部流场数值模拟 随着计算流体力学和计算机技术的快速发展，泵内部的流动特征成为热点研究方向，目前应用CFX 软件的科研人员还较少，所以将CFX使用的基本过程加以整理供初学者参考。如有不对之处敬请指教。 一、 CFX数值计算的完整流程 二、基于ICEM CFD的离心泵网格划分 2.1 导入几何模型 2.2 修整模型 2.3 创建实体 2.4 创建PRAT 2.5 设置全局参数 2.6 划分网格 2.7 检查网格质量并光顺网格 2.8 导出网格－选择求解器 2.9 导出网格 三、CFX-Pre 设置过程 3.1 基本步骤 3.2 新建文件 3.3 导入网格 3.4 定义模拟类型 3.5 创建计算域 3.6 指定边界条件 3.7 建立交界面 3.8 定义求解控制
3.10 写求解器输入文件 3.11 定义运行 3.12 计算过程 四、 CFX-Post后处理 4.1 计算泵的扬程和效率 4.2 云图 4.3 矢量图 4.4 流线图 2.1 导入几何模型 在ICEM CFD软件界面内，单击File→Imort Geometry→STEP/IGES（一般将离心泵装配文件保存成STEP格式），将离心泵造型导入ICEM，如图3所示。 图3 导入几何模型界面 2.2 修整模型 单击Geometry→Repair Geometry→Build Topology，设置Tolerence，然后单击Apply，如图4所示。拓扑分析后生成的曲线颜色指示邻近表面的关系：green = 自由边， yellow = 单边，red = 双边， blue =多边，线条
第一章 流体及其主要物理性质 例1： 已知油品的相对密度为0.85，求其重度。 解： 例2： 当压强增加5×104Pa 时，某种液体的密度增长0.02%，求该液体的弹性系数。 解： 例3： 已知：A ＝1200cm 2，V ＝0.5m/s μ1＝0.142Pa.s ，h 1＝1.0mm μ2＝0.235Pa.s ，h 2＝1.4mm 求：平板上所受的内摩擦力F 绘制：平板间流体的流速分布图 及应力分布图 解：（前提条件：牛顿流体、层流运 动） 因为 τ1＝τ2 所以 3 /980085.085.0m N ?=?=γδ0=+=?=dV Vd dM V M ρρρρρ d dV V -=Pa dp d dp V dV E p 84105.2105% 02.01111?=??==-==ρρβdy du μ τ=??????? -=-=?2221110 h u h u V μτμτs m h h V h u h u h u V /23.02 112212 2 11 =+= ?=-μμμμμN h u V A F 6.41 1=-==μ τ
第二章 流体静力学 例1： 如图，汽车上有一长方形水箱，高H ＝1.2m ，长L ＝4m ，水箱顶盖中心有一供加水用的通大气压孔，试计算当汽车以加速度为3m/s 2向前行驶时，水箱底面上前后两点A 、B 的静压强（装满水）。 解： 分析：水箱处于顶盖封闭状态，当加速时，液面不变化，但由于惯性力而引起的液体内部压力分布规律不变，等压面仍为一倾斜平面，符合 等压面与x 轴方向之间的夹角 例2： （1）装满液体容器在顶盖中心处开口的相对平衡 分析：容器内液体虽然借离心惯性力向外甩，但由于受容器顶限制，液面并不能形成旋转抛物面，但内部压强分布规律不变： 利用边界条件：r ＝0，z ＝0时，p ＝0 作用于顶盖上的压强： （表压） （2）装满液体容器在顶盖边缘处开口的相对平衡 压强分布规律： =+s gz ax g a tg = θPa L tg H h p A A 177552=??? ?? ?+==θγγPa L tg H h p B B 57602=??? ?? ?-==θγγC z g r p +-?=)2( 2 2ωγg r p 22 2ωγ =C z g r p +-?=)2( 2 2ω γ
借宝地写几个小短文，介绍CFD的一些实际的入门知识。主要是因为这里支持Latex，写起来比较方便。 CFD，计算流体力学，是一个挺难的学科，涉及流体力学、数值分析和计算机算法，还有计算机图形学的一些知识。尤其是有关偏微分方程数值分析的东西，不是那么容易入门。大多数图书，片中数学原理而不重实际动手，因为作者都把读者当做已经掌握基础知识的科班学生了。所以数学基础不那么好的读者往往看得很吃力，看了还不知道怎么实现。本人当年虽说是学航天工程的，但是那时本科教育已经退步，基础的流体力学课被砍得只剩下一维气体动力学了，因此自学CFD的时候也是头晕眼花。不知道怎么实现，也很难找到教学代码——那时候网络还不发达，只在教研室的故纸堆里搜罗到一些完全没有注释，编程风格也不好的冗长代码，硬着头皮分析。后来网上淘到一些代码研读，结合书籍论文才慢慢入门。可以说中间没有老师教，后来赌博士为了混学分上过CFD专门课程，不过那时候我已经都掌握课堂上那些了。 回想自己入门艰辛，不免有一个想法——写点通俗易懂的CFD入门短文给师弟师妹们。本人不打算搞得很系统，而是希望能结合实际，阐明一些最基本的概念和手段，其中一些复杂的道理只是点到为止。目前也没有具体的计划，想到哪里写到哪里，因此可能会很零散。但是我争取让初学CFD 的人能够了解一些基本的东西，看过之后，会知道一个CFD代码怎么炼成的（这“炼”字好像很流行啊）。欢迎大家提出意见，这样我尽可能的可以追加一些修改和解释。
言归正传，第一部分，我打算介绍一个最基本的算例，一维激波管问题。说白了就是一根两端封闭的管子，中间有个隔板，隔板左边和右边的气体状态（密度、速度、压力）不一样，突然把隔板抽去，管子内面的气体怎么运动。这是个一维问题，被称作黎曼间断问题，好像是黎曼最初研究双曲微分方程的时候提出的一个问题，用一维无粘可压缩Euler方程就可以描述了。 这里 这个方程就是描述的气体密度、动量和能量随时间的变化（）与它们各自的流量（密度流量，动量流量，能量流量 ）随空间变化（）的关系。 在CFD中通常把这个方程写成矢量形式 这里 进一步可以写成散度形式
汽车外部气体流动模拟 振动和噪声控制研究所 1.模型概述 在汽车外部建立一个较大的长方体几何空间，长度约为30m，宽度和高度约为5m，在空间内部挖出汽车形状的空腔，汽车尺寸参照本田CRV为4550mm*1820mm*1685mm。由于汽车向前开进，气体从车头流向车尾，因此将汽车前方空间设为气体入口，后方空间设为气体出口，模拟气体在车外的流动。另外为了节省计算成本将整个模型按1:100的比例缩小，考虑到模型和流体均是对称的，因此仅画出几何模型的一半区域，建立对称面以考虑生成包含理想气体的流体域。在Catia中建立的模型如图1.1所示。 图1.1几何模型 2.利用ICEM CFD进行网格划分 a)导入有Catia生成的stp格式的模型； b)模型修复，删除多余的点、线、面，允许公差设为0.1； c)生成体，由于本模型仅为流体区域，因此将全部区域划分为一个体，选取方法可以 使用整体模型选取； d)为了后面的设置边界方便，因此将具有相同特性的面设为一个part,共设置了in， out，FreeWalls，Symmetry和Body； e)网格划分，设置Max element=2，共划分了1333817个单元，有225390个节点； f)网格输出，设置求解器为ANSYS CFX，输出cfx5文件。 3.利用ANSYS CFX求解 a)生成域，物质选定Air Ideal Gas，参考压强设为1atm，浮力选项为无浮力模型，
域运动选项为静止，网格变形为无；流体模型设定中的热量传输设定为Isothermal，流体温度设定为288k，湍流模型设定为Shear Stress Transport模型，壁面函数 选择Automatic。 b)入口边界设定，类型为Inlet，位置选定在in，质量与栋梁选定Normal Speed，设 定为15m/s，湍流模型设定类型为Intensity and Length Scale=0.05，Eddy Len.Scale=0.1m。 c)出口边界设定，边界类型为Outlet，位置选out。质量与动量选项为Static Pressure，相对压强为0pa。 d)壁面边界设定，边界类型为Wall，位置选在FreeWalls。壁面边界详细信息中指定 WallInfluence On Flow为Free Slip。 e)对称边界设定，边界类型为Symmetry，位置选在Symmetry。 f)汽车外壁面设定，边界类型为Wall，位置设在Body，壁面详细信息选项中指定Wall Influence On Flow为No Slip，即汽车壁面为无滑移壁面。 g)初始条件设定，初始速度分量设为U方向为15m/s，其他两个方向的速度为零。 h)求解设置，残差类型选为RMS，残差目标设定为1e-5，当求解达到此目标时，求解 自动终止。求解之前的模型如图3.1所示。 图3.1求解之前的模型 4.结果后处理 从图4.1中可以看出计算收敛。
1.简述流体力学有哪些研究方法和优缺点? 实验方法就是运用模型实验理论设计试验装置和流程，直接观察流动现象，测量流体的流动参数并加以分析和处理，然后从中得到流动规律。实验研究方法的优点：能够直接解决工程实际中较为复杂的流动问题，能够根据观察到的流动现象，发现新问题和新的原理，所得的结果可以作为检验其他方法的正确性和准确性。实验研究方法的缺点主要是对于不同的流动需要进行不同的实验，实验结果的普遍性稍差。 理论方法就是根据流动的物理模型和物理定律建立描写流体运动规律的封闭方程组以及相应初始条件和边界条件，运 用数学方法准确或近似地求解流场，揭示流动规律。理论方法的优点是：所得到的流动方程的解是精确解，可以明确地给出各个流动参数之间的函数关系。解析方法的缺点是：数学上的困难比较大，只能对少数比较简单的流动给出解析解，所能得到的解析解的数目是非常有限的。 数值方法要将流场按照一定的规则离散成若干个计算点，即网格节点；然后，将流动方程转化为关于各个节点上流动 参数的代数方程；最后，求解出各个节点上的流动参数。数值方法的优点是：可以求解解析方法无能为力的复杂流动。数值方法的缺点是：对于复杂而又缺乏完整数学模型的流动仍然无能为力，其结果仍然需要与实验研究结果进行对比和验证。 2.写出静止流体中的应力张量,解释其中非0项的意义. 无粘流体或静止流场中，由于不存在切向应力，即p ij =0（i ≠j ），此时有 P =00000 0xx yy zz p p p ??????????=000000p p p -????-????-??=-p 00000011????1?????? = -p I 式中I 为单位张量，p 为流体静压力。 流体力学中，常将应力张量表示为 p =-+P I T （2-9） 式中p 为静压力或平均压力，由于其作用方向与应力定义的方向相反，所以取负值；T 称为偏应力张量，即 T =xx xy xz yx yy yz zx zy zz τττττττττ?????????? （2-10） 偏应力张量的分量与应力张量各分量的关系为：i ＝j 时，p ij 为法向应力，τii = p ij - p ；当i ≠j 时p ij 为粘性剪切应力，τij =p ij 。τii =0的流体称为非弹性流体或纯粘流体，τii ≠0的流体称为粘弹性流体。 3.分析可压缩(不可压缩)流体和可压缩(不可压缩)流动的关系. 当气体速度流动较小（马赫数小于0.3)时，其密度变化不大，或者说对气流速度的变化不十分敏感，气体的压缩性没有表现出来。因此，在处理工程实际问题时，可以把低速气流看成是不可压缩流动，把气体可以看作是不可压缩流体。而当气体以较大的速度流动时，其密度要发生明显的变化，则此时气体的流动必须看成是可压缩流动。 流场任一点处的流速v 与该点（当地)气体的声速c 的比值，叫做该点处气流的马赫数，用符号Ma 表示： Ma /v c v == （4－20） 当气流速度小于当地声速时，即Ma<1时，这种气流叫做亚声速气流；当气流速度大于当地声速时，即Ma>l 时，这种气流称为超声速气流；当气流速度等于当地声速时，即Ma=l 时，这种气流称为声速气流。以后将会看到，超声速气流和亚声速气流所遵循的规律有着本质的不同。 马赫数与气流的压缩性有着直接的联系。由式(4－11)可得 所以有 222Ma d ρv dv dv ρc v v =-=-。 （4－21） 当Ma≤0.3时，dρ/ρ≤0.09dv /v 。由此可见，当速度变化一倍时，气体的密度仅仅改变9%以下，一般可以不考虑密度的变化，即认为气流是不可压缩的。反之，当Ma>0.3时，气流必须看成是可压缩的。 4.试解释为什么有时候飞机飞过我们头顶之后才能听见飞机的声音. 5.试分析绝能等熵条件下截面积变化对气流参数(v ,p ,ρ,T )的影响.
第26卷第2期河北科技大学学报Vol.26,No.2 2005年6月Journal of Hebei University of Science and T echnology June2005 文章编号:100821542(2005)022******* 计算流体力学(CFD)的通用软件 翟建华 (河北科技大学国际交流与合作处,河北石家庄050018) 摘要:对化学工程领域中的通用CFD(Computational Fluid Dynamics)模拟软件Phoenics,Flu2 ent,CFX等的具体特点和应用情况进行了综述,指出了他们各自的结构特点、特有模块、包含的数学模型和成功应用领域;给出了选用CFD软件平台的7项准则,对今后CFD技术的发展进行了预测,指出,今后CFD研究的主要方向将集中在数学模型开发、工程改造和新设备开发及与工艺软件的匹配连用等方面。 关键词:计算流体力学;模拟软件;CFX;FLUENT;PH OENICS 中图分类号:T Q015.9文献标识码:A Review of commercial CFD software ZH AI Jian2hua (Department of Int ernation Exchange and Cooperation,H ebei University of Science and Technology,Shijiazhuang H ebei 050018,China) Abstr act:The paper summar izes the features and application of the CF D simulation software like Phoenics,F luent and CFX etc in chemical engineering,and discusses their str ucture features,special modules,mathematical models and successful application areas.It also puts forward seven r ules for the good choice of commercial CF D code for the CF D simulation resea rcher s.Based on t he predict ion of the technology development,it points out the possible r esear ch direction for CF D in the future will focus on the development of mathematical model,project transformat ion,new equipment and their matching application with technologi2 cal softwa re. Key words:CF D;simulation software;CF X;FLUENT;P HOENICS CFD(Computational Fluid Dynamics)软件是计算流体力学软件的简称,是用来进行流场分析、计算、预测的专用工具。通过CFD模拟,可以分析并且显示流体流动过程中发生的现象,及时预测流体在模拟区域的流动性能,并通过各种参数改变,得到相应过程的最佳设计参数。CFD的数值模拟,能使我们更加深刻地理解问题产生的机理,为实验提供指导,节省以往实验所需的人力、物力和时间,并对实验结果整理和规律发现起到指导作用。随着计算机软硬件技术的发展和数值计算方法的日趋成熟,出现了基于现有流动理论的商用CFD软件。这使许多不擅长CFD工作的其他专业研究人员能够轻松地进行流体数值计算,从而使研究人员从编制繁杂、重复性的程序中解放出来,以更多的精力投入到研究问题的物理本质、问题提法、边界(初值)条件和计算结果的合理解释等重要方面上,充分发挥商用CFD软件开发人员和其他专业研究人员各自的智力优势,为解决实际工程问题开辟了道路。 CFD研究走过了相当漫长的过程。早期数值模拟阶段,由于缺乏模拟工具,研究者一般根据自身工作性质和研究过程,自行编制模拟程序,其优点是针对性强,对具体问题的解决有一定精度,但是,带来的问题 收稿日期:2004208221;修回日期:2004211221;责任编辑:张军 作者简介:翟建华(19642),男,河北平乡人,教授,主要从事化工CFD、高效传质与分离和精细化工方面的研究。
高等流体力学 第一章 流体力学的基本概念 连续介质：流体是由一个紧挨着一个的连续的质点所组成的，没有任何空隙的连续体，即所 谓的连续介质。 流体质点：是指微小体积内所有流体分子的总和。 欧拉法质点加速度：时变加速度与位变加速度和 z u u y u u x u u t u dt du a x z x y x x x x x ??+??+??+??== 质点的随体导数：质点携带的物理量随时间的变化率称为质点的随体导数，用dt d 表示。在欧拉法描述中的任意物理量Q 的质点随体导数表述如下： x k k Q u t Q dt dQ ??+??= 式中Q 可以是标量、矢量、张量。质点的随体导数公式对任意物理量都成立，故将质点的 随体导数的运算符号表示如下： x k k u t dt d ??+??= 其中 t ?? 称为局部随体导数，x k k u ??称为对流随体导数，即在欧拉法描述的流动中，物理 量的质点随体导数等于局部随体导数与对流随体导数之和。 体积分的随体导数：质点携带的物理量随时间的变化率称为质点的随体导数。则在由流体质点组成的流动体积V 中标量函数Φ（x, t ）随时间的变化率就是体积分的随导函数。 由两部分组成①函数Φ 对时间的偏导数沿体积V 的积分，是由标量场的非恒定性引起的。②函数Φ通过表面S 的通量。由体积V 的改变引起的。 ()dV divv dt d dV v div t dS u dV t dV dt d v v n s v v ?? ? ???Φ+Φ=??????Φ+?Φ?=Φ+?Φ?=Φ??????????????()dV adivv dt da dV av div t a dS au dV t a adV dt d v v n s v v ?? ????+=??????+??=+??=?????????????? 变形率张量： 11ε 12ε13ε D ij = 21ε 22ε 23ε 31ε 32ε 33ε
CFD 实验报告一 姓名： 学号： 一、题目： 利用中心差分格式近似导数22/dx y d ，数值求解常微分方程 x dx y d 2sin 2 2= （10≤≤x ） 00==x y 4 2 s i n 11- ==x y 步长分别取x ?＝0.05, 0.01, 0.001，0.0001。 二、报告要求： 1）列出全部计算公式和步骤； 2）表列出程序中各主要符号和数组意义； 3）绘出数值计算结果的函数曲线，并与精确解比较； 4）比较不同差分格式和不同网格步长计算结果的精度和代价； 5）附源程序。 三、相关差分格式 二阶导数22/dx y d 的三点差分格式有向前差分、向后差分和中心差分，表达 式分别如下： ()()()22122 22122 211 222 222j j j j j j j j j u u u u O x x x u u u u O x x x u u u u O x x x ++--+--+?=+???-+?=+???-+?=+???一阶向前差分：一阶向后差分：二阶中心差分： 代入微分方程可以得到差分方程，表达式分别如下： 212 212 11 2 2=sin 22=sin 22=sin 2j j j j j j j j j j j j u u u x x u u u x x u u u x x ++--+--+?-+?-+?一阶向前差分：一阶向后差分：二阶中心差分： 对于三种差分格式，差分格式可以改写成AY b =的形式，其中A 是相同的，
非齐次项b 不同，如下所示： 2112112112A -????-?? ??=?? -?? ??-?? 系数矩阵 ()()()02112 3221sin 2sin 2sin 2k k y x x b x x x x y ---???? ??? ? ?=????? ?? ?-?? 一阶向前差分 ()()()()2202 322121sin 2sin 2sin 2sin 2k k x x y x x b x x x x y -???-????? ??=????? ???-?? 一阶向后差分 ()()()()2102 232 2211sin 2sin 2sin 2sin 2k k x x y x x b x x x x y --?? ?-????? ??=??????? ?-? ? 二阶中心差分 求解AY b =可以得到各节点y 的值[]T 1 22 1k k Y y y y y --= 。 四、计算公式和步骤； 1.关于精确解的推导： 已知22sin 2d y x dx =，对 x 进行两次积分，得到121 sin 24 y x C x C =-++，再结合 边界条件00 ==x y 和4 2 sin 11-==x y 得到相对应的1C 和2C ，确定最后精确解为： 1 sin 24y x x =-+。 2.关于数值求解方法： 对于方程组AY b =可直接求解，也可以使用追赶法求解，下面介绍简单追赶法求解三对角方程组的过程。
《高等流体力学》考试大纲 一、考试性质 《高等流体力学》是我校相关专业博士入学专业基础课考试科目。 二、考试形式与试卷结构 1、答卷方式：闭卷，笔试 2、答题时间；180分钟 3、题型比例 概念20% 计算与应用80% 4、参考书目 《高等流体力学》高学平，天津大学出版社，2005. 《高等工程流体力学》张鸣远等，西安交通大学出版社，2006. 三、考试要点 1、流体力学的基本概念 连续介质、欧拉法质点加速度、质点随体导数、体积分的随体导数、变形率张量、旋转角速度、判断有旋流与无旋流、涡量与速度环量的关系、应力张量的概念（包括切应力的特性、压应力的特性）、牛顿流体的本构方程（本构方程的概念、切应力和法向应力与变形的关系）。 2、流体运动的基本方程 微分形式的连续方程的表达形式、不可压缩流体的确切定义、理解其含义。N-S方程的各种表示形式、流体的能量包括哪几种形式，
并对各种形式进行解释，写出单位质量流体能量的表达式、流体运动微分形式的基本方程组有哪些方程组成，通常有几个未知量，方程组是否封闭、对于不可压缩流体，如何求解速度场、压强场以及温度场，说明其求解步骤。 3、势流运动 势流运动控制方程及求解步骤；势流求解常用的方法有哪些。速度势函数与流函数；复势与复速度；恒定平面势流的解析方法有哪几种途径；保角变换法的思路。 4、粘性流体运动 基本方程及求解途径；黏性流体运动的基本性质；黏性流体运动的解析解（如两平行板间的层流、普阿塞流的流速分布的推导）、小雷诺数流动近似解的思路；边界层的概念；边界层厚度（名义厚度、位移厚度）；边界层方程的相似性解的概念；边界层的分离现象。5、紊流运动 紊流的特征及分类；壁面剪切紊流的发生过程及紊流结构；时间平均法和系综平均法的概念。紊流运动方程—雷诺方程的推导思路，雷诺方程的形式及与N-S方程的区别，雷诺应力项的意义。紊流模型的用途，紊流模型通常有哪几类（零方程模型、一方程模型、二方程模型、其他模型）；紊流动能k、能量耗散率ε。 6、涡旋运动 涡旋的运动学性质、涡旋运动的基本方程；涡旋的形成。
流体力学计算题
水 水银 题型一：曲面上静水总压力的计算问题（注：千万注意方向，绘出压力体） 1、AB 曲面为一圆柱形的四分之一，半径R=0.2m ，宽度（垂直纸面）B=0.8m ，水深H=1.2m ，液体密度 3/850m kg =ρ，AB 曲面左侧受到液体压力。求作用在AB 曲面上的水平分力和铅直分力。（10分） 解：（1）水平分力： RB R H g A h P z c x ?- ==)2(ργ…….(3分) N 1.14668.02.0)2 2 .02.1(8.9850=??-??=，方向向右（2 分）。 （2）铅直分力：绘如图所示的压力体，则 B R R R H g V P z ??? ? ????+-==4)(2πργ……….（3分） 1.15428.04 2.014.32.0)2.02.1(8.98502=???? ? ?????+?-??=，方向向下（2分）。 2．有一圆滚门，长度l=10m ，直径D=4.2m ，上游水深H1=4.2m ，下游水深H2=2.1m ，求作用于圆滚门上的水平和铅直分压力。 l d Q h G B A 空 气 石 油 甘 油 7.623.66 1.52 9.14m 1 1
解题思路：（1）水平分力： l H H p p p x )(2 12 22121-=-=γ 方向水平向右。 （2）作压力体，如图，则 l D Al V p z 4 432 πγγγ? === 方向垂直向上。 3.如图示，一半球形闸门，已知球门的半径m R 1= ，上下游水位差m H 1= ，试求闸门受到的水平分力和竖直分力 的大小和方向。 解： (1)水平分力： ()2R R H A h P c πγγ?+===左,2R R A h P c πγγ?='=右 右左P P P x -= kN R H 79.30114.31807.92=???=?=πγ， 方向水平向右。 （２）垂直分力： V P z γ=，由于左、右两侧液体对曲面所形成的压力体均为半球面，且两侧方向相反，因而垂直方向总的压力为0。 4、密闭盛水容器，已知h 1=60cm,h 2=100cm ，水银测压计读值cm h 25=?。试求半径R=0.5m 的半球盖AB 所受总压力的水平分力和铅垂分力。
流体力学题库
B1 流体及其物理性质 1. 按连续介质的概念，流体质点是指( ) A.流体的分子 B.流体内的固体颗粒 C.几何的点 2.液体的粘性主要来源于液体（）。 A. 分子热运动； B. 分子间内聚力； C. 易变形性； D. 抗拒变形的能力。 3.当水的压强增加1个大气压时，水的密度约增大（）。 A. 1/20000； B. 1/10000； C. 1/4000； D. 1/2000 4.（）理想流体就是粘性为常数的流体； 5.（）当很小的切应力作用于流体时，流体是否流动还需要看其它条件； 6.（）当很小的切应力作用于流体时，流体不一定会流动； 7. 分别说明气体和液体的粘度与温度的关系及原因。 8. 写出牛顿粘性定律的表达式，指出各符号代表的变量及单位，说明什么是牛顿流体？ 9. 一圆锥体绕其中心轴作等角速度旋转ω=16/ rad s，锥体与固定壁面间的距离δ=1mm，用μ=0.1Pa s?的润滑油充满间隙，锥底半径R=0.3m，高H=0.5m。求作用于圆锥体的阻力矩。
10. 两无限大平行平板，保持两板的间距δ=0.2 mm 。板间充满锭子油，粘度为μ=0.01Pa ﹒s ，密度为ρ =800 kg/m 3。设下板固定，上板以U=0.5 m/s 的速度滑移，油内沿板垂直方向y 的速度u(y)为线性分布，试求： （1）锭子油的运动粘度υ； （2）上下板的粘性切应力τ1，τ2。 B2 流动分析基础 1. 非定常流动中,流线与迹线( )。 A.一定重合 B.一定不重合 C.特殊情况下可能重合 D.一定正交 2. 用欧拉法表示流体质点的加速度a r 等于（ ） A.2 2d r dt r B.u t ??r C.()u u ??r r D.u t ??r +()u u ??r r 3. 两根管径相同的圆管，以同样速度输送水和空气，不会出现（ ）情况。 A. 水管内为层流状态，气管内为湍流状态； B. 水管、气管内均为层流状态； C. 水管内为湍流状态，气管内为层流状态； D. 水管、气管内均为湍流状态。 4. 均匀流是（ ） A.当地加速度为零 B.迁移加速度为零 C.向心加速度为零 D.合加速度为零 5.（ ）若流体为稳定流动，则0dv dt =。 6.（ ）流动过程中，只要Re ＞2300，则流动为湍流。
《计算流体力学》教学大纲 课程编号：00002067 课程中文名称：计算流体力学 课程英文名称：Computational Fluid Mechanics 总学时：48实验学时：0 上机学时：0 学分：3 适用专业：工程力学 一、课程性质、目的和任务（300字内） 计算流体力学是工程力学专业的一门选修课。它的主要任务是通过教学环节，使学生理解和掌握计算流体力学的一些基本理论和基本计算方法。学会运用计算流体力学、计算机解决一些流体数值计算问题。为从事工程技术工作、科学研究及开拓新领域，打下坚实的基础。 通过本课程的学习，使学生掌握计算流体力学有关的基本概念与基本理论，学会将数学、计算机、流体力学知识有机地结合，对科学研究和工程应用中的流动问题进行数值模拟的方法。教学中应贯彻启发式，互动式。着重讲清基本概念、基本理论、分析问题的思路和方法，并配以适当讨论，逐步培养学生具有分析问题和解决问题的能力。指导学生阅读参考书、文献和资料，培养学生自学获取知识的能力。重视实践环节，要求学生上机计算，并分析计算结果。充分利用现代化的教学手段，使学生具有灵活运用知识、进行探索和开拓创新等方面的技能，并具有较高的综合素质。
二、课程教学内容及学时分配 1、理论讲授 3、教学环节的安排及学时分配 教材：1. 吴子牛主编的《计算流体力学基本原理》，科学出版社，2001年第一版。 2. 王福军主编《计算流体动力学》，清华大学出版社，2005年第一版。
参考资料： 1.R. Peyret & T. D. Taylor. Computational Methods for Floud Flow. 2.C.A.J.Fletcher. Computational Techniques for Fluid Dynamics. 3. H K Versteeg的《计算流体动力学导论》，世界图书出版公司，1995年第一版
高等工程流体力学 粘性流动 康顺 华北电力大学能源与动力工程系学院 Kangs@https://www.doczj.com/doc/777255862.html,
内容提纲 ?边界层及其方程 ?层流边界层流动转捩 ?湍流边界层结构 ?流动分离、二次流动与旋涡 能源动力领域流动问题的主要特征 ?全三维 ?非定常 ?粘性 ?高雷诺数，边界层 ?边界层：层流、转捩、湍流（紊流），分离流动，旋涡运动 叶轮机械（透平和压气机等）大多由单个或多个级组成。每个级含有一 排静子叶片列和一排转子叶片列。在级内的气流场中，一般至少有以下 几种流动现象发生：1、前缘马蹄涡；2、通道涡；3、顶部间隙涡；4、 边界层转捩；5、叶片尾迹；6、旋涡、尾迹等与叶片列周期性非定常相 互作用。 ?激波、激波与边界层相互作用
边界层流动 边界层 边界层概念：粘性很小的流体以大雷诺数运动时，在大部分流场上可以略去粘性的作用；但在物面附近的很薄的一层流体内必须考虑粘性作用。这一薄层流体称为边界层。 平板边界层示意图有边界的流动图谱 如右上图所示：流动分为三个区：边界层，尾迹区，位流区（外部势流区） 二维平板的边界层微分方程 设直匀流以零迎角平行流过一块长度为的平板,如左下图所示，人为规定，当某个y处的速度达到层外自由流的99％时，这一点到物体表面的距离（即y）称为边界层在改点的厚度，记为。显然，边界层的厚度是与X有关的，所以可以写成。 平板边界层 边界层的厚度很小，满足此关系式： 在忽略质量力的前提下，粘性平面不可压流的运动方程加上连续方程是: 用边界层条件式上式，y的数值限制在边界层之内，即 υ ∞l δδ（x） δ（x）l δ（x）<< 22 22 22 22 1 () 1 () u u u p u u u t x y x x y p u t x y y x y u x y υν ρ υυυυυ υν ρ υ ? ?????? ++=-++? ??????? ? ??????? ++=-++? ??????? ? ?? +=? ???? l δ（x）<<0yδ ≤≤