计算物理 ›› 2019, Vol. 36 ›› Issue (5): 517-532.DOI: 10.19596/j.cnki.1001-246x.7919

• • 上一篇    下一篇

求解含有高阶导数偏微分方程的局部间断Petrov-Galerkin方法

赵国忠1, 蔚喜军2, 郭虹平1, 董自明1   

  1. 1. 包头师范学院 数学科学学院, 包头 014030;
    2. 北京应用物理与计算数学研究所 计算物理实验室, 北京 100088
  • 收稿日期:2018-07-04 修回日期:2018-09-22 出版日期:2019-09-25 发布日期:2019-09-25
  • 通讯作者: Liu Haifeng (1968-),condensed matter physics,E-mail:liu_haifeng@iapcm.ac.cn
  • 作者简介:Zhao Guozhong (1977-),male,PhD,professor,research in computational fluid dynamics,E-mail:zhaoguozhongbttc@sina.com
  • 基金资助:
    National Natural Science Foundation of China (11761054, 11571002, 11261035), Program for Young Talents of Science and Technology in Universities of Inner Mongolia Autonomous Region (NJYT-15-A07), Natural Science Foundation of Inner Mongolia Autonomous Region, China (2015MS0108, 2012MS0102) and Science Research Foundation of Institute of Higher Education of Inner Mongolia Autonomous Region, China (NJZZ12198, NJZZ16234, NJZZ16235)

A Local Discontinuous Petrov-Galerkin Method for Partial Differential Equations with High Order Derivatives

ZHAO Guozhong1, YU Xijun2, GUO Hongping1, DONG Ziming1   

  1. 1. Faculty of Mathematics, Baotou Teachers'College, Baotou 014030, China;
    2. Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, Beijing 100088, China
  • Received:2018-07-04 Revised:2018-09-22 Online:2019-09-25 Published:2019-09-25
  • Supported by:
    National Natural Science Foundation of China (11761054, 11571002, 11261035), Program for Young Talents of Science and Technology in Universities of Inner Mongolia Autonomous Region (NJYT-15-A07), Natural Science Foundation of Inner Mongolia Autonomous Region, China (2015MS0108, 2012MS0102) and Science Research Foundation of Institute of Higher Education of Inner Mongolia Autonomous Region, China (NJZZ12198, NJZZ16234, NJZZ16235)

摘要: 构造一类求解三种类型偏微分方程的间断Petrov-Galerkin方法.求解的方程分别含有二阶、三阶和四阶偏导数,包括Burgers型方程、KdV型方程和双调和型方程.首先将高阶微分方程转化成为与之等价的一阶微分方程组,再将求解双曲守恒律的间断Petrov-Galerkin方法用于求解微分方程组.该方法具有四阶精度且具有间断Petrov-Galerkin方法的优点.数值实验表明该方法可以达到最优收敛阶而且可以模拟复杂波形相互作用,如孤立子的传播及相互碰撞等.

关键词: KdV型方程, 双调和方程, 局部间断Petrov-Galerkin方法, 孤立子演化

Abstract: A local discontinuous Petrov-Galerkin method is proposed for solving three types of partial differential equations with second, third and fourth order derivatives, respectively. They are Burgers type equations, KdV type equations and bi-harmonic type equations. The method extends discontinuous Petrov-Galerkin method for conservation laws by rewriting corresponding equations into a first order system and solving the system instead of the original equation. The method has a fourth order accuracy and maintains advantages of discontinuous Petrov-Galerkin method. Numerical simulations verify that the method reaches optimal convergence order and simulates well complex wave interaction such as soliton propagation and collision.

Key words: KdV type equation, bi-harmonic equation, local discontinuous Petrov-Galerkin method, soliton evolution

中图分类号: