找回密码
 快速注册
搜索
查看: 65|回复: 3

泊松方程

[复制链接]

3

主题

4

回帖

39

积分

积分
39

显示全部楼层

edward076923 发表于 2023-3-8 08:58 |阅读模式
最近在看一些泊松方程有关的材料,非常好奇如何推导电场和引力场的泊松方程,进而形象理解高斯定理?

48

主题

992

回帖

1万

积分

积分
14981
QQ

显示全部楼层

Czhang271828 发表于 2023-3-8 15:02
其实以下例子足以帮助形象地理解 Gauss 定理, 不需要知道电场磁场等预备知识, 只要知道哈士奇是一种有腿的狗即可.

一个不透光的袋子里有很多哈士奇, 每个哈士奇都会把若干条腿撑破袋子露在外面. Gauss 定理描述了这样一个恒等式
\[
\sum_{\text{🐕}\in\text{袋子}} \text{🐕露出的腿数}=\sum_{\text{袋子表面}}\text{单位表面的腿数}.
\]
假设哈士奇和猫一样是液体, 那么就没有''一只哈士奇''这样的说法; 而应该说哈士奇在袋内的密度函数, 哈士奇露在袋子表面的腿也是以密度函数的方式定义. 令
\begin{align*}
\rho_\text{🐕}(x)&:=\lim_{U_x\to \{x\}}\dfrac{U_x\text{ 中狗露出的腿数}}{U_x\text{ 的体积}}&&\text{$U_x$ 为 $x$ 的邻域(袋内三维拓扑空间)},\\
\tilde E_\text{🐕}(S)&:=\lim_{U_S\to \{S\}}\dfrac{U_S\text{ 上狗腿数}}{U_S\text{ 的面积}}&&\text{$U_S$ 为 $S$ 的邻域(袋面二维拓扑空间)}.
\end{align*}
数, 那么就有
\[
\int_{\text{袋子}}\rho_{\text{🐕}}(x)\mathrm{d}x=\int_{\text{袋面}}\tilde E_{\text{🐕}}(S)\mathrm dS.
\]
这样的式子确实正确, 但右式有一处挺别扭的地方: 袋面上单位 $\Delta S$ 处的狗腿密度 $\tilde E_\text{🐕}(S)$ 和该处袋子的倾斜角有关, 但 $\Delta S$ 上的狗腿数是由哈士奇决定的, 而非袋表面的倾斜角度. 因此我们用矢量 $\mathbf E(S)$ 表示哈士奇在袋表 $S$ 处贡献的腿矢量, 以及用矢量 $\mathbf S$ 表示袋表面 $S$ 处法相. 那么就有标准意义的 Gauss 公式
\[
\int_{\text{袋子}}\rho_{\text{🐕}}(x)\mathrm{d}x=\int_{\text{袋面}}\mathbf E_{\text{🐕}}(S)\mathbf{dS}.
\]

如果要具体地研究 Poisson 方程, 不同专业不同方向都有不同的记号. 如果不是为了做研究, 看看哈士奇就够了.

3

主题

4

回帖

39

积分

积分
39

显示全部楼层

 楼主| edward076923 发表于 2023-3-11 05:54
本帖最后由 edward076923 于 2023-3-11 07:13 编辑 感谢大佬指导,不过有一件事我有点没明白,就是您假设了每个在里面的狗狗一定会有固定的腿露出来,也就是Div f的空间积分正比于M或者Q。这是永远成立的吗还是只对inverse square relation成立?

PS:我不太理解为什么$\operatorname{div}\left(e_r\over r^2\right)=0$ 但是 $\operatorname{div}\left(e_r\over r\right)\neq 0$

48

主题

992

回帖

1万

积分

积分
14981
QQ

显示全部楼层

Czhang271828 发表于 2023-3-11 20:04
edward076923 发表于 2023-3-11 05:54
感谢大佬指导,不过有一件事我有点没明白,就是您假设了每个在里面的狗狗一定会有固定的腿露出来,也就是Di ...

第一个问题: 哈士奇模型的本质是"两种计数腿数量的方法结果相同", 推广就是连续化的情形. 如果某点处的狗在某个很小的空间角上伸出的腿连续均匀, 那确实蕴含三维空间中的平方反比定理(球表面积平方正比).
实际计算中不用刻意考虑这些, 只要套公式那就是对的. 如果某个连续可微的向量场在某个开集上处处不遵守平方反比定律, 那只能说明那个开集上有稠密的哈士奇在捣乱, 从而不能看作点电荷模型之类的.

