欢迎光临散文网 会员登陆 & 注册

数理方程来了||数理方法

2021-03-22 17:08 作者:湮灭的末影狐  | 我要投稿

//前面六章是复变函数与积分变换的主要内容。

//现在开始下半篇数学物理方程了

//开学之后生产力急剧下降

//但是这个学期应该是能更完的

7.1 数学物理方程的导出

有许多物理模型,描述它们的方程具有类似的形式。以下将导出一些经典物理模型中的微分方程,并将它们分为三大类:波动方程、输运方程与稳定场方程

7.1.1 波动方程

波动方程广泛出现在弹性力学、电磁学、力学等学科的相关模型中,它通常具有以下形式:

u_%7Btt%7D-a%5E2%5CDelta%20u%3D0

如果存在波源,则波动方程为:

u_%7Btt%7D-a%5E2%5CDelta%20u%3Df(%5Cvec%7Br%7D%2Ct)

其中,%5CDelta是拉普拉斯算符,根据模型的实际情况可能是1至3维。a为波速。一维波动方程通常写为

u_%7Btt%7D-a%5E2%20u_%7Bxx%7D%3Df(x%2Ct)

其实上面这种写法是数理方法书上的...我个人还是更习惯下面的写法:

%5Cnabla%5E2%20%5Cpsi%20-%5Cfrac%7B1%7D%7Bc%5E2%7D%5Cfrac%7B%5Cpartial%5E2%5Cpsi%7D%7B%5Cpartial%20t%5E2%7D%3Df(%5Cvec%7Br%7D%2Ct)

或者

%5Cfrac%7B%5Cpartial%5E2%5Cpsi%7D%7B%5Cpartial%20x%5E2%7D%20-%5Cfrac%7B1%7D%7Bc%5E2%7D%5Cfrac%7B%5Cpartial%5E2%5Cpsi%7D%7B%5Cpartial%20t%5E2%7D%3Df(x%2Ct)

这些也是大部分物理书籍、论文里波动方程的常见形式。波速为c.

例如,在电磁学中,有电磁势满足的波动方程:

%5Cnabla%5E2%20%5Cvec%7BA%7D%20-%20%5Cfrac1%7Bc%5E2%7D%20%5Cfrac%7B%5Cpartial%20%5E2%20%5Cvec%7BA%7D%7D%7B%5Cpartial%20t%5E2%7D%3D-%5Cmu_0%20%5Cvec%7Bj%7D

%5Cnabla%5E2%20%5Cphi%20%20-%20%5Cfrac1%7Bc%5E2%7D%20%5Cfrac%7B%5Cpartial%20%5E2%20%5Cphi%7D%7B%5Cpartial%20t%5E2%7D%3D-%5Cfrac%7B%5Crho%7D%7B%5Cepsilon_0%7D

弹性力学中,固体中的 纵波与横波都满足波动方程,波速分别是

c_1%3D%5Csqrt%7B%5Cfrac%20Y%5Crho%7D%2C%5C%2Cc_2%3D%5Csqrt%7B%5Cfrac%20G%5Crho%7D

Y%2CG分别是材料的杨氏模量和切变模量。

再如张力为T,线密度%5Clambda的弦上,微小的横振动也满足波动方程,波速为%5Csqrt%7BT%2F%5Clambda%7D

再如理想传输线(我们以此为例展示方程的推导过程):

设有平行导线,电阻和电磁波辐射忽略。二者之间单位长度的电容C,单位长度电感为L.

仅考虑两线上电荷线密度%5Clambda与电流I均等大反向的对称情况,则根据电流连续性有

%5Cfrac%7B%5Cpartial%20I%7D%7B%5Cpartial%20x%7D%20%3D-%5Cfrac%7B%5Cpartial%20%5Clambda%20%7D%7B%5Cpartial%20t%7D%20

再由电容电感的定义,可得

%5Clambda%3DC%5Cphi

%5Cfrac%7B%5Cpartial%20%5Cphi%20%7D%7B%5Cpartial%20x%7D%3D%20-L%5Cfrac%7B%5Cpartial%20I%7D%7B%5Cpartial%20t%7D

由上两式,得到

