第5章自适应时间步进与输出控制（ATO4）包

5.1 概述

在 MODFLOW 的代码体系中，针对时域问题（用于瞬态流模拟），采用后向差分公式进行离散化处理。时间步长的设定依赖于若干参数的预先指定，包括时间步数 (NSTP)、应力周期长度 (PERLEN)，以及用于按几何级数扩展时间步长的乘法因子 (TSMULT)。以下列举了该方法所存在的若干限制条件。

• 若矩阵求解器或非线性迭代方法（适用于无约束情形）未能在当前时间步内实现收敛，则整个计算过程将被迫终止。
• 即便最终得到了收敛的数值解，预设的时间步长可能并不具备较高的计算效率。
• 时间步长必须遵循几何级数增长规律，这导致用户在设定特定时间点以观察流动系统状态时缺乏足够的灵活性。因此，用户只能请求输出与时间序列相匹配的时间点，而这些点需在模拟开始前由用户明确指定。

为应对上述问题并提升整体解决方案的运行效率，开发了自适应时间步进与输出控制模块（ATO4）。该自适应方案依据系统预期非线性的程度动态选择合适的时间步长。当预计非线性程度较小时，将采用较大的时间步长以加快模拟进程；而当预计非线性较强时，则选择较小的时间步长以保障收敛性。若在当前设定的时间步内无法实现收敛，则进一步缩小时间步，并重复进行求解操作。

对于采用自动时间步进机制的情况，仿真结果对应的具体时间节点是未知的；因此，ATO4 模块还新增了专门用于输出控制的功能单元，该单元负责确定所