Bifilar pancake coil overunity experiment

[ayeaye]

One thing, time has to be in ms. The power calculation script.

#Resistor resistance in ohms
R = 47.0
#Frequency in Hz
F = 1000.0
#Voltages are in mV, time is in ms

def next(separator):
  global pos
  begin = pos + 1
  pos = csv.find(separator, begin)
  if (pos == -1): pos = len(csv)
  if (csv[begin : pos] == ""): return 0.0
  return float(csv[begin : pos])

csv = ""
pos = n = 0
xu = ei = eo = 0.0
f = open("input.txt")
while (True):
  s = f.readline()
  if (len(s) < 1): break
  csv = s.strip()
  if (len(csv) == 0): continue
  if (csv[0] != "-" and not csv[0].isdigit()): continue
  pos = -1
  t = next(",") * 1000
  if (n == 0): xu = t
  if (n == 1): xu = t - xu
  n += 1
  vs = next(",") * 1000
  vr = next(",") * 1000
  vl = vs - vr
  pl = vl * vr / R / 1000
  if (pl >= 0): ei += pl
  if (pl < 0): eo += pl
f.close()
ei *= xu
print("Input power was %.4f uW" % (F * ei))
eo *= -1 * xu
print("Output power was %.4f uW" % (F * eo))

Trinket.

https://trinket.io/python/2ee90868b9

Other thing, next() now returns 0 when there was no number. This enables to get all values, like the following calculates your energy per cycle. It's in J if your instantaneous power was in W. I don't know the frequency to calculate power. Your instantaneous power was always positive, why?

def next(separator):
  global pos
  begin = pos + 1
  pos = csv.find(separator, begin)
  if (pos == -1): pos = len(csv)
  if (csv[begin : pos] == ""): return 0.0
  return float(csv[begin : pos])

csv = ""
pos = n = 0
xu = ei = eo = 0.0
f = open("input.txt")
while (True):
  s = f.readline()
  if (len(s) < 1): break
  csv = s.strip()
  if (len(csv) == 0): continue
  if (csv[0] != "-" and not csv[0].isdigit()): continue
  pos = -1
  for j in range(6): next(",")
  t = next(",")
  if (n == 0): xu = t
  if (n == 1): xu = t - xu
  n += 1
  ei += next(",")
f.close()
print("Energy per cycle " + str(ei * xu))

Trinket.

https://trinket.io/python/66a825269d

MDO3034

[ayeaye]

It may look strange, but it's the right way to overcome what is missing in Python, and make it straightforward. Using list with an element string Instead of string, makes it modifiable by function. Inkscape script and its trinket.

#Scales in s and V, scale is 50 units
XU = 20 * 1e-9 / 50
YU1 = 2.0 / 50
YU2 = 0.2 / 50

def nextt(sl, separator):
  endp = sl[0].find(separator)
  if (endp == -1): endp = len(sl[0])
  tok = sl[0][ : endp]
  sl[0] = sl[0][endp + 1 : ]
  return int(round(float(tok)))

def tolist(path, yu):
  n = 0
  list = []
  segments = [[0, 0, 0, 0]]
  nextt(path, ",")
  nextt(path, " ")
  while (True):
    x = segments[-1][0] + abs(nextt(path, ","))
    y = segments[-1][1] - nextt(path, " ")
    segments[-1][2] = x
    segments[-1][3] = y
    if (path[0] == ""): break
    segments.append([x, y, 0, 0])
  for j in range(segments[-1][2]):
    while (j >= segments[n][2]): n += 1
    lx = segments[n][2] - segments[n][0]
    ly = segments[n][3] - segments[n][1]
    ratio = float(j - segments[n][0]) / lx
    v = segments[n][1] + int(round(ly * ratio))
    list.append(v * yu)
  return list

path1 = [""]
path2 = [""]
f = open("input.txt")
for s in f:
  startp = s.find("d=\"m ")
  if (startp == -1): continue
  startp = s.find(" ", startp)
  endp = s.find("\"", startp)
  if (path1[0] == ""): path1[0] = s[startp : endp]
  path2[0] = s[startp : endp]
f.close()
list1 = tolist(path1, YU1)
list2 = tolist(path2, YU2)
length = min(len(list1), len(list2))
for j in range(length):
  s = "%1.3E" % (j * XU) + ","
  s += str(list1[j]) + "," + str(list2[j])
  print(s)

https://trinket.io/python/eeffa85b19

Calculation script and its trinket.

# Resistor resistance in ohms
R = 47.0
#Frequency in Hz
F = 1000.0
#Voltages are in mV, time is in ms

