流体力学知识点大全

1、流体力学-笔记参考书籍：全美经典-流体动力学流体力学张兆顺、崔桂香流体力学吴望一一维不定常流流体力学课件清华大学王亮主讲目录：第一章绪论第二章流体静力学第三章流体运动的数学模型第四章量纲分析和相似性第五章粘性流体和边界层流动第六章不可压缩势流第七章一维可压缩流动第八章二维可压缩流动气体动力学第九章不可压缩湍流流动第十章高超声速边界层流动第十一章磁流体动力学第十二章非牛顿流体第十三章波动和稳定性第一章绪论1、牛顿流体：剪应力和速度梯度之间的关系式称为牛顿关系式，遵守牛顿关系式的流体是牛顿流体。2、理想流体：无粘流体，流体切应力为零，并且没有湍流？。此时，流体内部没有内摩擦，也就没有内耗散和损失。

2、层流：纯粘性流体，流体分层，流速比较小；湍流：随着流速增加，流线摆动，称过渡流，流速再增加，出现漩涡，混合。因为流速增加导致层流出现不稳定性。定常流：在空间的任何点，流动中的速度分量和热力学参量都不随时间改变，3、欧拉描述：空间点的坐标；拉格朗日：质点的坐标；4、流体的粘性引起剪切力，进而导致耗散。5、无黏流体无摩擦流动不分离无尾迹。6、流体的特性：连续性、易流动性、压缩性不可压缩流体：=0DtP=const是针对流体中的同一质点在不同时刻保持不变，即不可压缩流体的密度在任何时刻都保持不变。是一个过程方程。7、流体的几种线流线：是速度场的向量线，是指在欧拉速度场的描述；同一时刻、不同质点连接起

3、来的速度场向量线；drPU（x,t）ndrxU=0迹线：流体质点的运动轨迹，是流体质点运动的几何描述；同一质点在不同时刻的位移曲线；涡线：涡量场的向量线，舀=VxU,drP（n（x,t）ndrx=0涡线的切线和当地的涡量或准刚体角速度重合，所以，涡线是流体微团准刚体转动方向的连线，形象的说：涡线像一根柔性轴把微团穿在一起。第二章流体静力学AFdF1、压强：p=lim=AAt0AAdA静止流场中一点的应力状态只有压力。2、流体的平衡状态：1）、流体的每个质点都处于静止状态，=整个系统无加速度；2）、质点相互之间都没有相对运动，=整个系统都可以有加速度；由于流体质点之间都没有相对运动，导致剪应力处

4、处为零，故只有：体积力（重力、磁场力）和表面力（压强和剪切力）存在。3、表面张力：两种不可混合的流体之间的分界面是曲面，则在曲面两边存在一个压强差。4、正压流场：流体中的密度只是压力（压强）的单值函数。Pp5、涡量不生不灭定理拉格朗日定理：理想正压流体在势力场中运动时，如某一时刻连续流场无旋，则rrrr流场始终无旋。JndA=0,=VxU,有斯托克斯公式得：r=JUx=JndA=0,l0拉格朗日定理是判断理想正压流体在势力场中运动是否无旋的理论依据。涡量的产生原因：（A）流体的粘性；非理想流体；（B）非正压流体；大气和海洋中的密度分层（非正压）导致漩涡；（C）非有势力场；气流科氏力（非有势力）

