拉格朗日方法与欧拉方法

可以想象有一条河流，你站在岸边，将叠好的纸船放进了河里，你的目光跟随着小船的轨迹，以上情况用于描述小船运动轨迹的方法叫拉格朗日法，也是我们从小学到高中一直学习的方法，即跟踪一个物体的轨迹。小船飘走了，你走到桥上，站在最高点向下看小溪的潺潺流水，以上情况用于描述一个指定区域每时每刻变化的方法叫欧拉法。

生活中两种方法的比较

一开始接触这两种方法是在大三上流体力学课的时候，但是等理解过后发现，这两种方法更多的是给出了在现实生活中两种看问题的不同视角，不单单仅限于流体力学中使用。拉格朗日法是我们再熟悉不过的方法，从我们初中学习的动点问题，到高中的牛顿第二定律，都是对某一个质点进行研究，相当于拿一个摄像机跟踪这个质点。欧拉法应该是大家都陌生的一种方法，它涉及到一个概念，就是场。场可以理解为场地，这种描述的方法其实就是希望观察一个固定场景的变化。下面举一些好理解的例子供大家参考：

拍戏(我没拍过，脑补的)

家有儿女应该都看过(不怎么看电视剧),里面镜头跟随刘星动的拍摄可以理解为拉格朗日法。而在家里的拍摄可以理解为欧拉法。因为跟随刘星拍摄是需要镜头跟着刘星，对象不变，场景在变。而在家里拍摄，虽然人是来回走动，但是镜头不变，一直对着沙发和餐桌，这是场景不变，对象在变。

100m跑

我们在电视上观看100米跑决赛的时候会看到两种方式拍摄方式，第一种是摄影师拍摄影师跟着运动员在跑，这种就是拉格朗日法，比较费摄影师。还有一种就是一个镜头从上面俯瞰整个100m赛场，那就是欧拉法。

传记

看书的时候其实也有拉格朗日法和欧拉法，比如我们看的人物传记（普京传）就是拉格朗日手法描述的，而描述国家历史(明朝那些事儿)的就是欧拉法。

描述完上述内容，下面进行一些数学上公式的理解。

