描述流体运动地两种方法

实用文档

描述流体运动的两种方法

(姓名：张旺龙 学号：308081183 专业：流体力学)

引言:

描述流体运动的两种方法――拉各朗日方法和欧拉方法

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

一 拉格朗日方法

现在我们将上述描写运动的拉各朗日观点和方法用数学式子表达出来，为此首先必须用某种数学方法区别不同的流体质点。通常利用初始时刻流体质点的坐标作为区分不同流体质点的标志。设初始时刻0tt时，流体质点的坐标是（a,b,c），它可以是曲线坐标，也可以是直角坐标,,000xyz，重要的是给流体质点以标号而不在于采取什么具体的方式。我们约定采用a,b,c三个数的组合来区别流体质点，不同的a,b,c代表不同的质点。于是流体质点的运动规律数学上可表示为下列矢量形式：

,,,abctrr （1）

其中r是流体质点的失径。在直角坐标系中，有

,,,xxabct ,,,yyabct ,,,zzabct （2）

变数a,b,c,t称为拉各朗日变数。在式（2）中，如果固定a,b,c而令t改变，则得某一流体质点的运动规律。如果固定时间t而令a,b,c改变，则得同一时刻不同流体质点的位置分布。应该指出，在拉各朗日观点中，失径函数r的定义区域不是场，因为它不是空间坐标的函数，而是质点标号的函数。

现在从（1）式出发来求流体质点的速度和加速度。假设由（1）式确定的函数具有二阶连续偏导数。速度和加速度是对于同一质点而言的单位时间内位移变化率及速度变化率，设v，v分别表示速度矢量和加速度矢量，则

,,,rabcttv （3）

22,,,rabcttv （4）

既然对同一质点而言，a,b,c不变，因此上式写的是对时间t的偏导数。在直角坐标系中，速度和加速度的表达式是

,,,xabctut ,,,yabctvt ,,,zabctwt （5）

及 实用文档

22,,,uxabctt 22,,,vyabctt 22,,,wzabctt （6）

二 欧拉方法

现在来介绍描写流体运动的另一种观点和方法，即欧拉方法。和拉各朗日方法不同，欧拉方法不同，欧拉方法的着眼点不是流体质点而是空间点。设法在空间中的每一点上描述出流体运动随时间的变化状况。如果，每一点的流体运动都已知道，则整个流体的运动状况也就清楚了，那么应该用什么物理量来表现空间点上流体运动的变化情况呢？因为不同时刻将有不同流体质点经过空间某固定点，所以站在固定点上就无法观测和记录掠过的流体质点以前和以后的详细历史。也就是说我们无法象拉各朗日方法那样直接测量出每个质点的位置随时间的变化情况。虽然如此，不同时刻经过固定点的流体质点的速度是可以测出的，这样采用速度矢量来描写固定点上流体运动的变化状况就是十分自然的了。考虑到上面所说的情形，欧拉方法中流体质点的运动规律数学上可表示为下列矢量形式：

,tvvr （7）

在直角坐标系中有：

,,,uuxyzt ,,,vvxyzt ,,,wwxyzt （8）

要完全描述运动流体的状况还需要给定状态函数压力、密度、温度等

,,,ppxyzt ,,,xyzt ,,,TTxyzt （9）

变数,,,xyzt，称为欧拉变数，当,,xyz固定，t改变时，（7）式中的函数代表空间中固定点上速度随时间的变化规律，当t 固定，,,xyz改变时，它代表的是某一时刻中速度在空间的分布规律。应该指出，有（7）式确定的速度是定义在空间点上的，它们是空间点的坐标,,xyz的函数，所以我们研究的是场，如速度场，压力场、密度场等。因此当我们采用欧拉观点描述运动时，就可以广泛地利用场论的知识。若场内函数不依赖于失径r则称之为均匀场；反之称为不均匀场。若场内函数不依赖时间t 则称为定常场，反之称不定常场。

三 随体导数

3.1 定义求解

