新浪博客

[学习]GprMaxV2.0中GprMax3D输入文件的命令

2012-04-25 17:10阅读:
4.1 GPRMAX3D命令一般注释
大多数能够用在GPRMAX3D的命令几乎和GPRMAX2D是一样的。但是,有一些命令是GPRMAX3D独有的。为了简明扼要,这里只论述与GPRMAX2D不同的命令。
基本的空间和暂时离散步骤分别为t和△x,△y,△z

4.2 命令清淡
在版本2里面,GPRMAX3D一共有
42个命令用在3D GPR建模中。他们是:
#title
#domain:
#dx_dy_dz:
#time_step_stability_factor:
#messages:

#number_of_media:
#nips_number:
#media_file:
#geometry_file:
#medium:

#abc_type:
#abc_order:
#abc_stability_factors:
#abc_optimization_angles:
#abc__mixing_parameters:
#pml_layers:

#box:
#cylinder:
#cylinder_new:
#cylindrical_segment:
#sphere:
#plate:
#edge:
#triangle:
#bowtie:
#thin_wire:

#analysis:
#end_analysis:
#tx:
#rx:
#rx_box
#snapshot:
#tx_steps:
#rx_steps:

#huygens_surface:
#hertzian_dipole:
#voltage_source:
#transmissio_line:
#plane_wave:
#excitation_file:

4.3 GPRMAX3D命令参数
为了更好地介绍GprMax2D命令,我们将它们分成四类:
1. 一般命令:包括用于指定的大小和模型的离散
2. ABC相关命令: 允许定制和优化吸收边界条件
3. 介质和对象的构造命令: 用不同的参数来在模型中引入不同的介质和构造简单的几何形状
4. 激励和输出命令:用来放置源代码和模型的输出点

4.3.1 一般命令
#title
GprMax2D用法形同

#domain
用来指定模型的尺寸(单位:米)
#domain f1 f2 f3
f1f2f3分别是模型xyz方向上的大小

#dx_dy_dz:
指定xyz方向上的增量
#dx_dy_dz: f1 f2 f3
f1是空间步x方向上的增量,f2是空间步y方向上的增量,f3是空间步z方向上的增量。空间离散化控制最大允许时间步t和为了达到所需要仿真时间窗而提出的解决方案。△t和△x,△y,△z之间的关系是:
4.1
其中c为光速。在GPRMAX3D中,上式取等号。4.1中,一个小的x,△y,△z值导致△t得小值,这个△t小值意味着为了达到所给的仿真时间而需要更多的迭代次数。但是,需要指出更小的△x,△y,△z和△t值会让模型更精确。

#time_step_stability_factor
用法与GPR2D的相同。可以修改t的值。
#time_window
用法与GPR2D的相同。
#number_of_media:
用法与GPR2D的相同。
#geometry_file
用法与GPR2D相同。
#messages
用法与GPR2D相同。
#nips_number
用法与GPR2D相同。仅当GPRMAX3D被要求用时,才必须用到他。

4.3.2 ABC 相关命令
GPRMAX 3D这些影响Higdon ABCs配置和性能的命令与GPRMAX 2D用法相同。但是,GPRMAX3D用了更有力的PML ABC
#abc_type: pml
#abc_type: higdon 可以设置ABC的默认参数。

#pml_layers: i1
i1pml所占Yee单元的数目。默认值为8。所占用的Yee单元数越多,PML性能越好,但是所花的计算资源也越多。而且,有两点要注意:
a) PML是几何模型的一部分。但是,PML层里的域不参与计算并且用他们来计算就是错的。因此,不要把源和接收器设置在这个区域。PML的深度是用YEE 单元来度量而不是距离。
b) 注意:当前

我的更多文章

下载客户端阅读体验更佳

APP专享