 |
 |
|
| |
|
 |
|
| |
Polywell in Matlab
To get some idea of what electron look like traveling around in the cusp magnetic bottle of a polywell reactor I ploted the magnetic fields
from 6 square "coils" onto a 100x100x100 grid using a linear interpolation for field values in between grid points. For the integration of the
electron path I used a modified runge-kutta algorithm using a solution to the equations of motion in a uniform EM field, instead of the Euler approximation
used in a normal runge-kutta. The b-field is approximatly 0.7 T on face axis, which are 2m per side and 3m from opposing face.
No electric field or inter-electron effects are considered.
These first plots show 50 runs of electrons starting from the center at random directions and 10kev. The path is plotted for up to 10e-6 seconds or until the path escapes
the bounds of the box.
50 Random Runs.
(larger) |
Runs rotated about axial symetries.
(larger) |
Next are 15 runs shot almost directly at the line cusp between to faces, with slight deflection toward each face. Small deflections toward a face result in a
deflection toward a corner cusp.
Finally is a plot of projected electron density using a compilation of the 50 runs where equal time datapoints are slotted into
a grid. The plots are a slice in the X-Y directions. First the slice through Z=0m, and then slice at about Z=0.8m; Noting that the faces are at +- 1.5m.
The density is scaled so that the highest density value is 1.
|
|
| |
|
| |
|
|
|
|
| |
© 2001-2008 Andromeda Space. All Rights Reserved. |
|
