Difference between revisions of "Electromagnetic Potentials"

From S.H.O.
Jump to: navigation, search
(Scenario 2: Large drift velocities and negligible wire velocities)
(A New Idea: The Makerarc)
 
(99 intermediate revisions by the same user not shown)
Line 16: Line 16:
 
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.
 
causes transfers of energy ''within'' the magnetic rotor assembly (such as electron kinetic energy and inductive storage energy into wire kinetic energy (wired rotor example), or atomic electron kinetic energy to magnetic domain kinetic energy (permanent magnet rotor example)) that, in the low-frequency approximation, pretty much matches the amount of energy transfer between fields and the driving (stator) coils which is due to the transformer induction electric field <math>- \frac{\partial \mathbf{A}}{\partial t}</math> acting on currents in the coils from the relative motion of the magnetic rotor assembly.
  
As a result, the only way using conventional physics to explain various devices, such as the Marinov Motor<ref>http://redshift.vif.com/JournalFiles/Pre2001/V05NO3PDF/v05n3phi.pdf</ref>, the Distinti Paradox2<ref>http://www.distinti.com/docs/pdx/paradox2.pdf</ref>, and, especially, the Marinov Generator<ref name="Marinov Generator">http://overunity.com/14691/the-marinov-generator/</ref><ref name="Marinov Generator (paper)">http://www.overunityresearch.com/index.php?action=dlattach;topic=2470.0;attach=13897</ref> is to more accurately define the electric field <math>\mathbf{E}</math>, which is part of the full Lorentz Force equation, and then see if 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)
+
As a result, the only way using conventional physics to explain various devices, such as the Marinov Motor<ref name="Marinov Motor">http://redshift.vif.com/JournalFiles/Pre2001/V05NO3PDF/v05n3phi.pdf</ref>, the Distinti Paradox2<ref>http://www.distinti.com/docs/pdx/paradox2.pdf</ref>, and, especially, the Marinov Generator<ref name="Marinov Generator">http://overunity.com/14691/the-marinov-generator/</ref><ref name="Marinov Generator (paper)">http://www.overunityresearch.com/index.php?action=dlattach;topic=2470.0;attach=13897</ref> is to more accurately define the electric field <math>\mathbf{E}</math>, which is part of the full Lorentz Force equation, and then see if it leads to predictions that confirm observations, especially those found figure 9 of Cyril Smith's 2009 paper on the "Marinov Generator".<ref name="Marinov Generator (paper)"/> ''[[User:S.H.O.|S.H.O.]] <sup>[[User_talk:S.H.O.|talk]]</sup>'' 22:09, 5 March 2017 (PST)
  
 
'''Prior content in the "Comment Record" section:'''
 
'''Prior content in the "Comment Record" section:'''
Line 29: Line 29:
 
==Novel Force laws proposed by various researchers==
 
==Novel Force laws proposed by various researchers==
  
===Stefan Marinov's proposal===
+
===James Wesley's proposal===
  
Stefan Marinov proposed adding the "motional-transformer induction" on charge q:<ref>https://archive.org/stream/thornywayoftruthpart4maririch#page/104/mode/2up/search/motional-transformer+induction</ref>
+
James Wesley proposed adding the "motional induction" on charge q. In SI Units, this can be expressed as:<ref name="Marinov Motor"/>
  
: <math>(\mathbf{v_{wire}}\cdot\nabla)\mathbf{A}</math>
+
: <math>-q(\mathbf{v}\cdot\nabla)\mathbf{A}</math>
  
or rather:
+
to the Lorentz force.
  
: <math>-(\mathbf{v_{charge}}\cdot\nabla)\mathbf{A}</math>
+
The idea behind this was to explain an observation in an experiment involving a "Marinov Motor"<ref name="Marinov Motor"/> in which longitudinal induction forces were produced.
 
+
to the Lorentz force.
+
  
 
The extra term is equivalent to:<ref>http://www.overunityresearch.com/index.php?action=dlattach;topic=2470.0;attach=13908</ref>
 
The extra term is equivalent to:<ref>http://www.overunityresearch.com/index.php?action=dlattach;topic=2470.0;attach=13908</ref>
Line 55: Line 53:
 
<math>\mathbf{F} = q\left[-\nabla \varphi- \frac{\partial \mathbf{A}}{\partial t}+ \nabla_\mathbf{A}(\mathbf{v}\cdot\mathbf{A})-(\mathbf{v}\cdot\nabla)\mathbf{A} \right]</math>
 
<math>\mathbf{F} = q\left[-\nabla \varphi- \frac{\partial \mathbf{A}}{\partial t}+ \nabla_\mathbf{A}(\mathbf{v}\cdot\mathbf{A})-(\mathbf{v}\cdot\nabla)\mathbf{A} \right]</math>
  
Therefore, adding the extra term proposed by Stefan Marinov results in:
+
Therefore, adding the extra term proposed by Wesley results in:
  
 
<math>\mathbf{F} = q\left[-\nabla \varphi- \frac{\partial \mathbf{A}}{\partial t}+ \nabla_\mathbf{A}(\mathbf{v}\cdot\mathbf{A})-2(\mathbf{v}\cdot\nabla)\mathbf{A} \right]</math>
 
<math>\mathbf{F} = q\left[-\nabla \varphi- \frac{\partial \mathbf{A}}{\partial t}+ \nabla_\mathbf{A}(\mathbf{v}\cdot\mathbf{A})-2(\mathbf{v}\cdot\nabla)\mathbf{A} \right]</math>
Line 61: Line 59:
 
The problem with this modification:
 
The problem with this modification:
  
In the case of an electrical charge approaching a wire, this additional term proposed by Marinov '''doubles''' the force of deflection. This is not observed.
+
In the case of an electrical charge approaching a wire, this additional term proposed by Wesley would '''double''' the force of deflection. This is not observed.
  
Consider the vector potential due a current-carrying wire on the x-axis. Both the current and the vector potential of this current point in the +x direction. Now have a charge approaching this wire perpendicularly. Both the magnetic Lorentz force and the additional Marinov term predict the same force.
+
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.
 
+
Now consider the case where the wire is moving toward the charge. In this case, both the transformer induction <math>- \frac{\partial \mathbf{A}}{\partial t}</math> and the additional Marinov term predict the same force acting on the charge directed parallel to the vector potential of the current. The term proposed by Marinov '''doubles''' the observed force in these cases.
+
  
 
===Cyril Smith's proposal===
 
===Cyril Smith's proposal===
Line 103: Line 99:
 
