印度洋海啸波浪高度的地震成因仿真研究

在中美贸易战僵持的技术背景下,今天就不用美国的软件来讲这个案例了,我们来试试法国的开源软件Gerris Flow Solver (以下简称Gerris)。Gerris由Stéphane Popinet创建,并由Institut Jean le Rond d'Alembert(该研究院属于巴黎索邦大学)提供支持。现在Gerris的原始开发人员现在已经将开发转移到了Basilisk,它可以完成Gerris可以做的大部分工作,甚至更多。

640.webp.jpg


软件主要特点如下:

- 求解时间依赖的不可压缩变密度Euler,Stokes或Navier-Stokes方程

- 求解线性和非线性浅水波方程- 基于流场特征的网格自适应加密

- 复杂模型的网格自动化生成

- 空间和时间二阶精度

- 不限数量的对流、扩散粒子追踪- 可灵活加入源项

- MPI并行支持,动态负载平衡,并行可视化

- 基于VOF方法的多相流界面捕捉- 准确的表面张力模型

- 多相电-磁流体动力学      
       


案例介绍
   

2004年印度洋海啸是由于印度-澳大利亚和印度尼西亚的安达曼板块边界发生大规模断层破裂(> 1000公里)造成的。该案例运用格里利等人的断层模型作为圣维南原理分析海啸的初始条件。图1a中的动画展示了波高的演变。追踪波前中采用了自适应方法(图1.b),地形的动态重建依据ETOPO1数据集。

=640.gif

(图1a)波高动画图,图中峰值大于2m小于 -2m;
       


download.gif

(图1b)自适应动画图,图中波前峰值为海拔0.8海里和-101海里。



图2展示了断层破裂后在超过10小时的时间达到的最大波浪高度。642.webp.jpg

图2:超过10小时后最大波浪高度云图(以1m为显示间隔):

(a)孟加拉弯,最小值(蓝色)0m,最大值(红色)5m。

(b) 苏门答腊北部和泰国附近的细节,最小值(蓝色)0m,最大值(红色)大于8m。


最后,图3给出了在印度洋特定位置的不同时间观测波高(使用潮汐测量仪)和模拟波高的对比图。

640.webp (2).jpg

图3:在不同潮汐表位置观测到的波高和模拟波高的对比。

备注:水平轴是断层破裂后的时间(以小时为单位)。

地震导致印度洋海啸的波浪高度仿真的图8

640.webp (1).jpg

图4:观测波高(Jason-1卫星测高仪)和模拟波高对比。




关键参数设置:
           


# Segment 1

                Init {}

                {

                D = 0

              }

                InitOkada D

                {

                x = 94.57

                y = 3.83

                depth = 11.4857e3

                strike = 323

                dip = 12

                rake = 90

                length = 220e3

                width = 130e3

                U = 18

                }

                # Initial water level is at z = D

                Init

              {

                start = 0

              }

                {

                P = MAX (0., D - Zb)

                }


            # Segment 2

                EventList

                {

                start = 212

                step = 6

                end = 272

                }

                {

                Init {}

                {

                  D = 0

                }

                InitOkada D

                  {

                    x = 93.90

                    y = 5.22

                    depth = 11.4857e3

                    strike = 348

                    dip = 12

                    rake = 90

                    length = 150e3

                    width = 130e3

                    U = 23

                  }

                }

                # make sure the deformation is well resolved

                AdaptGradient

                {

                start = 212

                istep = 1

                end = 272

                }

              {

                cmax = 0.05

                cfactor = 2

                minlevel = 5

                maxlevel = LEVEL

                } D

                # add it to the current elevation field (only if wet)

                Init

                {

                start = 272

                }

                {

                P = (P < DRY ? P : MAX (0., P + D))

                }


            # Segment 3

                EventList

                {

                start = 528

                step = 6

                end = 588

                }

                {

                Init {}

                {

                  D = 0

                }

                InitOkada D

                {

                  x = 93.21

                  y = 7.41

                  depth = 12.525e3

                  strike = 338

                  dip = 12

                  rake = 90

                  length = 390e3

                  width = 120e3

                  U = 12

                }

                }

                # make sure the deformation is well resolved

                AdaptGradient

                {

                start = 528

                istep = 1

                end = 588

                }

              {

                cmax = 0.05

                cfactor = 2

                minlevel = 5

                maxlevel = LEVEL

                } D

                # add it to the current elevation field (only if wet)

                Init

                {

                  start = 588

                }

                {

                P = (P < DRY ? P : MAX (0., P + D))

                }


            # Segment 4

                EventList

                {

                start = 853

                step = 6

                end = 913

                }

                {

                Init {}

                    {

                  D = 0

                    }

                InitOkada D

                    {

                  x = 92.60

                  y = 9.70

                  depth = 15.12419e3

                  strike = 356

                  dip = 12

                  rake = 90

                  length = 150e3

                  width = 95e3

                  U = 12

                }

                }

                # make sure the deformation is well resolved

                AdaptGradient

                {

                start = 853

                istep = 1

                end = 913

                }

                {

                cmax = 0.05

                cfactor = 2

                minlevel = 5

                maxlevel = LEVEL

                } D  

                # add it to the current elevation field (only if wet)

                Init

                {

                start = 913

                }

                {

                P = (P < DRY ? P : MAX (0., P + D))

                }


            # Segment 5

                EventList

                {

                start = 1213

                step = 6

                end = 1273

                }

                {

                Init {}

                {

                  D = 0

                }

                InitOkada D

                {

                  x = 92.87

                  y = 11.70

                  depth = 15.12419e3

                  strike = 10

                  dip = 12

                  rake = 90

                  length = 350e3

                  width = 95e3

                  U = 12

                }

                }


免责声明:本文系网络转载或改编,未找到原创作者,版权归原作者所有。如涉及版权,请联系删

QR Code
微信扫一扫,欢迎咨询~

联系我们
武汉格发信息技术有限公司
湖北省武汉市经开区科技园西路6号103孵化器
电话:155-2731-8020 座机:027-59821821
邮件:tanzw@gofarlic.com
Copyright © 2023 Gofarsoft Co.,Ltd. 保留所有权利
遇到许可问题?该如何解决!?
评估许可证实际采购量? 
不清楚软件许可证使用数据? 
收到软件厂商律师函!?  
想要少购买点许可证,节省费用? 
收到软件厂商侵权通告!?  
有正版license,但许可证不够用,需要新购? 
联系方式 155-2731-8020
预留信息,一起解决您的问题
* 姓名:
* 手机:

* 公司名称:

姓名不为空

手机不正确

公司不为空