%5Cleft%20%5C%7B%5Cbegin%7Barray%7D%7Bc%7D%0A%5Cfrac%7B%5Cpartial%20I%7D%7B%5Cpartial%20x%7D%20%3D-%5Cfrac%7B%5Cpartial%20%5Clambda%20%7D%7B%5Cpartial%20t%7D%20%5C%5C%0A%5Cfrac%7B%5Cpartial%20%5Clambda%20%7D%7B%5Cpartial%20x%7D%3D%20-LC%5Cfrac%7B%5Cpartial%20I%7D%7B%5Cpartial%20t%7D%0A%5Cend%7Barray%7D%5Cright.

前式对x%0A微分,后式对t微分,得到

%0A%5Cfrac%7B%5Cpartial%5E2%20I%7D%7B%5Cpartial%20x%5E2%7D%20-LC%20%5Cfrac%7B%5Cpartial%5E2%20I%20%7D%7B%5Cpartial%20t%5E2%7D%0A%3D0

这是一维波动方程,波速为1%2F%5Csqrt%7BLC%7D.

有趣的是,电容、电感这样的参数应当是完全由传输线的形状与空间相对位置决定的,并且根据各自的定义有各自不同的求法,但我们在这里可以断言它们的关系一定是

1%2F%5Csqrt%7BLC%7D%3D1%2F%5Csqrt%7B%5Cmu_0%5Cepsilon_0%7D%3Dc

其中c是光速。这是因为,电线的传输速度当然应该是光速。

7.1.2 输运方程

输运方程的形式通常是

u_t%20-%20a%5E2%20%5CDelta%20u%3Df(%5Cvec%7Br%7D%2Ct%20)

它出现在热传导与扩散的模型中,u在两种模型中分别是温度、扩散物浓度。这里只简单解释各项的含义:对时间一阶微分项代表这一点强度随时间的变化速率;%5CDelta%20u一项描述这一点受到周围的影响;右边则是这一点的源强度。

7.1.3 稳定场方程

静电场、静磁场都满足稳定场方程。经典的泊松方程就属于稳定场方程,形如

%5Cnabla%5E2%5Cpsi%20%3D%20f(%5Cvec%20r)

可以描述静电场电势、稳定温度场、稳定浓度场、稳恒电流场等模型。

如静电场的泊松方程:

%5Cnabla%20%5E2%20%5Cphi%20%20%3D%20-%5Crho%20%2F%5Cepsilon%20_0

稳恒电流场中也有类似的方程:

%5Cnabla%20%5E2%20%5Cvec%20A%20%3D%20-%5Cmu_0%20%5Cvec%20j

对于存在热源的稳定温度场,也有

%5Cnabla%5E2%20T%20%3D%20f(%5Cvec%20r)

许多模型都有相近的数学结构。

7.2 定解条件

除了决定系统的微分方程,我们还需要定解条件来确定系统物理参量与时间、空间的函数关系。定解条件分为三种:初值条件、边界条件、衔接条件

以下将以实例说明几种定解条件:

场景1:弦张力T,线密度%5Clambda%20,一端固定于x%3D0,在x%3Da处固定了质量为m的质点。令%5Cpsi表示横向位移,x%3Dl处是自由端,输入入射波%5Cpsi(t)%3DA%5Ccos%5Comega%20t.

这个场景中出现了以下的定解条件:

① x%3D0处被固定,所以恒有%5Cpsi(0%2Ct)%3D0x%3Dl处有入射波,所以%5Cpsi(l%2Ct)%20%3D%20A%5Ccos%20%5Comega%20t.

这被称为第一类边界条件:直接规定所研究物理量在边界上的值。

② 初始时刻,弦的各点位置%5Cpsi(x%2C0)%3Df(x)被称为初值条件:规定t%3D0时整个系统的状态。

③ 在质点所在的位置,由于质需受力产生加速度,在这一点出现了%5Cpsi(x%2Ct)的偏导不连续的情况。在该模型中质点的运动方程可以写为:

%5Clim%20_%7Bx%20%5Crightarrow%20a%5E%2B%7D%20%20T%20%5Cfrac%7B%5Cpartial%20%5Cpsi%7D%7B%5Cpartial%20x%7D%20-%0A%5Clim%20_%7Bx%20%5Crightarrow%20a%5E-%7D%20%20T%20%5Cfrac%7B%5Cpartial%20%5Cpsi%7D%7B%5Cpartial%20x%7D%20%3D%20m%20%5Cfrac%7B%5Cpartial%5E2%20%5Cpsi%7D%7B%5Cpartial%20t%5E2%7D%20%7C_%7Bx%3Da%7D

