KRAKEN计算群速度(弹性海底)

KRAKEN输出的.prt文件中包含的群速度是根据液态海底的一个公式计算的,因此考虑弹性海底时群速度结果有误,此时可根据定义d\omega/dk_n计算群速度。特征值k_n可从.mod文件中读取,\omega=2\pi f,简单的数值差分就可以得到正确的群速度了。

下图为一个弹性海底的仿真示例,左图结果为.prt文件中的群速度,右图为真实群速度。可见无剪切波时左右两图的群速度曲线一致,有剪切波时左右两图的群速度曲线存在显著差异。

左图结果为.prt文件中的群速度,右图为真实群速度

KRAKEN作者Michael B. Porter的解释:

Are you taking the group speed out of the print file? If so, then the answers will be wrong because the group speed is calculated using a formula from our book that is only valid for purely acoustic problems. The alternative way you can do this is to read the wavenumbers out of the mode file and calculate the group speed as a numerical derivative dk/ d omega or something like that.

©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

相关阅读更多精彩内容

友情链接更多精彩内容