Difference between revisions of "Electromagnetic Potentials"

From S.H.O.
Jump to: navigation, search
(A New Idea: The Makerarc)
 
(149 intermediate revisions by the same user not shown)
Line 3: Line 3:
 
|titlemode=replace
 
|titlemode=replace
 
|keywords=electricity,magnetism,motor,generator
 
|keywords=electricity,magnetism,motor,generator
|description=The basic idea here is that the electromagnetic potentials <math>\phi</math> and <math>A</math> and their derivatives can be used to derive all electromagnetism.
+
|description=The basic idea here is that the electromagnetic potentials φ and A may be the underlying key to several inventions related to electromagnetic forces.
 
}}
 
}}
 
{{DISPLAYTITLE:Function Conjunction → Electromagnetic Potentials}}
 
{{DISPLAYTITLE:Function Conjunction → Electromagnetic Potentials}}
The basic idea here is that the electromagnetic potentials <math>\phi</math> and <math>A</math> and their derivatives can be used to derive all electromagnetism.
+
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.
  
==Comment Record==
+
==Introduction==
  
Beginning with the velocity-dependent electromagnetic potential, one normally may derive the Lorentz force from it. This occurs by substituting a gradient term via an approximate, rather than exact, vector identity <math>\nabla (\mathbf{A} \cdot \mathbf{v}) = \nabla_\mathbf{A}(\mathbf{A} \cdot \mathbf{v})</math>. This is "allowed" normally because velocities in Special Relativity are not an explicit function of the coordinates, a matter simply assumed to be fact.<ref>http://www.nhn.ou.edu/~gut/notes/cm/lect_09.pdf ("The triple cross product can be written in a more compact form where we use the fact that the velocity is not an explicit function of the coordinates.")</ref> However, by taking velocity to be an explicit function of the coordinates, as per the S.R.-like Lorentz Ether theory, the extra force term not seen in the Lorentz force appears, which is <math>\nabla_\mathbf{v}(\mathbf{A} \cdot \mathbf{v})</math>. ''[[User:S.H.O.|S.H.O.]] <sup>[[User_talk:S.H.O.|talk]]</sup>'' 03:18, 12 September 2016 (PDT)
+
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:
  
So herein, the actual vector identity for the gradient of a dot is employed, resulting in mathematical consistency, as opposed to "magically" waving away the velocity-gradient terms as is usually done to impose consistency of the Electromagnetic Lagrangian with the Lorentz Force. These extra terms are gauge-dependent, and so an appropriate gauge must be selected (by Nature itself) to render these (heretical) gauge-dependent forces meaningful. Applying the Lorenz gauge would make it consistent with the finite speed of light, while applying the Coulomb gauge would imply dependence of the force on the instaneous position of the sources of electromagnetic potential. ''[[User:S.H.O.|S.H.O.]] <sup>[[User_talk:S.H.O.|talk]]</sup>'' 13:37, 28 August 2016 (PDT)
+
: <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>
  
The extra force term <math>\nabla_\mathbf{v}(\mathbf{A} \cdot \mathbf{v})</math> consists of a changing effective mass-correction term (in units of kg/s) multiplied by a velocity. Such a force is not preserved under Galilean transformations, but neither is the derivative <math>\frac{d(m\mathbf{v})}{dt} = m \frac{d\mathbf{v}}{dt} + \mathbf{v} \frac{dm}{dt}</math>. The extra force term is essentially a <math>\mathbf{v} \frac{dm}{dt}</math> term and could conceivably substitute for it, unless other terms may also assume a similar role (more on this below). Such forces change the time-like component of the relativistic 4-momentum, and therefore they are related to changes of rest energy of a particle (charge) subject to potentials, as viewed by an arbitrary inertial observer. Whatever other such terms may be, after cancelling the terms on both sides, what remains on the left hand side is equal to <math>m \frac{d\mathbf{v}}{dt}</math> of the particle (charge) and the right hand side yields the standard Lorentz force on the particle (charge). Unlike the 4-momentum, which is Lorentz invariant, the 3-momentum (i.e. the set of the 3 space-like components of the 4-momentum) is not. The Lorentz force, as normally expressed, only deals with changes of the 3-momentum over time. An observer may observe the time-like component of its 4-momentum transform as it changes velocity resulting in an unaccounted for "thrust", particularly if the time-like component of the momentum <math>E/c</math> in some way depends on the ''velocity-dependent electromagnetic potentials'', which in turn would yield, in effect, a ''velocity-dependent rest mass''.
+
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.
  
