## thread for segments files and parameters for simulation runs

Discuss how polywell fusion works; share theoretical questions and answers.

Moderators: tonybarry, MSimon

happyjack27
Posts: 1438
Joined: Wed Jul 14, 2010 5:27 pm
cuboctahedron, electron plasma core:

happyjack27
Posts: 1438
Joined: Wed Jul 14, 2010 5:27 pm
rjaypeters wrote:6-coil circular six coils:

...

Regular polygons? As good as I can do.
3 decimal places (1 mm) should be plenty precise, i would think. any preferred scaling parameter? charge/current ratio? (on that note, anybody know how to convert "voltage on the grid" to charge density (columbs/meter))?

happyjack27
Posts: 1438
Joined: Wed Jul 14, 2010 5:27 pm
i hate to say it rj, but from the simulation it looks like you've got at least one of the current loops backwards. the electrons spiral and go in the same direction which suggests you have a pair of facing (e.g. opposite) coils with the same polarity.

(my first attempt had problems like this too, btw.)

rjaypeters
Posts: 869
Joined: Fri Aug 20, 2010 2:04 pm
Location: Summerville SC, USA
Huh. I'll check, but I thought I was careful to maintain counter-clockwise as seen from outside! Apologies and thanks.

Edit:

Code: Select all

``````16
0.0, 1.127, 1.273, -0.431, 1.041, 1.273
-0.431, 1.041, 1.273, -0.797, 0.797, 1.273
-0.797, 0.797, 1.273, -1.041, 0.431, 1.273
-1.041, 0.431, 1.273, -1.127, 0.0, 1.273
-1.127, 0.0, 1.273, -1.041, -0.431, 1.273
-1.041, -0.431, 1.273, -0.797, -0.797, 1.273
-0.797, -0.797, 1.273, -0.431, -1.041, 1.273
-0.431, -1.041, 1.273, 0.0, -1.127, 1.273,
0.0, -1.127, 1.273, 0.431, -1.041, 1.273
0.431, -1.041, 1.273, 0.797, -0.797, 1.273
0.797, -0.797, 1.273, 1.041, -0.431, 1.273
1.041, -0.431, 1.273, 1.127, 0.0, 1.273
1.127, 0.0, 1.273, 1.041, 0.431, 1.273
1.041, 0.431, 1.273, 0.797, 0.797, 1.273
0.797, 0.797, 1.273, 0.431, 1.041, 1.273
0.431, 1.041, 1.273, 0.0, 1.127, 1.273
16
-1.273, 1.127, 0.0, -1.273, 1.041, -0.431
-1.273, 1.041, -0.431, -1.273, 0.797, -0.797
-1.273, 0.797, -0.797, -1.273, 0.431, -1.041
-1.273, 0.431, -1.041, -1.273, 0.0, -1.127
-1.273, 0.0, -1.127, -1.273, -0.431, -1.041
-1.273, -0.431, -1.041, -1.273, -0.797, -0.797
-1.273, -0.797, -0.797, -1.273, -1.041, -0.431
-1.273, -1.041, -0.431, -1.273, -1.127, 0.0
-1.273, -1.127, 0.0, -1.273, -1.041, 0.431
-1.273, -1.041, 0.431, -1.273, -0.797, 0.797
-1.273, -0.797, 0.797, -1.273, -0.431, 1.041
-1.273, -0.431, 1.041, -1.273, 0.0, 1.127
-1.273, 0.0, 1.127, -1.273, 0.431, 1.041,
-1.273, 0.431, 1.041, -1.273, 0.797, 0.797
-1.273, 0.797, 0.797, -1.273, 1.041, 0.431
-1.273, 1.041, 0.431, -1.273, 1.127, 0.0
16
0.0, 1.127, -1.273, 0.431, 1.041, -1.273
0.431, 1.041, -1.273, 0.797, 0.797, -1.273
0.797, 0.797, -1.273, 1.041, 0.431, -1.273
1.041, 0.431, -1.273, 1.127, 0.0, -1.273
1.127, 0.0, -1.273, 1.041, -0.431, -1.273
1.041, -0.431, -1.273, 0.797, -0.797, -1.273
0.797, -0.797, -1.273, 0.431, -1.041, -1.273
0.431, -1.041, -1.273, 0.0, -1.127, -1.273
0.0, -1.127, -1.273, -0.431, -1.041, -1.273
-0.431, -1.041, -1.273, -0.797, -0.797, -1.273
-0.797, -0.797, -1.273, -1.041, -0.431, -1.273
-1.041, -0.431, -1.273, -1.127, 0.0, -1.273
-1.127, 0.0, -1.273, -1.041, 0.431, -1.273
-1.041, 0.431, -1.273, -0.797, 0.797, -1.273
-0.797, 0.797, -1.273, -0.431, 1.041, -1.273
-0.431, 1.041, -1.273, 0.0, 1.127, -1.273
16
1.273, 1.127, 0.0, 1.273, 1.041, 0.431
1.273, 1.041, 0.431, 1.273, 0.797, 0.797
1.273, 0.797, 0.797, 1.273, 0.431, 1.041
1.273, 0.431, 1.041, 1.273, 0.0, 1.127,
1.273, 0.0, 1.127, 1.273, -0.431, 1.041
1.273, -0.431, 1.041, 1.273, -0.797, 0.797
1.273, -0.797, 0.797, 1.273, -1.041, 0.431
1.273, -1.041, 0.431, 1.273, -1.127, 0.0
1.273, -1.127, 0.0, 1.273, -1.041, -0.431
1.273, -1.041, -0.431, 1.273, -0.797, -0.797
1.273, -0.797, -0.797, 1.273, -0.431, -1.041
1.273, -0.431, -1.041, 1.273, 0.0, -1.127
1.273, 0.0, -1.127, 1.273, 0.431, -1.041
1.273, 0.431, -1.041, 1.273, 0.797, -0.797
1.273, 0.797, -0.797, 1.273, 1.041, -0.431
1.273, 1.041, -0.431, 1.273, 1.127, 0.0
16
0.0, 1.273, 1.127, 0.431, 1.273, 1.041
0.431, 1.273, 1.041, 0.797, 1,273, 0.797
0.797, 1,273, 0.797, 1.041, 1.273, 0.431
1.041, 1.273, 0.431, 1.127, 1.273, 0.0
1.127, 1.273, 0.0, 1.041, 1.273, -0.431
1.041, 1.273, -0.431, 0.797, 1.273, -0.797
0.797, 1.273, -0.797, 0.431, 1.273, -1.041
0.431, 1.273, -1.041, 0.0, 1.273, -1.127
0.0, 1.273, -1.127, -0.431, 1.273, -1.041
-0.431, 1.273, -1.041, -0.797, 1.273, -0.797
-0.797, 1.273, -0.797, -1.041, 1.273, -0.431
-1.041, 1.273, -0.431, -1.127, 1.273, 0.0
-1.127, 1.273, 0.0, -1.041, 1.273, 0.431
-1.041, 1.273, 0.431, -0.797, 1.273, 0.797
-0.797, 1.273, 0.797, -0.431, 1.273, 1.041
-0.431, 1.273, 1.041, 0.0, 1.273, 1.127
16
0.0, -1.273, 1.127, -0.431, -1.273, 1.041
-0.431, -1.273, 1.041, -0.797, -1.273, 0.797,
-0.797, -1.273, 0.797, -1.041, -1.273, 0.431
-1.041, -1.273, 0.431, -1.127, -1.273, 0.0
-1.127, -1.273, 0.0, -1.041, -1.273, -0.431
-1.041, -1.273, -0.431, -0.797, -1.273, -0.797
-0.797, -1.273, -0.797, -0.431, -1.273, -1.041
-0.431, -1.273, -1.041, 0.0, -1.273, -1.127
0.0, -1.273, -1.127, 0.431, -1.273, -1.041
0.431, -1.273, -1.041, 0.797, -1.273, -0.797
0.797, -1.273, -0.797, 1.041, -1.273, -0.431
1.041, -1.273, -0.431, 1.127, -1.273, 0.0
1.127, -1.273, 0.0, 1.041, -1.273, 0.431
1.041, -1.273, 0.431, 0.797, -1.273, 0.797
0.797, -1.273, 0.797, 0.431, -1.273, 1.041
0.431, -1.273, 1.041, 0.0, -1.273, 1.127``````
100% checked. Found one error in third(?) ring. Lesson learned: I won't post until I do a 100% check.

