描述流体运动的两种方法

1、描述流体运动的两种方法（姓名：张旺龙学号：3专业：流体力学）引言:描述流体运动的两种方法拉各朗日方法和欧拉方法设流体质点在空间中运动，我们的任务就是确定描写流体运动的方法并且将它用数学式子表达出来。在流体力学中描写运动的观点和方法有两种，即拉各朗日方法和欧拉方法。拉各朗日方法，着眼于流体质点。设法描述出每个流体质点自始至终的运动过程，即它们的位置随时间变化的规律。如果知道了所有流体质点的运动规律，那么整个流体运动的状况也就清楚了。欧拉方法的着眼点不是流体质点而是空间点。设法在空间中的每一点上描述出流体运动随时间的变化状况。如果，每一点的流体运动都已知道，则整个流体的运动状况也就清楚了。一拉格朗

2、日方法现在我们将上述描写运动的拉各朗日观点和方法用数学式子表达出来，为此首先必须用某种数学方法区别不同的流体质点。通常利用初始时刻流体质点的坐标作为区分不同流体质点的标志。设初始时刻t=to时，流体质点的坐标是（a,b,c）,它可以是曲线坐标，也可以是直角坐标（x0,y0,z0），重要的是给流体质点以标号而不在于采取什么具体的方式。我们约定采用a,b,c三个数的组合来区别流体质点，不同的a,b,c代表不同的质点。于是流体质点的运动规律数学上可表示为下列矢量形式：r=r（a,b,c,t）（1）其中r是流体质点的失径。在直角坐标系中，有x=x（a,b,c,t）y=y（a,b,c,t）z=z（a,b

3、,c,t）（2）变数a,b,c,t称为拉各朗日变数。在式（2）中，如果固定a,b,c而令t改变，则得某一流体质点的运动规律。如果固定时间t而令a,b,c改变，则得同一时刻不同流体质点的位置分布。应该指出，在拉各朗日观点中，失径函数的定义区域不是场，因为它不是空间坐标的函数，而是质点标号的函数。现在从（1）式出发来求流体质点的速度和加速度。假设由（1）式确定的函数具有二阶连续偏导数。速度和加速度是对于同一质点而言的单位时间内位移变化率及速度变化率，设v，v分别表示速度矢量和加速度矢量，则dr(a,b,c,t)v=vdt3)d2r(a,b,c,t)v=4)既然对同一质点而言，a,b,c不变，因此上

4、式写的是对时间t的偏导数。在直角坐标系中，速度和加速度的表达式是dx(a,b,c,t)u=Ldy(a,b,c,t)v=dtdz(a,b,c,t)w=dt5)d2x(a,b,c,t)u=d2y(a,b,c,t)v=62z(a,b,c,t)w=6)欧拉方法现在来介绍描写流体运动的另一种观点和方法，即欧拉方法。和拉各朗日方法不同，欧拉方法不同，欧拉方法的着眼点不是流体质点而是空间点。设法在空间中的每一点上描述出流体运动随时间的变化状况。如果，每一点的流体运动都已知道，则整个流体的运动状况也就清楚了，那么应该用什么物理量来表现空间点上流体运动的变化情况呢因为不同时刻将有不同流体质点经过空间某固定点，所

5、以站在固定点上就无法观测和记录掠过的流体质点以前和以后的详细历史。也就是说我们无法象拉各朗日方法那样直接测量出每个质点的位置随时间的变化情况。虽然如此，不同时刻经过固定点的流体质点的速度是可以测出的，这样采用速度矢量来描写固定点上流体运动的变化状况就是十分自然的了。考虑到上面所说的情形，欧拉方法中流体质点的运动规律数学上可表示为下列矢量形式：7)v=v(r,t)在直角坐标系中有：u=u(x,y,z,t)v=v(x,y,z,t)w=w(x,y,z,t)8)要完全描述运动流体的状况还需要给定状态函数压力、密度、温度等9)p=p(x,y,z,t)p=p(x,y,z,t)T=T(x,y,z,t)变数x

6、,y,z,t，称为欧拉变数，当x,y,z固定，t改变时，（7）式中的函数代表空间中固定点上速度随时间的变化规律，当t固定，x,y,z改变时，它代表的是某一时刻中速度在空间的分布规律。应该指出，有（7）式确定的速度是定义在空间点上的，它们是空间点的坐标x,y,z的函数，所以我们研究的是场，如速度场，压力场、密度场等。因此当我们采用欧拉观点描述运动时，就可以广泛地利用场论的知识。若场内函数不依赖于失径r则称之为均匀场；反之称为不均匀场。若场内函数不依赖时间t则称为定常场，反之称不定常场。三随体导数定义求解假定速度函数（7）具有一阶连续偏导数，现在从（7）式出发求质点的加速度f，设dt某质点在场内运