To allude this possibility, consider that in cases where there are changing charge densities due to divergent/convergent electrical currents, and yet where current densities are constant, the kinetic forces <math>m \frac{d\mathbf{v}}{dt}</math> between charges do not sum to zero. A key such example of another non-Galilean invariant terms involving the variation mass over time at a given velocity, is the hypothetical Longitudinal force density conceived by Koen J. van Vlaenderen in his 2015 paper "General Classical Electrodynamics"<ref>http://vixra.org/abs/1512.0297</ref>, which from equations (2.8) and (2.9) can be written as <math>\mathbf{f}_L(\mathbf{x}) = \mathbf{J} (- \nabla \cdot \mathbf{A})</math>. The extra force term, adapted to help generalize Whittaker's force law to cover field force densities, was intended to preserve Newton's Third Law of Motion ("For every action there is an equal and opposite reaction"). Under the Lorenz gauge condition, this equals <math>\mathbf{f}_L(\mathbf{x}) = \frac{dq}{dV} \mathbf{v} \left( \epsilon_0 \mu_0 \frac{∂\varphi}{∂t} \right) = \frac{d(E_{coulomb})}{c^2dtdV} \mathbf{v}</math>, which when integrated over volume elements <math>dV</math> gives <math>\frac{d(m_{coulomb})}{dt} \mathbf{v} = \frac{d(E_{coulomb})}{c^2dt} \mathbf{v}</math>. When this <math>\mathbf{v} \frac{dm}{dt}</math> type term is subtracted from both sides of the full force equation (a <math>\frac{d(m \mathbf{v})}{dt}</math> type equation), we once again return to the Lorentz 3-force, which is the standard electromagnetic force of type <math>\mathbf{F} = m \frac{d\mathbf{v}}{dt} = m\mathbf{a}</math>. ''[[User:S.H.O.|S.H.O.]] <sup>[[User_talk:S.H.O.|talk]]</sup>'' 03:18, 12 September 2016 (PDT)
+
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)
  
==Background==
+
'''Prior content in the "Comment Record" section:'''
  
