matlab由频域变时域的代码使用PETSc在均匀各向同性介质中的3D声波传播
PETSc-
在这里,您可以找到C
+
PETSc实施方案解决3D声波方程的示例。
使用Krylov方法来找到Ax
=
b的迭代近似解
分散细节:
时域有限差分(FDTD)
隐式时间步进
O(2,4)
从泰勒级数推导的方案:
在空间[-1:16:-30:16:-1]
/
12dx2
及时[2:-5:4:-1]
/
dt2
型号详情
各向同性
同质
Dirichlet边界条件
使用方法:
必须安装PETSc。
确保设置了PETSC_DIR和PETSC_LIB环境库
make
all
./run_O22.sh
或者./run_O24.sh
可以在Shell脚本中更改运行时选项和处理器数量。
在代码中或从运行时更改标记,可以保存并绘制整个波场或仅在接收器位置绘制地震图。
运行时选项-vel
float-传播速度[km
/
s]
-xmax
float-模型尺寸[km]
-ymax浮点-zmax浮动-dt
float-时间步长[秒]
-tmax
float-最大模拟时间-isrc
int-源位置[网格点]
-jsr
2022-03-06 13:20:16
4.25MB
系统开源
1