二维、三维 Laplace 算子的极坐标表示

简介: (1) 设 $(r,\theta)$ 是 $\bbR^2$ 的极坐标, 即 $$\bex x=r\cos\theta,\quad y=r\sin \theta. \eex$$ 证明 Laplace 算子 $\dps{\lap=\frac{\p^2}{\p x^2}+ \frac{\p^2}{\p y...

(1) 设 $(r,\theta)$ 是 $\bbR^2$ 的极坐标, 即 $$\bex x=r\cos\theta,\quad y=r\sin \theta. \eex$$ 证明 Laplace 算子 $\dps{\lap=\frac{\p^2}{\p x^2}+ \frac{\p^2}{\p y^2}}$ 可以表示为 $$\bex \lap u=u_{rr}+\frac{1}{r}u_r+\frac{1}{r^2}u_{\theta\theta}. \eex$$ (2) 设 $(r,\theta,\phi)$ 是 $\bbR^3$ 的极坐标, 即 $$\bex x=r\sin \theta\cos \phi,\quad y=r\sin \theta\sin \phi,\quad z=r\cos \theta. \eex$$ 证明 Laplace 算子 $\dps{\lap=\frac{\p^2}{\p x^2}+ \frac{\p^2}{\p y^2}+\frac{\p^2}{\p z^2}}$ 可以表示为 $$\bex \lap u=\frac{1}{r^2}\frac{\p}{\p r}\sex{r^2\frac{\p u}{\p r}} +\frac{1}{r^2\sin \theta} \frac{\p}{\p \theta}\sex{\sin \theta\frac{\p u}{\p \theta}} +\frac{1}{r^2\sin^2\theta}\frac{\p^2u}{\p \phi^2}. \eex$$

证明: (1) 由 $$\bex x_r=\cos \theta,\quad x_\theta=-r\sin \theta=-y,\quad y_r=\sin \theta,\quad y_\theta=r\cos \theta=x \eex$$ 知 $$\beex \bea u_r&=u_x\cos \theta+u_y\sin \theta,\\ u_{rr}&=u_{xx}\cos^2\theta +2u_{xy}\sin \theta\cos \theta +u_{yy}\sin^2\theta,\\ u_\theta&=-yu_x+xu_y,\\ u_{\theta\theta}&= -xu_x-y(-yu_{xx}+xu_{xy})\\ &\quad-yu_y+x(-yu_{xy}+xu_{yy})\\ &=y^2u_{xx}-2xyu_{xy} +x^2u_{yy}-(xu_x+yu_y). \eea \eeex$$ 而 $$\bex u_{rr}+\frac{1}{r}u_r+\frac{1}{r^2}u_{\theta\theta} =u_{xx}+u_{yy}=\lap u. \eex$$ (2) 设 $\rho=r\sin \theta$, 则 $$\bex \ba{ll} x=\rho \cos \phi,&y=\rho\sin \phi,\\ z=r\cos \theta,&\rho=r\sin \theta. \ea \eex$$ 而由 (1), $$\beex \bea u_{xx}+u_{yy}&=u_{\rho\rho} +\frac{1}{\rho}u_\rho+\frac{1}{\rho^2}u_{\phi\phi},\\ u_{zz}+u_{\rho\rho} &=u_{rr}+\frac{1}{r}u_r+\frac{1}{r^2}u_{\theta\theta}. \eea \eeex$$ 两式相加得 $$\bee\label{3.2:1} \lap u=u_{rr}+\frac{1}{r}u_r+\frac{1}{\rho}u_\rho +\frac{1}{r^2}u_{\theta\theta}+\frac{1}{\rho^2}u_{\phi\phi}. \eee$$ 我们再计算 $u_\rho$ 如下: $$\bee\label{3.2:2} \bea u_\rho&=u_rr_\rho+u_\theta\theta_\rho\quad\sex{z=r\cos\theta,\ \rho=r\sin \theta,\quad u(z,\rho)=u(r,\theta)}\\ &=\frac{\rho}{r}u_r+\frac{z}{r^2}u_\theta\\ &=u_r\sin \theta+u_\theta\frac{\cos\theta}{r}. \eea \eee$$ 把 \eqref{3.2:2} 代入 \eqref{3.2:1}, 得 $$\beex \bea \lap u&=u_{rr}+\frac{1}{r}u_r+\frac{1}{r\sin\theta} \sex{u_r\sin \theta+u_\theta\frac{\cos \theta}{r}} +\frac{1}{r^2}u_{\theta\theta} +\frac{1}{r^2\sin^2\theta}u_{\phi\phi}\\ &=u_{rr}+\frac{2}{r}u_r+ \frac{1}{r^2\sin\theta}\cdot u_\theta\cos \theta +\frac{1}{r^2}u_{\theta\theta} +\frac{1}{r^2\sin^2\theta}u_{\phi\phi}\\ &=\frac{1}{r^2}(r^2u_{rr}+2ru_r) +\frac{1}{r^2\sin\theta} (u_\theta\cos \theta+u_{\theta\theta}\sin \theta) +\frac{1}{r^2\sin^2\theta}u_{\phi\phi}\\ &=\frac{1}{r^2}(r^2u_r)_r +\frac{1}{r^2\sin\theta}(u_\theta\sin \theta)_\theta +\frac{1}{r^2\sin^2\theta}u_{\phi\phi}. \eea \eeex$$

目录
相关文章
|
存储 Shell Linux
安卓手机上的终端模拟器Termux,实现了在手机上运行python
安卓手机上的终端模拟器Termux,可以在手机上运行python
2751 2
|
图形学
Unity 3D游戏-消消乐(三消类)教程和源码
Unity 消消乐教程和源码 本文提供全流程,中文翻译。Chinar坚持将简单的生活方式,带给世人!(拥有更好的阅读体验 —— 高分辨率用户请根据需求调整网页缩放比例) ...
6352 0
|
7月前
|
运维 Kubernetes API
解决Kubernetes集群中master节点无法与node节点通信的策略。
这些策略不仅需要执行命令来获取信息,更要深入理解集群组件如何交互,以便进行准确的故障定位与修复。一条一条地排查,并适时回顾配置文件,证书有效性等,通常可以找到问题所在。给出的命令需要根据具体环境的配置进行适当的修改。故障排除往往是一个细致且需求反复验证的过程,但遵循上述策略可以高效定位大部分通信故障的原因。
574 12
|
11月前
|
搜索推荐
课时10:sublime的基本设置
今天,我们来聊聊如何对SublimeText进行简单的个性化配置。在使用SublimeText的过程中,很多人都会遇到一些问题,比如Sublime自带的字体不太好看,或者字体大小不符合个人需求,不是偏大就是偏小。接下来,我们就详细看看如何调整这些设置。 1.字体大小与样式调整 2.主题安装与配置
1430 1
|
存储 缓存 小程序
微信小程序数据缓存与本地存储:优化用户体验
本文深入探讨微信小程序的数据缓存与本地存储,介绍其意义、机制及应用场景。通过合理使用内存和本地缓存,可减少网络请求、提升加载速度和用户体验。文中详细讲解了常用缓存API的使用方法,并通过一个新闻列表案例展示了缓存的实际应用。最后提醒开发者注意缓存大小限制、时效性和清理,以确保最佳性能。
|
编译器 C语言 计算机视觉
【ARM汇编速成】零基础入门汇编语言之指令集(二)
【ARM汇编速成】零基础入门汇编语言之指令集(二)
1434 0
|
监控 Java 数据处理
Spring Batch 是如何工作的?
Spring Batch 是如何工作的?
622 2
|
机器学习/深度学习 人工智能 数据可视化
斯坦福博士图解AlphaFold 3:超多细节+可视化还原ML工程师眼中的AF3
【8月更文挑战第8天】AlphaFold 3作为AI领域的重大突破,革新了蛋白质结构预测。斯坦福博士通过图解详析了其内部机制,展示了多尺度建模与图神经网络技术如何提升预测精度。尽管存在数据依赖性和计算成本等挑战,AlphaFold 3仍极大地加速了生物学研究与药物开发进程。论文详情参见:https://www.nature.com/articles/s41586-024-07487-w
781 4
|
存储 Python
数据包络分析(Data Envelopment Analysis, DEA)详解与Python代码示例
数据包络分析(Data Envelopment Analysis, DEA)详解与Python代码示例
|
缓存 Linux
linux性能分析之内存分析(free,vmstat,top,ps,pmap等工具使用介绍)
这些工具可以帮助你监视系统的内存使用情况、识别内存泄漏、找到高内存消耗的进程等。根据具体的问题和需求,你可以选择使用其中一个或多个工具来进行内存性能分析。注意,内存分析通常需要综合考虑多个指标和工具的输出,以便更好地理解系统的行为并采取相应的优化措施。
993 6