Xmpl2d20, 20th 'example' data file for CPO2D

Field and potential penetration through a flat mesh of parallel round wires.

One of a set of files:

test2d18.dat: Field penetration through a flat mesh of flat strips.

xmpl2d19.dat: Field penetration through a cylindrical mesh of round wires.

xmpl2d20.dat: Field penetration through a flat mesh of round wires.

xmpl3d23.dat: Field and potential penetration through a flat mesh of crossed round wires.

xmpl3d29.dat: Field and potential penetration through a flat mesh of crossed flat wires.

For further details see the papers:

Short and long range penetration of fields and potentials through meshes, grids or gauzes, by F H Read, N J Bowring, P D Bullivant and R R A Ward Nucl. Instr. Meth. A427, 363-367 (1999).

Penetration of electrostatic fields and potentials through meshes, grids or gauzes, by F H Read, N J Bowring, P D Bullivant and R R A Ward, Rev. Sci. Instrum. 69, 2000-2006 (1998).

The following data were obtained when the memory and speed of PC's was much more limited than at present, so the available number of segments was small and the requested inaccuracies were fairly high to give a quick demonstration.

The mesh is composed of 20 parallel wires (in planar symmetry) of radius 0.0125, spaced by 1.0. The mesh is at potential 0 and is enclosed by a box that has end potentials of 0 at z=-10 and +10 at z=+10 (but see below). The sides of the box at z>0 and z>0 are given uniform potential gradients.

The following quantities are obtained:

(1) v1 = the average potential at z = -r, where r is the radius of the wires

(2) v2 = the average potential at z = +r

(3) vc = the potential at the centre of the mesh (that is, at x=0, z=0)

(4) ec = the field at the centre of the mesh.

All these values are extrapolated to an infinite number of segments (see below).

For further details see the paper by F H Read, N J Bowring, P D Bullivant and R R A Ward, submitted to Rev Sci Instrum, to appear (hopefully) in 1998.

The mesh is at zero potential and fields are created on both sides of it by means of continuous end plates at a different potential. An extrapolation technique is used (see below) to simulate a mesh that is infinite in extent.

To simulate end plates at infinity, vl is applied to the plate at z=-10, 10+vr is applied to the plate at z=10, and va (=average of vl and vr) is applied to the mesh surround at z=0.

The resulting values of the potential v1 averaged over z=-r, and v2 averaged over z=r, are noted and then vl and vr are adjusted iteratively until v1=vl and v2=vr. The effective fields are then 0 and 1V/mm, as required.

vc and ec are the potential and field at the origin.

Results:

(The initial values of vl and vr are taken from a previous, less accurate study)

r |
nseg |
vr |
vl |
va |
v1 |
v2 |
vc |
ec |

0.0125 |
336 |
0.415 |
0.403 |
0.409 |
0.415618 |
0.403654 |
0.515885 |
-0.499872 |

0.0125 |
672 |
0.415 |
0.403 |
0.409 |
0.415477 |
0.403514 |
0.515732 |
-0.499869 |

0.0125 |
996 |
0.415 |
0.403 |
0.409 |
0.415446 |
0.403482 |
0.515697 |
-0.499865 |

0.0125 |
inf |
0.415 |
0.403 |
0.409 |
0.40944 |
0.00598 |
0.51568 |
-0.49987 |

0.0250 |
336 |
0.317 |
0.294 |
0.306 |
0.317979 |
0.294992 |
0.406112 |
-0.498199 |

0.0250 |
672 |
0.317 |
0.294 |
0.306 |
0.317844 |
0.294861 |
0.405966 |
-0.498193 |

0.0250 |
inf |
0.317 |
0.294 |
0.306 |
0.31780 |
0.29482 |
0.40592 |
-0.49819 |

0.0250 |
inf |
0.317 |
0.294 |
0.306 |
0.30631 |
0.01149 |
0.40592 |
-0.49819 |

0.050 |
336 |
0.229 |
0.187 |
0.208 |
0.230638 |
0.188409 |
0.296401 |
-0.490302 |

0.050 |
672 |
0.229 |
0.187 |
0.208 |
0.230513 |
0.188296 |
0.296265 |
-0.490279 |

0.050 |
996 |
0.229 |
0.187 |
0.208 |
0.230486 |
0.188271 |
0.296237 |
-0.490272 |

0.050 |
inf |
0.229 |
0.187 |
0.208 |
0.23047 |
0.18825 |
0.29622 |
-0.49027 |

0.050 |
inf |
0.229 |
0.187 |
0.208 |
0.20936 |
0.02110 |
0.29622 |
-0.49027 |

0.100 |
336 |
0.088 |
0.158 |
0.123 |
0.159834 |
0.089805 |
0.186737 |
-0.457138 |

