基于重心插值的二维抛物型方程的局部DQ法.pdfVIP

３３４ 　西安理工大学学报 ＪｏｕｒｎａｌｏｆＸｉ’ａｎＵｎｉｖｅｒｓｉｔｙｏｆＴｅｃｈｎｏｌｏｇｙ（２０１３）Ｖｏｌ．２９Ｎｏ．３　 　　文章编号：１００６４７１０（２０１３）０３０３３４０４ 基于重心插值的二维抛物型方程的局部 ＤＱ法 路小平，赵凤群，蒋卓韵 （西安理工大学 理学院，陕西 西安７１００５４） 摘要：针对传统微分求积法的局限性，提出一种基于重心插值的局部微分求积法，并应用于二维微 分方程的求解。在该方法中选择重心插值函数作为基函数以保证方法具有很好的数值稳定性。此 外，局部微分求积法能够克服微分求积法中节点过多出现的弊病。因而，本研究方法除了具有传统 微分求积法计算量少、精度高等优点外，还具有数值稳定性好、节点可以取到很多的优点。以二维 Ｂｕｒｇｅｒｓ方程组为例，数值结果表明了该算法的有效性。 关键词：重心插值；局部微分求积法；Ｂｕｒｇｅｒｓ方程组 中图分类号：Ｏ２４１　　　文献标志码：Ａ ＡＬｏｃａｌＤｉｆｆｅｒｅｎｔｉａｌＱｕａｄｒａｔｕｒｅＭｅｔｈｏｄＢａｓｅｄｏｎＧｒａｖｉｔｙＩｎｔｅｒｐｏｌａｔｉｏｎｆｏｒ ＳｏｌｖｉｎｇＴｗｏＤｉｍｅｎｓｉｏｎａｌＰａｒａｂｏｌｉｃＤｉｆｆｅｒｅｎｔｉａｌＥｑｕａｔｉｏｎｓ ＬＵＸｉａｏｐｉｎｇ，ＺＨＡＯＦｅｎｇｑｕｎ，ＪＩＡＮＧＺｈｕｏｙｕｎ （ＦａｃｕｌｔｙｏｆＳｃｉｅｎｃｅ，Ｘｉ’ａｎＵｎｉｖｅｒｓｉｔｙｏｆＴｅｃｈｎｏｌｏｇｙ，Ｘｉ’ａｎ７１００５４） Ａｂｓｔｒａｃｔ：Ａｉｍｉｎｇａｔｔｈｅｄｉｓａｄｖａｎｔａｇｅｓｏｆｔｒａｄｉｔｉｏｎａｌｄｉｆｆｅｒｅｎｔｉａｌｑｕａｄｒａｔｕｒｅｍｅｔｈｏｄ，ａｌｏｃａｌｄｉｆｆｅｒｅｎｔｉａｌ ｑｕａｄｒａｔｕｒｅｍｅｔｈｏｄｂａｓｅｄｏｎｇｒａｖｉｔｙｉｎｔｅｒｐｏｌａｔｉｏｎｉｓｐｒｏｐｏｓｅｄｉｎｔｈｉｓｐａｐｅｒ，ａｎｄｔｈｅｔｗｏｄｉｍｅｎｓｉｏｎａｌｐａｒａ ｂｏｌｉｃｅｑｕａｔｉｏｎｓａｒｅｓｏｌｖｅｄｂｙｔｈｅｎｅｗｍｅｔｈｏｄ．Ｉｎｔｈｉｓａｐｐｒｏａｃｈ，ｔｈｅｇｒａｖｉｔｙｉｎｔｅｒｐｏｌａｔｉｏｎｆｕｎｃｔｉｏｎｉｓｃｈｏ ｓｅｎａｓｔｈｅｂａｓｉｓｆｕｎｃｔｉｏｎ，ｗｈｉｃｈｇｕａｒａｎｔｅｅｓｔｈａｔｔｈｅｍｅｔｈｏｄｈａｓａｖｅｒｙｇｏｏｄｎｕｍｅｒｉｃａｌｓｔａｂｉｌｉｔｙ．Ｉｎａｄｄｉ ｔｉｏｎ，ｔｈｅｌｏｃａｌｄｉｆｆｅｒｅｎｔｉａｌｑｕａｄｒａｔｕｒｅｍｅｔｈｏｄｃａｎｏｖｅｒｃｏｍｅｔｈｅｄｒａｗｂａｃｋｏｆｄｉｆｆｅｒｅｎｔｉａｌｑｕａｄｒａｔｕｒｅｍｅｔｈ ｏｄｗｈｅｎｉｔｈａｓｔｏｏｍａｎｙｎｏｄｅｓ．Ａｃｃｏｒｄｉｎｇｌｙ，ｔｈｅｐｒｏｐｏｓｅｄｍｅｔｈｏｄｎｏｔｏｎｌｙｉｎｈｅｒｉｔｓｔｈｅａｄｖａｎｔａｇｅｓｏｆｔｈｅ ｔｒａｄｉｔｉｏｎａｌｄｉｆｆｅｒｅｎｔｉａｌｑｕａｄｒａｔｕｒｅｍｅｔｈｏｄｓｕｃｈａｓｌｅｓｓｃａｌｃｕｌａｔｉｏｎａｎｄｈｉｇｈｐｒｅｃｉｓｉｏｎ，ｂｕｔａｌｓｏｏｗｎｓｔｈｅ ｍｅｒｉｔｓｓｕｃｈａｓｇｏｏｄｎｕｍｅｒｉｃａｌｓｔａｂｉｌｉｔｙａｎｄｎｏｄｅｓｅｌｅｃｔｅｄｆｌｅｘｉｂｉｌｉｔｙ．ＴｗｏｄｉｍｅｎｓｉｏｎａｌＢｕｒｇｅｒｓｅｑｕａｔｉｏｎｓ ａｒｅｔａｋｅｎａｓｅｘａｍｐｌｅｔｏｅｘｐｌａｉｎｔｈｅａｐｐｌｉｃａｔｉｏｎｏｆｔｈｅｍｅｔｈｏｄ，ａｎｄｔｈｅｎｕｍｅｒｉｃａｌｒｅｓｕｌｔｓｉｎｄｉｃａｔｅｔｈｅｅｆ ｆｅｃｔｉｖｅｎｅｓｓｏｆｔｈｉｓｍｅｔｈｏｄ． Ｋｅｙｗｏｒｄｓ：ｇｒａｖｉｔｙｉｎｔｅｒｐｏｌａｔｉｏｎ；ｌｏｃａｌｄｉｆｆｅｒｅｎｔｉａｌｑｕａｄｒａｔｕｒｅｍｅｔｈｏｄ；Ｂｕｒｇｅｒｓｅｑｕａｔｉｏｎｓ 　　微分求积法（ＤｉｆｆｅｒｅｎｔｉａｌＱｕａｄｒａｔｕｒｅＭｅｔｈｏｄ，简 行，因而ＤＱ法只适用于规则区域，对于不规则区域 ［１］ 称ＤＱ法）是１９７１年Ｂｅｌｌｍａｎ和Ｃａｓｔｉ 提出的一种 无法直接应用； 采用均匀节点分布时节点过多会