Also, if some wish to join the fun, I will happy to post my models.
"Aqaba! By Land!" T. E. Lawrence

R. Peters

happyjack27
Posts: 1438
Joined: Wed Jul 14, 2010 5:27 pm
still something wrong with your y = +/- 1.273 coils. i ran facing pairs of coils, and all the other pairs look correct. i tried inverting one or the other coil and it still didn't work right.

happyjack27
Posts: 1438
Joined: Wed Jul 14, 2010 5:27 pm
maybe start w/my cuboctahedron (every 4 lines is a coil), which i know is correct 'cause i've tested, and push the coils out from the center to get the desired coil separation, then slice each segment into 2, then 4.

i'm going to see if i can get particle sources and sinks working over the weekend. then we can actually see electrons with the right KE per radial displacement.

rjaypeters
Posts: 869
Joined: Fri Aug 20, 2010 2:04 pm
Location: Summerville SC, USA
Frustrating, I'll send some tests soon.
"Aqaba! By Land!" T. E. Lawrence

R. Peters

happyjack27
Posts: 1438
Joined: Wed Jul 14, 2010 5:27 pm
i just made a 16-gon approximation of a wb-6 with unit-radius coils and no coil separation. just tell me what the distance of a coil should be from the center to get the desired separation and i'll plug it in using "find and replace".