第二个问题: 我们知道散度的定义是
\[
\mathrm{div}(\vec F(x)):=\lim_{U_x\to \{x\}}\,\,|V(U_x)|^{-1}\cdot \bigcirc\!\!\!\!\!\!\!\! \int\!\!\!\!\!\int_{\partial U_x}\vec F(S)\cdot \vec n_S\,\mathrm dS.
\]
此处 $U_x$ 是 $x$ 的开邻域(等价地, 看作开球无妨), $|V(U_x)|$ 是 $U_x$ 的体积, $\vec F$ 是向量场, $\partial U_x$ 是开邻域表面, $\vec n_S$ 是 $\partial U_x$ 表面 $S$ 处的法向量. 有时把 $\vec n_S\,\mathrm dS$ 写作向量 $\mathbf{dS}$.

关于 $\mathrm{div}(e_r\cdot r^{-2})=0$ 的实际意义, 我们考虑物理意义: 假定 $\vec r\neq 0$, 介质均匀, 则点电荷在 $\vec r$ 处诱导的电场的散度就是 $\mathrm{div}(e_r\cdot r^{-2})=0$ 的常数倍. 结合物理意义考虑散度的定义式, 把 $U_x$ 想象成均匀带电体, $\mathrm{div}(e_r\cdot r^{-2})=0$ (在相差常数的意义下)即"电场在均匀带电体 $U_x$ 表面创造的撕裂/挤压力再除以 $|V(U_x)|$". 在 $U_x\to \{x\}$ 时自然趋向 $0$. 这是因为点电荷电场在原点外连续可微, 其局部为匀强电场.  

补充一个数学计算, 考虑向量场 $\vec F(x,y,z)=(F_{(x)},F_{(y)},F_{(z)})=e_r\cdot r^{-2}$, 那么
\begin{align*}
\mathrm{div}\vec F(r)&=\nabla\cdot \vec F(r)\\
&=\dfrac{\partial }{\partial x}F_{(x)}+\dfrac{\partial }{\partial y}F_{(y)}+\dfrac{\partial }{\partial z}F_{(z)}\\
&=\dfrac{\partial }{\partial x}\dfrac{x}{r^3}+\dfrac{\partial }{\partial y}\dfrac{y}{r^3}+\dfrac{\partial }{\partial z}\dfrac{z}{r^3}\\
&=\dfrac{1}{r^3}-\dfrac{3x^2}{r^5}+\dfrac{1}{r^3}-\dfrac{3y^2}{r^5}+\dfrac{1}{r^3}-\dfrac{3z^2}{r^5}\\
&=0.
\end{align*}
此处 $\dfrac{\partial r}{\partial x}=\dfrac{x}{r}$.
反过来, 如果 $3$ 维空间的径向对称函数决定向量场, 根据量纲分析原则, 其形式一定是 $\vec r\cdot f(r)=(xf(r),yf(r),zf(r))$. 作为能量守恒定律的推论, 如果向量场连续可微, 则一定局部匀强, 即 $\mathrm{div}(\vec r\cdot f(r))=0$. 列散度方程
\[
0=\nabla\cdot (\vec r\cdot f(r))=\sum_{\text{轮换}}\dfrac{\partial xf(r)}{\partial x}=\sum_{\text{轮换}}\left[f(r)+\dfrac{x^2 f'(r)}r\right].
\]
即微分方程 $3f(r)+rf'(r)=0$. 此处看出是 $(r^3f(r))'=0$, 即 $f(r)$ 是 $r^{-3}$ 的常数倍.

结论: 能量守恒定律(推论是局部匀强电场)确实表明平方反比定律. 根据上述式子推算, $n+1$ 维空间中是 $n$-次方反比, $1$-维空间就是对数正比.
无钱佮歹看、无样佮歹生、无汉草佮无文采、无学历佮无能力、无高度无速度无力度共闲无代志。(闽南话)
口号:疼惜生命,远离内卷。

手机版|悠闲数学娱乐论坛(第3版)

GMT+8, 2025-3-4 16:35

Powered by Discuz!

× 快速回复 返回顶部 返回列表