研究磁性物理小论文范文.docVIP

　　研究磁性物理小论文范文 研究磁性物理小论文范文 导读：X格；有限差分；完全匹配层1引言现在永久磁体在许多技术中都是重要的组成部分。经过近40年的研究,稀土永磁体经历了第一代稀土永磁体(SmCo5)、第二代稀土永磁体(Sm2(Co,Fe)17)和第三代稀土永磁体(NdFeB)。目前,第三代永磁体(NdFeB)的磁能积已经接近其理论极限值。现代高技术(信息与通讯技术、航空航天技术、磁 基于混合X格的有限差分地震波场正演模拟 摘要：为克服单一的交错X格或者规则X格在波场计算中存在的不足，本文提出了一种基于混合X格的时域有限差分算法。 该算法在有效计算区域内采用四阶中心差分算子的规则X 格有限差分法， 与混合X格方法相比较具有较高的计算效率并能节省存储空间； 在边界处理 中则引入基于交错X格的完全匹配层吸收边界条件， 能够有效地抑制波场模拟正演算法中存 在的频散问题以及人工边界造成的反射干扰。 在理论分析的基础上， 我们解决了规则X格与 交错X格计算区域边界X格匹配与变量计算的问题， 给出了算法实现的详细步骤， 最后通过 Marmousi 模型等模拟算例进行了正演波场模拟，从模拟结果中验证了该算法的有效性。 关键字：声波方程；混合X格；有限差分；完全匹配层 1 引言现在永久磁体在许多技术中都是重要的组成部分。 经过近 40 年的研究,稀土永磁体经历了第一代稀土永磁体(SmCo5) 、 第二代稀土永磁体(Sm2(Co,Fe)17)和第三代稀土永 磁体(NdFeB) 。 目前,第三代永磁体(NdFeB)的磁能积已经接近其理论极限值。 现代高技术(信 息与通讯技术、航空航天技术、磁悬浮交通技术、核磁共振技术、纳米技术等)的发展对永 磁体的性能与质量提出了更高的要求。 2 原理地震波的传播过程是一个复杂的问题。 波动方程的组成包括波动函数相对于 空间坐标和时间坐标的偏导数计算。 虽然弹性 动方程包含的信息量更充分， 使用弹性波方程计算的数据更精确，但是这种方法的计算量非常大，计算消耗昂 贵。 并且决定场数据的主要因素是 P 波，可以近似的不考虑 S 波的影响.因此在 实际的应用中， 往往应用声波方程来近似的描述地震波的传播规律目前数值模拟 中常用的波动方程可分为二阶波动方程和一阶位移-应力波动方程。 2.1 规则X格有限差分法 无损流体媒质中声波的二阶偏微分支配方程如下[18]： 2 p ( x, z , t ) 2 p ( x, z , t ) 1 2 p ( x, z , t ) x 2 z 2 c ( x, z ) 2 t 2 (1) 其中 p( x, z, t ) 表示在某一时刻 t 上计算区域内二维空间任一点 ( x, z ) 处的位 移， c( x, z) 表示纵波的传播速度。 对方程(1)进行有限差分离散化。 2.1.1 波动方程二阶差分近似 先取二阶差分精度格式，也即忽略小量 x 和 z 的二阶及其高阶指数项。 从 而得到离散化的声波方程的二阶差分格式如下：p n 1 (i, j ) 2 p n (i, j ) p n 1 (i, j ) c 2 t 2 p n (i 1, j ) 2 p n (i, j ) p n (i 1, j ) x 2 (2) c 2 t 2 p n (i, 3 4 5 6 7 8 9 10 研究磁性物理小论文范文 导读： j 1) 2 p n (i, j ) p n (i, j 1) 2 z 2.1.2 波动方程四阶差分近似 由于有限差分法的缺点是数值频散问题， 因离散化求解波动方程而产生的伪 波动， 是差分方程所固有的本质特征。 有限差分法中的数值频散包括空间频散和 时间频散两部分，而空间频散是主要方面 [19]。 且数值分析表明[20-21]，X格频散 随X格间距的增大而加剧， 而减小X格间距必然会增大计算量。 但只要空间X格 间距和时间X格间距均满足采样定理，不产生时间和空间假频，随着差分方程阶 数的提高，对物性参数改变就越小，进而数值频散越少。 所以，可以通过提高差 分方程阶数的方法来提高波场模拟精度，压制数值频散。 本文用四阶差分逼近声 动方程(1)式中的空间偏导数，时间仍采用二阶 差分精度格式。 从而得到声波方程的四阶差分精度离散化结果如下：p n 1 (i, j ) 2 p n (i, j ) p n 1 (i, j ) c 2 t 2 p n (i 2, j ) 16 p n (i 1, j ) 30 p n (i, j ) 16 p n (i 1, j ) p n (i 2, j ) 2 12x c 2 t 2 p n (i, j 2) 16 p n (i, j