简介:针对矩形网射线追踪存在的模型剖分灵活性差、速度界面描述精度差等问题,研究了复杂结构三角网最小走时射线追踪全局算法。(1)根据剖分区域点、线、面的结构关系,遵循Delaunay三角剖分的优化准则进行三角网格剖分;(2)定义三角单元射线追踪的拓扑关系;(3)波源点及某一时刻波到达的每一个节点点构成波行面,在波行面扩展过程中计算节点的最小走时和次级源位置,实际次级源检索采用双曲线近似算法;(4)利用各节点走时和次级源方向信息,通过最小走时搜索,拾取从接收点到源点的射线路径。数值模拟结果表明,三角网射线追踪方法模型剖分时灵活性强、速度间断面的描述精度高,追踪结果准确。
简介:有限差分方法广泛应用于求解许多科技领域所涉及的偏微分方程,高阶显式有限差分方法通常用来提高求解精度,已经提出的高阶隐式有限差分方法和截断高阶显式有限差分方法可用来进一步提高模拟精度而不增加计算量。本文首先计算了针对常规网格上的一阶导数和二阶导数、交错网格上的一阶导数的有限差分系数,发现高阶隐式有限差分系数中存在一些小的系数。频散分析结果表明:忽略这些小的差分系数能够近似维持有限差分的精度,但是显著减小了计算量。然后,引入镜像对称边界条件来提高隐式有限差分方法的精度和稳定性,采用混合吸收边界条件来减小来自模型边界所不需要的反射。最后,给出了针对均匀和非均匀介质模型的弹性波模拟例子,表明了本文方法的优点。