24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

ycy66ds

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 229.9
帖子: 178
在线: 80.4小时
虫号: 2133185
注册: 2012-11-18
性别: MM
专业: 地震学

[求助] 谱元法波动方程数值模拟-波场快照

针对波动方程数值模拟，我采用的谱元法（spectral element method）
目前研究的是：二维区域，各向同性，自由边界条件，集中外力。
波场快照如附件所示
问题：
1.在P波和S波波前之间，因为集中力的原因，位移分量Ux,Uz应该为0，但是我所得的结果从附件中可以看出是非零的。
2.在单元与单元边界处，出现的位移分量有明显跳跃，而理论的波场快照肯定是连续的。
ps~望各位同仁不吝赐教

波场快照

回复此楼

» 猜你喜欢

1楼2013-08-27 16:29:34

已阅关注TA 给TA发消息送TA红花 TA的回帖

nahua

铁杆木虫 (正式写手)

钻石级木虫

应助: 2 (幼儿园)
金币: 16237.2
散金: 100
帖子: 859
在线: 715.6小时
虫号: 601151
注册: 2008-09-12
性别: GG
专业: 计算数学与科学工程计算

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

引用回帖:

11楼: Originally posted by ycy66ds at 2013-08-29 20:10:47
都是慢慢摸索的，程序出问题两个月了。差分有开始并行么？

...

谱元法搞自由边界比较简单，但是不知道并行效率怎么样。
我用有限差分做的并行，速度还可以。

你用的是弹性波动方程？波场是x方向位移分量而不是速度分量？
我用过高斯子波球面震源，感觉效果还可以。
调整一下图像的显示，只显示能量级别较高的分量看看效果。

赞一下(1人)

回复此楼

12楼2013-08-30 09:01:43

已阅关注TA 给TA发消息送TA红花 TA的回帖

nahua

铁杆木虫 (正式写手)

钻石级木虫

应助: 2 (幼儿园)
金币: 16237.2
散金: 100
帖子: 859
在线: 715.6小时
虫号: 601151
注册: 2008-09-12
性别: GG
专业: 计算数学与科学工程计算

引用回帖:

13楼: Originally posted by ycy66ds at 2013-09-01 20:34:20
国际上，谱元法的并行已经做得非常强大，现在已经可以模拟全球尺度的地震。
我求解的是弹性波动方程。波场大家给的是位移分量，波场快照貌似大家普遍给的是速度。
我给的震源时间函数是rick子波。你的球面震 ...

呵呵，谱元法的原理我不懂，所以就不讨论了。

波场快照我用的是速度。
球面震源也可以简化到二维问题中的。

我的意思不是让你改结果，只是调整图像的显示。例如波场的主振幅的绝对值为10e-6左右，则可以只显示绝对值比10e-7大的结果，而不显示例如10e-12量级左右的结果。

赞一下(1人)

回复此楼

14楼2013-09-01 22:10:48

已阅关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 ycy66ds 的主题更新

返回列表