一个非线性偏微分方程边值问题的对称约化及其数值解 Symmetry Reduction and Its Numerical Solution to the Boundary Value Problem of a Nonlinear Partial Differential Equation
 

Advances in Applied Mathematics
Vol.05 No.03(2016), Article ID:18350,6 pages
10.12677/AAM.2016.53046

Symmetry Reduction and Its Numerical Solution to the Boundary Value Problem of a Nonlinear Partial Differential Equation

Yanqing Han, Bilige Sudao*

College of Science, Inner Mongolia University of Technology, Hohhot Inner Mongolia

Received: Jul. 27th, 2016; accepted: Aug. 15th, 2016; published: Aug. 18th, 2016

Copyright © 2016 by authors and Hans Publishers Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/

ABSTRACT

We study the applications of the symmetry method on the boundary value problem for nonlinear partial differential equation. Firstly, the multi-parameter symmetry of a given boundary value problem for nonlinear partial differential equation is determined based on differential characteristic set algorithm. Secondly, by using the symmetry, the boundary value problem for nonlinear partial differential equation is reduced to an initial value problem of the original differential equation. Finally, we numerically solve the initial value problem of the original differential equations by using Runge-Kutta method.

Keywords:Boundary Value Problem for Nonlinear Partial Differential Equation, Differential Characteristic Set Algorithm, Symmetry Method, Runge-Kutta Method

一个非线性偏微分方程边值问题的对称约化及其数值解

韩雁清,苏道毕力格*

内蒙古工业大学理学院,内蒙古 呼和浩特

收稿日期:2016年7月27日;录用日期:2016年8月15日;发布日期:2016年8月18日

摘 要

本文研究了微分方程对称方法在非线性偏微分方程边值问题中的应用。首先,基于微分特征列集算法确定了给定非线性偏微分方程边值问题的多参数对称;其次,利用对称将非线性偏微分方程边值问题化为常微分方程初值问题;最后,利用龙格-库塔法求解了常微分方程初值问题的数值解。

关键词 :非线性偏微分方程边值问题,微分特征列集算法,对称方法,龙格–库塔法

1. 引言

Lie对称是公认的普适性方法之一,它以诸多传统方法为其特例,如:分离变量法、行波变换、相似变换等 [1] 。对称理论在现代数学、物理和力学等科学中有重要的理论和实际意义,并且已有了广泛的应用 [2] [3] 。但是除了在文 [1] [4] [5] 中研究者做了一些对称群在边值问题的应用外,这方面的研究还很少,所以利用对称群研究偏微分方程(简记为PDEs)边值问题是对称理论应用新的研究领域。在诸多求解非线性的问题中,常常利用相似变换对PDEs进行约化或降阶。我们知道PDEs的Lie变换群可以产生更一般形式的相似变换,并且这些变换有更多的数学和物理意义,所以利用对称方法研究非线性PDEs边值问题比直接用相似变换更有其优越性。

最近我们基于微分特征列集算法 [6] [7] 研究了对称方法在非线性PDEs边值问题中的应用,并且将对称方法、对称分类与数值方法、同伦分析方法进行结合解决了非线性PDEs边值问题 [8] - [11] 。也有研究者基于微分特征列集算法,将对称方法与变分迭代法、同伦摄动法结合,求解了PDEs边值问题 [12] [13] 。本文中,我们将基于微分特征列集算法,有效结合对称方法和数值方法求解一个非线性PDEs边值问题。

2. 一个非线性PDEs边值问题的对称约化及其数值解

考虑一个如下非线性PDEs边值问题

(1)

其中为常数。相应的边界条件为

,若 (2)

,若 (3)

下面引入流函数

(4)

将上面的关系代入(1)~(3),得到

(5)

满足边界条件

,若 (6)

,若 (7)

其中函数在后面根据边界条件在对称作用下的不变性来确定。

2.1. 边值问题的对称约化

下面利用对称方法约化边值问题(5)~(7)。

第一步:确定PDE边值问题(5)~(7)的多参数对称。

假设非线性PDE(5)的对称的无穷小向量为

其中为该对称的无穷小生成函数。

1) 产生关于无穷小生成函数的确定方程组。

根据Lie算法,利用产生确定方程组的算法 [7] 可得到方程(5)的对称的确定方程组,但是难于人工求解。我们可以利用微分特征列集程序包进行计算得到与确定方程组等价的特征列集对应的方程组,即

(8)

2) 确定多参数对称。

通过求解方程组(8),得到无穷小生成函数

即可得到对称的无穷小向量为

(9)

其中为任意常数,为任意函数。

第二步:约化边值问题。利用对称(9)将边值问题(5)~(7)约化为常微分方程初值问题。

对称(9)的特征方程为

,得到不变量

同理可从特征方程得到

(10)

将(10)代入方程(5),可得到常微分方程

(11)

根据PDEs边值问题的边界条件在对称群作用下的不变性 [1] 知,边界条件(6),(7)在对称(9)的延拓作用下不变,故有

,当 (12)

,当 (13)

,当 (14)