这被称为衔接条件,联系了不连续的质点两侧。

场景2:质量不可忽略,长度l,劲度系数k线密度为%5Clambda的弹簧,右端有质点m,左端是自由端。令%5Cpsi表示弹簧上某质点相对初始位置的偏移。同时质点处于力场F(%5Cpsi)%3D-k%5Cpsi中。

我们可以证明弹簧上某点的张力可以写为F%20%3D%20-kl%5Cfrac%7B%5Cpartial%20%5Cpsi%7D%7B%5Cpartial%20x%7D

 ① 由于自由端不受外力,有%5Cfrac%7B%5Cpartial%20%5Cpsi%7D%7B%5Cpartial%20x%7D%7C_%7Bx%3D0%7D%3D0

这被称为第二类边界条件,给出所研究物理量在边界法线方向上的导数。

 ② 质点受力,可列出运动方程:-kl%5Cfrac%7B%5Cpartial%20%5Cpsi%7D%7B%5Cpartial%20x%7D%7C_%7Bx%3Dl%7D-k%5Cpsi(l%2Ct)%3D%20m%20%5Cfrac%7B%5Cpartial%5E2%20%5Cpsi%7D%7B%5Cpartial%20t%5E2%7D%20%7C_%7Bx%3Dl%7D

这被称为第三类边界条件,给出所研究物理量及其在边界法线方向导数的关系。

7.3 数理方程分类

首先,形如

%5Csum_%7Bj%3D1%7D%5E%7Bn%7D%20%5Csum_%7Bi%3D1%7D%5E%7Bn%7D%20a_%7Bi%20j%7D%20u_%7Bx_%7Bi%7D%20x_%7Bj%7D%7D%2B%5Csum_%7Bi%3D1%7D%5E%7Bn%7D%20b_%7Bi%7D%20u_%7Bx_%7Bi%7D%7D%2Bc%20u%2Bf%3D0

被称为线性二阶偏微分方程。这是本课程重点研究的对象。

我们考虑两个自变量x%2Cy的情况,分类并化为标准形式。

a_%7B11%7D%20u_%7Bx%20x%7D%2B2%20a_%7B12%7D%20u_%7Bx%20y%7D%2Ba_%7B22%7D%20u_%7By%20y%7D%2Bb_%7B1%7D%20u_%7Bx%7D%2Bb_%7B2%7D%20u_%7By%7D%2Bc%20u%2Bf%3D0

这是此类方程的一般形式。作自变量代换:

%5Cleft%5C%7B%5Cbegin%7Barray%7D%7Bl%7D%0A%5Cxi%3D%5Cxi(x%2C%20y)%20%5C%5C%0A%5Ceta%3D%5Ceta(x%2C%20y)%0A%5Cend%7Barray%7D%5Cright.

(代换的雅可比行列式%5Cfrac%7B%5Cpartial(%5Cxi%2C%20%5Ceta)%7D%7B%5Cpartial(x%2C%20y)%7D%20%5Cneq%200.)

我们可以证明,变量代换后微分方程将化为新变量的二阶线性偏微分方程

A_%7B11%7D%20u_%7B%5Cxi%20%5Cxi%7D%2B2%20A_%7B12%7D%20u_%7B%5Cxi%20%5Ceta%7D%2BA_%7B22%7D%20u_%7B%5Ceta%20%5Ceta%7D%2BB_%7B1%7D%20u_%7B%5Cxi%7D%2BB_%7B2%7D%20u_%7B%5Ceta%7D%2BC%20u%2BF%3D0

变换前后方程各项系数满足关系:

