Equations are not being displayed properly on some articles. We hope to have this fixed soon. Our apologies.

Turtur, C. (2011). Construction guidelines for a ZPE-Converter on the basis of realistic DFEM-Computations. PHILICA.COM Article number 233.

ISSN 1751-3030  
Log in  
  997 Articles and Observations available | Content last updated 19 February, 15:50  
Philica entries accessed 2 793 311 times  

NEWS: The SOAP Project, in collaboration with CERN, are conducting a survey on open-access publishing. Please take a moment to give them your views

Submit an Article or Observation

We aim to suit all browsers, but recommend Firefox particularly:

Construction guidelines for a ZPE-Converter on the basis of realistic DFEM-Computations

Claus W. Turturconfirmed userThis person has donated to Philica (Fachbereich Elektrotechnik, University of Applied Sciences Braunschweig-Wolfenbuettel)

Published in physic.philica.com

In [Tur 11] the theory of a powerful vacuum-energy converter was developed, and such converters have been simulated with a dynamic finite element method (DFEM). The result was a theoretical description of the machine which should be appropriate for technical applications.
Due to many questions from colleagues who read the mentioned article, the author decided to continue his development on the DFEM-algorithm in order to simulate a zero-point-energy (ZPE) motor on the computer, as close to reality as possible.
The theoretical background of the simulation is explained in detail here, so that every colleague should be able, to use the algorithm in the appendix of the publication and to adapted it to the setup of a vacuum-energy motor according to his own conception.

Article body


Table of contents
1. Physical fundament and preliminary work
2. Motion of the components of the ZPE-Converter
3. Evaluation of the results of a converter example
4. Computation example for a concrete ZPE-motor
5. A concret EMDR vacuum energy converter
6. The EMDR-Converter with mechanical power-extraction
7. Practical advice for experimenalists, who want to build an EMDR-Converter
8. Resumée
9. Literature References
10. Appendix: Source-Code of the DFEM-Algorithm

1. Physical fundament and preliminary work

The algorithm is designed to simulate electric and magnetic ZPE-motors by principle, and it is not restricted to one special design or setup. Thus it allows simulating ZPE-motors with arbitrary position and numbers of coils as well as arbitrary positions and numbers of magnets. Even electrostatic ZPE-motors can be simulated. Also interaction with external entities can be included into the simulation, such as connections of coils with electrical circuits.

The functioning principle of the DFEM-algorithm is the following: Any motion (of mechanical components as well as of electrical charge or even of electrical and magnetical fields) is being brought back to differential equations, which also contain the information of external electrical circuits. Permanent magnets have to be simulated by conductor loops containing electrical current without power supply (which is not an unusual point of view). This allows the algorithm to determine Lorentz-forces, with which external magnetic fields interact with permanent magnets.

In order to present a concrete result (for the suggestion of a prototype), the algorithm which is shown in the appendix, is designed to simulate two coils and one permanent magnet as being drawn in figure 1.


Fig.1: In a three-dimensional Cartesian coordinate system (blue colour), we see two coils orient­ed parallel to the yz-plain (red colour). The corners of the coils are located at the Cartesian coordinates as written in red colour. A cylindrical permanent magnet is being simulated by two conductor loops, one at its top end and the other one at its bottom end. The magnet can rotate around the z-axis.The number and the arrangement of the coils can be chosen arbitrarily in the DFEM-algorithm. The example as being shown here, corresponds to the source-code in the appendix of the publication.


In order to prepare the solution of the differential-equations within the DFEM-algorithm, we need the following:

(a.)       The computation of the induced voltage, which the rotating permanent magnet induces in the coils,


(b.)      The computation of the magnetic force Lorentz-force with which the coils act onto the permanent magnet.

These both computations have the purpose to realise the coupling of the mechanical and the electrical parts of the system to each other (see [Tur 11]). On the one hand, the mechanical rotation of the magnet is influenced by the electrical current in the coils due to the Lorentz-force (between these currents and the permanent magnet), and on the other hand the electrical current in the LC-oscillation-circuit is influenced by the induced voltage which the rotation of the permanent magnet brings into the coils.

We now want to turn our attention to those both calculations as given under (a.) and (b.):


Details of a:

The determination of the voltage induced in the coils (due to the motion of the magnet) is based on the time dependent alteration of the magnetic flux , which has its reason in the rotation of the permanent magnet (in our example around the z-axis), which is normally calculated as being seen in equation (1). [Jac 81]

This computation begins with a determination of the magnetic field of the permanent magnet as a vector field which has to be stored in a data-array. This can be done experimentally or theoretically, as for instance with one of the subroutines „Magnetfeld_zuweisen“, which are marked with different numbers in the source-code in order to allow all the emulation of different permanent magnets. The vector field is now fixed rigidly to the permanent magnets so that each of permanent magnets has its own vector field. As soon as we rotate a magnet, the field is rotating together with its magnet. The rotation around the z-axis is realised with a coordinate-transformation as usual (see equation 2). The angle φ describes the orientation of the length-symmetry-axis of the magnet relatively to the y-axis. (x,y) are the coordinates in the system without rotation (as shown in figure 1), and (x’,y’) are the coordinates in the system being rotated relatively to (x,y) by an angle of φ. The responsible subroutine in the algorithm has the name „Magnet_drehen“.

By this means, we calculate the magnetic field strength, which the magnet produces at any arbitrary position in the space, namely as a function of the angle φ. Responsible for the computation of the field strength is the subroutine with the name „Feldstaerke_am_Ort_suchen“. Furthermore, at the end of this subroutine, the magnetic flux through the coil is being determined, in which the induced voltage has to be calculated. Therefore we use a subdivision of the coils into finite area-element, so that the magnetic flux can be taken into account as a function of the position, where it passes the coils.

In our very example, this computation is being simplified by the fact, that all orthogonal-vectors on all area elements of each coil are orientated exactly into the direction of the x-axis, so that the scalar-product of the field with orthogonal-vectors of the coil-area-elements can be derived as simple as shown in equation (3), namely as the x-component of the flux.

The magnetic flux through the whole coil is then being calculated as the sum of all magnetic flux-elements through all finite area-elements forming the coil, so that the total flux follows equation (4).

Remark regarding the area-elements of the coil (index „SFE):

In order to formulate the possibilities for the variation of the geometry of the coils as flexible as possible, each coil has to be described as a polygonal line, connecting arbitrarily defined support points. Each coil is modulated connecting its support points with each other. This allows the definition of arbitrarily shaped areas to be surrounded by coils. Finite conductor-loop elements are defined by the geometrical connections between support point and support point. Finite area-elements of the coils fill up the area surrounded by the conductor-loops, so that the magnetic flux through each coil can be calculated as the sum of the finite magnetic flux elements through all area-elements of the individual coil.


From this result as described in equation (4), the calculation of the induced voltage is simply a derivation to time as shown in equation (1). Therefore we need the time-dependency of the orientation of the magnet (given by the angle φ), this is the angular velocity of the permanent magnet at each moment of time. If we check our calculation with a constant angular velocity of the magnet, we come to a result as shown in figure (2).



Fig.2: For the verification of the com­putation-method, a permanent magnet producing a homogeneous magnetic field was rotated with constant angular velocity, and the induced voltage in both red coils of figure 1 was plotted. The fact, that both voltages differ by a factor of 2 (from coil to coil) has its reason in the fact, that one coil has twice as many windings as the other one.

For additional computations (see b.), we need to determine the torque which the magnet experiences due to the electric current in the coils. This requirement makes it necessary to emulate the permanent magnet by a configuration of conductor loops. If the magnetic field is not just a simple homogeneous field (on which fig.2 is based), the computation of the magnetic flux depends on the spatial resolution of the computation of the magnetic field. Due to this reason, the magnetic flux of more complicated magnets always display some numerical noise (due to the fact, that the finite elements are not continuous but disctere in spatial resolution). And the problem is, that the numerical noise is enhanced remarkably, when we calculate the derivative to time. An example for such numerical noise is shown in figure 3, as it was calculated by the subroutine „Magnetfeld_zuweisen_02“. When we will emulate a real cylindrical bar-magnet later (with the subroutine “Magnetfeld_zuweisen_03”), the numerical noise will be even much worse.

Fig.3: Although the magnetic flux from a magnet rotating with constant angular-velocity has only moderate numerical noise, the computation of the induced voltage needs numerical smoothing urgently. Otherwise it would not be sufficient for any use in the DFEM-algorithm.

In order to smooth the numerical noise of the voltage-signal, a Fourier-series was developed (in the subroutine “Fourier_Entwicklung”). It is important to take only low order components (maximum up to fifth order), in order to assure, that high frequency components are excluded. The less high-order components we take, the smoother the signal we get.

An additional effect of the Fourier-series is, that it helps to save CPU-time remarkably, when we will have to calculate the magnetic flux very often for the solution of the differential equation later. The computation of the magnetic flux itself contains a sum (see equation 4), which contains time-consuming operations (see equations 2 and 3), which take much more CPU-time, then the calculation of the Fourier-series with only five rather simple expressions to be summed up. This is important, because the solution of the differential-equation will have to be done with very many time-steps of few nanoseconds, so that the computation of the magnetic flux has to be done several 108 or 109 during each run of the DFEM-algorithm. This forces us to speed up the very innerst loops in the program, as the computation of the magnetic flux is one of them.

The fact that the Fourier-series approximation is done in rather low order (N0≤5) allows us to determine the Fourier-coefficients very easy by the use of the Gauß’ian method of the least square fit, to compare the Fourier approximation with the original data. This method is used to determine the Fourier-coefficients Aν.


Details of b:

For the determination of the Lorentz-forces, by which the electrical currents in the coils accelerate the permanent magnets, the emulation of the permanent magnets have being realised (as mentioned above) by conductor loops containing electrical current. In our example for the demonstration, we want to emulate cylindrical bar-magnets, because they are easy to buy and not very expensive, with regard to experiments. For the sake of simplicity, the cylindrical bar-magnets are emulated by two circular conductor loops, located at each end of the cylindrical bar. This means, the magnetic field of the bar-magnets is emulated by the magnetic field of one pair of circular coils.

The parameters we need are only the length of the bar, the diameter of the bar and the magnetic field at each end of the bar. Especially the magnetic field-strength can be adjusted by the choice of the electric current in the conductor loops emulating the bar-magnet.

The calculation of the magnetic field of the emulation-conductor-loops is done with the use of Biot-Savart’s law. The approach is illustrated in figure 4.

Fig.4: Illustration of the geometry of a conductor loop (green) whose elements are parametrized by a position vector l.  According to Biot-Savart, we calculate the magnetic field, which the conductor loop produces at any arbitrary point of interest  s. The conductor loop shown here describes the top end of the cylindrical bar-magnet, which is orientated along the y-axis. Thus the loop is orientated parallel to the xz-plane, has the radius „MEro“, and is located at the y-position „MEyo“.

The parameterisation of the conductor loop can be realised rather simple according to equation 6.

In principle, we now can introduce the expressions of (8) and (9) into the outer product of (7), but in order to make (6) complete, we additionally need the electrical charge q1 in (8). This has to be determined from the current I in the coil and the propagation-velocity of the electrical charge, as being described by the angular velocity ω. Therefore I  and ω have to be combined in such a way, that the motion of the electrical charge is being described appropriately. As we know, the electrical current is defined as the amount of electrical charge flowing per time. Thus we can write:


Therefore either q1 or ω can be chosen arbitrarily, and the other one has to be adjusted adequately, so that the electrical current I is correct to produce of the magnetic field which has to be emulated. We decide to chose arbitrarily q1=1Ampere, and to adjust ω. This can be done by a calculating ω from equation 10 as being shown in equation (11):

For the summation of the infinitesimal field-elements of equation (7), we could in principle solve the integral of equation (12). But the algorithm is designed for arbitrarily shaped conductor loops, and we already have N discrete finite elements, so that we can solve equation (12) by an approximation of a discrete sum as also shown in the same equation (12). But we shall keep in mind that a discrete sum always makes numerical noise (similar as shown in figure 3).

