Manuals / IBM / Communications / PDAs & Smartphones

IBM Release 1.93 manual DEFINE matname= Load, type=dielectric n=nload, k=kload

Models: Release 1.93

1 93

Download 93 pages 50.46 Kb

Page 58

Image 58

DEFINE matname= Load

type=dielectric n=n_load

k=k_load

END

WORLD

SUBDOMAIN ALL

XRANGE 0 xsize

YRANGE 0 ysize

ZRANGE 0 zsize

BASICSTEP step

BOUNDARY ymax illum

BOUNDARY xmin illum

BOUNDARY zmax illum

END

OBJECT

/* Absorber around the outside of the domain */

HOLLOWBOX matname=AirPML xlo=0 xhi=xsize ylo=0 yhi=ysize zlo=0 zhi=zsize thickness=Tpml

/* Air filling all space inside PML */

BLOCK matname=Air xlo=Tpml xhi=xsize-Tpml ylo=Tpml yhi=ysize-Tpml , zlo=Tpml zhi=zsize-Tpml

/* Dipole

FAN matname=PCL taper=linear taperpar=0,

xlo1 = (xsize-DipoleLengthUM*um)/2 xhi1=(xsize-DipoleLengthUM*um)/2, ylo1 = (ysize-DipoleHeight)/2 yhi1=(ysize+DipoleHeight)/2,

zlo1 = Zvertex-ZRakeUM*um-DipoleThickness/2 zhi1= Zvertex-ZRakeUM*um+DipoleThickness/2, xlo2 = (xsize-LoadLength)/2 xhi2=(xsize-LoadLength)/2,

ylo2 = (ysize-DipoleHeight)/2 yhi2=(ysize+DipoleHeight)/2,

Page 58

Image 58

IBM Release 1.93 manual DEFINE matname= Load, type=dielectric n=nload, k=kload, BOUNDARY zmax illum END OBJECT

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 EMPOST2.6. Cluster Control 2.5. The Visualization System VIS5Dempost ordersfile resultsfile 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= BlackGlass type=dielectric n=1.52 k=1e-3 Parameters epsReal epsImag muReal muImagDEFINE matname=copper type=conductor conductivity=5.8e7 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 timestep3.1.9. POSTPROCESS Group phase stateParameters variable xlo xhi ylo yhi zlo zhi x y z file 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 Syntax GUESS var1 = expression var2 =, expression not be updated at postprocessing time. That means that thearbitrary boundary. Not implemented in this release 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.6.4. Worked Example Glass Ridge Waveguide to Free Space Coupler 3.6.3. Worked Example Doped Silica Waveguide Mode3.7. Predefined Constants 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 r3.8.1. Arithmetic Operators 3.8. Predefined Mathematical Functionsa. Differences from C/C++ 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 Functions1. Startup and Steady State A. TEMPEST and General FDTD Information2. Time step 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 =Reflecting through centroid. Contracting by 0.5 in 1D Reflecting through centroid. Contracting by 0.5 in 1DConverged in 56 function evaluations B.1. FDTD and TEMPEST Appendix B. ConfigurationB.1.1. tempest patches b. Understanding convergence B.1.3. Advice common to all or most FDTD programsB.1.2. tempest limitations c. Hints and KinksB.3.1. Sample X11 Configuration B.3. X Window System ConfigurationB.2. REXX B.4. Release Notes B.3.2. Running Vis5DB.4.1. Wish list B.4.2. Beta Release Limitations YZ plane X perpendicular +Y right, +Z up, +X out of screen XY plane Z perpendicular +X right, +Y up, +Z out of screenZX plane Y perpendicular +Z right, +X up, +Y 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