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: C:\LANL\Examples\RadioFrequency\PillboxCavities\805ROMA.AF 11-17-2007 14:57:34 SF.INI file: C:\LANL\SF.INI 11-12-2007 14:55:22 Problem description: 0.805-GHz TM010 Modified Pillbox Cavity Pillbox cavity with a 1.63-cm-radius bore tube [Originally appeared in 1987 User's Guide 1.4] 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 = 7.55 DY = 0.1732051 [=DX*sin(60)] YMIN = 0.0 YMAX = 14.52 DX1 = 0.1986842 KMAX = 39 DY1 = 0.1728571 LMAX = 85 ITOT = 3567 Memory used for the solution file: 128.412 K Memory used for Automesh setup data: 135.710 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.52 3 1 2.680 14.52 4 2 2.680 12.48 2.040 0.000 2.040 2.040 5 1 4.720 3.270 6 2 6.350 3.270 1.630 270.0 1.630 1.630 7 1 7.550 1.640 8 1 7.550 0.000 9 1 0.000 0.000 Number of user-supplied fixed points = 9 Fixed points added at line regions = 0 Fixed points added on overlapping regions = 0 Total number of fixed points = 9 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.52 85 85 3 1 2.680 14.52 14 14 4 2 4.720 12.48 2.040 2.040 19 19 5 1 4.720 3.270 54 54 6 2 6.350 1.640 1.630 1.630 16 16 7 1 7.550 1.640 7 7 8 1 7.550 0.000 10 10 9 1 0.000 0.000 39 39 Region 1 mesh points K L X Y 1 1 0.00000000 0.00000000 1 2 0.00000000 0.172857143 1 3 0.00000000 0.345714286 1 4 0.00000000 0.518571429 1 5 0.00000000 0.691428571 1 6 0.00000000 0.864285714 1 7 0.00000000 1.03714286 1 8 0.00000000 1.21000000 1 9 0.00000000 1.38285714 1 10 0.00000000 1.55571429 1 11 0.00000000 1.72857143 1 12 0.00000000 1.90142857 1 13 0.00000000 2.07428571 1 14 0.00000000 2.24714286 1 15 0.00000000 2.42000000 1 16 0.00000000 2.59285714 1 17 0.00000000 2.76571429 1 18 0.00000000 2.93857143 1 19 0.00000000 3.11142857 1 20 0.00000000 3.28428571 1 21 0.00000000 3.45714286 1 22 0.00000000 3.63000000 1 23 0.00000000 3.80285714 1 24 0.00000000 3.97571429 1 25 0.00000000 4.14857143 1 26 0.00000000 4.32142857 1 27 0.00000000 4.49428571 1 28 0.00000000 4.66714286 1 29 0.00000000 4.84000000 1 30 0.00000000 5.01285714 1 31 0.00000000 5.18571429 1 32 0.00000000 5.35857143 1 33 0.00000000 5.53142857 1 34 0.00000000 5.70428571 1 35 0.00000000 5.87714286 1 36 0.00000000 6.05000000 1 37 0.00000000 6.22285714 1 38 0.00000000 6.39571429 1 39 0.00000000 6.56857143 1 40 0.00000000 6.74142857 1 41 0.00000000 6.91428571 1 42 0.00000000 7.08714286 1 43 0.00000000 7.26000000 1 44 0.00000000 7.43285714 1 45 0.00000000 7.60571429 1 46 0.00000000 7.77857143 1 47 0.00000000 7.95142857 1 48 0.00000000 8.12428571 1 49 0.00000000 8.29714286 1 50 0.00000000 8.47000000 1 51 0.00000000 8.64285714 1 52 0.00000000 8.81571429 1 53 0.00000000 8.98857143 1 54 0.00000000 9.16142857 1 55 0.00000000 9.33428571 1 56 0.00000000 9.50714286 1 57 0.00000000 9.68000000 1 58 0.00000000 9.85285714 1 59 0.00000000 10.0257143 1 60 0.00000000 10.1985714 1 61 0.00000000 10.3714286 1 62 0.00000000 10.5442857 1 63 0.00000000 10.7171429 1 64 0.00000000 10.8900000 1 65 0.00000000 11.0628571 1 66 0.00000000 11.2357143 1 67 0.00000000 11.4085714 1 68 0.00000000 11.5814286 1 69 0.00000000 11.7542857 1 70 0.00000000 11.9271429 1 71 0.00000000 12.1000000 1 72 0.00000000 12.2728571 1 73 0.00000000 12.4457143 1 74 0.00000000 12.6185714 1 75 0.00000000 12.7914286 1 76 0.00000000 12.9642857 1 77 0.00000000 13.1371429 1 78 0.00000000 13.3100000 1 79 0.00000000 13.4828571 1 80 0.00000000 13.6557143 1 81 0.00000000 13.8285714 1 82 0.00000000 14.0014286 1 83 0.00000000 14.1742857 1 84 0.00000000 14.3471429 1 85 0.00000000 14.5200000 2 85 0.206153846 14.5200000 3 85 0.412307692 14.5200000 4 85 0.618461538 14.5200000 5 85 0.824615385 14.5200000 6 85 1.03076923 14.5200000 7 85 1.23692308 14.5200000 8 85 1.44307692 14.5200000 9 85 1.64923077 14.5200000 10 85 1.85538462 14.5200000 11 85 2.06153846 14.5200000 12 85 2.26769231 14.5200000 13 85 2.47384615 14.5200000 14 85 2.68000000 14.5200000 15 85 2.87118293 14.5110217 16 85 3.06015962 14.4842651 17 85 3.24482119 14.4402492 18 84 3.42323488 14.3797900 19 84 3.59370659 14.3039354 19 83 3.75480747 14.2138941 20 83 3.90537506 14.1109678 21 82 4.04449598 13.9964929 22 82 4.17148245 13.8717903 22 81 4.28582441 13.7381447 23 80 4.38716540 13.5967750 23 79 4.47525828 13.4488383 24 78 4.54994479 13.2954180 24 77 4.61112643 13.1375338 24 76 4.65874593 12.9761497 24 75 4.69277289 12.8121826 25 74 4.71319286 12.6465136 25 73 4.72000000 12.4800000 25 72 4.72000000 12.3062264 25 71 4.72000000 12.1324528 25 70 4.72000000 11.9586792 25 69 4.72000000 11.7849057 25 68 4.72000000 11.6111321 25 67 4.72000000 11.4373585 25 66 4.72000000 11.2635849 25 65 4.72000000 11.0898113 25 64 4.72000000 10.9160377 25 63 4.72000000 10.7422642 25 62 4.72000000 10.5684906 25 61 4.72000000 10.3947170 25 60 4.72000000 10.2209434 25 59 4.72000000 10.0471698 25 58 4.72000000 9.87339623 25 57 4.72000000 9.69962264 25 56 4.72000000 9.52584906 25 55 4.72000000 9.35207547 25 54 4.72000000 9.17830189 25 53 4.72000000 9.00452830 25 52 4.72000000 8.83075472 25 51 4.72000000 8.65698113 25 50 4.72000000 8.48320755 25 49 4.72000000 8.30943396 25 48 4.72000000 8.13566038 25 47 4.72000000 7.96188679 25 46 4.72000000 7.78811321 25 45 4.72000000 7.61433962 25 44 4.72000000 7.44056604 25 43 4.72000000 7.26679245 25 42 4.72000000 7.09301887 25 41 4.72000000 6.91924528 25 40 4.72000000 6.74547170 25 39 4.72000000 6.57169811 25 38 4.72000000 6.39792453 25 37 4.72000000 6.22415094 25 36 4.72000000 6.05037736 25 35 4.72000000 5.87660377 25 34 4.72000000 5.70283019 25 33 4.72000000 5.52905660 25 32 4.72000000 5.35528302 25 31 4.72000000 5.18150943 25 30 4.72000000 5.00773585 25 29 4.72000000 4.83396226 25 28 4.72000000 4.66018868 25 27 4.72000000 4.48641509 25 26 4.72000000 4.31264151 25 25 4.72000000 4.13886792 25 24 4.72000000 3.96509434 25 23 4.72000000 3.79132075 25 22 4.72000000 3.61754717 25 21 4.72000000 3.44377358 25 20 4.72000000 3.27000000 24 19 4.72783214 3.11040225 25 18 4.75132594 2.95197290 25 17 4.79047027 2.79590401 25 16 4.84523524 2.64343555 25 15 4.91555396 2.49587820 26 15 5.00129494 2.35463414 27 14 5.10222640 2.22120973 27 13 5.21796569 2.09722623 28 12 5.34791896 1.98441313 29 12 5.49120908 1.88458376 29 11 5.64660802 1.79957838 30 11 5.81247360 1.73118053 31 11 5.98672153 1.68099756 32 11 6.16684551 1.65032274 33 10 6.35000000 1.64000000 34 10 6.55000000 1.64000000 35 10 6.75000000 1.64000000 36 10 6.95000000 1.64000000 37 10 7.15000000 1.64000000 38 10 7.35000000 1.64000000 39 10 7.55000000 1.64000000 39 9 7.55000000 1.45777778 39 8 7.55000000 1.27555556 39 7 7.55000000 1.09333333 39 6 7.55000000 0.911111111 39 5 7.55000000 0.728888889 39 4 7.55000000 0.546666667 39 3 7.55000000 0.364444444 39 2 7.55000000 0.182222222 39 1 7.55000000 0.00000000 38 1 7.35131579 0.00000000 37 1 7.15263158 0.00000000 36 1 6.95394737 0.00000000 35 1 6.75526316 0.00000000 34 1 6.55657895 0.00000000 33 1 6.35789474 0.00000000 32 1 6.15921053 0.00000000 31 1 5.96052632 0.00000000 30 1 5.76184211 0.00000000 29 1 5.56315789 0.00000000 28 1 5.36447368 0.00000000 27 1 5.16578947 0.00000000 26 1 4.96710526 0.00000000 25 1 4.76842105 0.00000000 24 1 4.56973684 0.00000000 23 1 4.37105263 0.00000000 22 1 4.17236842 0.00000000 21 1 3.97368421 0.00000000 20 1 3.77500000 0.00000000 19 1 3.57631579 0.00000000 18 1 3.37763158 0.00000000 17 1 3.17894737 0.00000000 16 1 2.98026316 0.00000000 15 1 2.78157895 0.00000000 14 1 2.58289474 0.00000000 13 1 2.38421053 0.00000000 12 1 2.18552632 0.00000000 11 1 1.98684211 0.00000000 10 1 1.78815789 0.00000000 9 1 1.58947368 0.00000000 8 1 1.39078947 0.00000000 7 1 1.19210526 0.00000000 6 1 0.993421053 0.00000000 5 1 0.794736842 0.00000000 4 1 0.596052632 0.00000000 3 1 0.397368421 0.00000000 2 1 0.198684211 0.00000000 1 1 0.00000000 0.00000000 Region 1 done, 9 fixed points, 237 total boundary points. Logical boundary segment end points for region 1: Segment Starting Point Ending Point K L K L 1 1 1 1 85 2 1 85 14 85 3 14 85 25 73 4 25 73 25 20 5 25 20 33 10 6 33 10 39 10 7 39 10 39 1 8 39 1 1 1 Drive point at X =00.0, K = 1, Y = 14.52, L = 85. Region 2 mesh points K L X Y 1 85 0.00000000 14.5200000 Processing boundary data... Memory used for additional setup arrays: 102.856 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 85 1 85 2 85 1 85 98 14 85 3 98 14 85 116 25 73 4 116 25 73 169 25 20 5 169 25 20 184 33 10 6 184 33 10 190 39 10 7 190 39 10 199 39 1 8 199 39 1 237 1 1 Relaxation parameters, 2011 unknown points. Optimizing mesh by successive over relaxation... -------X coordinate-------- -------Y coordinate-------- Cycle Residual ETA RHO Residual ETA RHO 1 2.685E-02 1.0000 1.6000 1.310E-03 1.0000 1.6000 2 3.488E-02 0.6724 1.6000 1.599E-03 0.6490 1.6000 3 2.374E-02 0.6805 1.6000 1.086E-03 0.6796 1.6000 4 1.604E-02 0.6764 1.6000 7.507E-04 0.6910 1.6000 5 1.094E-02 0.6817 1.6000 5.252E-04 0.6995 1.6000 6 7.407E-03 0.6777 1.6000 3.709E-04 0.7062 1.6000 7 5.065E-03 0.6839 1.6000 2.644E-04 0.7130 1.6000 8 3.448E-03 0.6809 1.6000 1.908E-04 0.7217 1.6000 9 2.370E-03 0.6875 1.6000 1.393E-04 0.7300 1.6000 10 1.629E-03 0.6873 1.6000 1.035E-04 0.7433 1.6000 11 1.154E-03 0.7085 1.6000 8.054E-05 0.7778 1.6000 12 8.143E-04 0.7058 1.6000 6.033E-05 0.7491 1.6000 13 5.970E-04 0.7332 1.6000 4.824E-05 0.7995 1.6000 14 4.477E-04 0.7500 1.6000 3.816E-05 0.7910 1.6000 15 3.574E-04 0.7984 1.6000 3.225E-05 0.8452 1.6000 16 2.819E-04 0.7888 1.6000 2.781E-05 0.8624 1.6000 17 2.388E-04 0.8472 1.6000 2.421E-05 0.8705 1.6000 18 2.046E-04 0.8567 1.6000 2.115E-05 0.8737 1.6000 19 1.812E-04 0.8857 1.6000 1.852E-05 0.8757 1.6000 20 1.585E-04 0.8749 1.6000 1.630E-05 0.8801 1.6000 21 1.427E-04 0.9000 1.6000 1.440E-05 0.8836 1.6000 22 1.272E-04 0.8920 1.6000 1.278E-05 0.8870 1.6000 23 1.152E-04 0.9051 1.6000 1.137E-05 0.8901 1.6000 24 1.038E-04 0.9013 1.6000 1.016E-05 0.8931 1.6000 25 9.427E-05 0.9082 1.6000 9.100E-06 0.8960 1.6000 26 1.314E-04 1.3939 1.7445 1.184E-05 1.3009 1.7313 27 1.168E-04 0.8887 1.7445 1.043E-05 0.8811 1.7313 28 1.033E-04 0.8844 1.7445 9.190E-06 0.8810 1.7313 29 9.138E-05 0.8849 1.7445 8.107E-06 0.8822 1.7313 30 8.067E-05 0.8828 1.7445 7.169E-06 0.8842 1.7313 31 7.129E-05 0.8838 1.7445 6.358E-06 0.8869 1.7313 32 6.296E-05 0.8831 1.7445 5.660E-06 0.8902 1.7313 33 5.569E-05 0.8845 1.7445 5.060E-06 0.8939 1.7313 34 4.926E-05 0.8846 1.7445 4.541E-06 0.8974 1.7313 35 4.364E-05 0.8859 1.7445 4.090E-06 0.9006 1.7313 36 3.869E-05 0.8866 1.7445 3.696E-06 0.9039 1.7313 37 3.435E-05 0.8878 1.7445 3.353E-06 0.9070 1.7313 38 3.051E-05 0.8882 1.7445 3.050E-06 0.9097 1.7313 39 2.714E-05 0.8895 1.7445 2.783E-06 0.9123 1.7313 40 2.417E-05 0.8907 1.7445 2.545E-06 0.9146 1.7313 41 2.155E-05 0.8916 1.7445 2.333E-06 0.9166 1.7313 42 1.924E-05 0.8926 1.7445 2.142E-06 0.9183 1.7313 43 1.719E-05 0.8939 1.7445 1.971E-06 0.9199 1.7313 44 1.539E-05 0.8948 1.7445 1.815E-06 0.9211 1.7313 45 1.379E-05 0.8960 1.7445 1.674E-06 0.9222 1.7313 46 1.237E-05 0.8972 1.7445 1.545E-06 0.9230 1.7313 47 1.111E-05 0.8982 1.7445 1.427E-06 0.9237 1.7313 48 9.992E-06 0.8994 1.7445 1.319E-06 0.9243 1.7313 Iteration converged. Generation completed. Problem description: 0.805-GHz TM010 Modified Pillbox Cavity Pillbox cavity with a 1.63-cm-radius bore tube [Originally appeared in 1987 User's Guide 1.4] Problem file: C:\LANL\Examples\RadioFrequency\PillboxCavities\805ROMA.AF 11-17-2007 14:57:34 Problem file length: 820 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 A 0.950000000 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.198684211 First X mesh interval (at XMIN) DXMIN 0.198684211 Minimum X mesh interval (found by Automesh) DYMIN 0.172857143 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 A 1 0 for X,Y problems, 1 for Z,R problems IMAX 41 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 3567 (KMAX+2)*(LMAX+2) KDRI 1 K coordinate of the drive point KMAX 39 Number of horizontal logical mesh points KMETHOD A 1 Wavenumber computation method (1= use beta) KPROB A 1 Problem type indicator (Superfish) LAST35 0 Code for last program to update T35 file LCYCLE 48 Iteration number in mesh optimization LDRI 85 L coordinate of the drive point LINT 1 Logical-mesh coordinate for Ez integration LMAX 85 Number of vertical logical mesh points MAXCY -1 Maximum number of cycles (-1: use default) MAXPPR 238 Maximum points per region NAIR 2207 Number of air points NBND 39 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 3487 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 237 Total number of boundary points in the mesh NPINP 2247 Total points in problem NPONTS 2011 Number of unknown relaxation points NREG 2 Number of regions NRMSEG 1 Normalization segment number for NORM=2 NSEG 8 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.711756368E-02 Average area of all triangles TRIMAX 2.376455740E-02 Area of the largest positive-area triangle TRIMIN 8.329824235E-03 Area of the smallest positive-area triangle VOLUME 3078.12075 Cavity volume (cylindrical symmetry only) XDRI A 0.0 X coordinate of the drive point XMAXG 7.55000000 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 72.8523510 Total cross sectional area YDRI A 14.5200000 Y coordinate of the drive point YMAXG 14.5200000 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