The summation has been realised within the subroutine „Magnetfeld_zuweisen_03“, with the variable of summation being I=0…N, with the aim to make the argument of the parametrization run from t=0…T=2π/ω  (for detailed understanding, please see subroutine in the appendix).

The results have been checked by the classical formula for the calculation of the field-strength (see equation 13) along the axis of the coil (which here is the y-axis), and the check confirms our results.


Now our simulation of the magnetic field of a cylindrical bar-magnet by the use of two conductor loops at each end of the cylinder is complete.

The determination of the magnetic field of a permanent magnet is not our goal, but it is one important step on our way towards the goal. Our goal finally is the determination of the torque, with which the coil accelerates or decelerates the rotating permanent magnet. Therefore we have to determine the Lorentz-force with which the currents in the coils (red colour in figure 1) does act onto the conductor loops emulating the permanent magnet. With other words: We have to calculate the Lorentz-force, which is the fundament for the calculation of the torque, which the permanent magnet experiences.

Therefore we again use Biot-Savart’s law. Now we apply it in a way that we calculate the magnetic field produced by the red coils at the position of the permanent magnet emulation conductor loops, which have been used to emulate the permanent magnet. This is necessary that we can calculate the Lorentz-force, which the permanent magnet emulation conductor loops experience within this field of the red coils

This means, that the conductor loops producing the field to be calculated now, are described by the polynomial line of the red coil, and the points of interest at which the field has to be calculated is the position of each conductor loop-element, which experiences a Lorentz-force. The situation is illustrated in figure 5.

Due to the rectangular shape of the red loops, we could in principle use a polynomial line with only for support points. But in reality this is not sensible, because we need to have several (many) finite area-elements within the red coils, so that the magnetic flux through the coils will be calculated properly (see equation 4).

Fig.5: Illustration of the vectors as being used for the calculation of the magnetic fields of the red coils at the position of the conductor loops emulating the permanent magnet as being used for the application of Biot-Savart’s law.  s=(sx,sy,sz)=conductor-loop element, r=(rx,ry,rz)=point of interest,  ds=(dsx,dsy,dsz)=motion of the electrical charges in the conductor-loop element.


With regard to the parameters according to figure 5, we can write Biot-Savart’s law according to equation (14):


As usual, we put the expression of (15) into equation (14) to calculate the finite field-elements dH, which are produced by all finite conductor elements of the red coils. This is the way how we come to the field which the red coils produce at the position of the magnet-emulation-coils. This calculation has to be done individually for each loop-elements of the magnet-emulation-coils. From there we come to Lorentz-force acting onto each magnet-emulation-coil (see equation 16).

From all finite Lorentz-force-elements dF we now calculate the finite torque-elements with which they act onto the rotation of the permanent magnet. The summation of these torque elements (as a discrete sum, see equation (17)) delivers the torque on the rotating magnet as being used in the DFEM-algorithm.

The calculations are realised in the subroutine „Drehmoment“.


Due to the spatial discretization, we also have numerical noise which has to be smoothed by a Fourier-series. Therefore we develop the torque as a function of the angle φ of the orientation of the magnet. We again restrict ourselves to fifth order or less in order to exclude high-frequency components for sure (compare (5)).

Again, the elapsed CPU-time to calculate the torque is rather large, so that the explicit torque-calculation is not recommendable within the solution of the differential-equation. Here we have again the advantage to save computer time due to the Fourier-series.

Now we reach the point, that the preparations are complete, as they are the following two steps:

-                     the calculation of the induced voltage, which the rotation of the permanent magnet induces into the coils,                                    and

-                     the computation of the torque, with which the electrical current in the coils act onto the permanent magnet.

The subroutines which do this calculations in fast manner (due to the quick Fourier-series) have the names „Schnell_Drehmoment“, „Fluss_T“ and „Fluss_I“. They can be seen in the source code in the appendix.

We summarise the results of section 1 in figure 6 and figure 7, which corresponds to the geometry of figure 5.



Fig.6:  The magnetic flux, which the permanent magnet according to figure 5 produces in the red coils depends on the orientation of the permanent magnet (angle).

Due to the symmetry of the setup, the magnetic flux has the same value for both coils.

The blue signal allows us to estimate the numerical noise, which has its reason in the spatial discretization of the geometry. The maximum of the numerical noise occurs in the moment, when the windings of the magnet emulation coils, come most close to the wires of the red coils. The purple signal is a smoothing of the blue signal as being calculated by the approximation of a 5th order Fourier series.

Fig.7:  If we later derive the magnetic flux to time, we calculate the induced voltage in the red coils (later, not printed here). Therefore we need to know the angular velocity of the rotating magnet. The curve of the torque is made to check the computation. Therefore the red coils are driven with a constant current, so that the torque depends mainly on the orientation of the magnet. The values for the parameters are:

- permanent magnet, cylindrically, 4 cm thick, 8 cm long, field strength 1 Tesla at ist end.

- „Red“ coils, rectengular, 6 cm wide, 12 cm height, located at x=–2 cm and at x=+2cm.


We explicitely emphasize, that the number of coils can be chosen arbitrarily same as the number of permanent magnets within the DFEM-algorithm presented here. This means that the algorithm can be used for the computation of every ZPE-converter working electrically or magnetically:

- If the DFEM-algorithm is used to emulate several permanent magnets, their interaction can be analysed. This is an interesting application of the algorithm for the simulation of self-running ZPE-magnetmotors, as they can for instance be found at [Hoh 11], [Jeb 06].

- If the DFEM-algorithm is used to emulate several coils, which might be connected with each other by a yoke if requested (in order to conduct the magnetic flux in an appropriate way), motionless ZPE-converters can be simulated, for instance as can be seen at [Mar 88-98], [Bea 02].

- If the DFEM-algorithm is used to emulate one coil and one permanent magnet, the “Electro-Mechanic Double Resonance” converter (EMDR-converter) as proposed by the author of this article can be simulated [Tur 11].

- If the DFEM-algorithm is used to emulate two coils and several permanent magnets fixed to each other in appropriate manner, the Keppe-motor can be simulated [Kep 10].

- If 6 cylindrical bar-magnets are mounted within 6 coils and are arranged with each other in a hexagon, the Coler-apparatus can be simulated. The behaviour of these elements will lead us to a rather complicated differential-equation, and some electrical elements forming an electric circuit are introduced as boundary conditions into the system of differential-equations. (For differential-equations, please see also section 2.) Perhaps, a simulation of the Coler-apparatus might help to decide whether this motionless-converter can work or not. [Hur 40], [Mie 84], [Nie 83]

- Also dynamic input and output of energy is no problem, because the differential-equations describing the motion can be expanded with some input- voltages, load-resistors, an so on… Such elements have to be taken into the differential-equations additionally. Also mechanical force of torque can be applied as boundary conditions in the differential-equations. In order to illustrate this, the source-code in the appendix contains two coils (see fig.5). In the source code, the left coil has of the name input-coil on the right coil has the name turbo-coil (see subroutine “U7”). Additionally there is a load-resistor („RLast“) being connected with the turbo-coil. The rotating axis of the permanent magnet is being supported initially with a given angular velocity, being applied as initial-condition for the solution of the differential-equation. This initial rotation brings energy into the system once at the very beginning of the motion and can then be disconnected. This means that the EMDR-converter is a self running motor, which needs energy support only at the very beginning of the motion to initialise the rotation. Furthermore the input-coil is not active in the source code as printed in the appendix, because of EMDR-converter does not need permanent input-energy, for it is a self running engine. Nevertheless the input-coil can be used if somebody wants to do this, as for instance for the purpose to control the “rounds per minute” of the motor.

- Furthermore the algorithm allows mechanical extraction of power. In our differential-equation this is simulated by a decelerating torque proportional to the angular velocity of the rotation (of the magnet). Such type of energy-extraction can be used to simulate friction as well as to simulate energy-extraction for technical application.

- In order to develop an exemplary ZPE-converter and bring it into a stable permanent mode of operation, it is a convenient method, to control some input-power or to control the extraction of power. With regard to a self running engine, there is no input power, so that we decided to control the output-power in the DFEM-algorithm shown in the appendix. Therefore we define a special value for the speed of revolution (rounds per minute). If the rotation is faster, the energy extraction is enhanced, and of the rotation is less fast, the energy extraction has to be reduced (respecting a given hysteresis, necessary for proper switching). Later we will find out, that our exemplary ZPE-motor can also work without a regulation, but with constand extraction of energy.

The variability of the DFEM-method is large, so that it is not restricted only to magnetic ZPE-converters, but it is also applicable to electrostatic ZPE-converters. The only necessary change is, to replace the Lorentz-force from equation (16) by the Coulomb-force as seen in equation (18).

The precision of the calculation mainly depends on the precision of the input-data, namely the data of the mechanical and electrical components as well as of the data of the interacting fields with which the components act onto each other.


2. Motion of the components of the ZPE-Converter

In section 1 we did the preparation of the necessary fundamental equations of physics. This is now done, and we can turn our attention towards the solution of the differential-equations describing the motion in the ZPE-converter. Hereby we speak about motions of the mechanical components as well as about motions of the electrical components (such as electrical charges and fields or magnetic fields).

The functioning principle of every ZPE-converter can be described by the motion of its components. The adequate means for this description are the differential-equations based on the interactions of the components with each other.

For ZPE-converters of course do not consist only of one single component, but of several components, which have to interact with each other, we always have two put up and solve coupled systems of differential-equations of higher order. If converters need input-energy (such types which do not work as self running engines, but only as over-unity systems), the energy-input has to be introduced as perturbation-function in the differential-equations. For all types of interaction with some external elements components, the appropriate means is the introduction of perturbation-functions into the differential-equations. This makes the higher order differential-equation systems inhomogeneous.

Mathematically, this has the consequence, that we cannot simply derive an analytical solution, valid for each type of differential-equation system. Consequently, the central core of computation of the DFEM-algorithm is a numerical iterative solver of the differential equations. This topic is, to what we want to focus our attention in section 2.

On this background it is clear, that the DFEM-algorithm needs some certain amount of CPU-time, if you want to have the solution with sufficient precision. You cannot expect the DFEM-algorithm to come to good convergence within few seconds.

The understanding, how the differential-equations of the motions have to be formulated, is the central point of understanding, which every user of the DFEM-algorithm has to work on. Only on the basis of this understanding, the user can apply the DFEM-algorithm onto his own system.


Let us now begin to develop the differential-equations for the example of the EMDR-converter, which define the core of computation of the algorithm. We will do this in analogy to [Tur 11].

The differential-equations by principle operate completely dynamically, so that the DFEM-algorithm has to do the same (taking the finite speed of propagation of the interacting fields into account). This has the consequence, that all physical sizes, entities and values under consideration have to be drawn back to the oscillating electrical charges or to the rotating magnet, using the fundamental entities of q, dq/dt,d2q/dt2 and o φ, dφ/dt,d2φ/dt2   , where q if the electrical charge and φ is the angle of rotation of the permanent magnet.

From [Tur 11] we learned, that it is necessary to have very fine time-steps in order to get reliable results. Consequently it will not be possible to save all data of all time-steps within special data-arrays. So the program now allows the calculation of as many time steps as required, but the data-storage will be done only for maximum 35,000 points which is a sensible upper limit for the data-export to Excel.

Values as for instance the inductivity of the cylindrical coil (see equation 19) or the momentum of inertia of the rotating magnet as a massive cylinder (see equation 20) are taken from standard textbooks of physics or engineering disciplines.


Fig.8: This is the setup of a ZPE-converter for whose simulation we want to develop the DFEM differential-equations. It consists of two coils (red colour) and one rotating permanent magnet (black colour) and additionally an input-voltage (blue colour) and one capacitor (green colour).

For our EMDR-converter as our very example of a ZPE-converter can be operated as a self running motor, the input-voltage is not applied at all, so that the left coil together with voltage-supply is taken away completely. Within the source-code, those lines for the simulation of the input-voltage-components are left away; they are marked as comment so that they do not take part at the calculation at all. In order to make the setup symmetrically, the right coil was moved symmetrically around the axis of rotation of the permanent magnet. In order to make clear to everybody how the setup is now looking, figure 9 was drawn. Those elements which are not necessary for the simulation of our exemplary EMDR-converter are only printed as comments in the source-code of the algorithm for those who want to use the algorithm for other ZPE-motors. These comments have only the purpose to help colleagues to simulate their ZPE-machines.