Code: Select all

``````//============ top, back, and right coils
-1.0000000, 1.000000a, 0.0000000,-0.9238795, 1.000000a, 0.3826834,-1
-0.9238795, 1.000000a, 0.3826834,-0.7071068, 1.000000a, 0.7071068,-1
-0.7071068, 1.000000a, 0.7071068,-0.3826834, 1.000000a, 0.9238795,-1
-0.3826834, 1.000000a, 0.9238795, 0.0000000, 1.000000a, 1.0000000,-1
0.0000000, 1.000000a, 1.0000000, 0.3826834, 1.000000a, 0.9238795,-1
0.3826834, 1.000000a, 0.9238795, 0.7071068, 1.000000a, 0.7071068,-1
0.7071068, 1.000000a, 0.7071068, 0.9238795, 1.000000a, 0.3826834,-1
0.9238795, 1.000000a, 0.3826834, 1.0000000, 1.000000a, 0.0000000,-1
1.0000000, 1.000000a, 0.0000000, 0.9238795, 1.000000a,-0.3826834,-1
0.9238795, 1.000000a,-0.3826834, 0.7071068, 1.000000a,-0.7071068,-1
0.7071068, 1.000000a,-0.7071068, 0.3826834, 1.000000a,-0.9238795,-1
0.3826834, 1.000000a,-0.9238795, 0.0000000, 1.000000a,-1.0000000,-1
0.0000000, 1.000000a,-1.0000000,-0.3826834, 1.000000a,-0.9238795,-1
-0.3826834, 1.000000a,-0.9238795,-0.7071068, 1.000000a,-0.7071068,-1
-0.7071068, 1.000000a,-0.7071068,-0.9238795, 1.000000a,-0.3826834,-1
-0.9238795, 1.000000a,-0.3826834,-1.0000000, 1.000000a, 0.0000000,-1

-1.0000000, 0.0000000, 1.000000a,-0.9238795, 0.3826834, 1.000000a
-0.9238795, 0.3826834, 1.000000a,-0.7071068, 0.7071068, 1.000000a
-0.7071068, 0.7071068, 1.000000a,-0.3826834, 0.9238795, 1.000000a
-0.3826834, 0.9238795, 1.000000a, 0.0000000, 1.0000000, 1.000000a
0.0000000, 1.0000000, 1.000000a, 0.3826834, 0.9238795, 1.000000a
0.3826834, 0.9238795, 1.000000a, 0.7071068, 0.7071068, 1.000000a
0.7071068, 0.7071068, 1.000000a, 0.9238795, 0.3826834, 1.000000a
0.9238795, 0.3826834, 1.000000a, 1.0000000, 0.0000000, 1.000000a
1.0000000, 0.0000000, 1.000000a, 0.9238795,-0.3826834, 1.000000a
0.9238795,-0.3826834, 1.000000a, 0.7071068,-0.7071068, 1.000000a
0.7071068,-0.7071068, 1.000000a, 0.3826834,-0.9238795, 1.000000a
0.3826834,-0.9238795, 1.000000a, 0.0000000,-1.0000000, 1.000000a
0.0000000,-1.0000000, 1.000000a,-0.3826834,-0.9238795, 1.000000a
-0.3826834,-0.9238795, 1.000000a,-0.7071068,-0.7071068, 1.000000a
-0.7071068,-0.7071068, 1.000000a,-0.9238795,-0.3826834, 1.000000a
-0.9238795,-0.3826834, 1.000000a,-1.0000000, 0.0000000, 1.000000a

