레이블이 debye shielding인 게시물을 표시합니다. 모든 게시물 표시
레이블이 debye shielding인 게시물을 표시합니다. 모든 게시물 표시

2018년 4월 15일 일요일

[플라즈마 물리][Plasma Physics]CH1 Introduction - Debye Shielding 디바이 차폐

"배워서 남주자"라는 가치를 가진 산동일크무크를 응원해주세요!

[Youtube Link]

Debye Shielding

Consider a test charge \(e\) and electron cloud around it in a singly charged plasma. Assume that

\(\bullet\) The ions are fixed; \(\frac{m_i}{m_e} \rightarrow \infty\)

\[\begin{aligned} \label{eq_6} n_i(r)=n \end{aligned}\]

\(\bullet\) The electrons obey Boltzman relation. In the presence of a potential energy \(q\phi\), the electron distribution function is

\[f(u)=A e^{-\frac{\frac{1}{2}mu^2+q\phi}{KT_e}}\]

Integrate this over \(du\)

\[\begin{aligned} \label{eq_7} n_e(r)=ne^{\frac{e\phi(r)}{K_BT}} \end{aligned}\]

Poisson’s equation in one dimension is

\[\begin{aligned} \epsilon_0 \nabla^2 \phi = \epsilon_0 \frac{d^2 \phi}{dx^2}=-e(n_i - n_e) \end{aligned}\]

substitution in to the Poisson’s equation,

\[\begin{aligned} \nabla^2 \phi(r) = \frac{e}{\epsilon_0}n(e^{\frac{e\phi(r)}{KT}}-1)\end{aligned}\]

assume that \(|\frac{e\phi}{KT}|<<1\), then by using taylor expansion; \(e^{\frac{e\phi(r)}{KT}}=1+\frac{e\phi(r)}{KT} +...\)

\[\begin{aligned} \nabla^2 \phi(r) &= \frac{e^2n}{\epsilon_0 KT}\phi(r)\\ &=\frac{\phi(r)}{\lambda_D^2}\end{aligned}\]

here we define Debye length \(\lambda_D\) as

\[\lambda_D = \sqrt{\frac{\epsilon_0 K T }{ne^2}}\]

In spherical coordinate

\[\begin{aligned} \frac{1}{r^2}\frac{d}{dr}\left( r^2 \frac{d\phi}{dr}\right) -\frac{\phi}{\lambda_D^2}&=0\\ \phi^{\prime \prime} + \frac{2}{r}\phi^{\prime} - \frac{1}{\lambda_D^2}\phi&=0\\ r\phi^{\prime \prime} + 2\phi^{\prime} - \frac{1}{\lambda_D^2}r\phi&=0\end{aligned}\]

Let \(\psi(r)=r\phi\), then \(\psi^{\prime}=\phi+r\phi^{\prime}\) and \(\psi^{\prime \prime}=2\phi^{\prime}+r\phi^{\prime \prime}\) so that

\[\psi^{\prime \prime}-\frac{1}{\lambda_D^2}\psi =0\]

\[\begin{aligned} \psi(r) &= C_1 e^{-\frac{r}{\lambda_D}}+C_2 e^{\frac{r}{\lambda_D}}\\ \phi(r) &= \frac{C_1}{r} e^{-\frac{r}{\lambda_D}}+\frac{C_2}{r} e^{\frac{r}{\lambda_D}}\end{aligned}\]

Applying the boundary conditions,

\[\begin{aligned} \phi &\rightarrow 0 \quad \quad \quad \quad ,\quad \quad r \rightarrow \infty \\ \phi &\rightarrow \frac{e}{4\pi\epsilon_0 r} \quad \quad , \quad \quad r \rightarrow 0\end{aligned}\]

we can find the constants.

\[C_1 = \frac{e}{4\pi \epsilon_0} \quad \quad \quad C_2 = 0\]

Hence, the solution is given by

\[\begin{aligned} \label{eq_debyelength} \boxed{\phi(r) = \frac{e}{4\pi\epsilon_0 r}e^{-\frac{r}{\lambda_D}}}\end{aligned}\]

The quantity \(\lambda_D\), called the Debye length, is a measure of the shielding distance or thickness of the sheath over which the influence of an individual charged particle is dominant.

\(\bullet\) \(\lambda_D\) = how long the shielding is effective.

\(\bullet\) Notice that electron temperature is used to define Debye length because it is more mobile; most of time this is true.

\(\bullet\) Useful forms of Eq.([eq_debyelength]) are

  • \(\lambda_D = 69\sqrt{\frac{T_e}{n}}[m]\) \(T_e\) in \(^{\circ}K\)

  • \(\lambda_D = 7430\sqrt{\frac{KT_e}{n}}[m]\) \(KT_e\) in \(eV\)

\(\bullet\) Trend

  • effective shielding \(\quad n \uparrow\) \(\lambda_D \downarrow\)

  • poor shielding \(\quad T \uparrow\) \(\lambda_D \uparrow\)

4차 산업혁명에 걸맞는 인터넷 기반 고등교육기관

이메일: ilkmooc@ilkmooc.kr

홈페이지주소: http://ilkmooc.kr

산동일크무크란? https://goo.gl/FnvqXd

2018년 4월 11일 수요일

[플라즈마 물리][Plasma Physics]Debye Shielding 디바이 차폐

