24小时热门版块排行榜    

查看: 1705  |  回复: 9
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

Ego悄吟

新虫 (初入文坛)

[求助] 基于有限差分的波场模拟,pml边界,震源为单力源已有2人参与

本人最近在做波场模拟的程序设计,在波场分离的时候遇到下图问题,震源附近发生突变,有没有大佬能够指出问题所在?
模型大小500x500x500,vp为3200km/s,vs为vp/1.732,密度为2250kg/m^3,震源位于中心,采用雷克子波,中心频率为50hz。

基于有限差分的波场模拟,pml边界,震源为单力源
水平.jpg


基于有限差分的波场模拟,pml边界,震源为单力源-1
3D.jpg
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Ego悄吟

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by diwushuxue at 2019-04-17 11:39:28
网格多大,你这跟波场分离没关系,是你正演没做好。有可能是频散过于严重

您好,看到您的回帖之后,我回去划细网格了,之前网格大小为10m,计算步长为1.60e-3s,之后网格大小为2.5m,计算步长为0.8e-3s,可是依然出现这个情况。
基于有限差分的波场模拟,pml边界,震源为单力源-2
3D.jpg


基于有限差分的波场模拟,pml边界,震源为单力源-3
水平.jpg

3楼2019-04-18 09:21:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Ego悄吟

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by diwushuxue at 2019-04-18 17:34:15
还是正演问题,你空间网格10m,时间采样1ms,主频30Hz。再试一下!...

您好,根据您给的建议,我回去试了试,结果还是出现之前的问题。
基于有限差分的波场模拟,pml边界,震源为单力源-4
3D.jpg


基于有限差分的波场模拟,pml边界,震源为单力源-5
水平.jpg

5楼2019-04-19 09:57:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Ego悄吟

新虫 (初入文坛)

引用回帖:
6楼: Originally posted by baobiao007 at 2019-04-19 19:21:44
明显是正演的递推过程没写对。。

您好,您的意思是差分出错了么?我检查了我的程序,貌似并没有问题。

发自小木虫IOS客户端
7楼2019-04-20 08:39:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Ego悄吟

新虫 (初入文坛)

引用回帖:
8楼: Originally posted by baobiao007 at 2019-04-20 10:36:59
送你段例子代码吧,直接运行即可


% Finite-differences in time domain(FDTD) acoustic wave propagation in 2D
% medium

close all;
% Output every ... time steps
IT_DISPLAY = 10;

%% MODEL
% M ...

您好,谢谢您的代码,我的模型的问题确实出现在差分上。

发自小木虫IOS客户端
10楼2019-04-21 09:28:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Ego悄吟 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 双非院士组还是211普通老师 +7 timergoi 2024-05-14 8/400 2024-05-16 01:18 by timergoi
[教师之家] 加上“青年”两个字,意义就变了 +10 zylfront 2024-05-13 13/650 2024-05-15 23:36 by flasheagle
[考博] 有博导缺学生吗,我可以当牛做马 +26 刘秋彤 2024-05-13 44/2200 2024-05-15 22:58 by 驴肉火没驴
[基金申请] 面上项目代表作必须有通讯作者文章吗? +27 zlxnnx111 2024-05-09 46/2300 2024-05-15 22:33 by 星火12
[基金申请] 南林海归博士首聘期 没拿到国自然,被降级后自杀 +30 babu2015 2024-05-13 32/1600 2024-05-15 22:11 by Howard28
[考博] 深度学习申博 +3 酸菜鱼. 2024-05-15 3/150 2024-05-15 21:29 by Youjun02
[基金申请] 河南工业大学过年开学后,就不让用国基金报账,这合法合规吗? 100+8 知己而已 2024-05-09 38/1900 2024-05-15 21:26 by msjy
[论文投稿] 当前scientific reports还值得投稿么? +4 lizhengke06 2024-05-14 5/250 2024-05-15 20:16 by vmuch520
[教师之家] 宋凯 +9 yexuqing 2024-05-14 9/450 2024-05-15 19:22 by otani
[教师之家] 学生家长私下联系老师修改成绩不成,唆使19名学生联名要求复核成绩 +22 sjtu2012 2024-05-11 25/1250 2024-05-15 16:45 by newfuzzy1
[教师之家] 另一半是老师会加分哦 +7 acogw9004 2024-05-14 7/350 2024-05-15 16:09 by feng6531
[基金申请] 基金委治打招呼顽疾越治越严重 +36 zzahkj 2024-05-10 65/3250 2024-05-15 16:06 by eseayl
[基金申请] 评审规则突发奇想 +18 平凡冰雪花 2024-05-13 23/1150 2024-05-15 15:54 by 平凡冰雪花
[职场人生] 硕士毕业四年工作不顺心 +4 noora2018 2024-05-14 4/200 2024-05-15 15:19 by 1158057902
[论文投稿] 投稿竟然倒在了Scientific Reports上。 +12 lizhengke06 2024-05-10 14/700 2024-05-15 09:34 by zqwcr
[基金申请] 2024工材口的国自然开始涵评了吗? +3 Thanks1989 2024-05-13 3/150 2024-05-14 22:27 by weiyongscut
[教师之家] 今天才知道硕士考进大专比博士进本科工资高得多 +3 ZHONGWU_U 2024-05-14 4/200 2024-05-14 18:40 by 起个新名字好难
[论文投稿] 中文投稿 +5 liujinyu80 2024-05-13 5/250 2024-05-14 17:13 by icm639
[基金申请] 再谈一作和二作 +19 jklily 2024-05-10 33/1650 2024-05-13 17:12 by mierbushiwo
[考博] 24/25申博,科研助理找导师 +9 cc毕业了吗 2024-05-12 9/450 2024-05-13 14:31 by zxl_1105
信息提示
请填处理意见