1.000000a,-1.0000000, 0.0000000, 1.000000a,-0.9238795, 0.3826834
1.000000a,-0.9238795, 0.3826834, 1.000000a,-0.7071068, 0.7071068
1.000000a,-0.7071068, 0.7071068, 1.000000a,-0.3826834, 0.9238795
1.000000a,-0.3826834, 0.9238795, 1.000000a, 0.0000000, 1.0000000
1.000000a, 0.0000000, 1.0000000, 1.000000a, 0.3826834, 0.9238795
1.000000a, 0.3826834, 0.9238795, 1.000000a, 0.7071068, 0.7071068
1.000000a, 0.7071068, 0.7071068, 1.000000a, 0.9238795, 0.3826834
1.000000a, 0.9238795, 0.3826834, 1.000000a, 1.0000000, 0.0000000
1.000000a, 1.0000000, 0.0000000, 1.000000a, 0.9238795,-0.3826834
1.000000a, 0.9238795,-0.3826834, 1.000000a, 0.7071068,-0.7071068
1.000000a, 0.7071068,-0.7071068, 1.000000a, 0.3826834,-0.9238795
1.000000a, 0.3826834,-0.9238795, 1.000000a, 0.0000000,-1.0000000
1.000000a, 0.0000000,-1.0000000, 1.000000a,-0.3826834,-0.9238795
1.000000a,-0.3826834,-0.9238795, 1.000000a,-0.7071068,-0.7071068
1.000000a,-0.7071068,-0.7071068, 1.000000a,-0.9238795,-0.3826834
1.000000a,-0.9238795,-0.3826834, 1.000000a,-1.0000000, 0.0000000

//============ bottom, front, and left coils
-1.0000000,-1.000000a, 0.0000000,-0.9238795,-1.000000a, 0.3826834
-0.9238795,-1.000000a, 0.3826834,-0.7071068,-1.000000a, 0.7071068
-0.7071068,-1.000000a, 0.7071068,-0.3826834,-1.000000a, 0.9238795
-0.3826834,-1.000000a, 0.9238795, 0.0000000,-1.000000a, 1.0000000
0.0000000,-1.000000a, 1.0000000, 0.3826834,-1.000000a, 0.9238795
0.3826834,-1.000000a, 0.9238795, 0.7071068,-1.000000a, 0.7071068
0.7071068,-1.000000a, 0.7071068, 0.9238795,-1.000000a, 0.3826834
0.9238795,-1.000000a, 0.3826834, 1.0000000,-1.000000a, 0.0000000
1.0000000,-1.000000a, 0.0000000, 0.9238795,-1.000000a,-0.3826834
0.9238795,-1.000000a,-0.3826834, 0.7071068,-1.000000a,-0.7071068
0.7071068,-1.000000a,-0.7071068, 0.3826834,-1.000000a,-0.9238795
0.3826834,-1.000000a,-0.9238795, 0.0000000,-1.000000a,-1.0000000
0.0000000,-1.000000a,-1.0000000,-0.3826834,-1.000000a,-0.9238795
-0.3826834,-1.000000a,-0.9238795,-0.7071068,-1.000000a,-0.7071068
-0.7071068,-1.000000a,-0.7071068,-0.9238795,-1.000000a,-0.3826834
-0.9238795,-1.000000a,-0.3826834,-1.0000000,-1.000000a, 0.0000000

