Difference between revisions of "Solution details"

From New IAC Wiki
Jump to navigation Jump to search
Line 45: Line 45:
  
  
The modified Bessel functions, first and second, are the solution for the previous equation but hte boundary conditions determines which one to use, in this case  r'-> 0, <math> n-> \inifty </math>
+
The modified Bessel functions, first and second, are the solution for the previous equation but hte boundary conditions determines which one to use, in this case  r'-> 0, <math> n-> \infty </math>

Revision as of 22:54, 25 October 2013

asymptotic solution details for Boltzmann equation for a hole has a uniform electric field

[math] (\frac {\partial^2{}}{\partial{x^2}} +\frac {\partial^2{}}{\partial{x^2}})[/math]n + [math] D_L \frac {\partial^2{}}{\partial{z^2}}[/math] - [math] W \frac {\partial{}}{\partial{z}}[/math] n = 0

Steps to solve Boltzmann equation

for the previous equation let consider the asymptotic solution has the form:

[math] n(x', y', z') = e^{\lambda_L z'} V(x,y,z) [/math]

so

[math] \nabla'^2 V = \lambda_L^2 V [/math]

where

[math] \nabla'^2 V = \frac {\partial^2{}}{\partial{x'^2}} + \frac {\partial^2{}}{\partial{y'^2}} + \frac {\partial^2{}}{\partial{z^2}}[/math]

and

[math] x' = \frac {D_L}{D} x [/math] [math] y' = \frac {D_L}{D} y [/math]

In spherical coordinates:

[math] \frac {1}{r'^2} \frac{\partial{}}{\partial{r'}}r'^2 \frac{\partial{V}}{\partial{r'}} + \frac {1}{r'^2 sin\theta'} \frac{\partial{}}{\partial{\theta}} sin\theta \frac{\partial{V}}{\partial{\theta}} = \lambda_L^2 V [/math] which is symmetric in [math]\phi[/math] direction.

Assuming [math]V(r',\theta) = R_k(r')P_k(\mu) [/math]the solution of the zenith angle direction is the Legendre polynomial, and can be written as:

[math]\frac {1}{r'sin\theta} \frac{\partial{}}{\partial{\theta}} sin\theta\frac{\partial{V}}{\partial{\theta}} = R_k(r') \frac{d}{d \mu} \left [ (1- \mu^2) \frac{d{P_k(\mu)}}{d{\mu}} \right] [/math]

and


[math] \frac{d}{d \mu} \left [ (1- \mu^2) \frac{d{P_k(\mu)}}{d{\mu}} \right]= -k(k+1) P_k(\mu) [/math]

so,

[math] \frac {1}{r'^2} \frac{d}{dr'}\left (r'^2 \frac{dR_k}{dr'}\right) - \left [ \frac{k(k+1)}{r'^2} +\lambda_L^2 \right]R_k = \frac{d^2 R_k}{dr'^2} +\frac{2}{r'} \frac{dR_k}{dr'}-\left [ \frac{k(k+1)}{r'^2} +\lambda_L^2 \right]R_k = 0 [/math]


The modified Bessel functions, first and second, are the solution for the previous equation but hte boundary conditions determines which one to use, in this case r'-> 0, [math] n-\gt \infty [/math]