Los Alamos National Laboratory Poisson Superfish Program Automesh written by James H. Billen and Lloyd M. Young The original Poisson Superfish codes were developed by Ron F. Holsinger in collaboration with Klaus Halbach. These programs are provided as a service to the accelerator community by the Los Alamos Accelerator Code Group (LAACG). (c) Copyright 1985-2005, by the Regents of the University of California. This software was produced under U. S. Government contract W-7405-ENG-36 by Los Alamos National Laboratory, which is operated by the University of California for the U. S. Department of Energy. Neither the Government nor the University makes any warranty, express or implied, or assumes any liability or responsibility for its use, or represents that use of this software would not infringe privately owned rights. Unpublished - rights reserved under Copyright Laws of the United States. Program Automesh 7.17 released 1-13-2006 Program file: C:\LANL\AUTOFISH.EXE Problem file: J:\.WINE\PILLBOX-SLAC\PILLBOX-SLAC-805.AF 4-28-2008 11:58:54 SF.INI file: J:\\SF.INI 3-31-2008 17:23:32 Problem description: 0.805-GHz TM010 Short Pillbox Cavity In this problem, Kmax < Lmax Coordinates and lengths have dimensions of centimeters. Region data: IREG = 1 MAT = 1 CUR = 0.0 DEN = 0.0 ITRI = 0 [equal weight triangles] IBOUND = 1 [Neumann boundary] IPDIAG = 0 [no extra Automesh diagnostics] DX = 0.2 XMIN = 0.0 XMAX = 5 DY = 0.1732051 [=DX*sin(60)] YMIN = 0.0 YMAX = 14.25 DX1 = 0.2 KMAX = 26 DY1 = 0.1737805 LMAX = 83 ITOT = 2380 Memory used for the solution file: 85.680 K Memory used for Automesh setup data: 110.800 K Region 1 boundary points from the PO namelist: Point NT X0 Y0 X Y R THETA A B 1 1 0.000 0.000 2 1 0.000 14.25 3 1 5.000 14.25 4 1 5.000 0.000 5 1 0.000 0.000 Number of user-supplied fixed points = 5 Fixed points added at line regions = 0 Fixed points added on overlapping regions = 0 Total number of fixed points = 5 Logical path finding for region 1 (includes extra boundary points at intersections of line regions plus points from other user-supplied segments) Point NT X Y R A B Fwd Bkw 2 1 0.000 14.25 83 83 3 1 5.000 14.25 26 26 4 1 5.000 0.000 83 83 5 1 0.000 0.000 26 26 Region 1 mesh points K L X Y 1 1 0.00000000 0.00000000 1 2 0.00000000 0.173780488 1 3 0.00000000 0.347560976 1 4 0.00000000 0.521341463 1 5 0.00000000 0.695121951 1 6 0.00000000 0.868902439 1 7 0.00000000 1.04268293 1 8 0.00000000 1.21646341 1 9 0.00000000 1.39024390 1 10 0.00000000 1.56402439 1 11 0.00000000 1.73780488 1 12 0.00000000 1.91158537 1 13 0.00000000 2.08536585 1 14 0.00000000 2.25914634 1 15 0.00000000 2.43292683 1 16 0.00000000 2.60670732 1 17 0.00000000 2.78048780 1 18 0.00000000 2.95426829 1 19 0.00000000 3.12804878 1 20 0.00000000 3.30182927 1 21 0.00000000 3.47560976 1 22 0.00000000 3.64939024 1 23 0.00000000 3.82317073 1 24 0.00000000 3.99695122 1 25 0.00000000 4.17073171 1 26 0.00000000 4.34451220 1 27 0.00000000 4.51829268 1 28 0.00000000 4.69207317 1 29 0.00000000 4.86585366 1 30 0.00000000 5.03963415 1 31 0.00000000 5.21341463 1 32 0.00000000 5.38719512 1 33 0.00000000 5.56097561 1 34 0.00000000 5.73475610 1 35 0.00000000 5.90853659 1 36 0.00000000 6.08231707 1 37 0.00000000 6.25609756 1 38 0.00000000 6.42987805 1 39 0.00000000 6.60365854 1 40 0.00000000 6.77743902 1 41 0.00000000 6.95121951 1 42 0.00000000 7.12500000 1 43 0.00000000 7.29878049 1 44 0.00000000 7.47256098 1 45 0.00000000 7.64634146 1 46 0.00000000 7.82012195 1 47 0.00000000 7.99390244 1 48 0.00000000 8.16768293 1 49 0.00000000 8.34146341 1 50 0.00000000 8.51524390 1 51 0.00000000 8.68902439 1 52 0.00000000 8.86280488 1 53 0.00000000 9.03658537 1 54 0.00000000 9.21036585 1 55 0.00000000 9.38414634 1 56 0.00000000 9.55792683 1 57 0.00000000 9.73170732 1 58 0.00000000 9.90548780 1 59 0.00000000 10.0792683 1 60 0.00000000 10.2530488 1 61 0.00000000 10.4268293 1 62 0.00000000 10.6006098 1 63 0.00000000 10.7743902 1 64 0.00000000 10.9481707 1 65 0.00000000 11.1219512 1 66 0.00000000 11.2957317 1 67 0.00000000 11.4695122 1 68 0.00000000 11.6432927 1 69 0.00000000 11.8170732 1 70 0.00000000 11.9908537 1 71 0.00000000 12.1646341 1 72 0.00000000 12.3384146 1 73 0.00000000 12.5121951 1 74 0.00000000 12.6859756 1 75 0.00000000 12.8597561 1 76 0.00000000 13.0335366 1 77 0.00000000 13.2073171 1 78 0.00000000 13.3810976 1 79 0.00000000 13.5548780 1 80 0.00000000 13.7286585 1 81 0.00000000 13.9024390 1 82 0.00000000 14.0762195 1 83 0.00000000 14.2500000 2 83 0.200000000 14.2500000 3 83 0.400000000 14.2500000 4 83 0.600000000 14.2500000 5 83 0.800000000 14.2500000 6 83 1.00000000 14.2500000 7 83 1.20000000 14.2500000 8 83 1.40000000 14.2500000 9 83 1.60000000 14.2500000 10 83 1.80000000 14.2500000 11 83 2.00000000 14.2500000 12 83 2.20000000 14.2500000 13 83 2.40000000 14.2500000 14 83 2.60000000 14.2500000 15 83 2.80000000 14.2500000 16 83 3.00000000 14.2500000 17 83 3.20000000 14.2500000 18 83 3.40000000 14.2500000 19 83 3.60000000 14.2500000 20 83 3.80000000 14.2500000 21 83 4.00000000 14.2500000 22 83 4.20000000 14.2500000 23 83 4.40000000 14.2500000 24 83 4.60000000 14.2500000 25 83 4.80000000 14.2500000 26 83 5.00000000 14.2500000 26 82 5.00000000 14.0762195 26 81 5.00000000 13.9024390 26 80 5.00000000 13.7286585 26 79 5.00000000 13.5548780 26 78 5.00000000 13.3810976 26 77 5.00000000 13.2073171 26 76 5.00000000 13.0335366 26 75 5.00000000 12.8597561 26 74 5.00000000 12.6859756 26 73 5.00000000 12.5121951 26 72 5.00000000 12.3384146 26 71 5.00000000 12.1646341 26 70 5.00000000 11.9908537 26 69 5.00000000 11.8170732 26 68 5.00000000 11.6432927 26 67 5.00000000 11.4695122 26 66 5.00000000 11.2957317 26 65 5.00000000 11.1219512 26 64 5.00000000 10.9481707 26 63 5.00000000 10.7743902 26 62 5.00000000 10.6006098 26 61 5.00000000 10.4268293 26 60 5.00000000 10.2530488 26 59 5.00000000 10.0792683 26 58 5.00000000 9.90548780 26 57 5.00000000 9.73170732 26 56 5.00000000 9.55792683 26 55 5.00000000 9.38414634 26 54 5.00000000 9.21036585 26 53 5.00000000 9.03658537 26 52 5.00000000 8.86280488 26 51 5.00000000 8.68902439 26 50 5.00000000 8.51524390 26 49 5.00000000 8.34146341 26 48 5.00000000 8.16768293 26 47 5.00000000 7.99390244 26 46 5.00000000 7.82012195 26 45 5.00000000 7.64634146 26 44 5.00000000 7.47256098 26 43 5.00000000 7.29878049 26 42 5.00000000 7.12500000 26 41 5.00000000 6.95121951 26 40 5.00000000 6.77743902 26 39 5.00000000 6.60365854 26 38 5.00000000 6.42987805 26 37 5.00000000 6.25609756 26 36 5.00000000 6.08231707 26 35 5.00000000 5.90853659 26 34 5.00000000 5.73475610 26 33 5.00000000 5.56097561 26 32 5.00000000 5.38719512 26 31 5.00000000 5.21341463 26 30 5.00000000 5.03963415 26 29 5.00000000 4.86585366 26 28 5.00000000 4.69207317 26 27 5.00000000 4.51829268 26 26 5.00000000 4.34451220 26 25 5.00000000 4.17073171 26 24 5.00000000 3.99695122 26 23 5.00000000 3.82317073 26 22 5.00000000 3.64939024 26 21 5.00000000 3.47560976 26 20 5.00000000 3.30182927 26 19 5.00000000 3.12804878 26 18 5.00000000 2.95426829 26 17 5.00000000 2.78048780 26 16 5.00000000 2.60670732 26 15 5.00000000 2.43292683 26 14 5.00000000 2.25914634 26 13 5.00000000 2.08536585 26 12 5.00000000 1.91158537 26 11 5.00000000 1.73780488 26 10 5.00000000 1.56402439 26 9 5.00000000 1.39024390 26 8 5.00000000 1.21646341 26 7 5.00000000 1.04268293 26 6 5.00000000 0.868902439 26 5 5.00000000 0.695121951 26 4 5.00000000 0.521341463 26 3 5.00000000 0.347560976 26 2 5.00000000 0.173780488 26 1 5.00000000 0.00000000 25 1 4.80000000 0.00000000 24 1 4.60000000 0.00000000 23 1 4.40000000 0.00000000 22 1 4.20000000 0.00000000 21 1 4.00000000 0.00000000 20 1 3.80000000 0.00000000 19 1 3.60000000 0.00000000 18 1 3.40000000 0.00000000 17 1 3.20000000 0.00000000 16 1 3.00000000 0.00000000 15 1 2.80000000 0.00000000 14 1 2.60000000 0.00000000 13 1 2.40000000 0.00000000 12 1 2.20000000 0.00000000 11 1 2.00000000 0.00000000 10 1 1.80000000 0.00000000 9 1 1.60000000 0.00000000 8 1 1.40000000 0.00000000 7 1 1.20000000 0.00000000 6 1 1.00000000 0.00000000 5 1 0.800000000 0.00000000 4 1 0.600000000 0.00000000 3 1 0.400000000 0.00000000 2 1 0.200000000 0.00000000 1 1 0.00000000 0.00000000 Region 1 done, 5 fixed points, 215 total boundary points. Logical boundary segment end points for region 1: Segment Starting Point Ending Point K L K L 1 1 1 1 83 2 1 83 26 83 3 26 83 26 1 4 26 1 1 1 Drive point at X = 1, K = 6, Y = 14.25, L = 83. Region 2 mesh points K L X Y 6 83 1.00000000 14.2500000 Processing boundary data... Memory used for additional setup arrays: 72.784 K Region 1 Material = 1 Region boundary indicator = 1 Equal-weight triangles Region 1 is closed. Region 2 Material = 1 Region boundary indicator = -1 Driving field = 1 Equal-weight triangles Region 2 contains a single point. Logical boundary segment end points for all regions: Segment Starting Point Ending Point Index K L Index K L 1 1 1 1 83 1 83 2 83 1 83 108 26 83 3 108 26 83 190 26 1 4 190 26 1 215 1 1 Relaxation parameters, 1944 unknown points. Optimizing mesh by successive over relaxation... -------X coordinate-------- -------Y coordinate-------- Cycle Residual ETA RHO Residual ETA RHO 1 2.992E-02 1.0000 1.6000 1.583E-09 1.0000 1.6000 2 4.156E-02 0.6765 1.6000 2.219E-09 0.6943 1.6000 3 2.830E-02 0.6817 1.6000 1.545E-09 0.6964 1.6000 4 1.922E-02 0.6797 1.6000 1.081E-09 0.6994 1.6000 5 1.311E-02 0.6825 1.6000 7.583E-10 0.7016 1.6000 6 8.919E-03 0.6806 1.6000 5.347E-10 0.7052 1.6000 7 6.100E-03 0.6842 1.6000 3.780E-10 0.7069 1.6000 8 4.161E-03 0.6823 1.6000 2.685E-10 0.7104 1.6000 9 2.864E-03 0.6884 1.6000 1.916E-10 0.7133 1.6000 10 1.967E-03 0.6871 1.6000 1.376E-10 0.7184 1.6000 11 1.373E-03 0.6981 1.6000 9.931E-11 0.7216 1.6000 12 9.608E-04 0.7000 1.6000 7.216E-11 0.7266 1.6000 13 6.919E-04 0.7203 1.6000 5.277E-11 0.7313 1.6000 14 5.057E-04 0.7309 1.6000 3.900E-11 0.7390 1.6000 15 3.858E-04 0.7631 1.6000 2.909E-11 0.7459 1.6000 16 3.026E-04 0.7844 1.6000 2.197E-11 0.7553 1.6000 17 2.480E-04 0.8196 1.6000 1.684E-11 0.7664 1.6000 18 2.086E-04 0.8411 1.6000 1.314E-11 0.7805 1.6000 19 1.804E-04 0.8652 1.6000 1.045E-11 0.7953 1.6000 20 1.584E-04 0.8781 1.6000 8.493E-12 0.8127 1.6000 21 1.411E-04 0.8905 1.6000 7.053E-12 0.8304 1.6000 22 1.265E-04 0.8967 1.6000 5.980E-12 0.8479 1.6000 23 1.142E-04 0.9030 1.6000 5.165E-12 0.8637 1.6000 24 1.035E-04 0.9064 1.6000 4.531E-12 0.8773 1.6000 25 9.423E-05 0.9102 1.6000 4.025E-12 0.8884 1.6000 26 1.344E-04 1.4266 1.7468 4.892E-12 1.2154 1.7236 27 1.191E-04 0.8862 1.7468 4.400E-12 0.8992 1.7236 28 1.056E-04 0.8866 1.7468 3.963E-12 0.9008 1.7236 29 9.388E-05 0.8890 1.7468 3.577E-12 0.9025 1.7236 30 8.370E-05 0.8915 1.7468 3.234E-12 0.9041 1.7236 31 7.488E-05 0.8946 1.7468 2.928E-12 0.9055 1.7236 32 6.721E-05 0.8976 1.7468 2.655E-12 0.9067 1.7236 33 6.053E-05 0.9006 1.7468 2.410E-12 0.9079 1.7236 34 5.468E-05 0.9034 1.7468 2.191E-12 0.9089 1.7236 35 4.953E-05 0.9060 1.7468 1.993E-12 0.9099 1.7236 36 4.499E-05 0.9083 1.7468 1.816E-12 0.9108 1.7236 37 4.096E-05 0.9104 1.7468 1.655E-12 0.9116 1.7236 38 3.736E-05 0.9122 1.7468 1.510E-12 0.9125 1.7236 39 3.414E-05 0.9138 1.7468 1.379E-12 0.9132 1.7236 40 3.124E-05 0.9152 1.7468 1.261E-12 0.9141 1.7236 41 2.864E-05 0.9165 1.7468 1.153E-12 0.9147 1.7236 42 2.628E-05 0.9176 1.7468 1.056E-12 0.9154 1.7236 43 2.414E-05 0.9186 1.7468 9.672E-13 0.9162 1.7236 44 2.219E-05 0.9195 1.7468 8.867E-13 0.9168 1.7236 45 2.042E-05 0.9203 1.7468 8.137E-13 0.9177 1.7236 46 1.881E-05 0.9210 1.7468 7.470E-13 0.9180 1.7236 47 1.734E-05 0.9217 1.7468 6.866E-13 0.9192 1.7236 48 1.599E-05 0.9223 1.7468 6.314E-13 0.9196 1.7236 49 1.476E-05 0.9229 1.7468 5.810E-13 0.9202 1.7236 50 1.363E-05 0.9235 1.7468 5.351E-13 0.9210 1.7236 51 1.679E-05 1.2318 1.8197 6.648E-13 1.2424 1.8075 52 1.529E-05 0.9105 1.8197 6.039E-13 0.9084 1.8075 53 1.392E-05 0.9104 1.8197 5.490E-13 0.9090 1.8075 54 1.267E-05 0.9107 1.8197 4.994E-13 0.9097 1.8075 55 1.154E-05 0.9108 1.8197 4.543E-13 0.9098 1.8075 56 1.052E-05 0.9110 1.8197 4.135E-13 0.9101 1.8075 57 9.580E-06 0.9110 1.8197 3.765E-13 0.9104 1.8075 Iteration converged. Generation completed. Problem description: 0.805-GHz TM010 Short Pillbox Cavity In this problem, Kmax < Lmax Problem file: J:\.WINE\PILLBOX-SLAC\PILLBOX-SLAC-805.AF 4-28-2008 11:58:54 Problem file length: 501 bytes Originating program: Autofish Problem type: Unknown RF Cavity Problem constants and variables. Letter A in the code column indicates a value supplied in the Automesh input file. Variable Code Value Description ALPHAT 3.930000000E-03 Temperature coefficient of resistance ASCALE 3767.30313 Scaling factor for H at drive point BETA 0.0 Particle velocity BETA1 0.100000000 Starting BETA in transit-time table BETA2 0.950000000 Ending BETA in transit-time table CCLDELK 1.00000000 Increment for coupling for table in SFO CCLMAXK 6.00000000 Highest coupling for table in SFO CCLMINK 1.00000000 Lowest coupling for table in SFO CLENGTH 0.0 Cavity length for normalization in SFO CLIGHT 2.997924580E+10 Exact speed of light in cm/sec CONV 1.00000000 Length conversion (number of units per cm) DBETA 5.000000000E-02 BETA increment in transit-time table DELFR 0.0 Frequency step size for a resonance search DIAGDLL 0 If 1, DLL writes diagnostics to DiagDLL.txt DPHI 180.000000 Phase length of the problem geometry DSLOPE 1.00000000 Slope of D(k^2) function DSTOLER 2.000000000E-02 Tolerance required on D(k^2) DX1 0.200000000 First X mesh interval (at XMIN) DXMIN 0.200000000 Minimum X mesh interval (found by Automesh) DYMIN 0.173780488 Minimum Y mesh interval (found by Automesh) ENORM 1000000.00 Field normalization for NORM=4 option EPS0 8.854187818E-12 Permittivity of free space EPSIK 1.000000000E-04 Frequency convergence parameter EPSO 1.000000000E-05 Convergence parameter in mesh optimization EZERO 1000000.00 E0 for normalization in SFO when NORM=0 EZEROT 1000000.00 E0*T for normalization in SFO when NORM=1 FMU0 1.256637061E-06 Permeability of free space FREQ A 805.000000 RF cavity resonant frequency FREQD 0.0 Design frequency for a cavity (MHz) HPHI 5000.00000 Normalization magnetic field for NORM=2 IBETA 0 If >0, SFO writes transit-time vs BETA ICCC 1 1 for real arrays, 2 for complex arrays ICCP 1 If 1, compute material power loss ICORNER1 0 First corner segment for computing average H ICORNER2 0 Last corner segment for computing average H ICYLIN 1 0 for X,Y problems, 1 for Z,R problems IMAX 28 KMAX+2 INFODATA 0 Number of tuning-code parameters IOBSEG -1 First segment of the CCL outer boundary IPIVOT 1 Pivoting in matrix inversion routines IRESID 0 If 1, calculate potential residuals IRMAX 25 Used in optimization of RHOXY IRTYPE 0 Surface resistance option indicator ISLOT 0 If 1, SFO computes coupling-slot power loss ITFILE 0 If 1, SFO writes transit-time plot file ITOT 2380 (KMAX+2)*(LMAX+2) KDRI 6 K coordinate of the drive point KMAX 26 Number of horizontal logical mesh points KMETHOD 0 Wavenumber computation method (1= use beta) KPROB A 1 Problem type indicator (Superfish) LAST35 0 Code for last program to update T35 file LCYCLE 57 Iteration number in mesh optimization LDRI 83 L coordinate of the drive point LINT 1 Logical-mesh coordinate for Ez integration LMAX 83 Number of vertical logical mesh points MAXCY -1 Maximum number of cycles (-1: use default) MAXPPR 216 Maximum points per region NAIR 2131 Number of air points NBND 26 Number of Dirichlet boundary points NBSLF 1 Left-side boundary condition NBSLO 0 Lower boundary condition NBSRT 1 Right-side boundary condition NBSUP 1 Upper boundary condition NDRI 2331 Drive point index = IRLAX(NPINP) NEGAT 0 Zero-area triangle indicator NFE 0 Number of iron points NHSTEM 1 Number of half stems NINTER 0 Number of interface points NMATR 0 Number of material records in T35 file NORM 0 Normalization method in SFO NPBOUND 215 Total number of boundary points in the mesh NPINP 2158 Total points in problem NPONTS 1944 Number of unknown relaxation points NREG 2 Number of regions NRMSEG 1 Normalization segment number for NORM=2 NSEG 4 Number of boundary segments NSPL 1 Number of special-potential points NSTEP 0 Number of steps for a resonance search OMEGAM 1.000000000E-03 Used in optimization of RHOXY PI 3.14159265 The number pi to machine precision PLCELL 360.000000 Phase length per cell for multicell problems RESIDR 1.000000000E-08 Residual resistance of a superconductor RFMU 1.00000000 Permeability for rf surface resistance RHO 1.724100000E-06 Material resistivity (Ohm-m) RHOR 1.724100000E-06 Reference resistivity (Ohm-cm) at TEMPR RHOXY 1.60000000 Over-relaxation factor in mesh optimization RMASS -2.00000000 Rest mass energy of particle in SFO RS 0.0 RF surface resistance (Ohms) RSTEM 1.00000000 Stem radius in cm SLOSS 3.000000000E-02 Coupling-slot power factor per % coupling TC 9.20000000 Critical temperature of a superconductor TEMPC 20.0000000 Normal conductor operating temperature TEMPK 2.00000000 Operating temperature of a superconductor TEMPR 20.0000000 Reference temperature for IRTYPE=3 TRIAVG 1.737804878E-02 Average area of all triangles TRIMAX 2.291078426E-02 Area of the largest positive-area triangle TRIMIN 1.044111782E-02 Area of the smallest positive-area triangle VOLUME 3189.69829 Cavity volume (cylindrical symmetry only) XDRI A 1.00000000 X coordinate of the drive point XMAXG 5.00000000 Upper X bound of the problem geometry XMING 0.0 Lower X bound of the problem geometry XNORM1 0.0 Starting X for NORM=4 integration path XNORM2 0.0 Ending X for NORM=4 integration path XYAREA 71.2500000 Total cross sectional area YDRI A 14.2500000 Y coordinate of the drive point YMAXG 14.2500000 Upper Y bound of the problem geometry YMING 0.0 Lower Y bound of the problem geometry YNORM1 0.0 Starting Y for NORM=4 integration path YNORM2 0.0 Ending Y for NORM=4 integration path ZCTR 0.0 Reference Z in transit-time integrals