==Explaining the Marinov Motor and Cyril Smith's "Marinov Generator" using Conventional Physics==
 
==Explaining the Marinov Motor and Cyril Smith's "Marinov Generator" using Conventional Physics==
  
In the sub-sections below, [https://en.wikipedia.org/wiki/Centimetre%E2%80%93gram%E2%80%93second_system_of_units CGS units] are used unless otherwise noted.
+
{| 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 E-field derived from the Liénard–Wiechert Potentials===
+
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>.
  
From the paper titled "Onoochin's Paradox" by Kirk T. McDonald, we have following statement:
+
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>.
  
<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>http://web.archive.org/web/20170318210550/http://puhep1.princeton.edu/~mcdonald/examples/ph501/ph501lecture24.pdf</ref>),
+
===The Liénard–Wiechert electric fields for electrically-neutral currents===
  
:<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>
+
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:
  
''where <math>\mathbf{a}</math> is the acceleration <math>\mathbf{a}</math> of the charge <math>q</math> at the present time.</blockquote>
+
<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>),
  
Let's consider the situation  where the acceleration of charge <math>q</math> is negligible:
+
:<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{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 <math>\mathbf{a}</math> is the acceleration <math>\mathbf{a}</math> of the charge <math>q</math> at the present time.</blockquote>
  
Next, let's clarify that <math>q</math> is the source of E field by attaching a subscript ''s'' to it:
+
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_s\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>
+
:<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:
 
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_s\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{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 receiver charge <math>q_r</math> located at <math>\mathbf{r}</math> at rest in the observer's inertial frame. The inertial observer and the charge <math>q_r</math> agree on what the electric field <math>E</math> is, they agree that that there is no magnetic force on <math>q_r</math>, and finally, they agree on the acceleration of <math>q_r</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>.
  
The above equation can be broken up into two parts, one based on the ''transverse relative velocity'' of the source <math>\mathbf{v}_{s\bot}</math>, and one based on the ''longitudinal relative velocity'' of the source<math>\mathbf{v}_{s\parallel}</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_s\frac{\mathbf{\hat{r}}}{r^2} \left(\frac{\mathbf{v}_{s\bot}^2}{2c^2} - 2\frac{\mathbf{v}_{s\parallel}^2}{2c^2} \right)</math>
+
:<math>\mathbf{E} = q\ \frac{\mathbf{\hat{r}}}{r^2} \frac{v^2}{2c^2} \left(1 - 3\ cos^2\theta \right)</math>
  
===Scenario 1: All charges confined to a line===
+
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:
  
We can consider in the surrounding environment a system of charges with arbitrary position and velocity which observe an electric field <math>\mathbf{E}</math> that is unique to each charge's rest frame.
+
:<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>
  
Let's start off by considering the simplest case where all the charges are distributed along a straight line and all of their velocities are tangent to this line. Therefore, the velocity-dependent electric field experienced by a charge <math>q_r</math> in its rest frame due a charge <math>q_s</math> is:
+
Since the '''relative ''azimuthal'' velocity''' and the '''relative ''radial'' velocity''' are orthogonal, we can express the following:
  
:<math>\mathbf{E} = - q_s\frac{\mathbf{\hat{r}}}{r^2} \frac{\left(\mathbf{v}_{s}-\mathbf{v}_{r}\right)^2}{c^2}</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>\mathbf{E} = - q_s\frac{\mathbf{\hat{r}}}{r^2} \frac{\left(\mathbf{v}_{s \to r}\right)^2}{c^2}</math>
+
  
Let's consider that along this line there is an element of a wire at velocity <math>\mathbf{v}_{wire}</math> relative to charge <math>q_r</math> carrying loose electrons with drift velocity <math>\mathbf{v}_{drift}</math> relative to that wire. And let's say the total electrical charge of that wire <math>Q_{wire}=+Q</math> + loose electrons <math>Q_{drift}=-Q</math> is zero. What is the electric field created by this current element according to a receiver charge <math>q_{r}</math> moving at velocity <math>\mathbf{v}_r</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{E} = - \frac{\mathbf{\hat{r}}}{r^2} \left( \frac{Q \mathbf{v}_{wire \to r}^2}{c^2} + \frac{-Q \left(\mathbf{v}_{wire \to r}+\mathbf{v}_{drift}\right)^2}{c^2}\right)</math>
+
:<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>
  
:<math>\mathbf{E} = - \frac{\mathbf{\hat{r}}}{r^2} \left( \frac{Q \mathbf{v}_{wire \to r}^2}{c^2} - \frac{Q \left(\mathbf{v}_{wire \to r}+\mathbf{v}_{drift}\right)^2}{c^2}\right)</math>
+
Therefore:
  
:<math>\mathbf{E} = - \frac{\mathbf{\hat{r}}}{r^2} \left( \frac{2 Q \mathbf{v}_{wire \to r}\mathbf{v}_{drift} + \mathbf{v}_{drift}^2}{c^2}\right)</math>
+
:<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>
  
To simplify this result, let's consider the average velocity of charges comprising the wire+current:
+
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{v}_{rel} = \mathbf{v}_{wire \to r} + \mathbf{v}_{drift}/2</math>
+
:<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>
  
This is effectively the relative velocity of a current element <math>\mathbf{I}_s dl = Q_{drift} \mathbf{v}_{drift}</math> with respect to receiver charge <math>q_r</math>
+
Consider the existence of four charges:
  
From here, we can find the alternate expression for electric field at charge <math>q_{r}</math> due to the current element <math>\mathbf{I}_s dl</math>:
+
:<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
  
:<math>\mathbf{E} = - \frac{\mathbf{\hat{r}}}{r^2} \left( \frac{Q \left(\mathbf{v}_{avg} - \mathbf{v}_{drift}/2\right)^2}{c^2} + \frac{-Q \left(\mathbf{v}_{avg} + \mathbf{v}_{drift}/2\right)^2}{c^2}\right)</math>
+
Consider the existence of their corresponding current elements:
  
:<math>\mathbf{E} = - \frac{\mathbf{\hat{r}}}{r^2} \left( \frac{Q \left(\mathbf{v}_{avg} - \mathbf{v}_{drift}/2\right)^2}{c^2} - \frac{Q \left(\mathbf{v}_{avg} + \mathbf{v}_{drift}/2\right)^2}{c^2}\right)</math>
+
:<math>id\mathbf{l}</math> : the source current element
 +
:<math>Id\mathbf{L}</math> : the target current element
  
:<math>\mathbf{E} = - \frac{\mathbf{\hat{r}}}{r^2} \left( \frac{-Q \mathbf{v}_{avg} \mathbf{v}_{drift}}{c^2}\right)</math>
+
Where:
  
:<math>\mathbf{E} = - \frac{\mathbf{\hat{r}}}{r^2} \left( \frac{Q_{drift} \mathbf{v}_{avg} \mathbf{v}_{drift}}{c^2}\right)</math>
+
:<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.
  
:<math>\mathbf{E} = - \frac{\mathbf{\hat{r}}}{r^2} \left( \frac{\mathbf{I}_s \mathbf{v}_{avg}}{c^2}\right)</math>
+
The current elements are equal to:
  
====Special Scenario 1.1: Forces between co-linear current elements confined to a straight line====
+
:<math>id\mathbf{l} = q_- \mathbf{v}_d</math>
 +
:<math>Id\mathbf{L} = Q_- \mathbf{V}_d</math>
  
If we wish to consider the force on charge <math>q_s</math> due to the electric field of this current element, we simply multiply by target charge <math>q_r</math>.
+
Where:
  
:<math>\mathbf{F} = - \frac{\mathbf{\hat{r}}}{r^2} \left( \frac{\mathbf{I}_s q_r \mathbf{v}_{avg}}{c^2}\right)</math>
+
:<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.
  
If <math>q_s</math> was a loose electron embedded in another electrically-neutral wire element, let <math>\mathbf{I}_r dl</math> stand for this target current element, analogous to the source current element <math>\mathbf{I}_s dl</math>. Therefore the force on current element <math>\mathbf{I}_r dl</math> due to current element <math>\mathbf{I}_s dl</math> is:
+
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:
  
:<math>\mathbf{F} = - \frac{\mathbf{\hat{r}}}{r^2} \left( \frac{\mathbf{I}_s \mathbf{I}_r}{c^2}\right)</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>
 +
:Force <math>\mathbf{F}_{++}</math> on <math>Q_+</math> by field <math>\mathbf{E}_{++}</math> of <math>q_+</math>
  
In '''SI units''', this becomes:
+
These four forces are dependent on four different relative velocities (source velocity w.r.t. target velocity):
  
:<math>\mathbf{F} = - \frac{\mathbf{\hat{r}}}{4 \pi \epsilon r^2} \left( \frac{\mathbf{I}_s \mathbf{I}_r}{c^2}\right)</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>
 +
:Relative velocity <math>\mathbf{v}_{++} = \mathbf{v}_+ - \mathbf{V}_+</math>
  
This is equivalent to:
+
The corresponding relative speeds are:
  
:<math>\mathbf{F} = - \left( \mathbf{I}_r \dot\ \nabla \right)\mathbf{A_s}</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_{x++} = \left(v_+ - V_+\right)_x</math>
  
This is essentially the "Marinov Force", but with restrictions. Like the Marinov Force, it exists without necessitating a magnetic field <math>\mathbf{B}</math> and so may exist in regions where the curl of <math>\mathbf{A_s}</math> is zero. However, this scenario is constricted to the case that "all charges are confined to a line" with the restriction that the current elements are electrically neutral.
+
: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>
  
===Scenario 2: Large drift velocities and negligible wire velocities===
+
: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 magnetization currents embodied inside permanent magnets consist of electrons with very high effective velocities per the Virial theorem which relates the kinetic energy to the potential energy of the system<ref>https://en.wikipedia.org/wiki/Virial_theorem</ref>. Therefore, in the case of permanent magnets interacting with other permanent magnets, one can simply consider the effective velocity of the charges comprising the effective magnetization currents and disregard the velocities of the atoms which confine these currents in place.
+
The two drift velocities are:
  
Let <math>\mathbf{v}_{rel}</math> be the relative velocity of charge <math>q_s</math> respect to <math>q_r</math>. In order to account for the forces due to the relative motion with respect to the positive charges which screen the Coulomb electric field of <math>q_s</math>, the velocity <math>\mathbf{v}_r</math> of charge <math>q_r</math> is added to <math>\mathbf{v}_{rel}</math>. Therefore the effective relative velocity <math>v_{net}</math> is equal to:
+
:<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
  
:<math>v_{net} = \mathbf{v}_{rel} + \mathbf{v}_r = \mathbf{v}_s - \mathbf{v}_r + \mathbf{v}_r = \mathbf{v}_s</math>
+
Let the effective velocities of the current elements be:
  
As indicated in the parent section, the electric field due to charge <math>q_s</math> moving at velocity <math>\mathbf{v}_s</math> is equal to:
+
:<math>\mathbf{v} = \left( \mathbf{v}_- + \mathbf{v}_+ \right)/2</math>
 +
:<math>\mathbf{V} = \left( \mathbf{V}_- + \mathbf{V}_+ \right)/2</math>
  
:<math>\mathbf{E} = q_s\frac{\mathbf{\hat{r}}}{r^2} \left(\frac{\mathbf{v}_{s\bot}^2}{2c^2} - 2\frac{\mathbf{v}_{s\parallel}^2}{2c^2} \right)</math>
+
So the effective relative velocity between the current elements (source velocity w.r.t. target velocity) is:
  
The electric force on <math>q_r</math> by <math>q_s</math> is therefore:
+
:<math>\mathbf{v}_{rel} = \mathbf{v} - \mathbf{V}</math>
  
:<math>\mathbf{F} = q_r q_s\frac{\mathbf{\hat{r}}}{r^2} \left(\frac{\mathbf{v}_{s\bot}^2}{2c^2} - 2\frac{\mathbf{v}_{s\parallel}^2}{2c^2} \right)</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:
  
====Special Scenario 2.1: Forces between magnetization currents from magnets of the same grade====
+
:<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.
  
In the case where the magnets are of the same grade (e.g. N52), the effective drift speed of the magnetization currents are identical. This makes it possible to calculate the above using current elements <math>\mathbf{I}_s</math> and <math>\mathbf{I}_r</math>.
+
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>.
  
The differential version of this is:
+
: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>
  
:<math>d^2\mathbf{F} = \frac{\mathbf{\hat{r}}}{r^2} \left(\frac{\mathbf{I}_{s\bot}^2}{2c^2} - 2\frac{\mathbf{I}_{s\parallel}^2}{2c^2} \right)</math>
+
The corresponding relative speeds are:
  
In '''SI units''', this becomes:
+
: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>
  
:<math>d^2\mathbf{F} = \frac{\mathbf{\hat{r}}}{4 \pi \epsilon r^2} \left(\frac{\mathbf{I}_{s\bot}^2}{2c^2} - 2\frac{\mathbf{I}_{s\parallel}^2}{2c^2} \right)</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>
  
In this special scenario:
+
: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>
  
* <math>\mathbf{I}_{s\parallel}</math> is the component of <math>\mathbf{I}_{s}</math> which is parallel to <math>\mathbf{r}</math>.
+
The equation for the electric field contains squared values of the speed. As noted before:
* <math>\mathbf{I}_{s\bot}</math> is the component of <math>\mathbf{I}_{s}</math> which is perpendicular to <math>\mathbf{r}</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>
 +
 
 +
The electric fields on charges <math>Q_+</math> and <math>Q_-</math>, in their own respective and distinct rest frames, are as follows:
 +
 
 +
:<math>\mathbf{E}_- = \frac{\mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( q_- \left(v_{x--}^2 + v_{y--}^2 - 2 v_{r--}^2\right) + q_+ \left(v_{x+-}^2 + v_{y+-}^2 - 2 v_{r+-}^2\right) \right)</math>
 +
:<math>\mathbf{E}_+ = \frac{\mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( q_- \left(v_{x-+}^2 + v_{y-+}^2 - 2 v_{r-+}^2\right) + q_+ \left(v_{x++}^2 + v_{y++}^2 - 2 v_{r++}^2\right) \right)</math>
 +
 
 +
===The Forces due to the Liénard–Wiechert electric fields for electrically-neutral currents===
 +
 
 +
It may be more helpful to calculate the forces due to relative radial velocities separately from the forces due to relative azimuthal velocities. Therefore:
 +
 
 +
:<math>\mathbf{F}_x = \frac{\mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( Q_- q_- v_{x--}^2 + Q_- q_+ v_{x+-}^2 + Q_+ q_- v_{x-+}^2 + Q_+ q_+ v_{x++}^2 \right)</math>
 +
:<math>\mathbf{F}_y = \frac{\mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( Q_- q_- v_{y--}^2 + Q_- q_+ v_{y+-}^2 + Q_+ q_- v_{y-+}^2 + Q_+ q_+ v_{y++}^2 \right)</math>
 +
:<math>\mathbf{F}_r = \frac{\mathbf{\hat{r}}}{r^2} \frac{- 2}{2c^2} \left( Q_- q_- v_{r--}^2 + Q_- q_+ v_{r+-}^2 + Q_+ q_- v_{r-+}^2 + Q_+ q_+ v_{r++}^2 \right)</math>
 +
 
 +
It follows that <math>\mathbf{F}_x</math> and <math>\mathbf{F}_y</math> are functions of currents perpendicular to radial vector <math>\mathbf{r}</math> while <math>\mathbf{F}_r</math> is a function of currents co-linear with radial vector <math>\mathbf{r}</math>
 +
 
 +
All charges (<math>Q_-</math>, <math>Q_+</math>, <math>q_-</math>, and <math>q_+</math>) may contribute simultaneously to the azimuthally-directed ('''transverse''') currents in <math>x</math> and <math>y</math> and the radially-directed ('''longitudinal''') currents in <math>r</math>.
 +
 
 +
It will be very advantageous to simplify these formulas. For electrically-neutral currents, we can recognize the following:
 +
 
 +
:<math>q_+ = - q_-</math>
 +
:<math>Q_+ = - Q_-</math>
 +
 
 +
This allows us to factor out <math>Q_+ q_+</math> with the following result:
 +
 
 +
:<math>\mathbf{F}_x = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( v_{x--}^2 - v_{x+-}^2 - v_{x-+}^2 + v_{x++}^2 \right)</math>
 +
:<math>\mathbf{F}_y = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( v_{y--}^2 - v_{y+-}^2 - v_{y-+}^2 + v_{y++}^2 \right)</math>
 +
:<math>\mathbf{F}_r = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{- 2}{2c^2} \left( v_{r--}^2 - v_{r+-}^2 - v_{r-+}^2 + v_{r++}^2 \right)</math>
 +
 
 +
Next, we will work on simplifying the contents within the parentheses. To make matters simpler, we will move the subscript to the lower right corner of the parentheses.
 +
 
 +
:<math>\mathbf{F}_x = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 \right)_x</math>
 +
:<math>\mathbf{F}_y = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 \right)_y</math>
 +
:<math>\mathbf{F}_r = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{- 2}{2c^2} \left( v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 \right)_r</math>
 +
 
 +
These can be distributed back to the contents within the parentheses after they are substituted for a different expression in terms of the currents. The terms in the parentheses (disregarding their sign) are as follows:
 +
 
 +
:<math>v_{--}^2 = \left(u + v_{rel} - U\right)^2 = u^2 + v_{rel}^2 + U^2 + (2u - 2U) v_{rel} - 2 u U</math>
 +
:<math>v_{+-}^2 = \left(-u + v_{rel} - U\right)^2 = u^2 + v_{rel}^2 + U^2 + (-2u - 2U) v_{rel} + 2 u U</math>
 +
:<math>v_{-+}^2 = \left(u + v_{rel} + U\right)^2 = u^2 + v_{rel}^2 + U^2 + (2u + 2U) v_{rel} + 2 u U</math>
 +
:<math>v_{++}^2 = \left(-u + v_{rel} + U\right)^2 = u^2 + v_{rel}^2 + U^2 + (-2u + 2U) v_{rel} - 2 u U</math>
 +
 
 +
Therefore sum of the forces on <math>Q_-</math> and <math>Q_+</math> depends on:
 +
 
 +
:<math>v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 = ((2u - 2U) - (-2u - 2U) - (2u + 2U) + (-2u + 2U))v_{rel} + ((-2) - (2) - (2) + (-2)) u U</math>
 +
:<math>v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 = ((2u - 2U) + (2u + 2U) - (2u + 2U) - (2u - 2U))v_{rel} - 8 u U</math>
 +
:<math>v_{--}^2 - v_{+-}^2 - v_{-+}^2 - v_{+-}^2 + v_{++}^2 = ((-2U) + (2U) - (2U) - (- 2U))v_{rel} - 8 u U</math>
 +
:<math>v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 = (0)v_{rel} - 8 u U</math>
 +
:<math>v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 = - 8 u U</math>
 +
:<math>v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 = - 8 (v_d/2)(V_d/2)</math>
 +
:<math>v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 = - 2 (v_d)(V_d)</math>
 +
 
 +
The force on <math>Q_-</math> depends on:
 +
 
 +
:<math>v_{--}^2 - v_{+-}^2 = ((2u - 2U) - (-2u - 2U))v_{rel} + ((-2) - (2)) u U</math>
 +
:<math>v_{--}^2 - v_{+-}^2 = 4 u v_{rel} - 4 u U</math>
 +
 
 +
The force on <math>Q_+</math> depends on:
 +
 
 +
:<math>v_{++}^2 - v_{-+}^2 = ((-2u + 2U) - (2u + 2U))v_{rel} + (- (2) + (-2)) u U</math>
 +
:<math>v_{++}^2 - v_{-+}^2 = -4 u v_{rel} - 4 u U</math>
 +
 
 +
Recalling that:
 +
 
 +
:<math>\mathbf{F}_x = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 \right)_x</math>
 +
:<math>\mathbf{F}_y = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 \right)_y</math>
 +
:<math>\mathbf{F}_r = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{- 2}{2c^2} \left( v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 \right)_r</math>
 +
:<math>v_{--}^2 - v_{+-}^2 - v_{-+}^2 + v_{++}^2 = - 2 (v_d)(V_d)</math>
 +
 
 +
Substitution yields:
 +
 
 +
:<math>\mathbf{F}_x = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( - 2 (v_d)(V_d) \right)_x</math>
 +
:<math>\mathbf{F}_y = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{2c^2} \left( - 2 (v_d)(V_d) \right)_y</math>
 +
:<math>\mathbf{F}_r = \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{- 2}{2c^2} \left( - 2 (v_d)(V_d) \right)_r</math>
 +
 
 +
Simplified:
 +
 
 +
:<math>\mathbf{F}_x = - \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( (v_d)(V_d) \right)_x</math>
 +
:<math>\mathbf{F}_y = - \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( (v_d)(V_d) \right)_y</math>
 +
:<math>\mathbf{F}_r = + \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{2}{c^2} \left( (v_d)(V_d) \right)_r</math>
 +
 
 +
As stated earlier, the current elements are equal to:
 +
 
 +
:<math>id\mathbf{l} = q_- \mathbf{v}_d</math>
 +
:<math>Id\mathbf{L} = Q_- \mathbf{V}_d</math>
 +
 
 +
We can now substitute the currents into the equation. First we substitute <math>Q_- q_-</math> for <math>Q_+ q_+</math>:
 +
 
 +
:<math>\mathbf{F}_x = - \frac{Q_- q_- \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( (v_d)(V_d) \right)_x</math>
 +
:<math>\mathbf{F}_y = - \frac{Q_- q_- \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( (v_d)(V_d) \right)_y</math>
 +
:<math>\mathbf{F}_r = + \frac{Q_- q_- \mathbf{\hat{r}}}{r^2} \frac{2}{c^2} \left( (v_d)(V_d) \right)_r</math>
 +
 
 +
Next, we assign each charge with their corresponding drift velocities:
 +
 
 +
:<math>\mathbf{F}_x = - \frac{\mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( (Q_- V_d) (q_- v_d) \right)_x</math>
 +
:<math>\mathbf{F}_y = - \frac{\mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( (Q_- V_d) (q_- v_d) \right)_y</math>
 +
:<math>\mathbf{F}_r = + \frac{\mathbf{\hat{r}}}{r^2} \frac{2}{c^2} \left( (Q_- V_d) (q_- v_d) \right)_r</math>
 +
 
 +
Next, we make a substitution for the current elements:
 +
 
 +
:<math>\mathbf{F}_x = - \frac{\mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( (Id\mathbf{L}) \cdot (id\mathbf{l}) \right)_x</math>
 +
:<math>\mathbf{F}_y = - \frac{\mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( (Id\mathbf{L}) \cdot (id\mathbf{l}) \right)_y</math>
 +
:<math>\mathbf{F}_r = + \frac{\mathbf{\hat{r}}}{r^2} \frac{2}{c^2} \left( (Id\mathbf{L}) \cdot (id\mathbf{l}) \right)_r</math>
 +
 
 +
This can be written as:
 +
 
 +
:<math>\mathbf{F}_x = - \frac{\mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( Ii (d\mathbf{L} \cdot d\mathbf{l}) \right)_x</math>
 +
:<math>\mathbf{F}_y = - \frac{\mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( Ii (d\mathbf{L} \cdot d\mathbf{l}) \right)_y</math>
 +
:<math>\mathbf{F}_r = + \frac{\mathbf{\hat{r}}}{r^2} \frac{2}{c^2} \left( Ii (d\mathbf{L} \cdot d\mathbf{l}) \right)_r</math>
 +
 
 +
Or:
 +
 
 +
:<math>\mathbf{F}_x = - \frac{\mathbf{\hat{r}}}{r^2} \frac{Ii}{c^2} \left( d\mathbf{L} \cdot d\mathbf{l} \right)_x</math>
 +
:<math>\mathbf{F}_y = - \frac{\mathbf{\hat{r}}}{r^2} \frac{Ii}{c^2} \left( d\mathbf{L} \cdot d\mathbf{l} \right)_y</math>
 +
:<math>\mathbf{F}_r = + \frac{\mathbf{\hat{r}}}{r^2} \frac{2Ii}{c^2} \left( d\mathbf{L} \cdot d\mathbf{l} \right)_r</math>
 +
 
 +
Or:
 +
 
 +
:<math>\mathbf{F}_x = - \frac{\mathbf{\hat{r}}}{r^2} \frac{Ii}{c^2} \left( d\mathbf{L}_x \cdot d\mathbf{l}_x \right)</math>
 +
:<math>\mathbf{F}_y = - \frac{\mathbf{\hat{r}}}{r^2} \frac{Ii}{c^2} \left( d\mathbf{L}_y \cdot d\mathbf{l}_y \right)</math>
 +
:<math>\mathbf{F}_r = + \frac{\mathbf{\hat{r}}}{r^2} \frac{2Ii}{c^2} \left( d\mathbf{L}_r \cdot d\mathbf{l}_r \right)</math>
 +
 
 +
Adding the forces together, their sum is:
 +
 
 +
:<math>\mathbf{F} = \frac{\mathbf{\hat{r}}}{r^2} \frac{Ii}{c^2} \left( 2 d\mathbf{L}_r \cdot d\mathbf{l}_r - d\mathbf{L}_x \cdot d\mathbf{l}_x - d\mathbf{L}_y \cdot d\mathbf{l}_y \right)</math>
 +
:<math>\mathbf{F} = \frac{\mathbf{\hat{r}}}{r^2} \frac{Ii}{c^2} \left( 2 d\mathbf{L}_r \cdot d\mathbf{l}_r - d\mathbf{L}_\theta \cdot d\mathbf{l}_\theta \right)</math>
 +
 
 +
In S.I. Units:
 +
 
 +
:<math>\mathbf{F} = \frac{\mathbf{\hat{r}}}{4\pi\epsilon_0r^2} \frac{Ii}{c^2} \left( 2 d\mathbf{L}_r \cdot d\mathbf{l}_r - d\mathbf{L}_\theta \cdot d\mathbf{l}_\theta \right)</math>
 +
 
 +
Or:
 +
 
 +
:<math>\mathbf{F} = \frac{\mu_0 I i \mathbf{\hat{r}}}{4\pi r^2} \left( 2 d\mathbf{L}_r \cdot d\mathbf{l}_r - d\mathbf{L}_\theta \cdot d\mathbf{l}_\theta \right)</math>
 +
 
 +
The second differential is:
 +
 
 +
:<math>d^2\mathbf{F} = \frac{\mu_0 \mathbf{\hat{r}}}{4\pi r^2} \left( 2 \mathbf{I}_r \cdot \mathbf{i}_r - \mathbf{I}_\theta \cdot \mathbf{i}_\theta \right)</math>
 +
 
 +
===Special Scenario: No relative motion between positive charges===
 +
 
 +
In the case that the positive charges <math>q_+</math> and <math>Q_+</math> are essentially stationary, '''we can simplify the derivation of the field and force equations in the rest frame of the positive charges'''.
 +
 
 +
As stated in the parent section, the two drift velocities are:
 +
 
 +
:<math>\mathbf{v}_d = \mathbf{v}_- - \mathbf{v}_+</math> : The drift velocity of the loose electrons of the source current element
 +
:<math>\mathbf{V}_d = \mathbf{V}_- - \mathbf{V}_+</math> : The drift velocity of the loose electrons of the target current element
 +
 
 +
These become:
 +
 
 +
:<math>\mathbf{v}_d = \mathbf{v}_-</math> : The drift velocity of the loose electrons of the source current element
 +
:<math>\mathbf{V}_d = \mathbf{V}_-</math> : The drift velocity of the loose electrons of the target current element
 +
 
 +
As stated in the parent section, the effective velocities of the current elements are:
 +
 
 +
:<math>\mathbf{v} = \left( \mathbf{v}_- + \mathbf{v}_+ \right)/2</math>
 +
:<math>\mathbf{V} = \left( \mathbf{V}_- + \mathbf{V}_+ \right)/2</math>
 +
 
 +
These become:
 +
 
 +
:<math>\mathbf{v} = \mathbf{v}_-  /2</math>
 +
:<math>\mathbf{V} = \mathbf{V}_- / 2</math>
 +
 
 +
As stated in the parent section, a new variable, the deviation velocity, is defined as one-half the drift velocity of the electrons:
 +
 
 +
:<math>\mathbf{u} = \mathbf{v}_d / 2</math> is the deviation velocity of the source current element.
 +
:<math>\mathbf{U} = \mathbf{V}_d / 2</math> is the deviation velocity of the target current element.
 +
 
 +
Therefore:
 +
 
 +
:<math>\mathbf{v} = \mathbf{u} </math>
 +
:<math>\mathbf{V} = \mathbf{U} </math>
 +
 
 +
The effective relative velocity between the current elements (source velocity w.r.t. target velocity) is:
 +
 
 +
:<math>\mathbf{v}_{rel} = \mathbf{v} - \mathbf{V}</math>
 +
 
 +
Therefore:
 +
 
 +
:<math>\mathbf{v}_{rel} = \mathbf{u} - \mathbf{U}</math>
 +
 
 +
As derived in the parent section, all four relative velocities can be expressed in terms of the deviation velocities <math>\mathbf{u}</math> and <math>\mathbf{U}</math> together with the relative velocity <math>\mathbf{v}_{rel}</math>.
 +
 
 +
:Relative velocity <math>\mathbf{v}_{--} = \mathbf{v}_- - \mathbf{V}_- = \mathbf{u} + \mathbf{v}_{rel} - \mathbf{U}</math>
 +
:Relative velocity <math>\mathbf{v}_{-+} = \mathbf{v}_- - \mathbf{V}_+ = \mathbf{u} + \mathbf{v}_{rel} + \mathbf{U}</math>
 +
:Relative velocity <math>\mathbf{v}_{+-} = \mathbf{v}_+ - \mathbf{V}_- = -\mathbf{u} + \mathbf{v}_{rel} - \mathbf{U}</math>
 +
:Relative velocity <math>\mathbf{v}_{++} = \mathbf{v}_+ - \mathbf{V}_+ = -\mathbf{u} + \mathbf{v}_{rel} + \mathbf{U}</math>
 +
 
 +
By substituting for <math>\mathbf{v}_{rel}</math>, we get:
 +
 
 +
:Relative velocity <math>\mathbf{v}_{--} = \mathbf{v}_- - \mathbf{V}_- = 2\mathbf{u} - 2\mathbf{U}</math>
 +
:Relative velocity <math>\mathbf{v}_{-+} = \mathbf{v}_- - \mathbf{V}_+ = 2\mathbf{u}</math>
 +
:Relative velocity <math>\mathbf{v}_{+-} = \mathbf{v}_+ - \mathbf{V}_- = -2\mathbf{U}</math>
 +
:Relative velocity <math>\mathbf{v}_{++} = \mathbf{v}_+ - \mathbf{V}_+ = 0</math>
 +
 
 +
As derived in the parent section, the following is a series sum of terms (a function of relative speeds) which will be used to help calculate the forces between currents:
 +
 
 +
:<math>v_{--}^2 - v_{-+}^2 - v_{+-}^2 + v_{++}^2</math>
 +
 
 +
Based on the above results, we have:
 +
 
 +
:<math>v_{--}^2 - v_{-+}^2 - v_{+-}^2 + v_{++}^2 = (2u - 2U)^2 - (2u)^2 - (-2U)^2 + 0</math>
 +
:<math>v_{--}^2 - v_{-+}^2 - v_{+-}^2 + v_{++}^2 = (4u^2 + 4U^2 - 8uU) - 4u^2 - 4U^2 + 0</math>
 +
:<math>v_{--}^2 - v_{-+}^2 - v_{+-}^2 + v_{++}^2 = - 8uU</math>
 +
:<math>v_{--}^2 - v_{-+}^2 - v_{+-}^2 + v_{++}^2 = - 8 (v_{d}/2)(V_{d}/2)</math>
 +
:<math>v_{--}^2 - v_{-+}^2 - v_{+-}^2 + v_{++}^2 =  -2 (v_{d})(V_{d})</math>
 +
 
 +
This produces the same results as the parent section titled "Explaining the Marinov Motor and Cyril Smith's 'Marinov Generator' using Conventional Physics" and the previous section titled "The Forces due to the Liénard–Wiechert electric fields for electrically-neutral currents":
 +
 
 +
In S.I. Units:
 +
 
 +
:<math>\mathbf{F} = \frac{\mathbf{\hat{r}}}{4\pi\epsilon_0r^2} \frac{Ii}{c^2} \left( 2 d\mathbf{L}_r \cdot d\mathbf{l}_r - d\mathbf{L}_\theta \cdot d\mathbf{l}_\theta \right)</math>
 +
 
 +
Or:
 +
 
 +
:<math>\mathbf{F} = \frac{\mu_0 I i \mathbf{\hat{r}}}{4\pi r^2} \left( 2 d\mathbf{L}_r \cdot d\mathbf{l}_r - d\mathbf{L}_\theta \cdot d\mathbf{l}_\theta \right)</math>
 +
 
 +
The second differential is:
 +
 
 +
:<math>d^2\mathbf{F} = \frac{\mu_0 \mathbf{\hat{r}}}{4\pi r^2} \left( 2 \mathbf{I}_r \cdot \mathbf{i}_r - \mathbf{I}_\theta \cdot \mathbf{i}_\theta \right)</math>
 +
 
 +
===Deriving James Wesley's additional force term in the case of co-linear current elements===
 +
 
 +
As pointed out above in the sub-section titled "James Wesley's proposal", James Wesley proposed adding the "motional induction" on charge <math>q</math>. In SI Units, this can be expressed as:<ref name="Marinov Motor"/>
 +
 
 +
: <math>-q(\mathbf{v}\cdot\nabla)\mathbf{A}</math>
 +
 
 +
How may we derive the same result for current elements located and oriented along the same line? Equation 2b of the paper titled "Observations of the Marinov Motor"<ref name="Marinov Motor"/> can be adapted to the form above, resulting in:
 +
 
 +
: <math>\mathbf{F} = -Q_- \left(\mathbf{V}_d \cdot \nabla \right) \mathbf{A} = Q_- \left(-\mathbf{V}_d \cdot \nabla \right) \mathbf{A} = - I \left( d\mathbf{L} \cdot \nabla \right) \mathbf{A} </math>
 +
 
 +
In the case that we are dealing only with co-linear current elements, the result can be expressed as:
 +
 
 +
: <math>\mathbf{F}_r = \left( -Q_- \left(\mathbf{V}_d \cdot \nabla \right) \mathbf{A} = Q_- \left(-\mathbf{V}_d \cdot \nabla \right) \mathbf{A} = - I \left( d\mathbf{L} \cdot \nabla \right) \mathbf{A} \right)_r</math>
 +
 
 +
As stated in the parent section, the force on <math>Q_-</math> depends on:
 +
 
 +
:<math>v_{--}^2 - v_{+-}^2 = ((2u - 2U) - (-2u - 2U))v_{rel} + ((-2) - (2)) u U</math>
 +
:<math>v_{--}^2 - v_{+-}^2 = 4 u v_{rel} - 4 u U</math>
 +
 
 +
As stated in the parent section, The force on <math>Q_+</math> depends on:
 +
 
 +
:<math>v_{++}^2 - v_{-+}^2 = ((-2u + 2U) - (2u + 2U))v_{rel} + (- (2) + (-2)) u U</math>
 +
:<math>v_{++}^2 - v_{-+}^2 = -4 u v_{rel} - 4 u U</math>
 +
 
 +
In the case that that there is no relative velocity <math>v_{rel}</math> between the current elements:
 +
 
 +
:<math>v_{--}^2 - v_{+-}^2 = - 4 u U = -4 (v_d/2)(V_d/2) = - (v_d)(V_d)</math>
 +
:<math>v_{++}^2 - v_{-+}^2 = - 4 u U = -4 (v_d/2)(V_d/2) = - (v_d)(V_d)</math>
 +
 
 +
As stated in the section "The Forces due to the Liénard–Wiechert electric fields for electrically-neutral currents", the electric forces on charges <math>Q_-</math> and <math>Q_+</math> are:
 +
 
 +
:<math>\mathbf{F}_{r-} = + \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( v_{--}^2 - v_{+-}^2  \right)_r</math>
 +
:<math>\mathbf{F}_{r+} = + \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( v_{++}^2 - v_{-+}^2  \right)_r</math>
 +
 
 +
From this and above, it follows that:
 +
 
 +
:<math>\mathbf{F}_{r-} = + \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( - (v_d)(V_d)  \right)_r</math>
 +
:<math>\mathbf{F}_{r+} = + \frac{Q_+ q_+ \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( - (v_d)(V_d)  \right)_r</math>
 +
 
 +
Therefore:
 +
 
 +
:<math>\mathbf{F}_{r-} = + \frac{Q_- q_- \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( - (v_d)(V_d) \right)_r</math>
 +
:<math>\mathbf{F}_{r+} = + \frac{Q_+ q_- \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( + (v_d)(V_d) \right)_r</math>
 +
 
 +
The effective electric fields experience by each charge is:
 +
 
 +
:<math>\mathbf{E}_{r-} = + \frac{q_- \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( - (v_d)(V_d) \right)_r</math>
 +
:<math>\mathbf{E}_{r+} = + \frac{q_- \mathbf{\hat{r}}}{r^2} \frac{1}{c^2} \left( + (v_d)(V_d) \right)_r</math>
 +
 
 +
Given that:
 +
 
 +
:<math>Q_+ = - Q_-</math>
 +
:<math>\mathbf{U} = \mathbf{v}_d / 2</math>
 +
:<math>\mathbf{v}_d = \mathbf{U} - (-\mathbf{U})</math>
 +
:<math>Id\mathbf{L} = Q_- \mathbf{V}_d = Q_- (+ \mathbf{U}) + Q_+ (-\mathbf{U})</math>
 +
 
 +
It can be shown that:
 +
 
 +
:<math>\mathbf{F}_{r+} = \left( - Q_+ \left((-\mathbf{U}) \cdot \nabla \right) \mathbf{A} \right)_r</math>
 +
:<math>\mathbf{F}_{r-} = \left( - Q_- \left((+\mathbf{U}) \cdot \nabla \right) \mathbf{A} \right)_r</math>
 +
:<math>\mathbf{F}_r = \left( - Q_- \left(\mathbf{V}_d \cdot \nabla \right) \mathbf{A} \right)_r = \mathbf{F}_{r-} + \mathbf{F}_{r+}</math>
 +
 
 +
This is based on the directional derivative of vector potential <math>\mathbf{A}</math> along <math>Id\mathbf{L}</math>.
 +
 
 +
It is important to note that the total force <math>\mathbf{F}_r</math> in <math>r</math> remains unaffected by the relative motion between the target current element <math>I_rd\mathbf{L}_r</math> and the source current element <math>i_rd\mathbf{l}_r</math> responsible for vector potential <math>\mathbf{A}_r</math>. However, forces <math>\mathbf{F}_{r-}</math> and <math>\mathbf{F}_{r+}</math> are affected by the relative speed <math>(v_{rel})_r</math> between the current elements in <math>r</math>, which leads to a kind of electromagnetic induction acting between co-linear current elements in relative motion in <math>r</math>. This feature may help to explain Cyril Smith's "Marinov Generator"<ref name="Marinov Generator"/>.
  
 
==A New Idea: The Makerarc==
 
==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:
+
The above result supports a development beyond the S.H.O. Drive. A much more powerful and compact permanent magnet system is now envisioned.
 +
 
 +
The preliminary name is Makerarc (MAY-KERR-ARK), which stands for:
  
 
* '''M'''agnetic
 
* '''M'''agnetic
Line 235: Line 565:
  
 
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)
 
Details pending. Stay tuned. Sincerely,  ''[[User:S.H.O.|S.H.O.]] <sup>[[User_talk:S.H.O.|talk]]</sup>'' 22:09, 5 March 2017 (PST)
 +
 +
What makes a magnet tick? Clues can be found in the video titled "MAGNETS: How Do They Work?" by Veritasium<ref>https://www.youtube.com/watch?v=hFAOXdXZ5TM</ref>. ''[[User:S.H.O.|S.H.O.]] <sup>[[User_talk:S.H.O.|talk]]</sup>'' 22:08, 23 March 2017 (PDT)
 +
 +
<youtube>https://www.youtube.com/watch?v=hFAOXdXZ5TM</youtube>
 +
 +
An attempt to produce the Makerarc back in March 2017<ref>https://www.facebook.com/Sho.Drives/videos/790278744462849/</ref> failed to demonstrate the predictions of the term <math>\mathbf{F}_r = \left( - Q_- \left(\mathbf{V}_d \cdot \nabla \right) \mathbf{A} \right)_r</math><ref>https://www.facebook.com/Sho.Drives/videos/791514291005961/</ref> where I assumed that an external vector potential <math>\mathbf{A}</math> would lead to a force co-linear to current elements. This may be due to the supercurrent nature of electrons bounded to atoms, allowing atoms to be capable of preserving or "fixing" the amount of magnetic flux that passes through them, just as "macroscopic" superconducting currents do for superconductors, and therefore, by extension be capable of resisting changes in the magnetic vector potential. This does not negate the possibility of a Marinov Generator, as designed by Cyril Smith, because in his case the currents receiving power were inside a conductor where an externally applied vector potential is not fully shielded against, permitting <math>\mathbf{F}_r = \left( - Q_- \left(\mathbf{V}_d \cdot \nabla \right) \mathbf{A} \right)_r</math> to yield a force. At the moment, I am still trying to decide on the correct model, but hopefully the surviving concept is similar to what I have developed so far in this page. Sincerely,  ''[[User:S.H.O.|S.H.O.]] <sup>[[User_talk:S.H.O.|talk]]</sup>'' 19:55, 14 July 2017 (PDT)
  
 
==Explaining "Altered" Lenz' Law Devices==
 
==Explaining "Altered" Lenz' Law Devices==
  
It is anticipated that the longitudinal force described in the previous section 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:
+
It is anticipated that the longitudinal force described above in the section titled "Explaining the Marinov Motor and Cyril Smith's 'Marinov Generator' using Conventional Physics" may explain some types of purported Reduced-Lenz devices. A good example can be found in the video below ("Ray's No Back EMF Generator"), although in this example, the longitudinal force ''increases'' the drag, mainly in positions where pancake generator coil is mostly outside the cylindrical boundary of the permanent magnet. This creates an illusion of a "Reduced" Lenz' Law effect when the magnet is mostly within the cylindrical boundary of the permanent magnet:
  
 
<youtube>x3TFALmHtMw</youtube>
 
<youtube>x3TFALmHtMw</youtube>
Line 245: Line 581:
  
 
==Explaining the Newman Motor==
 
==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.
 
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.

Latest revision as of 18:55, 14 July 2017

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

Introduction

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

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

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

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

Prior content in the "Comment Record" section:

Prior content in the "Background" section:

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

Novel Force laws proposed by various researchers

James Wesley's proposal

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

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

to the Lorentz force.

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

The extra term is equivalent to:[5]

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

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

The Lorentz force is:

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

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

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

The problem with this modification:

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

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

Cyril Smith's proposal

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

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

This is equal to:

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

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

The Lorentz force is:

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

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

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

The problem with this modification:

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Therefore:

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

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

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

Consider the existence of four charges:

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

Consider the existence of their corresponding current elements:

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

Where:

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

The current elements are equal to:

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

Where:

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

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

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

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

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

The corresponding relative speeds are:

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

The two drift velocities are:

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

Let the effective velocities of the current elements be:

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

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

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

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

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

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

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

The corresponding relative speeds are:

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Recalling that:

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

Substitution yields:

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

Simplified:

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

As stated earlier, the current elements are equal to:

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

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

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

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

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

Next, we make a substitution for the current elements:

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

This can be written as:

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

Or:

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

Or:

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

Adding the forces together, their sum is:

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

In S.I. Units:

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

Or:

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

The second differential is:

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

Special Scenario: No relative motion between positive charges

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

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

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

These become:

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

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

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

These become:

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

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

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

Therefore:

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

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

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

Therefore:

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

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

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

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

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

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

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

Based on the above results, we have:

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

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

In S.I. Units:

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

Or:

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

The second differential is:

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

From this and above, it follows that:

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

Therefore:

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

The effective electric fields experience by each charge is:

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

Given that:

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

It can be shown that:

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

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

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

A New Idea: The Makerarc

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

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

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

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

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

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

Explaining "Altered" Lenz' Law Devices

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

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

Explaining the Newman Motor

See also: Memory Lane

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

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

References

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

See also

Site map

HQGlossaryApril 2016 Presentation