Sommerfeld radiation condition

Arnold Johannes Wilhelm Sommerfeld (1868–1951)

复习

• curl 一个向量场的Curl 描述了一个三维的向量场$F$的无穷小旋转（infinitesimal rotation），有时也用$\nabla \times F,rot F$来表示。直观上，如果向量场表示流体的速度，那么$curl F$描述的就是循环密度（circulation density of the fluid）,$curl F=0$的场被称为无旋的。分析上的定义 $$(\nabla \times F)\cdot \hat{n}:=\lim_{A\to 0}{\frac{1}{|A|}\oint_{\partial A}F\cdot dr}$$ 直观上，举个简单例子:取 $\displaystyle F=y\hat{x}-x\hat{y}$，它的旋只有$\hat{z}$方向有分量，可以轻易算出来$\displaystyle \frac{2\pi r^2}{\pi r^2}=2$,直观上就是单位面积上的动量。在笛卡尔坐标系下 ，设$\displaystyle F=[F_x,F_y,F_z]$则
  $$\nabla \times F=\begin{vmatrix}\hat{i}&\hat{j}&\hat{k}\\ \partial_x&\partial_y&\partial_z\\ F_x&F_y&F_z\end{vmatrix}$$ 在一般坐标下面 $$\displaystyle (\nabla \times F)^k=\varepsilon^{klm}\nabla_l F_m$$ 其中，$\displaystyle \varepsilon$是Levi-Civita tensor,$\nabla$是协变导数。
• 常见计算关系 $$\nabla \times (v \times F)=(\nabla \cdot F+f \cdot F)v-(\nabla \cdot v +v\cdot \nabla)F \tag{K1}$$ $$\nabla \times (\nabla \times F)=\nabla(\nabla \cdot F)-\nabla^2 F\tag{K2}$$ $$\nabla \times(\nabla \varphi)=0 \tag{K3}$$ $$\nabla \times (\varphi F)=\nabla \varphi \times F+ \varphi \nabla \times F\tag{K4}$$

辐射条件定义

• 
  Helmholtz方程

  (∇2+k2)u=−f  in Rn(1)
  $$(\nabla^2+k^2)u=-f \ \ in\ \mathbb{R^n}\tag{1}$$ 其中$n=2,3$代表空间的维数;$f$是具有紧致支集的函数，代表有界能量源;$k>0$是常数，称为波数。

• $Sommerfeld$ 辐射条件
  

  lim|x|→∞|x|n−12(∂∂|x|−ik)u(x)=0(2)
  $$\lim_{|x|\rightarrow \infty} |x|^{\frac{n-1}{2}}(\frac{\partial}{\partial |x|}-ik)u(x)=0\tag{2}$$ 这里假设时间调和场是 $e^{-iwt}u$ ;辐射条件就是用来找到方程的唯一解。

• 辐射解radiating
  如果方程（1）的解$u$满足条件（2），就称$u$是辐射的（radiating）

简单例子

考虑

f=δ(x−x0)

$$f=\delta(x-x_0)$$ 此时，方程有无数解，比如像： $$u=cu_{+}+(1-c)u_{-}$$ 这里的 $$u(x)_{\pm}=\frac{e^{\pm ik|x-x_0|}}{4\pi|x-x_0|}$$ 练习：推导出上述解; 但只有$u_{+}$符合Sommerfeld 边界条件，他代表从点$x_0$发出的辐射场，而其他解都是非物理的。比如$u_{-}$可以解释成从无穷远处来的能量，最后Sinking 在$x_0$处。

关于赫姆霍兹方程

• 动机
  考虑波动方程

  （∇2−1c2∂2∂t2）u(r,t)=0
  $$（\nabla^2-\frac{1}{c^2}\frac{\partial^2}{\partial t^2}）u(r,t)=0$$ 用分离变量法$u(r,t)=A(r)T(t)$ 代入原方程可得：
  $$\frac{\nabla^2A(r)}{A(r)}=\frac{1}{T(t)c^2}\frac{\partial^2}{\partial t^2}T(t)$$ 观察方程两边，左边只依赖于$r$,右边只依赖于$t$,所以只有在左右两边都等于常数的才会成立
  $$\frac{\nabla^2A(r)}{A(r)}=\frac{1}{T(t)c^2}\frac{d^2}{d t^2}T(t)=-k^2$$ 所以会得到赫姆霍兹方程 $$(\Delta +k^2)A(r)=0\tag{1}$$ $$(\frac{d^2}{dt^2}+c^2k^2)T(t)=0\tag{2}$$
  注：$\displaystyle w:=ck$称之为角频率。