5、作用导致漩涡；（D）流场的间断，高速气流中的曲面激波后，产生有旋流流场；第三章流体运动的数学模型1、积分型的流体方程a）、质量守恒定律：物理意义：流出控制体表面的净质量流量等于控制体内质量对时间的减少率。rraJpV-dA=-JpdTatC.SC.Vb）、动量守恒：牛顿第二定律F(表面力)+JB-dT(体积力)sC.Vrarrrr=F=JpVdT+JpVV-dAatC.VC.SC）、角动量JrXdF+JrXB-dT=J?x(pVr)dT+J?x(pVvdAsatC.SC.VC.VC.S每一项物理意义：JrrxdFr：控制面上的力对原点的力矩sC.SrJ；xBdT:体积力对原点的力矩，C.V0j

6、rxCPdT:质量元的角动量，控制体内流体的总角动量,atC.VjrxCVv）dA:通过控制面的角动量流出率C.Sd）、能量守恒（热力学第一定律）Q-W=AEdQ-dWdtdt=JpedT+J(e+p/p)pVdAatC.VC.SDJpEdV=D*(t)rrrpfUdV+JJjTnrUdA+D*(t)Z*(t)JpqdV+JJ/VTdAD*(t)Z*(t)质量体内的总能量增长率：DfpEdV,DtD*t)rrTUdAn边界面上因热传导输入的热量：0九n-VTdAz*(t)S是系统的熵体积力所作的功率：fpfr-UrdV；表面力所作的功率：0ffD*(t)质量体内的生成热：fpq&dVe)、热力

7、学第二定律D*(t)dST-0,2、有积分形式到微分形势的方程，有三种方法、应用矢量的微积分；、积分应用于体积元，有体积元趋于零，取极限推得；、将系统的方程直接应用体积元，再将积分表达式取极限欧拉坐标，即：笛卡尔坐标，V二V(r，t)=V(x,y,z,t)；拉格朗日，刚体描述，速度、加速度分别为：r,r3、微分型的流体方程1)、连续性方程：单位时间流入控制体的质量等于控制体内质量的增加dp-dt+V-(pVr)=0定常流dpidt=0nV-(pVr)=0r不可压缩：DpfDt=0nV-V=0一维定常流：pAV=pAV1112222)、动量方程：单位时间流入控制体的动量以及作用于控制体上的外力之

8、和，等于控制体动量的增加。应力张量：代表剪应力和正应力；应力张量一定是对称的；否则，当体积元收缩成无限小时，必将以无限大的角速度旋转。因此，应力张量只能有六个分量。局部加速度：非定常流动，对流加速度：面积的变化；欧拉坐标系和拉格朗日中的速度和加速度其大小和方向都不会改变rDV/Dt欧拉=&拉格朗日rr涡量：速度矢量的旋度，=VxV角速度：占=-r=-VxVQ=0无旋流动22dtrt=pB+V-Fpb:体积力，F面积力;、能量方程：单位时间流入流体的能量、外界传入的热量、外力做功的总和等于控制体内能量的增加。V-q+pqR热传导非传导热Q(pE)(rrrftr+V-VpEV)=pB-V+V-P-

9、V-dt增加量流入量体积力做功表面力做功E二e+1V2,2V-q=V（九VT）,q=XVT,Fourier热传导定律九：热传导系数，pq：非热传导热，即：热辐射、化学生成热，R几种特殊情况：、定常流体：=0;dt、绝热过程：V-q=q=0，没有外界热传入Rr、质量力有势：B=-VG；(4)、理想流体：trrP=p=-np。n本构方程：求解方程组，流体微团的应力状态和微团变形运动状态间的物性关系式本构方程是张量方程；使得控制方程得以封闭，可以求解方程；控制方程+热力学状态方程+本构方程边界条件：固体壁面的不可穿透条件；垂直于壁面的法向速度连续(U-n)=(U-n),XbXrrU为固壁的速度，U为

10、同一点的流体质点的速度b无穷远条件无穷远处，流体保持静止状态；|x|T8,U=0,p=p,p=p,绕流条件参考系固结在运动物体上，无穷远处的来流条件：x|T8,U=U,p=p,P=P,ggg4、求解物理问题的基本步骤：1）、特定的物理问题；2）、物理模型描述；3）、数学模型的建立4）、求解数学方程；5）、实验验证结果；5、理想流体动力学无粘性，亦即无热传导，压力分布rDVDt欧拉方程：纳维-斯托克斯方程兰姆（Lamb）方程：rQVrrr+V-VV=f-VpTOC o 1-5 h zQtPrrDVQVrrr1rrrV-VV=VrV2rrVX。,r。=VXV=0,V2、rrr1QVV2、VxQ=f

