Electro-Gravi-Magnetics (EGM) Practical Modelling Methods of The Polarizable Vacuum - VI
Electro-Gravi-Magnetics (EGM) Practical Modelling Methods of The Polarizable Vacuum - VI
Electro-Gravi-Magnetics (EGM) Practical Modelling Methods of The Polarizable Vacuum - VI
Abstract
A number of tools to facilitate the experimental design process are presented. These include
the development of a design matrix based on the unit amplitude spectrum, the derivation of Harmonic
and Spectral Similarity Equations, Critical Phase Variance, Critical Field Strengths and Critical
Frequency.
1
[email protected], [email protected].
1
1 ITRODUCTIO
1.1 GENERAL
In [1-10], a number of practical engineering tools for application to the Polarizable Vacuum
(PV) model of gravity was derived by application of Buckingham Π Theory (BPT) and dimensional
analysis techniques. BPT is a well-tested and experimentally verified method that relates a
mathematical model to an Experimental Prototype2 (EP).
The power of BPT to facilitate and articulate the derivation of the mathematical constructs
defined in [6-10] has advanced theoretical boundaries in particle physics to a higher level. The tools
derived have volunteered precise calculations leading to the determination of the mass-energy and radii
of Leptons, Quarks and Bosons. Most notably, the radius prediction for the Proton has been
experimentally verified by the SELEX Collaboration to astounding precision. [7]
Two important aspects of the works derived are Critical Ratio “KR” and General Similarity
Equations (GSE). “KR” is defined as the ratio of the sum of the magnitudes of the applied
electromagnetic (EM) fields, to the magnitude of the background gravitational field. The General
Modelling Equations (GME) exploit this definition to produce the “GSE” derived in [5].
1.3 OBJECTVES
The formulations derived herein assist in the qualitative and quantitative experimental design
process as follows,
i. Harmonic representation of |KR| = 1 in the Fourier domain over an elemental displacement “∆r”,
termed the Critical Harmonic Operator “KR H”.
ii. Utilisation of “KR H” to formulate harmonic representations of various other physical variables
for consideration in the experimental design process.
iii. Utilisation of “KR H” to simplify “GSEx”, on a modal basis, to Harmonic Similarity Equations
(HSEx).
iv. Graphical visualisation of “HSEx” based on Complex Phasor Forms of the magnitude of the
applied Electric Field Intensity “EA” and Magnetic Flux Density “BA”.
v. Reduction of “HSEx” into simplified electromagnetic (EM) design consideration forms “HSEx R”.
vi. Spectral Similarity Equations (SSE): these qualify and quantify the similarity of a singularly
applied experimental EM source to the frequencies that inhabit the ambient Electro-Gravi-
Magnetic (EGM) spectrum. [4]
vii. Determination of the applied EM phase requirements with respect to the background
gravitational field utilising “SSEx”.
viii. Assess suitability of Maxwells Equations (ME) to experimental investigations utilising “SSEx”.
1.4 RESULTS
2
The experimental prototype represents the PV at the surface of the Earth by which all similarity
conditions are referenced.
2
2 THEORETICAL MODELLIG3
Assuming complete dynamic, kinematic and geometric similarity between the EP and the
mathematical model (|KR| =1) where the harmonic mode of the PV “nPV” approaches the harmonic cut-
off mode “nΩ” [|nPV| → nΩ and |nΩ| < ∞], “KR” has many representations. One such representation
incorporating the change in harmonic frequency modes across “∆r” shall be derived.
The spectral characteristics of the EP may be articulated at the surface of the Earth assuming4,
i. The Zero-Point-Field (ZPF) physically exists as a spectrum of frequencies and wave vectors.
ii. The sum of all ZPF wave vectors at the surface of the Earth is coplanar with the gravitational
acceleration vector. This represents the only vector of practical experimental consequence.
iii. A modified Complex Fourier Series representation of “g” is physically real and is representative
of the magnitude of the resultant ZPF wave vector.
iv. A physical relationship exists between gravity, electricity and magnetism such that the physical
interaction of applied EM fields with the PV, in accordance with the hypotheses to be tested [2]
may be investigated and potentially modified.
It was illustrated in [3] that, for an engineered change in “g” by application of BPT and the
Equivalence Principle, a change in the PV may be described [as |nPV| → nΩ] by,
2
E PV k PV, n PV, t
∆K 0( ω , X ) n PV, k PV
∆g ≡ ∆a PV .
r B PV k PV, n PV, t
2
n PV, k PV
i (1)
Where,
Variable Description Units
∆g Change of gravitational acceleration vector m/s2
∆aPV Change in PV acceleration vector
EPV(kPV,nPV,t) Magnitude of PV electric field vector V/m
BPV(kPV,nPV,t) Magnitude of PV magnetic field vector T
ω Field frequency Hz
kPV Harmonic wave vector of PV m-1
i Denotes initial conditions of PV None
∆K0(ω,X) Engineered relationship function
r Magnitude of position vector from centre of mass m
Table 1,
Subsequently, considering only the resultant ZPF wave vector relating to “g” in a practical laboratory
experiment, equation (1) may be usefully simplified by removing “kPV” notation and relating it to a
generalised Fourier representation5 of constant “g” over “∆r” as |nPV| → nΩ by,
2
E PV n PV, t
∆ K 0( ω , X ) n PV G. M . 2 π .n PV .∆ω δr( 1 , r , ∆r , M ) .t .i
. i . .e
r B
2 2 π . n PV
PV n PV, t r n PV
n PV
i (2)
Where, “∆ωδr” denotes the beat frequency across “∆r” as defined in [5].
3
All mathematical modelling and output was formed using “MathCad 8 Professional” and appears in
standard product notation.
4
A spherical object with uniform mass distribution.
5
The Right Hand Side (RHS) of equation (2) is analogous to the form utilised in [4].
3
2
E PV n PV, t
2
n PV
c
2
B PV n PV, t
n PV
(3)
G.M .
∆K 0( ω , X ) KR
2
r .c (4)6
Substituting equation (3) and (4) into (2) yields the PV EM harmonic representation of the ideal value
of the magnitude of “KR” for the complete reduction of “g” over “∆r” in a laboratory at the surface of
the Earth7 as |nPV| → nΩ,
2. 1 . π .n PV .∆ω δr( 1 , r , ∆r , M ) .t .i
K R( r , ∆ r , M ) i . e
π n PV
n PV
(5)
Where the maximum amplitude occurs at time index,
1
t n PV, r , ∆r , M
2 . n PV .∆ω δr( 1 , r , ∆r , M )
(6)
Yields the unit amplitude spectrum analogous to the result previously found [4] as |nPV| → nΩ(r+∆r),
2
K R n PV
H .
π n PV
(7)
3 MATHEMATICAL MODELLIG
Utilising equation (7) a table of expressions for the magnitude of the amplitude spectrum of
various experimental design considerate relationships may be formulated for complete dynamic,
kinematic and geometric similarity between the EP and the mathematical model (|KR| =1) as |nPV| →
nΩ ,
Description Primitive Form Harmonic Form Result
Engineered G.M . G.M . Σ∆K0 H
∆K 0( ω , X ) KR ∆K 0 n PV, r , M , X K R n PV
Relationship 2 H
r .c
2 H →∆K0
r .c
Function
Refractive Index G .M K PV n PV, r , M K PV( r , M ) .K R n PV ΣKPV H
2.
2 H H
r .c →KPV
K PV e
Engineered 2 . ∆K 0( ω , X ) K EGM n PV, r , M , K R K EGM r , M , K R .K R n PV ΣKEGM H
K EGM K PV. e H H
Refractive Index →KEGM
Gravitational G. M . m U g n PV, r , M U g ( r , M ) .K R n PV ΣUg H
H H
Potential Energy r G. M →Ug
Ug
(per unit mass) m r
Critical Factor8 ∆K C( r , ∆r , M ) ε0 µ0 Σ∆KC H
. ∆K C n PV, r , ∆r , M ∆U PV( r , ∆r , M ) .K R n PV .
KR H H ε0 →∆KC
∆U PV( r , ∆r , M ) µ0
Table 1,
6
“∆K0(ω,X)” is defined as the engineered relationship function. [3]
7
Where, “i” on the RHS of equation (2) represents complex number notation.
8
Representing the applied experimental field. [2,3]
4
Where, the permittivity and permeability of free space (“ε0” and “µ0” respectively) act as the
impedance function9 such that,
Variable Description Units
∆UPV(r,∆r,M) Change in energy density of PV Pa
∆KC(r,∆r,M) Change in critical factor PaΩ
G Gravitational constant m3kg-1s-2
Table 2,
4 PHYSICAL MODELLIG
A family of “HSEx” may be defined by relating the EP to the mathematical model on a modal
basis, termed discrete similarity for |∆r| << ∞. Utilising “GSE1,2” derived in [5] yields “HSE1,2” as
follows,
2 2 2 π .n PV .∆ω δr( 1 , r , ∆r , M ) .t .i
i . EA k A, n A, t ± c .B A k A , n A , t .e
(8)
HSE E A , B A , k A , n A , n PV, r , ∆r , M , t
1, 2 2 2
π . n PV. c . B A k A , n A , t
9
The impedance function “ Z = µ 0 /ε 0 ” is independent of KPV in the PV representation.
5
Where, [5]
µ0
St α ( r , ∆r , M ) ∆U PV( r , ∆r , M ) .
ε0
(15)
2
3.M .c . 1 1
∆ U PV( r , ∆ r , M )
4.π (r ∆ r)
3 3
r (16)
Variable Description Units
EA(kA,nA,t) Magnitude of applied electric field vector in complex form V/m
BA(kA,nA,t) Magnitude of applied magnetic field vector in complex form T
c Velocity of light in a vacuum m/s
M Mass kg
Table 3,
2.
B rms B A B 0 , n B, φ , r , ∆ r , M , t B0
2 (20)
Variable Description Units
EA(E0,nE,r,∆r,M,t) Applied electric field vector V/m
BA(B0,nB,φ,r,∆r,M,t) Applied magnetic field vector T
E0 Amplitude of electric field vector V/m
B0 Amplitude of magnetic field vector T
nE Harmonic mode number of the ZPF with respect to “EA” None
nB Harmonic mode number of the ZPF with respect to “BA” None
φ Relative phase variance between “EA” and “BA” None
ωE Harmonic frequency of the ZPF with respect to “EA” Hz
ωB Harmonic frequency of the ZPF with respect to “BA” Hz
Erms Root-mean-square of “EA” V/m
Brms Root-mean-square of “BA” T
Table 4,
10
The expression inside the magnitude notation on the RHS of equation (8-10).
6
Equations (17,18) are functions in Complex Form and contain Real and Imaginary
components. For visualisation purposes, only the Real component is required. Figure (1) includes a
graphical representation of “EA” and “BA” for arbitrary illustrational values.
The representations for “Re(HSE1,2)” have been accentuated for illustrational purposes by a
large value of “φ” (180°). Typically, values of “0°” would be expected in accordance with classical EM
propagation, or “90°” in accordance with Maxwell’s Equations11.
V
Re E A 1 . , 1 , R E , ∆r , M E , t
m
Re B A 1 .( T ) , 1 , 180 .( deg ) , R E , ∆r , M E , t
V
Re HSE 1 1 . , 1 .( T ) , 1 , 1 , 180 .( deg ) , 3 , R E , ∆r , M E , t
m
V
Re HSE 2 1 . , 1 .( T ) , 1 , 1 , 180 .( deg ) , 3 , R E , ∆r , M E , t
m
t
Electric Forcing Function
Magnetic Forcing Function
HSE 1
HSE 2
Figure 1,
Figure (2) includes arbitrary illustrational values but also contains important information
regarding “φ”. Exploratory graphical analysis demonstrates that “Re(HSE3)” is in-phase with
“Re(HSE4)” and out-of-phase with “Re(HSE5)” for key values (0° and 90°) of “φ”. The significance of
this being that “Re(HSE3)” is analogous to the Poynting Vector and implies that “Re(HSE4)” is
representative of constructive EGM interference and “Re(HSE5)” is representative of destructive EGM
interference.
“HSE4,5” were formed from General Modelling Equation “1” and “2” (GME1,2) as described in
[5]. “GME1” is proportional to a solution of the Poisson equation applied to Newtonian gravity, where
the resulting acceleration is a function of the geometry of the energy densities. “GME2” is proportional
to a solution of the Lagrange equation where the resulting acceleration is a function of the Lagrangian
densities of the EM field harmonics in a vacuum.
Therefore, experimental investigations with the objective of reducing the local gravitational
acceleration on a test bench, by means of EGM interference, should bias engineering designs governed
by “HSE5”. However, designs favouring “HSE4” should not be completely discounted and should form
part of any complete design process.
V
Re HSE 3 1 . , 1 .( T ) , 1 , 1 , 0 .( deg ) , 3 , R E , ∆r , M E , t
m
V
Re HSE 4 1 . , 1 .( T ) , 1 , 1 , 0 .( deg ) , 3 , R E , ∆r , M E , t
m
V
Re HSE 5 1 . , 1 .( T ) , 1 , 1 , 0 .( deg ) , 3 , R E , ∆r , M E , t
m
t
HSE 3
HSE 4
HSE 5
Figure 2,
11
This shall be discussed in greater detail in proceeding sections.
7
4.3 REDUCTION OF HSEx
2 . ( cos ( 2 . φ ) 1)
HSE 1 φ , n PV
R π . n PV
(21)
2 . ( cos ( 2 . φ ) 1)
HSE 2 φ , n PV
R π . n PV
(22)
K PV( r , M ) . St α ( r , ∆ r , M )
HSE 3 E rms , B rms , n PV, r , ∆ r , M
R π . n PV. E rms . B rms
(23)
12
See appendix (A).
13
Of the form “HSEx R”.
14
Alteration in notation is required for graphing purposes. It is a limitation of “MathCad 8
Professional” that axial arguments may not be written precisely in the form “HSE1 R(φ,1)” etc.
8
π
π
2
HSE 1_R ( φ , 1 )
Harmonic Similarity
HSE 1_R ( φ , 2 )
HSE 2_R ( φ , 1 )
HSE 2_R ( φ , 2 )
φ
Phase Variance
Figure 3,
Figure (4) analysis [Y-Axis is logarithmic]:
iv. |HSE3 R| → 1 as |nPV| → nΩ ZPF: This is consistent with Poynting Vector characteristics described
in [2,3].
Where,
Variable Description Units
nΩ(r,∆r,M)ZPF ZPF beat cut-off mode across “∆r” at “r” [5] None
RE Radius of the Earth m
ME Mass of the Earth kg
Table 5,
Harmonic Similarity
n PV
Harmonic Mode
Figure 4,
Reducing the number of independent variables where possible can make obvious gains. To
achieve this, we shall express “HSE4,5 R” as time averaged functions “HSE4A,5A R” explicitly in terms of
“φ” and “HSE3 R” as follows15,
1 . HSE
HSE 4A E rms , B rms , φ , n PV, r , ∆ r , M 3 E rms , B rms , n PV, r , ∆ r , M
R cos ( φ ) R
(26)
1 . HSE
HSE 5A E rms , B rms , φ , n PV, r , ∆ r , M 3 E rms , B rms , n PV, r , ∆ r , M
R sin ( φ ) R
(27)
Hence,
v. |HSE4,5 R| → 1 at “φ = 0°” and “φ = 90°” respectively.
vi. |HSE4| = |HSE5| at “φ = 45°” and “φ = 135°”.
15
Refer to appendix (A) for derivation.
9
4.6 SPECTRAL SIMILARITY EQUATIONS
The preceding sections define the requirements for complete dynamic, kinematic and
geometric similarity with any specific mode in the background EGM field. However, reproduction of
only one specific mode for experimental investigations is extremely limiting. Alternatively, it is highly
advantageous to consider the reproduction of a harmonically averaged distribution for each “HSE”,
termed Spectral Similarity Equations (SSE).
“SSE” are defined as a family of equations that quantify and qualify the similarity of a single
field source defined by “HSE” with respect to the spectrum of frequencies that inhabit the background
EGM field. “SSE” differs from “GSE” in that “GSE” represents similarity of multiple EM sources with
respect to the background EGM field.
Therefore, utilising “HSE” above, the magnitude of the average spectral similarity per
frequency mode with respect to the applied forcing function may be generalised as follows
1 .
SSE HSE
n Ω ( r, ∆ r, M ) 1
ZPF n PV
(28)
Where, “nPV” has the odd harmonic distribution: -nΩ ZPF, 2 - nΩ ZPF …. nΩ ZPF
Recognising that16,
ln 2 . n Ω ( r , ∆ r , M ) γ
1 . 1 ZPF
n Ω ( r, ∆ r, M ) 1 n PV n Ω ( r, ∆ r, M ) 1
ZPF n PV ZPF
(29)
As nPV → nΩ ZPF and nΩ ZPF >>1: Substituting equation (29) into the “HSE” yields,
ln 2 . n Ω ( r , ∆ r , M ) γ
2 . ( cos ( 2 . φ ) 1) . ZPF
SSE 1( φ , r , ∆ r , M )
π n Ω ( r, ∆ r, M ) 1
ZPF (30)
ln 2 . n Ω ( r , ∆ r , M ) γ
2 . ( cos ( 2 . φ ) 1) . ZPF
SSE 2( φ , r , ∆ r , M )
π n Ω ( r, ∆ r, M ) 1
ZPF (31)
.
K PV( r , M ) . St α ( r , ∆ r , M ) ln 2 n Ω ( r , ∆ r , M ) ZPF γ
SSE 3 E rms , B rms , r , ∆ r , M .
π . E rms . B rms n Ω ( r, ∆ r, M ) 1
ZPF (32)
1 . SSE
SSE 4 E rms , B rms , φ , r , ∆ r , M 3 E rms , B rms , r , ∆ r , M
cos ( φ ) (33)
1 . SSE
SSE 5 E rms , B rms , φ , r , ∆ r , M 3 E rms , B rms , r , ∆ r , M
sin ( φ ) (34)
“φC” is defined as the phase difference between “EA” and “BA” for complete dynamic,
kinematic and geometric similarity with the background EGM field “|SSEx| = 1”. Therefore, by
analyses of the preceding figures or the appropriate transformation of equation (31-34), “φC” may be
easily determined.
For proportional solutions to the Poisson equation applied to Newtonian gravity where the
resulting acceleration is a function of the geometry of the energy densities, “φC = 0°”. For proportional
solutions to the Lagrange equation where the resulting acceleration is a function of the Lagrangian
densities, “φC = 90°”.
16
Error “< 6.7x10-6(%)” at “nΩ ZPF > 106”.
10
4.8 CRITICAL FIELD STRENGTH
“EC” and “BC” are derived utilising the reciprocal harmonic distribution describing the EGM
amplitude spectrum. Solutions to “|SSE4,5| = 1” represent conditions of complete dynamic, kinematic
and geometric similarity with the amplitude of the background EGM spectrum. “EC” and “BC” denote
RMS values satisfying the proceeding equation,
π
SSE 4 E rms , B rms , 0 , r , ∆ r , M SSE 5 E rms , B rms , , r, ∆ r, M 1
2 (35)
Where,
E rms E C( r , ∆ r , M )
(36)
E C( r , ∆ r , M )
B rms
c (37)
4.9 DC-OFFESTS
The value of “EC” and “BC” may be decreased by the application of an offset function “DC”.
This denotes a percentage offset of the forcing function and may be applied to facilitate a specific
experimental configuration. For example, if “DC = 100(%)” the value of “EC” and “BC” computed
above yield,
π 1
SSE 4 ( 1 DC) . E rms, B rms, 0 , r , ∆ r , M SSE 5 E rms, ( 1 DC) . B rms, , r, ∆ r, M
2 2 (38)
π 1
SSE 4 ( 1 DC) . E rms, ( 1 DC) . B rms, 0 , r , ∆ r , M SSE 5 ( 1 DC) . E rms, ( 1 DC) . B rms, , r, ∆ r, M
2 4 (39)
Therefore, by re-computing the value of “EC” and “BC” at “|SSE4,5| = 1” a decrease in Critical
Amplitude shall be observed.
5 MAXWELLS EQUATIOS
5.1 GENERAL
By considering “ME” in relation to the applied EM fields and the requirements of similarity, it
is possible to deduce important design characteristics for further consideration. “ME” (in MKS units)
for time-varying fields are as follows (“ρ” is the charge density and “J” is the vector current density),
[11]
ρ ∂B . ∂E
∇ .E ,∇ E , ∇ B 0, ∇ B µ 0.J ε 0.µ 0.
ε 0 ∂t ∂t
(40)
Consequently as |SSE5| → 1, optimal similarity occurs when:
1. The divergence of “EA” is maximised.
2. The curl of “BA” is maximised.
“ωC” is defined as a half wavelength over “∆r” by applied fields and represents a minimum
frequency threshold for the application of “ME” within this experimental context,
c
ω C( ∆ r )
2.∆ r (41)
11
6 COCLUSIOS
A number of tools that facilitate the experimental design process are presented. These include
the development of a design matrix based on the unit amplitude spectrum, the derivation of Harmonic
and Spectral Similarity Equations, Critical Phase Variance, Critical Field Strengths and Critical
Frequency.
References
12
APPEDIX A
π .
2 .π .n E .t i
2
substitute , E A E 0 .e
π
2 .π .n B .t φ .i
2
substitute , B A B 0 .e
2 2. 2
i . EA c BA i
substitute , n E n B . ( exp( 2 . i .φ ) 1)
2
π . n PV. c . B A
2 π . n PV
substitute , B 0 2 . B rms
substitute , E 0 2 . E rms
simplify (A1)
1
2
i . ( exp( 2 . i . φ ) 1. ( cos ( 2 . φ ) 1)
1 ) simplify 2.
π . n PV π n PV
2
(A2)
2 . ( cos ( 2 . φ ) 1)
HSE 1 φ , n PV
R π . n PV
(A3)
π .
2 .π .n E .t i
2
substitute , E A E 0 .e
π
2 .π .n B .t φ .i
2
substitute , B A B 0 .e
2 2 2
i . EA c .B A i
substitute , n E n B . ( exp( 2 . i .φ ) 1)
π . n PV. c B A
2 2. π . n PV
substitute , B 0 2 . B rms
substitute , E 0 2 . E rms
simplify (A4)
1
2
( exp( 2 .i . φ ) 1) 1. ( 2 . cos ( 2 . φ ) 2)
i . simplify
π . n PV π n PV
2
(A5)
2 . ( cos ( 2. φ ) 1)
HSE 2 φ , n PV
R π . n PV
(A6)
π .
2 .π .n E .t i
2
substitute , E A E 0 .e
π
2 .π .n B .t φ .i
2
substitute , B A B 0 .e
2 . i . K PV( r , M ) . St α ( r , ∆ r , M ) St α ( r , ∆ r , M )
substitute , n E n B i . K PV( r , M ) . . exp i . 4 . π . n .
EM t π φ
π . n PV. E A . B A π . n PV. E rms . B rms
substitute , n B n EM
substitute , B 0 2 . B rms
substitute , E 0 2 . E rms
simplify (A7)
St α ( r , ∆ r , M ) 1. St α ( r , ∆ r , M )
i . K PV( r , M ) . . exp i . 4. π . n .
EM t π φ simplify K PV( r , M ) .
π . n PV. E rms. B rms π n PV. E rms. B rms
(A8)
K PV( r , M ) . St α ( r , ∆ r , M )
HSE 3 E rms , B rms , n PV, r , ∆ r , M
R π . n PV. E rms . B rms
(A9)
π .
2 .π .n E .t i
2
substitute , E A E 0 .e
π
2 .π .n B .t φ .i
2
substitute , B A B 0 .e
2 substitute , n E n B
4 .i . St α ( r , ∆ r , M ) . K PV( r , M ) . c . B A K PV( r , M )
2 . i . St α ( r , ∆ r , M ) . . exp( i . φ )
2 2 2 substitute , n B n EM E rms . B rms .π . n PV. exp i . π . 4 .n EM . t 1 exp i . 4 . π . n EM . t π 2.φ
π . n PV. E A . E A c .B A
substitute , B 0 2 . B rms
substitute , E 0 2 . E rms
E rms
substitute , c
B rms
13
2 .HSE 3 E rms , B rms , n PV, r , ∆ r , M
R
HSE 4 E rms , B rms , n EM , φ , n PV, r , ∆ r , M , t
R
cos 4 . π . n EM .t . cos 4 . π . n EM . t 2 . φ ...
+ sin 4 . π . n EM . t .sin 4 .π . n EM . t 2 . φ ...
+1 (A12)
π .
2 .π .n E .t i
.e 2
substitute , E A E0
π
2 .π .n B .t φ .i
2
substitute , B A B 0 .e
2 substitute , n E n B
4 .i . St α ( r , ∆ r , M ) . K PV( r , M ) . c . B A exp( i . φ )
2 .i . St α ( r , ∆ r , M ) . K PV( r , M ) .
2 2 2 substitute , n B n EM E rms . B rms . π . n PV. exp i .π . 4 . n EM . t 1 exp i . 4 .π . n EM . t π 2.φ
π . n PV.E A . E A c .B A
substitute , B 0 2 . B rms
substitute , E 0 2 . E rms
E rms
substitute , c
B rms
i .φ 2 2
2.i . e . K PV( r , M ) . St α ( r , ∆ r , M ) 1. 2
K PV( r , M )
simplify 2 . St α ( r , ∆ r , M ) .
i .π . 4 .n EM .t 1 i . 4 .π .n EM .t π 2 .φ π 2 2 2
B rms . E rms . n PV . cos 4 . π . n EM . t . cos 4 . π . n EM . t 2. φ 1 sin 4 . π . n EM . t . sin 4 . π . n EM . t 2.φ
E rms . B rms . π . n PV. e e
(A14)
2 . HSE 3 E rms , B rms , n PV, r , ∆ r , M
R
HSE 5 E rms , B rms , n EM , φ , n PV, r , ∆ r , M , t
R
cos 4 . π . n EM .t . cos 4 . π . n EM . t 2 . φ ...
+ sin 4 . π . n EM . t .sin 4 .π . n EM . t 2 . φ ...
+ 1 (A15)
T δr n PV, r , ∆ r , M
HSE 4A E rms, B rms, φ , n PV, r , ∆ r , M ∆ω δr n PV, r , ∆ r , M . HSE 4 E rms, B rms, n EM , φ , n PV, r , ∆ r , M , t dt
R R
0.( s ) (A16)
T δr n PV, r , ∆ r , M
2 . HSE 3 E rms , B rms , n PV, r , ∆ r , M T n , r, ∆ r, M
∆ω δr n PV, r , ∆ r , M . R
dt ∆ω δr n PV, r , ∆ r , M . δr PV . HSE
1
3 E rms , B rms , n PV, r , ∆ r , M R
cos 4 . π . n EM . t . cos 4 . π . n EM . t 2 . φ ... 2
2
+ sin 4 . π . n EM . t . sin 4 . π . n EM . t 2 . φ ... cos ( φ )
+1
0.( s ) (A17)
Recognising,
∆ω δr n PV, r , ∆ r , M .T
δr n PV, r , ∆ r , M 1
(A18)
Simplifying the integral yields the Time Averaged HSE4 at each frequency mode “nPV” as follows:
1 .HSE
HSE 4A E rms, B rms, φ , n PV, r , ∆ r , M 3 E rms, B rms, n PV, r , ∆ r , M
R cos ( φ ) R
(A19)
Hence, it follows that:
T δr n PV, r , ∆ r , M
HSE 5A E rms , B rms , φ , n PV, r , ∆ r , M ∆ω δr n PV, r , ∆ r , M . HSE 5 E rms , B rms , n EM , φ , n PV, r , ∆ r , M , t dt
R R
0.( s ) (A20)
T δr n PV, r , ∆ r , M
K PV( r , M ) . St α ( r , ∆ r , M )
2.
π . n PV. E rms . B rms i
∆ω δr n PV, r , ∆ r , M . dt . HSE
3 E rms , B rms , n PV, r , ∆ r , M R
cos 4 . π . n EM . t . cos 4 . π . n EM . t 2 . φ ... sin ( φ )
2
14