According Emil John Konopinski, a nuclear scientist<ref name="nyt obit">[http://query.nytimes.com/gst/fullpage.html?res=9C0CE6D81739F93BA15756C0A966958260  Emil Konopinski, 78, Atomic Bomb Scientist], New York Times</ref> and professor of Mathematics<ref>[http://www.genealogy.ams.org/id.php?id=104307 Eugene Greuling] at the Mathematics Genealogy Project</ref> who worked on the Manhattan Project<ref>{{cite journal|title=Obituary: Emil J. Konopinski|journal=Physics Today|date=October 1991|volume=44|issue=10|pages=144|url=http://www.physicstoday.org/resource/1/phtoad/v44/i10/p144_s1?bypassSSO=1|doi=10.1063/1.2810306|bibcode = 1991PhT....44j.144E }}</ref>, the electromagnetic fields <math>\mathbf{E}</math> and <math>\mathbf{B}</math> can be re-expressed in terms of the electromagnetic potentials <math>\varphi</math> and <math>A</math> through substitutions. From his article on "What the electromagnetic vector potential describes"<ref>http://exvacuo.free.fr/div/Sciences/Dossiers/EM/ScalarEM/J%20Konopinski%20-%20What%20the%20Electromagnetic%20Vector%20Potential%20Describes%20-%20ajp_46_499_78.pdf</ref>, he presents the equation of motion for a localized point charge:
+
* http://www.sho.wiki/index.php?title=Electromagnetic_Potentials&diff=1036&oldid=951
 +
* http://www.sho.wiki/index.php?title=Electromagnetic_Potentials&diff=next&oldid=1165
  
: <math>d(M\mathbf{v})/dt = q \left[ \mathbf{E} + \mathbf{v} \times \mathbf{B}/c \right]\ (Gaussian\ units)</math>:
+
'''Prior content in the "Background" section:'''
: <math>d(M\mathbf{v})/dt = q \left[ \mathbf{E} + \mathbf{v} \times \mathbf{B} \right]\ (SI\ units)</math>
+
  
With the standard substitutions for the fields in terms of the potentials, which were taken to be:
+
http://www.sho.wiki/index.php?title=Electromagnetic_Potentials&diff=next&oldid=1162
  
: <math>\mathbf{E} = - \nabla \varphi - \frac{∂\mathbf{A}}{c∂t}\ (Gaussian\ units)</math>
+
==Novel Force laws proposed by various researchers==
: <math>\mathbf{E} = - \nabla \varphi - \frac{∂\mathbf{A}}{∂t}\ (SI\ units)</math>
+
: <math>\mathbf{B} = \nabla \times \mathbf{A}\ </math>
+
  
Konopinski determined the function relating the time derivative of the "total" momentum with the velocity-dependent potential, as evaluated by an inertial observer who sees the localized charge <math>q</math> of mass <math>m</math> travelling with velocity <math>\mathbf{v}</math> which is subject to a magnetic vector potential <math>\mathbf{A}</math> and an electric scalar potential <math>\varphi</math>.
+
===James Wesley's proposal===
  
: <math>\frac{d}{dt} \left[ M\mathbf{v} + (q/c)\mathbf{A} \right] = - q \left[ \nabla \varphi - (\mathbf{v}/c) \cdot \mathbf{A} \right]\ (Gaussian\ units)</math>
+
James Wesley proposed adding the "motional induction" on charge q. In SI Units, this can be expressed as:<ref name="Marinov Motor"/>
: <math>\frac{d}{dt} \left[ M\mathbf{v} + q\mathbf{A} \right] = - q \left[ \nabla \varphi - \mathbf{v} \cdot \mathbf{A} \right]\ (SI\ units)</math>
+
  
The terms on the right can be separated as follows:
+
: <math>-q(\mathbf{v}\cdot\nabla)\mathbf{A}</math>
  
: <math>\frac{d}{dt} \left[ M\mathbf{v} + q\mathbf{A} \right] = - q \nabla \varphi + q \nabla (\mathbf{v} \cdot \mathbf{A})\ (SI\ units)</math>
+
to the Lorentz force.
  
Using Feynman subscript notation, we can separate the last term on the right into two separate terms:
+
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.
  
: <math>q\nabla (\mathbf{v} \cdot \mathbf{A}) = q\nabla_\mathbf{v}(\mathbf{v} \cdot \mathbf{A}) + q\nabla_\mathbf{A}(\mathbf{v} \cdot \mathbf{A})</math>
+
The extra term is equivalent to:<ref>http://www.overunityresearch.com/index.php?action=dlattach;topic=2470.0;attach=13908</ref>
  
In an article titled "A Discussion on the Magnetic Vector Potential"<ref>http://www.overunityresearch.com/index.php?action=dlattach;topic=2470.0;attach=13908</ref>, Cyril W. Smith (Professor Ph.D of Electronic and Electrical Engineering from 1964-1989<ref>http://www.positivehealth.com/author/cyril-smith-ph-d</ref>), the last term (without the charge <math>q</math>) can be expressed as:
+
: <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 \\
: <math>\nabla_\mathbf{A}(\mathbf{v} \cdot \mathbf{A}) = \begin{matrix}
+
- \left[ v_x \frac{∂A_y}{∂x} + v_y \frac{∂A_y}{∂y} + v_z \frac{∂A_y}{∂z} \right] \mathbf{e}_y \\
\mathbf{a}_x \left[ v_x \frac{∂A_x}{∂x} + v_y \frac{∂A_y}{∂x} + v_z \frac{∂A_z}{∂x} \right] \\
+
- \left[ v_x \frac{∂A_z}{∂x} + v_y \frac{∂A_z}{∂y} + v_z \frac{∂A_z}{∂z} \right] \mathbf{e}_z
\mathbf{a}_y \left[ v_x \frac{∂A_x}{∂y} + v_y \frac{∂A_y}{∂y} + v_z \frac{∂A_z}{∂y} \right] \\
+
\mathbf{a}_z \left[ v_x \frac{∂A_x}{∂z} + v_y \frac{∂A_y}{∂z} + v_z \frac{∂A_z}{∂z} \right]  
+
 
\end{matrix}</math>
 
\end{matrix}</math>
  
Written in this form, <math>\mathbf{a}_x</math>, <math>\mathbf{a}_y</math>, and <math>\mathbf{a}_z</math> are the unit basis vectors for x, y, and z, respectively. However, for the purposes of the S.H.O. Drive Wiki Site, <math>\mathbf{a}</math> will stand for the vector for acceleration, so it doesn't hurt that the standard variables for the basis vectors are really <math>\mathbf{e}_x</math>, <math>\mathbf{e}_y</math>, and <math>\mathbf{e}_z</math>.
+
Where <math>\mathbf{v}</math> is the velocity of the charge.
  
Predictably, the values for first term on the right (without the charge <math>q</math>) can be expressed as:
+
The Lorentz force is:
  
: <math>\nabla_\mathbf{v}(\mathbf{v} \cdot \mathbf{A}) = \begin{matrix}
+
<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>
\left[ A_x \frac{∂v_x}{∂x} + A_y \frac{∂v_y}{∂x} + A_z \frac{∂v_z}{∂x} \right] \mathbf{e}_x \\
+
 
\left[ A_x \frac{∂v_x}{∂y} + A_y \frac{∂v_y}{∂y} + A_z \frac{∂v_z}{∂y} \right] \mathbf{e}_y \\
+
Therefore, adding the extra term proposed by Wesley results in:
\left[ A_x \frac{∂v_x}{∂z} + A_y \frac{∂v_y}{∂z} + A_z \frac{∂v_z}{∂z} \right] \mathbf{e}_z
+
 
\end{matrix}</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>
 +
 
 +
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:<ref name="Marinov Generator"/>
 +
 
 +
<math>- \nabla_\mathbf{A}(\mathbf{v}\cdot\mathbf{A})</math>
  
The sum of the first and last terms gives:
+
This is equal to:
  
:<math>\nabla (\mathbf{v} \cdot \mathbf{A}) = \begin{matrix}
+
: <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} + A_x \frac{∂v_x}{∂x} + A_y \frac{∂v_y}{∂x} + A_z \frac{∂v_z}{∂x} \right] \mathbf{e}_x \\
+
- \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} + A_x \frac{∂v_x}{∂y} + A_y \frac{∂v_y}{∂y} + A_z \frac{∂v_z}{∂y} \right] \mathbf{e}_y \\
+
- \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} + A_x \frac{∂v_x}{∂z} + A_y \frac{∂v_y}{∂z} + A_z \frac{∂v_z}{∂z} \right] \mathbf{e}_z
+
- \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>
 