%5Cleft%5C%7B%5Cbegin%7Barray%7D%7Bl%7D%0AA_%7B11%7D%3Da_%7B11%7D%20%5Cxi_%7Bx%7D%5E%7B2%7D%2B2%20a_%7B12%7D%20%5Cxi_%7Bx%7D%20%5Cxi_%7By%7D%2Ba_%7B22%7D%20%5Cxi_%7By%7D%5E%7B2%7D%20%5C%5C%0AA_%7B12%7D%3Da_%7B11%7D%20%5Cxi_%7Bx%7D%20%5Ceta_%7Bx%7D%2Ba_%7B12%7D%5Cleft(%5Cxi_%7Bx%7D%20%5Ceta_%7By%7D%2B%5Cxi_%7By%7D%20%5Ceta_%7Bx%7D%5Cright)%2Ba_%7B22%7D%20%5Cxi_%7By%7D%20%5Ceta_%7By%7D%20%5C%5C%0AA_%7B22%7D%3Da_%7B11%7D%20%5Ceta_%7Bx%7D%5E%7B2%7D%2B2%20a_%7B12%7D%20%5Ceta_%7Bx%7D%20%5Ceta_%7By%7D%2Ba_%7B22%7D%20%5Ceta_%7By%7D%5E%7B2%7D%20%5C%5C%0AB_%7B1%7D%3Da_%7B11%7D%20%5Cxi_%7Bx%20x%7D%2B2%20a_%7B12%7D%20%5Cxi_%7Bx%20y%7D%2Ba_%7B22%7D%20%5Cxi_%7By%20y%7D%2Bb_%7B1%7D%20%5Cxi_%7Bx%7D%2Bb_%7B2%7D%20%5Cxi_%7By%7D%20%5C%5C%0AB_%7B2%7D%3Da_%7B11%7D%20%5Ceta_%7Bx%20x%7D%2B2%20a_%7B12%7D%20%5Ceta_%7Bx%20y%7D%2Ba_%7B22%7D%20%5Ceta_%7Br%20y%7D%2Bb_%7B1%7D%20%5Ceta_%7Bx%7D%2Bb_%7B2%7D%20%5Ceta_%7By%7D%20%5C%5C%0AC%3Dc%20%5C%5C%0AF%3Df%20%0A%5Cend%7Barray%7D%5Cright.

(书上这一段的完整推导过程看起来真的非常恐怖,但是我们只要知道其中变量代换的原理就行了,剩下都是体力活)

我们只要合理地选取%5Cxi%2C%20%5Ceta就可以简化变换后的方程。例如,我们注意到,如果令方程

a_%7B11%7D%20z_%7Bx%7D%5E%7B2%7D%2B2%20a_%7B12%7D%20z_%7Bx%7D%20z_%7By%7D%2Ba_%7B22%7D%20z_%7By%7D%5E%7B2%7D%3D0

的特解为%5Cxi(x%2Cy)%2C%5Ceta(x%2Cy),就可以使A_%7B11%7D%3DA_%7B22%7D%3D0,从而完成简化。

而上面这个方程可以化为常微方程。设有曲线族z(x%2Cy)%3DC,则有

%5Cfrac%7B%5Cmathrm%7Bd%7D%20y%7D%7B%5Cmathrm%7Bd%7D%20x%7D%20%3D-%5Cfrac%7Bz_x%7D%7Bz_y%7D%20

再根据前面方程化为

a_%7B11%7D%5Cleft(-%5Cfrac%7Bz_%7Bx%7D%7D%7Bz_%7By%7D%7D%5Cright)%5E%7B2%7D-2%20a_%7B12%7D%5Cleft(-%5Cfrac%7Bz_%7Bx%7D%7D%7Bz_%7By%7D%7D%5Cright)%2Ba_%7B22%7D%3D0

所以曲线族满足

%5Cfrac%7B%5Cmathrm%7Bd%7D%20y%7D%7B%5Cmathrm%7B~d%7D%20x%7D%3D%5Cfrac%7Ba_%7B12%7D%5Cpm%5Csqrt%7Ba_%7B12%7D%5E%7B2%7D-a_%7B11%7D%20a_%7B22%7D%7D%7D%7Ba_%7B11%7D%7D

所以方程给出2个曲线族,分别对应%5Cxi(x%2Cy)%3DC.

再根据上式根号内对偏微分方程分类:%5Cleft%5C%7B%5Cbegin%7Barray%7D%7Bl%7D%0Aa_%7B12%7D%5E%7B2%7D-a_%7B11%7D%20a_%7B22%7D%3E0%EF%BC%8C%20%E5%8F%8C%E6%9B%B2%E5%9E%8B%20%5C%5C%0Aa_%7B12%7D%5E%7B2%7D-a_%7B11%7D%20a_%7B22%7D%3D0%EF%BC%8C%20%E6%8A%9B%E7%89%A9%E5%9E%8B%5C%5C%0Aa_%7B12%7D%5E%7B2%7D-a_%7B11%7D%20a_%7B22%7D%3C0%EF%BC%8C%20%E6%A4%AD%E5%9C%86%E5%9E%8B%0A%5Cend%7Barray%7D%5Cright.

