土壤水动力学SWD土壤水运动的数值解法和软件应用教案.pptVIP

SUTRA由USGS开发，采用有限元法进行二维/三维、饱和/非饱和、变密度情况下地下水流动和溶质/能量运移的模拟]。我国已有报道介绍该软件，但应用尚不多见。国外多采用该模型进行海水入侵方面的模拟研究。 土壤水模拟软件HYDRUS HYDRUS1D、HYDRUS2D/3D系列软件 模型采用水头形式的Richards方程，可以对饱和-非饱和区进行统一模拟，可以处理诸如有压入渗、地下水位线波动情况下的区域水流运动问题。 在给定充分供水条件下的地面蒸发和作物蒸腾量条件下，该模型可以模拟土壤水盐胁迫下的实际根系吸水（作物蒸腾）和地面实际蒸发过程。 可以处理多种灵活的边界条件（大气边界、自由排水、定水位、动水位等），具有良好的用户界面和灵活的输入输出功能。模型采用有限元方法进行数值离散求解。 9.3 HYDRUS应用举例 （1） 应用软件模拟前要考虑： 计算区域和计算时段？ 土壤水分运动参数？ 如何给定初始条件？ 设起初为零通量稳定剖面，下边界含水率为0.3，各处初始含水率？ 边界条件 设左右边界为零通量边界 上边界--？ 下边界--？ 9.3 HYDRUS应用举例(2) 9.3应用举例 设定计算区域和计算时段 给定土壤水分运动参数 给定初始条件 设为地面零通量的稳定剖面，地下水埋深10m，各处初始基质势？ 边界条件 渠道边界？-定压力水头VS定土水势？ 其他边界-零通量边界 2 天后土壤含水率 2 天后土壤基质势 入渗强度变化 累积入渗量 * * 土壤水动力学 第9章土壤水动力学的数值求解与软件应用 土壤水动力学 (Soil Water Dynamics)第9章 土壤水动力学的数值求解与软件应用 (Numerical simulation application of Soil Water dynamics) 毛 晓 敏 中国农业大学水利与土木工程学院 Tel:Email: maoxiaomin@tsinghua.org.cn 第9章 土壤水动力学的数值求解与软件应用 土壤水分运动的数值计算方法 常用土壤水模拟软件 应用举例 9.1 土壤水分运动的数值计算方法 定解问题及常用数值计算方法 Richards方程的差分离散 边界条件处理 土壤水分运动参数取值 差分方程的求解 9.1.1 定解问题及常用数值计算方法 定解问题： Richards方程 边界条件 初始条件 求解方法： 解析法、半解析法： 数值方法： 解析法 ：求出解的函数表达式: 精度高，但有很大局限性:简化、特殊情况。 数值法：求用数值表示的有限个离散点和离散时段上的近似解: 具有广泛适用性；容易程序化，数值模拟软件多；随着计算机技术的提高和数值算法的改进得以迅速发展。 常用数值计算方法： 全离散化方法：时间和空间进行离散，利用微分方程推导出物理量在节点值的代数方程组，求解方程组得到微分方程在节点上的近似解 有限差分法(FDM)：用差商替代导数 控制容积法(CVM)：积分方程离散，守恒 有限元(FEM)：借助基函数，直接求函数的近似解 半离散方法：时间或空间离散 有限解析法（FAM） 线法(MOL) 有限差分法－利用Taylor 级数用差商代替导数： 前向差分： 后向差分： 中心差分： 一阶导数的近似表达式： 前差 后差 有限差分法的离散 ： i＝1 2 。。。。。。。。 I j＝1 2 。 。 。 。 。 。 。 J f i,j Δx 中心差 若采用中心差的形式，二阶导数如何表示？ 有限元法－用简单的插值函数代替每个单元上未知函数分布，形成一系列代数方程组的方法 有限元法的离散 ： 有限元法变化形式较多，如： 1）从权函数的选择来说，有配置法、矩量法、最小二乘法和伽辽金法； 2）从计算单元网格的形状来划分，有三角形网格、四边形网格和多边形网格； 3）从插值函数的精度来划分，又分为线性插值函数和高次插值函数等。 f j Galerkin法引入基函数ψi 并采用ψi作权函数，使得： 举例—Galerkin法： 对于： 可用矩阵表示为： 式中H、dH/dt分别为由h1，h2,…, hn和dh1/dt, dh2/dt, … dhn/dt组成的列矢量；D为导水矩阵，P为贮水矩阵。 若采用隐式方法： 通过求解大型稀疏矩阵方程，得到Hk+1 数值计算的收敛性与稳定性 收敛性：Δt、Δx→ 0时，差分方程解→微分方程的解 稳定性：计算误差有界 数值弥散问题： 对流占优情况下容易出现数值弥散，利用特