跳到主要内容

東北大学 理学研究科 物理学専攻 2018年8月実施 問題3(電磁気学)

Author

Miyake

Description

原題に基づく要約(日本語)

公式原題

[1] 帯電円環

半径 aa の細い絶縁円環を原点中心の xyxy 平面に置き、正の線電荷密度 ρ\rho で一様に帯電させる。真空の定数は ε0,μ0\varepsilon_0,\mu_0、観測点は P=(0,0,z)P=(0,0,z)

  1. 長さ dldl の微小部分が PP に作る電場の大きさを求める。
  2. 円の中心に対して反対側の同じ長さの部分と組み合わせ、合成電場の大きさと z>0,z<0z>0,z<0 での向きを求める。
  3. 円環全体が作る電場を求める。
  4. a=ρ=ε0=1a=\rho=\varepsilon_0=1 として Ez(z)E_z(z) の概形を描く。
  5. 円環を正の角速度 ω\omegazz 軸のまわりに回転させる。回転方向は右ねじが +z+z に進む方向とする。PP の磁束密度を求める。

[2] 空洞共振器

完全導体で囲まれた 0<x<a, 0<y<a, 0<z<b0<x<a,\ 0<y<a,\ 0<z<b の空洞を考える。a>b>0a>b>0 とし、内部は無源の真空で Maxwell 方程式を満たす。Ex=Ey=0, Ez=X(x)Y(y)Z(z)sinωtE_x=E_y=0,\ E_z=X(x)Y(y)Z(z)\sin\omega t とする。

  1. E=0\nabla\cdot\boldsymbol E=0 から、非零モードの ZZ が定数となることを示す。
  2. 波動方程式 2E=ε0μ0t2E\nabla^2\boldsymbol E=\varepsilon_0\mu_0\partial_t^2\boldsymbol E から X/X+Y/Y=ε0μ0ω2X''/X+Y''/Y=-\varepsilon_0\mu_0\omega^2 を導く。
  3. 左辺の各項を kx2,ky2-k_x^2,-k_y^2 と置き、X,YX,Y の一般解を求める。
  4. 導体表面の接線方向電場が0である条件を用い、EzE_z の正弦関数による表示と許される kx,kyk_x,k_y を求める。
  5. このモードの最小共振角振動数 ω0\omega_0 を求める。
  6. ω=ω0\omega=\omega_0 における磁束密度 B\boldsymbol B を求める。

円環と直方体空洞の幾何条件。

题目描述

【带电圆环】

半径为 aa 的细绝缘圆环位于 xyxy 平面,圆心为原点,均匀正线电荷密度为 ρ\rho,观察点为 P=(0,0,z)P=(0,0,z),真空常数为 ε0,μ0\varepsilon_0,\mu_0

  1. 求长度 dldl 的圆环微元在 PP 产生的电场大小。
  2. 将它与圆心对侧的等长微元配对,求合电场,并分别说明 z>0,z<0z>0,z<0 时的方向。
  3. 求整个圆环的电场。
  4. a=ρ=ε0=1a=\rho=\varepsilon_0=1,画出 Ez(z)E_z(z) 的概形。
  5. 令圆环绕 zz 轴以正角速度 ω\omega 转动,方向满足右手定则指向 +z+z,求 PP 处的磁感应强度。

【空腔谐振器】

完全导体包围的真空空腔为 0<x<a, 0<y<a, 0<z<b0<x<a,\ 0<y<a,\ 0<z<b,其中 a>b>0a>b>0。内部没有电荷和电流,设 Ex=Ey=0E_x=E_y=0Ez=X(x)Y(y)Z(z)sinωtE_z=X(x)Y(y)Z(z)\sin\omega t

  1. 用电场无散条件证明非零模态的 ZZ 为常数。
  2. 由波动方程导出 X/X+Y/Y=ε0μ0ω2X''/X+Y''/Y=-\varepsilon_0\mu_0\omega^2
  3. 令两项分别为 kx2,ky2-k_x^2,-k_y^2,求 X,YX,Y 的通解。
  4. 用导体表面切向电场为零的条件,写出正弦形式的 EzE_z 及允许的 kx,kyk_x,k_y
  5. 求该类模态的最低共振角频率 ω0\omega_0
  6. ω=ω0\omega=\omega_0 时的磁感应强度 B\boldsymbol B

Kai

[1]

1)

微小区間にある電荷は ρdl\rho dl で、APの長さは a2+z2\sqrt{a^2+z^2} であるから、

求める電場の大きさは

14πε0ρdla2+z2\begin{aligned} \frac{1}{4 \pi \varepsilon_0} \frac{\rho dl}{a^2+z^2} \end{aligned}

である。

2)

z軸に垂直な成分は打ち消し合うことを考慮して、 求める電場の大きさは

