Claus W. Turtur (Fachbereich Elektrotechnik, University of Applied Sciences BraunschweigWolfenbuettel)
Published in physic.philica.com Abstract In [Tur 11] the theory of a powerful vacuumenergy 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 DFEMalgorithm in order to simulate a zeropointenergy (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 vacuumenergy motor according to his own conception. Article body
Table of contents 1. Physical fundament and preliminary work 2. Motion of the components of the ZPEConverter 3. Evaluation of the results of a converter example 4. Computation example for a concrete ZPEmotor 5. A concret EMDR vacuum energy converter 6. The EMDRConverter with mechanical powerextraction 7. Practical advice for experimenalists, who want to build an EMDRConverter 8. Resumée 9. Literature References 10. Appendix: SourceCode of the DFEMAlgorithm 1. Physical fundament and preliminary work The algorithm is designed to simulate electric and magnetic ZPEmotors by principle, and it is not restricted to one special design or setup. Thus it allows simulating ZPEmotors with arbitrary position and numbers of coils as well as arbitrary positions and numbers of magnets. Even electrostatic ZPEmotors 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 DFEMalgorithm 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 Lorentzforces, 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 threedimensional Cartesian coordinate system (blue colour), we see two coils oriented parallel to the yzplain (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 zaxis.The number and the arrangement of the coils can be chosen arbitrarily in the DFEMalgorithm. The example as being shown here, corresponds to the sourcecode in the appendix of the publication. In order to prepare the solution of the differentialequations within the DFEMalgorithm, we need the following: (a.) The computation of the induced voltage, which the rotating permanent magnet induces in the coils, and (b.) The computation of the magnetic force Lorentzforce 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 Lorentzforce (between these currents and the permanent magnet), and on the other hand the electrical current in the LCoscillationcircuit 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 zaxis), 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 dataarray. 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 sourcecode 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 zaxis is realised with a coordinatetransformation as usual (see equation 2). The angle φ describes the orientation of the lengthsymmetryaxis of the magnet relatively to the yaxis. (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 areaelement, 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 orthogonalvectors on all area elements of each coil are orientated exactly into the direction of the xaxis, so that the scalarproduct of the field with orthogonalvectors of the coilareaelements can be derived as simple as shown in equation (3), namely as the xcomponent of the flux.
The magnetic flux through the whole coil is then being calculated as the sum of all magnetic fluxelements through all finite areaelements forming the coil, so that the total flux follows equation (4).
Remark regarding the areaelements 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 conductorloop elements are defined by the geometrical connections between support point and support point. Finite areaelements of the coils fill up the area surrounded by the conductorloops, so that the magnetic flux through each coil can be calculated as the sum of the finite magnetic flux elements through all areaelements 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 timedependency 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 computationmethod, 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 barmagnet 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 angularvelocity 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 DFEMalgorithm. In order to smooth the numerical noise of the voltagesignal, a Fourierseries 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 highorder components we take, the smoother the signal we get. An additional effect of the Fourierseries is, that it helps to save CPUtime 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 timeconsuming operations (see equations 2 and 3), which take much more CPUtime, then the calculation of the Fourierseries with only five rather simple expressions to be summed up. This is important, because the solution of the differentialequation will have to be done with very many timesteps of few nanoseconds, so that the computation of the magnetic flux has to be done several 10^{8} or 10^{9} during each run of the DFEMalgorithm. 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 Fourierseries approximation is done in rather low order (N_{0}≤5) allows us to determine the Fouriercoefficients 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 Fouriercoefficients A_{ν}. Details of b: For the determination of the Lorentzforces, 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 barmagnets, because they are easy to buy and not very expensive, with regard to experiments. For the sake of simplicity, the cylindrical barmagnets are emulated by two circular conductor loops, located at each end of the cylindrical bar. This means, the magnetic field of the barmagnets 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 fieldstrength can be adjusted by the choice of the electric current in the conductor loops emulating the barmagnet. The calculation of the magnetic field of the emulationconductorloops is done with the use of BiotSavart’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 BiotSavart, 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 barmagnet, which is orientated along the yaxis. Thus the loop is orientated parallel to the xzplane, has the radius „MEro“, and is located at the yposition „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 q_{1} in (8). This has to be determined from the current I in the coil and the propagationvelocity 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 q_{1} 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 q_{1}=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 fieldelements 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 fieldstrength (see equation 13) along the axis of the coil (which here is the yaxis), and the check confirms our results. Now our simulation of the magnetic field of a cylindrical barmagnet 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 Lorentzforce 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 Lorentzforce, which is the fundament for the calculation of the torque, which the permanent magnet experiences. Therefore we again use BiotSavart’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 Lorentzforce, 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 loopelement, which experiences a Lorentzforce. 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 areaelements 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 BiotSavart’s law. s=(s_{x},s_{y},s_{z})=conductorloop element, r=(r_{x},r_{y},r_{z})=point of interest, ds=(ds_{x},ds_{y},ds_{z})=motion of the electrical charges in the conductorloop element. With regard to the parameters according to figure 5, we can write BiotSavart’s law according to equation (14): As usual, we put the expression of (15) into equation (14) to calculate the finite fieldelements 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 magnetemulationcoils. This calculation has to be done individually for each loopelements of the magnetemulationcoils. From there we come to Lorentzforce acting onto each magnetemulationcoil (see equation 16).
From all finite Lorentzforceelements dF we now calculate the finite torqueelements 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 DFEMalgorithm.
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 Fourierseries. 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 highfrequency components for sure (compare (5)). Again, the elapsed CPUtime to calculate the torque is rather large, so that the explicit torquecalculation is not recommendable within the solution of the differentialequation. Here we have again the advantage to save computer time due to the Fourierseries. 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 Fourierseries) 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 DFEMalgorithm presented here. This means that the algorithm can be used for the computation of every ZPEconverter working electrically or magnetically:  If the DFEMalgorithm is used to emulate several permanent magnets, their interaction can be analysed. This is an interesting application of the algorithm for the simulation of selfrunning ZPEmagnetmotors, as they can for instance be found at [Hoh 11], [Jeb 06].  If the DFEMalgorithm 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 ZPEconverters can be simulated, for instance as can be seen at [Mar 8898], [Bea 02].  If the DFEMalgorithm is used to emulate one coil and one permanent magnet, the “ElectroMechanic Double Resonance” converter (EMDRconverter) as proposed by the author of this article can be simulated [Tur 11].  If the DFEMalgorithm is used to emulate two coils and several permanent magnets fixed to each other in appropriate manner, the Keppemotor can be simulated [Kep 10].  If 6 cylindrical barmagnets are mounted within 6 coils and are arranged with each other in a hexagon, the Colerapparatus can be simulated. The behaviour of these elements will lead us to a rather complicated differentialequation, and some electrical elements forming an electric circuit are introduced as boundary conditions into the system of differentialequations. (For differentialequations, please see also section 2.) Perhaps, a simulation of the Colerapparatus might help to decide whether this motionlessconverter can work or not. [Hur 40], [Mie 84], [Nie 83]  Also dynamic input and output of energy is no problem, because the differentialequations describing the motion can be expanded with some input voltages, loadresistors, an so on… Such elements have to be taken into the differentialequations additionally. Also mechanical force of torque can be applied as boundary conditions in the differentialequations. In order to illustrate this, the sourcecode in the appendix contains two coils (see fig.5). In the source code, the left coil has of the name inputcoil on the right coil has the name turbocoil (see subroutine “U7”). Additionally there is a loadresistor („R_{Last}“) being connected with the turbocoil. The rotating axis of the permanent magnet is being supported initially with a given angular velocity, being applied as initialcondition for the solution of the differentialequation. 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 EMDRconverter is a self running motor, which needs energy support only at the very beginning of the motion to initialise the rotation. Furthermore the inputcoil is not active in the source code as printed in the appendix, because of EMDRconverter does not need permanent inputenergy, for it is a self running engine. Nevertheless the inputcoil 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 differentialequation this is simulated by a decelerating torque proportional to the angular velocity of the rotation (of the magnet). Such type of energyextraction can be used to simulate friction as well as to simulate energyextraction for technical application.  In order to develop an exemplary ZPEconverter and bring it into a stable permanent mode of operation, it is a convenient method, to control some inputpower 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 outputpower in the DFEMalgorithm 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 ZPEmotor can also work without a regulation, but with constand extraction of energy. The variability of the DFEMmethod is large, so that it is not restricted only to magnetic ZPEconverters, but it is also applicable to electrostatic ZPEconverters. The only necessary change is, to replace the Lorentzforce from equation (16) by the Coulombforce as seen in equation (18).
The precision of the calculation mainly depends on the precision of the inputdata, 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 ZPEConverter 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 differentialequations describing the motion in the ZPEconverter. 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 ZPEconverter can be described by the motion of its components. The adequate means for this description are the differentialequations based on the interactions of the components with each other. For ZPEconverters 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 differentialequations of higher order. If converters need inputenergy (such types which do not work as self running engines, but only as overunity systems), the energyinput has to be introduced as perturbationfunction in the differentialequations. For all types of interaction with some external elements components, the appropriate means is the introduction of perturbationfunctions into the differentialequations. This makes the higher order differentialequation systems inhomogeneous. Mathematically, this has the consequence, that we cannot simply derive an analytical solution, valid for each type of differentialequation system. Consequently, the central core of computation of the DFEMalgorithm 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 DFEMalgorithm needs some certain amount of CPUtime, if you want to have the solution with sufficient precision. You cannot expect the DFEMalgorithm to come to good convergence within few seconds. The understanding, how the differentialequations of the motions have to be formulated, is the central point of understanding, which every user of the DFEMalgorithm has to work on. Only on the basis of this understanding, the user can apply the DFEMalgorithm onto his own system. Let us now begin to develop the differentialequations for the example of the EMDRconverter, which define the core of computation of the algorithm. We will do this in analogy to [Tur 11]. The differentialequations by principle operate completely dynamically, so that the DFEMalgorithm 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,d^{2}q/dt^{2} and of φ, dφ/dt,d^{2}φ/dt^{2} , 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 timesteps in order to get reliable results. Consequently it will not be possible to save all data of all timesteps within special dataarrays. So the program now allows the calculation of as many time steps as required, but the datastorage will be done only for maximum 35,000 points which is a sensible upper limit for the dataexport 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 ZPEconverter for whose simulation we want to develop the DFEM differentialequations. It consists of two coils (red colour) and one rotating permanent magnet (black colour) and additionally an inputvoltage (blue colour) and one capacitor (green colour). For our EMDRconverter as our very example of a ZPEconverter can be operated as a self running motor, the inputvoltage is not applied at all, so that the left coil together with voltagesupply is taken away completely. Within the sourcecode, those lines for the simulation of the inputvoltagecomponents 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 EMDRconverter are only printed as comments in the sourcecode of the algorithm for those who want to use the algorithm for other ZPEmotors. These comments have only the purpose to help colleagues to simulate their ZPEmachines.
Fig.9: This very simple setup of an EMDRconverter already allows the powerful conversion of ZPEenergy 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 ZPEmotor. The development and formulation of the differential equation system is done rather similar as in [Tur 11]. (a.) The differential equation of a harmonic oscillation of an electric LCoscillation circuit is described by equation (21). (b.) The differential equation of an attenuated oscillation of an electric LCRcircuit is described by equation (22). Its numerical iterative solution, as being achieved with the solver of the differential equations in our DFEMalgorithm, 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 LCRoscillation circuit for the check of the differential equation (22) and the solution (23). (c.) 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 outputpower of our ZPEmotor. 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 wellknown 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. (d.) 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). (e.) On this basis, the numericaliterative 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 DFEMalgorithm performs the calculation of the motions of the components of the ZPEconverter, 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 signalshape are not allowed, because they do not follow the full dynamics of the differential equations. Those who want to adopt the DFEMalgorithm 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 LCoscillationcircuit (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 LCcircuit, 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 chargecarriers (such as for instance electrons or electron holes). There are two reasons explaining this argument: The first reason is, that the chargecarriers 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 chargecarrier 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 LCoscillation 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 capacitorplate 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 chargecarriers within the wire of the coil, but these displacements are rather small  by several orders of magnitude too small to transport chargefrom 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 onedimensional consideration makes it easier, and it is in good agreement with the onedimensional 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 airpressure, 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 “driftspeed” 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 LCoscillationcircuit. 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 DFEMcomputation is traced back to the conversion of ZPEenergy, as explained in detail in [Tur 10a] and [Tur 10b]. It shall help the readers of this publication to understand, how the LCoscillationcircuit is indeed used, to control the speed of propagation of the interacting fields, as being necessary for ZPEconversion. Now the principle our calculation is explained. The execution of the calculation can be seen in details in the algorithm in the sourcecode. 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 convertersystem 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 ZPEconverter. The table is of course made for the example of our EMDRconverter, 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 DFEMalgorithm is developed so far, that we can perform realistic computations of arbitrary electric and/or magnetic ZPEconverters. An example therefore shall be calculated in the following sections  namely with the EMDRconverter suggested by the author of this publication. 4. Computation example for a concrete ZPEmotor The namely reason for the development of the DFEMalgorithm presented here is the fact, that the author wants to show realistic computations on his EMDRConverter („ElectroMechanical DoubleResonance“ 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 parameterset 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 DFEMalgorithm 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 Inputparameters:)  epo:=8.854187817E12(As/Vm); (electric field constant)  muo:=4*pi*1E7(Vs/Am); (magnetic field constant)  LiGe:=Sqrt(1/muo/epo)(m/s); (speed of light) (For the solution of the differentialequations and for the display of the results:)  AnzP:=5000000; (number of time steps of the numerical iteration)  dt:=1E6; (Sec. ) (Duration of each single time step)  Abstd:=1; (only for preparation, do not alter the value)  PlotAnfang:=0000; (For Dataexport to Excel: First PlotPunkt)  PlotEnde:=5000000; (For Dataexport to Excel: Last PlotPunkt)  PlotStep:=200; (For Dataexport to Excel: step width of the data being exported) (Remark: Excel is restricted to maximal 32.767 Datagroups. 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 (DFEMmeshing 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 Inputcoil, left coil in figure .1)  Nturbo:=12; (number of windings of the Turbocoil, right coil in figure .1)  nebeninput:=8; (windings sidebyside in Inputcoil)  ueberinput:=10; (windings on top of each other Inputcoil)  nebenturbo:=3; (windings sidebyside in der Turbocoil)  ueberturbo:=3; (windings on top of each other Turbocoil) (Remark: Here the parameters are used to define rectangular coils according to Fig. 1. The crosssection of the Inputcoil consists of 8 windings sidebyside and 10 such layers on top of each other. The crosssection of the Turbocoil consists of 3 windings sidebyside 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:=1E2; (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 barmagnet according to Fig. 8 and Fig.9. The shape can be altered if required, but the meshingsubroutines have to be altered also.) (Remark regarding the datastorage 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 dataarray at finite geometrical steps. When the magnet is moving, the field is moving together with the magnet. The step length for the datastorage of the field is „Bsw“.) (Remark: The weird value for the electrical current in the magnetemulationcoils 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.35E8; (Ohm*m) (Specific electrical resistance of copper, [Koh 96])  rhoMag:=7.8E3; (kg/m^3) (Density of the magnetmaterial, Iron, [Koh 96])  CT:=36.61E6; (Farad) (Capacity within the oscillation circuit of the Turbocoil)  CI:=100E6; (Farad) (Capacity within the oscillation circuit of the Inputcoil) (Remark: In our example, the inputcoil has been modelled in order to prepare it for everybody who wants to have an additional coil in the DFEMalgorithm. But the inputcoil is not used for the computation of the converter, and the complete oscillation circuit containing the inputcoil is left away for the computation of the self running EMDRconverter.)  Rlast:=0.0111; (Ohm) (Ohm'ian load resistor in the Turbocircuit 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 – capacitorvoltage, current in the coil) (Remark: Even a self running ZPEmotor 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 ZPEconverter is not a self running engine, but only an overunity engine, it permanently needs some classical inputenergy for operation. This can be supplied mechanically (at the axis of rotation) or electrically with some inputvoltage. The last version is displayed in the sourcecode printed in the appendix, by the use of a subroutine with the name “U7”. Nevertheless this inputvoltage is not used for the solution of the differentialequations, because the EMDRconverter 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 Inputcoil)  DLT:=4*(yo+zo)*Spsw*Nturbo; (Meter) (Length of the wire of the Turbocoil )  RI:=rho*(DLI)/(pi/4*DD*DD); (Ohm) (Ohm`ian resistance of the wire of the Inputcoil)  RT:=rho*(DLT)/(pi/4*DD*DD); (Ohm) (Ohm`ian resistance of the wire of the Turbocoil)  BreiteI:=nebeninput*DD; HoeheI:=ueberinput*DD; (Width and height of the Inputcoil)  BreiteT:=nebenturbo*DD; HoeheT:=ueberturbo*DD; (Width and height of the Turbocoil)  fkI:=Sqrt(HoeheI*HoeheI+4/pi*2*yo*2*zo)/HoeheI; (Induktivitycorrection for short coil)  fkT:=Sqrt(HoeheT*HoeheT+4/pi*2*yo*2*zo)/HoeheT; (Induktivitycorrection for short coil)  LI:=muo*(2*yo+BreiteI)*(2*zo+BreiteI)*Ninput*Ninput/(HoeheI*fkI); (Geometrical average => Induktivity of the Inputcoil)  LT:=muo*(2*yo+BreiteT)*(2*zo+BreiteT)*Nturbo*Nturbo/(HoeheT*fkT); (Geometrical average => Induktivity of the Turbocoil)  omT:=1/Sqrt(LT*CT); (Resonanceangularfrequency of Turbocicuit of LT & CT)  TT:=2*pi/omT; (classical duration of oscillation of Turbocircuit 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 DFEMAlgorithm.) With these parameters, the ZPEmotor is modelled as being displayed in figure 11. As explained above, there is only the Turbocoil and no Inputcoil.
Fig.11: EMDRconverter 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 Turbocoil 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 crosssection of the coil has an inside area of 5 × 6 cm^{2}, and the thickness of the wire 10 mm. The time steps for the numerical iterative the solution of the differentialequations 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 DFEMalgorithm, 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 systemparameters 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 dataoutput 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 differentialequations 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 analysistime. If the energygain plus the energy being extracted from the system is positive, the system converts ZPEenergy into classical energy. If this energygain plus the energy being extracted is negative, the system converts in the opposite manner classical energy into ZPEenergy. 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 energyextraction.)  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 inputvoltage and the inputpower in connection with the inputsupply is printed on the computer screen, but it is ZERO, because there is no classical energy brought into the system during operation. The EMDRconverter is a self running system. The values are only displayed for those, who want to modify the DFEMalgorithm to simulate an overunity machine.  The total duration of the observation, which is the sum of all timesteps “dt”. We do not speak about the elapsed CPUtime, but we speak about the simulated time of operation during which the converter is running. The simple online dataevaluation 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 systemparameters 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 datafile for Excel, where they can be displayed graphically. The listing of the contents of the different Excelcolumns 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 ZPEenergy, 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 LCoscillation 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 ZPEconverter, 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 CPUtime 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 DFEMprogram in the appendix can only run, if there is a datafile in the same directory with the name 'schonda', which can be downloaded together with sourcecode of the DFEMprogram for free from the Internetpage of the author of the publication presented here. The datafile 'schonda' only has the purpose to save CPUtime, as following: During the phase of initialisation of the main program, all preliminary work takes some CPUtime, 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 geometryparameters 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 datafile '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 sourcecode in the appendix. The typical construction of a EMDRconverter 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 Hallprobe), its dimensions are put into the inputpata lines of the sourcecode. 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 Turbocoil and to put all the other requested data into the DFEMalgorithm. 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 Turbooscillationcircuit. 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. Highspeed 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 ZPEconverter, 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 ZPEpower 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 powergain is much larger than the electrical powergain. 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 LCcircuit 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 adjustmentquality 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 LCoscillationcircuit 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 Turbocoil 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 Turborotor 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 ZPEconverter 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 inputdata 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. DFEMSimulation of an EMDRmotor (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.000E0007 (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 Dataexport to Excel ) PlotEnde = 10000000 (first point for the Dataexport to Excel ) PlotStep = 400 (step width for the Dataexport to Excel ) (Definition of the both coils (turbo and input):) Spsw = 0.010000 (Meters: step width for the automatic meshgeneration 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 inputcoil) Nturbo = 9 (number of windings of the turbocoil) nebeninput = 10 (windings sidebyside of the inputcoil) ueberinput = 10 (layers of windings of inputcoil) nebenturbo = 3 (windings sidebyside of the turbocoil) ueberturbo = 3 (layers of windings of the turbocoil) Bsw = 1.0E0002 (spatial resolution of the computation of the magnetic field) MEyo = 5.00000E0002 (ycoordinates of the conductorloops for the emulation of the permanent magnet) MEro = 1.00000E0002 (Radius of the conductorloops for the emulation of the permanent magnet) MEI = 1.58998E+0004 (Electrical current in the conductorloops for the emulation of the permanent magnet) further technical dimensions: DD = 0.0100000 (Meters) (diameter of the wire for the coils) rho = 1.35000000000000E0008 (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.83000E0005 (Farad) (capacitor in the turbocircuit) CI = 1.00000E0004 (Farad) (capacitor in the inputcircuit) additional necessary values: Rlast = 6.400000E0002 (Ohm) (Ohm´s load resistor in the LCTurbocircuit UmAn = 30000.00 (U/min) (mechanical initial conditions  rpm of the rotating magnet) Uc = 0.00 (Volt) (electrical initial conditions  voltage at the turbocapacitor) 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 Inputcoil) DLT:=4*(yo+zo)*Spsw*Nturbo = 3.96000 (Meter, length of the wire oft he turbocoil) RI:=rho*(DLI)/(pi/4*DD*DD) = 0.00756 (Ohm) (Ohm`resistance of the wire of the Inputcoil) RT:=rho*(DLT)/(pi/4*DD*DD) = 0.00068 (Ohm) (Ohm`resistance of the wire of the turbocoil) BreiteI:=nebeninput*DD = 0.10000 (Width of the inputcoil) HoeheI:=ueberinput*DD = 0.10000 (height of the inputcoil) BreiteT:=nebenturbo*DD = 0.03000 (Width of the Turbocoil) HoeheT:=ueberturbo*DD = 0.03000 (height of the Turbocoil) 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.24238647244960E0001 (Induktivity, Inputcoil) LT:=muo*(2*yo+BreiteT)*(2*zo+Breitet)*Nturbo*Nturbo/(HoeheT*fkT) = 9.93606632469255E0004 (Induktivity, Turbocoil) omT:=1/Sqrt(LT*CT) = 3.19974964955735E+0003 (classical angular frequency of the Turbocircuit) TT:=2*pi/omT = 1.96364903362012E0003 (classical period of the Turbocircuit ) Mmag:=rhoMag*(pi*MEro*MEro)*(2*MEyo) = 0.245 kg (Mass of the Magnet) J:=Mmag/4*(MEro*MEro+4*MEyo*MEyo/3) = 2.10329628157837E0004 (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.10329628157837E0004 kg*m^2 total duration of observation: 2.00000000000000E+0000 sec. ExcelExport: 0.00000E+0000… 2.00000E+0000 sec., Step 8.00000E0005 sec. These are 25000 datapoint 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 CPUtime is much larger on a normal computer). During this timeinterval, 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 nonclassical source which we regard as invisible for classical eyes. According to former explanations, it should be the energy of the electromagnetic zeropoint 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 Exceldata, 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 voltageamplitude is to be found at U_{c,max}=Q_{max}/C_{turbo}=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 loadresistor 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 Turbocoil) 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 (crosssection 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 Lorentzforces between the magnet and the coil of the LCcircuit. 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 observationtime, 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 highresolution 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 timezoom 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 computertime 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 timeresolution 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 equilibriumcondition of constant angularvelocity. 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 inputcircuit is not existing at all, for there is no input of energy. The machine is a selfrunning engine (not an overunity system). A longterm 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 shortterm 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 EMDRsystem is capable as a ZPEmotor at all ! The answer will be YES (!), as we will soon see in section 6. The EMDRmotor needs some mechanical resistance (such as for instance friction or power extraction) that it can work as a ZPEmotor.
Fig.19: Electrical power being extracted from the loadresistor in the Turbo circuit in a longtime analysis. Let us discuss this now: If we want to have the EMDRconverter in powerful operation, the machine needs a constant phasedifference between the electric oscillation in the LCcircuit and the mechanical rotation of the permanent magnet. There is an optimum value for this phaseshift, which has to be maintained during longterm operation, in order to maintain the conversion of ZPEenergy constantly. This phase shift is absolutely necessary, that the oscillating magnetic field made by the turbocoil can accelerate the rotating magnet exactly in this moment, in which the turbocoilfield is the most strong. But on the other hand the turbocoilfield 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 LCcircuit 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 phaseshift", 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 capacitorplates. 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 oscillationcircuit, 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 bluecolour. 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 t_{1}, 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 phasedifference between the fields in the LCcircuit 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 "southpole" to the position where we formally had the “northpole” 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 LCcircuit ? 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 dischargecurrent 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 LCcircuit is present.
Fig.20: Rotation of a magnet in the coil of a LCcircuit with constant “beneficial phaseshift” 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 LCcircuit. But this would have the consequence, that the phaseshift between the electric LCoscillation and rotation would disappear. And then the phaseshift 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 vacuumenergy). 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” phaseshift between the rotation of the magnet and the oscillation of the LCcircuit. From here we understand, that the EMDRconverter permanently needs a "beneficial phaseshift" between the electrical and magnetic motion. If this phase shift decreases to "zero", it is not further possible to convert any ZPEenergy. This means that we need the "beneficial phaseshift" for proper operation of the ZPEmotor. But the phase shift cannot be generated electrically (as for instance by a loadresistor). Thus we have to generate the "beneficial phaseshift" 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 LCcircuit without “beneficial phaseshift” between the orientation 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 EMDRconverter can be operated as a ZPEconverter only, if permanently some mechanical power is extracted. The fact that the rotating magnet can never overtake the oscillation of the LCcircuit (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 LCcircuit. In order to convert ZPEenergy, the EMDRmotor must have a positive "beneficial phaseshift", of which the ideal case is shown in figure 20, representing a phaseshift of 45°. For it is not possible to maintain this phaseshift 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 phaseshift". With regard to this aspect, our ZPEmotor 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 wellcontrolled 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 phaseshift" will always be regulated most close to its optimum value of 45°. This would allow longterm stable operation of the EMDRmotor 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 "upperlevel" and the "lowerlevel" for the angular velocity, and as soon as the rotating magnet becomes faster than the “upperlevel”, 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 "lowerlevel". This is the basic idea on which the following chapter 6 is constructed. 6. The EMDRConverter with mechanical powerextraction As we saw in the figures 13, 14, 15 and especially in figure 19, the conversion of zeropointenergy in our EMDRsystem 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 LCoscillationcircuit goes asymptotically down to zero, if we do not have some special technique maintain some “beneficial phaseshift” remarkably different from zero. We know the reason from section 5. Therefore we have two apply mechanical torque M_{mech}(t) to the shaft, in order to get a beneficial operationmode of the EMDRsystem. This torque M_{mech}(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 EMDRmotor longterm stable and extract remarkable power. The DFEMsimulation is very encouraging, and therefore it’s sourcecode 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 sourcecode 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 EMDRmotor shall be kept. In order to make the results convincing (with regard to the conversion of ZPEenergy), 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 ZPEenergy. As can be seen in the inputdata of the sourcecode, the solution of the differential equation is made by numerical iteration with 10^{8} steps, of which everyone has 43 nanoseconds. It was verified, that this is indeed a sufficient timeresolution, 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 EMDRconverter is initialised with an angular velocity of 30,000 rounds per minute, while there is no initial electrical energy within the LCoscillationcircuit. The rotation of the magnet induces a voltage into the coil and thus brings electrical energy into the LCoscillationcircuit. Of this reason, the initialisation of the operation extracts some energy from the rotation of the magnet and brings it into the LCoscillationcircuit. We see this in figure 22. But even during the first oscillation we see the energy gain from the ZPEenergy, 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 phaseshift".
Fig.22: Angular velocity of the rotating permanent magnet in a powerful EMDRconverter. The graphic is plotted from column "I" of the Excel dataexport. The control of the power extraction is being done via a timedependent regulation of the coefficient of friction c_{r}. The timedependent 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 coefficient of friction, as it is used for the control of mechanical power extraction from the EMDRmotor. The graphic is plotted from column "X" of the Excel dataexport. The amount of power actually being extracted is plotted in figure 24. With regard to the size of the setup, a mechanical poweroutput 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 poweroutput and the machine (at the load resistor), which is not yet fully optimised with regard to poweroutput.
Fig.24: Mechanical power output from the EMDRconverter. The graphic is plotted from column "Y" of the Excel dataexport. 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 LCoscillationcircuit (see figure 25).
Fig.25: Electrical charge Q in the capacitor of the EMDRconverter of our example. The graphic is plotted from column "A" of the Excel dataexport. The capacitorvoltage 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 capacitorvoltage 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 EMDRconverter is finetuned. 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 finetuning 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 finetuning 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 LCoscillationcircuit of our EMDRconverter. The graphic is plotted from column "B" of the Excel dataexport. 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. d^{2}q/dt^{2} in our EMDRconverter. The graphic is plotted from column "C" of the Excel dataexport. The angular acceleration of the rotating magnet (column "J" in Excel) with an amplitude of nearly 1000 rad/s^{2} makes us understand, that the components of the EMDRconverter 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·m^{2}. 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 EMDRsystem, because in the stable and durable longterm 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 ZPEenergy is immediately converted into mechanical energy per time, being extracted at the shaft. Fig.29: Total energy in the EMDRsystem. The graphic is plotted from column "T" of the Excel dataexport. Additional question: Is it allowed to drive the EMDRmotor without elaborated regulation of mechanical power extraction? This would make it much easier to build up a prototype. Answer: 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 EMDRmotor 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:=37E6”, 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 dataexport. Fig.31: Mechanical power extraction from the rotating shaft. The loadtorque is proportional to the angular velocity of the axis. The graphic is plotted from column „I“ of the Excel dataexport. Experimentalists, who build up an EMDRconverter with simple constant friction c_{r}, 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 longterm 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 loadcoefficient of figure 23 (this is cr:=crAnfang:=54E6), the EMDRconverter 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:=54E6" with good regulation of "cr", we can extract 537 Watts mechanically, but if we apply "cr:=54E6=const" without any regulation, EMDRconverter does work as a selfrunning ZPEengine. Here are some examples for loadcoefficients recommendable or not recommendable: Here are some examples for loadcoefficients 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 loadcoefficient, 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 EMDRengine will not have a chance to adjust the “beneficial phaseshift” 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 EMDRmotors build later, in order to maximise the extractable power in technical application. This is, what we got from regulation of the loadcoefficient: initial phase: crAnfang:=4.5·10^{5}, after regulation : cr:=5.4·10^{5}, => Pmech = 537 Watt (longterm 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 EMDRConverter With the end of section 6, the theory of the EMDRsystem 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 EMDRconverter is the magnet. But the algorithm allows to simulate almost every available magnet or configuration of magnets. For the sake of simplicity, our calculationexample was done with the simple cylindrical barmagnet, which should be available very easily. In [Tur 11], the calculationexample 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 calculationexample in [Tur 11] has been a rather rough approximation, but the accuracy now is better. Experimentalists who use a cylindrical barmagnet 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 phaseshift" 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 outputpower, and thus the aerodynamic drag should be minimized. This can be for instance done by inserting the cylindrical barmagnet 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 carrierdisc from nonconductive material in order to avoid eddycurrents carefully. Thus it is a good idea, to manufacture the round carrierdisc 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 barmagnet in order to minimise the aerodynamic drag due to the rotation of the magnet. Regarding eddycurrents: Up to now, there are no investigations about the question, whether they disturb the EMDRmotor or not. There is even no theoretical analysis of the role of eddycurrents up to now. There have been hints in discussions with colleagues, that eddycurrents would prevent the EMDRmotor to run at all, but perhaps eddycurrents might perhaps define a mechanical load and thus help to maintain the "beneficial phaseshift". It should be part of the experimental investigations to find out, whether eddycurrents disturb the EMDRsystem or not. But for the first experimental trials and approaches, I recommend to avoid eddycurrents, because the theory works fine without eddycurrents. Several colleagues discussed about the question, whether rareearth magnets (such as neodymium) are mechanically stable enough to withstand the centrifugalforces at 30,000 rounds per minute. It is wellknown 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 IronCobaltNickel alloy for the magnet (without any encapsulation). It is possible to make a cylindrical barmagnets of such alloy with a magnetic fieldstrength of 1 Tesla at each end of the bar. Such fieldstrength is absolutely sufficient for the operation in the EMDRsystem. 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 eddycurrents is not neglectable for any encapsulation. We can see this, when we regard the energytransfer from the rotating magnet into the coil as some special type of eddycurrentloss, which we need for the operation of the EMDRconverter. Nevertheless, eddycurrentlosses in the axis of rotation (we speak about the axis necessary to keep the magnet in its rotatingposition), should not be regarded as a very large problem. It was discussed to use special glassfibreplastic (for instance "GFK") in order to get an axis without any eddycurrentlosses. 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 eddycurrentlosses are the reason to use an axis made of glassfibreplastic, the bearings, which are rather close to the magnet and inside the coil, should also be of some isolating material, as for instance ceramicballbearings. 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 nonferromagnetic metallic axis, which is long enough, to mount the steelbearings 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 inputcoil, 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 “inputcoil” into “outputcoil”. But the algorithm in the appendix makes it rather easy to add such an outputcoil into our EMDRmodel. Such an outputcoil would have the advantage, that we can adjust the number of windings to the requirements of the voltagecurrentcharacteristics, which will be preferable for future applications. (The outputcoil can also be mounted perpendicular to the turbocoil 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 multipolemagnet (of higher order), as it can be seen for instance in figure 33. There we have 16 barmagnets (dipoles) mounted around a wheel, so that the whole magnetic setup has a multipolarity of 16, instead of 2, as we have it with simple cylindrical barmagnet. The consequence is, that such a multipolewheel will have 16 changes of the magnetic polarity within every turn instead of 2. An enhancement of the number of polaritychanges 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 barmagnet). A simple numerical example shows the advantage: If we for instance have an electrical LCoscillationcircuit with a resonance frequency of the 32,000 rounds per minute, the mechanical rotation of the multipolewheel 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: MultipoleMagnet, which is manufactured by mounting several magnetic dipoles around a wheel.
Nevertheless it should be pointed out again, that a dipolemagnet has to be mounted mechanically very stable, same as a multipolemagnet of higherorder in order to avoid accidents. But the situation is less critical for multipolemagnet 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 ZPEmotors, 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 ZPEconverter, which can be the EMDRdesign (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 selfrunning ZPEmotor (this is the EMDRengine). 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. (19702008). 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 9783817120178 [Dub 90] Beitz, W., Küttner, K.H. et. al. (1990), Dubbel  Taschenbuch für den Maschinenbau , 17.Auflage, SpringerVerlag. ISBN 3540523812 [Ger 95] Vogel, H. (1995), Gerthsen Physik, Springer Verlag. ISBN 3540592784 [Hoh 11] Hohl, D. (2011). http://www.magnetmotor.at/ In Internet many selfrunning magnetmotors 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 SubCommittee [Jac 81] Jackson, J. D. (1981), Klassische Elektrodynamik, Walter de Gruyter Verlag, ISBN 311007415X [Jeb 06] Jebens, K. (2006). Die Urkraft aus dem Universum, JupiterVerlag. ISBN 3906571238 [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, ISBN13: 9783519230014 [Mar 8898] Marinov, S. (19881998).  Der KugellagerMotor und der HuberEffekt, Stefan Marinov, raum&zeit 32 (1988) 8184  Repetition of Silvertooth's experiment for measuring the aether drift", Stefan Marinov, Speculations in Science and Technology 12 No.3 (1989) 187179  Propulsive and Rotating Ampère Bridges and the Principle of Relativity, Stefan Marinov, Physics Essays 4 No.1 (1991) 3036  The Missing Magnetic Force Law, Galilean Electrodynamics 9 No.2, Stefan Marinov, (March/April 1998) 3537 [Mie 84] Mielordt, S. (1984). Kompendium Hypertechnik. Tachyonenenergie, Hyperenergie, Antigravitation. Berlin, 1984, Reprint, 4. Auflage vom raum&zeit Verlag, ISBN 3890050050 [Nie 83] Nieper, H.A. (1983). Konversion von SchwerkraftFeldEnergie. Revolution in Technik, Medizin, Gesellschaft. MITVerlag, Oldenburg, 1983, 4. erw. Auflage, ISBN 3925188002 [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 9783817118601 [Tur 10a] Turtur, C.W. (2010). Example of a simple Algorithm for the Construction of Zeropointenergy Converters, PHILICA.COM, ISSN 17513030, Article number 207, (9. Okt. 2010) [Tur 10b] Turtur, C.W. (2010). DFEMComputation of a Zeropointenergy Converter with realistic Parameters for a practical setup, PHILICA.COM, ISSN 17513030, Article number 213, (7. Dez. 2010) [Tur 11] Turtur, C.W. (2011). DFEMSimulation of a Zeropointenergy Converter with realisable Dimensions and a Poweroutput in the Kilowattrange, PHILICA.COM, ISSN 17513030, Article number 219, (7. Feb. 2011) 10. Appendix: SourceCode of the DFEMAlgorithm (EMDR_009i.dpr) [Bor 99] Hint: If somebody wants to run the program, he or she will find the sourcecode for free download at the internetpage of the author.
Information about this Article This Article has not yet been peerreviewed This Article was published on 3rd April, 2011 at 07:34:00 and has been viewed 21482 times. 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 ZPEConverter on the basis of realistic DFEMComputations. PHILICA.COM Article number 233. 
