Manuals / IBM / Communications / PDAs & Smartphones

IBM Release 1.93 manual zhi = zsize-Tpml-step xinside=xsize/2 yinside=ysize/2 zinside=zsize/2

Models: Release 1.93

1 93

Download 93 pages 50.46 Kb

Page 62

Image 62

zhi = zsize-Tpml-step ,

xinside=xsize/2 yinside=ysize/2 zinside=zsize/2

FLUX name=MiddleFlux , /* crossing the middle (z=zsize/2) xLo = Tpml ,

xhi = xsize-Tpml, yLo = Tpml , yhi = ysize-Tpml , zlo = zsize/2 , zhi = zsize/2 ,

xinside=xsize/2 yinside=ysize/2 zinside=zsize

INTEGRAL name=diss variable=dissipation, xlo=(xsize-LoadLength)/2 xhi=(xsize+LoadLength)/2, ylo = (ysize-LoadHeight)/2 yhi=(ysize+LoadHeight)/2,

zlo = Zvertex-LoadThickness/2 zhi=Zvertex+LoadThickness/2,

SLICE orientation=xy variable=indexn , xlo=0 xhi=xsize,

ylo=0 yhi=ysize,

zlo=zsize/2+step zhi=zsize/2+step, file=Dipole2IIndexXY.bmp phase=0

SLICE orientation=zx variable=indexn , xlo=0 xhi=xsize, ylo=ysize/2+step yhi=ysize/2+step, zlo=0 zhi=zsize, file=Dipole2IIndexZX.bmp phase=0

SLICE orientation=xy variable=PoyntingZ , xlo=0 xhi=xsize,

ylo=0 yhi=ysize, zlo=zsize/2 zhi=zsize/2,

file=Dipole2IPZ_XY.bmp phase=0 palette=redblue

Page 62

Image 62

IBM Release 1.93 manual zhi = zsize-Tpml-step xinside=xsize/2 yinside=ysize/2 zinside=zsize/2, ylo=0 yhi=ysize

Contents

IBM T. J. Watson Research Center Yorktown Heights, NY January 22 Updated January 23 Phil HobbsA Programmable Optimizing Electro-Magnetic Simulator Release ?88,.d88b, d8888b d8888bPage IBM T. J. Watson Research Center Yorktown Heights, NY A Programmable Optimizing Electro-Magnetic Simulator ReleaseAugust 21 Phil Hobbs POEMSChapter Appendix A. V-Antenna Optimization RunChapter Introduction 1.1. Motivation1.1.1. Philosophy 1.1.2. Structure 1.1.3. OptimizationPage FDTD Engine FIDO + CLUSTER SCRIPT or TEMPEST Optimization LoopFront-End Script POEMS.CMD Postprocessor2.2. The Front-End Script poems.cmd 2. HOW POEMS WORKS2.1. Program Organization 2.2.1. Script Operation2.3. The FDTD Engine FIDO/TEMPEST 2.4. The Postprocessor EMPOSTempost ordersfile resultsfile 2.5. The Visualization System VIS5D2.6. Cluster Control 2.6.1. Parallel Processing 3.1.1. poems Command-Line Options 3. USING POEMS3.1. Command Reference 3.1.2. GLOBAL GroupFREQ 200*THz aliases in different SUBDOMAIN statements in the parameters file NB Hosts are identified by their hostname as specified by the$HOSTNAME environment variable, rather than by IP address. This which means that the host’s predefined hostname is not usedMACRO Syntax PRINT string expression 3.1.3. WORLD Group 3.1.4. MATERIAL Group Syntax XRANGE minexpression maxexpressionDEFINE matname=copper type=conductor conductivity=5.8e7 Parameters epsReal epsImag muReal muImagDEFINE matname= BlackGlass type=dielectric n=1.52 k=1e-3 3.1.5. OBJECT Group xhi1,yhi1,zhi1 and xlo2,ylo2,zlo2, xhi2,yhi2,zhi2. These must beParameters matname orientation xlo xhi ylo yhi zlo zhi widthfunc phasefunc intervalxhi2 ylo2 yhi2 zlo2 zhi2 Parameters mattype xlo xhi ylo yhi zlo zhi thicknessNot implemented in this release Parameters mattype xlo xhi ylo yhi zlo zhiParameters mattype slo shi perp hfunc vfunc Parameters mattype axis lo hi maxradius insidefunc3.1.6. SOURCE Group Parameters x y z width kx ky kz Ex Ey Ez mag phase n k orientationParameters x y z kx ky kz Ex Ey Ez mag phase orientation ParametersNot implemented in this release3.1.8. OUTPUT Group 3.1.7. COMMAND GroupParameters x y z Ex Ey Ez mag phase Parameters reltolerance maxcycles mincycles timestepParameters variable xlo xhi ylo yhi zlo zhi x y z file phase state3.1.9. POSTPROCESS Group towards from xInside, yInside, zInside supplied plane, xmin ≤ x ≤ xmax, ymin ≤ y≤ ymax, zmin ≤ z ≤ zmaxParameters name xlo xhi ylo yhi zlo zhi file direction Parameters name xinside yinside zinside xlo xhi ylo yhi zlo zhiParameters file orientation xlo xhi ylo yhi zlo zhi phase indexn The value of coordinates can be polar or rectangular, and arrangementvariable6 variable7 variable8 variable9 xlo xhi ylo yhi zlo zhi file Parametersfile orientation xlo xhi ylo yhi zlo zhi dt frames variableParameters frames variable1 variable2 variable3 variable4 variable5 MOVIEDISSIPATION Parameters name xlo xhi ylo yhi zlo zhi file xinside yinside zinsidelevel level3.1.10. OPTIMIZE Group arbitrary boundary. Not implemented in this release not be updated at postprocessing time. That means that theSyntax GUESS var1 = expression var2 =, expression Parameters iterations tolerance restarts 3.1.11. SCHEDULE GroupPENALTY expression Stepping is not implemented in this releaseUnfolded 2N-2 Rows 3.2. The Computational Domain3.2.1. Symmetry Actual Domain N rows Image Domain N-2 Rows Next Period3.3. OBJECTS 3.5. SOURCES 3.3.1. Perfectly-Matched Layers3.4. MATERIALS 3.5.1. POINT SOURCES3.5.2. PLANE WAVES Page 3.5.3. BEAM SOURCES 3.5.4. MODEFILE SOURCES3.6. Optimization 3.6.1. Merit Functionsd. RMS Voltage across terminal points c. Phase uniformity across a planee. Power dissipation in a load dissipation f. Reflection from an interface3.6.2. Worked Example Optimizing a V Antenna F i g u r e 2 . 1 2 Optimized V-antennaPage 3.7. Predefined Constants 3.6.3. Worked Example Doped Silica Waveguide Mode3.6.4. Worked Example Glass Ridge Waveguide to Free Space Coupler 3.7.1. Reserved Names i l l u m i n a t i o n p l a n e a f t e ra. Differences from C/C++ 3.8. Predefined Mathematical Functions3.8.1. Arithmetic Operators 3.8.4. ACOS 3.8.2. Logical Operators3.8.3. ABS 3.8.5. ACOSH3.8.12. COS 3.8.10. ATAN23.8.11. CEIL 3.8.13. ELINTK3.8.21. MAX 3.8.19. INTEGRAL3.8.20. LN 3.8.22. MIN3.8.28. ROUND 3.8.26. RANDOM3.8.27. ROOT1D 3.8.29. SIGN3.9. Analytical Pupil Functions 3.10. Material Parameter Functions2. Time step A. TEMPEST and General FDTD Information1. Startup and Steady State dt=2.0*PI/omega*FLOATperiodstep printftime step after adjustment %e\n,dtfflushstdoutperiodstep+=2 i=periodstep%4 if i!=0 periodstep=periodstep+4-i fprintfstderr,periodstep %d\n,periodstepAppendix A. V-Antenna Optimization Run The run will produce lots of bitmaps at appropriately chosen locationsA.1. POEMS INPUT DIPOLE2I.PAR set dipoleheight=0.2*um set dipolethickness=0.30*um Fundamental dimensionsset phi=0 /* phase of excitation */ set ky0=step/lambda set kx0=0 /* direction of plane wave set loadlength=2*stepk=kload DEFINE matname= Loadtype=dielectric n=nload END WORLD SUBDOMAIN ALL XRANGE 0 xsize YRANGE 0 ysize ZRANGE 0 zsizezlo2 = Zvertex-DipoleThickness/2 zhi2= Zvertex+DipoleThickness/2 END COMMAND COMPUTE maxcycles = 25 reltolerance = END OUTPUTzhi=zsize phase=0 state=steady file=Dipole2Iexi FAN matname=PCL taper=linear taperpar=0END OPTIMIZE STORE Dipole2I.simplex MERIT diss*1e12POSTPROCESS PHASEEXxlo=0 xhi=xsize ylo=0 yhi=ysize zlo=0 zhi=zsize file=Dipole2IPhaseEx AMPLEXzhi = zsize-Tpml-step xinside=xsize/2 yinside=ysize/2 zinside=zsize/2 zlo=zsize/2+step zhi=zsize/2+step, file=Dipole2IIndexXY.bmp phase=0file=Dipole2IPZXY.bmp phase=0 palette=redblue FLUX name=MiddleFlux , /* crossing the middle z=zsize/2 xLo = Tpmlfile=Dipole2IPXZX.bmp phase=0 palette=redblue file=Dipole2IPZZX.bmp phase=0 palette=redbluefile=Dipole2IPXXY.bmp phase=0 palette=redblue file=Dipole2IAmpExZX.bmp phase=0 palette=flamefile=Dipole2IExIZX.bmp phase=0 palette=redblue SLICE orientation=xy variable=PhaseEx , xlo=0 xhi=xsizeSLICE orientation=zx variable=PhaseEx , xlo=0 xhi=xsize file=Dipole2IDissXY.bmp palette=redbluefile=Dipole2IExQZX.bmp phase=pi/2 palette=redblue END ylo=ysize/2 yhi=ysize/2, zlo=0 zhi=zsize0 100 index 1.000000 A.2. tempest Input File DIPOLE2I.PAR.INrectangle node plot hy steady 0.250000 node 0 99 0 99 0 99 dipole2ihyqWritten by 56 pages of pointsource statements omitted A.3. Postprocessor orders DIPOLE2I.ORDERS empost dipole2i.orders outfile v5dfileDIPOLE2IEYI NULLEXI ARRAY Domain for comparison e.g. modematch 0.00000 SLICE POSTPROC.16.COMPARISONDOMAIN 0.00000 SLICE A.4. Run Results DIPOLE2I.SIMPLEX 0.00000069266248611 90.849532778 0.00000152214633334 Penalty = Contracting by 0.5 in 1D 1-D shrink unsuccessful0.00000064767135 98.83191 0.000001519205 Penalty = 0.000000883423566665 77.7563866665 0.0000017351705 Penalty =Reflecting through centroid. Contracting by 0.5 in 1D Reflecting through centroid. Contracting by 0.5 in 1DReflecting through centroid. Contracting by 0.5 in 1D Reflecting through centroid. Contracting by 0.5 in 1DReflecting through centroid. Contracting by 0.5 in 1D Reflecting through centroid. Contracting by 0.5 in 1DReflecting through centroid. Contracting by 0.5 in 1D Failed to converge in 36 function evaluations1-D shrink unsuccessful Reflecting through centroid. Contracting by 0.5 in 1DReflecting through centroid. Contracting by 0.5 in 1D 0.000000572297880735 102.392565535 0.00000148853488604 Penalty =Converged in 56 function evaluations Reflecting through centroid. Contracting by 0.5 in 1DReflecting through centroid. Contracting by 0.5 in 1D B.1.1. tempest patches Appendix B. ConfigurationB.1. FDTD and TEMPEST b. Understanding convergence B.1.3. Advice common to all or most FDTD programsB.1.2. tempest limitations c. Hints and KinksB.2. REXX B.3. X Window System ConfigurationB.3.1. Sample X11 Configuration B.4.1. Wish list B.3.2. Running Vis5DB.4. Release Notes B.4.2. Beta Release Limitations ZX plane Y perpendicular +Z right, +X up, +Y out of screen XY plane Z perpendicular +X right, +Y up, +Z out of screenYZ plane X perpendicular +Y right, +Z up, +X out of screen Appendix C. Multiprocessing BASICSTEP 16, 18, 24, 54, 62 Boundary 7, 15, 19, 29, 31 INDEXAiry pattern Antenna 1, 39, 40, 45, 51, 52 C++ 7Periodic boundary conditions 31, 34