\end{matrix}</math>
  
By choosing the <math>x</math> Cartesian axis so that it is aligned with the velocity <math>\mathbf{v}</math> of charge <math>q</math> such that <math>\mathbf{v}/|\mathbf{v}| = \mathbf{\hat{v}} = \mathbf{e}_x</math>, the above reduces to:
+
The idea behind this was to explain an observation in an experiment involving a "Marinov Generator"<ref name="Marinov Generator"/> in which longitudinal induction forces were produced.
  
:<math>\nabla (\mathbf{v} \cdot \mathbf{A}) = \begin{matrix}
+
The Lorentz force is:
\left[ v_x \frac{∂A_x}{∂x} + A_x \frac{∂v_x}{∂x} + A_y \frac{∂v_y}{∂x} + A_z \frac{∂v_z}{∂x} \right] \mathbf{e}_x \\
+
 
\left[ v_x \frac{∂A_x}{∂y} \right] \mathbf{e}_y \\
+
<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>
\left[ v_x \frac{∂A_x}{∂z} \right] \mathbf{e}_z
+
 
 +
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>
 
\end{matrix}</math>
  
<math>\nabla_\mathbf{v}(\mathbf{v} \cdot \mathbf{A})</math> simplifies to:
+
==Explaining the Marinov Motor and Cyril Smith's "Marinov Generator" using Conventional Physics==
  
:<math>\left[ \mathbf{A} \cdot \mathbf{a}/|\mathbf{v}| \right] \mathbf{\hat{v}} = \left[ A_x \frac{∂v_x}{∂x} + A_y \frac{∂v_y}{∂x} + A_z \frac{∂v_z}{∂x} \right] \mathbf{e}_x</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 above may also be expressed as:
+
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>.
  
:<math>\nabla_\mathbf{v}(\mathbf{A} \cdot \mathbf{v}) = (\mathbf{A} \cdot \mathbf{a})\mathbf{v} / |\mathbf{v}|^2</math>
+
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>.
  
==Draft==
+
===The Liénard–Wiechert electric fields for electrically-neutral currents===
  
The field experienced by a charge <math>q</math> viewed at rest in a static electromagnetic field 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>\mathbf{F}_{rest,static}/q = - \nabla \varphi</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>),
  
