Relativistic particle in non-uniform magnetic field (math)

In summary: F=mv^2/r^2).In summary, the code simulates a 3-D magnetic field based on a combination of wires and wire loops. The code produces a 3D vector field of 11x11x11. The code is attempting to send a virtual particle into the field so that the behavior can be analyzed. There are some problems with the equations even in the non-relativistic limit and the equation for determining the position in the field is nonsense.
  • #1
msailer1234
13
1
Hi everyone,
I am currently working on a project in Physics analyzing the possibility of magnetic fields protecting against radiation particles such as fully ionized nuclei or electrons. I wrote a code in Matlab that simulates a 3-D magnetic field based on a combination of wires and wire loops. The code produces a 3D vector field of 11x11x11. So I have Bx,By, and Bz for every point in the 11x11x11 box. I am now attempting to create a code which sends a virtual particle into the magnetic field so I can see what happens.

What I am doing is treating each point as a box with a uniform magnetic field which the particle travels through. I want to find an equation that can take the point and velocity the particle enters into a box and calculates the point and velocity it leaves the box. I am tripping over the calculations. I have attached my hand written calculations (I hope they are clearly written enough), but I don't know if I am incorporating relativity correctly. I also don't think I am approaching it correctly because the equation gives me Radius_x, Radius_y, and Radius_z which doesn't make sense to me. If anyone has any corrections or ideas, please let me know! Thanks!
 

Attachments

  • img20180117_20264237.jpg
    img20180117_20264237.jpg
    27.6 KB · Views: 577
Physics news on Phys.org
  • #2
Do you already have code for a virtual non-relativistic particle entering your box? If not, I would think that you need to write that first to gain the experience. Then you can attempt the more complicated relativistic case and check it by making sure that it reduces to the non-relativistic results for small speeds. I should also point out that you got the components of v wrong, if that's what you have under "given". You write the x-component as ##\frac{v_x}{\sqrt{1-v^2/c^2}}.## It should be just ##v_x.## There are other problems with your equations even in the non-relativistic limit. The equation ##r=mv^2/(q\vec{v}\times \vec B )## is nonsense.
 
  • Like
Likes msailer1234
  • #3
kuruman said:
Do you already have code for a virtual non-relativistic particle entering your box? If not, I would think that you need to write that first to gain the experience. Then you can attempt the more complicated relativistic case and check it by making sure that it reduces to the non-relativistic results for small speeds. I should also point out that you got the components of v wrong, if that's what you have under "given". You write the x-component as ##\frac{v_x}{\sqrt{1-v^2/c^2}}.## It should be just ##v_x.## There are other problems with your equations even in the non-relativistic limit. The equation ##r=mv^2/(q\vec{v}\times \vec B )## is nonsense.
Thank you for responding! Okay, so that looks like one of my problems. I got that equation by combining the force equations F=mv^2/r and F=qvxB. So am I correct in saying that that reasoning only works if used with the velocity being perpendicular to the B field making the equation r = mv/qB? Also, if this is a relativistic velocity, shouldn't it be relativistic in each spatial dimension?
 
  • #4
To begin with the symbol for v in relativistic equation stands for the relative velocity between two frames, O and O'. If Observer O sees observer O' move with velocity ##\vec{v}##, observer O' will see observer O move with velocity ##-\vec{v}##. That's all there is to that. You need the solve the following differential equation in 3D:
$$\frac{d\vec{v}}{dt}=\vec{\omega}\times \vec{v}$$ where $$\vec{\omega}_B=\frac{q\vec{B}}{\gamma m c};~~~\left( \gamma =\frac{1}{\sqrt{1-v^2/c^2}} \right)$$
Note that ##\omega_B## is velocity-dependent. Do you see why it is easier to do the non-relativistic calculation with ##\gamma = 1##?
 
  • #5
A radius with x,y,z coordinates is a perfectly valid thing. In addition to the curvature radius it also gives you the direction of the curvature.
You'll have to find a way to determine where your curved track (or circle) intersects the boundary of the cubes unless the deflection is very small.
 
  • Like
Likes msailer1234
  • #6
Since you are working with Cartesian coordinates, why are you trying to move to "r" (spherical?) ... it's only going to make things more complicated... I think you have 11x11x11 boxes (aka regions with a given Bx,By,Bz)...
if you assume that the particle in each box (i) it transcends gets a kick to its velocity in the x,y,z-direction (given by the relevant force), I think it's pretty simple.
for example if in box 0 it has a velocity [itex](v_x, v_y, v_z)[/itex], and use that to determine the box1 (next step), in box1 it would have [itex](v_x+ F_x/m , v_y + F_y/m , v_z+ F_z/m)[/itex]... ok that was for classical mechanics... for relativistic, the acceleration kick (F/m's) should be changed appropriately (see eg. https://en.wikipedia.org/wiki/Relativistic_mechanics#Force)
 
  • Like
Likes msailer1234
  • #7
Do you have to work with boxes? Following the track of particles could be easier.
Knowing some key parameters of the system would help to see which approximations work how well.
 
  • Like
Likes msailer1234
  • #8
mfb said:
Do you have to work with boxes? Following the track of particles could be easier.
Well by saying we have cubic spaces of known magnetic field, I think that's the best way to deal with the problem ?
The track of the particles is supposed to be emulated based on how the particle interacts with the field [at each point the velocity should be updated to get the next point the particle will be in].
Or maybe I haven't understood the problem.
 
  • Like
Likes msailer1234
  • #9
The approach so far looks more like OP wants to model the trajectory inside the cubes as circle and then let that intersect the cube boundaries to do the same again in the next cube.
 
  • Like
Likes msailer1234
  • #10
mfb said:
The approach so far looks more like OP wants to model the trajectory inside the cubes as circle and then let that intersect the cube boundaries to do the same again in the next cube.
From what my code gives me, I do think having uniform magnetic fields in several boxes is the best approach for me to model this. Thank you everyone for the replies. I am still looking into it, and I am consulting some professors at my college.
 
  • #11
Hi everyone,
I thought I would update you on how the project went. I apologize for taking so long to respond, but the box method did work great for an approximation. In hindsight, calculating the new magnetic field at each time step based on the particle's position would have been more efficient and practical, but it was a great learning experience. Thank you for your help!
 
  • Like
Likes mfb

1. What is a relativistic particle in a non-uniform magnetic field?

A relativistic particle in a non-uniform magnetic field is a theoretical concept in physics that describes the motion of a charged particle in a magnetic field that varies in strength and direction. This type of motion is characterized by the principles of relativity and can be described using mathematical equations.

2. How does a non-uniform magnetic field affect the trajectory of a relativistic particle?

A non-uniform magnetic field can cause a relativistic particle to experience a curved trajectory. This is because the particle will experience a varying magnetic force as it moves through the field, causing it to accelerate in different directions at different points in its path.

3. What is the role of the Lorentz factor in describing a relativistic particle in a non-uniform magnetic field?

The Lorentz factor is a mathematical term that is used to describe the relativistic effects on a particle's mass, length, and time as it moves at high speeds. In the case of a relativistic particle in a non-uniform magnetic field, the Lorentz factor is used to account for the particle's changing velocity and mass as it moves through the field.

4. How can the motion of a relativistic particle in a non-uniform magnetic field be calculated?

The motion of a relativistic particle in a non-uniform magnetic field can be calculated using a set of mathematical equations, such as the Lorentz force law and the equations of motion. These equations take into account the particle's charge, mass, velocity, and the strength and direction of the magnetic field to determine its trajectory.

5. What are some real-world applications of studying relativistic particles in non-uniform magnetic fields?

The study of relativistic particles in non-uniform magnetic fields has many practical applications, such as in the design and operation of particle accelerators and other high-energy physics experiments. It also helps us to better understand the behavior of charged particles in natural phenomena, such as cosmic rays and solar flares.

Similar threads

  • High Energy, Nuclear, Particle Physics
Replies
2
Views
913
  • High Energy, Nuclear, Particle Physics
Replies
8
Views
1K
  • Special and General Relativity
Replies
20
Views
919
  • Special and General Relativity
Replies
25
Views
1K
  • Introductory Physics Homework Help
Replies
16
Views
320
  • Introductory Physics Homework Help
Replies
14
Views
1K
  • Quantum Physics
Replies
10
Views
2K
  • Introductory Physics Homework Help
Replies
7
Views
768
  • Advanced Physics Homework Help
Replies
1
Views
842
Replies
26
Views
1K
Back
Top