理论教育 强解积分表达式-偏微分方程数值解法

强解积分表达式-偏微分方程数值解法

时间:2023-10-26 理论教育 版权反馈
【摘要】:若令方程余量和自然边界上余量加权积分的线性组合为零,则由此产生的积分表达式,称为强解积分表达式。在应用伽辽金法求解问题时,由式和式表示的两种积分表达式都可以采用,此时自然边界条件在表达中自动满足,只要求近似解满足强加边界条件,因弱解积分表达式比较简单,所以经常取弱解积分表达式进行计算。图2.5.1 柱坐标系一维边值问题先写出伽辽金积分表达式进行分部积分,并注意到在边界上δu=0,因此有这就是伽辽金弱解表达式。

强解积分表达式-偏微分方程数值解法

若令方程余量和自然边界上余量加权积分的线性组合为零,则由此产生的积分表达式,称为强解积分表达式。一般情况下,方程余量加权积分和边界余量的加权积分的线性组合都采取最简单的线性组合形式,即相加或相减,取相加还是相减取决于数学或物理上的方便与否。数学上的考虑是指方程余量的加权积分经分部积分运算以后,可与边界余量的加权积分的一部分相抵消,物理上则反映求解区域内部物理量只能按某种方式相结合使之达到平衡。仍以上面的二阶常微分方程的边值问题式(2.5.1)~式(2.5.3)为例。强解积分表达式为

978-7-111-44528-9-Chapter02-255.jpg

之所以将方程余量加权积分减去边界余量加权积分,是为了便于式(2.5.7)的简化。事实上,对式(2.5.7)进行分部积分,并注意到δux=0=0,有

978-7-111-44528-9-Chapter02-256.jpg

从形式上看,由强解表达式所导出的式(2.5.8)和弱解表达式(2.5.5)完全相同,它们的区别主要反映在强解积分表达式(2.5.8)是伽辽金法求解的出发点,对近似解

978-7-111-44528-9-Chapter02-257.jpg

中基函数的要求比弱解积分表达式的要强一些。按式(2.5.7)的要求,基函数φk应具有二阶导数可积的光滑程度,这比弱解所要求的一阶导数平方可积的要求要强一些。

在应用伽辽金法求解问题时,由式(2.5.5)和式(2.5.8)表示的两种积分表达式都可以采用,此时自然边界条件在表达中自动满足,只要求近似解满足强加边界条件,因弱解积分表达式比较简单,所以经常取弱解积分表达式进行计算。例如,对于平面渗流问题

978-7-111-44528-9-Chapter02-258.jpg

式(2.5.9)第二式为强加边界条件,第三式为自然边界条件。

为了得到弱解积分表达式,我们从伽辽金积分表达式出发,利用高斯公式,有

978-7-111-44528-9-Chapter02-259.jpg

注意到在Γ1978-7-111-44528-9-Chapter02-260.jpg,δφ=0;在Γ2978-7-111-44528-9-Chapter02-261.jpg,于是,有

978-7-111-44528-9-Chapter02-262.jpg

或写成

978-7-111-44528-9-Chapter02-263.jpg

此即伽辽金弱解的积分表达式。

强解的积分表达式可直接写出

978-7-111-44528-9-Chapter02-264.jpg

我们讨论下面的二阶常微分方程边值问题的求解

978-7-111-44528-9-Chapter02-265.jpg

伽辽金积分表达式为

978-7-111-44528-9-Chapter02-266.jpg

对式(2.5.13)进行分部积分,得到弱解积分表达式

978-7-111-44528-9-Chapter02-267.jpg

取基函数φ1=xφ2=x2,…,取二级近似,令

u2=α1x+α2x2

它显然满足强加边界条件式(2.5.11)。将近似函数代入弱解积分表达式,得

978-7-111-44528-9-Chapter02-268.jpg

积分后,得

978-7-111-44528-9-Chapter02-269.jpg

该二元一次方程组解为

978-7-111-44528-9-Chapter02-270.jpg

于是,得到近似解

978-7-111-44528-9-Chapter02-271.jpg

若直接采用伽辽金积分表达式,则此时对近似函数要求比较苛刻,它必须同时满足强加边界条件式(2.5.11)及自然边界条件式(2.5.12)。对于本例这个较简单的情况,可以较容易地选取近似解

978-7-111-44528-9-Chapter02-272.jpg

它显然符合两个给定的边界条件。代入积分表达式,有(www.daowen.com)

978-7-111-44528-9-Chapter02-273.jpg

积分后,得

978-7-111-44528-9-Chapter02-274.jpg

于是,得

978-7-111-44528-9-Chapter02-275.jpg

近似解为

978-7-111-44528-9-Chapter02-276.jpg

本定解问题的解析解为

978-7-111-44528-9-Chapter02-277.jpg

解析解以及强解和弱解的计算结果列于表2.5.1。

2.5.1 解析解以及强解和弱解的计算结果

978-7-111-44528-9-Chapter02-278.jpg

下面我们讨论柱坐标系中一维导热拉普拉斯方程边值问题,如图2.5.1所示。定解条件为

978-7-111-44528-9-Chapter02-279.jpg

这是一个轴对称圆管壁内稳态热传导方程,u为温度,内壁面和外壁面分别给定温度u1u2。一般地,边界上也可给定热流量,即给出978-7-111-44528-9-Chapter02-280.jpg的值。给定温度的边界条件是强加边界条件,给定热流量的是自然边界条件。

978-7-111-44528-9-Chapter02-281.jpg

图2.5.1 柱坐标系一维边值问题

先写出伽辽金积分表达式

978-7-111-44528-9-Chapter02-282.jpg

进行分部积分,并注意到在边界上δu=0,因此有

978-7-111-44528-9-Chapter02-283.jpg

这就是伽辽金弱解表达式。

设近似解为

u=u0+α1φ1

其中,u0应是满足非齐次边界条件的特解,假定它是一个线性函数

978-7-111-44528-9-Chapter02-284.jpg

φ1是基函数,它应满足齐次边界条件,即在r=r1r=r2处,φ1=0,符合这个条件的最简单的多项式函数应是

φ1=r1r2-(r1+r2r+r2

于是,近似解成为

978-7-111-44528-9-Chapter02-285.jpg

其中,α1是待定系数。将上式代入伽辽金弱解表达式中,可求得

978-7-111-44528-9-Chapter02-286.jpg

于是得到方程的一级近似解,即

978-7-111-44528-9-Chapter02-287.jpg

本问题是传热学中的经典问题,其解析解为

978-7-111-44528-9-Chapter02-288.jpg

请读者通过预设u1u2r1r2的值,比较一级近似解u与解析解u∗之间的误差。

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

我要反馈