调和方程(拉普拉斯方程)基本解和边界元方法的积分计算

前言

调和方程的基本解以及边界元方法中某积分的解析结果。

调和方程

设 $u(x_1,...,x_n)=f(r)$ （其中 $r=\sqrt{x_1^2+...+x_n^2}$ ）是 $n$ 维调和函数 （即满足方程 $\frac{{\partial }^{2}u}{\partial x_1^2}+...,+\frac{{\partial }^{2}u}{\partial x_n^2}=0$）,证明如下等式成立
$f(r)=c_1+\frac{c_2}{r^{n-2}}(n\ne 2)$
$f(r)=c_1+c_{2}ln\frac{1}{r}(n=2)$
其中$c_1,c_2$ 为任意常数。
证明：
$$u=f(r),\frac{\partial u}{\partial x_i}=f^{{}^{\prime }}(r)\cdot \frac{\partial r}{\partial x_i}=f^{{}^{\prime }}(r)\cdot \frac{x_i}{r}$$

$$\frac{{\partial }^{2}u}{\partial x_i^2}=f^{{}^{\prime \prime }}(r)\cdot \frac{x_i^2}{r^2}+f^{{}^{\prime }}(r)\cdot \frac{1}{r}-f^{{}^{\prime }}(r)\cdot \frac{x_i^2}{r^3}$$

$$\sum _{i=1}^n\frac{{\partial }^{2}u}{\partial x_i^2}=f^{{}^{\prime \prime }}(r)\cdot \frac{{\sum }_{i=1}^n{x}_i^2}{r^2}+f^{{}^{\prime }}(r)\cdot \frac{n}{r}-f^{{}^{\prime }}(r)\cdot \frac{{\sum }_{i=1}^n{x}_i^2}{r^3}=f^{{}^{\prime \prime }}(r)+\frac{n-1}{r}{f}^{{}^{\prime }}(r)$$
即方程 $\Delta u=0$ 化为 $f^{{}^{\prime \prime }}(r)+\frac{n-1}{r}{f}^{{}^{\prime }}(r)=0$

$$\frac{f^{{}^{\prime \prime }}(r)}{f^{{}^{\prime }}(r)}=-\frac{n-1}{r}(1)$$

所以

$$f^{{}^{\prime }}(r)=A_1{r}^{-(n-1)}(2)$$

若 $n\ne 2$，积分得

$$f(r)=\frac{A_1}{-n+2}{r}^{-n+2}+c_1$$

即 $n\ne 2$, 则 $f(r)=c_1+\frac{c_2}{r^{n-2}}$

若 $n=2$，则 $f^{{}^{\prime }}(r)=\frac{A_1}{r}$ 故 $f(r)=c_1+A_{1}lnr$

即$n=2$，则 $f(r)=c_1+c_{2}ln\frac{1}{r}$

上面的证明中 (1) 到 (2) 的推导如下：
可令 $z=f^{{}^{\prime }}(r)$， 则等式 (1) 变成 $\frac{z^{{}^{\prime }}}{z}=-\frac{n-1}{r}$

下面就类似求常微分方程 $\frac{dy}{dx}=P(x)y$ 的通解。
变量分离法：
$$\frac{dy}{y}=P(x)dx$$
两边积分，即得
$$ln\left|y\right|=\int P(x)dx+c_1$$
这里 $c_1$ 是任意常数，由对数定义，即有
$$\left|y\right|=e^{\int P(x)dx+c_1}$$
即
$$y=\pm e^{c_1}\cdot e^{\int P(x)dx}$$
令 $\pm e^{c_1}=c$ 得到
$$y=c{e}^{\int P(x)dx}$$

边界元方法的基本知识

二维空间的拉普拉斯方程的基本解为

$$u_s^{\ast }(r)=-\frac{1}{2\pi }lnr$$

$$q_s^{\ast }(r)=\frac{\partial u_s^{\ast }}{\partial \vec{n}}=-\frac{(\vec{r},\vec{n})}{2\pi r^2}$$

三维空间的拉普拉斯方程的基本解为

$$u_s^{\ast }(r)=\frac{1}{4\pi r}$$

$$q_s^{\ast }(r)=\frac{\partial u_s^{\ast }}{\partial \vec{n}}=-\frac{(\vec{r},\vec{n})}{4\pi r^3}$$

积分计算

如下图所示计算积分

$$I={\int }_L\frac{\vec{r}\cdot \vec{n}}{r^2}dl$$

其中 $\vec{r}=(x,y),L=\overline{AB}$

经计算可得
$$I=\beta ,其中\beta 为线段OA与线段OB的夹角的弧度值$$

注：可用如下网址验证 Integral Calculator