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 17 Current density = 2456.6870 Integration path, region of interest = 17 Point K L X Y Triangle leg to next point 1 94 424 101.88100 293.92900 space both sides 2 95 424 102.39364 293.92900 space both sides 3 96 424 102.90627 293.92900 space both sides 4 97 424 103.41891 293.92900 space both sides 5 98 424 103.93154 293.92900 space both sides 6 99 424 104.44418 293.92900 space both sides 7 100 424 104.95682 293.92900 space both sides 8 101 424 105.46945 293.92900 space both sides 9 102 424 105.98209 293.92900 space both sides 10 103 424 106.49473 293.92900 space both sides 11 104 424 107.00736 293.92900 space both sides 12 105 424 107.52000 293.92900 space both sides 13 106 424 108.03264 293.92900 space both sides 14 107 424 108.54527 293.92900 space both sides 15 108 424 109.05791 293.92900 space both sides 16 109 424 109.57054 293.92900 space both sides 17 110 424 110.08318 293.92900 space both sides 18 111 424 110.59582 293.92900 space both sides 19 112 424 111.10846 293.92900 space both sides 20 113 424 111.62109 293.92900 space both sides 21 114 424 112.13373 293.92900 space both sides 22 115 424 112.64636 293.92900 space both sides 23 116 424 113.15900 293.92900 space both sides 24 116 425 113.15900 294.90350 space both sides 25 116 426 113.15900 295.87800 space both sides 26 116 427 113.15900 296.85250 space both sides 27 116 428 113.15900 297.82700 space both sides 28 116 429 113.15900 298.80150 space both sides 29 116 430 113.15900 299.77600 space both sides 30 116 431 113.15900 300.75050 space both sides 31 116 432 113.15900 301.72500 space both sides 32 116 433 113.15900 302.69950 space both sides 33 116 434 113.15900 303.67400 space both sides 34 116 435 113.15900 304.64850 space both sides 35 116 436 113.15900 305.62300 space both sides 36 115 436 112.64636 305.62300 space both sides 37 114 436 112.13373 305.62300 space both sides 38 113 436 111.62109 305.62300 space both sides 39 112 436 111.10846 305.62300 space both sides 40 111 436 110.59582 305.62300 space both sides 41 110 436 110.08318 305.62300 space both sides 42 109 436 109.57054 305.62300 space both sides 43 108 436 109.05791 305.62300 space both sides 44 107 436 108.54527 305.62300 space both sides 45 106 436 108.03264 305.62300 space both sides 46 105 436 107.52000 305.62300 space both sides 47 104 436 107.00736 305.62300 space both sides 48 103 436 106.49473 305.62300 space both sides 49 102 436 105.98209 305.62300 space both sides 50 101 436 105.46945 305.62300 space both sides 51 100 436 104.95682 305.62300 space both sides 52 99 436 104.44418 305.62300 space both sides 53 98 436 103.93154 305.62300 space both sides 54 97 436 103.41891 305.62300 space both sides 55 96 436 102.90627 305.62300 space both sides 56 95 436 102.39364 305.62300 space both sides 57 94 436 101.88100 305.62300 space both sides 58 94 435 101.88100 304.64850 space both sides 59 94 434 101.88100 303.67400 space both sides 60 94 433 101.88100 302.69950 space both sides 61 94 432 101.88100 301.72500 space both sides 62 94 431 101.88100 300.75050 space both sides 63 94 430 101.88100 299.77600 space both sides 64 94 429 101.88100 298.80150 space both sides 65 94 428 101.88100 297.82700 space both sides 66 94 427 101.88100 296.85250 space both sides 67 94 426 101.88100 295.87800 space both sides 68 94 425 101.88100 294.90350 space both sides 69 94 424 101.88100 293.92900 Fr = 6.18208E+05 N/radian Fz = 3.75667E+05 N/radian 2piFz = 2.36039E+06 N Fr = 6.30397E+04 kg/radian Fz = 3.83074E+04 kg/radian 2piFz = 2.40693E+05 kg Fr = 1.39002E+05 lb/radian Fz = 8.44679E+04 lb/radian 2piFz = 5.30727E+05 lb