第三章 一元流体动力学基础 1.直径为150mm 的给水管道，输水量为h kN /7.980，试求断面平均流速。 解：由流量公式vA Q ρ= 注意：()vA Q s kg h kN ρ=?→// A Q v ρ= 得：s m v /57.1= 2.断面为300mm ×400mm 的矩形风道,风量为2700m 3/h,求平均流速.如风道出口处断面收缩为150mm ×400mm,求该断面的平均流速 解：由流量公式vA Q = 得：A Q v = 由连续性方程知2211A v A v = 得：s m v /5.122= 3.水从水箱流经直径d 1=10cm,d 2=5cm,d 3=2.5cm 的管道流入大气中. 当出口流速10m/ 时,求 (1)容积流量及质量流量;(2)1d 及2d 管段的流速 解：(1)由s m A v Q /0049.0333== 质量流量s kg Q /9.4=ρ (2)由连续性方程： 33223311,A v A v A v A v == 得：s m v s m v /5.2,/625.021== 4.设计输水量为h kg /294210的给水管道，流速限制在9.0∽s m /4.1之间。试确定管道直径，根据所选直径求流速。直径应是mm 50的倍数。 解：vA Q ρ= 将9.0=v ∽s m /4.1代入得343.0=d ∽m 275.0 ∵直径是mm 50的倍数，所以取m d 3.0= 代入vA Q ρ= 得m v 18.1= 5.圆形风道，流量是10000m 3/h,，流速不超过20 m/s 。试设计直径，根据所定直径求流速。直径规定为50 mm 的倍数。 解：vA Q = 将s m v /20≤代入得：mm d 5.420≥ 取mm d 450= 代入vA Q = 得：s m v /5.17= 6.在直径为d 圆形风道断面上，用下法选定五个点，以测局部风速。设想用和管轴同心但不同半径的圆周，将全部断面分为中间是圆，其他是圆环的五个面积相等的部分。测点即位于等分此部分面积的圆周上，这样测得的流速代表相应断面的平均流速。（1）试计算各测点到管心的距离，表为直径的倍数。（2）若各点流速为54321u u u u u ，，，，，空气密度为ρ，求质量流量G 。
计算流体动力学(简称CFD)是建立在经典流体动力学与数值计算方法基础之上的一门新型独立学科，通过计算机数值计算和图像显示的方法，在时间和空间上定量描述流场的数值解，从而达到对物理问题研究的目的。它兼有理论性和实践性的双重特点。 第一章节 流体流动现象大量存在于自然界及多种工程领域中，所有这些过程都受质量守恒、动量守恒和能量守恒等基本物理定律的支配。本章向读者介绍这些守恒定律的数学表达式，在此基础上提出数值求解这些基本方程的思想，阐述计算流体力学的任务及相关基础知识，最后简要介绍目前常用的计算流体动力学商用软件。 计算流体动力学((Computational Fluid Dynamics简称CFD)是通过计算机数值计算和图像显示，对包含有流体流动和热传导等相关物理现象的系统所做的分析。CFD的基本思想可以归结为:把原来在时间域及空间域上连续的物理量的场，如速度场和压力场，用一系列有限个离散点上的变量值的集合来代替，通过一定的原则和方式建立起关于这些离散点上场变量之间关系的代数方程组，然后求解代数方程组获得场变量的近似值。 CFD可以看做是在流动基本方程(质量守恒方程、动量守恒方程、能量守恒方程)控制卜对流动的数值模拟。通过这种数值模拟，我们可以得到极其复杂问题的流场内各个位置上的基本物理量(如速度、压力、温度、浓度等)的分布，以及这些物理量随时间的变化情况，确定旋涡分布特性、空化特性及脱流区等。还可据此算出相关的其他物理量，如旋转式流体机械的转矩、水力损失和效率等。此外，与CAD联合，还可进行结构优化设计等。 1.1.2计算流体动力学的工作步骤 采用CFD的方法对流体流动进行数值模拟，通常包括如下步骤: (1)建立反映工程问题或物理问题本质的数学模型。具体地说就是要建立反映问题各个量之间关系的微分方程及相应的定解条件，这是数值模拟的出发点。没有正确完善的数 学模型，数值模拟就毫无意义。流体的基本控制方程通常包括质量守恒方程、动量守恒方程、能量守恒方程，以及这些方程相应的定解条件。 (2}}寻求高效率、高准确度的计算方法，即建立针对控制方程的数值离散化方法，如有限差分法、有限元法、有限体积法等。这里的计算方法不仅包括微分方程的离散化方法及求解方法，还包括贴体坐标的建立，边界条件的处理等。这些内容，可以说是c}}的核心。 (3})编制程序和进行计算。这部分工作包括计算网格划分、初始条件和边界条件的输入、控制参数的设定等。这是整个工作中花时间最多的部分。由于求解的问题比较复杂，比如Na}ier-Stakes方程就是一个讨，分复杂的非线性方程，数值求解方法在理论上不是绝对完善的，所以需要通过实验加以验证。正是从这个意义上讲.数值模拟又叫数值试验。应该指出，这部分工作不是轻而易举就可以完成的。 4})显示计算结果。计算结果一般通过图表等方式显示，这对检查和判断分析质量和结果有重要参考意义。 以上这些步骤构成了CFD数值模拟的全过程。其中数学模型的建立是理论
工程流体力学闻德课后习题答案 第五章 实际流体动力学基础 5—1设在流场中的速度分布为u x =2ax ，u y =-2ay ，a 为实数，且a >0。试求切应力τxy 、τyx 和附加压应力p ′x 、p ′y 以及压应力p x 、p y 。 解：0y x xy yx u u x y ττμ??? ?==+= ????? 24x x u p a x μμ?'=-=-?，24y y u p a y μμ?'=-=?， 4x x p p p p a μ'=+=-，4y y p p p p a μ'=+=+ 5－2 设例５－1中的下平板固定不动，上平板以速度v 沿x 轴方向作等速运动（如图 所示），由于上平板运动而引起的这种流动，称柯埃梯（Couette ）流动。试求在这种流动情况下，两平板间的速度分布。（请将 d 0d p x =时的这一流动与在第一章中讨论流体粘性时的流动相比较） 解：将坐标系ox 轴移至下平板，则边界条件为 y ＝０，0X u u ==；y h =，u v =。 由例５－1中的（11）式可得 2d (1)2d h y p y y u v h x h h μ=- - （１） 当d 0d p x =时，y u v h =，速度ｕ为直线分布，这种特殊情况的流动称简单柯埃梯流动或简单剪切流动。它只是由于平板运动，由于流体的粘滞性带动流体发生的流动。 当 d 0d p x ≠时，即为一般的柯埃梯流动，它是由简单柯埃梯流动和泊萧叶流动叠加而成，速度分布为 (1)u y y y p v h h h =-- （２） 式中2d ()2d h p p v x μ= - （３） 当p ＞０时，沿着流动方向压强减小，速度在整个断面上的分布均为正值；当p ＜０时，沿流动方向压强增加，则可能在静止壁面附近产生倒流，这主要发生p ＜－１的情况． 5－3 设明渠二维均匀（层流）流动，如图所示。若忽略空气阻力，试用纳维—斯托克斯方程和连续性方程，证明过流断面上的速度分布为2sin (2)2 x g u zh z ，单宽流量 3 sin 3 gh q 。
《计算流体动力学分析》学习报告 计算流体力学基础： 本章主要讲解流体动力学的核心思想以及流体动力学的控制方程。 1、计算流体动力学（Computational Fluid Dynamic ）基本思想：把原来在时间和空间上的连续的物理量，用一系列离散点上的变量值来代替，通过一定的原则和方式建立变量之间的代数方程式，求解之后获得变量的近似值。 2、CFD 控制方程： 质量守恒方程 0)·=?+??u t ρρ（ 动量守恒方程（Navier-Stokes 方程） Fz z y x z u w div t w F z y x y u v div t v F z y x x u u div t u zz zx zx y zy yy xy x zx yx xx +??+??+??+??-=+??+??+??+??+??-=+??+??+??+??+??-=+??τττρρρτττρρρτττρρρ)()()()()()( 能量守恒方程 T p S gradT c k div T u div t +=+??)()(T ( ρρ） S T 为粘性耗散项。 方程含有u ，v ，w ，p ，T 和ρ六个未知量，所以还需要一个方程组，才能使其封闭，而这个方程组就是联系P 和ρ的状态方程组：P=（ρ，T ）。 组分质量守恒方程（在一个系统中，可能存在质的交换，或者存在化学组分时使用。） ()s s s s S c grad D div c u div t +=+??)()(c (s ρρρ ） 为便于对控制方程进行计算和分析，对CFD 控制方程写成通用格式： ()S z z y y x x z w y v x u t S grad div u div t +??Γ??+??Γ??+??Γ??=??+??+??+??+Γ=+??)()()()()()())()(φφφφρφρφρρφφφρρφ 依次为瞬态项，对流项，扩散项和源项。 3、湍流控制方程 三维的N-S 方程无论对于层流还是湍流都是是使用的，但由于直接求解三维瞬态的控制方程，对计算机的内存和速度要求很高，因此在工程上广为采用的方法是对瞬态的N-S 方程进行实践平均处理，同时补充反应湍流特性的其他方程，例如湍动能方程以及湍流耗散率方程
不可压缩流体动力学基础 1．已知平面流场的速度分布为xy x u x +=2，y xy u y 522+=。求在点（1，-1）处流体微团的线变形速度，角变形速度和旋转角速度。 解：（1）线变形速度：y x x u x x +=??=2θ 角变形速度：()x y y u x u x y z +=???? ????+??=222 121ε 旋转角速度：()x y x u x u x y z -=???? ????-??= 222121ω 将点（1，-1）代入可得流体微团的1=x θ，1=y θ；23/z =ε；21/z =ω 2．已知有旋流动的速度场为322+=y u x ，x z u y 32+=，y x u z 32+=。试求旋转角速度，角变形速度和涡线方程。 解：旋转角速度：2 121=???? ????-??=z u y u y z x ω 角变形速度：2521=???? ????+??= z u y u y z x ε 由z y x dz dy dx ωωω==积分得涡线的方程为： 1c x y +=，2c x z += 3．已知有旋流动的速度场为22z y c u x +=，0=y u ，0=z u ，式中c 为常数，试求流场的涡量及涡线方程。
解：流场的涡量为： 旋转角速度分别为：0=x ω 则涡线的方程为：c dz dy z y +=??ωω 即c y dz z dy +-=?? 可得涡线的方程为：c c y =+22 4．求沿封闭曲线2 22b y x =+,0=z 的速度环量。（1）Ax u x =，0=y u ；（2） Ay u x =，0=y u ；（3）0=y u ，r A u =θ。其中A 为常数。 解：（1）由封闭曲线方程可知该曲线时在z =0的平面上的圆周线。 在z =0的平面上速度分布为： Ax u x =，0=y u 涡量分布为：0=z Ω 根据斯托克斯定理得：0==?z A z s dA ΩΓ （2）涡量分布为：A z -=Ω 根据斯托克斯定理得：2b A dA z A z s πΩΓ-==? （3）由于0=r u ，r A u =θ
第三章流体动力学基础复习题 一、概念部分 1、描述流体运动的方法有和；前者以为研究对象，而后者以为研究对象。 2、流体运动的几何描述有：，，和。 3、流线有什么特点？流线、脉线和迹线有什么区别和联系？ 4、流体微团基本运动形式有，和变形运动等， 而变形运动又包括和两种。 5、描述有旋运动几何要素有、和。 6、判断正误：理想流体不存在有旋运动是否正确？为什么？试举例说明。 7、表征涡流的强弱的参数有和。 8、在无涡流空间画出的封闭周线上的速度环量为。 9、简述汤姆孙定理的内容 10、速度势函数?存在的条件是什么？流函数存在的条件是什么？ 11、简述流函数的物理意义的内容，并证明。 12、流网存在的条件是什么？简述流网的性质所包含的内容？ 13、无环量圆柱绕流运动由流、流和流叠加而成，有环量的圆柱绕流运动是无环量的圆柱绕流运动与流叠加而成。 14、是驻点。通过驻点的流线一定是零流线，是否正确？为什么？零流线是。轮廓线是。 15、描述流体运动的微分方程有、和。 写出它们的表达式。 16、纳维-斯托克斯方程中的速度只能是平均速度，是否正确？为什么？ 17、写出总水头和测压管水头的表达式，并说明各项的物理意义。 18、写出总压、全压和势压得表达式，并说明各项的物理意义。 19、简述系统和控制体的定义和特点 二、计算部分 1、已知拉格朗日描述：求速度与加速度的欧拉描述 2、试判断下列流场的描述方式：并转换成另一种描述方式 3、已知用欧拉法表示的流场速度分布规律为： 试求在t=0时刻位于点（a,b)的流体质点的运动轨迹及拉格朗日法表示的速度场 4、粘性流体在半径为R 的直圆管内做定常流动。设圆管截面（指垂直管轴的平面截面）上?????==-t t be y ae x ()()?????+-=+-=-t y t x e b u e a u 1111???+=+=t y u t x u y x
第七章不可压缩流体动力学基础在前面的章节中，我们学习了理想流体和粘性流体的流动分析，按照水力学的观点，求得平均量。但是，很多问题需要求得更加详细的信息，如流速、压强等流动参数在二个或三个坐标轴方向上的分布情况。本章的内容介绍流体运动的基本规律、基本方程、定解条件和解决流体问题的基本方法。 第一节流体微团的运动分析 运动方式：①移动或单纯的位移（平移）②旋转③线性变形④角变形。位移和旋转可以完全比拟于刚体运动，至于线性变形和脚变形有时统称为变形运动则是基于液体的易流动性而特有的运动形式，在刚体是没有的。 在直角坐标系中取微小立方体进行研究。
一、平移：如果图（a ）所示的基体各角点的质点速度向量完全相同时，则构成了液体基体的单纯位移，其移动速度为z y x u u u 、、。基体在运动中可能沿直线也可能沿曲线运动，但其方位与形状都和原来一样（立方基体各边的长度保持不变）。 二、线变形：从图（b ）中可以看出，由于沿y 轴的速度分量，B 点和C 点都比A 点和D 点大了 dy y u y ??，而 y u y ??就代表1=dy 时液体基体运动时，在单位时间内沿 y 轴方向的伸长率。 x u x ??，y u y ??，z u z ?? 三、角变形（角变形速度） d d d D C A B C D B A
dt y u dy dt dy y u d x x ??=???=α dt x u dx dt dx x u d y y ??=???=β θβθα+=-d d 2 βαθd d -= ∴ 角变形： ???? ????+??=+=-=x u y u d d d y x z 212βαθαθ ?? ? ????+??= x u z u z x y 21θ ???? ????+??=y u z u z y x 21θ 四、旋转（旋转角速度） ??? ? ????-??=-=y u x u x y z 21θω ??? ? ????-??=z u y u y z x 21ω 即， ?? ? ????-??=x u z u z x y 21ω z y x u u u z y x k j i ??????= 21ω 那么，代入欧拉加速度表达式，得： z x x x x x x z y y z z y y y y y y y x z z x x z z z z z z z y x x y y x x y du u u u u u u u dt t x u u u u u u u u dt t y u u u u u u u u dt t z αθθωωαθθωωαθθωω??? = =++++-???? ????==++++-???? ????==++++-? ??? 各项含义： （1） 平移速度 （2）线变形运动所引起的速度增量
不可压缩流体动力学基础 1．已知平面流场的速度分布为xy x u x +=2，y xy u y 522+=。求在点（1，-1）处流体微团的线变形速度，角变 形速度和旋转角速度。 解：（1）线变形速度：y x x u x x +=??=2θ 54+=??=xy y u y y θ 角变形速度：()x y y u x u x y z +=??? ? ????+??=222121ε 旋转角速度：()x y x u x u x y z -=???? ????-??=2221 21ω 将点（1，-1）代入可得流体微团的1=x θ，1=y θ；23/z =ε；21/z =ω 2．已知有旋流动的速度场为322+=y u x ，x z u y 32+=，y x u z 32+=。试求旋转角速度，角变形速度和 涡线方程。 解：旋转角速度：21 21=???? ????-??=z u y u y z x ω 2 121=??? ????-??=x u z u z x y ω 2121=???? ????-??=y u x u x y z ω 角变形速度：2 521=???? ????+??=z u y u y z x ε 2 521=??? ????-??=x u z u z x y ε 2521=??? ? ????-??=y u x u x y z ε 由z y x dz dy dx ωωω==积分得涡线的方程为： 1c x y +=，2c x z +=
3．已知有旋流动的速度场为22z y c u x +=，0=y u ，0=z u ，式中c 为常数，试求流场的涡量及涡线方程。 解：流场的涡量为： 0=??-??=z u y u y z x Ω 22z y cz x u z u z x y +=??-??=Ω 22z y cy y u x u x y z +-=??-??=Ω 旋转角速度分别为：0=x ω 222z y cz y +=ω 222z y cy z +-=ω 则涡线的方程为：c dz dy z y +=??ωω 即c y dz z dy +-=?? 可得涡线的方程为： c c y =+22 4．求沿封闭曲线2 22b y x =+,0=z 的速度环量。（1）Ax u x =，0=y u ；（2）Ay u x =，0=y u ；（3）0=y u ，r A u =θ。其中A 为常数。 解：（1）由封闭曲线方程可知该曲线时在z =0的平面上的圆周线。 在z =0的平面上速度分布为： Ax u x =，0=y u 涡量分布为：0=z Ω 根据斯托克斯定理得：0==?z A z s dA ΩΓ （2）涡量分布为：A z -=Ω 根据斯托克斯定理得：2b A dA z A z s πΩΓ-==?
CFD 流体动力学软件介绍 CFD—计算流体动力学，因历史原因，国一直称之为计算流体力学。其结构为： 提出问题—流动性质（流、外流；层流、湍流；单相流、多相流；可压、不可压等等），流体属性（牛顿流体：液体、单组分气体、多组分气体、化学反应气体；非牛顿流体） 分析问题—建模—N-S方程（连续性假设），Boltzmann方程（稀薄气体流动），各类本构方程与封闭模型。 解决问题—差分格式的构造/选择，程序的具体编写/软件的选用，后处理的完成。 成果说明—形成文字，提交报告，赚取应得的回报。 CFD实现过程： 1.建模——物理空间到计算空间的映射。 主要软件： 二维： AutoCAD： 大家不要小看它，非常有用。一般的网格生成软件建模都是它这个思路，很少有参数化建模的。相比之下AutoCAD的优点在于精度高，草图处理灵活。可以这样说，任何一个网格生成软件自带的建模工具都是非参数化的，而对于非参数化建模来说，AutoCAD应该说是最好的，毕竟它发展了很多很多年！ 三维： CATIA：航空航天界CAD的老大，法国人的东西，NB，实体建模厉害，曲面建模独步武林。本身可以生成有限元网格，前几天又发布了支持ICEM-CFD的插件ICEM-CFD CAA V5。有了它和ICEM-CFD，可以做任何建模与网格划分！ UG：总觉得EDS脑袋进水了，收了I-deas这么久了，也才发布个几百M的UG NX 2.0，还被大家争论来争论去说它如何的不好用！其实，软件本身不错，大公司用得也多，可是就这么打市场，早晚是走下坡路。按CAD建模的功能来说它排不上第一，也不能屈居第二，尤其是加上了I－DEAS更是如虎添翼。现
思考题及答案 一、选择 (1) 二、例题 (2) 三、问答 (14) 一、选择 问题：恒定流是： A、流动随时间按一定规律变化； B、流场中任意空间点的运动要素不随时间变化； C、各过流断面的速度分布相同； D、各过流断面的压强相同。 问题：非恒定流是： A、； B、； C、； D、。 问题：一元流动是： A、均匀流； B、速度分布按直线变化； C、运动参数是一个空间坐标和时间变量的函数； D、限于直线流动。 问题：均匀流是： A、当地加速度为零； B、迁移加速度为零； C、向心加速度为零； D、合加速度为零。 问题1：流速势函数存在的必要与充分条件是: A、平面无旋流动； B、理想流体平面流动； C、不可压缩流体平面流动； D、无旋流动。 问题2：设流速势函数j=xyz，则点B（1，2，1）处的速度u 为： B A、5； B、1； C、3； D、2。
判断：公式（3－14）与公式（3－16）两式形式完全相同，因此其应用条件也相同。 你的回答：对错 判断：土坝渗流中的流网网格一定是直线正方形网格。 你的回答：对错 二、例题 例1如图3-7，已知流速场为，其中C为常数，求流 线方程。 解：由式得 图3-7 积分得： 则： 此外，由得： 因此，流线为Oxy平面上的一簇通过原点的直线，这种流动称为平面点源流动（C＞0时）或平
面点汇流动（C＜0时） 例2已知平面流动 试求：（1）t=0时，过点M（-1，-1）的流线。 （2）求在t=0时刻位于x=-1，y=-1点处流体质点的迹线。解：（1）由式 （2）由式 得 得 得： 由t=0时，x=-1，y=-1得C 1=0, C 2 =0，则有： 将：t=0，x=-1，y=-1 代入得瞬时流线 xy=1 最后可得迹线为： 即流线是双曲线。 例3已知流动速度场为
课后习题 第一章 1.计算流体动力学的基本任务是什么 计算流体动力学是通过计算机数值计算和图像显示，对包含有流体流动和热传导等相关物理现象的系统所做的分析。 2.什么叫控制方程？常用的控制方程有哪几个？各用在什么场合？ 流体流动要受物理守恒定律的支配，基本的守恒定律包括：质量守恒定律、动量守恒定律、能量守恒定律。如果流动包含有不同组分的混合或相互作用，系统还要遵守组分守恒定律。如果流动处于湍流状态，系统还要遵守附加的湍流输运方程。控制方程是这些守恒定律的数学描述。 常用的控制方程有质量守恒方程、动量守恒方程、能量守恒方程、组分质量守恒方程。质量守恒方程和动量守恒方程任何流动问题都必须满足，能量守恒定律是包含有热交换的流动系统必须满足的基本定律。组分质量守恒方程，在一个特定的系统中，可能存在质的交换，或者存在多种化学组分，每种组分都需要遵守组分质量守恒定律。 4.研究控制方程通用形式的意义何在？请分析控制方程通用形式中各项的意义。 建立控制方程通用形式是为了便于对各控制方程进行分析，并用同一程序对各控制方程进行求解。
各项依次为瞬态项、对流项、扩散项、源项。 6.CFD商用软件与用户自行设计的CFD程序相比，各有何优势？常用的商用CFD软件有哪些？特点如何？ 由于CFD的复杂性及计算机软硬件条件的多样性，用户各自的应用程序往往缺乏通用性。 CFD商用软件的特点是 功能比较全面、适用性强。 具有比较易用的前后处理系统和其他CAD及CFD软件的接口能力，便于用户快速完成造型、网格划分等工作。 具有比较完备的容错机制和操作界面，稳定性高。 可在多种计算机、多种操作系统，包括并行环境下运行。 常用的商用CFD软件有PHOENICS、CFX、SRAR-CD、FIDAP、FLUENT。PHOENICS除了通用CFD软件应该拥有的功能外，PHOENICS软件有自己独特的功能：开放性、CAD接口、运动物体功能、多种模型选择、双重算法选择、多模块选择。 CFX除了可以使用有限体积法外，还采用基于有限元的有限体积法。用于模拟流体流动、传热、多相流、化学反应、燃烧问题。其优势在于处理流动物理现象简单而几何形状复杂的问题。 SRAR-CD基于有限体积法，适用于不可压流体和可压流的计算、热力学的计算及非牛顿流的计算。它具有前处理器、求解器、后处理器三大模块，以良好的可视化用户界面把建模、求解及后处理与全部的物理模型和算法结合在一个软件包中。
第五章 实际流体动力学基础 5—1设在流场中的速度分布为u x =2ax ，u y =-2ay ，a 为实数，且a >0。试求切应力τxy 、τyx 和附加压应力p ′x 、p ′y 以及压应力p x 、p y 。 解：0y x xy yx u u x y ττμ??? ?==+= ????? 24x x u p a x μμ?'=-=-?，24y y u p a y μμ?'=-=?， 4x x p p p p a μ'=+=-，4y y p p p p a μ'=+=+ 5－2 设例５－1中的下平板固定不动，上平板以速度v 沿x 轴方向作等速运动（如图所示），由于上平板运动而引起的这种流动，称柯埃梯（Couette ）流动。试求在这种流动情况下，两平板间的速度分布。（请将 d 0d p x =时的这一流动与在第一章中讨 论流体粘性时的流动相比较） 解：将坐标系ox 轴移至下平板，则边界条件为 y ＝０，0X u u ==；y h =，u v =。 由例５－1中的（11）式可得 2d (1)2d h y p y y u v h x h h μ=- - （１） 当d 0d p x =时，y u v h =，速度ｕ为直线分布，这种特殊情况的流动称简单柯埃梯流动或简单剪切流动。它只是由于平板运动，由于流体的粘滞性带动流体发生的流动。 当d 0d p x ≠时，即为一般的柯埃梯流动，它是由简单柯埃梯流动和泊萧叶流动叠加而成，速度分布为 (1)u y y y p v h h h =-- （２） 式中2d ()2d h p p v x μ= - （３） 当p ＞０时，沿着流动方向压强减小，速度在整个断面上的分布均为正值；当p ＜０时，沿流动方向压强增加，则可能在静止壁面附近产生倒流，这主要发生p ＜－１的情况． 5－3 设明渠二维均匀（层流）流动，如图所示。若忽略空气阻力，试用纳维—斯托克斯方程和连 续性方程，证明过流断面上的速度分布为2 sin (2)2x g u zh z r q m = -，单宽流量3sin 3gh q r q m =。 解：（1）因是恒定 二维流动， 0y x z u u u t t t ???===抖?，u u x =，0y u =， 0z u =，由纳维——斯托克 斯方程和连续性方程可 得
非流体、热动专业CFD新手入门 首先掌握流体力学基本原理，丁祖荣主编的流体力学这本教材，仔细看两天，这样就会知道gambit中为什么会有边界层设置，边界层厚度如何设置；雷诺系数如何确定来判断层流与湍流；马赫数如何确定来判断流体是可压还是不可压，这样就能解决Fluent，是基于压力还是基于密度求解。能够对实际中一些看似简单的流体现象有深刻的认识，能够准确判断是定常流还是非定常流。 CFD网格划分 网格划分对于初学者所接触案例，其实非常简单。但实际工程中，大项目，特别涉及到整套工程，如环保，飞机，网格质量与数量都要求非常高，往往服务器类的PC才能解决问题，所谓的内存128G，CPU四核主频3.0以上。初学者，简单的管道，一般的机器还是没问题。有机械三维软件基础的，对于gambit建模就非常容易了。往往大项目，复杂的结构gambit 建模显得力不从心，所以对于流体工作者来说，学习三维软件对于建模有莫大的帮助，如Proe。 1.1Gambit介绍 网格的划分使用Gambit软件，首先要启动Gambit，在Dos下输入Gambit
第三章 流体动力学基础 习 题 一、单选题 1、在稳定流动中，在任一点处速度矢量是恒定不变的，那么流体质点是 （ ） A ．加速运动 B ．减速运动 C ．匀速运动 D ．不能确定 2、血管中血液流动的流量受血管内径影响很大。如果血管内径减少一半，其血液的流量将变为原来的（ ）倍。 A ．21 B ．41 C ．81 D ．161 3、人在静息状态时，整个心动周期内主动脉血流平均速度为0.2 m/s ，其内径d =2×10-2 m ，已知血液的粘度η =×10-3 Pa·S，密度ρ=×103 kg/m 3 ，则此时主动脉中血液的流动形态处于（ ）状态。 A ．层流 B ．湍流 C ．层流或湍流 D ．无法确定 4、正常情况下，人的小动脉半径约为3mm ，血液的平均速度为20cm/s ，若小动脉某部分被一硬斑阻塞使之变窄，半径变为2mm ，则此段的平均流速为（ ）m/s 。 A ．30 B ．40 C ．45 D ．60 5、有水在同一水平管道中流动，已知A 处的横截面积为S A =10cm 2 ，B 处的横截面积为 S B =5cm 2，A 、B 两点压强差为1500Pa ，则A 处的流速为（ ）。 A ．1m/s B ．2m/s C ．3 m/s D ．4 m/s 6、有水在一水平管道中流动，已知A 处的横截面积为S A =10cm 2 ，B 处的横截面积为S B =5cm 2 ，A 、B 两点压强之差为1500Pa ，则管道中的体积流量为（ ）。 A ．1×10-3 m 3 /s B ．2×10-3 m 3 /s C ．1×10-4 m 3 /s D ．2×10-4 m 3 /s 7、通常情况下，人的小动脉内径约为6mm ，血流的平均流速为20cm/s ，若小动脉某处被一硬斑阻塞而变窄，测得此处血流的平均流速为80cm/s ，则小动脉此处的内径应为（ ）mm 。 A ．4 B ．3 C ．2 D ．1 8、正常情况下，人的血液密度为×103 kg/m 3 ，血液在内径为6mm 的小动脉中流动的平均速度为20cm/s ，若小动脉某处被一硬斑阻塞而变窄，此处内径为4mm ，则小动脉宽处与窄处压强之差（ ）Pa 。 二、判断题
借宝地写几个小短文，介绍CFD的一些实际的入门知识。主要是因为这里支持Latex，写起来比较便。 CFD，计算流体力学，是一个挺难的学科，涉及流体力学、数值分析和计算机算法，还有计算机图形学的一些知识。尤其是有关偏微分程数值分析的东西，不是那么容易入门。大多数图书，片中数学原理而不重实际动手，因为作者都把读者当做已经掌握基础知识的科班学生了。所以数学基础不那么好的读者往往看得很吃力，看了还不知道怎么实现。本人当年虽说是学航天工程的，但是那时本科教育已经退步，基础的流体力学课被砍得只剩下一维气体动力学了，因此自学CFD的时候也是头晕眼花。不知道怎么实现，也很难找到教学代码——那时候网络还不发达，只在教研室的故纸堆里搜罗到一些完全没有注释，编程风格也不好的冗长代码，硬着头皮分析。后来网上淘到一些代码研读，结合书籍论文才慢慢入门。可以说中间没有老师教，后来赌博士为了混学分上过CFD专门课程，不过那时候我已经都掌握课堂上那些了。 回想自己入门艰辛，不免有一个想法——写点通俗易懂的CFD入门短文给师弟师妹们。本人不打算搞得很系统，而是希望能结合实际，阐明一些最基本的概念和手段，其中一些复杂的道理只是点到为止。目前也没有具体的计划，想到哪里写到哪里，因此可能会很零散。但是我争取让初学CFD 的人能够了解一些基本的东西，看过之后，会知道一个CFD代码怎么炼成的（这“炼”字好像很流行啊）。欢迎大家提出意见，这样我尽可能的可以追加一些修改和解释。
言归正传，第一部分，我打算介绍一个最基本的算例，一维激波管问题。说白了就是一根两端封闭的管子，中间有个隔板，隔板左边和右边的气体状态（密度、速度、压力）不一样，突然把隔板抽去，管子面的气体怎么运动。这是个一维问题，被称作黎曼间断问题，好像是黎曼最初研究双曲微分程的时候提出的一个问题，用一维无粘可压缩Euler程就可以描述了。 这里 这个程就是描述的气体密度、动量和能量随时间的变化（）与它们各自的流量（密度流量，动量流量，能量流量）随空间变化（）的关系。 在CFD常把这个程写成矢量形式 这里 进一步可以写成散度形式
第八章管道不可压缩流体恒定流 有压管流是日常生活中最常见的输水方式，本章主要介绍了有压管流的水力特点，计算问题以及简单管道与串联、并联和管网的水力计算原理与应用。 概述 一、概念 有压管流（penstock）：管道中流体在压力差作用下的流动称为有压管流。 有压恒定管流：管流的所有运动要素均不随时间变化的有压管流。 有压非恒定管流：管流的运动要素随时间变化的有压管流。 观看录像 二、分类 1.有压管道根据布置的不同，可分为： 简单管路：是指管径、流速、流量沿程不变，且无分支的单线管道。 复杂管路：是指由两根以上管道所组成的管路系统。 2.按局部水头损失和流速水头之和在总水头损失中所占的比重，管道可分为 长管：指管道中以沿程水头损失为主，局部水头损失和流速水头所占比重小于（5%-10%）的沿程水头损失，从而可予以忽略的管道。 短管：局部水头损失和流速水头不能忽略的、需要同时计算的管道。 三、有压管道水力计算的主要问题 1.验算管道的输水能力：在给定作用水头、管线布置和断面尺寸的情况下，确定输送的流量。 2.确定水头：已知管线布置和必需输送的流量，确定相应的水头。 3.绘制测压管水头线和总水头线：确定了流量、作用水头和断面尺寸（或管线）后，计算沿管线各断面的压强、总比能，即绘制沿管线的测压管水头线和总水头线。 第一节简单管道的水力计算 一、基本公式 1.淹没出流 图8-1中，列断面１－１与２－２的能量方程(4-15)，
图8-1 令： 且w1>>w, w2>>w，则有 （8-1） 说明：简单管道在淹没出流的情况下，其作用水头H0完全被消耗于克服管道由于沿程阻力、局部阻力所作负功所产生的水头损失上。即： 管道中的流速与流量为： （8-2） （8-3） 式中： ——管系流量系数，，它反映了沿程阻力和局部阻力对管道输水能力的影响。H0——作用水头，指上、下游水位差加上游行进流速的流速水头。 ——局部阻力系数，包含出口损失。 问题：图示两根完全相同的长管道，只是安装高度不同，两管道的流量关系为：
不可压缩流体动力学基础 1.已知平面流场的速度分布为xy x u x +=2,y xy u y 522+=。求在点(1,-1)处流体微团的线变形速度,角变形速度与 旋转角速度。 解:(1)线变形速度:y x x u x x +=??=2θ 54+=??=xy y u y y θ 角变形速度:()x y y u x u x y z +=??? ? ????+??=222121ε 旋转角速度:()x y x u x u x y z -=???? ????-??=2221 21ω 将点(1,-1)代入可得流体微团的1=x θ,1=y θ;23/z =ε;21/z =ω 2.已知有旋流动的速度场为322+=y u x ,x z u y 32+=,y x u z 32+=。试求旋转角速度,角变形速度与涡线方程。 解:旋转角速度:21 21=???? ????-??=z u y u y z x ω 2 121=??? ????-??=x u z u z x y ω 2121=???? ????-??=y u x u x y z ω 角变形速度:2 521=???? ????+??=z u y u y z x ε 2 521=??? ????-??=x u z u z x y ε 2521=??? ? ????-??=y u x u x y z ε 由z y x dz dy dx ωωω==积分得涡线的方程为: 1c x y +=,2c x z +=
3.已知有旋流动的速度场为22z y c u x +=,0=y u ,0=z u ,式中c 为常数,试求流场的涡量及涡线方程。 解:流场的涡量为: 0=??-?? =z u y u y z x Ω 22z y cz x u z u z x y +=??-??=Ω 22z y cy y u x u x y z +-=??-??=Ω 旋转角速度分别为:0=x ω 2 22z y cz y +=ω 222z y cy z +-=ω 则涡线的方程为:c dz dy z y += ??ωω 即c y dz z dy +-=?? 可得涡线的方程为:c c y =+22 4.求沿封闭曲线2 22b y x =+,0=z 的速 度环量。(1)Ax u x =,0=y u ;(2)Ay u x =,0=y u ;(3)0=y u ,r A u =θ。其中A 为常数。 解:(1)由封闭曲线方程可知该曲线时在z =0的平面上的圆周线。 在z =0的平面上速度分布为: Ax u x =,0=y u 涡量分布为:0=z Ω 根据斯托克斯定理得:0==?z A z s dA ΩΓ (2)涡量分布为:A z -=Ω 根据斯托克斯定理得:2 b A dA z A z s πΩΓ-==?
1738年瑞士数学家：伯努利在名著《流体动力学》中提出了伯努利方程。 1755年欧拉在名著《流体运动的一般原理》中提出理想流体概念，并建立了理想流体基本方程和连续方程，从而提出了流体运动的解析方法，同时提出了速度势的概念。 1781年拉格朗日首先引进了流函数的概念。 1826年法国工程师纳维，1845年英国数学家、物理学家斯托克思提出了著名的N-S方程。 1876年雷诺发现了流体流动的两种流态：层流和紊流。 1858年亥姆霍兹指出了理想流体中旋涡的许多基本性质及旋涡运动理论，并于1887年提出了脱体绕流理论。 19世纪末，相似理论提出，实验和理论分析相结合。 1904年普朗特提出了边界层理论。 20世纪60年代以后，计算流体力学得到了迅速的发展。流体力学内涵不断地得到了充实与提高。 理想势流伯努利方程 （3-14） 或（3-15） 物理意义:在同一恒定不可压缩流体重力势流中，理想流体各点的总比能相等即在整个势流场中，伯努利常数C 均相等。 （应用条件：“”所示） 符号说明 二、沿流线的积分
1.只有重力作用的不可压缩恒定流，有 2.恒定流中流线与迹线重合： 沿流线（或元流）的能量方程： （3-16） 注意：积分常数C，在非粘性、不可压缩恒定流流动中，沿同一流线保持不变。一般不同流线各不相同（有旋流）。 （应用条件：“”所示，可以是有旋流） 流速势函数（势函数）观看录像>> ?存在条件：不可压缩无旋流，即或 必要条件存在全微分d 直角坐标 （3-19） 式中：——无旋运动的流速势函数，简称势函数。 ?势函数的拉普拉斯方程形式 对于不可压缩的平面流体流动中，将（3-19）式代入连续性微分方程（3-18），有： 或（3-20） 适用条件：不可压缩流体的有势流动。 点击这里练习一下 极坐标 （3-21） 流函数
第二章计算流体力学的基本知识 流体流动现象大量存在于自然界及多种工程领域中，所有这些工程都受质量守恒、动量守恒和能量守恒等基本物理定律的支配。这章将首先介绍流体动力学的发展和流体力学中几个重要守恒定律及其数学表达式，最后介绍几种常用的商业软件。 2.1 计算流体力学简介 2.1.1计算流体力学的发展 流体力学的基本方程组非常复杂，在考虑粘性作用时更是如此，如果不靠计算机，就只能对比较简单的情形或简化后的欧拉方程或N-S方程进行计算。20世纪30～40年代，对于复杂而又特别重要的流体力学问题，曾组织过人力用几个月甚至几年的时间做数值计算，比如圆锥做超声速飞行时周围的无粘流场就从1943年一直算到1947年。 数学的发展，计算机的不断进步，以及流体力学各种计算方法的发明，使许多原来无法用理论分析求解的复杂流体力学问题有了求得数值解的可能性，这又促进了流体力学计算方法的发展，并形成了"计算流体力学"。 从20世纪60年代起，在飞行器和其他涉及流体运动的课题中，经常采用电子计算机做数值模拟，这可以和物理实验相辅相成。数值模拟和实验模拟相互配合，使科学技术的研究和工程设计的速度加快，并节省开支。数值计算方法最近发展很快，其重要性与日俱增。 自然界存在着大量复杂的流动现象，随着人类认识的深入，人们开始利用流动规律来改造自然界。最典型的例子是人类利用空气对运动中的机翼产生升力的机理发明了飞机。航空技术的发展强烈推动了流体力学的迅速发展。 流体运动的规律由一组控制方程描述。计算机没有发明前，流体力学家们在对方程经过大量简化后能够得到一些线形问题解析解。但实际的流动问题大都是复杂的强非线形问题，无法求得精确的解析解。计算机的出现以及计算技术的迅速发展使人们直接求解控制方程组的梦想逐步得到实现，从而催生了计算流体力
第三章流体动力学基础 描述流体运动的两种方法： 拉格朗日法和欧拉法。除个别质点的运动问题外，都应用欧拉法。 拉格朗日法：是以个别质点为研究对象，观察该质点在空间的运动，然后将每个质点的运动情况汇总，得到整个流体的运动。质点的运动参数是起始坐标和时间变量t的连续函数。 欧拉法：是以整个流动空间为研究对象，观察不同时刻各空间点上流体质点的运动，然后将每个时刻的情况汇总起来，描述整个运动。空间点的物理量是空间坐标）和时间变量t的连续函数。 恒定流：各空间点上的运动参数都不随时间变化的流动。 非恒定流：各空间点上的运动参数随时间变化的流动。 一（二、三）元流：流体流动时各空间点上的运动参数是一（二、三）个空间坐标和时间变量的连续函数。 均匀流：流线是平行直线的流动。 非均匀流：流线不是平行直线的流动。 流线：表示某时刻流动方向的曲线，曲线上各质点的速度矢量都与该曲线相切。迹线：流体质点在一段时间内的运动轨迹。 流管：某时刻，在流场内任意做一封闭曲线，过曲线上各点做流线，所构成的管状曲面。 流束：充满流体的流管。 过流断面：与所有流线正交的横断面。 元流：过流断面无限小的流束，断面上各点的运动参数均相同。
总流：过流断面为有限大小的流束，断面上各点的运动参数不相同。流量：单位时间内通过某一过流断面的流体量。以体积计为体积流量，简称流量；以质量计为质量流量；以重量计为重量流量 非均匀渐变流：在非均匀流中流线近似于平行直线的流动。 水头线：总流或元流沿程能量变化的几何图示。 水力坡度：单位流程内的水头损失。 （简答）流线有哪些主要性质？流线和迹线有无重合的情况？答：流线性质：（1）在恒定流中，流线的形状和位置不随时间变化；（2）在同一时刻，一般情况下流线不能相交或转折。在恒定流中流线与迹线重合，非恒定流中一般情况下两者不重合，但当速度方向不随时间变化只是速度大小随时间变化时，两者仍重合。 试述流动分类：（1）根据运动参数是否随时间变化，分为恒定流和非恒定流；（2）根据运动参数与空间坐标的关系，分为一元流、二元流和三元流；（3）根据流线是否平行，分为均匀流和非均匀流。 不可压缩流体的连续性微分方程：不可压缩流体运动必须满足该方程。