基于波动方程中非结构化网络研究热点.pdfVIP

基于波动方程中非结构化网络研究热点

基于波动方程中非结构化网络研究热点

导读：数值模拟作为波动理论一个重要的研究手段得到了许多研究

人员的重视,采用有限差分方法对控制方程进行离散，计算区域的划分采

用矩形结构化网格。由本站硕士论文中心整理。

第1章绪论

课题研究背景

波的传播是经典的物理问题。随着科学技术的不断发展，波动理论的研究

在各类工程技术、军事技术等领域的实际问题中都有着重要的意义。如在

轮机工程中存在大量的振动问题，不仅会使动力装置的精度和其他性能降

低，同日寸会降低系统的隐身性能。例如，舰船轴系振动引起推进轴断裂

、汽轮机组强烈振动和管路系的振动噪声等等。因此，以波动理论为基础

，了解其振动传播的基本规律，从Ifu寻求实际解决的办法是很重要的一

个方面。十是求解波动方程的数值方法研究也越来越得到重视。

波动方程是重要的偏微分方程，它可以表述所有种类的波，例如声波、光

波、水波等等。通常，这些偏微分方程很难用解析方法求解，尤其在复杂

的非均匀介质中。因此必须借助数值计算方法或其它渐近解。目前，求解

波动问题所采用的数值方法主要分为这样几类，有限兀方法(Finite

ElementsMethod，边界兀方法(BoundaryElement

Method)，伪谱法(pseudospectracmethod)、有限差分方法(Finite

DifferenceMethod)，近年来还出现了界十有限差分法不I

I有限兀法之间的有限体积方法(FiniteVolumeMethod).

其中有限差分方法是比较经典的一种方法，它是将波动方程中波场函数的

空间导数和时间导数用相应的空间、时间的差分代替。有限差分方法计算

速度快，占用内存小，能准确模拟波在各种介质中的传播规律，它以其简

单高效的优点得到了很多研究人员的青睐。

然传统的有限差分方法也存在着一定的局限性。例如，因为其计算域网格

往往是基十单一的矩形网格来进行划分的，当面临计算域有弯曲界面或起

伏表面时，就会使得其计算受到限制;Ifn又因为其计算是基十固定步长的

网格进行的，因此使得结果稳定性变差。面对这些缺点，很多研究人员便

开始了对有限差分方法进行优化改进的工作。对十改进差分方法的研究也

逐渐成为波动理论的数值方法研究中的一个热点问题.

1.2.1国内外发展概况

求解波动问题常用的几种数值方法

求解波动问题所采用的数值方法主要有这样几种，有限兀方法，边界兀方

法，伪谱法、有限差分方法，还有近年来出现的有限体积法。

有限兀方法的基础是变分原理和剖分插值。把边值问题转化为变分问题，

使其对应的欧拉方程等十待解的微分方程，然后经过剖分插值方法求解其

泛函的极值函数，解出待解的边值问题。通过灵活的网格划分，使得其能

处理多种介质和边界条件，模拟精度较高，但是有限兀方法的问题在十其

不适肩{大规模的模型计算，b_计算量大.

边界兀方法对时间变量进行一维傅氏变换，把波动问题转换为Helmholtz

方程边值问题，再利用格林公式将新的边值问题转化为积分方程，最后用

边界单兀法求解此积分方程。与有限兀法相比，边界兀法把所讨论问题的

维数降低了一维，所用节点数比有限兀法少，但是边界兀法得到的代数方

程组所对应的矩阵一般是非稀疏的，有限兀法得到的代数方程组对应的矩

阵一般是稀疏的，因此对同一问题的处理所需的存储量和计算量未必比有

限兀法少.

伪谱法也称虚谱法。是在20世纪七十年代引入数值模拟计算领域的，它是

将波动方程中的空间微分变换成频域中的乘法运算，只需要较少的空间格

点就可以得到较高的计算精度。可以看成是一种无限阶的有限差分法，是

传统有限差分法的一个推广。伪谱法在粗网格上也能实现高精度的计算，

相对有限兀法实现起来较容易.

有限差分法，是一种最常用的数值模拟方法。它是对波动方程中的微分用

差分近似代替，把波动理论上的微分形式变成适用十数值计算的差分形式

，得到差分方程。在数值模拟计算时，给定初始条件，通过差分方程即可

由空间中过去时刻的波场值的分布和变化推导出当前时刻的波场值分布。

有限差分法有两种算法:显式算法和隐式算法。有限差分方法在模拟波传

播的计算中有两个主要的优点:(1)能够模拟介质的完全响应，描述波在介

质中传播产生的反射、透射和散射等所有振相，其精度仅受网格密度和网

格区域大小的影响;;(2)网格划分灵活、简单。与其它方法相比，更适用

十非均匀复杂介质模型的波场模拟.

有限体积法是上世纪60年代被引入数值计算领域的。但直到80年代