tough2-eos9模块练习

Eos9范例问题描述：假设模拟区域为一个二维XZ剖面的非均质各向同性介质，尺寸为10m长×10m厚。顶部为回填土，底部为风化土，左右边界为无通量边界，顶部边界为固定大气压力边界（=1.013E5Pa），底部边界固定压力为1.503E5Pa。目标：1.包气带的水流模型（使用eos9）2.建立二维网格使用meshmaker3.给每个网格赋上岩石属性4.给定边界条件5.添加源汇项（地表入渗）6.得到初始条件7.模拟入渗测试8.时间步长和输出参数9.考虑非等温流（使用eos3）1.网格生成Meshmaker：2.岩性参数分四个区：SAPRO（主岩性）、ATMOS（上边界）、FILLI（2m厚的回填土）、BOUND（下边界）还需要给ELEME中每个网格赋上岩性：1.AssignRock：给网格赋上岩性SAPRO、FILLI、BOUND2.AddBound：可以很方便的把某一行或某一列的网格变成边界，适合于与空气接触边界3.边界条件垂向边界：无通量边界，不需要做任何改变底边界：恒定压力边界（狄利克雷边界），将底边界单元格的体积给1E50，同时，将ROCKS里BOUND的SPHT从1000改为1.0E5（不输出底边界单元格的质量）上边界：恒定的入渗流量，通过GENER块给定4.初始条件PARAM4.0：默认初始条件Sa=0.9INDOM：默认边界压力条件INCON：无5.运行模型达到稳定状态运行6.入渗测试将稳定状态的SAVE作为入渗测试的初始条件INCON通过GENER给定入渗量