0.100 |
672 |
0.088 |
0.158 |
0.123 |
0.159704 |
0.089719 |
0.186616 |
-0.457062 |

0.100 |
996 |
0.088 |
0.158 |
0.123 |
0.159676 |
0.089699 |
0.186590 |
-0.457042 |

0.100 |
inf |
0.088 |
0.158 |
0.123 |
0.15965 |
0.08968 |
0.18657 |
-0.45703 |

0.100 |
inf |
0.088 |
0.158 |
0.123 |
0.12467 |
0.03499 |
0.18657 |
-0.45703 |

Summary of final results for potentials:

r |
vm |
vd |
vc |
ec |

0.0125 |
0.40944 |
0.00598 |
0.51568 |
-0.49987 |

0.0250 |
0.30631 |
0.01149 |
0.40592 |
-0.49819 |

0.050 |
0.20936 |
0.02110 |
0.29622 |
-0.49027 |

0.100 |
0.12467 |
0.03499 |
0.18657 |
-0.45703 |

Investigation of average and maximum values of field, and fourier components of potential:

e = exp(-2*pi*z)

r = 0.0125:

z |
av*e |
max*e |
f1*e |
f2*e**2 |
f3*e**3 |

0 |
1.269 |
10.61 |
0.1513 |
-0.0716 |
0.045 |

-0.2 |
1.01894 |
1.39279 |
0.1588 |
-0.0792 |
0.053 |

-0.4 |
0.99948 |
1.08368 |
0.1588 |
-0.0790 |
0.052 |

-0.6 |
0.99766 |
1.01371 |
0.1587 |
-0.0773 |
0.016 |

-0.8 |
0.99644 |
1.02022 |
0.1583 |
-0.0600 |
-1.248 |

-1 |
0.99549 |
1.1042 |
0.1574 |
0.1128 |
-45.72 |

average should be 2*pi*f1 at large z.

Taking f1*e=0.1588, 2*pi*f1*e=0.9978, which is consistent with the values found for av*e.

Take f2*e**2=-0.079, f3*e**3=0.053, r = 0.025:

z |
av*e |
max*e |
f1*e |
f2*e**2 |
f3*e**3 |

0.000E+00 |
9.926E-01 |
5.299E+00 |
1.437E-01 |
-6.402E-02 |
3.752E-02 |

-2.000E-01 |
2.876E-01 |
3.929E-01 |
1.576E-01 |
-7.783E-02 |
5.122E-02 |

-4.000E-01 |
8.031E-02 |
8.792E-02 |
1.576E-01 |
-7.761E-02 |
4.991E-02 |

-6.000E-01 |
2.282E-02 |
2.405E-02 |
1.574E-01 |
-7.522E-02 |
4.806E-04 |

-8.000E-01 |
6.498E-03 |
7.206E-03 |
1.569E-01 |
-4.975E-02 |
-1.837E+00 |

-1.000E+00 |
1.893E-03 |
2.500E-03 |
1.554E-01 |
2.225E-01 |
-6.965E+01 |

Take f1*e=0.1576, f2*e**2=-0.078, f3*e**3=0.051

r = 0.05:

z |
av*e |
max*e |
f1*e |
f2*e**2 |
f3*e**3 |

0 |
0.7851 |
2.868 |
0.12815 |
-0.0487 |
0.02287 |

-0.2 |
0.97467 |
1.31725 |
0.15235 |
-0.0724 |
0.04573 |

-0.4 |
0.95788 |
1.05565 |
0.15227 |
-0.0721 |
0.04408 |

-0.6 |
0.95689 |
1.04234 |
0.15205 |
-0.0690 |
-0.02027 |

-0.8 |
0.96611 |
1.18192 |
0.15138 |
-0.0343 |
-2.525 |

-1 |
1.10527 |
1.71413 |
0.14913 |
0.357 |
-101.675 |

Take f1*e=0.1523, f2*e**2=-0.072, f3*e**3=0.046

Summary of fourier coefficients:

r |
a10 |
a20 |
a30 |

0.0125 |
0.1588 |
-0.079 |
0.053 |

0.025 |
0.1576 |
-0.078 |
0.051 |

0.050 |
0.1523 |
-0.072 |
0.046 |

Compare with 3D results (xmpl3d23.dat):

r |
a10 |
a20 |
a30 |
a11 |

0.0125 |
0.070 |
-0.037 |
0.025 |
0.028 |

0.025 |
0.067 |
-0.035 |
0.024 |
0.033 |

0.050 |
0.061 |
-0.031 |
0.020 |
0.040 |

Appendix A

Program used to deduce v1, v2, vc and ec from the results in the output data.File temp20a.dat:

-see first part of program 16.

Appendix B

Program used to derive fourier components of the potential, and average and maximum values of the field, from the results in the output data file (temp20a.dat), when the above file has the lines:

-see second part of program 16.