求解自由边界问题的自适应投影方法.PDFVIP

第３４卷第５期 重庆工商大学学报（自然科学版） ２０１７年１０月 Ｖｏｌ．３４　ＮＯ．５　　　　　　　　ＪＣｈｏｎｇｑｉｎｇＴｅｃｈｎｏｌＢｕｓｉｎｅｓｓＵｎｉｖ（ＮａｔＳｃｉＥｄ）　 　　　　　　　 Ｏｃｔ．２０１７ － ｄｏｉ：１０．１６０５５／ｊ．ｉｓｓｎ．１６７２０５８Ｘ．２０１７．０００５．００２ 求解自由边界问题的自适应投影方法 钟艳丽，严月月，钟思超 （重庆师范大学 数学科学学院，重庆４０１３３１） 　　摘　要：对一类自由边界问题，提出了基于线性互补问题的自适应投影算法．采用有限差分格式将自由 边界问题离散为一个线性互补问题，然后用自适应投影迭代算法求其数值解，该方法在迭代过程中自动调 整参数，达到加快收敛速度的目的，每一步迭代只需要求解一个线性方程组．给出了具体算法过程，并利用投 影性质得到了它们的收敛性分析．最后用数值算例对算法验证，与已有的算法比较，结果表明：参数对自适应 投影算法影响较小，该方法收敛速度更快． 关键词：自由边界问题；有限差分；线性互补；自适应投影法 中图分类号：Ｏ２４１．１８２　　　文献标志码：Ａ　　　文章编号：１６７２０５８Ｘ（２０１７）０５０００７０６ ＝ ＝ ＝ 　　自由边界问题在所考虑区域内部的微分方程 ｎ １时，令 ［ａ，ｂ］，则边界 为ａ，ｂ两端点，Ｌ Ω Γ 是非线性的，通常不存在精确解或很难求其精确 ２  ＝ ＝ ；当ｎ ２时， 为平面有界闭区域，其边界 Ω Γ ２ ｘ 解，因此需要用数值方法求解．目前已有的方法有有  ２ ２ ［１５］   限差分法、有限元法等 ．文献［３］用有限差分法 ＝ ＝ ＋ Ω，ｘ （ｘ，ｘ），Ｌ ．如果已知函数 φ∈ １ ２ ２ ２ ｘ ｘ   将求解的问题离散为线性互补问题，提出了投影算 １ ２ ２ ２ Ｌ（）， Ｌ（ ），问题（１）等价于变分不等式问 Ω ψ∈ Γ 法．本文将障碍问题离散成