《海洋数值模拟》实验2指导书解析.docxVIP

PAGE PAGE # 目录 微分方程数值实验 错 微分方程数值实验 错.. 误！未定义书签 TOC \o 1-5 \h \z \o Current Document 1. 孤立波简介 2 \o Current Document 2. 有限差分格式 4 \o Current Document 蛙跳格式 5 \o Current Document 跳点格式 7 分裂格式 8 \o Current Document 3. 实例 9 4.附录：MATLAB源码 12 微分方程数值实验 KDV方程是历史上最著名的方程之一，其起源是水波中的孤立波。 孤立波简介 我们自己先设想一下：在一条窄河道中，当迅速前进的船突然停下时，在船 头形成的一个孤立的水波迅速离开船头，继续向前运动，如下图(1)所示。 图1.孤立波的运动图 孤立波是由是苏格兰一位优秀的造船工程师 Russell(John Scott Russell 1808-1882)在1834年研究船舶在运动中所受到的阻力时发现的。 在一次学术报告中，他是这样描述的： “我把注意力集中在船舶给予流体的运动上，立刻就观察到一个非同寻常而又非 常绚丽的现象，它是如此之重要，以致我将首先详细描述它所表现出来的外貌。 当我正在观察一只高速运动的船舶， 让它突然停止时，在船舶周围所形成的小波 浪中，一个紊乱的扰动现象吸引了我的注意。在船身长度的中部附近， 许多水聚 集在一起，形成一个廓线很清楚的水堆，最后还出现一尖峰，并以相当高的速度 开始向前运动，到船头后，继续保持它的形状不变，在静止流体的表面上，完全 孤立地向前运动，成为一孤立行进波，直到河道的转弯处才开始消失掉。” 纵观力学和物理学的发展史不难看到，每当开始引入一种新思想或新概念的 时候，往往会受到怀疑和非难，并常会引起激烈的争论，孤立波的命运也是如此。 当时科学界的权威们对Russell的这些结果，一开始时就表示了怀疑和反对。 甚至连英国流体力学家斯托克斯 (George Gabriel Stokes,1819-1903爵士也对此提 出质疑，怀疑在静止水面上能存在不变形的行波。他们的怀疑的问题主要有： 一 个完整的波动为什么会全部在水面上，而不是一部分在水面上，一部分在水面下； 波在传播的过程中，为什么波幅不会衰减；波的运动速度也与他们的研究结果不 符。 这一争论延续到19世纪70年代才初步得到解决。H.E.巴津(Bazin,H.E.)和英 国科学家瑞利(John William Strut Rayleigh，1842-1919)对孤立波进行了一系列的 研究，证明了 Russell的工作是正确的。 Russell与艾里，斯托克斯的争论，最终于 1895年由数学家D.J.科尔特弗 (Korteweg ，D.J.)和他的学生 G德.弗里斯(Vires ,Gde)所解决。他们在小振幅与 长波的假定下，从流体动力学导出了关于孤立波的方程 (后人称它为KdV方程)。 这一方程的行波解，在波长趋于无限的情况下，正是 Russell所发现的孤立波。 KdV方程的提出，从理论上阐明了孤立波的存在，给这场争论划上了句号。 从Russell的发现到KdV方程的提出，大约经历了 60年时间,孤立波才为学 术界普遍接受。 数学家Korteweg ,D J和G.de Vires从流体动力学中导出了关于孤立波的方程 (后人称它为KdV方程) 3 旦，启十0， -:t fx ：：x3 第二项空间采用中心差分 蛙跳大部分时间 (1)在数学物理方程中，可解得其精确解一行波解［9］，如下 u(x,t) =3Csech2 Ax - Bt D 1， (2) 其中C和D是常数；A = 1 IC,B八AC。 2 ■■；: 下面从物理学和数学角度理解孤立波的特点， 物理学：孤立波是物质非线性效应的一种特殊产物。 KdV方程中含有非线性项 -u —。 :X 数学：某些非线性偏微分方程的一类稳定的，能量有限的不弥散解。 孤立波可谓是孤立，如下 在传播过程中保持波形和速度不变，像是“透明地”穿过对方。 碰撞时两个孤立波重叠在一起，其高度低于碰撞前孤立波高度较高的一个(这 表明在非线性过程中，不存在线性叠加原理)两个孤立波碰撞时互相穿透且维 持原来的波形和速度；孤立波的波幅愈高，其传播速度愈快，如下图 (2)。 (a)(M(d) (a) (M (d) 图 2?两孤立波的传播的交会情形，(a)t = = ； (b)t =t2 t1 ； (a)t = t3 t2 ; (a)t = t3 t4。 碰撞后孤立波的轨道与碰撞前有些偏离(即发生了相移)。 在海洋中也存在孤立波，如下 涌浪向缓慢倾斜的海岸传播时，随着涌浪接近海岸，波峰变尖，到水