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

» 猜你喜欢

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

ycy66ds

铁虫 (小有名气)

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

1.收敛条件CFL，是满足的，所以应该不是数值频散；
2.对于有限元，谱元法，自由边界是自动满足的，不用参与计算，而有限差分则需特殊处理。
3.我给的波场快照是波已经传播到边界，波未传到边界的波场快照也有，这种现象仍然存在；

3楼2013-08-27 20:49:33

ycy66ds

铁虫 (小有名气)

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

问题补充：震源处，采用的是集中力，只是在某个时刻有作用力。雷克子波震源也会出现这种问题。

困扰2个月了，求大神关注帮忙

4楼2013-08-28 08:55:59

ycy66ds

铁虫 (小有名气)

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

熟悉SEM的冒个泡呗！！！！

5楼2013-08-29 10:05:36

ycy66ds

铁虫 (小有名气)

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

引用回帖:

6楼: Originally posted by baobiao007 at 2013-08-29 18:20:06
谱元法太高深了，膜拜

那你们平时用什么方法呢？

8楼2013-08-29 19:56:11

ycy66ds

铁虫 (小有名气)

引用回帖:

7楼: Originally posted by Smilelzx at 2013-08-29 19:26:13

9楼2013-08-29 19:56:38

ycy66ds

铁虫 (小有名气)

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

引用回帖:

10楼: Originally posted by baobiao007 at 2013-08-29 20:08:26
我就搞了点有限差分...

都是慢慢摸索的，程序出问题两个月了。差分有开始并行么？

11楼2013-08-29 20:10:47

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子波。你的球面震源是不是用来处理三维问题的？
  还有一点，我们能够对所得到的的图像进行额外人工的处理吗？，个人觉得，算出来的结果不能随便更改吧？

13楼2013-09-01 20:34:20

ycy66ds

铁虫 (小有名气)

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

引用回帖:

14楼: Originally posted by nahua at 2013-09-01 22:10:48
呵呵，谱元法的原理我不懂，所以就不讨论了。

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

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

球面震源是不是就是爆炸源，纯膨胀，没有剪切，只产生纵波，不产生横波？图形显示方面受教了，感谢~

15楼2013-09-02 12:02:17