14πε0ρdla2+z2za2+z22=12πε0zρdl(a2+z2)3/2\begin{aligned} \frac{1}{4 \pi \varepsilon_0} \frac{\rho dl}{a^2+z^2} \frac{|z|}{\sqrt{a^2+z^2}} \cdot 2 = \frac{1}{2 \pi \varepsilon_0} \frac{|z| \rho dl}{(a^2+z^2)^{3/2}} \end{aligned}

である。

また、電場の向きは、 z>0z \gt 0 のときは +z方向、 z<0z \lt 0 のときは -z方向である。

3)

z方向の単位ベクトルを z^\hat{z} として、

E=0πa12πε0zρdl(a2+z2)3/2z^=12πε0zρπa(a2+z2)3/2z^=12ε0azρ(a2+z2)3/2z^\begin{aligned} \vec{E} &= \int_0^{\pi a} \frac{1}{2 \pi \varepsilon_0} \frac{z \rho dl}{(a^2+z^2)^{3/2}} \hat{z} \\ &= \frac{1}{2 \pi \varepsilon_0} \frac{z \rho \cdot \pi a}{(a^2+z^2)^{3/2}} \hat{z} \\ &= \frac{1}{2 \varepsilon_0} \frac{az \rho}{(a^2+z^2)^{3/2}} \hat{z} \end{aligned}

である。

4)

a=ρ=ε0=1a=\rho=\varepsilon_0=1 とすると

Ez(z)=z2(1+z2)3/2,Ez(z)=12z22(1+z2)5/2.E_z(z)=\frac{z}{2(1+z^2)^{3/2}},\qquad E_z'(z)=\frac{1-2z^2}{2(1+z^2)^{5/2}}.

奇関数で、原点を通り、z±z\to\pm\infty00 に近づく。極大点は (1/2,1/(33))(1/\sqrt2,1/(3\sqrt3))、極小点は (1/2,1/(33))(-1/\sqrt2,-1/(3\sqrt3)) である。

円環の軸上電場の概形。

5)

z方向の単位ベクトルを z^\hat{z} として、

B=02πaμ04πaωρdla2+z2aa2+z2z^=02πaμ04πa2ωρdl(a2+z2)3/2z^=μ04πa2ωρ2πa(a2+z2)3/2z^=μ0a3ωρ2(a2+z2)3/2z^\begin{aligned} \vec{B} &= \int_0^{2 \pi a} \frac{\mu_0}{4 \pi} \frac{a \omega \rho dl}{a^2+z^2} \frac{a}{\sqrt{a^2+z^2}} \hat{z} \\ &= \int_0^{2 \pi a} \frac{\mu_0}{4 \pi} \frac{a^2 \omega \rho dl}{(a^2+z^2)^{3/2}} \hat{z} \\ &= \frac{\mu_0}{4 \pi} \frac{a^2 \omega \rho \cdot 2 \pi a}{(a^2+z^2)^{3/2}} \hat{z} \\ &= \frac{\mu_0 a^3 \omega \rho}{2 (a^2+z^2)^{3/2}} \hat{z} \end{aligned}

である。

[2]

1)

E=0\nabla\cdot\mathbf E=0 より

0=E=Exx+Eyy+Ezz=Ezz=X(x)Y(y)dZ(z)dzsinωt\begin{aligned} 0 = \vec{\nabla} \cdot \vec{E} = \frac{\partial E_x}{\partial x} + \frac{\partial E_y}{\partial y} + \frac{\partial E_z}{\partial z} = \frac{\partial E_z}{\partial z} = X(x) Y(y) \frac{dZ(z)}{dz} \sin \omega t \end{aligned}

であるから、

dZ(z)dz=0\begin{aligned} \frac{dZ(z)}{dz} = 0 \end{aligned}

であり、 Z(z)Z(z)zz に依存しない定数である。

2)

波動方程式の zz 成分の両辺に Ez(x,y,t)=X(x)Y(y)ZsinωtE_z(x,y,t) = X(x) Y(y) Z \sin \omega t を代入すると、

d2X(x)dx2Y(y)Zsinωt+X(x)d2Y(y)dy2Zsinωt=ε0μ0ω2X(x)Y(y)Zsinωt\begin{aligned} \frac{d^2 X(x)}{dx^2} Y(y) Z \sin \omega t + X(x) \frac{d^2 Y(y)}{dy^2} Z \sin \omega t = - \varepsilon_0 \mu_0 \omega^2 X(x) Y(y) Z \sin \omega t \end{aligned}

であるから、

1Xd2Xdx2+1Yd2Ydy2=ε0μ0ω2\begin{aligned} \frac{1}{X} \frac{d^2 X}{dx^2} + \frac{1}{Y} \frac{d^2 Y}{dy^2} = - \varepsilon_0 \mu_0 \omega^2 \end{aligned}

を得る。

3)

