- #1
- 1,798
- 33
A colleague and I are looking at modelling a hot electron beam hitting a initially charge neutral plasma. Initially we're looking at the 1D problem, the equations we're using are:
[tex]
\begin{array}{rcl}
\gamma^{3}(v/c)\left(\frac{\partial v}{\partial t}+v\frac{\partial v}{\partial x}\right) & = & -\frac{e}{m}E \\
\frac{\partial n}{\partial t}+\frac{\partial }{\partial x}(nv) & = & 0 \\
\frac{\partial E}{\partial t}+\frac{E}{\eta\epsilon_{0}} & = & \frac{nev}{\epsilon_{0}}
\end{array}
[/tex]
where e is the charge on the electron, m is the mass of the electron, [itex]\eta[/itex] is the resistivity of the charge neutral plasma, v is the speed of the electrons, n is the number density of the electrons and E is the electric field produced by the electron beam. It was brought to our attention that Ohms law is not relativitically invariant but it is possible to make it so.
So my question is this, "Are the equations we have correct for a simple 1D model or do we need to change something?"
[tex]
\begin{array}{rcl}
\gamma^{3}(v/c)\left(\frac{\partial v}{\partial t}+v\frac{\partial v}{\partial x}\right) & = & -\frac{e}{m}E \\
\frac{\partial n}{\partial t}+\frac{\partial }{\partial x}(nv) & = & 0 \\
\frac{\partial E}{\partial t}+\frac{E}{\eta\epsilon_{0}} & = & \frac{nev}{\epsilon_{0}}
\end{array}
[/tex]
where e is the charge on the electron, m is the mass of the electron, [itex]\eta[/itex] is the resistivity of the charge neutral plasma, v is the speed of the electrons, n is the number density of the electrons and E is the electric field produced by the electron beam. It was brought to our attention that Ohms law is not relativitically invariant but it is possible to make it so.
So my question is this, "Are the equations we have correct for a simple 1D model or do we need to change something?"