11、一Vp,+Vl2)PQtV2)2丿=f-PVp,Pr空+VQt将欧拉方程中的对流导数项换成旋量形式，即是Lamb型方程6、速度势r。=0,=一+V-VV=f-Vp+AU,不可压、粘性流DtQtPrr因为无旋，故有速度势存在；VxU=0,nU=VQ,静止不可压缩理想流体在瞬时脉冲压强作用下产生的流动是无旋的，它的速度势等于负压强冲量除以密度；通过欧拉方程，在短时间内进行积分处理，得出：TOC o 1-5 h z1$tVIIU=Jpdt=,=+C,PPP00物理意义：不可压缩流体的无旋流动可由瞬时压强的冲量产生。7、流函数rQuQv在不可压缩流体的二维运动中，V-V=Qu+QV=0,Q屮Q屮u=,

12、v=一QyQxQxQy满足上式的全微分函数：d屮=udyvdx=0,流函数的定义式子：屮=J（udy-vdx）,流函数的等值线是流线;流函数等值线和势函数等值线是正交的。因为流函数的切线表示速度，而速度一定垂直于势函数，故，二者正交。8、复势以速度势为实部，流函数为虚部组成的复函数，W(z)=(x,y)+i(x,y),复速度：以平面无旋流场的速度分量组成的复数U=u+iv,9、理想不可压缩流体的有旋流动理想不可压缩流体在非有势力作用下将产生有旋流动；有旋流动的流函数：有旋流动无速度势，但不可压缩流体存在流函数:屮(x,y)d屮=udy一vdx=0,河河u=,v=dydx苦=VxU,.不可压缩流

13、体，U=0，不可压缩流体法向应力等于热力学压强;.可压缩流体，流体微团体积发生变化，引起压强P变化，称为“容积_m粘性系数”或”第二粘性系数”，因此，反应由体积变化引起流体偏离热力学压强的粘性应力。描述不可压缩、粘性流体的动量方程(运动学方程)称为：纳维斯托克斯方程rDVDtr3Vrrr1r一+V凡V=fVp+AU，dtp2、粘性流体运动的基本特性（1）、粘性流体运动的有旋性无粘流体满足Euler方程，满足边界的不可穿透条件；而无旋条件只能使得N-S方程满足粘性的部分条件，故粘性流体有旋；（2）、粘性流体运动的耗散性在不可压缩牛顿流体流动的能量方程中有一粘性耗散项，它使得流体质点的熵增加，即：

14、绝热系统中牛顿流体运动是熵增的不可逆耗散系统；（3）、粘性流体运动的扩散性r方程中的vAU具有扩散性质，使得具有有旋性的流体有旋区域不断扩大；3、流体绕物体流动区域：One:邻近物体表面的薄层（边界层），摩擦起主要作用；Tw0：另一区域摩擦可以忽略；当粘性流体绕流的特征雷诺数很大时（即：粘性很小时），在物体表面形成粘性起主导作用的薄层，即：边界层。普朗特提出边界层理论定常绕流中流体粘性只在贴近物面极薄的一层内主宰流体运动，称这一层为边界层；边界层外的流动可近似为无粘的理想流动。研究内容：A：边界层的厚度；B：导致的速度分布；C：压强的分布；D：流体作用的固体表面的力的方法；边界层内的流动开始是

15、层流，但沿物体表面边界层增厚，如果表面足够长，会出现一个转区，边界层内的流可以转变为湍流。4、边界层的流动与分离第一阶段：流动方向压强减小，称为顺压梯度区。此时如dx0,（62u/Qy2）0,（a2uay2）0.0流体流动过程中受到两个力的作用，一个是粘性力，一个是压强梯度力。在第一阶段，粘性力减速，而压强梯度力加速，即阻碍粘性力的减速。在第三阶段，粘性力和逆压强梯度力共同减速流体，甚至导致壁面附近的流体质点出现倒流。5、内流：考虑粘性的N-S方程流向（X轴）和横向（Y轴）的无量纲化转换：流向尺度x:L，横向尺度5:L08连续性方程、动量方程Qp+v-(pVr)=0,QtrQPV+V.QtDP

