I have some questions relating movement of a particle calculating with mathematica.
A particle (mass: m, charge: q) moves on x-axis from left with speed $v_0$ to a magnetic field $\vec{B} = (0,0,B)$. Entry is at $\vec{r}(t=0)=(0,0,0)$.
- How would equation of movement look like for this problem ?
I tried this:
In[23]:= r[t_] := {x[t], y[t], z[t]}
In[24]:= electricField = 0, Efield, 0;
In[24]:= magneticField = 0 {0, 0, Bfield};
In[26]:= force = q (electricField + Cross[r'[t], magneticField]);
In[27]:= DSolve[{m r''[t] == force, r[0] == {0, 0, 0},
r[0] == {v1, v2, v3}}, r[t], t] //. {Bfield -> m \[Omega] / q,
Efield -> V Bfield} // ExpandAll // Simplify
During evaluation of In[27]:= DSolve::bvnul: For some branches of the general solution, the given boundary conditions lead to an empty solution.
Out[27]= {}
In[16]:= Clear[eqns]
In[29]:= eqns = Map[Thread, {m r''[t] == force, r[0] == {0, 0, 0},
r'[0] == {v1, v2, v3}}] //. {Bfield -> m \[Omega] /q,
Efield -> V Bfield} // ExpandAll // Simplify
Out[29]= {{electricField q == m (x^\[Prime]\[Prime])[t],
electricField q == m (y^\[Prime]\[Prime])[t],
electricField q == m (z^\[Prime]\[Prime])[t]}, {x[0] == 0,
y[0] == 0, z[0] == 0}, {v1 == Derivative[1][x][0],
v2 == Derivative[1][y][0], v3 == Derivative[1][z][0]}}
In[31]:= solution1 =
Dsolve[eqns, {x[t], y[t], z[t]}, t] // ExpandAll // Simplify
Out[31]= Dsolve[{{electricField q == m (x^\[Prime]\[Prime])[t],
electricField q == m (y^\[Prime]\[Prime])[t],
electricField q == m (z^\[Prime]\[Prime])[t]}, {x[0] == 0,
y[0] == 0, z[0] == 0}, {v1 == Derivative[1][x][0],
v2 == Derivative[1][y][0], v3 == Derivative[1][z][0]}}, {x[t],
y[t], z[t]}, t]
- When we start with $\vec{r}^{(0)}(t)=q(v_0t,0,0)$ as zeroth approximation and insert it in Lorentz-force, how would equation of movement look like ?
How can I write a routine that automatises this iteration till 6 iterations ? I read about "Nest" command but how to use it ?
- Charged particles describe circuits in homogenous magnetic field. How can I find exact solution with mathematica and see what angular frequency and radius they have ? There should be a possibility to develop exact solutions in power of $B$, how tovisualise it with mathematica ?