The field experienced by a charge <math>q</math> viewed at rest in a dynamic electromagnetic field is:
+
:<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{F}_{rest,dynamic}/q = - \nabla \varphi - ∂\mathbf{A}/∂t</math>
+
''where <math>\mathbf{a}</math> is the acceleration <math>\mathbf{a}</math> of the charge <math>q</math> at the present time.</blockquote>
  
The field experienced by a moving charge <math>q</math> in a dynamic electromagnetic field (ignoring dilation of proper time relative to coordinate time) is:
+
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}_{moving,dynamic}/q = - \nabla \varphi_q - \mathbf{A}_q/∂t</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>
  
Where:
+
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>\varphi_q = \varphi - \mathbf{v} \cdot A</math> is the scalar potential experienced by the moving charge.<ref>http://exvacuo.free.fr/div/Sciences/Dossiers/EM/ScalarEM/J%20Konopinski%20-%20What%20the%20Electromagnetic%20Vector%20Potential%20Describes%20-%20ajp_46_499_78.pdf</ref>
+
:<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>
  
* <math>∂ \mathbf{A}_q/∂t = ∂\mathbf{A}/∂t + (\mathbf{v} \cdot \nabla)\mathbf{A}</math> is the partial time derivative of the magnetic vector potential experienced by the moving charge.
+
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>.
  
Substituting per the above, the field experienced by the moving charge <math>q</math> is:
+
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{F}/q = - \nabla (\varphi-\mathbf{v} \cdot \mathbf{A}) - ∂\mathbf{A}/∂t - (\mathbf{v} \cdot \nabla)\mathbf{A}</math>
+
:<math>\mathbf{E} = q\ \frac{\mathbf{\hat{r}}}{r^2} \frac{v^2}{2c^2} \left(1 - 3\ cos^2\theta \right)</math>
: <math>\mathbf{F}/q = - \nabla \varphi + \nabla (\mathbf{v} \cdot \mathbf{A}) - ∂\mathbf{A}/∂t - (\mathbf{v} \cdot \nabla)\mathbf{A}</math>
+
  
Using Feynman subscript notation:
+
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>\nabla (\mathbf{v} \cdot \mathbf{A}) = \nabla_\mathbf{v}(\mathbf{v} \cdot \mathbf{A}) + \nabla_\mathbf{A}(\mathbf{v} \cdot \mathbf{A})</math>
+
:<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>
: <math>\nabla_\mathbf{A}(\mathbf{v} \cdot \mathbf{A}) = \mathbf{v} \times \nabla \times \mathbf{A} + (\mathbf{v} \cdot \nabla)\mathbf{A}</math>
+
: <math>\nabla_\mathbf{v}(\mathbf{v} \cdot \mathbf{A}) = \mathbf{A} \times \nabla \times \mathbf{v} + (\mathbf{A} \cdot \nabla)\mathbf{v}</math>
+
  
Substuting for the curl of the vector potential and the curl of the immediate velocity field for a moving charge, we have:
+
Since the '''relative ''azimuthal'' velocity''' and the '''relative ''radial'' velocity''' are orthogonal, we can express the following:
  
: <math>\nabla_\mathbf{A}(\mathbf{v} \cdot \mathbf{A}) = \mathbf{v} \times \mathbf{B} + (\mathbf{v} \cdot \nabla)\mathbf{A}</math>
+
:<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>
: <math>\nabla_\mathbf{v}(\mathbf{v} \cdot \mathbf{A}) = \mathbf{A} \times \mathbf{ω}_\mathbf{v} + (\mathbf{A} \cdot \nabla)\mathbf{v}</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:
 
Where:
  
: <math>\mathbf{B} = \nabla \times \mathbf{A}</math> is the magnetic field.
+
:<math>i</math> is the source current and <math>d\mathbf{l}</math> is its length element.
: <math>\mathbf{ω}_\mathbf{v} = \nabla \times \mathbf{v}</math> is the angular rate of deflection.
+
:<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>
  
Substituting per the above, the field experienced by the moving charge is:
+
In the case that that there is no relative velocity <math>v_{rel}</math> between the current elements:
  