以上几类方程经过(x%2Cy)%5Crightarrow(%5Cxi%2C%5Ceta)的变量代换之后,均可以化为标准形式。

双曲型的标准形式:

u_%7B%5Cxi%20%5Ceta%7D%3D-%5Cfrac%7B1%7D%7B2%20A_%7B12%7D%7D%5Cleft%5BB_%7B1%7D%20u_%7B%5Cxi%7D%2BB_%7B2%7D%20u_%7B%5Ceta%7D%2BC%20u%2BF%5Cright%5D

而如果再令%5Cxi%20%3D%20%5Calpha%2B%5Cbeta%2C%5C%3B%5Ceta%20%3D%20%5Calpha-%5Cbeta 则可以得到另一标准形式:

u_%7B%5Calpha%20%5Calpha%7D-u_%7B%5Cbeta%20%5Cbeta%7D%3D-%5Cfrac%7B1%7D%7BA_%7B12%7D%7D%5Cleft%5B%5Cleft(B_%7B1%7D%2BB_%7B2%7D%5Cright)%20u_%7B%5Calpha%7D%2B%5Cleft(B_%7B1%7D-B_%7B2%7D%5Cright)%20u_%7B%5Cbeta%7D%2B2%20C%20u%2B2%20F%5Cright%5D

例如,波动方程u_%7Btt%7D-c%5E2u_%7Bxx%7D%3D0就属于双曲型方程。根据前面理论,如果令%5Cxi%3Dx%2Bct%2C%5C%3B%5Ceta%3Dx-ct,则方程就可以化为:

u_%7B%5Cxi%5Ceta%7D%3D0

其通解是u%3Df(x%2Bct)%2Bg(x-ct)f%2Cg为任意函数。从这个解就很容易看出c作为“波速”的实际意义。

抛物型的标准形式:

%5Cxi%0A满足a_%7B11%7D%20z_%7Bx%7D%5E%7B2%7D%2B2%20a_%7B12%7D%20z_%7Bx%7D%20z_%7By%7D%2Ba_%7B22%7D%20z_%7By%7D%5E%7B2%7D%3D0,而%5Ceta不满足上式,可以证明抛物型的标准形式是

u_%7B%5Ceta%20%5Ceta%7D%3D-%5Cfrac%7B1%7D%7BA_%7B22%7D%7D%5Cleft%5BB_%7B1%7D%20u_%7B%5Cxi%7D%2BB_%7B2%7D%20u_%7B%5Ceta%7D%2BC%20u%2BF%5Cright%5D

椭圆型的标准形式:

除去一个负号外,椭圆型和双曲型应当是一样的:

u_%7B%5Cxi%20%5Ceta%7D%3D-%5Cfrac%7B1%7D%7B2%20A_%7B12%7D%7D%5Cleft%5BB_%7B1%7D%20u_%7B%5Cxi%7D%2BB_%7B2%7D%20u_%7B%5Ceta%7D%2BC%20u%2BF%5Cright%5D

但是这一次%5Cxi%2C%5Ceta变为复共轭的。为简化,取%5Cxi%20%3D%20%5Calpha%2B%5Cmathrm%20i%20%5Cbeta%2C%20%5Ceta%20%3D%20%5Calpha-%5Cmathrm%20i%20%5Cbeta%5C%3B(%5Calpha%2C%5Cbeta%5Cin%5Cmathbb%20R),得到

u_%7B%5Calpha%20%5Calpha%7D%2Bu_%7B%5Cbeta%20%5Cbeta%7D%3D-%5Cfrac%7B1%7D%7BA_%7B12%7D%7D%5Cleft%5B%5Cleft(B_%7B1%7D%2BB_%7B2%7D%5Cright)%20u_%7B%5Calpha%7D%2B%5Cleft(B_%7B2%7D-B_%7B1%7D%5Cright)%20u_%7B%5Cbeta%7D%2B2%20C%20u%2B2%20F%5Cright%5D

(先写这么多吧,这一章开始越来越麻烦了...以后不一定一篇笔记对应书上一章了)

参考文献

[1] 梁昆淼. 数学物理方法(第四版)[M]. 北京:高等教育出版社,2009.8,107~134.

数理方程来了||数理方法的评论 (共 条)

分享到微博请遵守国家法律