二维土壤溶质输运方程的有限体积元格式.pdfVIP

2012年2月 计算数学 第34卷第l期 Fbb．2012 MATHEMATICANUMERICAS】NICA V01．34．N01 二维土壤溶质输运方程的有限体积元格式+1) 刘 群 孙 萍 罗振东 (贵州师范大学数学与计算机科学学院，贵阳550001) 摘 要 本文导出二维的土壤溶质输运方程的有限体积元格式，并分析其误差通过数值例子说明，有 限体积元格式比有限元格式稳定 关键词：土壤溶质输运方程；有限体积元格式；有限元格式；误差估计 MR(2000)主题分类：65N30，35Q10，65M30 1．引 言 土壤所含的溶质不仅使得土壤发生盐碱化，而且会直接影响到环境问题．因此，对溶质在 土壤中的输运规律的研究，具有重要的实际意义土壤中的溶质输运问题往往归结为带有复 杂的源项及与土壤水分相关的系数的偏微分方程(参看[1—2D，因此，要求出其解析解是不容 易的，有效的方法是求其数值解(参看[1—3])[1】中给出了土壤中的溶质输运问题的有限差分 方法，[2，3】中给出了一些的有限元方法．由于有限体积元法(参见阻6】)，也称为盒子法(参见 【7])或广义差分方法(参见【8，9])，能保持局部的质量或能量守恒，比有限差分法精度高而且能 适应边界复杂的计算域，又与有限元方法有同阶精度但要比有限元方法便于计算，因此有限体 积元法被公认是最有效的数值计算方法之一，已经被广泛应用于求解各种类型的偏微分方程 (例如，二阶椭圆方程、抛物型方程、stokes方程等)的数值解(参见p14】) l7]提出了一些土壤溶质输运方程的有限体积元格式，但都没有对 据我们所知，虽然[15 格式做误差分析因此，本文导出二维的土壤溶质输运方程的有限体积元格式，并分析其误差． 通过数值例子说明，有限体积元格式比有限元格式稳定 本文的安排如下，第2节导出二维的土壤溶质输运方程的有限体积元格式；第3节给出 二维的土壤溶质输运方程的有限体积元解的误差分析；第4节给出一些数值例子，并通过与 有限元格式的计算结果比较说明有限体积法的优越性；最后，第5节是结论和展望 2．二维的土壤溶质输运方程的有限体积元格式 根据土壤溶质输运的对流和水动力扩散(即分子传播和水力学扩散)原理，二维的非饱和 土壤中溶质输运方程可以归结为下面的方程(参见[1])． ’2011年2月28日收到 1’基金项目： A2010001663)和贵州省科技计划项目(批准号黔科合J字∞11]2367) 万方数据 计算数学 2012年 f Q一(口1(”，日)％)。 {c(z，可，o)=co(。，可)，(z，掣)∈n， (21) 【c(z，Ⅳ，t)=cr(z，Ⅳ，t)，(z，Ⅳ)∈r，￡∈(o，T)， 其中n 见，不失一般性，我们不妨在下面的理论分析中假定cr(z，9，￡)=co(z，∥)=o 成 问题II．求c(￡)：[o，卅1÷u满足 ‰叭 』(c“叫佃(c，叫“(c，”)-㈣1叫， (22) I c(z，掣，o)=o，(z，掣)∈n， 积 设r是时间步长而且扩是c(t)在如=nr(n=o，1，，Ⅳ=[T／T】)处关于时间离散的 逼近，那么问题II关于时间的向后一步半离散化格式可写为： 问题III．求c”∈u(1≤n≤Ⅳ)满足