def nextt(sl, separator):
  endp = sl[0].find(separator)
  if (endp == -1): endp = len(sl[0])
  tok = sl[0][: endp]
  sl[0] = sl[0][endp + 1 : ]
  if (tok == ""): return 0.0
  return float(tok)

n = 0
csv = [""]
xu = ei = eo = 0.0
f = open("input.txt")
while (True):
  s = f.readline()
  if (len(s) < 1): break
  csv[0] = s.strip()
  if (len(csv[0]) == 0): continue
  if (csv[0][0] != "-" and not csv[0][0].isdigit()):
    continue
  t = nextt(csv, ",") * 1000
  n += 1
  if (n == 1): xu = t
  if (n == 2): xu = t - xu
  vs = nextt(csv, ",") * 1000
  vr = nextt(csv, ",") * 1000
  vl = vs - vr
  pl = vl * vr / R / 1000
  if (pl >= 0): ei += pl
  if (pl < 0): eo += pl
f.close()
ei *= xu
print("Input power was %.4f uW" % (F * ei))
eo *= -1 * xu
print("Output power was %.4f uW" % (F * eo))

https://trinket.io/python/839fcae20d

Calculation of power of the Tektronix sample data above, and its trinket.

def nextt(sl, separator):
  endp = sl[0].find(separator)
  if (endp == -1): endp = len(sl[0])
  tok = sl[0][: endp]
  sl[0] = sl[0][endp + 1 : ]
  if (tok == ""): return 0.0
  return float(tok)

n = 0
csv = [""]
xu = ei = eo = 0.0
f = open("input.txt")
while (True):
  s = f.readline()
  if (len(s) < 1): break
  csv[0] = s.strip()
  if (len(csv[0]) == 0): continue
  if (csv[0][0] != "-" and not csv[0][0].isdigit()):
    continue
  for j in range(6): nextt(csv, ",")
  t = nextt(csv, ",")
  n += 1
  if (n == 1): xu = t
  if (n == 2): xu = t - xu
  ei += nextt(csv, ",")
f.close()
print("Energy per cycle " + str(ei * xu))

https://trinket.io/python/126f63131a

[ayeaye]

Not really the topic here, but there are ways in Inkscape to convert oscilloscope traces to path, by not copying the path from the image completely by hand. I'm no expert in Inkscape, so the Inkscape experts forgive me, i can only tell the little that i know, please share your knowledge if you know some great things.

As much as i know, there is no way to convert a trace into a single line path automatically. If one knows any way to do that, please say, maybe there are some tools other than Inkscape that can do that.

In the Inkscape Path menu there is an item "Trace Bitmap". This and the "Stroke to Path" from the same menu, enable to convert traces into a circular paths. These will not be lines, they have some thickness.

At the right edge of the Inkscape screen there are buttons to select different modes of snapping. One is "Snap to paths". One can draw a zigzag path inside the circular path of the trace, so that the ends of the segments will be on that circular path. The other snap mode is "Snap from and to midpoints of line segments". This enables to draw path through the midpoints of the segments of that zigzag path.

This path will be a calculated approximation of the trace by a single path. This is not hand copied, but transformed from the original image. It may be too complicated though. There still has to be snapping to grid, and also considering all kind of imperfections of the image, a path copied by hand may much better match the trace than this semi-automatic transformation. But just to let everyone know, that there is such option.

A line from the beginning of the x axis may be drawn later. At the left side click on "Edit paths by nodes". Holding down the shift key, select the paths and the line segments to be joined. Then "Join selected nodes" at the upper edge of the screen, joins the paths into a single path.

[ayeaye]

Btw, what one thinks about the circuit on the figure below? What happens there, when the pulse is on, and after the pulse ends? Can one simulate it in any way? The caveat there is that these are not two separate coils, the two parallel lines indicate that these two coils are on the same core, so they cannot be simulated as two separate coils. That is to understand, what the capacitance between the windings really does.

[ayeaye]

This can be simulated with python, or c if one likes, if there is no existing simulator that can simulate this. Something like, solenoid  phi = mu0 * mur * n * i * a / l , what about toroid? The current is split, part goes to the left coil and part to the capacitor, magnetic flux is however the sum of the magnetic flux that both coil segments generate. Faraday's law of induction  e = - dphi / dt , which is voltage on the coil segment. Capacitor  i = c * dv / dt . One should start from induction, because this is the engine of change.

No, unity is not written into these equations, if there is unity or overunity, then it supposed to come from applying these equations.

Interesting? This question usually means, who will now do the dirty work. Why should i do?