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_B3.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_B3.T35 arrays: 15.763 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_B3.AM 7-29-2009 13:05:48 Problem file length: 6858 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.998127899 Rate of convergence in air, cycle N ETAFE 0.998121720 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 4390 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 172770 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 47069 Number of iron points NGMAX 48706 Number of points for recalculating gammas NGSAM 47555 Number of points for reluctivity when NM6=NM1 NINTER 2125 Number of interface points NMATR 1 Number of material records in T35 file NOTE 1 Determines relaxation order NPBOUND 8545 Total number of boundary points in the mesh NPINP 224000 Total points in problem NPONTS 213452 Number of unknown relaxation points NPTC 0 Number of arc points for interpolation NREG 33 Number of regions NSEG 155 Number of boundary segments NSPL 0 Number of special-potential points NTERM 0 Number of coefficients in harmonic analysis NWMAX 48956 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.406782273E-07 Residual for air RESIDI 4.979882317E-07 Residual for iron RHOAIR 1.99516843 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.118884108 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 13 Current density = 2449.8361 Integration path, region of interest = 13 Point K L X Y Triangle leg to next point 1 94 348 101.88100 217.44100 space both sides 2 95 348 102.39150 217.44100 space both sides 3 96 348 102.90200 217.44100 space both sides 4 97 348 103.41250 217.44100 space both sides 5 98 348 103.92300 217.44100 space both sides 6 99 348 104.43350 217.44100 space both sides 7 100 348 104.94400 217.44100 space both sides 8 101 348 105.45450 217.44100 space both sides 9 102 348 105.96500 217.44100 space both sides 10 103 348 106.47550 217.44100 space both sides 11 104 348 106.98600 217.44100 space both sides 12 105 348 107.49650 217.44100 space both sides 13 106 348 108.00700 217.44100 space both sides 14 107 348 108.51750 217.44100 space both sides 15 108 348 109.02800 217.44100 space both sides 16 109 348 109.53850 217.44100 space both sides 17 110 348 110.04900 217.44100 space both sides 18 111 348 110.55950 217.44100 space both sides 19 112 348 111.07000 217.44100 space both sides 20 113 348 111.58050 217.44100 space both sides 21 114 348 112.09100 217.44100 space both sides 22 115 348 112.60150 217.44100 space both sides 23 116 348 113.11200 217.44100 space both sides 24 117 348 113.62250 217.44100 space both sides 25 118 348 114.13300 217.44100 space both sides 26 119 348 114.64350 217.44100 space both sides 27 120 348 115.15400 217.44100 space both sides 28 121 348 115.66450 217.44100 space both sides 29 122 348 116.17500 217.44100 space both sides 30 123 348 116.68550 217.44100 space both sides 31 124 348 117.19600 217.44100 space both sides 32 125 348 117.70650 217.44100 space both sides 33 126 348 118.21700 217.44100 space both sides 34 126 349 118.21700 218.50409 space both sides 35 126 350 118.21700 219.56718 space both sides 36 126 351 118.21700 220.63027 space both sides 37 126 352 118.21700 221.69336 space both sides 38 126 353 118.21700 222.75645 space both sides 39 126 354 118.21700 223.81955 space both sides 40 126 355 118.21700 224.88264 space both sides 41 126 356 118.21700 225.94573 space both sides 42 126 357 118.21700 227.00882 space both sides 43 126 358 118.21700 228.07191 space both sides 44 126 359 118.21700 229.13500 space both sides 45 125 359 117.70650 229.13500 space both sides 46 124 359 117.19600 229.13500 space both sides 47 123 359 116.68550 229.13500 space both sides 48 122 359 116.17500 229.13500 space both sides 49 121 359 115.66450 229.13500 space both sides 50 120 359 115.15400 229.13500 space both sides 51 119 359 114.64350 229.13500 space both sides 52 118 359 114.13300 229.13500 space both sides 53 117 359 113.62250 229.13500 space both sides 54 116 359 113.11200 229.13500 space both sides 55 115 359 112.60150 229.13500 space both sides 56 114 359 112.09100 229.13500 space both sides 57 113 359 111.58050 229.13500 space both sides 58 112 359 111.07000 229.13500 space both sides 59 111 359 110.55950 229.13500 space both sides 60 110 359 110.04900 229.13500 space both sides 61 109 359 109.53850 229.13500 space both sides 62 108 359 109.02800 229.13500 space both sides 63 107 359 108.51750 229.13500 space both sides 64 106 359 108.00700 229.13500 space both sides 65 105 359 107.49650 229.13500 space both sides 66 104 359 106.98600 229.13500 space both sides 67 103 359 106.47550 229.13500 space both sides 68 102 359 105.96500 229.13500 space both sides 69 101 359 105.45450 229.13500 space both sides 70 100 359 104.94400 229.13500 space both sides 71 99 359 104.43350 229.13500 space both sides 72 98 359 103.92300 229.13500 space both sides 73 97 359 103.41250 229.13500 space both sides 74 96 359 102.90200 229.13500 space both sides 75 95 359 102.39150 229.13500 space both sides 76 94 359 101.88100 229.13500 space both sides 77 94 358 101.88100 228.07191 space both sides 78 94 357 101.88100 227.00882 space both sides 79 94 356 101.88100 225.94573 space both sides 80 94 355 101.88100 224.88264 space both sides 81 94 354 101.88100 223.81955 space both sides 82 94 353 101.88100 222.75645 space both sides 83 94 352 101.88100 221.69336 space both sides 84 94 351 101.88100 220.63027 space both sides 85 94 350 101.88100 219.56718 space both sides 86 94 349 101.88100 218.50409 space both sides 87 94 348 101.88100 217.44100 Fr = 6.00421E+05 N/radian Fz = 1.53619E+05 N/radian 2piFz = 9.65217E+05 N Fr = 6.12259E+04 kg/radian Fz = 1.56648E+04 kg/radian 2piFz = 9.84247E+04 kg Fr = 1.35003E+05 lb/radian Fz = 3.45408E+04 lb/radian 2piFz = 2.17027E+05 lb