理论教育 2D和3D问题的欧拉方程解法-偏微分方程数值解法:土建类

2D和3D问题的欧拉方程解法-偏微分方程数值解法:土建类

时间:2023-10-26 理论教育 版权反馈
【摘要】:下面我们直接从变分的角度讨论一般二维和三维边值问题的欧拉方程。于是,有考虑到δu在域D内和边界Γ2上都是任意的,要使泛函的变分δJ=0,必须使得下面的两式成立此即二维问题的欧拉方程。对于泛函应用推导二维问题欧拉方程相同的方法,可导出三维问题的欧拉方程。因此,在三维问题中,由式表示的泛函J,如果在u=u上取得极值,则函数u=u必须满足欧拉方程式。图2.2.3 三维边值问题图2.2.4 可变域变分

2D和3D问题的欧拉方程解法-偏微分方程数值解法:土建类

前面讨论的泛函J[yx)],其中y是只包括一个变量x的函数,本节先讨论的泛函J[uxy)],其中函数u是两个自变量xy的函数。对于泛函

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

式中,xy为两个独立变量,是平面两个坐标;uxy)是任意连续平面函数,该函数在区域D的边界Г上有已知值,即具有第一类边界条件,如图2.2.1所示的曲面,在这种情况下求泛函取得极值的曲面uxy)。

因已知uxy)|Г=f,因此这是一个固定边界的变分问题,即泛函J取极值时,极值曲面的一系列邻近曲面在区域D中的值可改变,但在边界上的值不变。与前面介绍的固定端点变分类似,作极值曲面uxy)的一系列邻近曲面

uxyε)=uxy)+εηxy) (2.2.13)

式中,ηxy)为在边界Г上函数值等于零的任意光滑连续函数,即ηxy)|Г=0。ε为一绝对值较小的参变量,并与xy无关。于是

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

图2.2.1 边界值固定的平面场

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

将这些邻近曲面代入式(2.2.12),得

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

泛函取极值的条件为978-7-111-44528-9-Chapter02-109.jpg,故

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

εxy无关,所以978-7-111-44528-9-Chapter02-111.jpg,上式简化为

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

为进一步简化式(2.2.14),可应用格林公式。格林公式的形式为

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

式中,Xxy)和Yxy)是定义在区域D和边界Γ上的两个任意连续函数。取出式(2.2.14)中的后两项

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

对上式右端的第一个积分式应用格林公式,令978-7-111-44528-9-Chapter02-115.jpg978-7-111-44528-9-Chapter02-116.jpg,得

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

代回式(2.2.14),得

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

对于固定边界ηxy)|Γ=0,式(2.2.16)成为

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

式(2.2.17)中积分值对任意函数ηxy)都为零,于是,极值曲面uxy)必须满足的条件为

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

式(2.2.18)即为欧拉方程。

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

由于xy为独立变量,因此上式可以简写为

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

同理得

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

式(2.2.18)是极值曲面必须满足的微分方程,函数在边界Γ上的已知值作为欧拉微分方程的定解条件。

下面考虑一个典型的二维边值问题,求泛函

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

的极值曲面为uxy),其中,曲面u在域D的边界Γ上的温度值为已知,即

uxy)|Γ=fxy

应用欧拉方程式(2.2.18),其中

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

注意到F中不包含xyu,代入式(2.2.19)和式(2.2.20),得

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

因此

978-7-111-44528-9-Chapter02-127.jpg(www.daowen.com)

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

于是

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

代入欧拉方程式(2.2.18),得

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

可见,极值曲面为满足拉普拉斯方程的一系列曲面,为了确定唯一的极值曲面,必须利用固定边界给出的定解条件,这里指定为第一类边界条件,即强加边界条件。

下面我们直接从变分的角度讨论一般二维和三维边值问题的欧拉方程。考虑二维边值问题,设有泛函

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

其中,uxy),Gxyu)是沿边界Γ2取值的函数。这里假定u在边界Γ1上的函数值是给定的,在边界Γ2上是未定的,如图2.2.2所示。

假定泛函J[u]存在,并在u=uxy)上实现极值,在极值曲面uxy)的邻域,函数可写成

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

在边界Γ1上,有δuxy)=0。于是,式(2.2.21)成为

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

图2.2.2 二维边值问题

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

上式的泛函为参数ε的函数。由积分符号下求微分的方法,得到

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

ε=0时,有

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

再由部分积分

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

于是

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

令上式右边第二项中的978-7-111-44528-9-Chapter02-139.jpg978-7-111-44528-9-Chapter02-140.jpg,根据格林公式

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

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

其中,Γ=Γ1+Γ2。考虑到边界Γ1上的δu=0,上式右边只剩下沿边界Γ2上的线积分。在边界上dy=nxdΓ,dx=-nydΓ。这里nx=cos(n,x),ny=cos(n,y),为边界外法线的方向余弦。于是,有

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

考虑到δu在域D内和边界Γ2上都是任意的,要使泛函的变分δJ=0,必须使得下面的两式成立

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

此即二维问题的欧拉方程。因此,由式(2.2.21)所给的泛函,如在u=uxy)上取得极值,函数uxy)必须满足欧拉方程式(2.2.22)。

用同样的方法可以推导出三维的欧拉方程。设函数u=uxyz)为空间域Dxyz)内的连续二阶可微函数,u在边界Γ1u=u1为已知的,在边界Γ2上为未知的,如图2.2.3所示。对于泛函

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

应用推导二维问题欧拉方程相同的方法,可导出三维问题的欧拉方程。

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

式中,

nx=cos(n,x),ny=cos(n,y),nz=cos(n,z)为边界外法线的方向余弦。因此,在三维问题中,由式(2.2.23)表示的泛函Ju),如果在u=uxyz)上取得极值,则函数u=uxyz)必须满足欧拉方程式(2.2.24)。

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

图2.2.3 三维边值问题

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

图2.2.4 可变域变分

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

我要反馈