: <math>\mathbf{F}/q = - \nabla \varphi - ∂\mathbf{A}/∂t + \mathbf{v} \times \mathbf{B} + (\mathbf{v} \cdot \nabla)\mathbf{A} + \mathbf{A} \times \mathbf{ω}_\mathbf{v} + (\mathbf{A} \cdot \nabla)\mathbf{v} - (\mathbf{v} \cdot \nabla)\mathbf{A}</math>
+
:<math>v_{--}^2 - v_{+-}^2 = - 4 u U = -4 (v_d/2)(V_d/2) = - (v_d)(V_d)</math>
: <math>\mathbf{F}/q = - \nabla \varphi - ∂\mathbf{A}/∂t + \mathbf{v} \times \mathbf{B} + \mathbf{A} \times \mathbf{ω}_\mathbf{v} + (\mathbf{A} \cdot \nabla)\mathbf{v}</math>
+
:<math>v_{++}^2 - v_{-+}^2 = - 4 u U = -4 (v_d/2)(V_d/2) = - (v_d)(V_d)</math>
  
This field includes the field from Lorentz plus two additional terms:
+
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}/q = \mathbf{F}_{Lorentz}/q + \mathbf{A} \times \mathbf{ω}_\mathbf{v} + (\mathbf{A} \cdot \nabla)\mathbf{v}</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>
: <math>(\mathbf{A} \cdot \nabla)\mathbf{v}</math> is the dot product of the magnetic vector potential with the gradient of the velocity field.
+
:<math>\mathbf{F}_{r+} = + \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( v_{++}^2 - v_{-+}^2  \right)_r</math>
  
For a velocity field defined in the immediate neighborhood of a moving charge <math>q</math> at point <math>p</math>, where the local <math>(\nabla_\mathbf{v} \mathbf{A})_p</math> is a tangent vector on <math>\mathbf{A}</math> (the Lie derivative of <math>\mathbf{v}</math> along <math>\mathbf{A}</math>), the above is equivalent to:
+
From this and above, it follows that:
  
: <math>(\mathbf{A} \cdot \nabla)\mathbf{v} = |\mathbf{A}|(\nabla_\mathbf{v} \mathbf{A})_p = (\mathbf{A} \cdot \mathbf{v})\mathbf{a}/|\mathbf{v}|^2</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>
 +
:<math>\mathbf{F}_{r+} = + \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( - (v_d)(V_d)  \right)_r</math>
  
Where <math>\mathbf{a}</math> is the convective acceleration of the charge, which equals:
+
Therefore:
  
: <math>\mathbf{a} = (∂\mathbf{v}/∂x)|(\mathbf{x}/∂t)|</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>
 +
