二维土壤溶质输运方程基于POD方法的降阶CN有限体积元外推算法.pdfVIP

2014年 8月 计 算 数 学 第 36卷第 3期 Aug．，2014 M ATHEMATICA NUMERICA SINICA Vo1．36，No．3 二维土壤溶质输运方程基于 POD方法的 降阶 CN 有限体积元外推算法术1) 腾 飞 (凯里学院数学科学学院 贵州凯里556011) 罗振东) (华北电力大学数理学院，北京102206) 摘 要 利用Crank—Nicolson(CN)有限体积元方法和特征投影分解方法建立二维土壤溶质输运方程 的一种维数很低、精度足够高的降阶 CN有限体积元外推算法，并给出这种外推算法的降阶 CN 有限体积元解的误差估计和算法的实现．最后用数值例子说明数值结果与理论结果相吻合，并阐 明这种降阶 CN有限体积元外推算法的优越性． 关键词：特征投影分解；土壤溶质输运方程；降阶 CN有限体积元外推算法 MR (2000)主题分类：65N12，65M15，65N30 1．引 言 土壤溶质输运方程 [1-3J是描述土壤中的溶质运动变化规律的方程．土壤的溶质含量对土 壤和环境变化有直接影响．因此，对土壤 中的溶质含量变化做预测预报有着很重要的实际意 义。由于土壤溶质输运方程往往含有复杂的源项和水分流通率及水动力扩散系数，因此，要想 求出其解析解是不容易的，有效的方法是求其数值解 [1-3]．文献 …给出了土壤溶质输运方程 的有限差分方法，文献 [2】给出了一种有限元方法．正如 [3，4】所说，与有限差分方法和有限元 方法相比，有限体积元方法有如下优点：有限体积元方法能保持局部的质量或能量守恒，比有 限差分法精度高而且能适应边界复杂的计算域，又与有限元方法有同阶精度但要比有限元方 法便于计算 (事实上，有限体积元方法最终可通过体积元转化差分格式计算，但可以借助有限 元方法做解的存在性、稳定性和收敛性的理论分析)，因此有限体积元方法被公认是最有效的 数值计算方法之一．这样，利用有限体积元方法求解土壤溶质输运方程比用有限元法和有限 差分方法更有优势．虽然文献 [33]已经给出了土壤溶质输运方程的一种时间一阶精度的经典 有限体积元方法，但是对于实际工程应用，经典的有限体积元格式也包含有很多 自由度 (等于 网格节点数 目)．这样，由于计算过程中截断误差的积累，使得经典的有限体积元格式计算若 干步后，就出现浮点溢出无法继续计算下去．因此，在保证有足够精度数值解的情况下，如何 降低有限体积元格式的自由度、简化计算、节省计算量和存储要求，减少计算过程中截断误差 的积累，是实际计算中需要解决重要的问题． 2013年8月 27日收到． )基金项目：国家自然科学基金 和贵州省教育厅自然科学研究项目(黔教合 KY字 。)通讯作者：罗振东，E—mail：zhdluo@ncepu．edu．cn． 258 计 算 数 学 大量实例已经证明，特征投影分解 (ProperOrthogonalDecomposition，简记为 POD)方 法 l5_是一种能够极大地减少数值模型 自由度即降低模型维数的有效逼近方法．该方法已经广 泛成功地应用于包括信号分析和模式识别 [6j、统计学及地球物理流体动力学和气象学f】等 各种领域．POD方法实质上是在最小二乘意义下寻找能代表已知数据的一组正交基，即一种 求已知数据的最优逼近方法．早年，POD方法主要用于统计中的主分量分析或寻找动力系统 的主要特征 [5-9J．直到最近，才有一些基于POD方法的降阶有限体积元格式 [1o-13】被建立． 虽然 1【3]已经提出了二维土壤溶质输运方程基于 POD方法的降阶有限体积元格式，但只有 一 阶精度时间．这样，为了得到足够高的时间精度，时间步长必须取做够小，从而增加了计算 迭代步数，也会增加计算过程中截断