【调剂】北京石油化工学院2024年16个专业接受调剂
查看: 1669  |  回复: 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悄吟 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 申博求助 +3 dskabdh 2024-04-24 8/400 2024-04-24 23:17 by dskabdh
[考研] 博士生招生!! +5 2458827884 2024-04-19 5/250 2024-04-24 22:28 by 不懂的风清
[基金申请] 国社科项目,你们学校都限额申报吗? +6 屡战屡败 2024-04-21 8/400 2024-04-24 22:04 by njuzmwu
[基金申请] 两类问题算是白选了~ +6 jurkat.1640 2024-04-23 10/500 2024-04-24 21:18 by lucky_my2010
[访问学者] CSC的访问学者申请,没有个评审意见,也不知道怎么改,还有必要申请吗 +4 flyingship 2024-04-20 4/200 2024-04-24 21:09 by 59038mute
[找工作] 普通院校药学硕士,做合成的,感觉找不到工作 +5 pom戴墨镜 2024-04-24 9/450 2024-04-24 20:19 by l419110028
[硕博家园] 聊天 +9 暮色恋伊人 2024-04-22 10/500 2024-04-24 19:34 by liubai000
[基金申请] 化学结清有情发出来了 +3 starboy7286 2024-04-20 3/150 2024-04-24 15:52 by 嘿,黑贝
[基金申请] 颜宁:基础研究应顶天立地,立地即从需求端作为问题的出发点,立项背景就是面向真实需求 +11 zju2000 2024-04-18 11/550 2024-04-24 14:55 by 1501311
[考研] 0854-0855调剂 +6 shangannum1 2024-04-21 9/450 2024-04-24 14:43 by z_longlong
[公派出国] 東京大学 应用化学专业柳田研究室 招收修士/博士生 +5 蝶梦123 2024-04-19 5/250 2024-04-23 22:40 by 90472jvvh
[硕博家园] +8 s150535912 2024-04-18 9/450 2024-04-23 15:48 by huixiong0627
[教师之家] 发现督导在听课,需要跟TA招呼示意吗?有同学看手机、课件有疏漏算大问题吗? +11 河西夜郎 2024-04-19 11/550 2024-04-23 15:01 by mylifecg
[电化学] 有偿收购废弃锂硫电池 +3 汤圆喵 2024-04-19 5/250 2024-04-23 10:41 by Happy C
[论文投稿] 一直找不到审稿人 +4 lizhengke06 2024-04-21 5/250 2024-04-22 23:37 by topedit
[考博] 25申博记录贴 +6 我属驴核动力驴 2024-04-18 12/600 2024-04-22 15:50 by 诩雾
[找工作] 化工专业 10+3 xiao太阳94 2024-04-18 7/350 2024-04-22 13:48 by mlnjut
[论文投稿] 研究光催化的,好中的三四区 20+3 sl.0117 2024-04-20 3/150 2024-04-22 09:53 by bobvan
[考研] 312求调剂 +3 Leroic000 2024-04-19 3/150 2024-04-20 18:19 by xinaishan
[基金申请] 国自然青年基金只有一篇能第一标注,其他都是第二标注能结题吗? 10+4 lice_1987 2024-04-18 6/300 2024-04-19 08:29 by kudofaye
信息提示
请填处理意见