Fig.9: This very simple setup of an EMDR-converter already allows the powerful conversion of ZPE-energy into classical electrical and classical mechanical energy. Its differential equation system was formulated and solved as explained on the following pages, together with the results describing the behaviour of this ZPE-motor.


The development and formulation of the differential equation system is done rather similar as in [Tur 11].


The differential equation of a harmonic oscillation of an electric LC-oscillation circuit is described by equation (21).



The differential equation of an attenuated oscillation of an electric LCR-circuit is described by equation (22).


Its numerical iterative solution, as being achieved with the solver of the differential equations in our DFEM-algorithm, is described in equation (23) in analogy with [Tur 11], which is in good agreement with the classical solution (see figure 10).


Fig.10: Verification of an attenuated oscillation in the electric LCR-oscillation circuit for the check of the differential equation (22) and the solution (23).


Different from the electrical oscillation, the mechanical motion (of the permanent magnet) is not an oscillation but a rotation and thus it has no restituive force. So we have to take two contributions into account to describe the torque M. The first contribution goes back to the magnetic field from the coil, which acts onto the permanent magnet. The second contribution goes back to the mechanical extraction of power, which will be the dominant part of the output-power of our ZPE-motor. This second contribution is made by the torque proportional to the angular velocity of the rotation. (Its mechanism could be friction or some other mechanism as well.) The first mentioned contribution is well-known from section 1. The last mentioned contribution will be discussed in detail in section 6.

Consequently the differential equation of the mechanical part of the system can be integrated rather simply as shown in equation (24). This also indicates that equation (24a) contains the coupling of the electrical part of the system into the mechanical part of the system.


Still missing in our system of differential equations, and thus to be introduced now, is the coupling of the mechanical part of the system into the electrical part of the system. The appropriate means therefore is the induced voltage (and this is the reason, why we did its calculation), which the rotation of the permanent magnet brings into the coil. Therefore the differential equation of the electrical system is expanded, from equation (22), so that we come to equation (25).



On this basis, the numerical-iterative solution of the electrical part of the differential equation system can be integrated according to equation (26) - and by the way in analogy to [Tur 11].

This was the explanation, how the main program of the DFEM-algorithm performs the calculation of the motions of the components of the ZPE-converter, analysing on the one hand the motion of the electrical charges and on the other hand the motion of the rotating magnet. This is a dynamical computation. Therefore the algorithm has the dynamical FEM (=DFEM).

General remark:

Within our calculations it is only allowed to use formulas which follow the dynamics of the system. Also the physical sizes used in these formulas have to respect that this criterion. This means that physical sizes as for instance the mean value of an electrical current, some effective values, and so on… are strictly forbidden. Even physical sizes, entities and values, which somehow refer to a special signal-shape are not allowed, because they do not follow the full dynamics of the differential equations. Those who want to adopt the DFEM-algorithm to their own machines and/or experiments have to be very careful, to respect this criterion without any exception. This is important.

Fundamental philosophical remark regarding the propagation of the fields (and first of all the speed of propagation of these fields) within the electrical circuit:

An electrical LC-oscillation-circuit (which for instance can be seen in figure 9) has an electrical oscillation within the coil and the capacitor.

The question is now: Which type of physical entities do oscillate in this circuit ?

Is it electrical charges which oscillate back and forth ?

No - this is for sure not the case !  This can be understood rather easy, when we follow the electrical charges within the LC-circuit, beginning at the moment of time, at which the coil is free from any electrical current. This is the moment, in which the capacitor is charged up to its maximum voltage and charge, so that the total energy of the oscillation circuit is now stored as electrostatic field energy of the electrostatic field between the capacitor plates. From now on, the capacitor begins to discharge, so that electrical fields (and voltages) propagate along the wire of the coil. From the positive capacitor plate one field propagates towards the negative capacitor plate, and in the opposite way another field propagates from the negative capacitor plate into the direction towards the positive capacitor plate. But the propagating entities are only fields, not charge-carriers (such as for instance electrons or electron holes). There are two reasons explaining this argument: The first reason is, that the charge-carriers can not propagate with the speed of v=1/√(LC), which we know to be the speed of the propagating electrical signal. The second reason is, that positive and negative charge-carrier would compensate each other as soon as they meet each other in the middle of the coil. If this would happen (as for instance electrons and electron holes would compensate each other), the complete oscillation would stop as soon as the charge carriers from the different capacitor plates meet each other. This would be the case after one quarter of a period of oscillation. Everybody knows that this is in contradiction with the real observation (at LC-oscillation circuits), so obviously the oscillating physical entities are not charge carriers.


But which are the oscillating physical entities - is it electrical fields ?

Yes - this is indeed the case !  Among the typical properties of electrical fields (as well as electrical waves) is the ability to superpose without disturbing each other, and even to cross the path of each other without taking notice of each other. So we can imagine one field as a wave crest and the other one as a way through, both of them moving from one capacitor plate to the opposite plate, not interfering with each other when they pass the coil at the same moment. Both of them follow their speed of propagation v=1/(LC), and thus the oscillation circuit behaves as we know it from wave theory.

Because of their ability to pass each other without disturbing each other, the positive wave crest and the negative way through reaches the opposite capacitor-plate after half an oscillation, and they take the other half of the oscillation to find their way back “home”, and so on… (cyclically…).

Of course the propagating fields cause small displacements of the charge-carriers within the wire of the coil, but these displacements are rather small - by several orders of magnitude too small to transport charge-from one capacitor plate the opposite one.

We can understand this rather easy, when we have a look to the acoustic analogon. If an acoustic signal propagates within a tube (the one-dimensional consideration makes it easier, and it is in good agreement with the one-dimensional propagation of the electrical signals in the wire), we can also send the wave crest and the wave trough from the opposite ends of the tube, and we will also see, that they pass each other without disturbing each other. What we regard here is the field of air-pressure, in which waves always superpose without disturbing each other. The conception is the following: The tube can be subdivided in small finite volumes, containing gas atoms of the air. And these volumes change their positions by a very small amount, as soon as the pressure field is coming. We face the same situation of the electrons in the wire of the coil. We can imagine small finite volumes, containing electrons. And these volumes alter their positions by a small amount as soon as they are exposed to an electrical field.

The gas atoms of the air are moving very fast, same as the electrons in the wire. But the small finite elements of volume filled by electrons resp. gas molecules only move with a very moderate speed, which we know as the “drift-speed” in the case of the electrons and as the “acoustical velocity” in the case of the gas molecules. Nevertheless signals and waves are propagating with typical signals speed, which we know to be the “speed of sound” in acoustics and v=1/(LC) in electrodynamics.

The signal speed defines the speed of propagation of the field and waves in the wire of the LC-oscillation-circuit. This means that we really control the speed of propagation of the interacting fields (of the electromagnetic interaction), when we adjust the inductivity “L” and the capacity “C” of the oscillation circuit.

This explanation demonstrates, how the DFEM-computation is traced back to the conversion of ZPE-energy, as explained in detail in [Tur 10a] and [Tur 10b]. It shall help the readers of this publication to understand, how the LC-oscillation-circuit is indeed used, to control the speed of propagation of the interacting fields, as being necessary for ZPE-conversion.

Now the principle our calculation is explained. The execution of the calculation can be seen in details in the algorithm in the source-code. And it is clear, that the calculation will give results. So, this is the time to begin to analyse the results. This is, to what we will turn our attention beginning with the next section.


3. Evaluation of the results of a converter example

A very first evaluation of the behaviour of the converter-system can be done with some first results, which the program displays on the screen, directly when the program is running. Table 1 is a list of the very first and important physical values for the evaluation of a ZPE-converter. The table is of course made for the example of our EMDR-converter, for which we later will have concrete hints how to build it up.


Several further physical values should be analysed due to their dynamic behaviour as a function of time. They are listed in table 2. The program exports these data into a file which can be read by Excel. There they are available for being displayed graphically. Thus table 2 also contains information about the column in which the data are to be found in Excel.

Now our DFEM-algorithm is developed so far, that we can perform realistic computations of arbitrary electric and/or magnetic ZPE-converters. An example therefore shall be calculated in the following sections - namely with the EMDR-converter suggested by the author of this publication.


4. Computation example for a concrete ZPE-motor

The namely reason for the development of the DFEM-algorithm presented here is the fact, that the author wants to show realistic computations on his EMDR-Converter („Electro-Mechanical Double-Resonance“ Converter), and to get reliable results with sensible precision on this machine. This has the purpose to prepare building up an experimental prototype. Even if the author of this article can not built up such a prototype by himself due to his very restricted possibilities, he hopes, that many colleagues will read this article and try to build up prototypes by themselves.

In order to make it most efficient for those colleagues, who already began to think about building up such prototypes, the author decided to develop the new calculations for a setup which is a rather similar to this one in his former calculation performed in [Tur 11]. The setup is, what we saw in figure 9.

The definition of the geometry in the computer program needs a set of the 32 input parameters. Additionally the program works with some constants of nature and some other parameters derived from the input parameters (additionally 18 such parameters), which are displayed on the screen during the runtime of the program. Sensible definition of the input parameters needs exact fine tuning of their values, and thus take several hours/days even when the author does it by himself.

Caution: The parameter-set displayed on the following pages corresponds directly to the geometry and setup of figure 9. Different design needs different parameters. Furthermore the program contains several subroutines, which are developed to perform automatic meshing for the finite element method. But this automatic meshing only works, if the setup is not altered by principle. If somebody wants to change the design of the converter, it will be necessary to change the subroutines as well.


(a.) Definition of the geometry of the setup

We start now with the explanation of the input parameters, which the realistic DFEM-algorithm requires. For didactical reasons we now do not want to take mechanical power extraction into consideration, because if we leave this aspect for later, it will be easier to understand the converter. Mechanical power extraction will be introduced later within this publication.

(Constants of nature, not Input-parameters:)

- epo:=8.854187817E-12(As/Vm);  (electric field constant)

- muo:=4*pi*1E-7(Vs/Am);       (magnetic field constant)

- LiGe:=Sqrt(1/muo/epo)(m/s);   (speed of light)

(For the solution of the differential-equations and for the display of the results:)

- AnzP:=5000000;   (number of time steps of the numerical iteration)

- dt:=1E-6; (Sec. )   (Duration of each single time step)

- Abstd:=1;    (only for preparation, do not alter the value)

- PlotAnfang:=0000;              (For Data-export to Excel: First Plot-Punkt)

- PlotEnde:=5000000;    (For Data-export to Excel: Last Plot-Punkt)

- PlotStep:=200;      (For Data-export to Excel: step width of the data being exported)

(Remark: Excel is restricted to maximal 32.767 Data-groups. If the number of time steps is larger than this value, not all computed data can be plotted graphically by Excel. In our example, only every 200th point is being  exported to Excel.)