-1.0000000, 0.0000000,-1.000000a,-0.9238795, 0.3826834,-1.000000a,-1
-0.9238795, 0.3826834,-1.000000a,-0.7071068, 0.7071068,-1.000000a,-1
-0.7071068, 0.7071068,-1.000000a,-0.3826834, 0.9238795,-1.000000a,-1
-0.3826834, 0.9238795,-1.000000a, 0.0000000, 1.0000000,-1.000000a,-1
0.0000000, 1.0000000,-1.000000a, 0.3826834, 0.9238795,-1.000000a,-1
0.3826834, 0.9238795,-1.000000a, 0.7071068, 0.7071068,-1.000000a,-1
0.7071068, 0.7071068,-1.000000a, 0.9238795, 0.3826834,-1.000000a,-1
0.9238795, 0.3826834,-1.000000a, 1.0000000, 0.0000000,-1.000000a,-1
1.0000000, 0.0000000,-1.000000a, 0.9238795,-0.3826834,-1.000000a,-1
0.9238795,-0.3826834,-1.000000a, 0.7071068,-0.7071068,-1.000000a,-1
0.7071068,-0.7071068,-1.000000a, 0.3826834,-0.9238795,-1.000000a,-1
0.3826834,-0.9238795,-1.000000a, 0.0000000,-1.0000000,-1.000000a,-1
0.0000000,-1.0000000,-1.000000a,-0.3826834,-0.9238795,-1.000000a,-1
-0.3826834,-0.9238795,-1.000000a,-0.7071068,-0.7071068,-1.000000a,-1
-0.7071068,-0.7071068,-1.000000a,-0.9238795,-0.3826834,-1.000000a,-1
-0.9238795,-0.3826834,-1.000000a,-1.0000000, 0.0000000,-1.000000a,-1

-1.000000a,-1.0000000, 0.0000000,-1.000000a,-0.9238795, 0.3826834,-1
-1.000000a,-0.9238795, 0.3826834,-1.000000a,-0.7071068, 0.7071068,-1
-1.000000a,-0.7071068, 0.7071068,-1.000000a,-0.3826834, 0.9238795,-1
-1.000000a,-0.3826834, 0.9238795,-1.000000a, 0.0000000, 1.0000000,-1
-1.000000a, 0.0000000, 1.0000000,-1.000000a, 0.3826834, 0.9238795,-1
-1.000000a, 0.3826834, 0.9238795,-1.000000a, 0.7071068, 0.7071068,-1
-1.000000a, 0.7071068, 0.7071068,-1.000000a, 0.9238795, 0.3826834,-1
-1.000000a, 0.9238795, 0.3826834,-1.000000a, 1.0000000, 0.0000000,-1
-1.000000a, 1.0000000, 0.0000000,-1.000000a, 0.9238795,-0.3826834,-1
-1.000000a, 0.9238795,-0.3826834,-1.000000a, 0.7071068,-0.7071068,-1
-1.000000a, 0.7071068,-0.7071068,-1.000000a, 0.3826834,-0.9238795,-1
-1.000000a, 0.3826834,-0.9238795,-1.000000a, 0.0000000,-1.0000000,-1
-1.000000a, 0.0000000,-1.0000000,-1.000000a,-0.3826834,-0.9238795,-1
-1.000000a,-0.3826834,-0.9238795,-1.000000a,-0.7071068,-0.7071068,-1
-1.000000a,-0.7071068,-0.7071068,-1.000000a,-0.9238795,-0.3826834,-1
-1.000000a,-0.9238795,-0.3826834,-1.000000a,-1.0000000, 0.0000000,-1
``````
and just for the hell of it,

Code: Select all

``````//vertices of a quadrant of a unit 32-gon:
0.0000000, 1.0000000, 0.0000000
0.1950903, 0.9807853, 0.0000000
0.3826834, 0.9238795, 0.0000000
0.5555702, 0.8314696, 0.0000000
0.7071068, 0.7071068, 0.0000000
0.8314696, 0.5555702, 0.0000000
0.9238795, 0.3826834, 0.0000000
0.9807853, 0.1950903, 0.0000000
1.0000000, 0.0000000, 0.0000000
``````
which, judging from this picture, is all the precision we could ever possibly need.

happyjack27
Posts: 1438
Joined: Wed Jul 14, 2010 5:27 pm
k, i made a 32-gon unit wb-6.

