24小时热门版块排行榜

南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知（长期有效）

返回列表

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

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的回帖

ycy66ds

铁虫 (小有名气)

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

引用回帖:

12楼: Originally posted by nahua at 2013-08-30 09:01:43
谱元法搞自由边界比较简单，但是不知道并行效率怎么样。
我用有限差分做的并行，速度还可以。

你用的是弹性波动方程？波场是x方向位移分量而不是速度分量？
我用过高斯子波球面震源，感觉效果还可以。
调整一 ...

国际上，谱元法的并行已经做得非常强大，现在已经可以模拟全球尺度的地震。
  我求解的是弹性波动方程。波场大家给的是位移分量，波场快照貌似大家普遍给的是速度。
  我给的震源时间函数是rick子波。你的球面震源是不是用来处理三维问题的？
  还有一点，我们能够对所得到的的图像进行额外人工的处理吗？，个人觉得，算出来的结果不能随便更改吧？