土壤溶质运移问题的最小二乘混合有限元法.pdfVIP

第35卷 第2期 肇 庆 学 院 学 报 Vo1．35，No．2 2014年3月 JOURNALOFZHAOQINGUNIVERSITY Mar．2014 土壤溶质运移问题的最小二乘混合有限元法 刘知雨，王桂霞，王 娟 (内蒙古师范大学 数学科学学院，内蒙古 呼和浩特 010022) 摘 要 ：给出了非粘性土壤溶质运移问题的最小二乘混合有限元法，利用该方法将方程降阶对方程进行离散，并 构造最小二乘混合有限元格式．最小二乘混合有限元方法可以避免标准混合有限元格式中的ladyzhenskaya。babus— ka—brezzi(LBB)F~$,]条件．通过分析此逼近格式的收敛性，给出相应的误差估计，误差表明此种数值方法具有最优的 收敛阶． 关键词：最小二乘；土壤溶质运移问题；混合有限元法：误差估计 中图分类号：O241．2 文献标志码：A 文章编号：1009．8445(2014)02，001204 0 引言 组成土壤液相的并不是纯水而是含有溶质的水溶质．土壤水中含有的溶质，不仅存在着土壤盐碱化的 问题，而且还会发展成为更为深远和广泛的水土环境问题Ⅲ．土壤中的溶质运移问题是非常复杂的，难 以求 出解析解．近年来，最 ／bz．乘混合有限元法的研究受到计算数学家们的重视 1．与传统的有限元法相比，最小 二乘混合有限元法在选择混合元空问时不需要满足LBB条件，空间选择更灵活，而且形成的系数矩阵对称 正定，在误差估计中，能得到最优的误差估计． 本文中，笔者将最小二乘混合有限元法推广到如下的非粘性土壤溶质运移问题中： 一 (，警+) ， c(z，o)=Co(z)，0zs ；c(L，力=C2(c)，￡∈(0，， (1) 一 击(，警+)一 ． 其巾：c(z，￡)为溶质浓度，S是源汇项，D ，为其水动力弥散系数，O(z，￡)，q(z，￡)和 ，0=q(z，t)／o(z，t)分别是土 壤含水率、土壤水分运动通量和空隙平均流速；c。(，c．和c2(t)分别为初始时刻边界溶质通量、上边界溶质 通量和下边界溶质浓度．为了便于理沦研究，假定Co(X，=cr(x，Y，=0． 1 方程离散及格式建立 本文用到的Sobolev空间都是标准的，Lz(O)中的内积用 (·，)·表示，令 m,p，1sp 。。是定义在 上的 收稿 日期：2014．01—02 基金项 目：国家自然科学基金资助项 目(1136l039)；教育部自然科学基金资助项 目(212029)；内蒙古自然科学基金资助项 目 (2013MS01l6)；内蒙古教育厅 自然科学基金资助项 目(NJZZI3037) 作者简介：刘知雨(199O一)，女(蒙古族)，内蒙古通辽人，内蒙古师范大学数学科学学院硕士研究生 通信作者：王桂霞(1968。)，女，内蒙古赤峰人，内蒙古师范大学数学科学学院教授 ． 一 一 囹簦 ±塑 堕垂整 墅 垦 三 坌查 Z垦 ． ． !三 Sobolev~．Ihq，相应的范数为I1 ，特别当p=2时，把W 记作 ，其范数简记为 ，m 0．定义内积与空间 如下 ： (“，=L删ddy，VM，veL )； (，∞)=∑(，∞)，Vo，weL(， H：{oJeL )；divweL )；=o，z--o}， 5={t，∈’(，t，=0在r上)． 给定时间步长At=T／(』、r)(Ⅳ是正整数)，令时间层 =nAt，n=0，1，2，…，N．HhCH．SCS是相应的有限元 空间．对于CeS，定义其椭圆投影尺c∈．s，使得 (V(Rc—c)，7v^)=A(Rc—C,U^)，Vt，^∈Js^． (2) 其中A是正常数．则有以下的误差估计 IIc-RclI㈣ …IIc)Il。㈣，Ilc,(t)一c))II