|Message: Trajectory in dipole field independent of energy !||Not Logged In (login)|
Click on the Forum title, e.g. on the "Forums by Category" page, to read a sequence of postings to the Forum and its threads all in one page. If you are only interested in one thread or the thread following a specific posting, click the thread or the posting, which takes you to a smaller page, which contains only the part you are interested in and may be easier to navigate.
Messages are "chained" if there are only replies at the first level, i.e. 1/1.html, 1/1/1.html etc. In case of "chained" messages the message number is replaced by the icon and there is no indentation.
Inline: Display the subject line only or also the text of the posting(s); for the choice "All" the "Outline" choices are switched off.
|1||0||1||no text / full text of posting|
|2||1||All||text for level 1 only / text for All postings|
Outline: Choose the depth of the posting thread, successive toggle controls provide increasing detail.
|1||2||1||2 levels / 1 level (original posting)|
|2||3||2||3 levels / 2 levels|
|3||3||All||3 levels / all levels (all postings)|
So I have my Earth with its dipole working, and running pretty speedily. Now comes the big problem: the paths of charged particles in the field are independent of energy! I have gone through the following steps, changing no code except the GetFieldValue() function... First, uniform B-field along z axis. As expected, the radius of curvature of the particle's path increases with energy. As a side note, I found that above a certain energy, about 8.5MeV for an electron in a 3 nT field, the particle would go absolutely straight, zero curvature. Then, a B-field along the z-axis where field strength is a function of the particle's z coordinate. Again, radius of curvature increases with energy. Now, a more complicated function, looking like Bfield = y Bfield = z Bfield = x And still the path depends on energy. Adding another layer of complexity: Bfield = y*z Bfield = z*x Bfield = x*y And as always, more energy makes the particle take larger loops with less curvature. So then I implemented my dipole field: Bfield = 3*coeff*z*x Bfield = 3*coeff*z*y Bfield = coeff*(3*z*z - r*r) and suddenly, the trajectory becomes independent of energy. For the same starting position and direction, an electron ranging anywhere from 10 eV to 10 TeV will take exactly the same path. The initial direction has a small effect; it determines whether the particle travels forward or backward along the field line trapping it but otherwise does not matter. For example, these two run macros: /gps/position -8000. -8000. -8000. km /gps/energy 100 eV /gps/direction 1. 0. 0. and /gps/position -8000. -8000. -8000. km /gps/energy 10 MeV /gps/direction 1. 0. 1. produce identical output. What is the cause of this problem?
|Inline Depth:||Outline Depth:||Add message:|