声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2243|回复: 1

[基础理论] 流体运动的拉格朗日描述和欧拉描述及声波方程

[复制链接]
发表于 2018-4-3 16:32 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
  摘要:拉格朗日描述与欧拉描述乃描述流体运动的两种体系。两种描述对于线性声学几无分别,但对非线性声学则不然。人们或惯用欧拉描述分析求解声学问题,但对有些非线性问题若采用拉格朗日描述更直接了当,更易求得解析解。本文概论概述两种描述,探讨拉格朗日描述下流体运动的基本方程,最后给出一维的非线性声波方程。

  何谓拉格朗日描述和欧拉描述?
  欲描述流体的运动,既可选择拉格朗日坐标体系,也可采用欧拉坐标体系。拉格朗日描述取初始时刻(如时间t=0)流体质点的三维位置矢量R=(a,b,c)作为空间变量,并标记该流体质点。随着流体的运动,t时刻的质点运动至新的三维空间位置r=(x,y,z)。显然,r与R具有函数关系:r=r(R,t),或用笛卡尔坐标分量表示
1.png
  该关系其实给出了流体质点在空间的运动轨迹。由此可求出质点运动的速度——流速v
2.png
  显然,不同的质点具有不同的速度,即流速v不但是时间的函数,也是初始位置R的函数。不仅如此,任何描述流体状态的物理量(例如密度ρ)均属流体质点所“携带”的。不同质点的运动状态当然不同;故流体的任何物理量不仅是时间t的函数,也是初始坐标矢量R的函数,例如,密度ρ=ρ(R,t)。此种以流体质点初始位置R和时间t描述流体运动的方法,即所谓的拉格朗日描述,其中R=(a, b, c)是拉格朗日坐标。其实,拉格朗日坐标相当于质点的“标签”。所以,拉格朗日描述记录了质点的时间演化,既可追溯质点的既往,也可预测质点的未来。

  与一般的流动不同,声学问题中的质点运动是围绕初始位置R=(a, b, c)的质点振动。记振动位移矢量为δr=(ξ,η,ζ),则r=R+δr,或者,
3.png
  当然,位移δr也是初始位置R=(a, b, c)和时间t的函数:δr=δr(R, t),或用分量表示:
4.png
  与拉格朗日描述不同,欧拉描述乃现实主义的描述手法。它仅关注“现在”——t 时刻处于空间r=(x,y,z)位置的质点。在此描述下,所有物理量皆表为空间位置坐标r(欧拉坐标)和时间t的函数,例如密度ρ=ρ(r,t)。所以,欧拉描述所描绘的是流场在不同时刻的瞬时空间分布,其空间坐标r=(x,y,z)是纯数学的独立变量。

  两种描述的关系
  虽然如此,t 时刻位于空间位置r=(x,y,z)的流体质点并非“空穴来风”,而是从别处迁移而来。运动规律r=r(R,t)或公式(1)可对质点“追根溯源”:原则上,根据当前的空间位置r反推初始位置R。据此,欧拉描述和拉格朗日描述可以相互转换。为明确起见,特以上标(L)(E)分别标记拉格朗日和欧拉体系下的物理量,例如密度ρ(L)(R,t)ρ(E)(r,t)。据上论述,
5.png
  诸如声波这般流体运动,位移δr的量值一般很小(即,δr的范数||δr||<<1),因此可以对上式右端作泰勒展开:
6.png
  式中,欧拉空间的梯度算符:
7.png
  可见,已知位移δr=(ξ,η,ζ),就可从欧拉体系的物理量求得拉格朗日体系的对应量。反之,若函数r=r(R,t)连续可逆,可求得逆函数R=R(r,t),从而可把拉格朗日体系的物理量转换为欧拉体系的【注1】。特别是,对于小位移振动的情形,有如下展开
8.png
  式中,拉格朗日空间的梯度算符
9.png
  可见,拉格朗日描述的物理量ρ(L)(R,t)与欧拉描述的ρ(E)(r,t)仅相差量级O(||δr||)
  数学上,函数关系r=r(R,t)或公式(1)所表达的是以时间t为参变数的空间映射:
10.png
  矢量r在拉格朗日空间的梯度:
11.png
  是一个“矢量的矢量”——张量,刻画了映射的性质,式中I是单位张量。J也可用等价的雅可比矩阵形式表示:
12.png
  式中的下标表示对相应变量(a,b,c)的导数,例如:
13.png

  连续性方程
  设某质点的初始体积ΔV0(→0)t 时刻的体积ΔV。根据数学分析,雅可比矩阵J的行列式(记为|J|)是在映射(4)下质点体积元的缩放因子:
14.png
  设流体初始的质量密度为ρ0,t时刻的质量密度为ρ。根据质量守恒定律,ρ0ΔV0=ρΔV。将上式代入,得到拉格朗日坐标下的连续性方程
15.png
  对于一维的情形,根据定义(6)知|J|=1+ξa,公式(7)遂简为:
16.png

  运动方程
  理想流体的运动遵循欧拉方程:
17.png
  式中,P是流体压强。在拉格朗日体系中,速度v也是拉格朗日坐标R=(a,b,c)和时间t的函数:v = v(R, t)。所以,加速度dv/dt是对v的时间偏导数。如此,上列欧拉方程可改写为
18.png
  在此方程的等式两端同时点乘矢量:
19.png
  并利用微分关系
20.png
  得到
21.png
  同理,可得
22.png
  三个方程合而为一,可表为
23.png
  此即拉格朗日体系的流体运动方程。式中,雅可比矩阵J由公式(6)给出。或者,利用公式(5)用位移矢量δr表示,
24.png
  如以矩阵表示,则可表为矩阵形式
25.png
  在一维情形,J = 1+ξa,方程(8)于是简化为:
26.png
  其中已经利用了连续性方程(7a)。可见,拉格朗日描述下质点的一维运动方程是线性的。当然,质点运动仍是非线性的,集中体现在连续性方程(7)以及(下述的)状态方程中。

  一维理想流体的非线性声波方程
  理想流体(如气体)满足绝热状态方程:P=P(ρ)。利用连续性方程 (7a),则:
27.png
  代入一维运动方程(10)的右端,得到非线性声波方程
28.png
  其中P'(ρ)=dP/dρ。如果|ξa|<<1,则可对方程(11)作线性化处理,从而得到熟知的线性波动方程。若|ξa|<1,则有泰勒展开:
29.png
  式中,c0是熟知的线性声速,而β是非线性系数,衡量媒质本身(二次)非线性的大小。若仅计及二阶非线性小量,则方程(11)可近似为
30.png
  方程的左端是标准的一维波动方程,而右端的非齐次项是二次非线性的贡献。

  【注1】在诸如湍流等复杂流动情形,可逆性或许不存在,两种描述难于转换。

  本文来源于网易聲之韻博客,作者王新龍,南京大学声学研究所。

回复
分享到:

使用道具 举报

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-3-29 09:49 , Processed in 0.071922 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表