土壤水流与溶质耦合运移问题的混合元迎风广义差分法研究.pdfVIP

2013年2月 计算数学 第35卷第1期 Feb．，2013 MATHEMATICANUMERICASINICA V01．35．No．1 土壤水流与溶质耦合运移问题的混合元 一迎风广义差分法研究q) 李焕荣 (重庆工商大学数学与统计学院，重庆400067； LASG，中国科学院大气物理研究所，北京100029) 摘 要 本文研究了一维非饱和土壤水流与溶质耦合运移问题的数学模型，建立了求其数值解的守恒 混合元一迎风广义差分格式．对非线性土壤水分入渗方程，采用守恒混合元法进行离散模拟，同时 得到了土壤含水量和水分通量；而对对流一扩散形式的溶质运移方程，利用迎风的广义差分法离 散求解．且分析了解的存在唯一性，并讨论了误差估计．最后给出数值算例，模拟结果表明利用本 文格式来求解非饱和土壤水流与溶质耦合运移问题是可靠的，且该格式具有稳定性和可实用性． 关键词：守恒混合元法；迎风广义差分法；数值模拟；水流与溶质耦合运移问题 MR(2000)主题分类：65N30，65M30 1．引 言 非饱和土壤水中溶质运移模型的研究对于精确地计算和预测土壤中养分和污染物运移规 律，改良盐渍土和防止土壤次生盐渍化等农业生产活动，水资源开发利用以及控制地下水污染 等都有着重大的现实意义． 溶质在非饱和土壤水中的运移问题可以简化描述为由一个非线性土壤水流方程和一个以 对流为主的对流一扩散方程组成的数学模型[1--4】．由于该模型的复杂性，一般难以求其解析 解【5】故常采用数值方法来求解．近几年，不少国内外学者已对非饱和土壤水中溶质运移问 元方法阳】进行了数值求解．但有限差分法对边界条件及土壤参数极为敏感，误差相对较大， 因此梅一等人【9]通过引入数值弥散因子来控制数值弥散减小了误差；而间断Galerkin方法 和有限元方法的模拟计算量又相对较大，因此本文作者发展了基于POD的降维有限元方法 【10】，在确保精度的前提下减小了计算量．然而以上研究工作只把水分流动考虑为一种稳定状 态，均没有考虑溶质随水流运移的动态耦合问题，因此本文作者用广义差分法111J模拟了该耦 合模型的数值解，同时求得了水分含量和溶质浓度随时空的分布；虽然该工作取得了较好的模 拟结果，但在该研究工作中，溶质方程中出现的水分通量只是通过含水量的差商得到的，并没 有考虑水分通量的直接求解． 因此本文将采用守恒混合元法和迎风广义差分法[12,a3]对非饱和土壤水流与溶质耦合运 移问题进行数值模拟，一方面对非线性土壤水分入渗方程，采用混合元法进行模拟，可以同时 +2011年12月5日收到． 课题资助． 万方数据 2 计算数学 得到土壤含水量和水分通量，避免了传统的方法通过差商得到水分通量时造成的数值震荡；另 一方面，利用迎风的广义差分法离散模拟对流．扩散形式的溶质运移方程，不仅具有较高的精 度且相对有效地克服了数值弥散现象． 2．水流与溶质耦合运移的守恒混合元一迎风广义差分格式 描述非饱和土壤一维溶质运移的基本方程[1j14】为： 掣=未(础，塞)一掣慨 江·， 方程(2．1)右端第一项描述溶质运移中的弥散作用，第二项描述溶质运移中的对流作用．式中， 是源汇项，c(x，t)为土壤溶液浓度(mL／L)，其溶质可为盐分、化肥、农药中的有害物质、 重金属等；O(x，t)为土壤容积含水量(cm3 为水分通量(cm／min)；t为时间坐标(min)． 土壤含水率除满足溶质运移基本方程外，还应满足非饱和土壤水分运动基本方程11516】： 塞=￡(即)赛)+百OK(O)懈 (2．2) 义同前．若z为水平坐标，则(2．2)式右端无第二项． 试验资料11 K(p)，百