(For the definition of the geometry of the coils (DFEM-meshing is done automatically):

- Spsw:=0.01;                (Meters of step width of the meshing)

- xo:=0; yo:=6; zo:=5;   (Geometrical parameters according Fig.1, steps of Spsw)

- Ninput:=80;        (number of windings of the Input-coil, left coil in figure .1)

- Nturbo:=12;           (number of windings of the Turbo-coil, right coil in figure .1)

- nebeninput:=8;      (windings side-by-side in Input-coil)

- ueberinput:=10;     (windings on top of each other Input-coil)

- nebenturbo:=3;       (windings side-by-side in der Turbo-coil)

- ueberturbo:=3;       (windings on top of each other Turbo-coil)

(Remark: Here the parameters are used to define rectangular coils according to Fig. 1. The cross-section of the Input-coil consists of 8 windings side-by-side and 10 such layers on top of each other. The cross-section of the Turbo-coil consists of 3 windings side-by-side and 3 such layers on top of each other. “On top of each other”means, that the layers are built up radially.)

(For the emulation of the permanent magnet:)

- Bsw:=1E-2;      (Meters) (The magnetic field shall be stored in steps of centimetres.)

- MEyo:=0.05;     (Half length of the cylindrical bar magnet)

- MEro:=0.01;     (Radius of the cylindrical bar magnet)

- MEI:=15899.87553475;      (Amperes, current in the coil is to emulate the permanent magnet)

(Remark: We here use a cylindrical bar-magnet according to Fig. 8 and Fig.9. The shape can be altered if required, but the meshing-subroutines have to be altered also.)

(Remark regarding the data-storage of the magnetic field: The magnetic field is fixed rigidly to the magnet within a sufficiently extended volume of space. The values of the field strength are stored in a data-array at finite geometrical steps. When the magnet is moving, the field is moving together with the magnet. The step length for the data-storage of the field is „Bsw“.)

(Remark: The weird value for the electrical current in the magnet-emulation-coils has its reason in the fact, that a given value of the magnetic field has to be emulated. The field strength of this field to be emulated is displayed on the screen during the runtime of the program. It can be read in order to adjust the current MEI in such way that the required field strength is achieved. In our example we have a magnet with 1 Tesla at its ends.)

(Further technical dimensions:)

- DD:=0.10; (Meter)    (Thickness of the wire from which the coil is made)

- rho:=1.35E-8; (Ohm*m)    (Specific electrical resistance of copper, [Koh 96])

- rhoMag:=7.8E3; (kg/m^3)    (Density of the magnet-material, Iron, [Koh 96])

- CT:=36.61E-6; (Farad)    (Capacity within the oscillation circuit of the Turbo-coil)

- CI:=100E-6; (Farad)                (Capacity within the oscillation circuit of the Input-coil)

(Remark: In our example, the input-coil has been modelled in order to prepare it for everybody who wants to have an additional coil in the DFEM-algorithm. But the input-coil is not used for the computation of the converter, and the complete oscillation circuit containing the input-coil is left away for the computation of the self running EMDR-converter.)

- Rlast:=0.0111; (Ohm)    (Ohm'ian load resistor in the Turbo-circuit for the extraction of energy)

- UmAn:=50000; (U/min) (Initial angular velocity: mechanical boundary condition - rotating magnet)

- Uc:=0;(Volt)  Il:=0; (Ampere) (electrical boundary conditions – capacitor-voltage, current in the coil)

(Remark: Even a self running ZPE-motor needs an initial energy to be started. It will not start just by alone. The initial energy can be supplied mechanically (as it is done in the example here), but it can be supplied electrically as well, for instance by charging the capacitor the coil in order to initialize the operational the motor.)

- U7(t)=0;

(Remark: If the ZPE-converter is not a self running engine, but only an over-unity engine, it permanently needs some classical input-energy for operation. This can be supplied mechanically (at the axis of rotation) or electrically with some input-voltage. The last version is displayed in the source-code printed in the appendix, by the use of a subroutine with the name “U7”. Nevertheless this input-voltage is not used for the solution of the differential-equations, because the EMDR-converter is a self running engine and does not need any classical input power.)

(Composed Parameters, for the purpose to control the input data. Do not use them for input.)

- DLI:=4*(yo+zo)*Spsw*Ninput; (Meter)   (Length of the wire of the Input-coil)

- DLT:=4*(yo+zo)*Spsw*Nturbo; (Meter)    (Length of the wire of the Turbo-coil )

- RI:=rho*(DLI)/(pi/4*DD*DD); (Ohm)    (Ohm`ian resistance of the wire of the Input-coil)

- RT:=rho*(DLT)/(pi/4*DD*DD); (Ohm)    (Ohm`ian resistance of the wire of the Turbo-coil)

- BreiteI:=nebeninput*DD; HoeheI:=ueberinput*DD;    (Width and height of the Input-coil)

- BreiteT:=nebenturbo*DD; HoeheT:=ueberturbo*DD;    (Width and height of the Turbo-coil)

- fkI:=Sqrt(HoeheI*HoeheI+4/pi*2*yo*2*zo)/HoeheI;    (Induktivity-correction for short coil)

- fkT:=Sqrt(HoeheT*HoeheT+4/pi*2*yo*2*zo)/HoeheT;   (Induktivity-correction for short coil)

- LI:=muo*(2*yo+BreiteI)*(2*zo+BreiteI)*Ninput*Ninput/(HoeheI*fkI);

(Geometrical average => Induktivity of the Input-coil)

- LT:=muo*(2*yo+BreiteT)*(2*zo+BreiteT)*Nturbo*Nturbo/(HoeheT*fkT);

(Geometrical average => Induktivity of the Turbo-coil)

- omT:=1/Sqrt(LT*CT);    (Resonance-angular-frequency of Turbo-cicuit of LT & CT)

- TT:=2*pi/omT;         (classical duration of oscillation of Turbo-circuit of LT & CT)

- Mmag:=rhoMag*(pi*MEro*MEro)*(2*MEyo);    (Mass of the Magnet)

- J:=Mmag/4*(MEro*MEro+4*MEyo*MEyo/3);         (moment of inertia of the magnet of rotation)

- omAn:=UmAn/60*2*pi;    (Start angular velocity (rad/sec.) of the rotating magnet)

- UmSec:=UmAn/60;           (Start angular velocity, rotating Magnet (rounds per second))

(Remark: Some of these values are not only necessary to control the input data but they are also necessary for the further computation in the DFEM-Algorithm.)


With these parameters, the ZPE-motor is modelled as being displayed in figure 11. As explained above, there is only the Turbo-coil and no Input-coil.

Fig.11: EMDR-converter with a rotating magnet (red colour: length 10 cm, thickness 2 cm). (Drawn by the use of ANSYS [Ans 08].)

Comments to Fig.11: Also in red colour we see a thin axis of rotation, around which the magnet can rotate. At each end of the axis, there is a bearing (also red) to keep the axis of the rotating magnet at its position.

The windings of the Turbo-coil are drawn in light blue colour. As we see, the field flux lines of the rotating magnet cut the windings of the coil exactly perpendicular.

Because of practical aspects (as can be seen later), the windings of the coil have to be made from rather thick material, and the coil should be made of not very many windings. In our example, the cross-section of the coil has an inside area of 5 × 6 cm2, and the thickness of the wire 10 mm.


The time steps for the numerical iterative the solution of the differential-equations have to be choosen fine enough, that every period of motion consists of sufficiently many steps of calculation. Thefore the parameters as displayed above have time steps of dT=1µsec. at an angular velocity of 50,000 rpm. This is for sure not sufficient. The time steps must be much smaller.

But during the practical application of the DFEM-algorithm, we start with rather rough time steps in order to get a first rough feeling for the behaviour of the machine. Rough time steps save computer (CPU)- time. We can use such rough times steps, to see what is happening, during the phase of definition of the geometry of the machine. We can also use such rough time steps for the maximisation of the converted power, when we want to alter the system-parameters by hand. As soon as this task is done, we have to reduce the time steps remarkably in order to check the convergence of the algorithm. Good convergence is achieved, when we see that any further reduction of the length of the time steps does not have a remarkable influence on the results. As soon as we reach the condition of good convergence, we have to begin the fine tuning of the system parameters, together with a renewal of the maximisation of the converted power.

In order to make the optimisation of the system parameters as efficient as possible, there is a data-output routine included into the program. This allows a fast evaluation of the behaviour of the system, when we alter the system parameters by hand (by trial and error). During this phase, we see the following data on the screen:

- The start energy within the system, which is brought into the differential-equations by initial conditions. This energy is computed among others from the parameters UmAn, Uc, Il.

- The energy within the system at the end of the time of analysis. Therefore all amount of energy within the system summed up, this is the energy of motion of the magnets, the energy within the coils and the energy within the capacitors.

- The amount of energy being gained within the system during the analyzed time of operation. This is the difference between the start energy and the energy at the end of the analysis-time. If the energy-gain plus the energy being extracted from the system is positive, the system converts ZPE-energy into classical energy. If this energy-gain plus the energy being extracted is negative, the system converts in the opposite manner classical energy into ZPE-energy. An ideal classical electromotor will contain the same start energy plus input energy, as it contains energy at the end of operation plus energy being extracted during operation. (Friction to be taken into account as energy-extraction.)

- The power corresponding to the energy gain, which is calculated by dividing the energy gain trough the time of observation.

- The energy being extracted from the Ohm’ian resistor during the time of operation, and the power corresponding with this energy gain. Not very interesting is the knowledge about the Ohm’ian losses in the wires of the coil. These losses are included into the calculation, but they are not printed on the computer screen, because we do not have any influence on them.

- The average of the mechanically extracted power and the sum of the mechanically extracted energy during the time of observation. Mechanical energy and power can be extracted from the rotating shaft (which is located in the middle of the permanent magnet). Mechanical energy is not yet regarded now in the sections 4 and 5, but soon in section 6 it will be regarded detailed.

- Additional to the mechanically extracted power and energy, there is some mechanical gain of energy, which remains inside the system.

- The input-voltage and the input-power in connection with the input-supply is printed on the computer screen, but it is ZERO, because there is no classical energy brought into the system during operation. The EMDR-converter is a self running system. The values are only displayed for those, who want to modify the DFEM-algorithm to simulate an over-unity machine.

- The total duration of the observation, which is the sum of all time-steps “dt”. We do not speak about the elapsed CPU-time, but we speak about the simulated time of operation during which the converter is running.


The simple online data-evaluation as described above helps the user to get a quick impression about the mode of operation of the converter, which allows to perform a variation and optimization of the system-parameters by hand. This means, that we can alter the values of the system parameters, check the behaviour of the system, alter the values again, check again, and so on… By this means we should be capable to develop a design which should work properly (at least theoretically).

As soon as this design is found, it is recommended to evaluate the system more precise. This can be done by observing those variables as a function of time, which follow the dynamics of the motion. Therefore the time dependent variables of the system are exported into a data-file for Excel, where they can be displayed graphically. The listing of the contents of the different Excel-columns have been printed above in Table 2. From the evaluation of these data, we have to answer questions such as the following:

- Does the system run into a stable mode of operation ?

This can be seen for instance, when we regard the angular velocity ω as a function of time, because a stable mode of operation can only be achieved if the angular velocity is constant, or at least if it oscillates around a constant value (within a well controlled hysteresis). If this is not the case, the converter is still in the phase of initialisation, or it does not run stable at all. In order to decide between these both possibilities, the length of the time steps "dt" should be reduced, and the total oberservation time has to be enhanced. Then the calculation shall be repeated with longer observation time.

The same observation can be done with regard to the electrical current in the coil or with the voltage in the capacitor.

- Does the machine convert enough ZPE-energy, so that it will not be brought to standstill by friction ?

Therefore we check the energy of the mechanical rotation remaining in the rotation of the permanent magnet, as a function of time. The gain should be large enough, that we can expect, that it is sufficient to surmount the energy loss due to friction. If the calculation is done taking mechanical power extraction into account (see section 6), we can define the coefficient of friction and check the extracted power (as a function of time).

- Also the time dependent behaviour of the electrical currents and the voltages in the LC-oscillation circuit should be observed graphic only, because this contains more reliable and detailed information than the simple estimation of the maximum, which is printed during the runtime on the screen. If the machine is running properly as a ZPE-converter, the values of the voltage and the current are nomaly rather large, so that we have to be careful not to overload the wire of the coil or the capacitor.

By the way it should be mentioned, that the elapsed CPU-time can take several minutes or several hours, especially when the time steps are very short (for instance in the range of few nanoseconds).


5. A concret EMDR vacuum energy converter

The DFEM-program in the appendix can only run, if there is a data-file in the same directory with the name 'schonda', which can be downloaded together with source-code of the DFEM-program for free from the Internet-page of the author of the publication presented here. The data-file 'schonda' only has the purpose to save CPU-time, as following: During the phase of initialisation of the main program, all preliminary work takes some CPU-time, which is only necessary, if the parameters describing the geometry of the setup have been altered since the last run of the program. (In this case, the automatic meshing has to be renewed.) If the geometry-parameters remain unchanged since the last run of the program, the results of the initialisation can be taken directly from the last run of the program (together with the existing mesh, stored in 'schonda'). This is exactly what the program does, when it reads the data-file 'schonda', in order to save the time for the same initialisation which already has been calculated before. This makes it more efficient to repeat the program several times during the phase of the optimisation of the system parameters.

We now (in section 5) want to discuss the system parameters, with the values as printed in the source-code in the appendix.


The typical construction of a EMDR-converter begins with a search for an appropriate permanent magnet. As soon as the magnet is found and its field strength is measured (for instance with a Hall-probe), its dimensions are put into the input-pata lines of the source-code. Then we have to construct the conductor loops for the emulation of the permanent magnet and to apply an appropriate electrical current within the conductor loops, in order to reproduce the measured values of the magnetic field.

Then we have to modulate the Turbo-coil and to put all the other requested data into the DFEM-algorithm. Finally the initial angular velocity of the rotating magnet has to be adopted, and the very last step is the adjustment of the capacitor in the Turbo-oscillation-circuit. The criterion of the initial angular velocity of the rotating magnet shall be decided from the stability of the bearing keeping the rotating axis. The bearing must withstand the angular velocity of the stable mode of operation, which can be remarkable larger than the initial angular velocity. Thus the value of the maximum possible angular velocity of the magnet plays an important role for the adjustment of the capacitor. High-speed rotation has the consequence to enhance the amount of power being converted from the ZPE of the quantum vacuum. If the system operators as a ZPE-converter, the angular velocity is increasing during the initial phase of the operation. And we can use the amount of the increase of the angular velocity as an indicator for the quality of adjustment of the system parameters. The more ZPE-power we convert, the more increase of the angular velocity we observe. At the very beginning of the adjustment procedure, we apply a very small load resistor and adjust the capacity to a maximum of increase of the angular velocity of the magnet.

Rather often (depending on the geometry of the system) we observe, that the mechanical power-gain is much larger than the electrical power-gain. This is also the reason, that we will soon (in section 6) have to take the real benefit of the converter mechanically from the rotating shaft.

The next step of the optimisation now consists in enhancing the load resistor in many small steps, and always readjusting the capacitor, while checking the mechanical and electrical power gain. Depending on the configuration of the system parameters, the mechanical and the electrical power can increase both or decrease both at the same time, but it is also possible that one is increasing and the other one is decreasing. The load resistor should not be enhanced to much, otherwise it will attenuate the LC-circuit rather strong, and thus prevent the engine to start properly.

If the mechanical power gain of the machine is rather large, this is absolutely no problem, because it simply indicates that we have enough mechanical power, to overcome friction without any problems. The way how to extract mechanical power will be the topic of section 6.


Even the theoretical adjustment procedure makes clear, that the load resistor as well as the capacitor have to be adjusted with very high precision. The consequence is, that we need a capacitor and a load resistor for the practical setup, which can be adjusted very precisely. The precision of the adjustment-quality should be at least somewhere between 1% and 0.1%. When we will soon see, how large the voltage and the electrical current in the LC-oscillation-circuit really are, we will understand that this defines a really serious requirement to the capacitor and to the resistor.

The example under discussion here is not optimized with regard to the converted power, because of this optimisation will be much better soon in section 6. Nevertheless it should be mentioned, that the converted power can be enhanced remarkably, when we enhance the number of windings in the Turbo-coil a little bit (accepting the disadvantage, that the voltage in the capacitor will be enhanced strongly if we do so).

The most effective way to enhance the converted power is the angular velocity in the stable operation. Enhancing the angular velocity has the consequence to enhance the converted power tremendously. 30,000 rounds per minute as used in the example here is not a really high angular velocity. Such a rather moderate value has the purpose to make it easier to find an adequate bearing for the rotating permanent magnet. But if we have the typically angular velocity of Turbo-rotor in the automotive industry in mind (which spins with about 100,000 rounds per minute or even more), we see that an enhancement of the angular velocity should not be very difficult - and with it a remarkable enhancement of the converted power. Even much higher angular velocity is known from dentistry and from turbomolecular vacuumpumps in ultra high vacuum technology. If we could use such speedy rotation, the restriction to the power density of the ZPE-converter should probably be the electrical current in the copper wire, which should not be too strong, so that the copper wire will not get too hot.

If the results of the algorithm shall be documented, there is a simple option to press the “D”-button before you leave the program with the last , and the program will write the input-data as well as the most important results into a file named „Auswertung“. This file can be read with a text program. An example for such a file is printed here.

DFEM-Simulation of an EMDR-motor (here still without mechanical power extraction)

Parameters for the solution of the differential equation and the output of the results:
AnzP =     10000000    (number of time steps for the observed operation)
dt   =  2.000E-0007 (seconds,   duration of each time steps for the solution of the differential equations)
Abstd= 1   (only for preparation, do not alter the value)
PlotAnfang =  0       (first point for the Data-export to Excel )
PlotEnde   =  10000000      (first point for the Data-export to Excel )
PlotStep   =   400       (step width for the Data-export to Excel )
(Definition of the both coils (turbo and input):)
Spsw   =     0.010000     (Meters: step width for the automatic mesh-generation of the coils)
xo = 0,     (number of steps of Spsw)
yo = 6,     (number of steps of Spsw)
zo = 5,     (number of steps of Spsw)
Ninput =   100       (number of windings of the input-coil)
Nturbo =         9      (number of windings of the turbo-coil)
nebeninput =        10   (windings side-by-side of the input-coil)
ueberinput =        10    (layers of windings of input-coil)
nebenturbo =         3    (windings side-by-side of the turbo-coil)
ueberturbo =         3    (layers of windings of the turbo-coil)
Bsw =  1.0E-0002      (spatial resolution of the computation of the magnetic field)
MEyo =  5.00000E-0002      (y-coordinates of the conductor-loops for the emulation of the permanent magnet)
MEro =  1.00000E-0002      (Radius of the conductor-loops for the emulation of the permanent magnet)
MEI  =  1.58998E+0004      (Electrical current in the conductor-loops for the emulation of the permanent magnet)
further technical dimensions:
DD =    0.0100000  (Meters)    (diameter of the wire for the coils)
rho =  1.35000000000000E-0008  (Ohm*m)     (Specific electr. resistance of copper, depending on temperature)
rhoMag =  7.80000000000000E+0003  (kg/m^3)     (density of the permanent magnet, Iron)
CT =  9.83000E-0005  (Farad)    (capacitor in the turbo-circuit)
CI =  1.00000E-0004  (Farad)     (capacitor in the input-circuit)
additional necessary values:
Rlast =  6.400000E-0002  (Ohm)     (Ohm´s load resistor in the LC-Turbo-circuit
UmAn =   30000.00  (U/min)    (mechanical initial conditions - rpm of the rotating magnet)
Uc =       0.00  (Volt)        (electrical initial conditions - voltage at the turbo-capacitor)
Il =       0.00  (Ampere)    (electrical initial conditions - electrical current in the turbo coil)
composed parameters. These values are calculated, no input possible:
DLI:=4*(yo+zo)*Spsw*Ninput =   44.00000 (Meter,   length of the wire oft he Input-coil)
DLT:=4*(yo+zo)*Spsw*Nturbo =    3.96000 (Meter,  length of the wire oft he turbo-coil)
RI:=rho*(DLI)/(pi/4*DD*DD) =    0.00756 (Ohm)   (Ohm`resistance of the wire of the Input-coil)
RT:=rho*(DLT)/(pi/4*DD*DD) =    0.00068 (Ohm) (Ohm`resistance of the wire of the turbo-coil)
BreiteI:=nebeninput*DD =    0.10000      (Width of the input-coil)
HoeheI:=ueberinput*DD  =    0.10000      (height of the input-coil)
BreiteT:=nebenturbo*DD =    0.03000      (Width of the Turbo-coil)
HoeheT:=ueberturbo*DD  =    0.03000      (height of the Turbo-coil)
fkI:=Sqrt(HoeheI*HoeheI+4/pi*2*yo*2*zo)/HoeheI =  123.61179      (factor of correction for inductivity)
fkT:=Sqrt(HoeheT*HoeheT+4/pi*2*yo*2*zo)/HoeheT =  412.02703    (factor of correction for inductivity)
LI:=muo*(2*yo+BreiteI)*(2*zo+BreiteI)*Ninput*Ninput/(HoeheI*fkI) =  1.24238647244960E-0001 (Induktivity, Input-coil)
LT:=muo*(2*yo+BreiteT)*(2*zo+Breitet)*Nturbo*Nturbo/(HoeheT*fkT) =  9.93606632469255E-0004 (Induktivity, Turbo-coil)
omT:=1/Sqrt(LT*CT) =  3.19974964955735E+0003    (classical angular frequency of the Turbo-circuit)
TT:=2*pi/omT =  1.96364903362012E-0003    (classical period of the Turbo-circuit )
Mmag:=rhoMag*(pi*MEro*MEro)*(2*MEyo) =    0.245 kg    (Mass of the Magnet)
J:=Mmag/4*(MEro*MEro+4*MEyo*MEyo/3) =  2.10329628157837E-0004 (moment of inertia of the rotating magnet)
Several parameters, to be calculated from the above values:
Magnet: start angular velocity.: omAn =     3141.592654 rad/sec
Magnet: start angular velocity, Umdr./sec.: UmSec =  500.0000000000 Hz
Mass of the Magnet =   0.245044 kg
moment of inertia of the rotating magnet:  2.10329628157837E-0004 kg*m^2
total duration of observation:  2.00000000000000E+0000 sec.
Excel-Export:  0.00000E+0000… 2.00000E+0000 sec., Step  8.00000E-0005 sec.
These are    25000 data-point tob e exported to Excel.

Some of the results of the computation:
initial energy within the system:     1037.93511187 Joule
final energy within the system:       1171.18167853 Joule
power gain in the system:   66.62328333 Watt
energy extracted at the load resistor =    13.29139772 Joule
power extracted at the load resistor=   6.64569885828765E+0000 Watt
inserted energy by voltage supply:  0.00000000000000E+0000 Joule
inserted power by voltage supply:   0.00000000000000E+0000 Watt
total duration of the observation      2.00000000000000E+0000 sec.


This is the documentation of a set of important data written automatically by the program. We now want to discuss these data and further more other data:

The duration of observation is 2.0 seconds (the elapsed CPU-time is much larger on a normal computer). During this time-interval, the system gains a mechanical power of 66.62 Watts, and additionally an electrical power of 6.645 Watts. For both of them, there is no classical source of energy, so the energy can only come from some non-classical source which we regard as invisible for classical eyes. According to former explanations, it should be the energy of the electromagnetic zero-point waves of the quantum vacuum.


If we enhance the duration of the observation, we see that the converted power is reduced. The reason is not an incomplete convergence of the algorithm (for instance due to the duration of the finite time steps “dt”), but the reason is the fact, that the engine has totally different behaviour within the initial phase then after reaching a stable equilibrium mode of operation. We see this when we regard figure 12, which displays the angular velocity of the rotating magnet as a function of time (this is column ”I” in the Excel-data, of the automatically written file with the name „test“.) Further explanation will follow later.

 Fig.12: Angular velocity of the rotating magnet, which increases as function of time.

At the begin of the initial phase, the energy gain is rather large. After some time the operation converges to saturation, but within the two seconds under analysis, it does not yet reach the saturated stable condition.

The voltage of the capacitor can be calculated from the electrical charge inside the capacitor as being plotted in figure 13. Please notice, that the calculation was done with 10.000.000 time steps, which is far too much to be resolved by any computer graphics. Consequently we only can see the envelope of the oscillation.

Fig.13: In green colour we see the envelope of the oscillation of the electrical charge, which is flowing into and out of the capacitor.

The maximum of the voltage-amplitude is to be found at Uc,max=Qmax/Cturbo=0.01C/98.2µF=102Volts.

It should not be a problem to get such capacitors. If the number of windings is enhanced, it is possible to work with higher voltage, which allows to enhance the electrical power in the circuit as well as the power ready for extraction (at load-resistor as well as mechanically).

The current in the coil, shown in figure 14, has a maximum of a bit more than 30 Amperes. If the number of windings (of the Turbo-coil) is enhanced and the load resistor is adjusted together with other systems parameters (such as the capacitor) at the same time, a strong enhancement of the current is possible. With our wire according to figure 11, with a diameter of 10 mm (cross-section of 78.5 mm²), a current of 30 Amperes should not be a problem at all.



Fig.14: Envelope of the oscillating current in the coil as a function of time.


The voltage over the coil, which graphically has a shape very similar to figure 13 and 14, should have an amplitude of the same order of magnitude as the voltage of the capacitor, as soon as the system is adjusted sensible. In our example we find a value of 99 Volts, which fulfils this criterion surprisingly well.

The angular acceleration of the rotating magnet, of which we see the envelope in figure of 15, displays a remarkable oscillation within each period of rotation. This makes it clear, that there are remarkable Lorentz-forces between the magnet and the coil of the LC-circuit. This indicates, that those colleagues who want to build up the design in a practical experiment, should take care of mechanical forces, which should be done by a stable fixation of the coil and the axis of the magnet.

A stable mode of operation is achieved, as soon as the oscillation of the torque (and the angular acceleration) is symmetrically around the abscissa - which is not the case here in our example, as can be expected from the discussion of figure 12. As we see, even at the end of the observation-time, the average of the angular acceleration is clearly positive, so that the rotor did not yet come to its maximum (final) angular velocity.

Fig.15: Envelope of the angular velocity of the rotating magnet as a function of time.


The question about the quality of the numerical iteration can be answered, when we focus the plot to a high-resolution interval of time, so that we do not only see the envelopes of the curves, but the curves themselves. This is shown in figure 16, showing that time-zoom of the total duration of about 5 milliseconds, within which we can resolve about five periods of oscillation. The plot displays the angular acceleration. If the curves are irregular or not smooth, we have the typical case of too low scanning frequency (too long time steps “dt”), so that we face the necessity to reduce the length of time steps. Of course such an enhancement of the precision of the calculation does enhance the elapsed computer-time for the computation. In our example of figure 16, we see satisfactory numerical scanning, because the curve is a regular and smooth.

Fig.16: The observation with high time-resolution allows us to evaluate the quality of the numerical iteration.

Although we permanently extract power from the system (see figure 17), the total energy within the system is increasing strictly monotonously during the whole time of our observation (see figure 18). This will be the case as long as the rotor does not achieve the stable equilibrium-condition of constant angular-velocity. Remark: The situation is totally different, as soon as we extract mechanical power from the system, but we will see this later in section 6.



Fig.17: Electrical power extraction at the load resistor.

Fig.18: Total energy within the system. This is the energy of the coils, the capacitors and the rotation of the magnets. In our example we do not have to analyse any power being supplied to the system, because the input-circuit is not existing at all, for there is no input of energy. The machine is a self-running engine (not an over-unity system).


A long-term analysis of the electrical power being extracted from the load resistor is shown in figure 19. The only aspect which was changed with regard to the short-term analysis reported before, is the total duration of the observation, which is now 40 seconds for the computation to figure 19. Obviously the extractable power at the load resistor converges to ZERO. This observation is confirmed very clearly, if we perform a further enhancement of the duration time of the observation. This arises of the question, whether the EMDR-system is capable as a ZPE-motor at all !

The answer will be YES (!), as we will soon see in section 6. The EMDR-motor needs some mechanical resistance (such as for instance friction or power extraction) that it can work as a ZPE-motor.

 Fig.19: Electrical power being extracted from the load-resistor in the Turbo circuit in a long-time analysis.


Let us discuss this now:

If we want to have the EMDR-converter in powerful operation, the machine needs a constant phase-difference between the electric oscillation in the LC-circuit and the mechanical rotation of the permanent magnet. There is an optimum value for this phaseshift, which has to be maintained during long-term operation, in order to maintain the conversion of ZPE-energy constantly. This phase shift is absolutely necessary, that the oscillating magnetic field made by the turbo-coil can accelerate the rotating magnet exactly in this moment, in which the turbo-coil-field is the most strong. But on the other hand the turbo-coil-field should decelerate the magnet at this moment at which it is the most weak.

Let us remind, that a constant field would accelerate and decelerate the magnet to the same amount, so that the angular velocity will alter during each period of rotation, but it will become the same after a full period of rotation. Differently from this, an oscillating field has the possibility to make the acceleration different from deceleration.

We can illustrate this explanation by a rather simple model, as following: The oscillation within the LC-circuit causes the energy within this circuit move periodically back and forth between the (energy of the) electrostatic field in the capacitor and (energy of the) the magnetic field in the coil. This is exactly happening twice per period, because the capacitor plates will be charged alternating "positive" and "negative", as well as the coil is alternating between "north" and "south" (at each side). At these moments, in which the field energy is inside the coil (i.e. the field energy is the energy of a magnetic field), the magnet has to be accelerated remarkably. But in the moment when the field is inside the capacitor (i.e. the field energy is the energy of an electrostatic field), the magnet is decelerated by the coil only very weak, because the coil has almost no magnetic field.

A graphical illustration can be seen in figure 20 and 21, which is an animation of eight pictures cyclically following one after each other as a function of time. Let us start our considerations with figure 20, where we see the case of a "beneficial phase-shift", driving a strong rotation. The capacitor is drawn with purple colour, and the electrical charge within the capacitor is illustrated by "positive" and "negative" algebraic signs noted at the capacitor-plates. The more charge we find inside the capacitor plates, the more algebraic signs are noted. When the field (and some electrical charge) moves into the oscillation-circuit, an arrow in black colour symbolises the direction of propagation of the field. Obviously the field has to pass the long wire, which is formed as a coil drawn in blue-colour. This has the consequence, that there is an electric current within the coil, producing a magnetic field. This magnetic field is noted next to the wire of the coil by the use of the symbols ("N" and "S"). Again the number of symbols represents the strength of the field.


Let us now begin our considerations with part "1" of figure 20. This is the moment t1, where the capacitor still contains some electrical charge, but it is far away from being charged to its maximum. This means that some of the field’s energy is in the coil, causing some magnetic fields (also away from its maximum). But we see the polarity of the magnetic field and the polarity of the permanent magnet, which are orientated relatively to each other in such a way, that the magnet will begin to rotate clockwise. In part "2", the capacitor is discharged completely, so that the electrical current within the coil reaches its maximum. Now the total field’s energy is the energy of a magnetic field, so that the attractive force accelerating the permanent magnet is rather strong. This means that the angular velocity of the rotation is enhanced remarkably (because the phase-difference between the fields in the LC-circuit and the position of the rotating magnet is "beneficial"). In part "3", the magnetic fields of the coil is reduced partially, but due to the orientation of the permanent magnet, there is no more torque and thus no angular acceleration on the magnet. In this situation, the magnet is simply continuing his rotation constantly. In part "4", there is no more magnetic field in the coil, because there is no electrical current and the coil. This means that there is no deceleration which might like to move the magnet back with its "south-pole" to the position where we formally had the “north-pole” of the coil. In this situation, the magnet is also continuing his rotation constantly. But where is the energy of the oscillating field in the LC-circuit ? It is inside the capacitor, so that it does not have the chance to have any influence on the motion of the magnet. This means, up to now we can say, that the magnet was accelerated but not decelerated. In part "4", the capacitor is charged up to its maximum. From now on the capacitor begins to discharge, which can be ovserved clearly in part "5". Of course the discharge-current has the opposite direction as in part "1", so that the magnetic field of the coil in part "5" has the opposite polarity as in part "1". And this is fine, because the magnet also has the opposite direction as in part "1". This means that the field of the coil begins to accelerate the rotation of the magnet also clockwise in part "5". In part "5" we have a magnetic field growing slowly, but in part "6", this magnetic field already reaches its maximum value, so that we have a good angular acceleration of the magnet which is again orientated clockwise. (The situation of "6" corresponds to the situation of "2".) Part "7" now corresponds to the situation of "3", and after the above explanation it is clear, that neither in "7" nor in "8" there is any deceleration disturbing the rotation of the magnet. The rotation is going clockwise, so that part "8" is followed by part "1", and so on…

The acceleration of the magnet is active as long as the phase shift between the orientation of the magnet and the magnetic field in the LC-circuit is present.

Fig.20: Rotation of a magnet in the coil of a LC-circuit with constant “beneficial phase-shift” between the orientation of the magnet and the orientation of the oscillating magnetic field in the circuit.


If the magnet would be accelerated permanently, due to the beneficial phase shift, the angular velocity of the rotation would increase until the rotation of the magnet catches the rotation of the LC-circuit. But this would have the consequence, that the phase-shift between the electric LC-oscillation and rotation would disappear. And then the phase-shift would no longer be “beneficial”, but it would come into a condition which we can see in figure 21. Under this condition there is no more acceleration of the magnet (i.e. no further conversion of vacuum-energy).

Let us begin our explanations with part "1" of figure 21. Here we still see a slight clockwise acceleration acting on the permanent magnet. The magnetic field of the coil is not extremely strong, but it is existing. Part "2" does not need very much explanation, because the orientation of the permanent magnet does not allow any acceleration at all. But interesting is now part "3" with counterclockwise acceleration of the permanent magnet. And as we see, the absolute value of the deceleration in "3" is of the same size as the absolute value of the acceleration in "1". This means that the acceleration in "1" is completely compensated by the deceleration in "3", indicating that the time from "1…3" does not cause any acceleration or deceleration in sum at all. Part "4" does not change our general situation, because here we again have no acceleration or deceleration. There is simply no magnetic field of the coil in part "4". The train of thoughts, which we applied from "1…4" can be repeated for "5…8" analogously, so that we come to the conclusion, that during one turn, the angular velocity of the rotating magnetic is increasing and decreasing periodically, but there is no sum acceleration or deceleration at all. There is no “beneficial” phase-shift between the rotation of the magnet and the oscillation of the LC-circuit.

From here we understand, that the EMDR-converter permanently needs a "beneficial phase-shift" between the electrical and magnetic motion. If this phase shift decreases to "zero", it is not further possible to convert any ZPE-energy. This means that we need the "beneficial phase-shift" for proper operation of the ZPE-motor. But the phase shift cannot be generated electrically (as for instance by a load-resistor). Thus we have to generate the "beneficial phase-shift" mechanically by applying some mechanical load to the rotating shaft at the axis of the magnet.

Fig.21: Rotation of a magnet in the coil of a LC-circuit without “beneficial phase-shift” between the orien­ta­tion of the magnet and the orientation of the oscillating field in the circuit.

Due to the absence of mechanical load, the rotating magnet can catch the position of the oscillating field, so that no further ZPE- energy will be converted. This way of operation brings us to the knowledge, that the EMDR-con­verter can be operated as a ZPE-con­verter only, if permanently some mechanical power is extrac­ted.


The fact that the rotating magnet can never overtake the oscillation of the LC-circuit (even not if it is free from any load), assures that there is an upper limit to the angular velocity of the mechanical rotation. This prevents us from the danger of too speedy rotation, which might cause a damage of the machine or an accident. Of course we have to make sure, that the rotating magnet can withstand the frequency of the LC-circuit.

In order to convert ZPE-energy, the EMDR-motor must have a positive "beneficial phase-shift", of which the ideal case is shown in figure 20, representing a phase-shift of 45°. For it is not possible to maintain this phase-shift only by electrical means, mechanical power extraction from the shaft is absolutely necessary. This means that the rotating axis of the permanent magnet has to be loaded permanently with some mechanical torque, in order to maintain the requested "beneficial phase-shift".

With regard to this aspect, our ZPE-motor shows completely different behaviour than a classical electrical motor. If we for instance decide to use good bearings in order to minimise friction and to avoid mechanical load to the shaft, we will even not be able to get any electrical power of the engine. But already if we begin to enhance mechanical friction in the bearings, this would help to enhance the amount of electrical power which can be extracted.

Of course, the application of bad bearings with much friction is not, what we recommend. Preferable is the well-controlled extraction of beneficial mechanical power from the rotating shaft. In the ideal case, we could have some active control, to influence the mechanical power extraction in such a way, that the "beneficial phase-shift" will always be regulated most close to its optimum value of 45°. This would allow long-term stable operation of the EMDR-motor with a maximum of extracted energy. This regulation can be constructed with regard to the angular velocity, or with regard to the phase shift. If we decide to choose the first alternative, we can define an "upper-level" and the "lower-level" for the angular velocity, and as soon as the rotating magnet becomes faster than the “upper-level”, the amount of the power being extracted can be enhanced. Analogously, the amount of power being extracted can be reduced, as soon as the angular velocity becomes slower than the "lower-level". This is the basic idea on which the following chapter 6 is constructed.


6. The EMDR-Converter with mechanical power-extraction

As we saw in the figures 13, 14, 15 and especially in figure 19, the conversion of zero-point-energy in our EMDR-system can work efficiently during the initialisation of the operation. But the converted power goes asymptotically down to zero, in the same way as the phase shift between the rotation of the magnet and the LC-oscillation-circuit goes asymptotically down to zero, if we do not have some special technique maintain some “beneficial phase-shift” remarkably different from zero. We know the reason from section 5.

Therefore we have two apply mechanical torque Mmech(t) to the shaft, in order to get a beneficial operation-mode of the EMDR-system. This torque Mmech(t) has to be brought into equation (24a), as additional contribution, additionally to the torque of the magnetic forces (between permanent magnet and coil). The power which is extracted from the shaft by this torque is written in equation (27), and we find it in the second last line of table 1.

We will see, that this will enable us to run the EMDR-motor long-term stable and extract remarkable power. The DFEM-simulation is very encouraging, and therefore it’s source-code is printed completely in the appendix of the publication here.

What we apply is a torque proportional to the angular velocity of the rotation (this is an arbitrary decision, it could also be made different), which we find in the source-code in those lines which begin with the comment "(GG)". The coefficient of friction has the name "cr" in the program, and it can be controlled with the subroutine "Reibung_nachregeln". This subroutine works with an angular velocity called "phipZiel", around which we have a small hysteresis, within which the angular velocity of the EMDR-motor shall be kept. In order to make the results convincing (with regard to the conversion of ZPE-energy), we allow the magnet to spin a little bit faster than with the initial frequency, so that everybody can see that also the angular velocity is enhanced during operation, fed by ZPE-energy.

As can be seen in the input-data of the source-code, the solution of the differential equation is made by numerical iteration with 108 steps, of which everyone has 43 nanoseconds. It was verified, that this is indeed a sufficient time-resolution, so that the algorithm has converged to the serious result.

When you run the algorithm, you get the following data onto the screen:

power gained within the system:    5.04845399 Watt
total extraction of energy on the load resistor = 223.50737922 Joule
corresponding to a power of:       5.19784602848813E+0001 Watt
energy being supplied from input:  0.00000000000000E+0000 Joule
corresponding to a power of:         0.00000000000000E+0000 Watt
totally extraction of mechanical energy =   2271.25431928806 Joule
corresponding to a power of =               528.19867890420 Watt
at a duration of observation of    4.30000000000000E+0000 sec. 


We interpret this as following:

The start of the EMDR-converter is initialised with an angular velocity of 30,000 rounds per minute, while there is no initial electrical energy within the LC-oscillation-circuit. The rotation of the magnet induces a voltage into the coil and thus brings electrical energy into the LC-oscillation-circuit. Of this reason, the initialisation of the operation extracts some energy from the rotation of the magnet and brings it into the LC-oscillation-circuit. We see this in figure 22. But even during the first oscillation we see the energy gain from the ZPE-energy, so that the angular velocity of the magnet will be soon faster than at the very beginning. There we see the transient behaviour of the system, during which the energy is distributed between the electrical and mechanical part of the system as the laws of physics require. The transient behavior finally leads us to the frequency "phipZiel" with precision of a small hysteresis due to the control mechanism of the "beneficial phase-shift".

Fig.22: Angular velocity of the rotating permanent magnet in a powerful EMDR-converter. The graphic is plotted from column "I" of the Excel data-export.


The control of the power extraction is being done via a time-dependent regulation of the coefficient of friction cr. The time-dependent dynamic behaviour of this coefficient is printed in figure 23. Obviously the transient motion at the very beginning of the operation requires some control of the coefficient. But is seems as if this control is not further necessary, when the stable mode of operation is achieved. (We will discuss this later in detail.)

Fig.23: Dynamic behaviour of the coeffi­cient of friction, as it is used for the control of mechanical power ex­tract­ion from the EMDR-motor. The graphic is plotted from column "X" of the Excel data-export.


The amount of power actually being extracted is plotted in figure 24. With regard to the size of the setup, a mechanical power-output of a bit more than 530 whites is okay - as the permanent magnet has a length of only 10 cm. And it shall be mentioned that there is additional electrical power-output and the machine (at the load resistor), which is not yet fully optimised with regard to power-output.

Fig.24: Mechanical power output from the EMDR-converter. The graphic is plotted from column "Y" of the Excel data-export.

In order to get an imagination about the mechanical and electrical dimensions, necessary for the design of the prototype, we now want to inspect the electrical current and the voltage in the LC-oscillation-circuit (see figure 25).

Fig.25: Electrical charge Q in the capacitor of the EMDR-con­verter of our example. The graphic is plotted from column "A" of the Excel data-export.

The capacitor-voltage goes rather quick from the transient behaviour to the stable operation. Figure 25 displays the envelope of the oscillating charge in the capacitor, from which we calculate the capacitor-voltage according to equation (28), as being nearly 200 Volts. This value is not a problem for a practical experiment.

But it is important to know, that the capacitor must be tuned extremely fine, because this is the device, with which the EMDR-converter is fine-tuned. The adjustment of its value with a precision of about 1% … 0.1% is necessary, so that it is a good advice to use a capacitor bank for instance as shown in figure 26. All single capacitors must be connected parallel and not in series. The capacitor should have almost the same internal resistance. And it is recommendable that they have the same time constants for being charged and discharged. Nevertheless they shall have different capacity in order to make the fine-tuning possible over a wide range of capacity.

Finally it should be emphasized, that the adjustment of the capacitor is the means, by which the uncertainties and the approximations of the theoretical calculations have to be compensated !

This means, that the computation can deliver a theoretical value of the capacity deviating from the real experimental value even by a factor of 1 … 2 … 3 (more or less ?).


Fig.26: Capacitor bank for fine-tuning of the capacity. All single capacitors may have different capacity, so that a wide range of capacity can be tuned to this very fine adjustment.

An electrical current of about 60 Amperes (see figure 27) is not very much for the coil we have used, which has a diameter of the wire of 10 mm (i.e. a cross section 78.5 mm²). The problem is the handling of such a thick wire, but a trained mechanician should be able to do this. The capacitor bank must withstand a current of 60 Amperes.

Fig.27: Electrical current dQ/dt in the LC-oscillation-circuit of our EMDR-con­verter. The graphic is plotted from column "B" of the Excel data-export.


The voltage over the coil can be found by equation (29) from the law of induction. It’s amplitude is also close to 200 Volts, and thus far away from causing difficulties.

Fig.28: Time derive of the current in the coil, i.e. d2q/dt2 in our EMDR-con­verter. The graphic is plotted from column "C" of the Excel data-export.


The angular acceleration of the rotating magnet (column "J" in Excel) with an amplitude of nearly 1000 rad/s2 makes us understand, that the components of the EMDR-con­verter have to be mounted with appropriate mechanical stability. This angular acceleration acts on the magnet with a ponderable mass of 245 Gramms and a moment of inertia of rotation of 2.1·10-4 kg·m2. If the coil and the axis of rotation is not mounted with appropriate stability, the experiment might be dangerous.

The total energy of the system remains rather close to the value to which it is regulated (keep hysreresis in mind), as we see in figure 29. Similar behaviour can be seen in all channels of the energy analysis of the EMDR-system, because in the stable and durable long-term mode of operation, the energy is just oscillating between the different components of the system. And the amount of power per time, which is gained from the ZPE-energy is immediately converted into mechanical energy per time, being extracted at the shaft.



Fig.29: Total energy in the EMDR-system. The graphic is plotted from column "T" of the Excel data-export.

Additional question:

Is it allowed to drive the EMDR-motor without elaborated regulation of mechanical power extraction?

This would make it much easier to build up a prototype.


Yes, it is allowed. It is even not the problem at all.

Also when the coefficient of friction is simply constant (in algorithm, const: „cr = crAnfang“), the EMDR-motor will run into a stable mode of operation by alone, as long as the coefficient of friction is within certain limits, which are not extremely narrow. If for instance we run the algorithm as shown in the appendix, but let us apply one single change, namely to use a constant coefficient of friction without any regulation, applying the value “crAnfang:=37E-6”, the motor will come rather quick to a stable mode of operation, as we can see it in figure 30. Also the mechanically extracted power will come to a constant value rather soon, as we see in figure 31.



Fig.30: Angular velocity of the rotating permanent magnet. The graphic is plotted from column „I“ of the Excel data-export.

Fig.31: Mechanical power extraction from the rotating shaft. The load-torque is proportional to the angular velocity of the axis. The graphic is plotted from column „I“ of the Excel data-export.

Experimentalists, who build up an EMDR-converter with simple constant friction cr, should be careful not to use to strong friction, because strong mechanical load during the transient (initialization-) phase of operation does extract too much energy, so that the system cannot come into the long-term stable mode of operation. The allowed mechanical load without regulation is smaller than the allowed mechanical load with regulation. This means, if we take the average value of the load-coefficient of figure 23 (this is cr:=crAnfang:=54E-6), the EMDR-converter would never come to a stable mode of operation, because it would be slowed down too much during the initial phase. If we apply "crAnfang:=54E-6" with good regulation of "cr", we can extract 537 Watts mechanically, but if we apply "cr:=54E-6=const" without any regulation, EMDR-converter does work as a self-running ZPE-engine.

Here are some examples for load-coefficients recommendable or not recommendable:

Here are some examples for load-coefficients recommendable or not recommendable:
cr:=crAnfang:=2.0·10-5  =>  Pmech = 201 Watt    (goes to stable operation)
cr:=crAnfang:=2.5·10-5  =>  Pmech = 251 Watt    (goes to stable operation)
cr:=crAnfang:=3.0·10-5  =>  Pmech = 300 Watt    (goes to stable operation)
cr:=crAnfang:=3.5·10-5  =>  Pmech = 349 Watt    (goes to stable operation)
cr:=crAnfang:=3.7·10-5  =>  Pmech = 369 Watt    (goes to stable operation)
cr:=crAnfang:=4.0·10-5  =>  Pmech = 247 Watt    (no stable operation, angular velocity goes asymptotically down to zero)
cr:=crAnfang:=5.3·10-5  =>  Pmech = 117 Watt    (average of Fig.23, no stable operation)


Obviously the extractable mechanical power increases linearly with increasing load-coefficient, as long as the mechanical load is not to strong (see the range of  cr:=2.0·10-5    3.7·10-5). But if the load is too strong, the EMDR-engine will not have a chance to adjust the “beneficial phase-shift” between its electrical and its mechanical motion in such a way, that it can find its stable operation. This is, what we observe at cr:=4.0·10-5, and of course much more at cr:=5.3·10-5.

So we can say, that the control of the load coefficient is not absolutely necessary, especially not for the first prototypes, but it is nice to have for EMDR-motors build later, in order to maximise the extractable power in technical application. This is, what we got from regulation of the load-coefficient:

initial phase:    crAnfang:=4.5·10-5,

after regulation :      cr:=5.4·10-5, =>  Pmech = 537 Watt    (long-term operation, figure.23)

The reason is, that the load is regulated down in those moments in which the torque is weak (during initialization), and the load is regulated up during those moments in which strong torque and power is available.


7. Practical advice for experimenalists, who want to build an EMDR-Converter

With the end of section 6, the theory of the EMDR-system is discussed. But we now want to speak about practical aspects, which are interesting for those, who want to try to verify the theory experimentally [PC 11].

Central part of the EMDR-converter is the magnet. But the algorithm allows to simulate almost every available magnet or configuration of magnets. For the sake of simplicity, our calculation-example was done with the simple cylindrical bar-magnet, which should be available very easily. In [Tur 11], the calculation-example was restricted to a homogeneous magnetic field of a disc magnet, with an orientation of the magnetisation "in plain" (in order to make the computations easy). This restriction is now not necessary any further, with the new algorithm presented here. Thus the calculation-example in [Tur 11] has been a rather rough approximation, but the accuracy now is better.

Experimentalists who use a cylindrical bar-magnet should not forget about the aerodynamic drag of the rotating magnet due to its rather high angular velocity. We can accept this aerodynamic drag as a mechanical load which helps to maintain the "beneficial phase-shift" according to figure 20, when we will try to build the first prototypes.

But this is not the version which shall be used for future technical applications. For technical applications, engineers will have to try to maximise the output-power, and thus the aerodynamic drag should be minimized. This can be for instance done by inserting the cylindrical bar-magnet into a round disk, as shown in figure 32. But the disc must not be made of ferromagnetic material, because such a disc would guide the magnetic field flux lines into the totally wrong direction, or even keep the magnetic field flux lines inside the disc. And probably it is recommendable, to manufacture the round carrier-disc from nonconductive material in order to avoid eddy-currents carefully. Thus it is a good idea, to manufacture the round carrier-disc from plastic, and to use a type of plastic which is mechanically stable, so that it can withstand the centrifugal forces.



Fig.32: Possible suggestion to encapsulate a cylindrical bar-magnet in order to minimise the aerodynamic drag due to the rotation of the magnet.


Regarding eddy-currents: Up to now, there are no investigations about the question, whether they disturb the EMDR-motor or not. There is even no theoretical analysis of the role of eddy-currents up to now. There have been hints in discussions with colleagues, that eddy-currents would prevent the EMDR-motor to run at all, but perhaps eddy-currents might perhaps define a mechanical load and thus help to maintain the "beneficial phase-shift". It should be part of the experimental investigations to find out, whether eddy-currents disturb the EMDR-system or not. But for the first experimental trials and approaches, I recommend to avoid eddy-currents, because the theory works fine without eddy-currents.


Several colleagues discussed about the question, whether rare-earth magnets (such as neodymium) are mechanically stable enough to withstand the centrifugal-forces at 30,000 rounds per minute. It is well-known that neodymium magnets are less mechanically stable than steel. Nevertheless an encapsulation with steel in order to enhance the mechanical stability is absolutely forbidden, due to its influence on the field flux lines. For those, who want to maximise the mechanical stability, it is recommended to use an Iron-Cobalt-Nickel alloy for the magnet (without any encapsulation). It is possible to make a cylindrical bar-magnets of such alloy with a magnetic field-strength of 1 Tesla at each end of the bar. Such field-strength is absolutely sufficient for the operation in the EMDR-system. And furthermore the application of magnets of such alloy allows the mechanical stability of steel without encapsulating the magnet at all.

For sure the problem of eddy-currents is not neglectable for any encapsulation. We can see this, when we regard the energy-transfer from the rotating magnet into the coil as some special type of eddy-current-loss, which we need for the operation of the EMDR-converter.

Nevertheless, eddy-current-losses in the axis of rotation (we speak about the axis necessary to keep the magnet in its rotating-position), should not be regarded as a very large problem. It was discussed to use special glassfibre-plastic (for instance "GFK") in order to get an axis without any eddy-current-losses. There is no argument against such an axis, as long as it is stable enough. But due to the stability of such a material, it might be recommendable, to mount the bearings rather close to the magnet (as drawn in figure 11), which makes it difficult, to apply the initial angular velocity as well as to extract mechanical power during operation. Furthermore, if eddy-current-losses are the reason to use an axis made of glassfibre-plastic, the bearings, which are rather close to the magnet and inside the coil, should also be of some isolating material, as for instance ceramic-ball-bearings. But such bearings are not stable for the same angular frequency as steel bearings. This is the reason, why I do not want to exclude the use of non-ferromagnetic metallic axis, which is long enough, to mount the steel-bearings far away from the magnet and the coil (i.e. outside the coil). Finally the experiment will have to decide between these different suggestions for the setup - some of them might work and some others not.


Also a second coil has been under discussion, which could be brought to the position in which we have the input-coil, but which has the purpose is to extract energy, namely electrically induced energy due to the rotation of the magnet. This application would change the name of the “input-coil” into “output-coil”. But the algorithm in the appendix makes it rather easy to add such an output-coil into our EMDR-model. Such an output-coil would have the advantage, that we can adjust the number of windings to the requirements of the voltage-current-characteristics, which will be preferable for future applications. (The output-coil can also be mounted perpendicular to the turbo-coil in order to allow both of them to come most close to the rotating magnet.)


A further suggestion regards the reduction of the angular velocity of the rotation of the magnet, namely the use of a multipole-magnet (of higher order), as it can be seen for instance in figure 33. There we have 16 bar-magnets (dipoles) mounted around a wheel, so that the whole magnetic setup has a multi-polarity of 16, instead of 2, as we have it with simple cylindrical bar-magnet. The consequence is, that such a multipole-wheel will have 16 changes of the magnetic polarity within every turn instead of 2. An enhancement of the number of polarity-changes per turn by a factor of 8 means, that we can reduce the angular velocity (of the magnet) by a factor of 8 (in comparison with the simple cylindrical bar-magnet). A simple numerical example shows the advantage: If we for instance have an electrical LC-oscillation-circuit with a resonance frequency of the 32,000 rounds per minute, the mechanical rotation of the multipole-wheel only requests 4000 rounds per minute. The numerical example can be changed arbitrarily, as long as the number of the multipoles is even (not odd).

Abb.33: Multipole-Magnet, which is manufactured by mounting several magnetic dipoles around a wheel.

Nevertheless it should be pointed out again, that a dipole-magnet has to be mounted mechanically very stable, same as a multipole-magnet of higher-order in order to avoid accidents. But the situation is less critical for multipole-magnet of higher order, because it runs more smoothly than dipole magnet (similar as every electrical engine).


8. Resumée

The result of the work presented here, is a computation method for ZPE-motors, which is much more close to reality than the very first development in [Tur 11], which explained the very principle more than the engineering aspects. The algorithm in the appendix is written in a way, that every trained researcher should be able to use it for the simulation of his own ZPE-converter, which can be the EMDR-design (invented by the author of this article) or which can also be some other design.

Important part of the work presented here, is the design of concrete example for a self-running ZPE-motor (this is the EMDR-engine). The design is calculated and suggested practicable enough, that every trained experimentalist should be able to build it up, as soon as he or she has the possibility to work in an appropriately equipped laboratory. Of course there are still several open experimental questions to be solved, but I am convinced, that this should not be an existential problem. (It might be a problem of time, and thus I do not dare any predictions.)

The author of this article would like to build up such a prototype by himself very much (especially on the background that he is educated/trained as an experimental physicist), but unfortunately he does not have any possibility to work in a laboratory in the moment now.


9. Literature References
[Ans 08] Swanson, J. (1970-2008). Finite Element Program ANSYS,  ANSYS, Inc. Software Products, http://www.ansys.com
[Bea 02] Bearden, T. et. al. (2002) Motionless Electromagnetic Generator, US Patent, 6,362,718 from 26. März 2002, Inventors: Patrick L. Stephen, Thomas E. Bearden, James C. Hayes, Kenneth D. Moore, James L. Kenny, To be found at http://www.cheniere.org/
[Bor 99] Borland Pascal (Delphi 5 from 1999 or newer version)
[Bro 08] Bronstein, I. N., Semendjajew, K. A. , Musiol, G. , Mühlig, H. (2008), Taschenbuch der Mathematik, Verlag Harri Deutsch, 7.Auflage, ISBN 978-3-8171-2017-8
[Dub 90] Beitz, W., Küttner, K.-H. et. al. (1990), Dubbel - Taschenbuch für den Maschinenbau , 17.Auflage, Springer-Verlag. ISBN 3-540-52381-2
[Ger 95] Vogel, H. (1995), Gerthsen Physik, Springer Verlag. ISBN 3-540-59278-4
[Hoh 11] Hohl, D. (2011). http://www.magnetmotor.at/   In Internet many self-running magnet-motors can be found. Here is one example.
[Hur 40] Hurst, R. (1940). The Invention of Hans Coler, Relating To An Alleged New Source Of Power., B.I.O.S. Final Report No. 1043, B.I.O.S.Trip No. 2394, B.I.O.S. Target Number: C31/4799, British Intelligence Objectives Sub-Committee
[Jac 81] Jackson, J. D. (1981), Klassische Elektrodynamik, Walter de Gruyter Verlag, ISBN 3-11-007415-X
[Jeb 06] Jebens, K. (2006). Die Urkraft aus dem Universum, Jupiter-Verlag. ISBN 3-906571-23-8
[Kep 10] Keppe, N. and Pacheco, C. B. S. (2008). Keppe Motor Manual 1.0.   http://www.keppemotor.com/ and http://www.keppemotor.com/manual1.php
[Koh 96] Kohlrausch, F., et.al. (1996). Praktische Physik, 3 Bde., Verlag B.G.Teubner Stuttgart, 1996, ISBN-13: 978-3519230014
[Mar 88-98] Marinov, S. (1988-1998).
  - Der Kugellager-Motor und der Huber-Effekt, Stefan Marinov, raum&zeit 32 (1988) 81-84
  - Repetition of Silvertooth's experiment for measuring the aether drift", Stefan Marinov, Speculations in Science and Technology 12 No.3 (1989) 187-179
  - Propulsive and Rotating Ampère Bridges and the Principle of Relativity, Stefan Marinov, Physics Essays 4 No.1 (1991) 30-36
  - The Missing Magnetic Force Law, Galilean Electrodynamics 9 No.2, Stefan Marinov, (March/April 1998) 35-37
[Mie 84] Mielordt, S. (1984). Kompendium Hypertechnik. Tachyonenenergie, Hyperenergie, Antigravitation.   Berlin, 1984, Reprint, 4. Auflage vom raum&zeit Verlag, ISBN 3-89005-005-0
[Nie 83] Nieper, H.A. (1983). Konversion von Schwerkraft-Feld-Energie. Revolution in Technik, Medizin, Gesellschaft. MIT-Verlag, Oldenburg, 1983, 4. erw. Auflage, ISBN 3-925188-00-2
[PC 11] private communication, with numerous colleagues – too many people to be enumerated here.
[Stö 10] Stöcker, H. (2010). Taschenbuch der Physik, Horst Stöcker, Verlag Harri Deutsch,  ISBN 978-3-8171-1860-1
[Tur 10a] Turtur, C.W. (2010). Example of a simple Algorithm for the Construction of Zero-point-energy Converters, PHILICA.COM, ISSN 1751-3030, Article number 207, (9. Okt. 2010)
[Tur 10b] Turtur, C.W. (2010). DFEM-Computation of a Zero-point-energy Converter with realistic Parameters for a practical setup, PHILICA.COM, ISSN 1751-3030, Article number 213, (7. Dez. 2010)
[Tur 11] Turtur, C.W. (2011). DFEM-Simulation of a Zero-point-energy Converter with realisable Dimensions and a Power-output in the Kilowatt-range, PHILICA.COM, ISSN 1751-3030, Article number 219, (7. Feb. 2011)

10. Appendix: Source-Code of the DFEM-Algorithm (EMDR_009i.dpr)   [Bor 99]

Hint: If somebody wants to run the program, he or she will find the source-code  for free download at the internet-page of the author.

Information about this Article
This Article has not yet been peer-reviewed
This Article was published on 3rd April, 2011 at 07:34:00 and has been viewed 20872 times.

Creative Commons License
This work is licensed under a Creative Commons Attribution 2.5 License.
The full citation for this Article is:
Turtur, C. (2011). Construction guidelines for a ZPE-Converter on the basis of realistic DFEM-Computations. PHILICA.COM Article number 233.

<< Go back Review this ArticlePrinter-friendlyReport this Article

Website copyright © 2006-07 Philica; authors retain the rights to their work under this Creative Commons License and reviews are copyleft under the GNU free documentation license.
Using this site indicates acceptance of our Terms and Conditions.

This page was generated in 0.3612 seconds.