:<math>\mathbf{F}_{r+} = + \frac{Q_+ q_- \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( + (v_d)(V_d) \right)_r</math>
  
If the charge is taken as a point particle, the convective acceleration is the same as the acceleration.
+
The effective electric fields experience by each charge is:
  
: <math>\mathbf{a} = ∂²\mathbf{x}/∂t²</math>
+
:<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{A} \times \mathbf{ω}_\mathbf{v}</math> is the cross product of the magnetic vector potential and the angular rate of deflection.
+
:<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{ω}_\mathbf{v} = (\mathbf{v} \times \mathbf{a})/|\mathbf{v}|^2</math>
+
  
When fields are static, the field experienced by a moving charge is:
+
Given that:
  
: <math>\mathbf{F}_{moving,static}/q = - \nabla \varphi + \mathbf{v} \times \mathbf{B} + \mathbf{A} \times \mathbf{ω}_\mathbf{v} + (\mathbf{A} \cdot \nabla)\mathbf{v}</math>
+
:<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>
  
So for the case of static fields, the field on an accelerating charge is:
+
It can be shown that:
  
: <math>\mathbf{F}_{moving,static}/q = - \nabla \varphi + \mathbf{v} \times \mathbf{B} + \mathbf{A} \times (\mathbf{v} \times \mathbf{a})/|\mathbf{v}|^2 + (\mathbf{A} \cdot \mathbf{v})\mathbf{a}/|\mathbf{v}|^2</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{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>
  
While the power on an accelerating charge q subject to a static field is:
+
This is based on the directional derivative of vector potential <math>\mathbf{A}</math> along <math>Id\mathbf{L}</math>.
  
: <math>P_{moving,static} = q \left[\ -\nabla \varphi \cdot \mathbf{v} + (\mathbf{v} \times B) \cdot \mathbf{v} + \mathbf{A} \times (\mathbf{v} \times \mathbf{a}) \cdot \mathbf{v}/|\mathbf{v}|^2 + (\mathbf{A} \cdot \mathbf{v})(\mathbf{a} \cdot \mathbf{v})/|\mathbf{v}|^2\ \right]</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"/>.
: <math>P_{moving,static} = q \left[\ -\nabla \varphi \cdot \mathbf{v} + \mathbf{A} \times (\hat{\mathbf{v}} \times \mathbf{a}) \cdot \hat{\mathbf{v}} + (\mathbf{A} \cdot \mathbf{\hat{v}})(\mathbf{a} \cdot \mathbf{\hat{v}})\ \right]</math>
+
  
The field on a moving charge in a changing electromagnetic field becomes:
+
==A New Idea: The Makerarc==
  
: <math>\mathbf{F}_{moving,dynamic}/q = - \nabla \varphi - ∂\mathbf{A}/∂t + \mathbf{v} \times \mathbf{B} + \mathbf{A} \times (\mathbf{v} \times \mathbf{a})/|\mathbf{v}|^2 + (\mathbf{A} \cdot \mathbf{v})\mathbf{a}/|\mathbf{v}|^2</math>
+
The above result supports a development beyond the S.H.O. Drive. A much more powerful and compact permanent magnet system is now envisioned.
  
Summarizing the derivation of the last two terms above, we have:
+
The preliminary name is Makerarc (MAY-KERR-ARK), which stands for:
  
: <math>\nabla_\mathbf{v}(\mathbf{v} \cdot \mathbf{A}) = \mathbf{A} \times \nabla \times \mathbf{v} + (\mathbf{A} \cdot \nabla)\mathbf{v}</math>
+
* '''M'''agnetic
: <math>\mathbf{ω}_\mathbf{v} = \nabla \times \mathbf{v}</math> is the angular rate of deflection.
+
* '''A'''tom
: <math>\nabla_\mathbf{v}(\mathbf{v} \cdot \mathbf{A}) = \mathbf{A} \times \mathbf{ω}_\mathbf{v} + (\mathbf{A} \cdot \nabla)\mathbf{v}</math>
+
* '''K'''inetic
: <math>\mathbf{A} \times \mathbf{ω}_\mathbf{v} = \mathbf{A} \times (\mathbf{v} \times \mathbf{a})/|\mathbf{v}|^2</math>
+
* '''E'''nergy
: <math>(\mathbf{A} \cdot \nabla)\mathbf{v} = (\mathbf{A} \cdot \mathbf{v})\mathbf{a}/|\mathbf{v}|^2</math>
+
* '''R'''eservoir
: <math>\nabla_\mathbf{v}(\mathbf{v} \cdot \mathbf{A}) = \mathbf{A} \times (\mathbf{v} \times \mathbf{a})/|\mathbf{v}|^2 + (\mathbf{A} \cdot \mathbf{v})\mathbf{a}/|\mathbf{v}|^2</math>
+
* '''A'''nd
 +
* '''R'''esource
 +
* '''C'''hannel
  
A concise alternative to the above is:
+
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)
  
: <math>\nabla_\mathbf{v}(\mathbf{v} \cdot \mathbf{A}) = (\mathbf{A} \cdot \mathbf{a})\mathbf{v}/|\mathbf{v}|^2 </math>
+
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)
  
The field on a moving charge <math>q</math> in a changing electromagnetic field becomes:
+
<youtube>https://www.youtube.com/watch?v=hFAOXdXZ5TM</youtube>
  
: <math>\mathbf{F}_{moving,dynamic}/q = - \nabla \varphi - ∂\mathbf{A}/∂t + \mathbf{v} \times \mathbf{B} + (\mathbf{A} \cdot \mathbf{a})\mathbf{v}/|\mathbf{v}|^2</math>
+
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)
  
The power on a moving charge <math>q</math> in a changing electromagnetic field becomes:
+
==Explaining "Altered" Lenz' Law Devices==
  
: <math>P_{moving,dynamic} = q \left[\ \left(-\nabla \varphi - ∂\mathbf{A}/∂t\right) \cdot \mathbf{v} + (\mathbf{A} \cdot \mathbf{a})(\mathbf{v} \cdot \mathbf{v})/|\mathbf{v}|^2\ \right]</math>
+
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:
: <math>P_{moving,dynamic} = q \left[\ \left(-\nabla \varphi - ∂\mathbf{A}/∂t\right) \cdot \mathbf{v} + \mathbf{A} \cdot \mathbf{a}\ \right]</math>
+
  
