地下水溶质运移模型的优化数值算法.pdfVIP

第23鲞第4．期 武汉工业学院学报 v01．23No．4 2004年12月 Journalof Wuhan Dec．2004 PolytechnicUni”versity 文章编号：1009-4881(2004)04—0101—03 地下水溶质运移模型的优化数值算法 王惠芸1，马士进2 (1．武汉工业学院电气信息工程系，湖北武汉430023；2．河北工程学院土木工程系，湖北邯郸056038) 摘要：基于有限分析方法基础之上构造一种优化数值算法，阐述了该方法的基本理论和构 造过程，并对该数值方法进行稳定性分析。同时结合两个经典算例进行了计算和分析。结果 表明：该方法克服了有限分析系数中包含级数等计算复杂的缺点，能够消除数值弥散和振荡。 关键词：溶质运移；差分格式；对流一扩散；优化数值算法 中图分类号：023 文献标识码：A 0 引言 1 差分格式的构造 目前环境工程、石油工程以及水利工程中常常 考虑如下地下水溶质运移的定解问题 涉及到污染物在地下水中的迁移问题，而描述地下 誓+y等：仇誓 咒，￡∈[0，L]×[0，T] 水中溶质运动的控制方程为对流一扩散方程，可以 函 缸 ‘a矿 用此方程进行数值模拟和预测，但最终要归结为对 G(髫，O)=0 菇∈[0，￡] 该方程进行数值求解。数值求解的主要困难在于高 c(o，t)=Ot，C(￡，t)=卢t∈[0，T] Pelect数(即对流占优)的情形。目前常用的有限差 C(戈，0)=Co(戈) 髫E[0，L] 分格式有MacCormack格式¨1、特征型的Garlerkin (1) 格式和混合的E．一L旧1格式等，但它们或精度底或 其中c一污染物的浓度；卜流体空隙的平均流速； 格式复杂，对于以对流占优的问题常出现数值弥散和 仇一纵向弥散系数；L一计算时间段。 数值振荡。尽管许多学者对其进行了改进，但仍然克 设z．h分别为时间步长和空间步长，y、D。= 服不了以上存在的困难。八十年代初期由美国Iowa 大学的c．J．Chen提出的有限分析法L3J，其基本思想 是用控制方程的局部解析解组成整体的数值解，这样 格式，可得 的解可以较好地反映定解问题物理特性，并且可以得 T02cn+l (2) 到单调无震荡、数值稳定性好的数值解。为此，作者 垡+y型=DL 基于有限分析法方法基础之上，提出一种优化数值算 化简后可得如下微分方程： 法，该方法较有限分析法而言，一方面可以有效地克 服有限分析系数中级数计算繁杂的缺点，另外可以对 出2 P(3)7 可d2Cn+l一瓦VDL出idCn+1=警DLr 拟和因子进行修正，得到一种隐式的数值效应好的差 令譬：日，则对于常系数微分方程(3)式 分格式，并对该优化差分格式进行稳定性分析，同时 上，￡7 利用MATLAB软件H1数据可视化和程序易于实现的