HydroGeoSphere/输运条件
- twc 输运解的时间加权因子。
为输运解分配一个新的时间加权因子值,默认值为 0.5。值可以从 0.0(显式)到 0.5(中心或 Crank-Nicolson)或 1.0(完全隐式)。完全隐式时间加权不太可能表现出振荡,但比中心时间加权更容易发生数值弥散。
- • • •
- • • •
- pectol 佩克莱特数。
为佩克莱特数分配一个新值,默认值为 1 × 1020。佩克莱特数不会以任何方式影响解,它仅仅用于生成警告信息。佩克莱特数可以计算为
其中 是流速, 是单元长度,在本例中是 x 方向,而 是弥散系数。它是网格细化程度的衡量指标,较大的数值表明空间离散化较差。
较大的默认值抑制了警告信息的生成。
- • • •
- • • •
此指令使 HydroGeoSphere 将每个单元的佩克莱特数写入文件 prefixo.peclet_number,仅在第一个时间步长。
- • • •
- • • •
- courtol 库朗数。
为库朗数分配一个新值,默认值为 1 × 1020。库朗数不会以任何方式影响解,它仅仅用于生成警告信息。库朗数可以计算为
其中 是流速, 是单元长度,在本例中是 x 方向,而 是时间步长。它是时间步长大小的衡量指标,较大的数值表明时间离散化较差。
- • • •
- • • •
- restolc 输运求解器收敛准则。
为输运求解器收敛准则分配一个新值,默认值为 1 × 10−10。当此准则小于最新输运解的残差(误差)的绝对最大值时,就会获得收敛。
- • • •
- • • •
- isolv_infoc 输运求解器详细程度。
为输运求解器详细程度分配一个新值,默认值为 0。此值控制打印到屏幕和列表文件中的求解器性能信息的详细程度,可以取 0(无信息)或 1(完整信息)。
- • • •
- • • •
- maxtit 输运求解器最大迭代次数。
为输运求解器最大迭代次数分配一个新值,默认值为 2000。
- • • •
- • • •
- almax,btmax,gammax 分别为x-、y- 和 z-方向的速度上游加权因子。
导致使用速度的上游加权,而不是默认的中央加权。数值范围从 0.0(无上游加权)到 1.0(完全上游加权)。请注意,这些变量不适用于控制体积情况,在这种情况下,当 upstrvel 为 .TRUE. 时,始终使用完全上游加权。
- • • •
- • • •
此指令使HydroGeoSphere 使用 van Leer 通量限制器,如第 3.9.1 节所述。
- • • •
- • • •
- maxiter_flim, resmax_flim, delmax_flim 分别为非线性通量限制器迭代过程的最大迭代次数、残差和绝对收敛准则。
为影响 van Leer 通量限制器行为的参数分配新值,这些参数默认分别为 15、1 × 10−5 和 1 × 10−8。
- • • •
- • • •
此指令仅影响耦合地表/地下系统中的传输,并使HydroGeoSphere 忽略两个区域之间的弥散/扩散交换。这样做是为了评估对流交换相对于弥散交换的重要性,或在将HydroGeoSphere 与忽略弥散交换的代码进行比较时使用。
- • • •
- • • •
此指令仅影响耦合裂隙/地下系统中的传输,并使HydroGeoSphere 忽略两个区域之间的弥散/扩散交换。这样做是为了评估对流交换相对于弥散交换的重要性,或在将HydroGeoSphere 与忽略弥散交换的代码进行比较时使用。
- • • •
- • • •
以下三个指令可用于定义用于标记输出和选择性停止运行的阈值浓度。
- detection_threshold_conc 检测阈值浓度 [M L−3]。此指令设置将用于控制以下两个指令行为的值。
- • • •
- • • •
导致HydroGeoSphere 使用字符串 > 检测阈值标记观测井输出,如果节点处的浓度超过上面定义的detection_threshold_conc。
- • • •
- • • •
如果任何观测井节点处的浓度超过上面定义的detection_threshold_conc,则导致HydroGeoSphere 停止执行。
- • • •
- • • •