电脑桌面
添加小米粒文库到电脑桌面
安装后可以在桌面快捷访问

声波加pml边界条件VIP免费

声波加pml边界条件_第1页
1/7
声波加pml边界条件_第2页
2/7
声波加pml边界条件_第3页
3/7
雷克子波 因计算机资源有限,不便在太大的区域做PML 的计算,一般只在边界上取有限宽度的区域作为PML 计算区域。根据相关文献的研究,PML 边界区域最少长度应为半个波长[6]。本文综合考虑了效果与开销等因素,选取了边界上 50 层作为 PML 的计算区域。常规计算区域与 PML 边界区域的如图 2-4 所示。 xxxxxPvaPKtx zzzzzPvaPKtz xzxxvvPaPKtxz 衰减系数 2331lo g2xzPvRL 式中 L 为 PML 层的厚度, 为层内的点距 PML 与非 PML 的边界的距离,Pv 为纵波速度, R 那么在 PML 边界区域内,对于式(2-13)1xxxvPvtx 即为理论反射系数,一般取 0.001 较为合适, 为 方向的空间步长。 1xxxvPvtx,可看作为在常规的计算方程基础上,减去一项进行 PML 的阻尼修正项。因本文中只考虑各项同性介质中的地震波传播规律,故可做xz假设。在此利用一下三个假设: 2121,2121,21,21kjixkjixkjixvvv 21,2121,21,2121kjizkjizkjizvvv 1,1,21,21kjikjikjiPPP 因为以上三个近似精度均为时间方向上的近似,且时间精度均为二阶精度,因交错网格技术的时间精度为二阶,故以上近似不影响本式的计算精度。 故可得: NmkmjikmjimkjikjixkjixPPaxtvttv11,,21,2121,2121,15.015.011 (2-18a) 同理: NmkjmikjmimkjikjizkjizPPaztvttv1,1,21,21,2121,2115.015.011 (2-18b) NmkmjizkmjixmNmkjmizkjmizmkjikjikjivvaxtvvaztKPttP12121,2121,121,2121,2121,,1,5.015.011 (2-18c) 此为在PML边界区域内的弹性声波应力-速度方程组。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%交错网格---非均匀介质二维声波方程(一阶压力--速度)、2阶时间差分、2阶空间差分精度 %%加上边界 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% close all;clea...

1、当您付费下载文档后,您只拥有了使用权限,并不意味着购买了版权,文档只能用于自身使用,不得用于其他商业用途(如 [转卖]进行直接盈利或[编辑后售卖]进行间接盈利)。
2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。
3、如文档内容存在违规,或者侵犯商业秘密、侵犯著作权等,请点击“违规举报”。

碎片内容

声波加pml边界条件

确认删除?
VIP
微信客服
  • 扫码咨询
会员Q群
  • 会员专属群点击这里加入QQ群
客服邮箱
回到顶部