基于二次插值的非饱和土壤水流问题的特征差.pdfVIP

维普资讯 应 用 数 学 MATHEMATICA APPLICATA 2006。19(i)：i59～168 基于二次插值的非饱和土壤水流 问题 的 特征差分方法及数值模拟 宋丽叶 (中国科学院大气物理研究所，北京 100029) 摘要 ：本文针对一类非饱和土壤水流 问题，提 出了基 于二 次插值 的特征一差分格 式，得 到了严谨的L。模误差估计．并作 了数值试验，指明方法的有效性． 关键词 ：非饱和土壤水流；特征差分法 ；二次插值 ；误差估计 ；数值模拟 中图分类号：0241．82：O241．3 AMS(2000)主题分类 ：65MO6；65M25 文献标识码 ：A 文章编号 ：1001—9847(2006)Ol一0159—10 1．引言 均质土壤中的地下水流动可归结为非饱和土壤水的流动．由于非饱和流动的数学模型归 结为非线性 的偏微分方程 ，除了一些很特殊 的情况外 ，很难得到解析解．罗振东[2采用混合有 限元法建立的非饱和土壤水流的守恒形式 ，可以同时求出地下水及其通量的分布，用以统一计 算剖面入渗、蒸发、蒸腾和再分配以及这些现象交替出现时的水流运动过程．由Douglas，Rus— sellc 等提出的特征差分方法 ，可以大幅度加大时间步长 ，提高计算精度．利用这一方法 ，本文 针对均质土壤 ，地下水埋藏很深 的第二类非线性边值条件下的非饱和土壤水流问题建立了数 值模型，得到了基于二次插值的L 模误差估计． 设 G一 (O，L)，H (G)表示在 G内直到一阶导数平方可积的Sobolev空间．令 H (G)一 { ∈H (G)；v(L)一0)，(·，·)表示 G上的L 内积．定义 ’ I1 ‘ ，一(o，T； ·(G))一I／I{ ：[o，刀一W，(G)，I『I『~(O,TIWm一(㈤+cx。，s一0，1，…，kJ1． II l_．II 其 中 lI lI c0ITt c一 ess 。 0 -简记X(O，T；z(G))兰x(z)·用e表示小 的正常数 ，M～一M～(s， ，…，5r)表示与量 s，S，…，，有关的正常数．在不同估计式 中，同一符 号取不 同数值． 2．基于二次插值的特征一差分格式 我们考虑一维非饱和流问题 ，含水率有不同的时空分布．设 2轴垂直向下 ，坐标原点取为 地面 ，Q(z，￡)为在t时刻离地面距离为z处 的土壤含水率．假设地面有随时间变化的入渗或 · 收稿 日期 ：2005—04—04 基金项 目：国家重点基础研究专项资助项 目(2005CB321703)，国家 自然科学基金项 目 作者简介：宋丽叶，女 ，汉 。山东人 ，博士，主要从事冻土水热耦合方面的研究． 维普资讯 16O 应 用 数 学 2006 蒸发翠，人渗为正，蒸发为负，在底邵含水翠确随时l司而变化的分 ·则根据Darcy足俸，放j王缤 性原理 ，非饱和土壤水流 问题可归结为下面的模型方程 ： 求使得对于任意的 T≥ 0满足 ： f一 (D(Q) )+ 一s 0z ∈(o，丁)； JQ(z’O)=Q0(， 。≤ (2．1) IQ(L，)一卢()， ￡∈ (O，T)； IK(Q)一D(Q)!一q(￡)， z—o，￡∈(o，T)． 其中各参量的物理意义见 [2]．为了将边界条件齐次化 ，令 (z，￡)一Q(z，)一fl(t)．