一维Boussinesq方程有限差分格式的讨论.pdfVIP

一 维 Boussinesq方程有限差分格式的讨论 王全荣 ，唐仲华 ，李星 ，翟莉娟 (1．生物地质与环境地质教育部重点实验室，中国地质大学环境学院，武汉 430074； 2．中国地质大学数学与物理学院，武汉 430074；3．中国地质大学外 国语学院，武汉 430074) 摘要：目前，采用 “transmissivity”(导水系数)法将 Boussinesq方程线性化后再求解是解决问题的主 要手段 ，其 中求解方法包括调和平均值法，几何平均值法和算术平均值法，为 了对 比它们 的优劣 性，结合水 头线方程的性质 (凹凸性)和平均值的性质 ，以均质各向同性的潜水稳 定流 BoussinesⅡ 方程为例 ，通过解析解，建立不同差分格式的误差方程 ，推导差分格式的收敛性条件 ，最后 ，通过 数值实验检验不同差分格式的误差，结果表明：若网格步长 (单元格长度)不趋于0，采用调和平 均值建立的差分方程误差比较大，不能收敛到精确解；采用几何平均值时，误差较 小，但是采用算 术平均值时，差分方程是无条件收敛性的。 关键词：Boussinesq方程 ；误差方程 ；平均值 中图分 类号 ：P641．2 文献标识码 ：A Thediscussion ofthefinitedifferenceschemeof1．D boussinesqequation WangQuanrong，TangZhonghua，LiXing~，ZhaiLijuan (1．KeyLaboratoryofBiogeologyandEnvironmentalGeologyofMinistryofEducation，Schoolof EnvironmentalStudies，ChinaUniversityofGeosciences，Wuhan 430074，China； 2．SchoolofMathematicsandPhysics，ChinaUniversityofGeosciences，Wuhan 430074，China； 3．TheCollegeofForeignLanguage，ChinaUniversityofGeosciences，Wuhan 430074，China ) Abstract：ThemainmethodtosolveBoussinesqequation istolinearizeitby “transmissivity”firstly，which is calculatedbythreeknownaveragemethods，such asharmonicmean，geometricmeanandarithmeticmean．In orderto contrasttheerrorscaused by differentaveragemethods，combiningthe property ofthefunction of hydraulichead(concave／convex)andtheaveragemethods，theerrorformulaofvaryingfinitedifferenceschemes ofBoussinesqequation isestablishedtoexaminetheconvergentconditions，andcomparetheerrors throughthe numericalexperiments．The resultsshow thatcontrolequation solved by the finite difference method in all isotropieandhomogeneousaquifer，ifthe d spacingdoesnotapproach zero，thefinite difference schemeof Boussinesqequationwhosethicknessoftheaquiferiscalculatedbytheharmonicmeandoesnotconverge，while theconvergenceoftheo