Manuals / IBM / Communications / PDAs & Smartphones

IBM Release 1.93 manual A Programmable Optimizing Electro-Magnetic Simulator Release, Poems

Models: Release 1.93

1 93

Download 93 pages 50.46 Kb

Page 3

Image 3

POEMS:

A Programmable Optimizing Electro-Magnetic Simulator

Release 1.70

August 21, 2007

Phil Hobbs

IBM T. J. Watson Research Center Yorktown Heights, NY

iii

Page 3

Image 3

IBM Release 1.93 manual A Programmable Optimizing Electro-Magnetic Simulator Release, August 21 Phil Hobbs, Poems

Contents

?88,.d88b, d8888b d8888b January 22 Updated January 23 Phil HobbsA Programmable Optimizing Electro-Magnetic Simulator Release IBM T. J. Watson Research Center Yorktown Heights, NYPage POEMS A Programmable Optimizing Electro-Magnetic Simulator ReleaseAugust 21 Phil Hobbs IBM T. J. Watson Research Center Yorktown Heights, NYAppendix A. V-Antenna Optimization Run Chapter1.1. Motivation Chapter Introduction1.1.1. Philosophy 1.1.3. Optimization 1.1.2. StructurePage Postprocessor Optimization LoopFront-End Script POEMS.CMD FDTD Engine FIDO + CLUSTER SCRIPT or TEMPEST2.2.1. Script Operation 2. HOW POEMS WORKS2.1. Program Organization 2.2. The Front-End Script poems.cmd2.4. The Postprocessor EMPOST 2.3. The FDTD Engine FIDO/TEMPEST2.5. The Visualization System VIS5D 2.6. Cluster Controlempost ordersfile resultsfile 2.6.1. Parallel Processing 3.1.2. GLOBAL Group 3. USING POEMS3.1. Command Reference 3.1.1. poems Command-Line OptionsFREQ 200*THz which means that the host’s predefined hostname is not used NB Hosts are identified by their hostname as specified by the$HOSTNAME environment variable, rather than by IP address. This aliases in different SUBDOMAIN statements in the parameters fileMACRO Syntax PRINT string expression 3.1.3. WORLD Group Syntax XRANGE minexpression maxexpression 3.1.4. MATERIAL GroupParameters epsReal epsImag muReal muImag DEFINE matname= BlackGlass type=dielectric n=1.52 k=1e-3DEFINE matname=copper type=conductor conductivity=5.8e7 xhi1,yhi1,zhi1 and xlo2,ylo2,zlo2, xhi2,yhi2,zhi2. These must be 3.1.5. OBJECT GroupParameters mattype xlo xhi ylo yhi zlo zhi thickness phasefunc intervalxhi2 ylo2 yhi2 zlo2 zhi2 Parameters matname orientation xlo xhi ylo yhi zlo zhi widthfuncParameters mattype axis lo hi maxradius insidefunc Parameters mattype xlo xhi ylo yhi zlo zhiParameters mattype slo shi perp hfunc vfunc Not implemented in this releaseParametersNot implemented in this release 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 3.1.6. SOURCE GroupParameters reltolerance maxcycles mincycles timestep 3.1.7. COMMAND GroupParameters x y z Ex Ey Ez mag phase 3.1.8. OUTPUT Groupphase state 3.1.9. POSTPROCESS GroupParameters variable xlo xhi ylo yhi zlo zhi x y z file Parameters name xinside yinside zinside xlo xhi ylo yhi zlo zhi supplied plane, xmin ≤ x ≤ xmax, ymin ≤ y≤ ymax, zmin ≤ z ≤ zmaxParameters name xlo xhi ylo yhi zlo zhi file direction towards from xInside, yInside, zInsideThe value of coordinates can be polar or rectangular, and arrangement Parameters file orientation xlo xhi ylo yhi zlo zhi phase indexnMOVIE Parametersfile orientation xlo xhi ylo yhi zlo zhi dt frames variableParameters frames variable1 variable2 variable3 variable4 variable5 variable6 variable7 variable8 variable9 xlo xhi ylo yhi zlo zhi filelevel Parameters name xlo xhi ylo yhi zlo zhi file xinside yinside zinsidelevel DISSIPATION3.1.10. OPTIMIZE Group not be updated at postprocessing time. That means that the Syntax GUESS var1 = expression var2 =, expressionarbitrary boundary. Not implemented in this release Stepping is not implemented in this release 3.1.11. SCHEDULE GroupPENALTY expression Parameters iterations tolerance restartsActual Domain N rows Image Domain N-2 Rows Next Period 3.2. The Computational Domain3.2.1. Symmetry Unfolded 2N-2 Rows3.3. OBJECTS 3.5.1. POINT SOURCES 3.3.1. Perfectly-Matched Layers3.4. MATERIALS 3.5. SOURCES3.5.2. PLANE WAVES Page 3.5.4. MODEFILE SOURCES 3.5.3. BEAM SOURCES3.6.1. Merit Functions 3.6. Optimizationf. Reflection from an interface c. Phase uniformity across a planee. Power dissipation in a load dissipation d. RMS Voltage across terminal pointsF i g u r e 2 . 1 2 Optimized V-antenna 3.6.2. Worked Example Optimizing a V AntennaPage 3.6.3. Worked Example Doped Silica Waveguide Mode 3.6.4. Worked Example Glass Ridge Waveguide to Free Space Coupler3.7. Predefined Constants i l l u m i n a t i o n p l a n e a f t e r 3.7.1. Reserved Names3.8. Predefined Mathematical Functions 3.8.1. Arithmetic Operatorsa. Differences from C/C++ 3.8.5. ACOSH 3.8.2. Logical Operators3.8.3. ABS 3.8.4. ACOS3.8.13. ELINTK 3.8.10. ATAN23.8.11. CEIL 3.8.12. COS3.8.22. MIN 3.8.19. INTEGRAL3.8.20. LN 3.8.21. MAX3.8.29. SIGN 3.8.26. RANDOM3.8.27. ROOT1D 3.8.28. ROUND3.10. Material Parameter Functions 3.9. Analytical Pupil FunctionsA. TEMPEST and General FDTD Information 1. Startup and Steady State2. Time step fprintfstderr,periodstep %d\n,periodstep printftime step after adjustment %e\n,dtfflushstdoutperiodstep+=2 i=periodstep%4 if i!=0 periodstep=periodstep+4-i dt=2.0*PI/omega*FLOATperiodstepThe run will produce lots of bitmaps at appropriately chosen locations Appendix A. V-Antenna Optimization RunA.1. POEMS INPUT DIPOLE2I.PAR set kx0=0 /* direction of plane wave set loadlength=2*step Fundamental dimensionsset phi=0 /* phase of excitation */ set ky0=step/lambda set dipoleheight=0.2*um set dipolethickness=0.30*umEND WORLD SUBDOMAIN ALL XRANGE 0 xsize YRANGE 0 ysize ZRANGE 0 zsize DEFINE matname= Loadtype=dielectric n=nload k=kloadFAN matname=PCL taper=linear taperpar=0 END COMMAND COMPUTE maxcycles = 25 reltolerance = END OUTPUTzhi=zsize phase=0 state=steady file=Dipole2Iexi zlo2 = Zvertex-DipoleThickness/2 zhi2= Zvertex+DipoleThickness/2STORE Dipole2I.simplex MERIT diss*1e12 END OPTIMIZEAMPLEX PHASEEXxlo=0 xhi=xsize ylo=0 yhi=ysize zlo=0 zhi=zsize file=Dipole2IPhaseEx POSTPROCESSFLUX name=MiddleFlux , /* crossing the middle z=zsize/2 xLo = Tpml zlo=zsize/2+step zhi=zsize/2+step, file=Dipole2IIndexXY.bmp phase=0file=Dipole2IPZXY.bmp phase=0 palette=redblue zhi = zsize-Tpml-step xinside=xsize/2 yinside=ysize/2 zinside=zsize/2file=Dipole2IAmpExZX.bmp phase=0 palette=flame file=Dipole2IPZZX.bmp phase=0 palette=redbluefile=Dipole2IPXXY.bmp phase=0 palette=redblue file=Dipole2IPXZX.bmp phase=0 palette=redbluefile=Dipole2IDissXY.bmp palette=redblue SLICE orientation=xy variable=PhaseEx , xlo=0 xhi=xsizeSLICE orientation=zx variable=PhaseEx , xlo=0 xhi=xsize file=Dipole2IExIZX.bmp phase=0 palette=redblueylo=ysize/2 yhi=ysize/2, zlo=0 zhi=zsize file=Dipole2IExQZX.bmp phase=pi/2 palette=redblue ENDplot hy steady 0.250000 node 0 99 0 99 0 99 dipole2ihyq A.2. tempest Input File DIPOLE2I.PAR.INrectangle node 0 100 index 1.000000Written by 56 pages of pointsource statements omitted empost dipole2i.orders outfile v5dfile A.3. Postprocessor orders DIPOLE2I.ORDERSDIPOLE2IEYI NULLEXI ARRAY Domain for comparison e.g. modematch 0.00000 SLICE POSTPROC.16.COMPARISONDOMAIN 0.00000 SLICE A.4. Run Results DIPOLE2I.SIMPLEX 0.000000883423566665 77.7563866665 0.0000017351705 Penalty = Contracting by 0.5 in 1D 1-D shrink unsuccessful0.00000064767135 98.83191 0.000001519205 Penalty = 0.00000069266248611 90.849532778 0.00000152214633334 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 1DFailed to converge in 36 function evaluations 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 1D0.000000572297880735 102.392565535 0.00000148853488604 Penalty = Reflecting through centroid. Contracting by 0.5 in 1DReflecting through centroid. Contracting by 0.5 in 1D 1-D shrink unsuccessfulReflecting through centroid. Contracting by 0.5 in 1D Reflecting through centroid. Contracting by 0.5 in 1DConverged in 56 function evaluations Appendix B. Configuration B.1. FDTD and TEMPESTB.1.1. tempest patches c. Hints and Kinks B.1.3. Advice common to all or most FDTD programsB.1.2. tempest limitations b. Understanding convergenceB.3. X Window System Configuration B.3.1. Sample X11 ConfigurationB.2. REXX B.3.2. Running Vis5D B.4. Release NotesB.4.1. Wish list B.4.2. Beta Release Limitations XY plane Z perpendicular +X right, +Y up, +Z out of screen YZ plane X perpendicular +Y right, +Z up, +X out of screenZX plane Y perpendicular +Z right, +X up, +Y out of screen Appendix C. Multiprocessing C++ 7 INDEXAiry pattern Antenna 1, 39, 40, 45, 51, 52 BASICSTEP 16, 18, 24, 54, 62 Boundary 7, 15, 19, 29, 31Periodic boundary conditions 31, 34