第五章--孤立波(共27页).doc
《第五章--孤立波(共27页).doc》由会员分享,可在线阅读,更多相关《第五章--孤立波(共27页).doc(27页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、精选优质文档-倾情为你奉上第五章 孤立波第一节 历史回顾1. 一个奇特的水波相传约170年以前,1834年的一天,在从爱丁堡到格拉斯哥的运河上,一位苏格兰海军工程师罗素(J.Scott Russell)观察到一种奇特的水波。在1844年发表的一份报告中,他描述了当年观察到的这种奇特水波,并称这种波为孤立波(Solitary wave)。他是这样描述的:“我看到两匹骏马拉着一条船沿运河迅速前进。当船突然停止时,随船一起运动的船头处的水堆并没有停止下来。它激烈地在船头翻动起来,随即突然离开船头,并以巨大的速度向前推进。一个轮廓清晰又光滑的水堆,犹如一个大鼓包,沿着运河一直向前推进,在行进过程中其形
2、状与速度没有明显变化。我骑马跟踪注视,发现它保持着起始时约30英尺(1英尺25.4厘米)长,1-1.5英尺高的浪头,约以每小时8-9英里的速度前进。后来它的高度逐渐减低,经过约一英里(1英里1.609千米)的追逐后,在运河的拐弯处消失了”。为了探究上述的水波鼓包到底是一种什么样的现象,随后,罗素在水槽的一端用一重锤垂落入水中,对重锤激起的水浪的运动情况进行了反复的观察,如图5-1所示。他发现这种水浪与运河中出现的奇特水波是本相同。通过实验,他还总结出水波的移动速度v、水的深度d及水波幅度A之间的关系:B为一某比例常数。这实验结果说明,水波的运动速度与波幅的高度有关,波幅高的速度较快,且波幅的宽
3、度对高度之比也相对较窄。图5-1 罗素在浅水槽中做的水波实验然而罗素当年未能从流体力学出发给孤立波以合理的理论解释,因此没有引起人们的充分重视。直到半个世纪以后,即1895年,两位荷兰科学家科特维格(Kortweg)与德弗雷斯(de Vries)才对浅水槽中单向运动的奇特波动现象用一波动方程进行理论分析,得到了比较满意的解释。他们认为,这种现象是波动过程中非线性效应与色散现象互相平衡的结果。他们建立了以他们姓名的首写字母命名的方程,即KdV方程。KdV方程的形式如下:式中,u为相对于静止水面的高度,即波幅。KdV方程是一个非线性的偏微分方程,求解很难。为了求解KdV方程的孤立波解,可以只找与相
4、关的行波解,这里v为波速,于是偏微分方程化为的一个常微分方程,得解:,其图象与观察到的孤立波形状相同。但是此后孤立波现象的研究与KdV方程又被默默地遗忘了几十年。掀起这一邻域研究热潮的应归功于乌莱姆(S.Ulam)。1955年,在乌莱姆领导的美国阿尔莫斯国家实验室,著名物理学家费米(E.Fermi)、帕斯塔(J.Pasta)和乌莱姆数值计算了用非线性弹簧联结的64个质点组成的弦的振动,目的是从数值实验上验证统计力学中的能量均分定理。他们对少数质点进行激发,按照能量均分原理,由于弱的非线性相互作用,经长时间以后,初始的激发能量应有涨落地均衡的分布到每个质点。然而计算结果令人意外,长时间以后能量几
5、乎全部回到了初始集中在少数质点上的状态。这个结果预示着这个非线性系统可以出现孤立波。这就是著名的FPU问题。1965年,美国数学家采布斯基(Zabusky)与克鲁思卡尔(Kruskal),把FPU的非线性振子系统的能量不均分问题与KdV方程联系了起来。此后人们发现,在许多物理体系中都存在KdV方程,说明孤立波是一种普遍存在的物理现象。于是KdV方程被看作为数学物理的一个基本方程。此后人们又进一步发现,除KdV方程外,其它的一些偏微分方程也有孤立波解,从此一个广大的孤立波研究领域展开来了。由此可见,孤立波是既是一种特殊的,又是不难见到的波动现象,它被称为自然界里的相干结构(coherent st
6、ructure,或称拟序结构),即一种有序结构。从美丽的木星上的巨型红斑到固体中的电荷密度波都属于这样的有序结构。从运动形态上讲,相干结构与混沌运动既是相互对立的,但又是具有内在联系的两种非线性现象。混沌运动所呈现的是非线性中奇妙的无序状态,相干结构则反映了非线性系统中的惊人有序性。从尺度上讲,相干结构可以小到原于尺寸,大到天文范围。木星上巨型红斑的尺寸大达4108米,相当于地球与月亮之间的距离,图5-2就是上世纪70年代由旅行者1号拍摄到的一幅木星上的巨型红斑照片。七十年代,从“阿波罗联盟号”宇宙飞船上拍摄到的出现在泰国安达曼海面上的孤立波照片,照片上显示出五条几乎平行的直线表面波的波包,每
7、一条大约150公里宽,两个波包相互间大约相距10公里;罗素所观察到的水面上孤立水波的尺寸在1米量级,而在二硫化钽晶体中所呈现的电菏密度孤立波,其尺度仅约米。图5-2 木星上的巨型红斑照片2. 孤立波与孤立子遵循着FPU问题的思路,采布斯基与克鲁思卡尔还采用数值模拟的方法,用计算计机又计算了两个具有不同速度孤立波前后追逐中发生的现象。设有同向行进两个孤立波,波幅较高在后的孤立波,逐渐赶上前面幅度较低的孤立波,于是两个孤立波的相遇了。令人惊奇的是两个孤立波相遇后,又能很好地分离开来继续前进,而原来的波包形状没有发生大的变化。图5-3给出了两个计算出的孤立子随时间演化的情况。图( i )是时间t =
8、 0.5,两个高低不等同向行进的孤立子,高的在后,但速度较快;图( ii )是时间t = 0.1,幅度较高的孤立子逐渐赶上幅度较低的孤立子;图( iii )是时间t = 0,两个孤立子到达同一位置;图( iv )是时间t = 0.1,幅度较大的孤立子开始超过幅度较小的孤立子;图( v )是时间t = 0.5,幅度较大的已超过幅度较小的,于是两个孤立子分离开来了。图5-3 两个KdV孤立子的碰撞过程这个计算结果说明两个孤立波的相遇具有类似于物质粒子之间的碰撞特性,称为孤立波碰撞。说明孤立波具有非常独特的稳定性。为了强调孤立波的这种碰撞特性,采布斯基和克鲁思卡尔将具有碰撞特性孤立波称为“孤立子”或
9、称“孤子”(soliton)。事实上,并非所有的孤立波都具有这种稳定的碰撞特性,人们只把具有这种碰撞稳定性的孤立波称之为“孤立子”。从物理本本质上讲,孤立子是由非线性场所激发的、能量不弥散的、形态上稳定的准粒子。这种准粒子具有一切粒予所具有的特性,如能量、动量、质量、电荷、自旋等等,它门也遵循一般的自然规律,如能量、动量、质量守恒定律。另一方面,它又有自身的特征波动性,在一切可以出现波动的介质里,在一定条件下都可存在。除上面介绍的浅水层外,在水层深处、固体介质、电磁场、等离子体、生物体、以及微观粒子的波动性中都可能有孤立波存在。它是一种行波,它既可以以速度v在空间传播。又可以处于静止状态,成为
10、非传播的孤粒子。与量子力学所描述的微观粒子相比,孤立子遵循经典运动规律,服从牛顿运动方程或哈密顿运动方程。所以孤立子是一种新型的准粒子,它是本世纪物理学中提出的一个重要的新概念。由于对孤立子的研究具有十分重要的理论与应用价值,因此越来越受到各国的科学界的重视。目前,在理论上和实验上已对孤立波巳作了大量的研究。并巳有了明确的定义。通常在数学上把具有下列性质的非线性方程的解称为孤立波解:向单方向传播的行波;分布在空间的一个小区域中;波动形状不随时间演变而发生变化;孤立波之间的相互作用具有类似粒子一样的弹性碰撞。现在已经知道一系列非线性偏微分方程存在孤立波解,其中最有代表性的有四类方程:(1)KdV
11、方程(2)正弦高登(Sine-Gordon)方程(3)户田(M.Toda)非线性晶格方程(4)非线性薛定谔方程(NLSE)。至于孤立波的形状,象水面上的鼓包形状只是其中一种,除此以外还有几种另外的形状。图5-4给出了有代表性的四种类型:( i )波包型,(ii)凹陷型,(iii)扭结型,(iv)反扭结型。( i )、(ii)两种是在,;(iii)、(iv)两种是在,趋近于不同的数值。图5-4 四种孤立子类型第二节 KdV方程1. 1. 波动中的非线性会聚效应其实,波动中的会聚现象在生活中并非罕见。微风吹拂,水面只掀起层层碎浪;劲风吹来,浪尖则卷起浪花。这就是一种非线性会聚效应。这是因为在劲风的
12、推动下,在水浪的不同高度上有不同的前进速度。同样的情况可以出现在海滩边。远处传来的海浪越近海岸,浪头越高,终于在离海岸不远处卷起了浪花。这是因为海滩对水浪运动产生某种阻滞力,海浪的较低部分受到的阻滞力较大,较高的部分阻滞力较小。由此可见,当一个水浪的不同部分有不同的行进速度时,将会出现会聚效应。特别是当水浪高处前进速度大,低处前进速度小,水浪会在前进中越来越前倾,于是在某一时刻波前出现坍塌,卷起了浪花,如图5-5所示。图5-5 水浪在行进中波前逐渐变陡,最终波形出现坍塌并卷起了浪花现在用一个简单的数学模型来说明波动中的会聚现象。已知介质中的波动是一种随时间传播的扰动。设波动介质是由相互间没有作
13、用的粒子组成的。为简单起见我们研究一维情况,设在时刻t,介质中x处的粒子密度为。由于粒子既不会产生,也不会消灭,所以有: (5-2-1)将式(5-2-1)写成全微分: (5-2-2)由于为粒子的移动速度,因此式(5-2-2)可以写为: (5-2-3)在一般情况下速度v是坐标与时间的复杂函数。然而,如果速度v = v0是常数,这是线性情况,则式(5-2-3)具有行波解:这时介质的移动速度v0也就是波速。在初始时刻介质中出现的扰动,即波动 (5-2-4)将在传播中保持不变。波动将以速度v0无畸变地沿x方向前进。然而在非线性的情况下,式(5-2-3)的解将是很复杂的。例如,波动的速度v可能与介质的密
14、度n有关,即:。在这种情况下,一次近似具有下面形式: (5-2-5)由此可见,方程(5-2-5所描写的行波,波的各部分将具有不同的运动速度,特别是当时,波速随密度增大而增大。因此可以想象,随着波的传播,波包的前沿会越来越陡,意味着形成了某种会聚效应。随着行进中波包前沿越来越前倾,于是在某个时刻会发生类似于图5-4所示的波形坍塌现象。2. 波动中的色散波动在传播中往往存在色散现象。我们知道,一束平面单色波与一列正弦波相对应。一个频率为沿x方向传播的平面波可以表示为: (5-2-6)式中,为波长,在多维的情况下应写成矢量k,矢量指向代表波的传播方向,故称为波矢,也称传播常数。平面波(5-2-6)的
15、等相位面由下式给出:由得等相位面的运动速度,即相速v: (5-2-7)它代表了一列平面波的传播速度。一个线性微分方程 (5-2-8)它具有平面波解。将式(5-2-6)代入方程(5-2-8)可得 (5-2-9)由此得:或 (5-2-10)因为m是可以任意设定的常数,所以满足(5-2-8)方程的平面波解有许多个。根据波的叠加原理,方程(5-2-8)的通解,可以把具有不同波矢的许多平面波(谐波)叠加起来构成。式(5-2-10)称色散关系,显然常数,说明了不同k值(从而值也不同)的平面波有不同的相速。一个波动(或称波包)的移动速度是由群速度决定的。由式(5-2-10)得波的群速度: (5-2-11)因
16、此,具有不同k值的波(群),将具有不同的群速度。不同波数的各个子平面波以不同的群速度传播,所以方程(5-2-8)通解所描述的波动,在它运动时将改变它的形状并弥散开来。于是,初始时刻出现的波包,会随时间的推移而发生变化,波包发生弥散,以至在某个时刻波包完全消失。因此,如果一个波动的所有谐波都以同一的速度行进,=常数,就是非色散波;反之,如果每个谐波都有不同的行进速度,常数,就是色散波。附带说一下,线性波动方程(5-2-8)与色散关系(5-2-9)之间存在着对应关系。实际上对于平面波解(5-2-6),我们发现有下述的对应于是便可以在波动方程与色散关系之间建立直接的对应。知道了这样的对应关系,我们可
17、以有色散关系直接构造出方程来,反之亦然。例如,如果有色散关系则波动方程为:3. KdV方程如上所述,一个线性波动由于在介质中传播时既存在色散,所以该波动是不种稳定的。只有当在波动中存在非线性的会聚时,如果色散与会聚两种作用出现某种平衡,才会出现波形稳定的孤立波。在KdV方程中正是同时存在了这两种效应。KdV方程的推导过程是很复杂的,全面地介绍过于冗长,这里采用一种简化模型以了解一下KdV方程的推导思路。首先,对于不可压缩介质,式(5-2-3)所表示的粒子随时间与坐标变化关系中,粒子数密度n应可以用粒子速度v来替代,即有 (5-2-13)现在考虑水面波动的色散关系,可以证明,忽略表面张力,在重力
18、作用下水波的色散关系为: (5-2-14)式中g为重力加速度,h为水的深度。对式(5-2-14)进行级数展开,并略去高次项后有: (5-2-15)式中,与色散关系(5-2-15)对应的方程为 (5-2-16)现在将导致波形坍塌的非线性效应的方程(5-2-13)与色散效应的方程包含到一个方程中,就得到KdV方程 (5-2-17)4. KdV方程的孤立波解设方程(5-2-17)有如下形式的解: (5-2-18)将式(5-2-18)代入式(5-2-17),得到以为时间变量的方程:整理后 (5-2-19)这里,。将方程(5-2-19)对进行一次积分,得:用微分的简写形式来表示,即对一次微分写为u,二次
19、微分写为u” 等等。对上式乘以u,再进行一次积分得 (5-2-20)积分常数可以看成为体系的哈密顿量。方程(5-2-20)进一步可以写成为 (5-2-21)其中第一项为动能,第二项为势能 (5-2-22)是一三次曲线,如图5-6所示,一般有三个零点,即,和,且,它们可用a,和C来表示。根据三次代数方程的解与系数的关系,有:=3a。设于是将式(5-2-21)可改写为: (5-2-23)方程(5-2-23)可写成下述积分 (5-2-24)是待定的积分常数。图5-6 曲线为了保证有实数解,要求(5-2-24)积分号中的根号内为正值。在区间内为负;在区间或内随u是单调变化的,也就是无界的;无界,u将随
20、(或)而趋于无界,因此我们只能考虑区间内的积分。引入代换:式(5-2-24)变为 (5-2-25)式中,为椭圆函数的模数。设所以 (5-2-26)这是KdV方程的椭圆余弦波解。波的周期为所以的周期为这样是空间周期为L的一列行波,如图5-7所示。由于它不具有局域性质,所以它不是孤立波。图5-7 KdV方程的椭圆余弦波解当椭圆函数的模数,即时(图5-8a),椭圆余弦函数演变为三角余弦函数,于是解(5-2-26)化为 (5-2-27)这是振幅十分小的余弦波解,它是由于幅度很小的非线性项可以忽略的结果。当椭圆函数的模数,即时(图5-8b),椭圆余弦函数演变为双曲正割函数,于是解(5-2-26)化为 (
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 第五 孤立 27
限制150内