非饱与土壤水流问题的混合有限元法及数值模拟.pdfVIP

非饱和土壤水流I可题的混合有限元法及其数值模拟4 罗振东 锄 (中国科学院计算数学所，科学与工程计算国家重点实验室，北京100080) 谢正辉朱江 (中科院大气所大气科学和地球流体力学数值模拟国家重点实验室．北京100080) 摘 要 本文利用混合有限元法研究非饱和土壤水流问题，建立了同时求解非饱 和土壤水流动含水量和通量的数值计算模型，讨论7其广义解、半离散和全 离散的混合有限元解的存在性及离散解的误差．该数值计算模型可以用来统 一计算剞面入渗，蒸发，蒸腾和再分配以及当这些现象交替出现时的土壤畲 水量和水分通量分布．利用混合有限元法，不仅能同时求出含水量和通量， 而且能提高通量精度及使计算稳定．该数值计算模型可用于高精度高分辨率 陆面模式． 1 非饱和土壤水流动问题及其广义解的存在性 非饱和土壤水的流动在大气科学、土壤学、农业工程、环境工程和地下水动力学等方面 具有重要意义，是大气科学的热门课题【1,2】．该问题可归结为非饱和流动问题[3—7】．由于 方程是非线性的，除了～些很特殊的情况外，不能得到解析解．因此，典型的作法是用数值 逼近求解非饱和流方程．在陆面过程研究中，不仅要模拟好土壤含水量的分布，而且还要使 其通量有较好的近似．传统的土壤含水量方程的离散方法是通过先求含水量的分布，由此 作差商得到通量分布．本文利用混合有限元方法建立的非饱和土壤水流动的数值模型，可同 时求出含水量和通量，使通量的精度提高，并使计算稳定，它可用于高精度高分辨率陆面模 式，并统一计算剖面入渗，蒸发，蒸腾和再分配以及当这些现象交替出现时的含水量和通量 分布．基于大气环流模式的水平分辨率(1-5个经纬度)，土壤水在水平方向上的流动可以忽 略．则根据Darcy定律及连续性原理，非饱和土壤水流问题可归结为下面的一维模型方程： ㈦-，甏一囊㈣塞)十警幽，№¨即，T， (1．2) e(z，0)=eo(0，O≤zSL； (13) e(L，t)=口(￡)，t∈(0，r)； (1 ●、 4)Ⅳ(目)一D(日)宴Oz=q(f)，。=o，t∈(o，即 其中口表示体积含水率，一sf是根系吸水率，K(e)是水力传导系数，D(o)是土壤水扩散 ． 率(参见f3，41)，z是正向向下的垂直坐标．土壤水力传导比价格系数Ⅳ和土壤水扩散率口 国后科学基金的合资助资助． 426 与 日 的 如 下 3 K 、，●／一 讣 圣 m D 也虬上 J— 期{系，．。。J、．。。。L ㈣㈣∽ ∞㈤垫以o一以曲一一虬 黼。，任，一‘ 其中以和冉分别表示土壤水分饱和含水率和残余含水率，而且饱和水传导率K，，土质参数 下边界含水量和初始时刻的含水量．条件(2)，(3)，(4)分别是初始条件、下边界条件和上边界 条件，见[3,41． 注：上边界条件还其它的情况． ．． 界条件齐次化，令 r F(z，t)=o(z，t)一口(￡)， 1 F(z，t)=v(z，t)一q(t)