16、=0,Dt(rr)rIpVV丿=-Vp+yV2V,r空=0,dtQu*Qv*+一Qx*Qy*Qu*u*+vQx*=0,1|d2U*Re(dx*2dv*u*+vdx*(2)式中：Qu*Qp*=+Qy*Qx*Qv*1Qp*Qy*2Qy*Re2:O(1),Re?1,QuQv+=0,QxQyQuQuQpu+v=+yQxQyQxQvQvu+v=QxQy+2Qy*2/*1Q2v*+Qy*2丿1fQ2V+Re(Qx*22生+丫QyQ2uQ2u、+Qx2Qy2丿Q2vQ2v1:、Qx2Qy2丿Qu*Qv*+=0,Qx*Qy*Qu*u*+vQx*並=0,Qy*(1)Q2u*QQ卜厂，(2)Qy*Qx*Qy*2(

17、3)Qu*Qpu=Uu*,v=Uv*,x=Lx*,y=Ly*,p=pU2p*8(x*,y*)=p*(x*,8)=p*(x*,0),(3)式:1=0,中压强在法向为常数，即:Qy*有方程得出结论：沿壁面方向压强等于外部流场的当地.边界层内压强在垂直壁面方向不变壁面压强；.流向的分子粘性扩散远小于法向扩散(方程(2)中最后一项)；致使不可压缩流体定常流动的边界层方程有椭圆型的定常N-S方程退化为抛物型偏微分方程；V3.当Re?1,时，边界层横向尺度1/vRe,即：边界层的横向尺度与Re数的平方根成反比6、边界层厚度8(5)/、排挤厚度：5=f(1uU)dy,1e0物理意义：厚度为(55)的理想位势

18、流进入边界层后，由于近壁流速减小，它的外边界外移，相当于物面增加厚度5，故5称为位移厚度或排挤厚度;11(2)、动量损失厚度：5=f)(u；U)(1uU)dy,2ee边界层内流体的通量：fpu2dy,0流量相同的理想位势流的厚度等于8-5，其动量通量：pU2（5-5）1e1由于粘性，使流入边界层的动量通量和位势流相比损失量：pU2e（5-8）-fpu2dy,10已知U（5-5）=fudy,，故，动量通量损失为:e10e1ee0000则流过厚度5的动量通量：p5U2，22efpu（U-u）dy=p8e20U2,e=5=J一（U-u2Ue0e第六章不可压缩势流1、讨论不可压缩二维势流理论，适用于马

19、赫数小于0.3左右的亚声速流动。势流理论：无旋流动oV=-Vrw=VxV流体的旋度（或称涡量）：2、伯努利方程不可压缩、无旋流动、非定常的伯努利方程-dr二-fVp-dr-pfV屮-drpf空dr+pfvdtrr守恒体力的无旋流动中：B=VV，速度势：V=-Vr帥V2p小-+屮二Condt2prrIV23、速度势和流函数速度势：无旋流动ov=-ve不可压缩：V-V=0;V2=0任意二维流场，均可用来流函数屮表征。在二维流动中，等V线是流线，它在两流线之间的数值差等于该两流线之间的容积流率。流函数物理意义：由下图可知，沿从屮至屮的路径，流动从右到左为正向，笛12r卡尔坐标中以定义的V为：u=-d

20、屮6y;v=Q屮:;屮至屮之间的容积率为:V屮7刘7dx+dy=(Qxdy丿12Q=f(vdx-udy)=f12f即斗-屮21111从物理上讲：流函数是单值的。除沿任意包围奇点，如源或汇的封闭积分轮廓线外，沿任何封闭轮廓的积分f山=0。14、复函数在二维空间中，定义复函数必须速度势和流函数必须为调和函数且满足柯西黎曼方程，F=+“第七章一维可压缩流动(P160)一维非定常流见第八章二维可压缩流动气体动力学1、可压缩空气动力学流动问题：无摩擦、无旋和等熵的流动；在超声速流动中，可能会出现激波，激波中是不等熵的。绝热连续的流动过程是等熵过程；理想可压缩流动的方程组：连续、动量以及状态方程(与时间相

21、关时)p可以引入速度势的概念，进行化简求解，得到关于速度势的方程。2、在能量方程中：若流动是绝热(&=o)且连续的，即过程是绝热可逆的，有热力学第二定律：Ds=DqT，可导出熵增&=DsDt=&T=0，故：绝热连续的流动过程是等熵过程。一维声波的传播是非色散性的双向波，因为声速c=dp：dp),.7pp，由此可知声速只与热力学状态有关，与扰动的运动学特性，(扰动的频率、波长等无关)。马赫数M：流体的速度与当地的声速之比；物理解释：单位质量流体的惯性力与压强合力的量级之比；气体质点的单位质量的动能与内能的量级之比；马赫锥：在超声速绕流运动中产生的圆锥面角度；超声速运动的点扰动只能在下游马赫锥内传

22、播，而不能传播到马赫锥外。3、理想气体等熵流动的性质、理想气体定常绝热连续性流动中沿流线熵不变；、理想气体绝热定常流动沿流线h+U22=const；、克鲁克定理(Croco定理)QxU=TVs-Vh0有此公式可以判断：均熵、均焓及旋度之间的关系；当均熵、均焓时，流体无旋；当均熵、无旋时，流体均焓；当均焓、无旋时，流体均熵，等等滞止参数：在定常流动中，气体流动等熵地减速到速度等于零的状态，称为滞止状态，滞止状态的气流参数为滞止参数。滞止温度：因为等熵，故有能量方程：h=cT=h+U2.；2=cT+U2.2nT=T+竺,opopo2cp=pRT,pp“C,p=RTrfn=C,*p=(RT血(T、十

23、T0Y-1理想气体定常等熵流动中的最大速度：h0=cpT+T，WhenT=0*Umax临界参数：在理想气体定常等熵流动中，流体质点速度等于当地声速的状态称为临界状态，临界状态下的气体状态参数，称为临界参数。u*=a*=：yp*p*=：yRT*=：(y-1)h*,IU*2.；2+h*=h,i0速度系数：流体速度与临界速度之比*h=(y+1)h*；20c*化简：.UUccU九=0= HYPERLINK l bookmark164 o Current Document c*ccc*c0U2T=T+0T2c，*Tp(y-1)1+Ma22T0U2(y-1)U2=1+=1+2cT2yRTp-11+O竺=1

24、+OMa2,2c22U*2TU*2 HYPERLINK l bookmark190 o Current Document ,*一=1+2cT*2cT*ppr(y-1)=Ma1+Ma22T=T*+,*o=1+0*22RT*(y1)U*2(y+1)2c*21=MaTTT*04、激波理论在强扰动下，流动的参数发生突变的现象，称为激波；激波厚度约为分子自由程的量级，在这一薄层中，物理量迅速地从波前值变Y+11为波后值，速度梯度、压强梯度和密度梯度都很大，因此，研究激波层内流动时必须考虑粘性和热传导的作用。当激波层中不发生离解、电离等物理、化学过程时，气体穿过激波可认为是绝热过程。正激波：和气流速度垂直

25、的物理量间断面；驻激波：将坐标系固结在激波上，正激波可以看成是静止的平面；分析激波两侧的参数，考虑：连续性、动量、能量和状态方程面积分别为A,A,因为选取的控制体非常窄(dxT0)，故体积(dVT0)忽略.质量动量能量彳B!pdv+ffpUr-ndA=0,DA血pUdV+ffpUr0r)dA=-ffp-ndA,DfffpEdV+dtDffpE(U-n)dA=-ffpUdA,AAnpU=pU,1122np+pU2=p+pU2,11122211“ncT+U2=cT+U2,p121p222状态：1=2RTpRTp1122n1=2RTpRTp1122有第1,2公式可以得到：pU=pU,nU+U=U1+

26、v1122121p+pU2=p+pU2,111222、1+匕p2Jp-p=pU(U1211=pU11r11)+9IppJ12-U),21(p-p)12r11)+IppJ12=U2-U2,21再有第3式子，可以的11yp1ypcT+-U2=cT+-U2,n匕+U2=pp121p222y1p21y1p12r11)+9Ipp丿12、pp12y-1Ipp丿12同时乘以匕，整理后:p1+2U22,(y+1)匕(y1)pp2=1p1(y+1)(y1)匕p1(y+1)厶+(y1)p=p2一1-p1(y-1)昌+(y+1)p1有状态方程：上述关系式就是(y+1)仝+(y-1)(昱TT=pU!,T1(y+1)昱

27、+(y-1)p1兰金-于格尼奥(Rankine-Hugoniot)关系式有连续性方程、动量方程和能量方程推导出压强、密度和温度的比值。激波过程与等熵过程：.激波压缩是有限压缩，正激波后的密度增高有极限lim匕二离二6,p1C(、丫等熵压缩是无限的，lim2=lim邑=s,28P也8P丿p1p1v2.激波绝热曲线和等熵曲线在p.P1时相切，这说明，弱激波压缩接近等21熵压缩；.相同的密度比p.p1下，激波压缩过程的压强比大于等熵过程的压强比;21v4.激波压缩过程熵增必大于零，是绝热不可逆过程;As=s一s21=clnvP/PLIPfPY丿、11丿=clnvpp(P2P)1二cv激波压缩p/p1

28、时,，则有激波曲线和等熵曲线：p/p(p/p)，可知12121sAs=s一s0。21.激波膨胀是不可能的若有p.p1，激波后的压强小于激波前压强：pip1，于是:2121pp(p-p)，则出现As=s-s0，这是不可能发生的。2121s215、普朗特关系有动量方程除以连续方程，c2十+Uc22YU2应用临界参数的定义及动量方程：YpU2ypU2yp*U*2TOC o 1-5 h z1+4=2+亠=+一1P2Y一1P2Y一1P*212RTU2YRTU2YRT*U*21+4=2+彗=+,Y-12y-12Y-12(y-1)U2(y-1)U2C2+4=C2+21222(y+1)c*2(y-1)U2C2

29、=-4,122(y-1)U*2(y+1)c*2=C*2+=22(y+1)c*2(y-1)U2C2=-牛222将C2、C2的结果带入前面式子,12c*2c*2+U=+U,U1U212UU=c*2n12c2c2L+UYUiYU12（U-U）11_-12IUU丿12九九=1126、运动激波及其反射运动激波，选择激波作为相对坐标系0=0,Z）=Q-M7、斜激波理论：与气流方向不垂直的平面激波质量方程：fffpdV+ffpUr-ndA=0,dtDA动量方程：fjfpUdV+J!p击OhDdA=-J!p-ndA,dtDAAnp+pU2=p+pU2,U=U=U,111n21t2ttfffpEdV+JJpEA=-J!pUdA,nh+1U2=h+1U2,dt121n222nDAAPP1=2,RTpRTp1122能量方程：状态方程：激波压缩、等熵压缩对比P1U1n5、小扰动理论和线化理论6、特征线方法一维不定常流中介绍第九章不可压缩湍流流动1、由易到难的流动：位势流（流速很低）9层流（流速较低）9湍流（高）流场中存在无限小的扰动，当雷诺数很低时，扰动逐渐衰减，流动保持层流状态当雷诺数增大时，小扰动会逐渐增长，流动