其中为对称的一阶延拓的无穷小向量,为对称的二阶延拓的无穷小向量,即

由关系式(12),(13),(14)可以确定函数,即

其中为任意常数。

为了对应边界条件(6)和(7),取,故有

时,

根据边界条件(6),(7)和关系(10),得到初值条件

(15)

2.2. 边值问题的数值解

下面我们利用四阶龙格–库塔法求解常微分方程初值问题(11),(15)的数值解。

为了使用龙格–库塔法,首先将(11),(15)化为一阶常微分方程初值问题。令

则(11),(15)可以化为

(16)

相应的有初值条件

(17)

其中都是的函数,那么可以建立龙格–库塔公式

其中为步长。

常微分方程(16)中的常数和初始条件(17)中的取其不同值时,我们也可以得到相应的数值解。下面我们取,并且时,借助Mathematica符号系统可以得到函数在区间上的数值解,如图1所示。

Figure 1. Nomerical solution of in

图1.上的数值解

3. 结论

本文研究了微分方程对称方法在非线性PDE边值问题中的应用,也探索了利用对称方法与数值方法的有效结合求解PDE边值问题的新途径。我们首先基于微分特征列集算法确定了给定非线性PDE边值问题的多参数对称,并利用对称将非线性PDEs边值问题约化为常微分方程初值问题。其次我们充分考虑到对称方法与数值方法的彼此互补性,借助于Mathematica符号系统利用四阶龙格-库塔法求解了上一步得到的常微分方程初值问题的数值解。得到的结果充分体现了对称方法在偏微分方程应用中的优越性,并且我们在本文中有效结合对称方法和数值方法求解了边值问题,这项研究也推广了对称方法的应用范围。

基金项目

本文由国家自然科学基金项目(11571008),内蒙古自治区自然科学基金(2014MS0114)资助。

文章引用

韩雁清,苏道毕力格. 一个非线性偏微分方程边值问题的对称约化及其数值解
Symmetry Reduction and Its Numerical Solution to the Boundary Value Problem of a Nonlinear Partial Differential Equation[J]. 应用数学进展, 2016, 05(03): 375-380. http://dx.doi.org/10.12677/AAM.2016.53046

参考文献 (References)

  1. 1. Bluman, G.W. and Kumei, S. (1989) Symmetries and Differential Equations. Spring-Verlag, New York, Berlin.

  2. 2. Bluman, G., Cheviakov, A.F. and Anco, S.C. (2010) Applications of Symmetry Methods to Partial Differential. Spring-Verlag, New York.

  3. 3. Ibragimov, N.H. and Ibragimov, R.N. (2012) Applications of Lie Group Analysis to Mathematical Modelling in Natural Sciences. Mathematical Modelling of Natural Phenomena, 7, 52-65. http://dx.doi.org/10.1051/mmnp/20127205

  4. 4. Seshadri, R. and Na, T.Y. (1985) Group Invariance in Engineering Boundary Value Problems. Springer-Verlag, New York-Berlin-Heidelberg Tokyo.

  5. 5. Yürüsoy M,Pakdemirli M,Noyan O F. (2001) Lie Group Analysis of Creeping Flow of a Second Grade Fluid. Interna-tional Journal of Non-Linear Mechanics, 36, 955-960. http://dx.doi.org/10.1016/S0020-7462(00)00060-3

  6. 6. 朝鲁. 微分方程(组)对称向量的吴–微分特征列集算法及其应用[J]. 数学物理学报, 1999, 19(3): 326-332.

  7. 7. 特木尔朝鲁, 白玉山. 基于吴方法的确定和分类(偏)微分方程古典和非古典对称新算法理论[J]. 中国科学: A辑, 2010, 40(4): 1-18.

  8. 8. 苏道毕力格. 一些求解偏微分方程解析解方法的研究[D]: [博士学位论文]. 呼和浩特: 内蒙古工业大学, 2011.

  9. 9. 王晓民, 苏道毕力格, 特木尔朝鲁. 对称方法在非线性偏微分方程边值问题中的应用[J]. 内蒙古大学学报(自然科学版), 2013,44(2): 129-132.

  10. 10. 苏道毕力格, 王晓民, 乌云莫日根. 对称分类在非线性偏微分方程组边值问题中的应用[J]. 物理学报, 2014, 63(4): 040201.

  11. 11. 苏道毕力格, 王晓民, 鲍春玲. 利用对称方法求解非线性偏微分方程组边值问题的数值解[J]. 应用数学, 2014, 27(4): 10-15.

  12. 12. Lu, L. and Temuer, C. (2011) A New Method for Solving Boundary Value Problems for Partial Differential Equations. Computers and Mathematics with Applications, 61, 2164-2167. http://dx.doi.org/10.1016/j.camwa.2010.09.002

  13. 13. Eerdunbuhe and Temuerchaolu (2012) Approximate Solution of the Magneto-Hydrodynamic Flow over a Nonlinear Stretching Sheet. Chinese Physics B, 21, 035201. http://dx.doi.org/10.1088/1674-1056/21/3/035201

*通讯作者。

期刊菜单