7、动，其运动轨迹为厶。在t时刻，给质点位于M点，速度为v（M,t），过了At时间后，该质点运动于M点，速度为v（Mt+At）。根据定义，加速度的表达式是dt=limAttOv（M:t+At）-v（M,t）At10）从（10）式可以看到，速度的变化亦即加速度的获得主要是下面两个原因引起的。一方面，当质点由M点运动M点时，时间过去了At，由于场的不定常性速度将发生变化。另一方面与此同时M点在场内沿迹线移动了MM距离，由于场的不均匀性亦将引起速度的变化。根据这样的考虑，将（10）的右边分成两部分dvdtv（M；t+At）-v（M；t）v（M；t）-v（M,t）TOC o 1-5 h zlim+limA

8、ttOAtAttOAtv（M；t+At）-v（M；t）=limAttOAt（1十MM十v（Mt）-v（M,t） HYPERLINK l bookmark34 o Current Document +limlimAttOAtMMtOMM1）右边第一项当At0时MMQt，这一项代表由于场的不定常性引起的速度变化，称为局部导数或就地导数；右边第二项是VQv（M,t）Qs，它代表由于场的不均匀性引起的速度变化，称为位变导数或对流导数，其中7代表沿s方向移动单位长os度引起的速度变化，而如今在单位时间内移动了V的距离，因此s方向上的速度变化是ovV。这样总的速度变化即加速度就是局部导数和位变导数之和，称

9、之为随体导数。于是os有dvovov=+V- HYPERLINK l bookmark28 o Current Document dtotos从场论中得知空=(s.V)vOs0其中s是曲线L的单位切向矢量。考虑到Vs=v，得00空=v+(vV)vdtOt这就是矢量形式的加速度的表达式。12)13)在直角坐标系中采取下列形式dudtOuOuOuOu+u+v+wOtOxOyOzdvdtOvOvOvOv+u+v+wOtOxOyOz14)dwOw+udtOtOwOwOw+v+wOxOyOz级数求解从级数展开角度来求解欧拉下的加速度的表达式，用欧拉方法描述流场时，一、某空间点上的流体质点的速度是时间的函

10、数，所以速度随时间变化，二、原来在某空间点上的流体质点经过了At后到达了另一空间点，若这两点的速度不同，那么由于迁移，它也会有速度的变化。设在t时刻，位于P(x,y,z)点的一个微团具有速度u,v,w。经At后，该微团移到(x+uAt,y+vAt,z+wAt)。令uf(x,y,z,t)经过At后，u变成了u+Au，即u+Au=/x+uAt,y+vAt,z+wAt,t+15)16)=/(x,y，z,T)+生uAt+生vAt+FwAt+FAt+(AT)的高阶项IQxBydzdt丿略去高阶项，仅保留一阶项，得ATBTBxByBzAUBUBUBUBU=+-U+-V+-AtBtBxByBz即w竺=BF+

11、茎U+BFv+BFw此式右侧第一项是微团在(x,y,z)处其速度随时间的变化率，即当地导数或局部导数。后三项是由于微团流向不同的领点是而出现的速度变化率，即迁移导数。总的称为流体质点的随体导数。同样，V,W也有这样的随体导数dvBvBvBvBv=+U+V+WdTBTBxByBzdwBwBwBwBw=+U+V+WdTBTBxByBz微分求解随体导数的求解还可以通过直接微分的方式得到。设与轨迹L相对应的运动方程是r=r(T)x=x(T)y=y(T)z=z(T)于是速度函数可写成17)v=v(x(T),y(T),z(T),T)对v做复合函数微分，并考虑到drdTdx=UdTdy=VdTdz=WdT于