Example Problem: Debye Shielding
Consider a positive point charge immersed in a plasma as we discussed in the class. Show that the net charge in the Debye shielding cloud exactly cancels the test charge. Assume that the ions are fixed and that \(e\phi << KT_e\). Note that we assumed that ion distribution is similar to that of the electrons in the class.

Answer:

While deriving, we assume potential to be

\[\begin{aligned} \frac{e\phi}{KT_e} << 1 \quad \quad \quad n_i \simeq n_e\\ \lambda_D = \sqrt{\frac{\epsilon KT}{ne^2}}\\ \phi(r) = \frac{e}{4\pi\epsilon_0 r }e^{-\frac{r}{\lambda_D}}\end{aligned}\]

Assume \(\frac{m_i}{m_e} \rightarrow \infty\) so that ions are fixed; \(n_i(r)=n\). Further assume electrons obey the Boltzmann distribution.

\[\begin{aligned} f(u) &= A e^{\frac{\frac{-1}{2}mv^2 + q\phi}{KT_e}}\\ n_e(r) &= ne^{\frac{e\phi(r)}{K_B T}}\end{aligned}\]

\[\begin{aligned} \rho &= -e(n_i - n_e) \\ &=-e \left(n - ne^{\frac{e\phi(r)}{K_BT}} \right)\\ &= ne \left(e^{\frac{e\phi(r)}{K_BT}} - 1 \right)\\\end{aligned}\]

for \(\frac{e\phi}{KT_e} << 1\),

\[\begin{aligned} e^{\frac{e\phi(r)}{K_B T}} = 1 + \frac{1}{2}\frac{e\phi(r)}{K_B T} + ...\end{aligned}\]

\[\begin{aligned} \rho &= ne(1+\frac{e\phi(r)}{K_B T} - 1) \\ &=\frac{ne^2 \phi(r)}{K_B T} \\ &= \frac{\epsilon_0}{\lambda_D^2}\phi(r)\end{aligned}\]

\[\begin{aligned} \int \rho dv &= Q \\ &= \int_V \frac{\epsilon_0}{\lambda_D^2}\frac{e}{4\pi\epsilon_0 r }e^{-\frac{r}{\lambda_D}} r^2 \sin\theta dr d\theta d\phi\\ &= \int \frac{4\pi \epsilon_0}{\lambda_D^2}\frac{e}{4\pi \epsilon_0} r e^{-\frac{r}{\lambda_D}}dr \\ &= \frac{e}{\lambda_D^2} \int_{0}^{\infty} r e^{-\frac{r}{\lambda_D}}dr \\ &= \frac{e}{\lambda_D^2} (\lambda_D^2 e^{-\frac{r}{\lambda_D}}-\lambda_D r e^{-\frac{r}{\lambda_D}}) =0\end{aligned}\]

On the last step, when \(r = \lambda_D\), charge becomes zero. Hence, the Debye sheilding cloud exactly cancels the test charge.

4차 산업혁명에 걸맞는 인터넷 기반 고등교육기관

이메일: ilkmooc@ilkmooc.kr

홈페이지주소: http://ilkmooc.kr

산동일크무크란? https://goo.gl/FnvqXd

2018년 4월 10일 화요일

[플라즈마 물리]F.F. Chen Problem 1.5 Solution Debye Shielding

Problem In a strictly steady state situation both the ions and the electrons will follow the Boltzmann relation

\[n_i = n_0 e^{\frac{-q_i \phi }{KT_i}}\]

For the case of an infinite, transparent grid charged to a potential \(\phi\), show that the shielding distance is then given approximately by

\[\lambda_D^{-2} = \frac{ne^2}{\epsilon_0} \left( \frac{1}{KT_e} + \frac{1}{KT_i} \right)\]

Answer

\[\begin{aligned} \nabla \cdot \vec{E} &=\frac{\rho}{\epsilon_0} = \frac{e}{\epsilon_0}(n_i - n_e) \nonumber \\ &= - \frac{e}{\epsilon_0} n_0 \left( e^{-\frac{e\phi}{KT_i}} - e^{\frac{e\phi}{KT_e}} \right) \nonumber \\ &= - \frac{e}{\epsilon_0} n_0 \left( 1 -\frac{e\phi}{KT_i} - 1 - \frac{e\phi}{KT_e} \right) \nonumber \\ \frac{d^2 \phi }{dx^2 }&=\frac{e^2 n_0}{\epsilon_0}\left( \frac{1}{KT_i} + \frac{1}{KT_e} \right) \phi \nonumber\end{aligned}\]

Hence,

\[\begin{aligned} \frac{1}{\lambda_D^2} = \frac{e^2 n_0}{\epsilon_0}\left( \frac{1}{KT_i} + \frac{1}{KT_e} \right) \nonumber\end{aligned}\]

\[\lambda_D = \sqrt{\left( \frac{\epsilon_0 K T_e T_i}{ne^2 (T_i + T_e)} \right)}\]

  • \(T_i << T_e \quad \quad \lambda_D \simeq \sqrt{\left( \frac{\epsilon_0 K T_i}{ne^2} \right)}\)

  • \(T_e << T_i \quad \quad \lambda_D \simeq \sqrt{\left( \frac{\epsilon_0 K T_e}{ne^2} \right)}\)

4차 산업혁명에 걸맞는 인터넷 기반 고등교육기관

이메일: ilkmooc@ilkmooc.kr

홈페이지주소: http://ilkmooc.kr

산동일크무크란? https://goo.gl/FnvqXd