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.16 released 5-9-2005 Starting from file C:\my_poisson\GLUEX_SOL_02_B5.T35. Program file: C:\LANL\FORCE.EXE SF.INI file: C:\LANL\SF.INI 7-21-2006 18:34:36 Memory used for GLUEX_SOL_02_B5.T35 arrays: 15.756 M Memory used for other FORCE arrays: 332.020 K Problem description: GlueX (Hall D) solenoid - basic, 1500A Problem file: C:\my_poisson\GLUEX_SOL_02_B5.AM 7-29-2009 15:21:50 Problem file length: 6864 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.998107426 Rate of convergence in air, cycle N ETAFE 0.998099797 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 4360 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 173650 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 46409 Number of iron points NGMAX 47826 Number of points for recalculating gammas NGSAM 46895 Number of points for reluctivity when NM6=NM1 NINTER 1905 Number of interface points NMATR 1 Number of material records in T35 file NOTE 1 Determines relaxation order NPBOUND 8310 Total number of boundary points in the mesh NPINP 224000 Total points in problem NPONTS 213678 Number of unknown relaxation points NPTC 0 Number of arc points for interpolation NREG 32 Number of regions NSEG 149 Number of boundary segments NSPL 0 Number of special-potential points NTERM 0 Number of coefficients in harmonic analysis NWMAX 48076 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.364034500E-07 Residual for air RESIDI 4.979442540E-07 Residual for iron RHOAIR 1.99514354 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.8161560 Area of the largest positive-area triangle TRIMIN 0.118884112 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 12 Current density = 2452.4244 Integration path, region of interest = 12 Point K L X Y Triangle leg to next point 1 94 172 101.88100 42.00700 space both sides 2 95 172 102.37346 42.00700 space both sides 3 96 172 102.86592 42.00700 space both sides 4 97 172 103.35837 42.00700 space both sides 5 98 172 103.85083 42.00700 space both sides 6 99 172 104.34329 42.00700 space both sides 7 100 172 104.83575 42.00700 space both sides 8 101 172 105.32821 42.00700 space both sides 9 102 172 105.82067 42.00700 space both sides 10 103 172 106.31312 42.00700 space both sides 11 104 172 106.80558 42.00700 space both sides 12 105 172 107.29804 42.00700 space both sides 13 106 172 107.79050 42.00700 space both sides 14 107 172 108.28296 42.00700 space both sides 15 108 172 108.77542 42.00700 space both sides 16 109 172 109.26788 42.00700 space both sides 17 110 172 109.76033 42.00700 space both sides 18 111 172 110.25279 42.00700 space both sides 19 112 172 110.74525 42.00700 space both sides 20 113 172 111.23771 42.00700 space both sides 21 114 172 111.73017 42.00700 space both sides 22 115 172 112.22262 42.00700 space both sides 23 116 172 112.71508 42.00700 space both sides 24 117 172 113.20754 42.00700 space both sides 25 118 172 113.70000 42.00700 space both sides 26 119 172 114.19700 42.00700 space both sides 27 120 172 114.69400 42.00700 space both sides 28 121 172 115.19100 42.00700 space both sides 29 122 172 115.68800 42.00700 space both sides 30 122 173 115.68800 42.98158 space both sides 31 122 174 115.68800 43.95617 space both sides 32 122 175 115.68800 44.93075 space both sides 33 122 176 115.68800 45.90533 space both sides 34 122 177 115.68800 46.87992 space both sides 35 122 178 115.68800 47.85450 space both sides 36 122 179 115.68800 48.82908 space both sides 37 122 180 115.68800 49.80367 space both sides 38 122 181 115.68800 50.77825 space both sides 39 122 182 115.68800 51.75283 space both sides 40 122 183 115.68800 52.72742 space both sides 41 122 184 115.68800 53.70200 space both sides 42 121 184 115.19489 53.70200 space both sides 43 120 184 114.70179 53.70200 space both sides 44 119 184 114.20868 53.70200 space both sides 45 118 184 113.71557 53.70200 space both sides 46 117 184 113.22246 53.70200 space both sides 47 116 184 112.72936 53.70200 space both sides 48 115 184 112.23625 53.70200 space both sides 49 114 184 111.74314 53.70200 space both sides 50 113 184 111.25004 53.70200 space both sides 51 112 184 110.75693 53.70200 space both sides 52 111 184 110.26382 53.70200 space both sides 53 110 184 109.77071 53.70200 space both sides 54 109 184 109.27761 53.70200 space both sides 55 108 184 108.78450 53.70200 space both sides 56 107 184 108.29139 53.70200 space both sides 57 106 184 107.79829 53.70200 space both sides 58 105 184 107.30518 53.70200 space both sides 59 104 184 106.81207 53.70200 space both sides 60 103 184 106.31896 53.70200 space both sides 61 102 184 105.82586 53.70200 space both sides 62 101 184 105.33275 53.70200 space both sides 63 100 184 104.83964 53.70200 space both sides 64 99 184 104.34654 53.70200 space both sides 65 98 184 103.85343 53.70200 space both sides 66 97 184 103.36032 53.70200 space both sides 67 96 184 102.86721 53.70200 space both sides 68 95 184 102.37411 53.70200 space both sides 69 94 184 101.88100 53.70200 space both sides 70 94 183 101.88100 52.72742 space both sides 71 94 182 101.88100 51.75283 space both sides 72 94 181 101.88100 50.77825 space both sides 73 94 180 101.88100 49.80367 space both sides 74 94 179 101.88100 48.82908 space both sides 75 94 178 101.88100 47.85450 space both sides 76 94 177 101.88100 46.87992 space both sides 77 94 176 101.88100 45.90533 space both sides 78 94 175 101.88100 44.93075 space both sides 79 94 174 101.88100 43.95617 space both sides 80 94 173 101.88100 42.98158 space both sides 81 94 172 101.88100 42.00700 Fr = 4.59786E+05 N/radian Fz =-2.57892E+05 N/radian 2piFz =-1.62038E+06 N Fr = 4.68851E+04 kg/radian Fz =-2.62977E+04 kg/radian 2piFz =-1.65233E+05 kg Fr = 1.03382E+05 lb/radian Fz =-5.79864E+04 lb/radian 2piFz =-3.64339E+05 lb