理论教育 二维偏微分方程的常微分方程求解法

二维偏微分方程的常微分方程求解法

时间:2023-10-26 理论教育 版权反馈
【摘要】:在某些情况下,采用加权余量法可将两个自变量的偏微分方程定解问题化成常微分方程定解问题求解,现以二维泊松方程定解问题来说明。设定解问题为二维定解问题的解u=u(x,y)依赖于变量x和y,若取基函数仅为一个变量x的函数φ0,φ1,…

二维偏微分方程的常微分方程求解法

在某些情况下,采用加权余量法可将两个变量偏微分方程定解问题化成常微分方程定解问题求解,现以二维泊松方程定解问题来说明。设定解问题为

二维定解问题的解u=uxy)依赖于变量xy,若取基函数仅为一个变量x的函数φ0x),φ1x),…,φNx),则在构造解的近似函数时,系数αi必然是y的函数,即

若取基函数为两个变量xy的函数φ0xy),φ1xy),…,φNxy),则在构造解的近似函数时,系数αi可以视为常数,即

以上两种近似函数取法都是可行的,我们先采用后一种方法选取近似函数。根据定解问题的边界条件,取近似函数为

uxy;α)=1001φ1xy

这里假定φ0=100,以满足问题给定的边界条件。在选取φ1xy)时,为了使之能满足齐次边界条件,即

φ1xy)|x=±51xy)|y=±5=0

可取φ1xy)为

φ1xy=(25-x2)(25-y2

于是,得近似函数为

uxy;α)=1001(25-x2)(25-y2

可得余量为

R=-2α1(50-x2-y2+2

采用伽辽金法,可得方程

解得

于是得到近似解

下面采用前一种方法选取近似函数。根据定解问题的边界条件,取近似函数为

uxy;α)=1001yφ1x

为了让φ1满足问题给定的边界条件。

φ1x)|x=±5=0(www.daowen.com)

可取

φ1x=(25-x2

于是,得近似函数为

uxyα=100+(25-x2α1y

为使其满足y=±5处的边界条件,可推得

α1(5)=0,α1(-5)=0

此时,余量为

R=-2α1+(25-x2α″1+2

采用伽辽金法,可得方程

解得

10α″1-α1+1=0

这是一个常微分方程,考虑到α1(5)=α1(-5)=0,不难解得

得到近似解为

为比较上述两种结果,取计算域中心,即x=y=0处的值,可得两个近似解分别为

u(0,0)=115.63

u(0,0)=115.13

该定解问题的精确解为

u(0,0)=114.65可见两个近似解的误差均在1%以内。

免责声明:以上内容源自网络,版权归原作者所有,如有侵犯您的原创版权请告知,我们将尽快删除相关内容。

我要反馈