12、是得到賞二色+空竺+空空+空主dtdtdxdtdydtdzdt+u+v+wdtdxdydz嗚+(vV)v(18)上述将随体导数分解为局部导数和位变导数之和的方法对于任何矢量a和任何标量申都是成立的，此时有=+(vV)adtdt空+(vV)(pdtdt四两种流动描述方法之间的关系(19)(20)欧拉方法在数学处理上的最大困难是方程式的非线性，而拉各朗日方法中的加速度项则为线性。但是直接应用拉各朗日型的基本方程解决流体力学问题是困难的，因此在处理流动问题是，常常必须用拉各朗日的观点而却应用欧拉观点的方法，这里就必须研究拉各朗日与欧拉两种系统之间的变化关系。为此引用雅克比行列式(Jacobian)。

13、J(t)=det空dgi1)拉各朗日变数g与欧拉变数x可以互换的唯一条件是:J(t)H0,g(2雅克比行列式的时间导数：J亠=(Vu)Jdtdx(22)例1讨论不可压缩流体的数学表示根据定义，质点的密度在运动过程中不变的流体的称为不可压缩流体。换而言之，对于不可压缩流体而言，密度的随体导数为零，即dt这就是不可压缩流体的数学表示。应该特别指出不可压缩流体的数学表示晋=0和不可压缩匀质流体的数学表示P=常数是不同的，不能把它们混为一谈。字=0表示每个质点的密度在它运动的全过程中不变。但是这个质点的密度和那个质点的密度可以不同，因此不可压缩流体的密度不一定是常数，只有既为不可压缩流体同时又是匀质时

14、密度才是处处时时都是同一常数。这个事实也可推导如下：dP=0（不可压缩），Vp=0（匀质），根dt据随体导数的表达式，可知：叟=0，于是P=常数。ct迹线、流线迹线是流体质点在空间中运动时所描绘出来的曲线。流星在夜空划出的一道光线，五光色的烟火图案是迹线的例子，迹线给出了同一质点在不同时刻的空间位置和速度方向。如果流体运动速度已经给出，即v=v（r,t），则迹线方程可通过求解下列微分方程组而得到dx=udt(x,y,z,t)dy=v(dtx,y,z,t)(23)dydx二d二dtu(x,y,z,t)v(x,y,z,t)w(x,y,z,t)式中t是自变量，x,y,z都是t的函数。积分后在所得到的

15、表达式中消去时间t后即得到迹线的方程。流线是某一瞬时流场内一条想象的曲线，该曲线上各点的速度方向和曲线在该点的切线方向重合。这里定义的流线类似于物理学中定义的电场线和磁场线，它们都是利用矢量线来几何地表示一个矢量场。流线是在同一时刻由不同流体质点所组成的曲线，它给出该时刻不同流体质点的运动方向。应该着重指出流线是某一瞬时的曲线，另一瞬时若流场改变了，通过同一点的流线也会改变。设d=dxi+dyj+dzk是流线上某点的线元，而v=ui+vj+wk是该点的速度矢量,根据流线定义，dl和v相互平行，于是dlxv=0由上式可得(24)dxdydzu(x,y,z,t)v(x,y,z,t)w(x,y,z,

16、t)这是t时刻的流线应该满足的微分方程，t在积分时硬挨看作常数处理。从上述分析可知，迹线和流线是两个具有不同内容和意义的曲线。迹线是同一质点在不同时刻的空间位置连成的曲线，它与拉格朗日观点相关联；而流线则是同一时刻不同质点所组成的曲线，它是与欧拉观点相联系的。在非定常运动中，迹线和流线一般说来是不重合的；但是在定常流动中，由于流动与时间t无关，流线不随时间而改变，流体质点沿着流线运动，因此，流线即是迹线。例设有一平面流场，其速度表达式是u=x+1，v=-y+1，w=0，求t=0时，过（1，1）点的流线和迹线。解：（1）流线的微分方程是dx_dyx+1y+1式中t是参数，积分得（x+1）（-y+

17、1）=c以t=0时，x=y=-l代入上式，可确定积分常数c=1，所以所求流线方程为xy=1这是一条双曲线，根据题意所求流线应是第三象限的分支曲线（2）迹线应满足的方程是dxdy=x+1，=y+1dtdt这里t是自变量，以上两个方程的解分别是x=cett1，y=cet+t112以H时，x=y-1代入得：=叮0，消去t后得X+y=2这是一条直线。同时，由此例可见，在不定常运动时，迹线和流线一般是不重合的物质积分的随体导数线段元，面积元和体积元的随体导数dpd(Qor=lrr丿=vv=ovdtdt1212(2St=v5tdt5)d一0s=0s(Vv)一(Vv)0sdt对物质体积分的随体导数J申St=Jd6t+f申St=J空+申divvStdtTtdtTdtTIdt丿又因为空+申divv=空+vV申+Vv=