Difference between revisions of "Electromagnetic Potentials"

From S.H.O.
Jump to: navigation, search
(Introduction)
(A New Idea: The Makerarc)
 
(146 intermediate revisions by the same user not shown)
Line 10: Line 10:
 
==Introduction==
 
==Introduction==
  
From December 2016 to March 2017, I ( ''[[User:S.H.O.|S.H.O.]] <sup>[[User_talk:S.H.O.|talk]]</sup>'') have been conducting electromagnetic simulations using JavaScript and the THREE.js script library (http://threejs.org). Based on these results, I have determined that the magnetic component of the Lorentz force:
+
From December 2016 to March 2017, I (''[[User:S.H.O.|S.H.O.]] <sup>[[User_talk:S.H.O.|talk]]</sup>'') have been conducting electromagnetic simulations using JavaScript and the THREE.js script library (http://threejs.org). Based on these results, I have determined that the magnetic component of the Lorentz force:
  
: <math>\mathbf{F_{mag}} = q \mathbf{v} \times \mathbf{B} = \nabla_\mathbf{A}(\mathbf{v}\cdot\mathbf{A}) - (\mathbf{v}\cdot\nabla)\mathbf{A}</math>
+
: <math>\mathbf{F_{mag}} = q\ \mathbf{v} \times \mathbf{B} = q\ \left[ \nabla_\mathbf{A}(\mathbf{v}\cdot\mathbf{A}) - (\mathbf{v}\cdot\nabla)\mathbf{A} \right] </math>
  
causes changes in kinetic energy that, in the low-frequency approximation, pretty much matches the amount of energy induced into surrounding coils.
+
causes transfers of energy ''within'' the magnetic rotor assembly (such as electron kinetic energy and inductive storage energy into wire kinetic energy (wired rotor example), or atomic electron kinetic energy to magnetic domain kinetic energy (permanent magnet rotor example)) that, in the low-frequency approximation, pretty much matches the amount of energy transfer between fields and the driving (stator) coils which is due to the transformer induction electric field <math>- \frac{\partial \mathbf{A}}{\partial t}</math> acting on currents in the coils from the relative motion of the magnetic rotor assembly.
  
As a result, the only way using conventional physics to explain various devices, such as the Marinov Motor,<ref>http://redshift.vif.com/JournalFiles/Pre2001/V05NO3PDF/v05n3phi.pdf</ref> the Distinti Paradox2<ref>http://www.distinti.com/docs/pdx/paradox2.pdf</ref>, and, especially, the Marinov Generator<ref name="Marinov Generator">http://overunity.com/14691/the-marinov-generator/</ref><ref name="Marinov Generator (paper)">http://www.overunityresearch.com/index.php?action=dlattach;topic=2470.0;attach=13897</ref> is to more accurately define the electric field <math>\mathbf{E}</math>, which is part of the full Lorentz Force equation, and then see if the leads to predictions that confirm observations, especially those found figure 9 of Cyril Smith's 2009 paper on the "Marinov Generator".<ref name="Marinov Generator (paper)"/> ''[[User:S.H.O.|S.H.O.]] <sup>[[User_talk:S.H.O.|talk]]</sup>'' 22:09, 5 March 2017 (PST)
+
As a result, the only way using conventional physics to explain various devices, such as the Marinov Motor<ref name="Marinov Motor">http://redshift.vif.com/JournalFiles/Pre2001/V05NO3PDF/v05n3phi.pdf</ref>, the Distinti Paradox2<ref>http://www.distinti.com/docs/pdx/paradox2.pdf</ref>, and, especially, the Marinov Generator<ref name="Marinov Generator">http://overunity.com/14691/the-marinov-generator/</ref><ref name="Marinov Generator (paper)">http://www.overunityresearch.com/index.php?action=dlattach;topic=2470.0;attach=13897</ref> is to more accurately define the electric field <math>\mathbf{E}</math>, which is part of the full Lorentz Force equation, and then see if it leads to predictions that confirm observations, especially those found figure 9 of Cyril Smith's 2009 paper on the "Marinov Generator".<ref name="Marinov Generator (paper)"/> ''[[User:S.H.O.|S.H.O.]] <sup>[[User_talk:S.H.O.|talk]]</sup>'' 22:09, 5 March 2017 (PST)
  
 
'''Prior content in the "Comment Record" section:'''
 
'''Prior content in the "Comment Record" section:'''
Line 29: Line 29:
 
==Novel Force laws proposed by various researchers==
 
==Novel Force laws proposed by various researchers==
  
===Stefan Marinov's proposal===
+
===James Wesley's proposal===
  
Stefan Marinov proposed adding the "motional-transformer induction" on charge q:<ref>https://archive.org/stream/thornywayoftruthpart4maririch#page/104/mode/2up/search/motional-transformer+induction</ref>
+
James Wesley proposed adding the "motional induction" on charge q. In SI Units, this can be expressed as:<ref name="Marinov Motor"/>
  
: <math>(\mathbf{v_{wire}}\cdot\nabla)\mathbf{A}</math>
+
: <math>-q(\mathbf{v}\cdot\nabla)\mathbf{A}</math>
  
or rather:
+
to the Lorentz force.
  
: <math>-(\mathbf{v_{charge}}\cdot\nabla)\mathbf{A}</math>
+
The idea behind this was to explain an observation in an experiment involving a "Marinov Motor"<ref name="Marinov Motor"/> in which longitudinal induction forces were produced.
 
+
to the Lorentz force.
+
  
 
The extra term is equivalent to:<ref>http://www.overunityresearch.com/index.php?action=dlattach;topic=2470.0;attach=13908</ref>
 
The extra term is equivalent to:<ref>http://www.overunityresearch.com/index.php?action=dlattach;topic=2470.0;attach=13908</ref>
Line 55: Line 53:
 
<math>\mathbf{F} = q\left[-\nabla \varphi- \frac{\partial \mathbf{A}}{\partial t}+ \nabla_\mathbf{A}(\mathbf{v}\cdot\mathbf{A})-(\mathbf{v}\cdot\nabla)\mathbf{A} \right]</math>
 
<math>\mathbf{F} = q\left[-\nabla \varphi- \frac{\partial \mathbf{A}}{\partial t}+ \nabla_\mathbf{A}(\mathbf{v}\cdot\mathbf{A})-(\mathbf{v}\cdot\nabla)\mathbf{A} \right]</math>
  
Therefore, adding the extra term proposed by Stefan Marinov results in:
+
Therefore, adding the extra term proposed by Wesley results in:
  
 
<math>\mathbf{F} = q\left[-\nabla \varphi- \frac{\partial \mathbf{A}}{\partial t}+ \nabla_\mathbf{A}(\mathbf{v}\cdot\mathbf{A})-2(\mathbf{v}\cdot\nabla)\mathbf{A} \right]</math>
 
<math>\mathbf{F} = q\left[-\nabla \varphi- \frac{\partial \mathbf{A}}{\partial t}+ \nabla_\mathbf{A}(\mathbf{v}\cdot\mathbf{A})-2(\mathbf{v}\cdot\nabla)\mathbf{A} \right]</math>
Line 61: Line 59:
 
The problem with this modification:
 
The problem with this modification:
  
In the case of an electrical charge approaching a wire, this additional term proposed by Marinov '''doubles''' the force of deflection. This is not observed.
+
In the case of an electrical charge approaching a wire, this additional term proposed by Wesley would '''double''' the force of deflection. This is not observed.
  
Consider the vector potential due a current-carrying wire on the x-axis. Both the current and the vector potential of this current point in the +x direction. Now have a charge approaching this wire perpendicularly. Both the magnetic Lorentz force and the additional Marinov term predict the same force. Now consider the case where the wire is moving toward the charge. In this case, both the transformer induction <math>- \frac{\partial \mathbf{A}}{\partial t}</math> and the additional Marinov term predict the same force acting on the charge directed parallel to the vector potential of the current. The term proposed by Marinov '''doubles''' the observed force in these cases.
+
Consider the vector potential due a current-carrying wire on the x-axis. Both the current and the vector potential of this current point in the <math>+x</math> direction. Now have a charge approaching this wire perpendicularly. Both the magnetic Lorentz force and the additional Wesley term predict the same force. The proposed additional term is superfluous and would double the force of deflection if added.
  
 
===Cyril Smith's proposal===
 
===Cyril Smith's proposal===
Line 91: Line 89:
 
The problem with this modification:
 
The problem with this modification:
  
In the case of two parallel current-carrying wires, this additional term proposed by Cyril Smith '''negates''' the magnetic forces between the currents. It turns out that the extra term may yield forces ''perpendicular'' to the velocity. These are:
+
In the case of two parallel current-carrying wires, this additional term proposed by Cyril Smith '''negates''' the magnetic forces between the currents. It turns out that the extra term may yield forces ''perpendicular'' to the velocity. The relevant field components are:
  
 
: <math>- \nabla_\mathbf{A}(\mathbf{v} \cdot \mathbf{A})_\bot = \begin{matrix}
 
: <math>- \nabla_\mathbf{A}(\mathbf{v} \cdot \mathbf{A})_\bot = \begin{matrix}
 
- \left[ v_y \frac{∂A_y}{∂x} + v_z \frac{∂A_z}{∂x} \right] \mathbf{e}_x \\
 
- \left[ v_y \frac{∂A_y}{∂x} + v_z \frac{∂A_z}{∂x} \right] \mathbf{e}_x \\
 
- \left[ v_x \frac{∂A_x}{∂y} + v_z \frac{∂A_z}{∂y} \right] \mathbf{e}_y \\
 
- \left[ v_x \frac{∂A_x}{∂y} + v_z \frac{∂A_z}{∂y} \right] \mathbf{e}_y \\
- \left[ v_x \frac{∂A_x}{∂z} + v_y \frac{∂A_y}{∂z} \right] \mathbf{e}_z  
+
- \left[ v_x \frac{∂A_x}{∂z} + v_y \frac{∂A_y}{∂z} \right] \mathbf{e}_z
 
\end{matrix}</math>
 
\end{matrix}</math>
  
==Explaining Cyril Smith's "Marinov Generator" using standard Lorentz Transforms==
+
==Explaining the Marinov Motor and Cyril Smith's "Marinov Generator" using Conventional Physics==
  
As it turns out, under both cases, the role of the electric scalar potential <math>\varphi</math> was ignored. Under a rather simple restriction however, the longitudinal force observed in the Marinov Generator experiment can emerge without any additional perpendicular components. The term we are looking for corresponds to <math>-v_x \frac{∂A_x}{∂x} \mathbf{e}_x</math> in the case that the Cartesian x-axis is defined to be aligned with the velocity <math>\mathbf{v}</math>.
+
{| class=wikitable width=200 style="float: right"
 +
|-
 +
|<youtube>https://www.youtube.com/watch?v=1TKSfAkWWN0</youtube>
 +
|-
 +
|
 +
See the online bulletin thread titled "Electromagnetism and relativity" for details on this video<ref>https://www.physicsforums.com/threads/electromagnetism-and-relativity.747811/</ref>.
 +
|}
  
The solution is to simply require that the gradient of the electric scalar potential (in the rest frame of the charge) <math>\varphi'</math> is zero.
+
The Liénard–Wiechert electric field was derived from the Liénard–Wiechert potentials <math>\varphi</math> and <math>A</math><ref name="Liénard"/><ref name="Wiechert"/> by Kirk T. McDonald<ref name="McDonald notes"/>, Professor Emeritus of Princeton University in New Jersey<ref>https://dof.princeton.edu/about/clerk-faculty/emeritus/kirk-t-mcdonald</ref>.
  
The Lorentz transforms of the electric scalar potential <math>\varphi</math> and magnetic vector potential <math>\mathbf{A}</math> are as follows:<ref name="Physics Formulas 2010">The Cambridge Handbook of Physics Formulas, G. Woan, Cambridge University Press, 2010, ISBN 978-0-521-57507-2.</ref>
+
In the sub-sections below, [https://en.wikipedia.org/wiki/Gaussian_units Gaussian units] are used unless otherwise noted. Also, do note that all electric fields due to '''''source charges''''' are, in the following sub-sections, evaluated in the <u>'''rest frame'''</u> of '''<u>each</u> ''target charge''''' subject to them <u>'''separately'''</u>. This procedure has roots in an approach to electromagnetism introduced by Edward M. Purcell<ref>http://physics.weber.edu/schroeder/mrr/mrrtalk.html</ref> in Section 5.6 of the Berkeley Physics Course (Volume II) titled ''Electricity and Magnetism''<ref>https://www.scribd.com/doc/128728926/Electricity-and-Magnetism-Berkeley-Physics-Course-Purcell</ref><ref>https://en.wikipedia.org/wiki/Berkeley_Physics_Course</ref> and is explained by a video by Veritasium titled "How Special Relativity Makes Magnets Work"<ref name="Veritasium">https://www.youtube.com/watch?v=1TKSfAkWWN0</ref>. This avoids having to perform calculations based on the magnetic field viewed by an arbitrary inertial observer. This procedure relies on the relative velocities between the charges. The calculations in the sub-sections below are valid for <math>v \ll c</math>.
  
: <math>\varphi' = \gamma (\varphi - v A_\parallel)</math>
+
===The Liénard–Wiechert electric fields for electrically-neutral currents===
: <math>A_\parallel' = \gamma (A_\parallel - v \varphi /c^2)</math>
+
: <math>A_\bot' = A_\bot</math>
+
  
By setting the gradient of the scalar potential in the charge's rest frame to zero, the result is:
+
From the paper titled "Onoochin's Paradox" by Kirk T. McDonald<ref name="McDonald">http://citeseerx.ist.psu.edu/viewdoc/summary?doi=10.1.1.299.8534</ref><ref>http://freeweb.siol.net/markoor/onoochin.pdf</ref>, we have following statement:
  
: <math>\nabla \varphi' = 0 = \gamma \nabla (\varphi - v A_\parallel)</math>
+
<blockquote>''For calculations of the Lorentz force to be accurate to order <math>\frac{1}{c^2}</math>, it suffices to use eq. (4) for the magnetic field. However, to maintain the desired accuracy the electric field of a moving charge must also include effects of retardation, as can be obtained from an expansion of the Liénard–Wiechert fields <ref name="Liénard">https://docs.google.com/file/d/0B817m31MAj0wZTZjZmMwMjgtY2Y5YS00YTQ5LThjM2EtNzhjYTYzNzFlZDY0/edit?hl=en_GB&pli=1</ref><ref name="Wiechert">https://docs.google.com/file/d/0B817m31MAj0wMDI1YjllYjctY2NhOS00M2M2LWFlMTUtYjVmYTkyZmVlY2M2/edit?hl=en_GB</ref> (for details, see the appendix of <ref name="McDonald notes">http://web.archive.org/web/20170318210550/http://puhep1.princeton.edu/~mcdonald/examples/ph501/ph501lecture24.pdf</ref>),
: <math>- \nabla \varphi = - v \nabla A_\parallel</math>
+
  
This can also be expressed as:
+
:<math>\mathbf{E} \approx q\ \frac{\mathbf{\hat{r}}}{r^2} \left(1 + \frac{v^2}{2c^2} - 3\frac{\left(\mathbf{\hat{r}} \cdot \mathbf{v}\right)^2}{2c^2} \right) - \frac{q}{2c^2r}\left[\mathbf{a}+\left(\mathbf{a} \cdot \mathbf{\hat{r}} \right)\mathbf{\hat{r}}\right]</math>
  
<math> - |\mathbf{v}| \nabla A_\parallel </math>
+
''where <math>\mathbf{a}</math> is the acceleration <math>\mathbf{a}</math> of the charge <math>q</math> at the present time.</blockquote>
  
Substituting this result in the Lorentz force yields:
+
Let's consider the situation  where the acceleration <math>\mathbf{a}</math> of charge <math>q</math> is negligible. The electric field at <math>\mathbf{r}</math> due to '''''source charge''''' <math>q</math> located at the origin is:
  
<math>\mathbf{F} = q \left[ - |\mathbf{v}| \nabla A_\parallel - \frac{\partial \mathbf{A}}{\partial t} + \nabla_{A}(\mathbf{v}\cdot\mathbf{A}) - (\mathbf{v}\cdot\nabla)\mathbf{A} \right]</math>
+
:<math>\mathbf{E} = q\ \frac{\mathbf{\hat{r}}}{r^2} \left(1 + \frac{v^2}{2c^2} - 3\frac{\left(\mathbf{\hat{r}} \cdot \mathbf{v}\right)^2}{2c^2} \right)</math>
  
In the case that the magnetic field is steady <math> - \frac{\partial \mathbf{A}}{\partial t} = 0</math> and magnetic Lorentz force on the wire is opposed by the physical restriction to the apparatus, the observed longitudinal induction field is:
+
In the Coulomb gauge, the first term in the parentheses comes from the electric scalar potential of a charge at rest in the observer's inertial frame. In event that the charge is contained within an electrically-neutral body, the electric field reduces to:
  
<math>\mathbf{E} = - |\mathbf{v}| \nabla A_\parallel</math>
+
:<math>\mathbf{E} = q\ \frac{\mathbf{\hat{r}}}{r^2} \left(\frac{v^2}{2c^2} - 3\frac{\left(\mathbf{\hat{r}} \cdot \mathbf{v}\right)^2}{2c^2} \right)</math>
  
Integrated over the path taken by the charge yields the same voltages that were predicted by Cyril Smith,<ref name="Marinov Generator (paper)"/> but here we use a correct expression and without modifying Einstein's theory of Special Relativity.
+
Let's consider a '''''target charge''''' <math>Q</math> located at <math>\mathbf{r}</math> at rest in the observer's inertial frame. The inertial observer and the charge <math>q</math> agree on what the electric field <math>\mathbf{E}</math> is, they agree that there is no magnetic force on <math>Q</math>, and finally, they agree on the acceleration of <math>Q</math>.
  
===The New Idea===
+
Another way to express this result is in terms of the angle <math>\theta</math> between <math>\mathbf{v}</math> and <math>\mathbf{r}</math>:
  
This new result has inspired a development beyond the S.H.O. Drive. A much more powerful and compact permanent magnet system is now envisioned, the viability of which is rests on the condition that <math>\langle \nabla \varphi' \rangle = 0</math> each charge <math>q</math>. The preliminary name is Makerarc (MAY-KERR-ARK), which stands for:
+
:<math>\mathbf{E} = q\ \frac{\mathbf{\hat{r}}}{r^2} \frac{v^2}{2c^2} \left(1 - 3\ cos^2\theta \right)</math>
 +
 
 +
The above equation can be broken up into two parts, one based on the '''relative ''azimuthal'' velocity''' <math>\mathbf{v}_{\theta}</math> of the source, and one based on the '''relative ''radial'' velocity''' <math>\mathbf{v}_r</math> of the source. First we rearrange the equation:
 +
 
 +
:<math>\mathbf{E} = q\ \frac{\mathbf{\hat{r}}}{r^2} \left(\frac{v^2 - \left(\mathbf{\hat{r}} \cdot \mathbf{v}\right)^2}{2c^2} - 2\frac{\left(\mathbf{\hat{r}} \cdot \mathbf{v}\right)^2}{2c^2} \right)</math>
 +
 
 +
Since the '''relative ''azimuthal'' velocity''' and the '''relative ''radial'' velocity''' are orthogonal, we can express the following:
 +
 
 +
:<math>\mathbf{E} = q\ \frac{\mathbf{\hat{r}}}{r^2} \left(\frac{v_{\theta}^2}{2c^2} - \frac{v_r^2}{c^2} \right)</math>
 +
 
 +
The '''relative ''azimuthal'' velocity''' can be split into orthogonal components in <math>x</math> and <math>y</math> according to the following Pythagorean relation:
 +
 
 +
:<math>|\mathbf{v}_{\theta}|^2 = |\mathbf{v}_x|^2 + |\mathbf{v}_y|^2</math>
 +
:<math>v_{\theta}^2 = v_x^2 + v_y^2</math>
 +
 
 +
Therefore:
 +
 
 +
:<math>\mathbf{E} = q\ \frac{\mathbf{\hat{r}}}{r^2} \left(\frac{v_x^2}{2c^2} + \frac{v_y^2}{2c^2} - \frac{v_r^2}{c^2} \right)</math>
 +
 
 +
The variables <math>v_x</math>, <math>v_y</math>, and <math>v_r</math> are mutually independent from each other. Therefore, the electric field <math>\mathbf{E}</math> can be split into three co-radial contributions:
 +
 
 +
:<math>\mathbf{E}_x = + q\ \frac{\mathbf{\hat{r}}}{r^2} \frac{v_x^2}{2c^2}</math>
 +
:<math>\mathbf{E}_y = + q\ \frac{\mathbf{\hat{r}}}{r^2} \frac{v_y^2}{2c^2}</math>
 +
:<math>\mathbf{E}_r = - q\ \frac{\mathbf{\hat{r}}}{r^2} \frac{v_r^2}{c^2}</math>
 +
 
 +
Consider the existence of four charges:
 +
 
 +
:<math>q_-</math> : the negative charge of loose electrons of the source current element
 +
:<math>Q_-</math> : the negative charge of loose electrons of the target current element
 +
:<math>q_+</math> : the positive charge of metallic atoms (excluding loose electrons) of the source current element
 +
:<math>Q_+</math> : the positive charge of metallic atoms (excluding loose electrons) of the target current element
 +
 
 +
Consider the existence of their corresponding current elements:
 +
 
 +
:<math>id\mathbf{l}</math> : the source current element
 +
:<math>Id\mathbf{L}</math> : the target current element
 +
 
 +
Where:
 +
 
 +
:<math>i</math> is the source current and <math>d\mathbf{l}</math> is its length element.
 +
:<math>I</math> is the target current and <math>d\mathbf{L}</math> is its length element.
 +
 
 +
The current elements are equal to:
 +
 
 +
:<math>id\mathbf{l} = q_- \mathbf{v}_d</math>
 +
:<math>Id\mathbf{L} = Q_- \mathbf{V}_d</math>
 +
 
 +
Where:
 +
 
 +
:<math>\mathbf{v}_d</math> is the drift velocity of the electrons of the source current element.
 +
:<math>\mathbf{V}_d</math> is the drift velocity of the electrons of the target current element.
 +
 
 +
Say we want to calculate the force on the target current element due to the source current element. This requires us to analyze four different forces:
 +
 
 +
:Force <math>\mathbf{F}_{--}</math> on <math>Q_-</math> by field <math>\mathbf{E}_{--}</math> of <math>q_-</math>
 +
:Force <math>\mathbf{F}_{+-}</math> on <math>Q_-</math> by field <math>\mathbf{E}_{+-}</math> of <math>q_+</math>
 +
:Force <math>\mathbf{F}_{-+}</math> on <math>Q_+</math> by field <math>\mathbf{E}_{-+}</math> of <math>q_-</math>
 +
:Force <math>\mathbf{F}_{++}</math> on <math>Q_+</math> by field <math>\mathbf{E}_{++}</math> of <math>q_+</math>
 +
 
 +
These four forces are dependent on four different relative velocities (source velocity w.r.t. target velocity):
 +
 
 +
:Relative velocity <math>\mathbf{v}_{--} = \mathbf{v}_- - \mathbf{V}_-</math>
 +
:Relative velocity <math>\mathbf{v}_{+-} = \mathbf{v}_+ - \mathbf{V}_-</math>
 +
:Relative velocity <math>\mathbf{v}_{-+} = \mathbf{v}_- - \mathbf{V}_+</math>
 +
:Relative velocity <math>\mathbf{v}_{++} = \mathbf{v}_+ - \mathbf{V}_+</math>
 +
 
 +
The corresponding relative speeds are:
 +
 
 +
:Relative speed <math>v_{x--} = \left(v_- - V_-\right)_x</math>
 +
:Relative speed <math>v_{x+-} = \left(v_+ - V_-\right)_x</math>
 +
:Relative speed <math>v_{x-+} = \left(v_- - V_+\right)_x</math>
 +
:Relative speed <math>v_{x++} = \left(v_+ - V_+\right)_x</math>
 +
 
 +
:Relative speed <math>v_{y--} = \left(v_- - V_-\right)_y</math>
 +
:Relative speed <math>v_{y+-} = \left(v_+ - V_-\right)_y</math>
 +
:Relative speed <math>v_{y-+} = \left(v_- - V_+\right)_y</math>
 +
:Relative speed <math>v_{y++} = \left(v_+ - V_+\right)_y</math>
 +
 
 +
:Relative speed <math>v_{r--} = \left(v_- - V_-\right)_r</math>
 +
:Relative speed <math>v_{r+-} = \left(v_+ - V_-\right)_r</math>
 +
:Relative speed <math>v_{r-+} = \left(v_- - V_+\right)_r</math>
 +
:Relative speed <math>v_{r++} = \left(v_+ - V_+\right)_r</math>
 +
 
 +
The two drift velocities are:
 +
 
 +
:<math>\mathbf{v}_d = \mathbf{v}_- - \mathbf{v}_+</math> : The drift velocity of the loose electrons of the source current element
 +
:<math>\mathbf{V}_d = \mathbf{V}_- - \mathbf{V}_+</math> : The drift velocity of the loose electrons of the target current element
 +
 
 +
Let the effective velocities of the current elements be:
 +
 
 +
:<math>\mathbf{v} = \left( \mathbf{v}_- + \mathbf{v}_+ \right)/2</math>
 +
:<math>\mathbf{V} = \left( \mathbf{V}_- + \mathbf{V}_+ \right)/2</math>
 +
 
 +
So the effective relative velocity between the current elements (source velocity w.r.t. target velocity) is:
 +
 
 +
:<math>\mathbf{v}_{rel} = \mathbf{v} - \mathbf{V}</math>
 +
 
 +
The effective velocity of each current element is halfway between the velocity of the negative charges and the velocity of the positive charges, so one may rather define a new variable, the deviation velocity, to be one-half the drift velocity of the electrons:
 +
 
 +
:<math>\mathbf{u} = \mathbf{v}_d / 2</math> is the deviation velocity of the source current element.
 +
:<math>\mathbf{U} = \mathbf{V}_d / 2</math> is the deviation velocity of the target current element.
 +
 
 +
All four relative velocities can be expressed in terms of the deviation velocities <math>\mathbf{u}</math> and <math>\mathbf{U}</math> together with the relative velocity <math>\mathbf{v}_{rel}</math>.
 +
 
 +
:Relative velocity <math>\mathbf{v}_{--} = \mathbf{v}_- - \mathbf{V}_- = \mathbf{u} + \mathbf{v}_{rel} - \mathbf{U}</math>
 +
:Relative velocity <math>\mathbf{v}_{+-} = \mathbf{v}_+ - \mathbf{V}_- = -\mathbf{u} + \mathbf{v}_{rel} - \mathbf{U}</math>
 +
:Relative velocity <math>\mathbf{v}_{-+} = \mathbf{v}_- - \mathbf{V}_+ = \mathbf{u} + \mathbf{v}_{rel} + \mathbf{U}</math>
 +
:Relative velocity <math>\mathbf{v}_{++} = \mathbf{v}_+ - \mathbf{V}_+ = -\mathbf{u} + \mathbf{v}_{rel} + \mathbf{U}</math>
 +
 
 +
The corresponding relative speeds are:
 +
 
 +
:Relative speed <math>v_{x--} = \left(v_- - V_- = u + v_{rel} - U\right)_x</math>
 +
:Relative speed <math>v_{x+-} = \left(v_+ - V_- = -u + v_{rel} - U\right)_x</math>
 +
:Relative speed <math>v_{x-+} = \left(v_- - V_+ = u + v_{rel} + U\right)_x</math>
 +
:Relative speed <math>v_{x++} = \left(v_+ - V_+ = -u + v_{rel} + U\right)_x</math>
 +
 
 +
:Relative speed <math>v_{y--} = \left(v_- - V_- = u + v_{rel} - U\right)_y</math>
 +
:Relative speed <math>v_{y+-} = \left(v_+ - V_- = -u + v_{rel} - U\right)_y</math>
 +
:Relative speed <math>v_{y-+} = \left(v_- - V_+ = u + v_{rel} + U\right)_y</math>
 +
:Relative speed <math>v_{y++} = \left(v_+ - V_+ = -u + v_{rel} + U\right)_y</math>
 +
 
 +
:Relative speed <math>v_{r--} = \left(v_- - V_- = u + v_{rel} - U\right)_r</math>
 +
:Relative speed <math>v_{r+-} = \left(v_+ - V_- = -u + v_{rel} - U\right)_r</math>
 +
:Relative speed <math>v_{r-+} = \left(v_- - V_+ = u + v_{rel} + U\right)_r</math>
 +
:Relative speed <math>v_{r++} = \left(v_+ - V_+ = -u + v_{rel} + U\right)_r</math>
 +
 
 +
The equation for the electric field contains squared values of the speed. As noted before:
 +
 
 +
:<math>\mathbf{E} = q\ \frac{\mathbf{\hat{r}}}{r^2} \left(\frac{v^2 - \left(\mathbf{\hat{r}} \cdot \mathbf{v}\right)^2}{2c^2} - 2\frac{\left(\mathbf{\hat{r}} \cdot \mathbf{v}\right)^2}{2c^2} \right)</math>
 +
 
 +
The electric fields on charges <math>Q_+</math> and <math>Q_-</math>, in their own respective and distinct rest frames, are as follows:
 +
 
 +
:<math>\mathbf{E}_- = \frac{\mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( q_- \left(v_{x--}^2 + v_{y--}^2 - 2 v_{r--}^2\right) + q_+ \left(v_{x+-}^2 + v_{y+-}^2 - 2 v_{r+-}^2\right) \right)</math>
 +
:<math>\mathbf{E}_+ = \frac{\mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( q_- \left(v_{x-+}^2 + v_{y-+}^2 - 2 v_{r-+}^2\right) + q_+ \left(v_{x++}^2 + v_{y++}^2 - 2 v_{r++}^2\right) \right)</math>
 +
 
 +
===The Forces due to the Liénard–Wiechert electric fields for electrically-neutral currents===
 +
 
 +
It may be more helpful to calculate the forces due to relative radial velocities separately from the forces due to relative azimuthal velocities. Therefore:
 +
 
 +
:<math>\mathbf{F}_x = \frac{\mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( Q_- q_- v_{x--}^2 + Q_- q_+ v_{x+-}^2 + Q_+ q_- v_{x-+}^2 + Q_+ q_+ v_{x++}^2 \right)</math>
 +
:<math>\mathbf{F}_y = \frac{\mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( Q_- q_- v_{y--}^2 + Q_- q_+ v_{y+-}^2 + Q_+ q_- v_{y-+}^2 + Q_+ q_+ v_{y++}^2 \right)</math>
 +
:<math>\mathbf{F}_r = \frac{\mathbf{\hat{r}}}{r^2} \frac{- 2}{2c^2} \left( Q_- q_- v_{r--}^2 + Q_- q_+ v_{r+-}^2 + Q_+ q_- v_{r-+}^2 + Q_+ q_+ v_{r++}^2 \right)</math>
 +
 
 +
It follows that <math>\mathbf{F}_x</math> and <math>\mathbf{F}_y</math> are functions of currents perpendicular to radial vector <math>\mathbf{r}</math> while <math>\mathbf{F}_r</math> is a function of currents co-linear with radial vector <math>\mathbf{r}</math>
 +
 
 +
All charges (<math>Q_-</math>, <math>Q_+</math>, <math>q_-</math>, and <math>q_+</math>) may contribute simultaneously to the azimuthally-directed ('''transverse''') currents in <math>x</math> and <math>y</math> and the radially-directed ('''longitudinal''') currents in <math>r</math>.
 +
 
 +
It will be very advantageous to simplify these formulas. For electrically-neutral currents, we can recognize the following:
 +
 
 +
:<math>q_+ = - q_-</math>
 +
:<math>Q_+ = - Q_-</math>
 +
 
 +
This allows us to factor out <math>Q_+ q_+</math> with the following result:
 +
 
 +
:<math>\mathbf{F}_x = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( v_{x--}^2 - v_{x+-}^2 - v_{x-+}^2 + v_{x++}^2 \right)</math>
 +
:<math>\mathbf{F}_y = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( v_{y--}^2 - v_{y+-}^2 - v_{y-+}^2 + v_{y++}^2 \right)</math>
 +
:<math>\mathbf{F}_r = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{- 2}{2c^2} \left( v_{r--}^2 - v_{r+-}^2 - v_{r-+}^2 + v_{r++}^2 \right)</math>
 +
 
 +
Next, we will work on simplifying the contents within the parentheses. To make matters simpler, we will move the subscript to the lower right corner of the parentheses.
 +
 
 +
:<math>\mathbf{F}_x = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 \right)_x</math>
 +
:<math>\mathbf{F}_y = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 \right)_y</math>
 +
:<math>\mathbf{F}_r = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{- 2}{2c^2} \left( v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 \right)_r</math>
 +
 
 +
These can be distributed back to the contents within the parentheses after they are substituted for a different expression in terms of the currents. The terms in the parentheses (disregarding their sign) are as follows:
 +
 
 +
:<math>v_{--}^2 = \left(u + v_{rel} - U\right)^2 = u^2 + v_{rel}^2 + U^2 + (2u - 2U) v_{rel} - 2 u U</math>
 +
:<math>v_{+-}^2 = \left(-u + v_{rel} - U\right)^2 = u^2 + v_{rel}^2 + U^2 + (-2u - 2U) v_{rel} + 2 u U</math>
 +
:<math>v_{-+}^2 = \left(u + v_{rel} + U\right)^2 = u^2 + v_{rel}^2 + U^2 + (2u + 2U) v_{rel} + 2 u U</math>
 +
:<math>v_{++}^2 = \left(-u + v_{rel} + U\right)^2 = u^2 + v_{rel}^2 + U^2 + (-2u + 2U) v_{rel} - 2 u U</math>
 +
 
 +
Therefore sum of the forces on <math>Q_-</math> and <math>Q_+</math> depends on:
 +
 
 +
:<math>v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 = ((2u - 2U) - (-2u - 2U) - (2u + 2U) + (-2u + 2U))v_{rel} + ((-2) - (2) - (2) + (-2)) u U</math>
 +
:<math>v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 = ((2u - 2U) + (2u + 2U) - (2u + 2U) - (2u - 2U))v_{rel} - 8 u U</math>
 +
:<math>v_{--}^2 - v_{+-}^2 - v_{-+}^2 - v_{+-}^2 + v_{++}^2 = ((-2U) + (2U) - (2U) - (- 2U))v_{rel} - 8 u U</math>
 +
:<math>v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 = (0)v_{rel} - 8 u U</math>
 +
:<math>v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 = - 8 u U</math>
 +
:<math>v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 = - 8 (v_d/2)(V_d/2)</math>
 +
:<math>v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 = - 2 (v_d)(V_d)</math>
 +
 
 +
The force on <math>Q_-</math> depends on:
 +
 
 +
:<math>v_{--}^2 - v_{+-}^2 = ((2u - 2U) - (-2u - 2U))v_{rel} + ((-2) - (2)) u U</math>
 +
:<math>v_{--}^2 - v_{+-}^2 = 4 u v_{rel} - 4 u U</math>
 +
 
 +
The force on <math>Q_+</math> depends on:
 +
 
 +
:<math>v_{++}^2 - v_{-+}^2 = ((-2u + 2U) - (2u + 2U))v_{rel} + (- (2) + (-2)) u U</math>
 +
:<math>v_{++}^2 - v_{-+}^2 = -4 u v_{rel} - 4 u U</math>
 +
 
 +
Recalling that:
 +
 
 +
:<math>\mathbf{F}_x = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 \right)_x</math>
 +
:<math>\mathbf{F}_y = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 \right)_y</math>
 +
:<math>\mathbf{F}_r = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{- 2}{2c^2} \left( v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 \right)_r</math>
 +
:<math>v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 = - 2 (v_d)(V_d)</math>
 +
 
 +
Substitution yields:
 +
 
 +
:<math>\mathbf{F}_x = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( - 2 (v_d)(V_d) \right)_x</math>
 +
:<math>\mathbf{F}_y = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( - 2 (v_d)(V_d) \right)_y</math>
 +
:<math>\mathbf{F}_r = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{- 2}{2c^2} \left( - 2 (v_d)(V_d) \right)_r</math>
 +
 
 +
Simplified:
 +
 
 +
:<math>\mathbf{F}_x = - \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( (v_d)(V_d) \right)_x</math>
 +
:<math>\mathbf{F}_y = - \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( (v_d)(V_d) \right)_y</math>
 +
:<math>\mathbf{F}_r = + \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{2}{c^2} \left( (v_d)(V_d) \right)_r</math>
 +
 
 +
As stated earlier, the current elements are equal to:
 +
 
 +
:<math>id\mathbf{l} = q_- \mathbf{v}_d</math>
 +
:<math>Id\mathbf{L} = Q_- \mathbf{V}_d</math>
 +
 
 +
We can now substitute the currents into the equation. First we substitute <math>Q_- q_-</math> for <math>Q_+ q_+</math>:
 +
 
 +
:<math>\mathbf{F}_x = - \frac{Q_- q_- \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( (v_d)(V_d) \right)_x</math>
 +
:<math>\mathbf{F}_y = - \frac{Q_- q_- \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( (v_d)(V_d) \right)_y</math>
 +
:<math>\mathbf{F}_r = + \frac{Q_- q_- \mathbf{\hat{r}}}{r^2} \frac{2}{c^2} \left( (v_d)(V_d) \right)_r</math>
 +
 
 +
Next, we assign each charge with their corresponding drift velocities:
 +
 
 +
:<math>\mathbf{F}_x = - \frac{\mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( (Q_- V_d) (q_- v_d) \right)_x</math>
 +
:<math>\mathbf{F}_y = - \frac{\mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( (Q_- V_d) (q_- v_d) \right)_y</math>
 +
:<math>\mathbf{F}_r = + \frac{\mathbf{\hat{r}}}{r^2} \frac{2}{c^2} \left( (Q_- V_d) (q_- v_d) \right)_r</math>
 +
 
 +
Next, we make a substitution for the current elements:
 +
 
 +
:<math>\mathbf{F}_x = - \frac{\mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( (Id\mathbf{L}) \cdot (id\mathbf{l}) \right)_x</math>
 +
:<math>\mathbf{F}_y = - \frac{\mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( (Id\mathbf{L}) \cdot (id\mathbf{l}) \right)_y</math>
 +
:<math>\mathbf{F}_r = + \frac{\mathbf{\hat{r}}}{r^2} \frac{2}{c^2} \left( (Id\mathbf{L}) \cdot (id\mathbf{l}) \right)_r</math>
 +
 
 +
This can be written as:
 +
 
 +
:<math>\mathbf{F}_x = - \frac{\mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( Ii (d\mathbf{L} \cdot d\mathbf{l}) \right)_x</math>
 +
:<math>\mathbf{F}_y = - \frac{\mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( Ii (d\mathbf{L} \cdot d\mathbf{l}) \right)_y</math>
 +
:<math>\mathbf{F}_r = + \frac{\mathbf{\hat{r}}}{r^2} \frac{2}{c^2} \left( Ii (d\mathbf{L} \cdot d\mathbf{l}) \right)_r</math>
 +
 
 +
Or:
 +
 
 +
:<math>\mathbf{F}_x = - \frac{\mathbf{\hat{r}}}{r^2} \frac{Ii}{c^2} \left( d\mathbf{L} \cdot d\mathbf{l} \right)_x</math>
 +
:<math>\mathbf{F}_y = - \frac{\mathbf{\hat{r}}}{r^2} \frac{Ii}{c^2} \left( d\mathbf{L} \cdot d\mathbf{l} \right)_y</math>
 +
:<math>\mathbf{F}_r = + \frac{\mathbf{\hat{r}}}{r^2} \frac{2Ii}{c^2} \left( d\mathbf{L} \cdot d\mathbf{l} \right)_r</math>
 +
 
 +
Or:
 +
 
 +
:<math>\mathbf{F}_x = - \frac{\mathbf{\hat{r}}}{r^2} \frac{Ii}{c^2} \left( d\mathbf{L}_x \cdot d\mathbf{l}_x \right)</math>
 +
:<math>\mathbf{F}_y = - \frac{\mathbf{\hat{r}}}{r^2} \frac{Ii}{c^2} \left( d\mathbf{L}_y \cdot d\mathbf{l}_y \right)</math>
 +
:<math>\mathbf{F}_r = + \frac{\mathbf{\hat{r}}}{r^2} \frac{2Ii}{c^2} \left( d\mathbf{L}_r \cdot d\mathbf{l}_r \right)</math>
 +
 
 +
Adding the forces together, their sum is:
 +
 
 +
:<math>\mathbf{F} = \frac{\mathbf{\hat{r}}}{r^2} \frac{Ii}{c^2} \left( 2 d\mathbf{L}_r \cdot d\mathbf{l}_r - d\mathbf{L}_x \cdot d\mathbf{l}_x - d\mathbf{L}_y \cdot d\mathbf{l}_y \right)</math>
 +
:<math>\mathbf{F} = \frac{\mathbf{\hat{r}}}{r^2} \frac{Ii}{c^2} \left( 2 d\mathbf{L}_r \cdot d\mathbf{l}_r - d\mathbf{L}_\theta \cdot d\mathbf{l}_\theta \right)</math>
 +
 
 +
In S.I. Units:
 +
 
 +
:<math>\mathbf{F} = \frac{\mathbf{\hat{r}}}{4\pi\epsilon_0r^2} \frac{Ii}{c^2} \left( 2 d\mathbf{L}_r \cdot d\mathbf{l}_r - d\mathbf{L}_\theta \cdot d\mathbf{l}_\theta \right)</math>
 +
 
 +
Or:
 +
 
 +
:<math>\mathbf{F} = \frac{\mu_0 I i \mathbf{\hat{r}}}{4\pi r^2} \left( 2 d\mathbf{L}_r \cdot d\mathbf{l}_r - d\mathbf{L}_\theta \cdot d\mathbf{l}_\theta \right)</math>
 +
 
 +
The second differential is:
 +
 
 +
:<math>d^2\mathbf{F} = \frac{\mu_0 \mathbf{\hat{r}}}{4\pi r^2} \left( 2 \mathbf{I}_r \cdot \mathbf{i}_r - \mathbf{I}_\theta \cdot \mathbf{i}_\theta \right)</math>
 +
 
 +
===Special Scenario: No relative motion between positive charges===
 +
 
 +
In the case that the positive charges <math>q_+</math> and <math>Q_+</math> are essentially stationary, '''we can simplify the derivation of the field and force equations in the rest frame of the positive charges'''.
 +
 
 +
As stated in the parent section, the two drift velocities are:
 +
 
 +
:<math>\mathbf{v}_d = \mathbf{v}_- - \mathbf{v}_+</math> : The drift velocity of the loose electrons of the source current element
 +
:<math>\mathbf{V}_d = \mathbf{V}_- - \mathbf{V}_+</math> : The drift velocity of the loose electrons of the target current element
 +
 
 +
These become:
 +
 
 +
:<math>\mathbf{v}_d = \mathbf{v}_-</math> : The drift velocity of the loose electrons of the source current element
 +
:<math>\mathbf{V}_d = \mathbf{V}_-</math> : The drift velocity of the loose electrons of the target current element
 +
 
 +
As stated in the parent section, the effective velocities of the current elements are:
 +
 
 +
:<math>\mathbf{v} = \left( \mathbf{v}_- + \mathbf{v}_+ \right)/2</math>
 +
:<math>\mathbf{V} = \left( \mathbf{V}_- + \mathbf{V}_+ \right)/2</math>
 +
 
 +
These become:
 +
 
 +
:<math>\mathbf{v} = \mathbf{v}_-  /2</math>
 +
:<math>\mathbf{V} = \mathbf{V}_- / 2</math>
 +
 
 +
As stated in the parent section, a new variable, the deviation velocity, is defined as one-half the drift velocity of the electrons:
 +
 
 +
:<math>\mathbf{u} = \mathbf{v}_d / 2</math> is the deviation velocity of the source current element.
 +
:<math>\mathbf{U} = \mathbf{V}_d / 2</math> is the deviation velocity of the target current element.
 +
 
 +
Therefore:
 +
 
 +
:<math>\mathbf{v} = \mathbf{u} </math>
 +
:<math>\mathbf{V} = \mathbf{U} </math>
 +
 
 +
The effective relative velocity between the current elements (source velocity w.r.t. target velocity) is:
 +
 
 +
:<math>\mathbf{v}_{rel} = \mathbf{v} - \mathbf{V}</math>
 +
 
 +
Therefore:
 +
 
 +
:<math>\mathbf{v}_{rel} = \mathbf{u} - \mathbf{U}</math>
 +
 
 +
As derived in the parent section, all four relative velocities can be expressed in terms of the deviation velocities <math>\mathbf{u}</math> and <math>\mathbf{U}</math> together with the relative velocity <math>\mathbf{v}_{rel}</math>.
 +
 
 +
:Relative velocity <math>\mathbf{v}_{--} = \mathbf{v}_- - \mathbf{V}_- = \mathbf{u} + \mathbf{v}_{rel} - \mathbf{U}</math>
 +
:Relative velocity <math>\mathbf{v}_{-+} = \mathbf{v}_- - \mathbf{V}_+ = \mathbf{u} + \mathbf{v}_{rel} + \mathbf{U}</math>
 +
:Relative velocity <math>\mathbf{v}_{+-} = \mathbf{v}_+ - \mathbf{V}_- = -\mathbf{u} + \mathbf{v}_{rel} - \mathbf{U}</math>
 +
:Relative velocity <math>\mathbf{v}_{++} = \mathbf{v}_+ - \mathbf{V}_+ = -\mathbf{u} + \mathbf{v}_{rel} + \mathbf{U}</math>
 +
 
 +
By substituting for <math>\mathbf{v}_{rel}</math>, we get:
 +
 
 +
:Relative velocity <math>\mathbf{v}_{--} = \mathbf{v}_- - \mathbf{V}_- = 2\mathbf{u} - 2\mathbf{U}</math>
 +
:Relative velocity <math>\mathbf{v}_{-+} = \mathbf{v}_- - \mathbf{V}_+ = 2\mathbf{u}</math>
 +
:Relative velocity <math>\mathbf{v}_{+-} = \mathbf{v}_+ - \mathbf{V}_- = -2\mathbf{U}</math>
 +
:Relative velocity <math>\mathbf{v}_{++} = \mathbf{v}_+ - \mathbf{V}_+ = 0</math>
 +
 
 +
As derived in the parent section, the following is a series sum of terms (a function of relative speeds) which will be used to help calculate the forces between currents:
 +
 
 +
:<math>v_{--}^2 - v_{-+}^2 - v_{+-}^2 + v_{++}^2</math>
 +
 
 +
Based on the above results, we have:
 +
 
 +
:<math>v_{--}^2 - v_{-+}^2 - v_{+-}^2 + v_{++}^2 = (2u - 2U)^2 - (2u)^2 - (-2U)^2 + 0</math>
 +
:<math>v_{--}^2 - v_{-+}^2 - v_{+-}^2 + v_{++}^2 = (4u^2 + 4U^2 - 8uU) - 4u^2 - 4U^2 + 0</math>
 +
:<math>v_{--}^2 - v_{-+}^2 - v_{+-}^2 + v_{++}^2 = - 8uU</math>
 +
:<math>v_{--}^2 - v_{-+}^2 - v_{+-}^2 + v_{++}^2 = - 8 (v_{d}/2)(V_{d}/2)</math>
 +
:<math>v_{--}^2 - v_{-+}^2 - v_{+-}^2 + v_{++}^2 =  -2 (v_{d})(V_{d})</math>
 +
 
 +
This produces the same results as the parent section titled "Explaining the Marinov Motor and Cyril Smith's 'Marinov Generator' using Conventional Physics" and the previous section titled "The Forces due to the Liénard–Wiechert electric fields for electrically-neutral currents":
 +
 
 +
In S.I. Units:
 +
 
 +
:<math>\mathbf{F} = \frac{\mathbf{\hat{r}}}{4\pi\epsilon_0r^2} \frac{Ii}{c^2} \left( 2 d\mathbf{L}_r \cdot d\mathbf{l}_r - d\mathbf{L}_\theta \cdot d\mathbf{l}_\theta \right)</math>
 +
 
 +
Or:
 +
 
 +
:<math>\mathbf{F} = \frac{\mu_0 I i \mathbf{\hat{r}}}{4\pi r^2} \left( 2 d\mathbf{L}_r \cdot d\mathbf{l}_r - d\mathbf{L}_\theta \cdot d\mathbf{l}_\theta \right)</math>
 +
 
 +
The second differential is:
 +
 
 +
:<math>d^2\mathbf{F} = \frac{\mu_0 \mathbf{\hat{r}}}{4\pi r^2} \left( 2 \mathbf{I}_r \cdot \mathbf{i}_r - \mathbf{I}_\theta \cdot \mathbf{i}_\theta \right)</math>
 +
 
 +
===Deriving James Wesley's additional force term in the case of co-linear current elements===
 +
 
 +
As pointed out above in the sub-section titled "James Wesley's proposal", James Wesley proposed adding the "motional induction" on charge <math>q</math>. In SI Units, this can be expressed as:<ref name="Marinov Motor"/>
 +
 
 +
: <math>-q(\mathbf{v}\cdot\nabla)\mathbf{A}</math>
 +
 
 +
How may we derive the same result for current elements located and oriented along the same line? Equation 2b of the paper titled "Observations of the Marinov Motor"<ref name="Marinov Motor"/> can be adapted to the form above, resulting in:
 +
 
 +
: <math>\mathbf{F} = -Q_- \left(\mathbf{V}_d \cdot \nabla \right) \mathbf{A} = Q_- \left(-\mathbf{V}_d \cdot \nabla \right) \mathbf{A} = - I \left( d\mathbf{L} \cdot \nabla \right) \mathbf{A} </math>
 +
 
 +
In the case that we are dealing only with co-linear current elements, the result can be expressed as:
 +
 
 +
: <math>\mathbf{F}_r = \left( -Q_- \left(\mathbf{V}_d \cdot \nabla \right) \mathbf{A} = Q_- \left(-\mathbf{V}_d \cdot \nabla \right) \mathbf{A} = - I \left( d\mathbf{L} \cdot \nabla \right) \mathbf{A} \right)_r</math>
 +
 
 +
As stated in the parent section, the force on <math>Q_-</math> depends on:
 +
 
 +
:<math>v_{--}^2 - v_{+-}^2 = ((2u - 2U) - (-2u - 2U))v_{rel} + ((-2) - (2)) u U</math>
 +
:<math>v_{--}^2 - v_{+-}^2 = 4 u v_{rel} - 4 u U</math>
 +
 
 +
As stated in the parent section, The force on <math>Q_+</math> depends on:
 +
 
 +
:<math>v_{++}^2 - v_{-+}^2 = ((-2u + 2U) - (2u + 2U))v_{rel} + (- (2) + (-2)) u U</math>
 +
:<math>v_{++}^2 - v_{-+}^2 = -4 u v_{rel} - 4 u U</math>
 +
 
 +
In the case that that there is no relative velocity <math>v_{rel}</math> between the current elements:
 +
 
 +
:<math>v_{--}^2 - v_{+-}^2 = - 4 u U = -4 (v_d/2)(V_d/2) = - (v_d)(V_d)</math>
 +
:<math>v_{++}^2 - v_{-+}^2 = - 4 u U = -4 (v_d/2)(V_d/2) = - (v_d)(V_d)</math>
 +
 
 +
As stated in the section "The Forces due to the Liénard–Wiechert electric fields for electrically-neutral currents", the electric forces on charges <math>Q_-</math> and <math>Q_+</math> are:
 +
 
 +
:<math>\mathbf{F}_{r-} = + \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( v_{--}^2 - v_{+-}^2  \right)_r</math>
 +
:<math>\mathbf{F}_{r+} = + \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( v_{++}^2 - v_{-+}^2  \right)_r</math>
 +
 
 +
From this and above, it follows that:
 +
 
 +
:<math>\mathbf{F}_{r-} = + \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( - (v_d)(V_d)  \right)_r</math>
 +
:<math>\mathbf{F}_{r+} = + \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( - (v_d)(V_d)  \right)_r</math>
 +
 
 +
Therefore:
 +
 
 +
:<math>\mathbf{F}_{r-} = + \frac{Q_- q_- \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( - (v_d)(V_d) \right)_r</math>
 +
:<math>\mathbf{F}_{r+} = + \frac{Q_+ q_- \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( + (v_d)(V_d) \right)_r</math>
 +
 
 +
The effective electric fields experience by each charge is:
 +
 
 +
:<math>\mathbf{E}_{r-} = + \frac{q_- \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( - (v_d)(V_d) \right)_r</math>
 +
:<math>\mathbf{E}_{r+} = + \frac{q_- \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( + (v_d)(V_d) \right)_r</math>
 +
 
 +
Given that:
 +
 
 +
:<math>Q_+ = - Q_-</math>
 +
:<math>\mathbf{U} = \mathbf{v}_d / 2</math>
 +
:<math>\mathbf{v}_d = \mathbf{U} - (-\mathbf{U})</math>
 +
:<math>Id\mathbf{L} = Q_- \mathbf{V}_d = Q_- (+ \mathbf{U}) + Q_+ (-\mathbf{U})</math>
 +
 
 +
It can be shown that:
 +
 
 +
:<math>\mathbf{F}_{r+} = \left( - Q_+ \left((-\mathbf{U}) \cdot \nabla \right) \mathbf{A} \right)_r</math>
 +
:<math>\mathbf{F}_{r-} = \left( - Q_- \left((+\mathbf{U}) \cdot \nabla \right) \mathbf{A} \right)_r</math>
 +
:<math>\mathbf{F}_r = \left( - Q_- \left(\mathbf{V}_d \cdot \nabla \right) \mathbf{A} \right)_r = \mathbf{F}_{r-} + \mathbf{F}_{r+}</math>
 +
 
 +
This is based on the directional derivative of vector potential <math>\mathbf{A}</math> along <math>Id\mathbf{L}</math>.
 +
 
 +
It is important to note that the total force <math>\mathbf{F}_r</math> in <math>r</math> remains unaffected by the relative motion between the target current element <math>I_rd\mathbf{L}_r</math> and the source current element <math>i_rd\mathbf{l}_r</math> responsible for vector potential <math>\mathbf{A}_r</math>. However, forces <math>\mathbf{F}_{r-}</math> and <math>\mathbf{F}_{r+}</math> are affected by the relative speed <math>(v_{rel})_r</math> between the current elements in <math>r</math>, which leads to a kind of electromagnetic induction acting between co-linear current elements in relative motion in <math>r</math>. This feature may help to explain Cyril Smith's "Marinov Generator"<ref name="Marinov Generator"/>.
 +
 
 +
==A New Idea: The Makerarc==
 +
 
 +
The above result supports a development beyond the S.H.O. Drive. A much more powerful and compact permanent magnet system is now envisioned.
 +
 
 +
The preliminary name is Makerarc (MAY-KERR-ARK), which stands for:
  
 
* '''M'''agnetic
 
* '''M'''agnetic
Line 143: Line 564:
 
* '''C'''hannel
 
* '''C'''hannel
  
Details pending. Stay turned. Sincerely,  ''[[User:S.H.O.|S.H.O.]] <sup>[[User_talk:S.H.O.|talk]]</sup>'' 22:09, 5 March 2017 (PST)
+
Details pending. Stay tuned. Sincerely,  ''[[User:S.H.O.|S.H.O.]] <sup>[[User_talk:S.H.O.|talk]]</sup>'' 22:09, 5 March 2017 (PST)
 +
 
 +
What makes a magnet tick? Clues can be found in the video titled "MAGNETS: How Do They Work?" by Veritasium<ref>https://www.youtube.com/watch?v=hFAOXdXZ5TM</ref>. ''[[User:S.H.O.|S.H.O.]] <sup>[[User_talk:S.H.O.|talk]]</sup>'' 22:08, 23 March 2017 (PDT)
 +
 
 +
<youtube>https://www.youtube.com/watch?v=hFAOXdXZ5TM</youtube>
 +
 
 +
An attempt to produce the Makerarc back in March 2017<ref>https://www.facebook.com/Sho.Drives/videos/790278744462849/</ref> failed to demonstrate the predictions of the term <math>\mathbf{F}_r = \left( - Q_- \left(\mathbf{V}_d \cdot \nabla \right) \mathbf{A} \right)_r</math><ref>https://www.facebook.com/Sho.Drives/videos/791514291005961/</ref> where I assumed that an external vector potential <math>\mathbf{A}</math> would lead to a force co-linear to current elements. This may be due to the supercurrent nature of electrons bounded to atoms, allowing atoms to be capable of preserving or "fixing" the amount of magnetic flux that passes through them, just as "macroscopic" superconducting currents do for superconductors, and therefore, by extension be capable of resisting changes in the magnetic vector potential. This does not negate the possibility of a Marinov Generator, as designed by Cyril Smith, because in his case the currents receiving power were inside a conductor where an externally applied vector potential is not fully shielded against, permitting <math>\mathbf{F}_r = \left( - Q_- \left(\mathbf{V}_d \cdot \nabla \right) \mathbf{A} \right)_r</math> to yield a force. At the moment, I am still trying to decide on the correct model, but hopefully the surviving concept is similar to what I have developed so far in this page. Sincerely,  ''[[User:S.H.O.|S.H.O.]] <sup>[[User_talk:S.H.O.|talk]]</sup>'' 19:55, 14 July 2017 (PDT)
 +
 
 +
==Explaining "Altered" Lenz' Law Devices==
 +
 
 +
It is anticipated that the longitudinal force described above in the section titled "Explaining the Marinov Motor and Cyril Smith's 'Marinov Generator' using Conventional Physics" may explain some types of purported Reduced-Lenz devices. A good example can be found in the video below ("Ray's No Back EMF Generator"), although in this example, the longitudinal force ''increases'' the drag, mainly in positions where pancake generator coil is mostly outside the cylindrical boundary of the permanent magnet. This creates an illusion of a "Reduced" Lenz' Law effect when the magnet is mostly within the cylindrical boundary of the permanent magnet:
 +
 
 +
<youtube>x3TFALmHtMw</youtube>
 +
 
 +
Simulations in JavaScript and THREE.js have determined that in many other configurations of currents and magnets, the magnetic Lorentz forces <math>q\ \mathbf{v} \times \mathbf{B}</math> will be ''opposed'' in part by the additional force. This makes certain types of magnetic circuit arrangements more efficient at electrical power generation (in terms of output power vs. input power) but less efficient as an electrical motor. Simulations of the S.H.O. Drive indicated that the additional force may vary between assisting or opposing the magnetic Lorentz force, which is helpful depending on whether the goal is to develop an efficient motor or an efficient generator, respectively. The simulated S.H.O. Drive designs appeared to switch between these two extremes every quarter cycle. The only way to make it work would have been to increase the inductive reactance (so that the current would be delayed by about a quarter cycle), and even then, simulations indicated that the extra term was typically insufficient to completely reverse the net force. However simulations of various rotor and S.H.O. coil curve modifications for the S.H.O. Drive showed that it was possible for the additional force to be a significant percentage of the magnetic Lorentz force. Per more recent simulations (early March), the Makerarc design (previous section) will improve upon this many fold. ''[[User:S.H.O.|S.H.O.]] <sup>[[User_talk:S.H.O.|talk]]</sup>'' 23:39, 5 March 2017 (PST)
 +
 
 +
==Explaining the Newman Motor==
 +
 
 +
{{See also|Memory Lane}}
 +
 
 +
A Newman Motor-style coil and magnet arrangement, like that shown in the video below, have been simulated by me using JavaScript and THREE.js.
 +
 
 +
<youtube>APq8HaeJpQ8</youtube>
 +
 
 +
The extra electric field term predicts a significant opposition to the magnetic Lorentz force at angles slightly straying from the "top-dead-vertical" position, making it a better generator than a motor. However, when energy is discharged from the "generator coil" to the "motor coil" of Newman's motor, the rotor will have often changed position to the point where the magnetic Lorentz force becomes increasingly significant, helpful for motive purposes. Newman's motor operated at a high Q, which facilitated energy recovery. ''[[User:S.H.O.|S.H.O.]] <sup>[[User_talk:S.H.O.|talk]]</sup>'' 00:00, 6 March 2017 (PST)
  
 
==References==
 
==References==

Latest revision as of 18:55, 14 July 2017

The basic idea here is that the electromagnetic potentials [math]\varphi[/math] and [math]A[/math] may be the underlying key to several inventions related to electromagnetic forces.

Introduction

From December 2016 to March 2017, I (S.H.O. talk) have been conducting electromagnetic simulations using JavaScript and the THREE.js script library (http://threejs.org). Based on these results, I have determined that the magnetic component of the Lorentz force:

[math]\mathbf{F_{mag}} = q\ \mathbf{v} \times \mathbf{B} = q\ \left[ \nabla_\mathbf{A}(\mathbf{v}\cdot\mathbf{A}) - (\mathbf{v}\cdot\nabla)\mathbf{A} \right] [/math]

causes transfers of energy within the magnetic rotor assembly (such as electron kinetic energy and inductive storage energy into wire kinetic energy (wired rotor example), or atomic electron kinetic energy to magnetic domain kinetic energy (permanent magnet rotor example)) that, in the low-frequency approximation, pretty much matches the amount of energy transfer between fields and the driving (stator) coils which is due to the transformer induction electric field [math]- \frac{\partial \mathbf{A}}{\partial t}[/math] acting on currents in the coils from the relative motion of the magnetic rotor assembly.

As a result, the only way using conventional physics to explain various devices, such as the Marinov Motor[1], the Distinti Paradox2[2], and, especially, the Marinov Generator[3][4] is to more accurately define the electric field [math]\mathbf{E}[/math], which is part of the full Lorentz Force equation, and then see if it leads to predictions that confirm observations, especially those found figure 9 of Cyril Smith's 2009 paper on the "Marinov Generator".[4] S.H.O. talk 22:09, 5 March 2017 (PST)

Prior content in the "Comment Record" section:

Prior content in the "Background" section:

http://www.sho.wiki/index.php?title=Electromagnetic_Potentials&diff=next&oldid=1162

Novel Force laws proposed by various researchers

James Wesley's proposal

James Wesley proposed adding the "motional induction" on charge q. In SI Units, this can be expressed as:[1]

[math]-q(\mathbf{v}\cdot\nabla)\mathbf{A}[/math]

to the Lorentz force.

The idea behind this was to explain an observation in an experiment involving a "Marinov Motor"[1] in which longitudinal induction forces were produced.

The extra term is equivalent to:[5]

[math]-(\mathbf{v}\cdot\nabla)\mathbf{A} = \begin{matrix} - \left[ v_x \frac{∂A_x}{∂x} + v_y \frac{∂A_x}{∂y} + v_z \frac{∂A_x}{∂z} \right] \mathbf{e}_x \\ - \left[ v_x \frac{∂A_y}{∂x} + v_y \frac{∂A_y}{∂y} + v_z \frac{∂A_y}{∂z} \right] \mathbf{e}_y \\ - \left[ v_x \frac{∂A_z}{∂x} + v_y \frac{∂A_z}{∂y} + v_z \frac{∂A_z}{∂z} \right] \mathbf{e}_z \end{matrix}[/math]

Where [math]\mathbf{v}[/math] is the velocity of the charge.

The Lorentz force is:

[math]\mathbf{F} = q\left[-\nabla \varphi- \frac{\partial \mathbf{A}}{\partial t}+ \nabla_\mathbf{A}(\mathbf{v}\cdot\mathbf{A})-(\mathbf{v}\cdot\nabla)\mathbf{A} \right][/math]

Therefore, adding the extra term proposed by Wesley results in:

[math]\mathbf{F} = q\left[-\nabla \varphi- \frac{\partial \mathbf{A}}{\partial t}+ \nabla_\mathbf{A}(\mathbf{v}\cdot\mathbf{A})-2(\mathbf{v}\cdot\nabla)\mathbf{A} \right][/math]

The problem with this modification:

In the case of an electrical charge approaching a wire, this additional term proposed by Wesley would double the force of deflection. This is not observed.

Consider the vector potential due a current-carrying wire on the x-axis. Both the current and the vector potential of this current point in the [math]+x[/math] direction. Now have a charge approaching this wire perpendicularly. Both the magnetic Lorentz force and the additional Wesley term predict the same force. The proposed additional term is superfluous and would double the force of deflection if added.

Cyril Smith's proposal

Cyril Smith proposed adding the following gradient to the Lorentz force:[3]

[math]- \nabla_\mathbf{A}(\mathbf{v}\cdot\mathbf{A})[/math]

This is equal to:

[math]- \nabla_\mathbf{A}(\mathbf{v} \cdot \mathbf{A}) = \begin{matrix} - \left[ v_x \frac{∂A_x}{∂x} + v_y \frac{∂A_y}{∂x} + v_z \frac{∂A_z}{∂x} \right] \mathbf{e}_x \\ - \left[ v_x \frac{∂A_x}{∂y} + v_y \frac{∂A_y}{∂y} + v_z \frac{∂A_z}{∂y} \right] \mathbf{e}_y \\ - \left[ v_x \frac{∂A_x}{∂z} + v_y \frac{∂A_y}{∂z} + v_z \frac{∂A_z}{∂z} \right] \mathbf{e}_z \end{matrix}[/math]

The idea behind this was to explain an observation in an experiment involving a "Marinov Generator"[3] in which longitudinal induction forces were produced.

The Lorentz force is:

[math]\mathbf{F} = q \left[ -\nabla \varphi - \frac{\partial \mathbf{A}}{\partial t} + \nabla_\mathbf{A}(\mathbf{v}\cdot\mathbf{A}) - (\mathbf{v}\cdot\nabla)\mathbf{A} \right][/math]

Therefore, adding the extra term proposed by Cyril Smith results in:

[math]\mathbf{F} = q \left[ -\nabla \varphi - \frac{\partial \mathbf{A}}{\partial t} + 0\nabla_\mathbf{A}(\mathbf{v}\cdot\mathbf{A}) - (\mathbf{v}\cdot\nabla)\mathbf{A} \right][/math]

The problem with this modification:

In the case of two parallel current-carrying wires, this additional term proposed by Cyril Smith negates the magnetic forces between the currents. It turns out that the extra term may yield forces perpendicular to the velocity. The relevant field components are:

[math]- \nabla_\mathbf{A}(\mathbf{v} \cdot \mathbf{A})_\bot = \begin{matrix} - \left[ v_y \frac{∂A_y}{∂x} + v_z \frac{∂A_z}{∂x} \right] \mathbf{e}_x \\ - \left[ v_x \frac{∂A_x}{∂y} + v_z \frac{∂A_z}{∂y} \right] \mathbf{e}_y \\ - \left[ v_x \frac{∂A_x}{∂z} + v_y \frac{∂A_y}{∂z} \right] \mathbf{e}_z \end{matrix}[/math]

Explaining the Marinov Motor and Cyril Smith's "Marinov Generator" using Conventional Physics

See the online bulletin thread titled "Electromagnetism and relativity" for details on this video[6].

The Liénard–Wiechert electric field was derived from the Liénard–Wiechert potentials [math]\varphi[/math] and [math]A[/math][7][8] by Kirk T. McDonald[9], Professor Emeritus of Princeton University in New Jersey[10].

In the sub-sections below, Gaussian units are used unless otherwise noted. Also, do note that all electric fields due to source charges are, in the following sub-sections, evaluated in the rest frame of each target charge subject to them separately. This procedure has roots in an approach to electromagnetism introduced by Edward M. Purcell[11] in Section 5.6 of the Berkeley Physics Course (Volume II) titled Electricity and Magnetism[12][13] and is explained by a video by Veritasium titled "How Special Relativity Makes Magnets Work"[14]. This avoids having to perform calculations based on the magnetic field viewed by an arbitrary inertial observer. This procedure relies on the relative velocities between the charges. The calculations in the sub-sections below are valid for [math]v \ll c[/math].

The Liénard–Wiechert electric fields for electrically-neutral currents

From the paper titled "Onoochin's Paradox" by Kirk T. McDonald[15][16], we have following statement:

For calculations of the Lorentz force to be accurate to order [math]\frac{1}{c^2}[/math], it suffices to use eq. (4) for the magnetic field. However, to maintain the desired accuracy the electric field of a moving charge must also include effects of retardation, as can be obtained from an expansion of the Liénard–Wiechert fields [7][8] (for details, see the appendix of [9]),
[math]\mathbf{E} \approx q\ \frac{\mathbf{\hat{r}}}{r^2} \left(1 + \frac{v^2}{2c^2} - 3\frac{\left(\mathbf{\hat{r}} \cdot \mathbf{v}\right)^2}{2c^2} \right) - \frac{q}{2c^2r}\left[\mathbf{a}+\left(\mathbf{a} \cdot \mathbf{\hat{r}} \right)\mathbf{\hat{r}}\right][/math]
where [math]\mathbf{a}[/math] is the acceleration [math]\mathbf{a}[/math] of the charge [math]q[/math] at the present time.

Let's consider the situation where the acceleration [math]\mathbf{a}[/math] of charge [math]q[/math] is negligible. The electric field at [math]\mathbf{r}[/math] due to source charge [math]q[/math] located at the origin is:

[math]\mathbf{E} = q\ \frac{\mathbf{\hat{r}}}{r^2} \left(1 + \frac{v^2}{2c^2} - 3\frac{\left(\mathbf{\hat{r}} \cdot \mathbf{v}\right)^2}{2c^2} \right)[/math]

In the Coulomb gauge, the first term in the parentheses comes from the electric scalar potential of a charge at rest in the observer's inertial frame. In event that the charge is contained within an electrically-neutral body, the electric field reduces to:

[math]\mathbf{E} = q\ \frac{\mathbf{\hat{r}}}{r^2} \left(\frac{v^2}{2c^2} - 3\frac{\left(\mathbf{\hat{r}} \cdot \mathbf{v}\right)^2}{2c^2} \right)[/math]

Let's consider a target charge [math]Q[/math] located at [math]\mathbf{r}[/math] at rest in the observer's inertial frame. The inertial observer and the charge [math]q[/math] agree on what the electric field [math]\mathbf{E}[/math] is, they agree that there is no magnetic force on [math]Q[/math], and finally, they agree on the acceleration of [math]Q[/math].

Another way to express this result is in terms of the angle [math]\theta[/math] between [math]\mathbf{v}[/math] and [math]\mathbf{r}[/math]:

[math]\mathbf{E} = q\ \frac{\mathbf{\hat{r}}}{r^2} \frac{v^2}{2c^2} \left(1 - 3\ cos^2\theta \right)[/math]

The above equation can be broken up into two parts, one based on the relative azimuthal velocity [math]\mathbf{v}_{\theta}[/math] of the source, and one based on the relative radial velocity [math]\mathbf{v}_r[/math] of the source. First we rearrange the equation:

[math]\mathbf{E} = q\ \frac{\mathbf{\hat{r}}}{r^2} \left(\frac{v^2 - \left(\mathbf{\hat{r}} \cdot \mathbf{v}\right)^2}{2c^2} - 2\frac{\left(\mathbf{\hat{r}} \cdot \mathbf{v}\right)^2}{2c^2} \right)[/math]

Since the relative azimuthal velocity and the relative radial velocity are orthogonal, we can express the following:

[math]\mathbf{E} = q\ \frac{\mathbf{\hat{r}}}{r^2} \left(\frac{v_{\theta}^2}{2c^2} - \frac{v_r^2}{c^2} \right)[/math]

The relative azimuthal velocity can be split into orthogonal components in [math]x[/math] and [math]y[/math] according to the following Pythagorean relation:

[math]|\mathbf{v}_{\theta}|^2 = |\mathbf{v}_x|^2 + |\mathbf{v}_y|^2[/math]
[math]v_{\theta}^2 = v_x^2 + v_y^2[/math]

Therefore:

[math]\mathbf{E} = q\ \frac{\mathbf{\hat{r}}}{r^2} \left(\frac{v_x^2}{2c^2} + \frac{v_y^2}{2c^2} - \frac{v_r^2}{c^2} \right)[/math]

The variables [math]v_x[/math], [math]v_y[/math], and [math]v_r[/math] are mutually independent from each other. Therefore, the electric field [math]\mathbf{E}[/math] can be split into three co-radial contributions:

[math]\mathbf{E}_x = + q\ \frac{\mathbf{\hat{r}}}{r^2} \frac{v_x^2}{2c^2}[/math]
[math]\mathbf{E}_y = + q\ \frac{\mathbf{\hat{r}}}{r^2} \frac{v_y^2}{2c^2}[/math]
[math]\mathbf{E}_r = - q\ \frac{\mathbf{\hat{r}}}{r^2} \frac{v_r^2}{c^2}[/math]

Consider the existence of four charges:

[math]q_-[/math] : the negative charge of loose electrons of the source current element
[math]Q_-[/math] : the negative charge of loose electrons of the target current element
[math]q_+[/math] : the positive charge of metallic atoms (excluding loose electrons) of the source current element
[math]Q_+[/math] : the positive charge of metallic atoms (excluding loose electrons) of the target current element

Consider the existence of their corresponding current elements:

[math]id\mathbf{l}[/math] : the source current element
[math]Id\mathbf{L}[/math] : the target current element

Where:

[math]i[/math] is the source current and [math]d\mathbf{l}[/math] is its length element.
[math]I[/math] is the target current and [math]d\mathbf{L}[/math] is its length element.

The current elements are equal to:

[math]id\mathbf{l} = q_- \mathbf{v}_d[/math]
[math]Id\mathbf{L} = Q_- \mathbf{V}_d[/math]

Where:

[math]\mathbf{v}_d[/math] is the drift velocity of the electrons of the source current element.
[math]\mathbf{V}_d[/math] is the drift velocity of the electrons of the target current element.

Say we want to calculate the force on the target current element due to the source current element. This requires us to analyze four different forces:

Force [math]\mathbf{F}_{--}[/math] on [math]Q_-[/math] by field [math]\mathbf{E}_{--}[/math] of [math]q_-[/math]
Force [math]\mathbf{F}_{+-}[/math] on [math]Q_-[/math] by field [math]\mathbf{E}_{+-}[/math] of [math]q_+[/math]
Force [math]\mathbf{F}_{-+}[/math] on [math]Q_+[/math] by field [math]\mathbf{E}_{-+}[/math] of [math]q_-[/math]
Force [math]\mathbf{F}_{++}[/math] on [math]Q_+[/math] by field [math]\mathbf{E}_{++}[/math] of [math]q_+[/math]

These four forces are dependent on four different relative velocities (source velocity w.r.t. target velocity):

Relative velocity [math]\mathbf{v}_{--} = \mathbf{v}_- - \mathbf{V}_-[/math]
Relative velocity [math]\mathbf{v}_{+-} = \mathbf{v}_+ - \mathbf{V}_-[/math]
Relative velocity [math]\mathbf{v}_{-+} = \mathbf{v}_- - \mathbf{V}_+[/math]
Relative velocity [math]\mathbf{v}_{++} = \mathbf{v}_+ - \mathbf{V}_+[/math]

The corresponding relative speeds are:

Relative speed [math]v_{x--} = \left(v_- - V_-\right)_x[/math]
Relative speed [math]v_{x+-} = \left(v_+ - V_-\right)_x[/math]
Relative speed [math]v_{x-+} = \left(v_- - V_+\right)_x[/math]
Relative speed [math]v_{x++} = \left(v_+ - V_+\right)_x[/math]
Relative speed [math]v_{y--} = \left(v_- - V_-\right)_y[/math]
Relative speed [math]v_{y+-} = \left(v_+ - V_-\right)_y[/math]
Relative speed [math]v_{y-+} = \left(v_- - V_+\right)_y[/math]
Relative speed [math]v_{y++} = \left(v_+ - V_+\right)_y[/math]
Relative speed [math]v_{r--} = \left(v_- - V_-\right)_r[/math]
Relative speed [math]v_{r+-} = \left(v_+ - V_-\right)_r[/math]
Relative speed [math]v_{r-+} = \left(v_- - V_+\right)_r[/math]
Relative speed [math]v_{r++} = \left(v_+ - V_+\right)_r[/math]

The two drift velocities are:

[math]\mathbf{v}_d = \mathbf{v}_- - \mathbf{v}_+[/math] : The drift velocity of the loose electrons of the source current element
[math]\mathbf{V}_d = \mathbf{V}_- - \mathbf{V}_+[/math] : The drift velocity of the loose electrons of the target current element

Let the effective velocities of the current elements be:

[math]\mathbf{v} = \left( \mathbf{v}_- + \mathbf{v}_+ \right)/2[/math]
[math]\mathbf{V} = \left( \mathbf{V}_- + \mathbf{V}_+ \right)/2[/math]

So the effective relative velocity between the current elements (source velocity w.r.t. target velocity) is:

[math]\mathbf{v}_{rel} = \mathbf{v} - \mathbf{V}[/math]

The effective velocity of each current element is halfway between the velocity of the negative charges and the velocity of the positive charges, so one may rather define a new variable, the deviation velocity, to be one-half the drift velocity of the electrons:

[math]\mathbf{u} = \mathbf{v}_d / 2[/math] is the deviation velocity of the source current element.
[math]\mathbf{U} = \mathbf{V}_d / 2[/math] is the deviation velocity of the target current element.

All four relative velocities can be expressed in terms of the deviation velocities [math]\mathbf{u}[/math] and [math]\mathbf{U}[/math] together with the relative velocity [math]\mathbf{v}_{rel}[/math].

Relative velocity [math]\mathbf{v}_{--} = \mathbf{v}_- - \mathbf{V}_- = \mathbf{u} + \mathbf{v}_{rel} - \mathbf{U}[/math]
Relative velocity [math]\mathbf{v}_{+-} = \mathbf{v}_+ - \mathbf{V}_- = -\mathbf{u} + \mathbf{v}_{rel} - \mathbf{U}[/math]
Relative velocity [math]\mathbf{v}_{-+} = \mathbf{v}_- - \mathbf{V}_+ = \mathbf{u} + \mathbf{v}_{rel} + \mathbf{U}[/math]
Relative velocity [math]\mathbf{v}_{++} = \mathbf{v}_+ - \mathbf{V}_+ = -\mathbf{u} + \mathbf{v}_{rel} + \mathbf{U}[/math]

The corresponding relative speeds are:

Relative speed [math]v_{x--} = \left(v_- - V_- = u + v_{rel} - U\right)_x[/math]
Relative speed [math]v_{x+-} = \left(v_+ - V_- = -u + v_{rel} - U\right)_x[/math]
Relative speed [math]v_{x-+} = \left(v_- - V_+ = u + v_{rel} + U\right)_x[/math]
Relative speed [math]v_{x++} = \left(v_+ - V_+ = -u + v_{rel} + U\right)_x[/math]
Relative speed [math]v_{y--} = \left(v_- - V_- = u + v_{rel} - U\right)_y[/math]
Relative speed [math]v_{y+-} = \left(v_+ - V_- = -u + v_{rel} - U\right)_y[/math]
Relative speed [math]v_{y-+} = \left(v_- - V_+ = u + v_{rel} + U\right)_y[/math]
Relative speed [math]v_{y++} = \left(v_+ - V_+ = -u + v_{rel} + U\right)_y[/math]
Relative speed [math]v_{r--} = \left(v_- - V_- = u + v_{rel} - U\right)_r[/math]
Relative speed [math]v_{r+-} = \left(v_+ - V_- = -u + v_{rel} - U\right)_r[/math]
Relative speed [math]v_{r-+} = \left(v_- - V_+ = u + v_{rel} + U\right)_r[/math]
Relative speed [math]v_{r++} = \left(v_+ - V_+ = -u + v_{rel} + U\right)_r[/math]

The equation for the electric field contains squared values of the speed. As noted before:

[math]\mathbf{E} = q\ \frac{\mathbf{\hat{r}}}{r^2} \left(\frac{v^2 - \left(\mathbf{\hat{r}} \cdot \mathbf{v}\right)^2}{2c^2} - 2\frac{\left(\mathbf{\hat{r}} \cdot \mathbf{v}\right)^2}{2c^2} \right)[/math]

The electric fields on charges [math]Q_+[/math] and [math]Q_-[/math], in their own respective and distinct rest frames, are as follows:

[math]\mathbf{E}_- = \frac{\mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( q_- \left(v_{x--}^2 + v_{y--}^2 - 2 v_{r--}^2\right) + q_+ \left(v_{x+-}^2 + v_{y+-}^2 - 2 v_{r+-}^2\right) \right)[/math]
[math]\mathbf{E}_+ = \frac{\mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( q_- \left(v_{x-+}^2 + v_{y-+}^2 - 2 v_{r-+}^2\right) + q_+ \left(v_{x++}^2 + v_{y++}^2 - 2 v_{r++}^2\right) \right)[/math]

The Forces due to the Liénard–Wiechert electric fields for electrically-neutral currents

It may be more helpful to calculate the forces due to relative radial velocities separately from the forces due to relative azimuthal velocities. Therefore:

[math]\mathbf{F}_x = \frac{\mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( Q_- q_- v_{x--}^2 + Q_- q_+ v_{x+-}^2 + Q_+ q_- v_{x-+}^2 + Q_+ q_+ v_{x++}^2 \right)[/math]
[math]\mathbf{F}_y = \frac{\mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( Q_- q_- v_{y--}^2 + Q_- q_+ v_{y+-}^2 + Q_+ q_- v_{y-+}^2 + Q_+ q_+ v_{y++}^2 \right)[/math]
[math]\mathbf{F}_r = \frac{\mathbf{\hat{r}}}{r^2} \frac{- 2}{2c^2} \left( Q_- q_- v_{r--}^2 + Q_- q_+ v_{r+-}^2 + Q_+ q_- v_{r-+}^2 + Q_+ q_+ v_{r++}^2 \right)[/math]

It follows that [math]\mathbf{F}_x[/math] and [math]\mathbf{F}_y[/math] are functions of currents perpendicular to radial vector [math]\mathbf{r}[/math] while [math]\mathbf{F}_r[/math] is a function of currents co-linear with radial vector [math]\mathbf{r}[/math]

All charges ([math]Q_-[/math], [math]Q_+[/math], [math]q_-[/math], and [math]q_+[/math]) may contribute simultaneously to the azimuthally-directed (transverse) currents in [math]x[/math] and [math]y[/math] and the radially-directed (longitudinal) currents in [math]r[/math].

It will be very advantageous to simplify these formulas. For electrically-neutral currents, we can recognize the following:

[math]q_+ = - q_-[/math]
[math]Q_+ = - Q_-[/math]

This allows us to factor out [math]Q_+ q_+[/math] with the following result:

[math]\mathbf{F}_x = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( v_{x--}^2 - v_{x+-}^2 - v_{x-+}^2 + v_{x++}^2 \right)[/math]
[math]\mathbf{F}_y = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( v_{y--}^2 - v_{y+-}^2 - v_{y-+}^2 + v_{y++}^2 \right)[/math]
[math]\mathbf{F}_r = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{- 2}{2c^2} \left( v_{r--}^2 - v_{r+-}^2 - v_{r-+}^2 + v_{r++}^2 \right)[/math]

Next, we will work on simplifying the contents within the parentheses. To make matters simpler, we will move the subscript to the lower right corner of the parentheses.

[math]\mathbf{F}_x = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 \right)_x[/math]
[math]\mathbf{F}_y = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 \right)_y[/math]
[math]\mathbf{F}_r = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{- 2}{2c^2} \left( v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 \right)_r[/math]

These can be distributed back to the contents within the parentheses after they are substituted for a different expression in terms of the currents. The terms in the parentheses (disregarding their sign) are as follows:

[math]v_{--}^2 = \left(u + v_{rel} - U\right)^2 = u^2 + v_{rel}^2 + U^2 + (2u - 2U) v_{rel} - 2 u U[/math]
[math]v_{+-}^2 = \left(-u + v_{rel} - U\right)^2 = u^2 + v_{rel}^2 + U^2 + (-2u - 2U) v_{rel} + 2 u U[/math]
[math]v_{-+}^2 = \left(u + v_{rel} + U\right)^2 = u^2 + v_{rel}^2 + U^2 + (2u + 2U) v_{rel} + 2 u U[/math]
[math]v_{++}^2 = \left(-u + v_{rel} + U\right)^2 = u^2 + v_{rel}^2 + U^2 + (-2u + 2U) v_{rel} - 2 u U[/math]

Therefore sum of the forces on [math]Q_-[/math] and [math]Q_+[/math] depends on:

[math]v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 = ((2u - 2U) - (-2u - 2U) - (2u + 2U) + (-2u + 2U))v_{rel} + ((-2) - (2) - (2) + (-2)) u U[/math]
[math]v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 = ((2u - 2U) + (2u + 2U) - (2u + 2U) - (2u - 2U))v_{rel} - 8 u U[/math]
[math]v_{--}^2 - v_{+-}^2 - v_{-+}^2 - v_{+-}^2 + v_{++}^2 = ((-2U) + (2U) - (2U) - (- 2U))v_{rel} - 8 u U[/math]
[math]v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 = (0)v_{rel} - 8 u U[/math]
[math]v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 = - 8 u U[/math]
[math]v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 = - 8 (v_d/2)(V_d/2)[/math]
[math]v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 = - 2 (v_d)(V_d)[/math]

The force on [math]Q_-[/math] depends on:

[math]v_{--}^2 - v_{+-}^2 = ((2u - 2U) - (-2u - 2U))v_{rel} + ((-2) - (2)) u U[/math]
[math]v_{--}^2 - v_{+-}^2 = 4 u v_{rel} - 4 u U[/math]

The force on [math]Q_+[/math] depends on:

[math]v_{++}^2 - v_{-+}^2 = ((-2u + 2U) - (2u + 2U))v_{rel} + (- (2) + (-2)) u U[/math]
[math]v_{++}^2 - v_{-+}^2 = -4 u v_{rel} - 4 u U[/math]

Recalling that:

[math]\mathbf{F}_x = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 \right)_x[/math]
[math]\mathbf{F}_y = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 \right)_y[/math]
[math]\mathbf{F}_r = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{- 2}{2c^2} \left( v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 \right)_r[/math]
[math]v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 = - 2 (v_d)(V_d)[/math]

Substitution yields:

[math]\mathbf{F}_x = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( - 2 (v_d)(V_d) \right)_x[/math]
[math]\mathbf{F}_y = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( - 2 (v_d)(V_d) \right)_y[/math]
[math]\mathbf{F}_r = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{- 2}{2c^2} \left( - 2 (v_d)(V_d) \right)_r[/math]

Simplified:

[math]\mathbf{F}_x = - \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( (v_d)(V_d) \right)_x[/math]
[math]\mathbf{F}_y = - \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( (v_d)(V_d) \right)_y[/math]
[math]\mathbf{F}_r = + \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{2}{c^2} \left( (v_d)(V_d) \right)_r[/math]

As stated earlier, the current elements are equal to:

[math]id\mathbf{l} = q_- \mathbf{v}_d[/math]
[math]Id\mathbf{L} = Q_- \mathbf{V}_d[/math]

We can now substitute the currents into the equation. First we substitute [math]Q_- q_-[/math] for [math]Q_+ q_+[/math]:

[math]\mathbf{F}_x = - \frac{Q_- q_- \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( (v_d)(V_d) \right)_x[/math]
[math]\mathbf{F}_y = - \frac{Q_- q_- \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( (v_d)(V_d) \right)_y[/math]
[math]\mathbf{F}_r = + \frac{Q_- q_- \mathbf{\hat{r}}}{r^2} \frac{2}{c^2} \left( (v_d)(V_d) \right)_r[/math]

Next, we assign each charge with their corresponding drift velocities:

[math]\mathbf{F}_x = - \frac{\mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( (Q_- V_d) (q_- v_d) \right)_x[/math]
[math]\mathbf{F}_y = - \frac{\mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( (Q_- V_d) (q_- v_d) \right)_y[/math]
[math]\mathbf{F}_r = + \frac{\mathbf{\hat{r}}}{r^2} \frac{2}{c^2} \left( (Q_- V_d) (q_- v_d) \right)_r[/math]

Next, we make a substitution for the current elements:

[math]\mathbf{F}_x = - \frac{\mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( (Id\mathbf{L}) \cdot (id\mathbf{l}) \right)_x[/math]
[math]\mathbf{F}_y = - \frac{\mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( (Id\mathbf{L}) \cdot (id\mathbf{l}) \right)_y[/math]
[math]\mathbf{F}_r = + \frac{\mathbf{\hat{r}}}{r^2} \frac{2}{c^2} \left( (Id\mathbf{L}) \cdot (id\mathbf{l}) \right)_r[/math]

This can be written as:

[math]\mathbf{F}_x = - \frac{\mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( Ii (d\mathbf{L} \cdot d\mathbf{l}) \right)_x[/math]
[math]\mathbf{F}_y = - \frac{\mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( Ii (d\mathbf{L} \cdot d\mathbf{l}) \right)_y[/math]
[math]\mathbf{F}_r = + \frac{\mathbf{\hat{r}}}{r^2} \frac{2}{c^2} \left( Ii (d\mathbf{L} \cdot d\mathbf{l}) \right)_r[/math]

Or:

[math]\mathbf{F}_x = - \frac{\mathbf{\hat{r}}}{r^2} \frac{Ii}{c^2} \left( d\mathbf{L} \cdot d\mathbf{l} \right)_x[/math]
[math]\mathbf{F}_y = - \frac{\mathbf{\hat{r}}}{r^2} \frac{Ii}{c^2} \left( d\mathbf{L} \cdot d\mathbf{l} \right)_y[/math]
[math]\mathbf{F}_r = + \frac{\mathbf{\hat{r}}}{r^2} \frac{2Ii}{c^2} \left( d\mathbf{L} \cdot d\mathbf{l} \right)_r[/math]

Or:

[math]\mathbf{F}_x = - \frac{\mathbf{\hat{r}}}{r^2} \frac{Ii}{c^2} \left( d\mathbf{L}_x \cdot d\mathbf{l}_x \right)[/math]
[math]\mathbf{F}_y = - \frac{\mathbf{\hat{r}}}{r^2} \frac{Ii}{c^2} \left( d\mathbf{L}_y \cdot d\mathbf{l}_y \right)[/math]
[math]\mathbf{F}_r = + \frac{\mathbf{\hat{r}}}{r^2} \frac{2Ii}{c^2} \left( d\mathbf{L}_r \cdot d\mathbf{l}_r \right)[/math]

Adding the forces together, their sum is:

[math]\mathbf{F} = \frac{\mathbf{\hat{r}}}{r^2} \frac{Ii}{c^2} \left( 2 d\mathbf{L}_r \cdot d\mathbf{l}_r - d\mathbf{L}_x \cdot d\mathbf{l}_x - d\mathbf{L}_y \cdot d\mathbf{l}_y \right)[/math]
[math]\mathbf{F} = \frac{\mathbf{\hat{r}}}{r^2} \frac{Ii}{c^2} \left( 2 d\mathbf{L}_r \cdot d\mathbf{l}_r - d\mathbf{L}_\theta \cdot d\mathbf{l}_\theta \right)[/math]

In S.I. Units:

[math]\mathbf{F} = \frac{\mathbf{\hat{r}}}{4\pi\epsilon_0r^2} \frac{Ii}{c^2} \left( 2 d\mathbf{L}_r \cdot d\mathbf{l}_r - d\mathbf{L}_\theta \cdot d\mathbf{l}_\theta \right)[/math]

Or:

[math]\mathbf{F} = \frac{\mu_0 I i \mathbf{\hat{r}}}{4\pi r^2} \left( 2 d\mathbf{L}_r \cdot d\mathbf{l}_r - d\mathbf{L}_\theta \cdot d\mathbf{l}_\theta \right)[/math]

The second differential is:

[math]d^2\mathbf{F} = \frac{\mu_0 \mathbf{\hat{r}}}{4\pi r^2} \left( 2 \mathbf{I}_r \cdot \mathbf{i}_r - \mathbf{I}_\theta \cdot \mathbf{i}_\theta \right)[/math]

Special Scenario: No relative motion between positive charges

In the case that the positive charges [math]q_+[/math] and [math]Q_+[/math] are essentially stationary, we can simplify the derivation of the field and force equations in the rest frame of the positive charges.

As stated in the parent section, the two drift velocities are:

[math]\mathbf{v}_d = \mathbf{v}_- - \mathbf{v}_+[/math] : The drift velocity of the loose electrons of the source current element
[math]\mathbf{V}_d = \mathbf{V}_- - \mathbf{V}_+[/math] : The drift velocity of the loose electrons of the target current element

These become:

[math]\mathbf{v}_d = \mathbf{v}_-[/math] : The drift velocity of the loose electrons of the source current element
[math]\mathbf{V}_d = \mathbf{V}_-[/math] : The drift velocity of the loose electrons of the target current element

As stated in the parent section, the effective velocities of the current elements are:

[math]\mathbf{v} = \left( \mathbf{v}_- + \mathbf{v}_+ \right)/2[/math]
[math]\mathbf{V} = \left( \mathbf{V}_- + \mathbf{V}_+ \right)/2[/math]

These become:

[math]\mathbf{v} = \mathbf{v}_- /2[/math]
[math]\mathbf{V} = \mathbf{V}_- / 2[/math]

As stated in the parent section, a new variable, the deviation velocity, is defined as one-half the drift velocity of the electrons:

[math]\mathbf{u} = \mathbf{v}_d / 2[/math] is the deviation velocity of the source current element.
[math]\mathbf{U} = \mathbf{V}_d / 2[/math] is the deviation velocity of the target current element.

Therefore:

[math]\mathbf{v} = \mathbf{u} [/math]
[math]\mathbf{V} = \mathbf{U} [/math]

The effective relative velocity between the current elements (source velocity w.r.t. target velocity) is:

[math]\mathbf{v}_{rel} = \mathbf{v} - \mathbf{V}[/math]

Therefore:

[math]\mathbf{v}_{rel} = \mathbf{u} - \mathbf{U}[/math]

As derived in the parent section, all four relative velocities can be expressed in terms of the deviation velocities [math]\mathbf{u}[/math] and [math]\mathbf{U}[/math] together with the relative velocity [math]\mathbf{v}_{rel}[/math].

Relative velocity [math]\mathbf{v}_{--} = \mathbf{v}_- - \mathbf{V}_- = \mathbf{u} + \mathbf{v}_{rel} - \mathbf{U}[/math]
Relative velocity [math]\mathbf{v}_{-+} = \mathbf{v}_- - \mathbf{V}_+ = \mathbf{u} + \mathbf{v}_{rel} + \mathbf{U}[/math]
Relative velocity [math]\mathbf{v}_{+-} = \mathbf{v}_+ - \mathbf{V}_- = -\mathbf{u} + \mathbf{v}_{rel} - \mathbf{U}[/math]
Relative velocity [math]\mathbf{v}_{++} = \mathbf{v}_+ - \mathbf{V}_+ = -\mathbf{u} + \mathbf{v}_{rel} + \mathbf{U}[/math]

By substituting for [math]\mathbf{v}_{rel}[/math], we get:

Relative velocity [math]\mathbf{v}_{--} = \mathbf{v}_- - \mathbf{V}_- = 2\mathbf{u} - 2\mathbf{U}[/math]
Relative velocity [math]\mathbf{v}_{-+} = \mathbf{v}_- - \mathbf{V}_+ = 2\mathbf{u}[/math]
Relative velocity [math]\mathbf{v}_{+-} = \mathbf{v}_+ - \mathbf{V}_- = -2\mathbf{U}[/math]
Relative velocity [math]\mathbf{v}_{++} = \mathbf{v}_+ - \mathbf{V}_+ = 0[/math]

As derived in the parent section, the following is a series sum of terms (a function of relative speeds) which will be used to help calculate the forces between currents:

[math]v_{--}^2 - v_{-+}^2 - v_{+-}^2 + v_{++}^2[/math]

Based on the above results, we have:

[math]v_{--}^2 - v_{-+}^2 - v_{+-}^2 + v_{++}^2 = (2u - 2U)^2 - (2u)^2 - (-2U)^2 + 0[/math]
[math]v_{--}^2 - v_{-+}^2 - v_{+-}^2 + v_{++}^2 = (4u^2 + 4U^2 - 8uU) - 4u^2 - 4U^2 + 0[/math]
[math]v_{--}^2 - v_{-+}^2 - v_{+-}^2 + v_{++}^2 = - 8uU[/math]
[math]v_{--}^2 - v_{-+}^2 - v_{+-}^2 + v_{++}^2 = - 8 (v_{d}/2)(V_{d}/2)[/math]
[math]v_{--}^2 - v_{-+}^2 - v_{+-}^2 + v_{++}^2 = -2 (v_{d})(V_{d})[/math]

This produces the same results as the parent section titled "Explaining the Marinov Motor and Cyril Smith's 'Marinov Generator' using Conventional Physics" and the previous section titled "The Forces due to the Liénard–Wiechert electric fields for electrically-neutral currents":

In S.I. Units:

[math]\mathbf{F} = \frac{\mathbf{\hat{r}}}{4\pi\epsilon_0r^2} \frac{Ii}{c^2} \left( 2 d\mathbf{L}_r \cdot d\mathbf{l}_r - d\mathbf{L}_\theta \cdot d\mathbf{l}_\theta \right)[/math]

Or:

[math]\mathbf{F} = \frac{\mu_0 I i \mathbf{\hat{r}}}{4\pi r^2} \left( 2 d\mathbf{L}_r \cdot d\mathbf{l}_r - d\mathbf{L}_\theta \cdot d\mathbf{l}_\theta \right)[/math]

The second differential is:

[math]d^2\mathbf{F} = \frac{\mu_0 \mathbf{\hat{r}}}{4\pi r^2} \left( 2 \mathbf{I}_r \cdot \mathbf{i}_r - \mathbf{I}_\theta \cdot \mathbf{i}_\theta \right)[/math]

Deriving James Wesley's additional force term in the case of co-linear current elements

As pointed out above in the sub-section titled "James Wesley's proposal", James Wesley proposed adding the "motional induction" on charge [math]q[/math]. In SI Units, this can be expressed as:[1]

[math]-q(\mathbf{v}\cdot\nabla)\mathbf{A}[/math]

How may we derive the same result for current elements located and oriented along the same line? Equation 2b of the paper titled "Observations of the Marinov Motor"[1] can be adapted to the form above, resulting in:

[math]\mathbf{F} = -Q_- \left(\mathbf{V}_d \cdot \nabla \right) \mathbf{A} = Q_- \left(-\mathbf{V}_d \cdot \nabla \right) \mathbf{A} = - I \left( d\mathbf{L} \cdot \nabla \right) \mathbf{A} [/math]

In the case that we are dealing only with co-linear current elements, the result can be expressed as:

[math]\mathbf{F}_r = \left( -Q_- \left(\mathbf{V}_d \cdot \nabla \right) \mathbf{A} = Q_- \left(-\mathbf{V}_d \cdot \nabla \right) \mathbf{A} = - I \left( d\mathbf{L} \cdot \nabla \right) \mathbf{A} \right)_r[/math]

As stated in the parent section, the force on [math]Q_-[/math] depends on:

[math]v_{--}^2 - v_{+-}^2 = ((2u - 2U) - (-2u - 2U))v_{rel} + ((-2) - (2)) u U[/math]
[math]v_{--}^2 - v_{+-}^2 = 4 u v_{rel} - 4 u U[/math]

As stated in the parent section, The force on [math]Q_+[/math] depends on:

[math]v_{++}^2 - v_{-+}^2 = ((-2u + 2U) - (2u + 2U))v_{rel} + (- (2) + (-2)) u U[/math]
[math]v_{++}^2 - v_{-+}^2 = -4 u v_{rel} - 4 u U[/math]

In the case that that there is no relative velocity [math]v_{rel}[/math] between the current elements:

[math]v_{--}^2 - v_{+-}^2 = - 4 u U = -4 (v_d/2)(V_d/2) = - (v_d)(V_d)[/math]
[math]v_{++}^2 - v_{-+}^2 = - 4 u U = -4 (v_d/2)(V_d/2) = - (v_d)(V_d)[/math]

As stated in the section "The Forces due to the Liénard–Wiechert electric fields for electrically-neutral currents", the electric forces on charges [math]Q_-[/math] and [math]Q_+[/math] are:

[math]\mathbf{F}_{r-} = + \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( v_{--}^2 - v_{+-}^2 \right)_r[/math]
[math]\mathbf{F}_{r+} = + \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( v_{++}^2 - v_{-+}^2 \right)_r[/math]

From this and above, it follows that:

[math]\mathbf{F}_{r-} = + \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( - (v_d)(V_d) \right)_r[/math]
[math]\mathbf{F}_{r+} = + \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( - (v_d)(V_d) \right)_r[/math]

Therefore:

[math]\mathbf{F}_{r-} = + \frac{Q_- q_- \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( - (v_d)(V_d) \right)_r[/math]
[math]\mathbf{F}_{r+} = + \frac{Q_+ q_- \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( + (v_d)(V_d) \right)_r[/math]

The effective electric fields experience by each charge is:

[math]\mathbf{E}_{r-} = + \frac{q_- \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( - (v_d)(V_d) \right)_r[/math]
[math]\mathbf{E}_{r+} = + \frac{q_- \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( + (v_d)(V_d) \right)_r[/math]

Given that:

[math]Q_+ = - Q_-[/math]
[math]\mathbf{U} = \mathbf{v}_d / 2[/math]
[math]\mathbf{v}_d = \mathbf{U} - (-\mathbf{U})[/math]
[math]Id\mathbf{L} = Q_- \mathbf{V}_d = Q_- (+ \mathbf{U}) + Q_+ (-\mathbf{U})[/math]

It can be shown that:

[math]\mathbf{F}_{r+} = \left( - Q_+ \left((-\mathbf{U}) \cdot \nabla \right) \mathbf{A} \right)_r[/math]
[math]\mathbf{F}_{r-} = \left( - Q_- \left((+\mathbf{U}) \cdot \nabla \right) \mathbf{A} \right)_r[/math]
[math]\mathbf{F}_r = \left( - Q_- \left(\mathbf{V}_d \cdot \nabla \right) \mathbf{A} \right)_r = \mathbf{F}_{r-} + \mathbf{F}_{r+}[/math]

This is based on the directional derivative of vector potential [math]\mathbf{A}[/math] along [math]Id\mathbf{L}[/math].

It is important to note that the total force [math]\mathbf{F}_r[/math] in [math]r[/math] remains unaffected by the relative motion between the target current element [math]I_rd\mathbf{L}_r[/math] and the source current element [math]i_rd\mathbf{l}_r[/math] responsible for vector potential [math]\mathbf{A}_r[/math]. However, forces [math]\mathbf{F}_{r-}[/math] and [math]\mathbf{F}_{r+}[/math] are affected by the relative speed [math](v_{rel})_r[/math] between the current elements in [math]r[/math], which leads to a kind of electromagnetic induction acting between co-linear current elements in relative motion in [math]r[/math]. This feature may help to explain Cyril Smith's "Marinov Generator"[3].

A New Idea: The Makerarc

The above result supports a development beyond the S.H.O. Drive. A much more powerful and compact permanent magnet system is now envisioned.

The preliminary name is Makerarc (MAY-KERR-ARK), which stands for:

  • Magnetic
  • Atom
  • Kinetic
  • Energy
  • Reservoir
  • And
  • Resource
  • Channel

Details pending. Stay tuned. Sincerely, S.H.O. talk 22:09, 5 March 2017 (PST)

What makes a magnet tick? Clues can be found in the video titled "MAGNETS: How Do They Work?" by Veritasium[17]. S.H.O. talk 22:08, 23 March 2017 (PDT)

An attempt to produce the Makerarc back in March 2017[18] failed to demonstrate the predictions of the term [math]\mathbf{F}_r = \left( - Q_- \left(\mathbf{V}_d \cdot \nabla \right) \mathbf{A} \right)_r[/math][19] where I assumed that an external vector potential [math]\mathbf{A}[/math] would lead to a force co-linear to current elements. This may be due to the supercurrent nature of electrons bounded to atoms, allowing atoms to be capable of preserving or "fixing" the amount of magnetic flux that passes through them, just as "macroscopic" superconducting currents do for superconductors, and therefore, by extension be capable of resisting changes in the magnetic vector potential. This does not negate the possibility of a Marinov Generator, as designed by Cyril Smith, because in his case the currents receiving power were inside a conductor where an externally applied vector potential is not fully shielded against, permitting [math]\mathbf{F}_r = \left( - Q_- \left(\mathbf{V}_d \cdot \nabla \right) \mathbf{A} \right)_r[/math] to yield a force. At the moment, I am still trying to decide on the correct model, but hopefully the surviving concept is similar to what I have developed so far in this page. Sincerely, S.H.O. talk 19:55, 14 July 2017 (PDT)

Explaining "Altered" Lenz' Law Devices

It is anticipated that the longitudinal force described above in the section titled "Explaining the Marinov Motor and Cyril Smith's 'Marinov Generator' using Conventional Physics" may explain some types of purported Reduced-Lenz devices. A good example can be found in the video below ("Ray's No Back EMF Generator"), although in this example, the longitudinal force increases the drag, mainly in positions where pancake generator coil is mostly outside the cylindrical boundary of the permanent magnet. This creates an illusion of a "Reduced" Lenz' Law effect when the magnet is mostly within the cylindrical boundary of the permanent magnet:

Simulations in JavaScript and THREE.js have determined that in many other configurations of currents and magnets, the magnetic Lorentz forces [math]q\ \mathbf{v} \times \mathbf{B}[/math] will be opposed in part by the additional force. This makes certain types of magnetic circuit arrangements more efficient at electrical power generation (in terms of output power vs. input power) but less efficient as an electrical motor. Simulations of the S.H.O. Drive indicated that the additional force may vary between assisting or opposing the magnetic Lorentz force, which is helpful depending on whether the goal is to develop an efficient motor or an efficient generator, respectively. The simulated S.H.O. Drive designs appeared to switch between these two extremes every quarter cycle. The only way to make it work would have been to increase the inductive reactance (so that the current would be delayed by about a quarter cycle), and even then, simulations indicated that the extra term was typically insufficient to completely reverse the net force. However simulations of various rotor and S.H.O. coil curve modifications for the S.H.O. Drive showed that it was possible for the additional force to be a significant percentage of the magnetic Lorentz force. Per more recent simulations (early March), the Makerarc design (previous section) will improve upon this many fold. S.H.O. talk 23:39, 5 March 2017 (PST)

Explaining the Newman Motor

See also: Memory Lane

A Newman Motor-style coil and magnet arrangement, like that shown in the video below, have been simulated by me using JavaScript and THREE.js.

The extra electric field term predicts a significant opposition to the magnetic Lorentz force at angles slightly straying from the "top-dead-vertical" position, making it a better generator than a motor. However, when energy is discharged from the "generator coil" to the "motor coil" of Newman's motor, the rotor will have often changed position to the point where the magnetic Lorentz force becomes increasingly significant, helpful for motive purposes. Newman's motor operated at a high Q, which facilitated energy recovery. S.H.O. talk 00:00, 6 March 2017 (PST)

References

  1. 1.0 1.1 1.2 1.3 1.4 http://redshift.vif.com/JournalFiles/Pre2001/V05NO3PDF/v05n3phi.pdf
  2. http://www.distinti.com/docs/pdx/paradox2.pdf
  3. 3.0 3.1 3.2 3.3 http://overunity.com/14691/the-marinov-generator/
  4. 4.0 4.1 http://www.overunityresearch.com/index.php?action=dlattach;topic=2470.0;attach=13897
  5. http://www.overunityresearch.com/index.php?action=dlattach;topic=2470.0;attach=13908
  6. https://www.physicsforums.com/threads/electromagnetism-and-relativity.747811/
  7. 7.0 7.1 https://docs.google.com/file/d/0B817m31MAj0wZTZjZmMwMjgtY2Y5YS00YTQ5LThjM2EtNzhjYTYzNzFlZDY0/edit?hl=en_GB&pli=1
  8. 8.0 8.1 https://docs.google.com/file/d/0B817m31MAj0wMDI1YjllYjctY2NhOS00M2M2LWFlMTUtYjVmYTkyZmVlY2M2/edit?hl=en_GB
  9. 9.0 9.1 http://web.archive.org/web/20170318210550/http://puhep1.princeton.edu/~mcdonald/examples/ph501/ph501lecture24.pdf
  10. https://dof.princeton.edu/about/clerk-faculty/emeritus/kirk-t-mcdonald
  11. http://physics.weber.edu/schroeder/mrr/mrrtalk.html
  12. https://www.scribd.com/doc/128728926/Electricity-and-Magnetism-Berkeley-Physics-Course-Purcell
  13. https://en.wikipedia.org/wiki/Berkeley_Physics_Course
  14. https://www.youtube.com/watch?v=1TKSfAkWWN0
  15. http://citeseerx.ist.psu.edu/viewdoc/summary?doi=10.1.1.299.8534
  16. http://freeweb.siol.net/markoor/onoochin.pdf
  17. https://www.youtube.com/watch?v=hFAOXdXZ5TM
  18. https://www.facebook.com/Sho.Drives/videos/790278744462849/
  19. https://www.facebook.com/Sho.Drives/videos/791514291005961/

See also

Site map

HQGlossaryApril 2016 Presentation