1Xd2Xdx2=kx2,    1Yd2Ydy2=ky2\begin{aligned} \frac{1}{X} \frac{d^2 X}{dx^2} = -k_x^2 , \ \ \ \ \frac{1}{Y} \frac{d^2 Y}{dy^2} = -k_y^2 \end{aligned}

より、

d2Xdx2=kx2X,    d2Ydy2=ky2Y\begin{aligned} \frac{d^2 X}{dx^2} = -k_x^2 X , \ \ \ \ \frac{d^2 Y}{dy^2} = -k_y^2 Y \end{aligned}

であるから、 X,YX, Y の一般解は、次のようになる:

{X(x)=Asinkxx+BcoskxxY(y)=Csinkyy+Dcoskyy.\begin{aligned} \begin{cases} X(x) = A \sin k_x x + B \cos k_x x \\ Y(y) = C \sin k_y y + D \cos k_y y \end{cases} . \end{aligned}

ここで、 A,B,C,DA, B, C, D は積分定数である。

4)

(i) X(0)=0X(0)=0 より、 B=0B=0 であり、

X(x)=Asinkxx\begin{aligned} X(x) = A \sin k_x x \end{aligned}

を得る。

Y(0)=0Y(0)=0 より、 D=0D=0 であり、

Y(y)=Csinkyy\begin{aligned} Y(y) = C \sin k_y y \end{aligned}

を得る。

よって、

Ez=AsinkxxCsinkyyZsinωt=Ez0sinkxxsinkyysinωt\begin{aligned} E_z &= A \sin k_x x \cdot C \sin k_y y \cdot Z \sin \omega t \\ &= E_{z0} \sin k_x x \sin k_y y \sin \omega t \end{aligned}

を得る。 ここで、 Ez0=ACZE_{z0} = A C Z とした。

(ii) X(a)=0X(a)=0 より、

kx=πanx    (nx=1,2,)\begin{aligned} k_x = \frac{\pi}{a} n_x \ \ \ \ (n_x = 1, 2, \cdots) \end{aligned}

を得る。

Y(a)=0Y(a)=0 より、

ky=πany    (ny=1,2,)\begin{aligned} k_y = \frac{\pi}{a} n_y \ \ \ \ (n_y = 1, 2, \cdots) \end{aligned}

を得る。

5)

問2), 3)より、

kx2ky2=ε0μ0ω2  ω2=kx2+ky2ε0μ0\begin{aligned} -k_x^2 - k_y^2 = - \varepsilon_0 \mu_0 \omega^2 \end{aligned} \begin{aligned} \therefore \ \ \omega^2 = \frac{k_x^2 + k_y^2}{\varepsilon_0 \mu_0} \end{aligned}

である。

これに、問4)で得た kx=πnx/a,ky=πny/ak_x = \pi n_x / a, k_y = \pi n_y / a を代入すると、

ω2=π2(nx2+ny2)ε0μ0a2\begin{aligned} \omega^2 = \frac{\pi^2 (n_x^2 + n_y^2)}{\varepsilon_0 \mu_0 a^2} \end{aligned}

となる。

ω=ω0\omega = \omega_0 となるのは nx=ny=1n_x=n_y=1 のときなので、

ω02=2π2ε0μ0a2  ω0=2π2ε0μ0a2\begin{aligned} \omega_0^2 &= \frac{2 \pi^2}{\varepsilon_0 \mu_0 a^2} \\ \therefore \ \ \omega_0 &= \sqrt{\frac{2 \pi^2}{\varepsilon_0 \mu_0 a^2}} \end{aligned}

を得る。

6)

k=π/ak=\pi/a とおけば、基本モードの電場は

E=Ez0sin(kx)sin(ky)sin(ω0t)z^,ω02=2k2ε0μ0.\mathbf E=E_{z0}\sin(kx)\sin(ky)\sin(\omega_0t)\,\hat{\mathbf z},\qquad \omega_0^2=\frac{2k^2}{\varepsilon_0\mu_0}.

Faraday の法則 tB=×E\partial_t\mathbf B=-\nabla\times\mathbf E を時間積分すると、共振に伴う振動磁場は

B=Ez0kω0[sin(kx)cos(ky)x^cos(kx)sin(ky)y^]cos(ω0t).\boxed{\mathbf B=\frac{E_{z0}k}{\omega_0}\left[\sin(kx)\cos(ky)\,\hat{\mathbf x}-\cos(kx)\sin(ky)\,\hat{\mathbf y}\right]\cos(\omega_0t)}.

Bz=0B_z=0 であり、直接微分すると B=0\nabla\cdot\mathbf B=0×B=ε0μ0tE\nabla\times\mathbf B=\varepsilon_0\mu_0\partial_t\mathbf E も満たす。共振とは無関係な静磁場の加算はここでは含めていない。

Reference