假定速度函数（7）具有一阶连续偏导数，现在从（7）式出发求质点的加速度ddtv，设某质点在场内运动，其运动轨迹为L。在t 时刻，给质点位于M点，速度为,Mtv，过了t时间后，该质点运动于M点，速度为,Mttv。根据定义，加速度的表达式是

0,,limtMttMtddttvvv （10）

从（10）式可以看到，速度的变化亦即加速度的获得主要是下面两个原因引起的。一方面，当质点由M点运动M点时，时间过去了t，由于场的不定常性速度将发生变化。另一方面与此同时M点在场内沿迹线移动了MM距离，由于场的不均匀性亦将引起速度的实用文档

变化。根据这样的考虑，将（10）的右边分成两部分

ddtv0,,limtMttMttvv0,,limtMtMttvv

0,,limtMttMttvv00,,limlimtMMMtMtMMtMMvv (11)

右边第一项当0t时MM，因此它是,Mttv，这一项代表由于场的不定常性引起的速度变化，称为局部导数或就地导数；右边第二项是,MtVsv，它代表由于场的不均匀性引起的速度变化，称为位变导数或对流导数，其中sv代表沿s方向移动单位长度引起的速度变化，而如今在单位时间内移动了V的距离，因此s方向上的速度变化是Vsv。这样总的速度变化即加速度就是局部导数和位变导数之和，称之为随体导数。于是有

dVdttsvvv （12）

从场论中得知

0ssvv

其中0s是曲线L的单位切向矢量。考虑到0Vsv，得

ddttvvvv （13）

这就是矢量形式的加速度的表达式。

在直角坐标系中采取下列形式

duuuuuuvwdttxyz

dvvvvvuvwdttxyz （14）

dwwwwwuvwdttxyz

3.2 级数求解

从级数展开角度来求解欧拉下的加速度的表达式，用欧拉方法描述流场时，一、某空间点上的流体质点的速度是时间的函数，所以速度随时间变化，二、原来在某空间点上的流体质点经过了t后到达了另一空间点，若这两点的速度不同，那么由于迁移，它也会有速度实用文档

的变化。

设在t时刻，位于,,Pxyz点的一个微团具有速度,,uvw。经t后，该微团移到,,xutyvtzwt。令,,,ufxyzt

经过t后，u变成了uu，即

uu,,,fxutyvtzwttt

,,,fxyztffffutvtwttxyztt的高阶项 （15）

略去高阶项，仅保留一阶项，得

uffffuvwttxyz

即 uuuuuuvwttxyz （16）

此式右侧第一项是微团在,,xyz处其速度随时间的变化率，即当地导数或局部导数。后三项是由于微团流向不同的领点是而出现的速度变化率，即迁移导数。总的称为流体质点的随体导数。

同样，,vw也有这样的随体导数

dvvvvvuvwdttxyz

dwwwwwuvwdttxyz

3.3 微分求解

随体导数的求解还可以通过直接微分的方式得到。设与轨迹L相对应的运动方程是

trr

或

xxt yyt zzt

于是速度函数可写成

,,,xtytzttvv （17）

对v做复合函数微分，并考虑到

ddtrv

即 dxudt ， dyvdt ， dzwdt 实用文档

于是得到

ddxdydzdttxdtydtzdtvvvvv

uvwtxyzvvvv

tvvv （18）

上述将随体导数分解为局部导数和位变导数之和的方法对于任何矢量a和任何标量都是成立的，此时有

ddttaava (19)

ddttv (20)

四 两种流动描述方法之间的关系

欧拉方法在数学处理上的最大困难是方程式的非线性，而拉各朗日方法中的加速度项则为线性。但是直接应用拉各朗日型的基本方程解决流体力学问题是困难的，因此在处理流动问题是，常常必须用拉各朗日的观点而却应用欧拉观点的方法，这里就必须研究拉各朗日与欧拉两种系统之间的变化关系。为此引用雅克比行列式（Jacobian）。

detiixJt (21)