无单元Galerkin方法中周期边界条件的处理.pdfVIP

第28卷 第1期 工 程 数 学 学 报 Vo1．28No．1 2011年02月 CHINESE JOURNAL OFENGINEERING MATHEMATICS Feb．2011 文章编~-：1005—3085(2011)01—0067-08 无单元Galerkin方法中周期边界条件的处理木 王晓东，欧阳洁，苏 进 (西北工业大学应用数学系，西安 710129) 摘 要：本文研究了无单元Galerkin方法中周期边界条件的处理技术 ，将Lagrange乘子法用于周期边界 条件的处理 ．数值计算结果表明，该方法具有较高的计算精度．另外，它与无单元Galerkin方 法中本质边界条件处理的Lagrange乘子法具有统一性，对于周期、本质混合型边界条件的处理 尤为方便． 关键词：周期边界；无单元Galerkin方法；无网格方法；Lagrange乘子法 分类号：AMS(2000)65M60；65N30 中图分类号：O241．8 文献标识码：A 1 引言 一 些微分方程的解在整个求解区域内会具有某种周期性，例如：电磁场问题、某些流动问 题、聚合物分子的取向分布问题等．针对这类问题，我们可以利用周期性将问题限定在一个周 期内进行求解，从而降低 问题规模并大大简化计算．如何在限定区域上给定边界条件，将直接 关系到微分方程的求解 ．通常可 以给定两种边界条件，一种是 自由边界条件，另一种是周期边 界条件．当解同时具有某种对称性，并且可以确定对称线或对称面的位置时，可以将求解区域 取为两个相邻对称线或对称面之间的区域，并且 由于解在对称线或对称面上导数为零，所 以可 以给定 自然边界条件．而当解不具有对称性或无法确定对称线或对称面的位置时，则一般只能 给定周期边界条件． 无单元Galerkin(elementfreegalerkin，EFG)方法 【卜3l作为一种新型的微分方程求解方法， 近年来受到很多学者的关注．该方法不依赖于网格或单元，具有易于进行 自适应分析、容易构 造高阶可导的近似 函数、计算精度高等众多优点．由于存在这些优点，EFG方法非常适合求 解诸如复杂流体多尺度模拟4【—61中具有周期性边界条件的Fokker—Planck方程 7【-10]等 问题．它 的高精度可 以很好地避免FokkerPlanck方程求解过程 中由于误差积累导致的取向分布概率密 度函数失去守恒性的问题．要实现EFG方法在这类 问题中的应用 ，首先就必须要解决周期边 界条件的处理问题． EFG方法使用移动最小二乘法[，2】构造形 函数 ，而移动最小二乘法不具有Kronecker delta函数的性质，所以有限元等方法中周期边界条件的处理方法不能直接用于EFG方法．本 文借鉴在 EFG方法 中使用 Lagrange乘子法 (11,12】处理本质边界条件的思想，提出了周期边界 条件处理的Lagrange乘子法 ．该方法精度高，并且与本质边界条件处理的Lagrange乘子法具 有统一性 ，易于实现 ． 收稿日期：2009—04—09．作者简介：王晓东(1986年5月生)，男，硕士．研究方向：无网格方法及复杂流体数值模拟 基金项目：国家自然科学基金 10590353)；国家重点基础研究发展计划 (973)项目(2005CB321704)． 68 工 程 数 学 学 报 第28卷 2 EFG方法的函数逼近 EFG方法使用移动最小二乘法进行未知函数的逼近和检验函数 的构造 ．移动最小二乘 法[，。】(movingleastsquares，MLS)是局部近似与最小二乘相结合的产物，该方法直接基于空 间离散点来进行近似，不需要对区域进行网格划分． 在区域Q中，未知标量函数u()的移动最小二乘逼近可以表示为 ∞