拉格朗日方法(质点法）

位置与时间

这种方法如前面所说，要求观察者随时随地跟踪每个流体质点，记录该质点的运动历程从而获得整体流体的运动规律。设在某一初始时刻$t_0$，流体质点的初始位置坐标为$(a,b,c)$，那么任意时刻$t$，质点的空间位置可以表示为：

$$\begin{gathered} x=x(a,b,c,t) \\ y=y(a,b,c,t) \\ z=z(a,b,c,t) \end{gathered}$$
如果给定$t$那么上式表示在给定时间不同质点的空间位置；如果$(a,b,c)$给定了，那么就表示给定质点在不同时刻的位置变化，即质点的轨迹。

速度与加速度

如果知道给定质点的位置，则可以根据定义确定质点的速度：
$$\begin{gathered} u=\underset{\Delta t\to 0}{\lim }\frac{x(a,b,c,t+\Delta t)-x(a,b,c,t)}{\Delta t}=\frac{\partial x(a,b,c,t)}{\partial t} \\ v=\underset{\Delta t\to 0}{\lim }\frac{y(a,b,c,t+\Delta t)-y(a,b,c,t)}{\Delta t}=\frac{\partial y(a,b,c,t)}{\partial t} \\ w=\underset{\Delta t\to 0}{\lim }\frac{z(a,b,c,t+\Delta t)-z(a,b,c,t)}{\Delta t}=\frac{\partial z(a,b,c,t)}{\partial t} \end{gathered}$$
这里用时间偏导数，是因为流体质点是给定的。同样可以得到加速度表达式：
$$\begin{gathered} a_x=\underset{\Delta t\to 0}{\lim }\frac{u(a,b,c,t+\Delta t)-u(a,b,c,t)}{\Delta t}=\frac{\partial u(a,b,c,t)}{\partial t}=\frac{{\partial }^{2}x}{\partial t^2} \\ {a}_y=\underset{\Delta t\to 0}{\lim }\frac{v(a,b,c,t+\Delta t)-v(a,b,c,t)}{\Delta t}=\frac{\partial v(a,b,c,t)}{\partial t}=\frac{{\partial }^{2}y}{\partial t^2} \\ {a}_z=\underset{\Delta t\to 0}{\lim }\frac{w(a,b,c,t+\Delta t)-w(a,b,c,t)}{\Delta t}=\frac{\partial w(a,b,c,t)}{\partial t}=\frac{{\partial }^{2}z}{\partial t^2} \end{gathered}$$
这种跟踪质点的方法实际上是理论力学中质点系法的直接延伸，也就是我们日常最长接触的记录方法。这种方法概念清晰，便于物理定律的推广，缺点是记录了太多的资料。例如每年到汛期，人们仅想了解武汉段的长江水情，但用这种方法描述的时候，必须把通过武汉段长江内的所有水的质点的来龙去脉搞清楚，全程跟踪并记录各个质点的流动过程，才能刻画武汉段长江的流动特征。这样做实际上对许多不在武汉段长江内的水质点记录是无用的。

欧拉方法（流场法）

对于速度的表示

为了避免拉格朗日法不必要的资料记录，我们不标识质点，相反标识空间点，观察者相对于空间点不动，记录不同时刻，不同质点通过固定空间点的快慢，观察者直接记录的量是不同时刻通过给定空间点的质点速度值，如将每个空间点上布置一个观察者，记录各个空间点在各时刻的质点速度值，这样就可以对所考察的流动区域特征给出全面了解。这种方法虽然标识的是空间点，但是研究的仍然是流体质点。任一空间点的位置坐标是$(x,y,z)$，在$t$时刻由位于该处的观察者直接记录的通过该空间点流体质点速度为：
$$\begin{gathered} u=u(x,y,z,t) \\ v=v(x,y,z,t) \\ w=w(x,y,z,t) \end{gathered}$$
如果给定$(x,y,z)$，$t$变化，表示不同时刻不同流体质点通过同一空间点的速度(快慢)；若给定$t$，$(x,y,z)$在变化，表示给定时间，占据不同空间点流体质点的速度。空间点速度本质上指的是$t$瞬时恰好占据该空间点流体质点的速度。物理上，把一个布满了某种物理量的空间成为场。流体流动所占据的空间即为流场。如果物理量是速度，那就是速度场。

对加速度的推导

如果速度是具有一阶连续偏导数的函数，那么根据定义可以确定质点的加速度。在欧拉方法中，局部跟踪流体质点运动的跟随性条件是：
$$\frac{dx}{dt}=u，\frac{dy}{dt}=v，\frac{dz}{dt}=w$$
在研究流体质点加速度时，就需要局部跟踪所考察的质点速度变化情况。设某一流体质点在$t$时刻位于流场中$M(x,y,z)$，经过时段$\Delta t$移动到了$N(x+\Delta x,y+\Delta y,z+\Delta z)$,根据质点加速度定义，有：
$$\begin{gathered} a=\frac{dV}{dt}=\underset{\Delta t\to 0}{\lim }\frac{\Delta V}{\Delta t}=\underset{\Delta t\to 0}{\lim }\frac{V(N,t+\Delta t)-V(M,t)}{\Delta t} \\ =\underset{\Delta t\to 0}{\lim }\frac{V(N,t+\Delta t)-V(N,t)}{\Delta t}+\underset{\Delta t\to 0}{\lim }\frac{V(N,t)-V(M,t)}{\Delta t} \end{gathered}$$
此式表明，质点的加速度由两部分组成，第一部分时在$\Delta t$时间段内通过固定空间点的速度变化引起的加速度，第二部分是在给定时间质点通过相邻空间点速度变化引起的加速度。
对第一部分用泰勒展开：
$$V(N,t+\Delta t)=V(N,t)+\frac{\partial V(N,t)}{\partial t}\Delta t+O(\Delta t^2)$$
由此可得：
$$\underset{\Delta t\to 0}{\lim }\frac{V(N,t+\Delta t)-V(N,t)}{\Delta t}=\underset{\Delta t\to 0}{\lim }\frac{\frac{\partial V(N,t)}{\partial t}\Delta t}{\Delta t}=\frac{\partial V}{\partial t}$$
对于第二部分相当于固定时间由相邻不同空间位置流体质点速度变化引起的加速度，同样可以根据泰勒展开获得：
$$\begin{gathered} V(N,t)=V(x+\Delta x,y+\Delta y,z+\Delta z,t) \\ =V(x,y,z,t)+\frac{\partial V(x,y,z,t)}{\partial x}\Delta x+\frac{\partial V(x,y,z,t)}{\partial y}\Delta y+\frac{\partial V(x,y,z,t)}{\partial z}\Delta z+O(\Delta x^2,...) \end{gathered}$$
同理，带入第二部分，有：
$$\underset{\Delta t\to 0}{\lim }\frac{V(N,t)-V(M,t)}{\Delta t}=\underset{\Delta t\to 0}{\lim }\frac{\Delta x}{\Delta t}\frac{\partial V(M,t)}{\partial x}+\underset{\Delta t\to 0}{\lim }\frac{\Delta y}{\Delta t}\frac{\partial V(M,t)}{\partial y}+\underset{\Delta t\to 0}{\lim }\frac{\Delta z}{\Delta t}\frac{\partial V(M,t)}{\partial z}$$
根据流体的跟随性(跟踪的是同一个质点)
$$\underset{\Delta t\to 0}{\lim }\frac{V(N,t)-V(M,t)}{\Delta t}=u\frac{\partial V(M,t)}{\partial x}+v\frac{\partial V(M,t)}{\partial y}+w\frac{\partial V(M,t)}{\partial z}$$
综合起来，可以得到在欧拉框架下，流体质点的全加速度：
$$\begin{gathered} a=\frac{dV}{dt}=\frac{\partial V}{\partial t}+u\frac{\partial V(M,t)}{\partial x}+v\frac{\partial V(M,t)}{\partial y}+w\frac{\partial V(M,t)}{\partial z} \\ =\frac{\partial V}{\partial t}+（V\cdot \nabla ）V \end{gathered}$$
在全加速度表达式中，等式右边第一项表示速度对时间的偏导数，是由流场的非定常性引起的，称为局部加速度或当地加速度；第二项表示因流体质点位置变化引起的加速度，称为迁移加速度或对流加速度。两者相加合成全加速度或者随体加速度(跟随流体质点的加速度)。

转换公式实例

拉格朗日描述式变换成欧拉描述式(L-E变换)

给定拉格朗日位移描述式：
$$x_1=a_1{e}^{-\frac{2t}{k}},x_2=a_2{e}^{-\frac{t}{k}},x_3=a_3{e}^{-\frac{t}{k}}$$
求欧拉速度场。
解：先求出拉格朗日速度表达式：
$$\begin{gathered} u_1=\frac{\partial x_1}{\partial t}=-\frac{2a_1}{k}{e}^{-\frac{2t}{k}} \\ {u}_2=\frac{\partial x_2}{\partial t}=-\frac{a_2}{k}{e}^{-\frac{t}{k}} \\ {u}_3=\frac{\partial x_3}{\partial t}=-\frac{a_3}{k}{e}^{-\frac{t}{k}} \end{gathered}$$
目的是消去$a$,从位移解出$a$：
$$a_1=x_1{e}^{\frac{2t}{k}},a_2=x_2{e}^{\frac{t}{k}},a_3=x_3{e}^{\frac{t}{k}}$$
带入速度表达式得到：
$$u_1=-\frac{2}{k}{x}_1，u_2=-\frac{1}{k}{x}_1，u_3=-\frac{1}{k}{x}_1$$

欧拉描述式转换成拉格朗日描述式(E-L转换)

已知欧拉速度场:
$$u_1=x_1+t，u_2=x_2+t，u_3=0$$
初始条件为：
$$t=0,x_1=a_1,x_2=a_2,x_3=a_3$$
求质点位移和速度的拉格朗日表达式。
解：
$$\begin{gathered} u_1=\frac{\partial x_1}{\partial t}=x_1+t \\ {u}_2=\frac{\partial x_2}{\partial t}=x_2+t \\ {u}_3=\frac{\partial x_3}{\partial t}=0 \end{gathered}$$
该方程的一般解为：
$$\begin{gathered} x_1=c_1{e}^t-t-1 \\ {x}_2=c_2{e}^t-t-1 \\ {x}_3=c_3 \end{gathered}$$
带入初始条件，得到积分常数：$c_1=a_1+1,c_2=a_2+1,c_3=a_3$
得到拉格朗日表达式为：
$$\begin{gathered} x_1=(a_1+1)e^t-t-1 \\ {x}_2=(a_2+1)e^t-t-1 \\ {x}_3=a_3 \end{gathered}$$
速度的拉格朗日表达式为：
$$\begin{gathered} u_1=(a_1+1)e^t-1 \\ {u}_2=(a_2+1)e^t-1 \\ {u}_3=0 \end{gathered}$$