Los Alamos National Laboratory Poisson Superfish Program Force written by Ron F. Holsinger in collaboration with Klaus Halbach 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 Force 7.17 released 1-13-2006 Starting from file Z:\home\gen\.wine\my_poisson\GLUEX_SOL_02.T35. Program file: Z:\HOME\GEN\.WINE\LANL\FORCE.EXE SF.INI file: Z:\HOME\GEN\.WINE\LANL\SF.INI 7-21-2006 18:34:36 Memory used for GLUEX_SOL_02.T35 arrays: 15.759 M Memory used for other FORCE arrays: 332.020 K Problem description: GlueX (Hall D) solenoid - basic, 1500A Problem file: Z:\home\gen\.wine\my_poisson\GLUEX_SOL_02.AM 7-29-2009 9:49:50 Problem file length: 6691 bytes Originating program: Automesh Problem type: Magnetostatic Coordinates and lengths have dimensions of centimeters. Problem constants and variables. Letter A in the code column indicates a value supplied in the Automesh input file. Letter S in the code column indicates a value supplied in the SFO input file. Variable Code Value Description ANGLE 0.0 Extent of arc for interpolating potential ANGLZ 0.0 Initial point on arc for interpolation CLIGHT 2.997924580E+10 Exact speed of light in cm/sec CONV 1.00000000 Length conversion (number of units per cm) DIAGDLL 0 If 1, DLL writes diagnostics to DiagDLL.txt DX1 2.06896552 First X mesh interval (at XMIN) DXMIN 0.250000000 Minimum X mesh interval (found by Automesh) DYMIN 1.00000000 Minimum Y mesh interval (found by Automesh) EPS0 8.854187818E-12 Permittivity of free space EPSILA 5.000000000E-07 Convergence PARAMETER: air, interface points EPSILI 5.000000000E-07 Convergence parameter for iron points EPSO 1.000000000E-05 Convergence parameter in mesh optimization ETAAIR 0.998128892 Rate of convergence in air, cycle N ETAFE 0.998121870 Rate of convergence in iron, cycle N FIXEPS 9.00000000 Permittivity when MODE=-1 FIXGAM 4.000000000E-03 Reluctivity when MODE=-1 FMU0 1.256637061E-06 Permeability of free space ICAL 0 Type of formula for currents near coils ICCC 1 1 in Poisson/Pandira problems (all real arrays) ICYCLE 4400 Present iteration number ICYLIN A 1 0 for X,Y problems, 1 for Z,R problems ICYSEN 0 If zero, do not print boundary integrals IENERGY 0 If 1, calculate stored energy IHDL 100000 Cycles between H*dl integrals IMAX 322 KMAX+2 IPERM 0 0 = real currents, 1 = no real currents IPIVOT 1 Pivoting in matrix inversion routines IPRFQ 0 Print frequency during Poisson iterations IRMAX 25 Used in optimization of RHOXY ISKIP 1 Cycles between calculation of reluctance ITOT 226044 (KMAX+2)*(LMAX+2) IVERG 10 Number of cycles between convergence tests KBZERO 1 Logical K coordinate where B = BDES KMAX 320 Number of horizontal logical mesh points KMIN 1 Lower K bound for computing field KPROB A 0 Problem type indicator (Poisson) KTOP 320 Upper K bound for computing field KTYPE 2 Symmetry indicator for harmonic analysis LAST35 5 Code for last program to update T35 file LBZERO 1 Logical L coordinate where B = BDES LCYCLE 47 Iteration number LMAX 700 Number of vertical logical mesh points LMIN 1 Lower L bound for computing field LTOP 1 Upper L bound for computing field MAP 1 Conformal mapping parameter MAXCY 100000 Maximum number of cycles (-1: use default) MAXPPR 2038 Maximum points per region MODE A 0 Material property indicator NAIR 173267 Number of air points NAMAX 0 Number of points for recalculating couplings NBND 2036 Number of Dirichlet boundary points NBSLF A 0 Left-side boundary condition NBSLO A 0 Lower boundary condition NBSRT A 0 Right-side boundary condition NBSUP A 0 Upper boundary condition NEGAT 0 Zero-area triangle indicator NFE 46591 Number of iron points NGMAX 48222 Number of points for recalculating gammas NGSAM 47064 Number of points for reluctivity when NM6=NM1 NINTER 2106 Number of interface points NMATR 1 Number of material records in T35 file NOTE 1 Determines relaxation order NPBOUND 8375 Total number of boundary points in the mesh NPINP 224000 Total points in problem NPONTS 213542 Number of unknown relaxation points NPTC 0 Number of arc points for interpolation NREG 32 Number of regions NSEG 151 Number of boundary segments NSPL 0 Number of special-potential points NTERM 0 Number of coefficients in harmonic analysis NWMAX 48459 Number of points for recalculating couplings OMEGAM 1.000000000E-03 Used in optimization of RHOXY OMEGAP 1.000000000E-03 Optimization parameter for RHOAIR in Poisson PI 3.14159265 The number pi to machine precision RATIO 0.0 BZERO/XJFACT for air solution RESIDA 1.385090612E-07 Residual for air RESIDI 4.907815868E-07 Residual for iron RHOAIR 1.99516927 Over-relaxation factor: air, interface points RHOFE 1.00000000 Over-relaxation factor for iron points RHOGAM 8.000000000E-02 Under-relaxation parameter for reluctivity RHOPT1 1.90000000 If = RHOAIR, causes optimization of RHOAIR RHOXY 1.60000000 Over-relaxation factor for mesh optimization RINT 0.0 Radius of arc for interpolating potential RNORM 0.0 Aperture radius in harmonic analysis RZERO 1.00000000 Scaling factor for conformal transformation SNEGA 0.0 Total negative current, this solution SNEGG 0.0 Total negative current at generation SPOSA 6912000.00 total positive current, this solution SPOSG 6912000.00 Total positive current at generation STACK 1.00000000 Stacking factor STOTA 6912000.00 Total current in problem, this solution STOTG 6912000.00 Total current in problem at generation TNEGC 0.0 Negative current after conformal transformation TPOSC 0.0 Positive current after conformal transformation TRIAVG 0.896937407 Average area of all triangles TRIMAX 10.8161565 Area of the largest positive-area triangle TRIMIN 0.118884109 Area of the smallest positive-area triangle VOLUME 502654825. Cavity volume (cylindrical symmetry only) XJFACT 1.00000000 Factor multiplying currents XJFEND 1.00000000 Ending XJFACT for a range of multipliers XJSTEPS 0 Number of steps from XJFACT to XJFEND XJTOL 5.000000000E-07 Tolerance on XJFACT XMAXF 0.0 Upper X bound for computing field XMAXG 400.000000 Upper X bound of the problem geometry XMINF 0.0 Lower X bound for computing field XMING 0.0 Lower X bound of the problem geometry XORG 0.0 X coordinate center of harmonic analysis circle XYAREA 400000.000 Total cross sectional area YMAXF 0.0 Upper Y bound for computing field YMAXG 700.000000 Upper Y bound of the problem geometry YMINF 0.0 Lower Y bound for computing field YMING -300.000000 Lower Y bound of the problem geometry YORG 0.0 Y coordinate center of harmonic analysis circle Region 18 Current density = 2414.8556 Integration path, region of interest = 18 Point K L X Y Triangle leg to next point 1 94 441 101.88100 310.70900 space both sides 2 95 441 102.38006 310.70900 space both sides 3 96 441 102.87911 310.70900 space both sides 4 97 441 103.37817 310.70900 space both sides 5 98 441 103.87723 310.70900 space both sides 6 99 441 104.37628 310.70900 space both sides 7 100 441 104.87534 310.70900 space both sides 8 101 441 105.37440 310.70900 space both sides 9 102 441 105.87345 310.70900 space both sides 10 103 441 106.37251 310.70900 space both sides 11 104 441 106.87157 310.70900 space both sides 12 105 441 107.37062 310.70900 space both sides 13 106 441 107.86968 310.70900 space both sides 14 107 441 108.36874 310.70900 space both sides 15 108 441 108.86779 310.70900 space both sides 16 109 441 109.36685 310.70900 space both sides 17 110 441 109.86591 310.70900 space both sides 18 111 441 110.36496 310.70900 space both sides 19 112 441 110.86402 310.70900 space both sides 20 113 441 111.36307 310.70900 space both sides 21 114 441 111.86213 310.70900 space both sides 22 115 441 112.36119 310.70900 space both sides 23 116 441 112.86025 310.70900 space both sides 24 117 441 113.35930 310.70900 space both sides 25 118 441 113.85836 310.70900 space both sides 26 119 441 114.35742 310.70900 space both sides 27 120 441 114.85647 310.70900 space both sides 28 121 441 115.35553 310.70900 space both sides 29 122 441 115.85458 310.70900 space both sides 30 123 441 116.35364 310.70900 space both sides 31 124 441 116.85270 310.70900 space both sides 32 125 441 117.35175 310.70900 space both sides 33 126 441 117.85081 310.70900 space both sides 34 127 441 118.34987 310.70900 space both sides 35 128 441 118.84892 310.70900 space both sides 36 129 441 119.34798 310.70900 space both sides 37 130 441 119.84704 310.70900 space both sides 38 131 441 120.34609 310.70900 space both sides 39 132 441 120.84515 310.70900 space both sides 40 133 441 121.34421 310.70900 space both sides 41 134 441 121.84326 310.70900 space both sides 42 135 441 122.34232 310.70900 space both sides 43 136 441 122.84138 310.70900 space both sides 44 137 441 123.34043 310.70900 space both sides 45 138 441 123.83949 310.70900 space both sides 46 139 441 124.33855 310.70900 space both sides 47 140 441 124.83760 310.70900 space both sides 48 141 441 125.33666 310.70900 space both sides 49 142 441 125.83572 310.70900 space both sides 50 143 441 126.33477 310.70900 space both sides 51 144 441 126.83383 310.70900 space both sides 52 145 441 127.33289 310.70900 space both sides 53 146 441 127.83194 310.70900 space both sides 54 147 441 128.33100 310.70900 space both sides 55 147 442 128.33100 311.72935 space both sides 56 147 443 128.33100 312.74969 space both sides 57 147 444 128.33100 313.77003 space both sides 58 147 445 128.33100 314.79038 space both sides 59 147 446 128.33100 315.81072 space both sides 60 147 447 128.33100 316.83107 space both sides 61 147 448 128.33100 317.85141 space both sides 62 147 449 128.33100 318.87176 space both sides 63 147 450 128.33100 319.89210 space both sides 64 147 451 128.33100 320.91245 space both sides 65 147 452 128.33100 321.93279 space both sides 66 147 453 128.33100 322.95314 space both sides 67 147 454 128.33100 323.97348 space both sides 68 147 455 128.33100 324.99383 space both sides 69 147 456 128.33100 326.01417 space both sides 70 147 457 128.33100 327.03452 space both sides 71 147 458 128.33100 328.05486 space both sides 72 147 459 128.33100 329.07521 space both sides 73 147 460 128.33100 330.09555 space both sides 74 147 461 128.33100 331.11590 space both sides 75 147 462 128.33100 332.13624 space both sides 76 147 463 128.33100 333.15659 space both sides 77 147 464 128.33100 334.17693 space both sides 78 147 465 128.33100 335.19728 space both sides 79 147 466 128.33100 336.21762 space both sides 80 147 467 128.33100 337.23797 space both sides 81 147 468 128.33100 338.25831 space both sides 82 147 469 128.33100 339.27866 space both sides 83 147 470 128.33100 340.29900 space both sides 84 146 470 127.83194 340.29900 space both sides 85 145 470 127.33289 340.29900 space both sides 86 144 470 126.83383 340.29900 space both sides 87 143 470 126.33477 340.29900 space both sides 88 142 470 125.83572 340.29900 space both sides 89 141 470 125.33666 340.29900 space both sides 90 140 470 124.83760 340.29900 space both sides 91 139 470 124.33855 340.29900 space both sides 92 138 470 123.83949 340.29900 space both sides 93 137 470 123.34043 340.29900 space both sides 94 136 470 122.84138 340.29900 space both sides 95 135 470 122.34232 340.29900 space both sides 96 134 470 121.84326 340.29900 space both sides 97 133 470 121.34421 340.29900 space both sides 98 132 470 120.84515 340.29900 space both sides 99 131 470 120.34609 340.29900 space both sides 100 130 470 119.84704 340.29900 space both sides 101 129 470 119.34798 340.29900 space both sides 102 128 470 118.84892 340.29900 space both sides 103 127 470 118.34987 340.29900 space both sides 104 126 470 117.85081 340.29900 space both sides 105 125 470 117.35175 340.29900 space both sides 106 124 470 116.85270 340.29900 space both sides 107 123 470 116.35364 340.29900 space both sides 108 122 470 115.85458 340.29900 space both sides 109 121 470 115.35553 340.29900 space both sides 110 120 470 114.85647 340.29900 space both sides 111 119 470 114.35742 340.29900 space both sides 112 118 470 113.85836 340.29900 space both sides 113 117 470 113.35930 340.29900 space both sides 114 116 470 112.86025 340.29900 space both sides 115 115 470 112.36119 340.29900 space both sides 116 114 470 111.86213 340.29900 space both sides 117 113 470 111.36307 340.29900 space both sides 118 112 470 110.86402 340.29900 space both sides 119 111 470 110.36496 340.29900 space both sides 120 110 470 109.86591 340.29900 space both sides 121 109 470 109.36685 340.29900 space both sides 122 108 470 108.86779 340.29900 space both sides 123 107 470 108.36874 340.29900 space both sides 124 106 470 107.86968 340.29900 space both sides 125 105 470 107.37062 340.29900 space both sides 126 104 470 106.87157 340.29900 space both sides 127 103 470 106.37251 340.29900 space both sides 128 102 470 105.87345 340.29900 space both sides 129 101 470 105.37440 340.29900 space both sides 130 100 470 104.87534 340.29900 space both sides 131 99 470 104.37628 340.29900 space both sides 132 98 470 103.87723 340.29900 space both sides 133 97 470 103.37817 340.29900 space both sides 134 96 470 102.87911 340.29900 space both sides 135 95 470 102.38006 340.29900 space both sides 136 94 470 101.88100 340.29900 space both sides 137 94 469 101.88100 339.27866 space both sides 138 94 468 101.88100 338.25831 space both sides 139 94 467 101.88100 337.23797 space both sides 140 94 466 101.88100 336.21762 space both sides 141 94 465 101.88100 335.19728 space both sides 142 94 464 101.88100 334.17693 space both sides 143 94 463 101.88100 333.15659 space both sides 144 94 462 101.88100 332.13624 space both sides 145 94 461 101.88100 331.11590 space both sides 146 94 460 101.88100 330.09555 space both sides 147 94 459 101.88100 329.07521 space both sides 148 94 458 101.88100 328.05486 space both sides 149 94 457 101.88100 327.03452 space both sides 150 94 456 101.88100 326.01417 space both sides 151 94 455 101.88100 324.99383 space both sides 152 94 454 101.88100 323.97348 space both sides 153 94 453 101.88100 322.95314 space both sides 154 94 452 101.88100 321.93279 space both sides 155 94 451 101.88100 320.91245 space both sides 156 94 450 101.88100 319.89210 space both sides 157 94 449 101.88100 318.87176 space both sides 158 94 448 101.88100 317.85141 space both sides 159 94 447 101.88100 316.83107 space both sides 160 94 446 101.88100 315.81072 space both sides 161 94 445 101.88100 314.79038 space both sides 162 94 444 101.88100 313.77003 space both sides 163 94 443 101.88100 312.74969 space both sides 164 94 442 101.88100 311.72935 space both sides 165 94 441 101.88100 310.70900 Fr = 3.16213E+06 N/radian Fz =-3.89200E+05 N/radian 2piFz =-2.44542E+06 N Fr = 3.22448E+05 kg/radian Fz =-3.96874E+04 kg/radian 2piFz =-2.49363E+05 kg Fr = 7.10998E+05 lb/radian Fz =-8.75107E+04 lb/radian 2piFz =-5.49846E+05 lb