《海洋数值模拟》实验2指导书.doc

PAGE PAGE 13 目 录 TOC \o 1-2 \h \z \u HYPERLINK \l _Toc278465087 微分方程数值实验 PAGEREF _Toc278465087 \h 2 HYPERLINK \l _Toc278465088 1. 孤立波简介 PAGEREF _Toc278465088 \h 2 HYPERLINK \l _Toc278465089 2. 有限差分格式 PAGEREF _Toc278465089 \h 4 HYPERLINK \l _Toc278465090 2.1 蛙跳格式 PAGEREF _Toc278465090 \h 5 HYPERLINK \l _Toc278465091 2.2 跳点格式 PAGEREF _Toc278465091 \h 7 HYPERLINK \l _Toc278465092 2.3 分裂格式 PAGEREF _Toc278465092 \h 8 HYPERLINK \l _Toc278465093 3. 实例 PAGEREF _Toc278465093 \h 9 HYPERLINK \l _Toc278465094 4. 附录：Matlab源码 PAGEREF _Toc278465094 \h 12 微分方程数值实验 KDV方程是历史上最著名的方程之一，其起源是水波中的孤立波。 1. 孤立波简介 我们自己先设想一下：在一条窄河道中，当迅速前进的船突然停下时，在船头形成的一个孤立的水波迅速离开船头,继续向前运动，如下图(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 ,G.de)所解决。他们在小振幅与长波的假定下，从流体动力学导出了关于孤立波的方程(后人称它为KdV方程)。这一方程的行波解，在波长趋于无限的情况下，正是Russell所发现的孤立波。KdV方程的提出，从理论上阐明了孤立波的存在，给这场争论划上了句号。 从Russell的发现到KdV方程的提出，大约经历了60年时间,孤立波才为学术界普遍接受。 数学家Korteweg ,D.J和G.de Vires从流体动力学中导出了关于孤立波的方程(后人称它为KdV方程) ， 第二项空间采用中心差分 蛙跳大部分时间 (1)在数学物理方程中，可解得其精确解—行波解[9]，如下 ， (2) 其中和是常数；,。 下面从物理学和数学角度理解孤立波的特点, 物理学：孤立波是物质非线性效应的一种特殊产物。KdV方程中含有非线性项。 数学：某些非线性偏微分方程的一类稳定的，能量有限的不弥散解。 孤立波可谓是孤立，如下 在传播过程中