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 8 Current density = 2400.0150 Integration path, region of interest = 8 Point K L X Y Triangle leg to next point 1 94 310 101.88100 180.13400 space both sides 2 95 310 102.38006 180.13400 space both sides 3 96 310 102.87912 180.13400 space both sides 4 97 310 103.37818 180.13400 space both sides 5 98 310 103.87723 180.13400 space both sides 6 99 310 104.37629 180.13400 space both sides 7 100 310 104.87535 180.13400 space both sides 8 101 310 105.37441 180.13400 space both sides 9 102 310 105.87347 180.13400 space both sides 10 103 310 106.37253 180.13400 space both sides 11 104 310 106.87159 180.13400 space both sides 12 105 310 107.37065 180.13400 space both sides 13 106 310 107.86971 180.13400 space both sides 14 107 310 108.36876 180.13400 space both sides 15 108 310 108.86782 180.13400 space both sides 16 109 310 109.36688 180.13400 space both sides 17 110 310 109.86594 180.13400 space both sides 18 111 310 110.36500 180.13400 space both sides 19 112 310 110.86406 180.13400 space both sides 20 113 310 111.36312 180.13400 space both sides 21 114 310 111.86218 180.13400 space both sides 22 115 310 112.36123 180.13400 space both sides 23 116 310 112.86029 180.13400 space both sides 24 117 310 113.35935 180.13400 space both sides 25 118 310 113.85841 180.13400 space both sides 26 119 310 114.35747 180.13400 space both sides 27 120 310 114.85653 180.13400 space both sides 28 121 310 115.35559 180.13400 space both sides 29 122 310 115.85465 180.13400 space both sides 30 123 310 116.35371 180.13400 space both sides 31 124 310 116.85277 180.13400 space both sides 32 125 310 117.35182 180.13400 space both sides 33 126 310 117.85088 180.13400 space both sides 34 127 310 118.34994 180.13400 space both sides 35 128 310 118.84900 180.13400 space both sides 36 129 310 119.48100 180.13400 space both sides 37 130 310 120.11300 180.13400 space both sides 38 130 311 120.11300 181.12813 space both sides 39 130 312 120.11300 182.12225 space both sides 40 130 313 120.11300 183.11638 space both sides 41 130 314 120.11300 184.11050 space both sides 42 130 315 120.11300 185.10462 space both sides 43 130 316 120.11300 186.09875 space both sides 44 130 317 120.11300 187.09287 space both sides 45 130 318 120.11300 188.08700 space both sides 46 129 318 119.60656 188.08700 space both sides 47 128 318 119.10011 188.08700 space both sides 48 127 318 118.59367 188.08700 space both sides 49 126 318 118.08722 188.08700 space both sides 50 125 318 117.58078 188.08700 space both sides 51 124 318 117.07433 188.08700 space both sides 52 123 318 116.56789 188.08700 space both sides 53 122 318 116.06144 188.08700 space both sides 54 121 318 115.55500 188.08700 space both sides 55 120 318 115.04856 188.08700 space both sides 56 119 318 114.54211 188.08700 space both sides 57 118 318 114.03567 188.08700 space both sides 58 117 318 113.52922 188.08700 space both sides 59 116 318 113.02278 188.08700 space both sides 60 115 318 112.51633 188.08700 space both sides 61 114 318 112.00989 188.08700 space both sides 62 113 318 111.50344 188.08700 space both sides 63 112 318 110.99700 188.08700 space both sides 64 111 318 110.49056 188.08700 space both sides 65 110 318 109.98411 188.08700 space both sides 66 109 318 109.47767 188.08700 space both sides 67 108 318 108.97122 188.08700 space both sides 68 107 318 108.46478 188.08700 space both sides 69 106 318 107.95833 188.08700 space both sides 70 105 318 107.45189 188.08700 space both sides 71 104 318 106.94544 188.08700 space both sides 72 103 318 106.43900 188.08700 space both sides 73 102 318 105.93256 188.08700 space both sides 74 101 318 105.42611 188.08700 space both sides 75 100 318 104.91967 188.08700 space both sides 76 99 318 104.41322 188.08700 space both sides 77 98 318 103.90678 188.08700 space both sides 78 97 318 103.40033 188.08700 space both sides 79 96 318 102.89389 188.08700 space both sides 80 95 318 102.38744 188.08700 space both sides 81 94 318 101.88100 188.08700 space both sides 82 94 317 101.88100 187.09287 space both sides 83 94 316 101.88100 186.09875 space both sides 84 94 315 101.88100 185.10462 space both sides 85 94 314 101.88100 184.11050 space both sides 86 94 313 101.88100 183.11638 space both sides 87 94 312 101.88100 182.12225 space both sides 88 94 311 101.88100 181.12813 space both sides 89 94 310 101.88100 180.13400 Fr = 3.75585E+05 N/radian Fz =-2.47267E+05 N/radian 2piFz =-1.55363E+06 N Fr = 3.82990E+04 kg/radian Fz =-2.52142E+04 kg/radian 2piFz =-1.58426E+05 kg Fr = 8.44493E+04 lb/radian Fz =-5.55974E+04 lb/radian 2piFz =-3.49329E+05 lb