===Caveat===
+
<youtube>x3TFALmHtMw</youtube>
  
If the parent section to this subsection is correct, that the field on a moving charge <math>q</math> in a changing electromagnetic field becomes:
+
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)
  
: <math>\mathbf{F}_{moving,dynamic}/q = - \nabla \varphi - ∂\mathbf{A}/∂t + \mathbf{v} \times \mathbf{B} + (\mathbf{A} \cdot \mathbf{a})\mathbf{v}/|\mathbf{v}|^2</math>
+
==Explaining the Newman Motor==
  
The next question one may have is if the force generated by multiplying the field by the charge value (i.e. <math>\mathbf{F}_{moving,dynamic}</math>) is a <math>\frac{d(m \mathbf{v})}{dt}</math> type force or a <math>m \frac{d\mathbf{v}}{dt} = m\mathbf{a}</math> type force. If the latter is true, then the presence of a permanent magnet having a magnetic moment <math>\mathbf{m}</math> in the vicinity of a copper loop, which provides the vector potential <math>\mathbf{A}</math>, would cause the effective inductance to change via the last term on the right <math>(\mathbf{A} \cdot \mathbf{a})\mathbf{v}/|\mathbf{v}|^2</math>. However, the last term on the right is a <math>\mathbf{v} \frac{dm}{dt}</math> type term. If this is indeed the case, then such a term cannot contribute to acceleration and therefore may not contribute to any changes to velocity. Rather this force would be seen as changing the "effective-inertia" of the charge, implying that the term in essence carries with it, in part, the mass of other objects which provide the vector potential <math>\mathbf{A}</math>, in addition to the mass of the charge <math>q</math>. It could very be argued that mass is not a property of a particle-in-itself but rather a property of a particle coupled to the potentials and/or fields of other particles, such that the mass is the result of exchange forces between particles. If that is true, then to derive the force in terms of the common definition of <math>\mathbf{F} = m \frac{d\mathbf{v}}{dt} = m\mathbf{a}</math>, the last term should be subtracted from both sides. The result of doing so is the usual Lorentz Force:
+
{{See also|Memory Lane}}
  
: <math>\mathbf{F}_{Lorentz} = m \mathbf{a} = q \left[\ - \nabla \varphi - ∂\mathbf{A}/∂t + \mathbf{v} \times \mathbf{B}\ \right]</math>
+
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.
  
According to Kirk T. McDonald, professor of High-Energy Experimental Physics at Princeton,<ref>http://www.physics.princeton.edu/~mcdonald/</ref> the presence of a stationary permanent magnet does not affect the inductance of a stationary solenoid.<ref>Self Inductance of a Solenoid
+
<youtube>APq8HaeJpQ8</youtube>
with a Permanent-Magnet Core http://www.physics.princeton.edu/~mcdonald/examples/magsol.pdf</ref> This indicates that the extra term above cannot be allowed to affect in anyway the acceleration of a charge moving through the field of a magnet. This is great news for the S.H.O. Drive project because it means that it will not take additional energy to drive current to repel the central permanent magnet. The above predicted additional force term (depending on the gauge-dependent Magnetic Vector Potential) may be used to help enforce the rule of conservation of momentum (in 3-space) in the case of time-varying charge densities produced by non-radiating, discontinuous, yet constant current densities. We can add to this contribution the force term suggested by Koen J. van Vlaenderen in his 2015 paper "General Classical Electrodynamics"<ref>http://vixra.org/abs/1512.0297</ref> (as mentioned in the Comment Record section above) which is dependent on the gauge-dependent Electric Scalar Potential. Perhaps both of theses "forces" only change "mass" not velocity, or perhaps more accurately, they change the inertia that a charge carries by virtue of interactions via the Magnetic Vector Potential and Electric Scalar Potential.
+
  
So in the end, the actual accelerations of charges in 3-space would appear to follow exactly from what you would expect Special Relativity / Lorentz force to predict, but perhaps there are indeed underlying and non-apparent energy and force exchanges that may prove more apparent in future discoveries. Sincerely, ''[[User:S.H.O.|S.H.O.]] <sup>[[User_talk:S.H.O.|talk]]</sup>'' 14:53, 15 September 2016 (PDT)
+
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==
 
<references />
 
<references />
  

Latest revision as of 19: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