Code: Select all

`````` 0.000000a, 0.0000000, 1.0000000, 0.000000a, 0.1950903, 0.9807853,-1
0.000000a, 0.1950903, 0.9807853, 0.000000a, 0.3826834, 0.9238795,-1
0.000000a, 0.3826834, 0.9238795, 0.000000a, 0.5555702, 0.8314696,-1
...
``````
what should 0.000000a = ? (the major radius of the containment sphere) what should the spatial scale be? other parameters (electron gun distance from center, initial electron velocity, etc., coil current density, coil charge density, etc.)

rjaypeters
Posts: 869
Joined: Fri Aug 20, 2010 2:04 pm
Location: Summerville SC, USA
+- Y coils at cardinal points only.

Code: Select all

``````4
0.0, -1.273, 1.127, -1.127, -1.273, 0.0
-1.127, -1.273, 0.0, 0.0, -1.273, -1.127
0.0, -1.273, -1.127, 1.127, -1.273, 0.0
1.127, -1.273, 0.0, 0.0, -1.273, 1.127
4
0.0, 1.273, 1.127, 1.127, 1.273, 0.0
1.127, 1.273, 0.0, 0.0, 1.273, -1.127
0.0, 1.273, -1.127, -1.127, 1.273, 0.0
-1.127, 1.273, 0.0, 0.0, 1.273, 1.127``````
Cardinal points at a different distance from origin.

-y coil:
0.0, -2, 1.127 top of compass
-1.127, -2, 0.0 west
0.0, -2, -1.127 south
1.127, -2, 0.0 east

+y coil:
0.0, 2, 1.127 north
1.127, 2, 0.0 west
0.0, 2, -1.127 south
-1.127, 2, 0.0 east

Code: Select all

``````4
0.0, -2, 1.127, -1.127, -2, 0.0
-1.127, -2, 0.0, 0.0, -2, -1.127
0.0, -2, -1.127, 1.127, -2, 0.0
1.127, -2, 0.0, 0.0, -2, 1.127
4
0.0, 2, 1.127, 1.127, 2, 0.0
1.127, 2, 0.0, 0.0, 2, -1.127
0.0, 2, -1.127, -1.127, 2, 0.0
-1.127, 2, 0.0, 0.0, 2, 1.127``````
Triangles:

Code: Select all

``````3
0.0, -1.273, 1.127, -1.041, -1.273, -0.431
-1.041, -1.273, -0.431, 1.041, -1.273, -0.431
1.041, -1.273, -0.431, 0.0, -1.273, 1.127
3
0.0, 1.273, 1.127, 1.041, 1.273, -0.431
1.041, 1.273, -0.431, -1.041, 1.273, -0.431
-1.041, 1.273, -0.431, 0.0, 1.273, 1.127``````
I am using a sphere 1.5m in radius and the coil centers are 1.273m from the origin (the number that doesn't change in the examples above).
"Aqaba! By Land!" T. E. Lawrence

R. Peters

happyjack27
Posts: 1438
Joined: Wed Jul 14, 2010 5:27 pm
and the coil radius is 1.127 this for a unit coil radius, their centers should be 1.273/1.127 from the origin?

happyjack27
Posts: 1438
Joined: Wed Jul 14, 2010 5:27 pm
wb-6 style, distance from coil midplane to origin 1.1 times coil radius. (giving a decent amount of coil separation) electrons start as a hollow sphere. you can make out some weak line cusps due to the coil separation and the fact that the coils are circular (actually 32-gons).

rjaypeters
Posts: 869
Joined: Fri Aug 20, 2010 2:04 pm
Location: Summerville SC, USA
happyjack27 wrote:and the coil radius is 1.127 this for a unit coil radius, their centers should be 1.273/1.127 from the origin?
Yes, but you knew that already.
"Aqaba! By Land!" T. E. Lawrence

R. Peters

Posts: 6204
Joined: Thu Sep 17, 2009 11:18 pm
Location: North East Coast
happyjack27 wrote:wb-6 style, distance from coil midplane to origin 1.1 times coil radius. (giving a decent amount of coil separation) electrons start as a hollow sphere. you can make out some weak line cusps due to the coil separation and the fact that the coils are circular (actually 32-gons).