• 分离变量法解赫姆霍兹方程 考虑半径为$r=a$的圆盘上的情况， 引进极坐标$(r,\theta)$
  ,施加边界条件$\displaystyle A|_{r=a}=0$ 可得
  

  {Arr+1rAr+1r2Aθθ+k2AA(a,θ)=0=0
  $$\begin{cases}A_{rr}+\frac{1}{r}A_r+\frac{1}{r^2}A_{\theta\theta}+k^2A&=0\\ A(a,\theta)&=0 \end{cases}$$ 分离变量 $$A(r,\theta)=R(r)\Theta(\theta)$$
  得到 $$\begin{cases}\Theta'' +n^2\Theta&=0\\ r^2R''+rR'+r^2k^2R-n^2R&=0\end{cases}$$

• 近轴近似 将复的振幅写成$\displaystyle A(r)=u(r)e^{ikz}$, 定义$\Delta_{\bot}:=\frac{\partial^2}{\partial x^2}+\frac{\partial^2}{\partial y^2}$, 在满足条件 $$\bbox[18px,border:2px solid red]{\displaystyle \vert \frac{\partial^2 u}{\partial z^2} \vert \ll |k\frac{\partial u}{\partial z}|}\tag{轴不等式}$$ 时，$u$接近解以下方程的解 $$\bbox[18px,border:2px solid red]{\Delta_{\bot}u +2i k \frac{\partial u}{\partial z}=0}\tag{*}$$ 将$A(r)=u(r)e^{ikz}$带入Helmholtz 方程可得： $$(\frac{\partial^2}{\partial x^2}+\frac{\partial^2}{\partial y^2})ue^{ikz}+\frac{\partial^2}{\partial z^2}ue^{ikz}+2(\frac{\partial}{\partial z} u)ike^{ikz}=0$$ 由于轴不等式,可以将$\frac{\partial^2}{\partial z^2}ue^{ikz}$忽略不计，再将$u=Ae^{-ikz}$带回（*）可得关于$A(r)$的近轴方程 $$\Delta_{\bot}A +2i k \frac{\partial A}{\partial z}=0\tag{**}$$
• 非齐次赫姆霍兹方程 $$(\nabla^2+k^2)u=-f \ \ in\ \mathbb{R^n}$$ 为了使方程有唯一解，我们需要在无穷远处施加边界条件——$Sommerfeld$ 辐射条件： $$\lim_{r\rightarrow \infty} r^{\frac{n-1}{2}}(\frac{\partial}{\partial r}-ik)A(r\hat{x})=0$$ 这里的$|\hat{x}|=1$；在$Sommerfeld$边界条件下，我们可以用卷积的办法给出方程的解： $$\displaystyle A(x)=(G*f)(x)=\int_{\mathbb{R^n}} G(x-y)f(y) dy$$ 这里的$G$是方程的Green 函数。
  练习：考虑$f=\delta(x)$,推导出$G(x)=\frac{ie^{ik|x|}}{2k}$

Maxwell 方程的辐射条件 ^{Inverse Acoustic and ElectromagneticScattering Theory , Chapter6,p157}

• 本节考虑真空Maxwell方程 $$\displaystyle \begin{cases}curl\ E-ikH&=0\\ curl \ H +ikE&=0 \end{cases}\tag{1}$$ 这里$k$是波数，$\displaystyle k^2=(\varepsilon +\frac{i \sigma}{w})\mu w^2$
• 定理6.2 Stratton-Chu公式 设$D$ 是$C^2$的有界区域，$v$是$\partial D$的单位外法向量，设$E,H\in C^1(D)\cap C(\bar{D})$是Maxwell 方程（1）在区域$D$里面的解，则它们满足$Stratton-Chu$公式 $$\displaystyle E(x)=-curl \int_{\partial D}v(y)\times E(y)\Phi(x,y)ds(y) +\frac{1}{ik}curl\ curl \int_{\partial D}v(y)\times H(y)\Phi(x,y)ds(y)$$ $$\displaystyle H(x)=-curl \int_{\partial D}v(y)\times H(y)\Phi(x,y)ds(y) -\frac{1}{ik}curl\ curl \int_{\partial D}v(y)\times E(y)\Phi(x,y)ds(y)$$
• 定理6.4 设$E,H$是（1）的解，则$E,H$是无源场（$div E=0,div H=0$,比如说磁场，进多少出多少），且满足Helmholtz 方程 $$\Delta E+k^2E=0$$ $$\Delta H+k^2H=0$$ 反过来，设$E(or\ H)$是Helmholtz 方程的解，且$divE=0(or\ div H=0)$,那么$\displaystyle E,H:=\frac{curlE}{ik}(or\ H,E:=\frac{-curlH}{ik})$ 满足Maxwell 方程。
• Silver-Muller辐射条件 $$\lim_{r\to \infty} (H \times x-rE)=0\tag{2}$$ 或者 $$\lim_{r\to \infty} (E \times x-rH)=0$$ 这里的$\displaystyle r=|x|$
• 辐射解 满足条件（2）的Maxwell方程解被称为辐射。
• 定理6.6 Stratton-Chu公式 设有界区域$D$是某个无界$C^2$区域的开补集，$v$是$\partial D$单位外法向量。$E,H\in C^1(\mathbb{R^3}\setminus\bar{D})\cap C(\mathbb{R^3}\setminus D)$ 是Maxwell 方程在$\mathbb{R^3}\setminus \bar{D}$ 上的辐射解，则有Stratton-Chu公式： $$\displaystyle E(x)=curl \int_{\partial D}v(y)\times E(y)\Phi(x,y)ds(y) -\frac{1}{ik}curl\ curl \int_{\partial D}v(y)\times H(y)\Phi(x,y)ds(y)$$ $$\displaystyle H(x)=curl \int_{\partial D}v(y)\times H(y)\Phi(x,y)ds(y) +\frac{1}{ik}curl\ curl \int_{\partial D}v(y)\times E(y)\Phi(x,y)ds(y)$$