Manuals / IBM / Communications / PDAs & Smartphones

IBM Release 1.93 manual A.3. Postprocessor orders DIPOLE2I.ORDERS

Models: Release 1.93

1 93

Download 93 pages 50.46 Kb

Page 69

...
err_tol 2.000000E-04
max_cycle 8
plot refractive	node 0 97 0 97 0 119				dipole2i_refractive
plot ex steady	0.000000	node 0 97 0			97 0	119	dipole2i_exi
plot ex steady	0.250000	node 0 97 0			97 0	119	dipole2i_exq
plot ey steady	0.000000	node 0 97 0			97 0	119	dipole2i_eyi
plot ey steady	0.250000	node 0 97 0			97 0	119	dipole2i_eyq
plot ez steady	0.000000	node 0 97 0			97 0	119	dipole2i_ezi
plot ez steady	0.250000	node 0 97 0			97 0	119	dipole2i_ezq
plot hx steady	0.000000	node 0 97 0			97 0	119	dipole2i_hxi
plot hx steady	0.250000	node 0	97	0	97 0	119	dipole2i_hxq
plot hy steady	0.000000	node 0	97	0	97 0	119	dipole2i_hyi
plot hy steady	0.250000	node 0	97	0	97 0	119	dipole2i_hyq
plot hz steady	0.000000	node 0	97	0	97 0	119	dipole2i_hzi

A.3. Postprocessor orders: DIPOLE2I.ORDERS

This file is machine-generated. It is the input for EMPOST, and is regenerated on each iteration of the optimizer. Sometimes it’s useful to hack this file up a little, e.g. to change a bitmap’s scale linear to logarithmic scale, or to change a colour palette. To do this, make the change and call empost directly from the command line:

empost dipole2i.orders outfile v5dfile

Files outfile and v5dfile are, respectively, the file containing the calculated scalar results (OUTPUTFACETFLUX, INPUTFACETFLUX, MIDDLEFLUX, and DISS), and a dummy filename for when the Vis5D support is integrated into empost.

Alternatively, you can change the input file (here dipole2i.par), run poems.cmd to generate a new orders file, stop the simulation by pressing ctrl-break, and then run empost manually as above.

/* dipole2i: */
/* Optimization of a V antenna for maximum power absorption from free space */
/* Parameters optimized: length and rake of		antenna, load impedance */
Dipole2I	/*	Run title string (no whitespace)
	/*	World symmetry (X, Y, Z, XY, XZ, YZ, XYZ)
1.58000000E-006	/*	Free-space wavelength (metres)
5	/*	Number of subdomains
ALL	/*	Subdomain name (no whitespace)
0	/*	Type: Computational partition (0) or comparison (1)(e.g. for modematch)
DIPOLE2I_EXI	/*	ExI: file name for the real	part of Ex
DIPOLE2I_EXQ	/*	ExQ: file name for the imaginary part of Ex

65

Image 69

IBM Release 1.93 manual A.3. Postprocessor orders DIPOLE2I.ORDERS, empost dipole2i.orders outfile v5dfile

Contents

A Programmable Optimizing Electro-Magnetic Simulator Release January 22 Updated January 23 Phil HobbsIBM T. J. Watson Research Center Yorktown Heights, NY ?88,.d88b, d8888b d8888bPage August 21 Phil Hobbs A Programmable Optimizing Electro-Magnetic Simulator ReleaseIBM T. J. Watson Research Center Yorktown Heights, NY POEMSAppendix A. V-Antenna Optimization Run Chapter1.1. Motivation Chapter Introduction1.1.1. Philosophy 1.1.3. Optimization 1.1.2. StructurePage Front-End Script POEMS.CMD Optimization LoopFDTD Engine FIDO + CLUSTER SCRIPT or TEMPEST Postprocessor2.1. Program Organization 2. HOW POEMS WORKS2.2. The Front-End Script poems.cmd 2.2.1. Script Operation2.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. Command Reference 3. USING POEMS3.1.1. poems Command-Line Options 3.1.2. GLOBAL GroupFREQ 200*THz $HOSTNAME environment variable, rather than by IP address. This NB Hosts are identified by their hostname as specified by thealiases in different SUBDOMAIN statements in the parameters file which means that the host’s predefined hostname is not usedMACRO 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 Groupxhi2 ylo2 yhi2 zlo2 zhi2 phasefunc intervalParameters matname orientation xlo xhi ylo yhi zlo zhi widthfunc Parameters mattype xlo xhi ylo yhi zlo zhi thicknessParameters mattype slo shi perp hfunc vfunc Parameters mattype xlo xhi ylo yhi zlo zhiNot implemented in this release Parameters mattype axis lo hi maxradius insidefuncParameters x y z kx ky kz Ex Ey Ez mag phase orientation Parameters x y z width kx ky kz Ex Ey Ez mag phase n k orientation3.1.6. SOURCE Group ParametersNot implemented in this releaseParameters x y z Ex Ey Ez mag phase 3.1.7. COMMAND Group3.1.8. OUTPUT Group Parameters reltolerance maxcycles mincycles timestepphase state 3.1.9. POSTPROCESS GroupParameters variable xlo xhi ylo yhi zlo zhi x y z file Parameters name xlo xhi ylo yhi zlo zhi file direction supplied plane, xmin ≤ x ≤ xmax, ymin ≤ y≤ ymax, zmin ≤ z ≤ zmaxtowards from xInside, yInside, zInside Parameters name xinside yinside zinside xlo xhi ylo yhi zlo zhiThe value of coordinates can be polar or rectangular, and arrangement Parameters file orientation xlo xhi ylo yhi zlo zhi phase indexnParameters frames variable1 variable2 variable3 variable4 variable5 Parametersfile orientation xlo xhi ylo yhi zlo zhi dt frames variablevariable6 variable7 variable8 variable9 xlo xhi ylo yhi zlo zhi file MOVIElevel Parameters name xlo xhi ylo yhi zlo zhi file xinside yinside zinsideDISSIPATION level3.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 PENALTY expression 3.1.11. SCHEDULE GroupParameters iterations tolerance restarts Stepping is not implemented in this release3.2.1. Symmetry 3.2. The Computational DomainUnfolded 2N-2 Rows Actual Domain N rows Image Domain N-2 Rows Next Period3.3. OBJECTS 3.4. MATERIALS 3.3.1. Perfectly-Matched Layers3.5. SOURCES 3.5.1. POINT SOURCES3.5.2. PLANE WAVES Page 3.5.4. MODEFILE SOURCES 3.5.3. BEAM SOURCES3.6.1. Merit Functions 3.6. Optimizatione. Power dissipation in a load dissipation c. Phase uniformity across a planed. RMS Voltage across terminal points f. Reflection from an interfaceF 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.3. ABS 3.8.2. Logical Operators3.8.4. ACOS 3.8.5. ACOSH3.8.11. CEIL 3.8.10. ATAN23.8.12. COS 3.8.13. ELINTK3.8.20. LN 3.8.19. INTEGRAL3.8.21. MAX 3.8.22. MIN3.8.27. ROOT1D 3.8.26. RANDOM3.8.28. ROUND 3.8.29. SIGN3.10. Material Parameter Functions 3.9. Analytical Pupil FunctionsA. TEMPEST and General FDTD Information 1. Startup and Steady State2. Time step periodstep+=2 i=periodstep%4 if i!=0 periodstep=periodstep+4-i printftime step after adjustment %e\n,dtfflushstdoutdt=2.0*PI/omega*FLOATperiodstep fprintfstderr,periodstep %d\n,periodstepThe run will produce lots of bitmaps at appropriately chosen locations Appendix A. V-Antenna Optimization RunA.1. POEMS INPUT DIPOLE2I.PAR set phi=0 /* phase of excitation */ set ky0=step/lambda Fundamental dimensionsset dipoleheight=0.2*um set dipolethickness=0.30*um set kx0=0 /* direction of plane wave set loadlength=2*steptype=dielectric n=nload DEFINE matname= Loadk=kload END WORLD SUBDOMAIN ALL XRANGE 0 xsize YRANGE 0 ysize ZRANGE 0 zsizezhi=zsize phase=0 state=steady file=Dipole2Iexi END COMMAND COMPUTE maxcycles = 25 reltolerance = END OUTPUTzlo2 = Zvertex-DipoleThickness/2 zhi2= Zvertex+DipoleThickness/2 FAN matname=PCL taper=linear taperpar=0STORE Dipole2I.simplex MERIT diss*1e12 END OPTIMIZExlo=0 xhi=xsize ylo=0 yhi=ysize zlo=0 zhi=zsize file=Dipole2IPhaseEx PHASEEXPOSTPROCESS AMPLEXfile=Dipole2IPZXY.bmp phase=0 palette=redblue zlo=zsize/2+step zhi=zsize/2+step, file=Dipole2IIndexXY.bmp phase=0zhi = zsize-Tpml-step xinside=xsize/2 yinside=ysize/2 zinside=zsize/2 FLUX name=MiddleFlux , /* crossing the middle z=zsize/2 xLo = Tpmlfile=Dipole2IPXXY.bmp phase=0 palette=redblue file=Dipole2IPZZX.bmp phase=0 palette=redbluefile=Dipole2IPXZX.bmp phase=0 palette=redblue file=Dipole2IAmpExZX.bmp phase=0 palette=flameSLICE orientation=zx variable=PhaseEx , xlo=0 xhi=xsize SLICE orientation=xy variable=PhaseEx , xlo=0 xhi=xsizefile=Dipole2IExIZX.bmp phase=0 palette=redblue file=Dipole2IDissXY.bmp palette=redblueylo=ysize/2 yhi=ysize/2, zlo=0 zhi=zsize file=Dipole2IExQZX.bmp phase=pi/2 palette=redblue ENDrectangle node A.2. tempest Input File DIPOLE2I.PAR.IN0 100 index 1.000000 plot hy steady 0.250000 node 0 99 0 99 0 99 dipole2ihyqWritten 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.00000064767135 98.83191 0.000001519205 Penalty = Contracting by 0.5 in 1D 1-D shrink unsuccessful0.00000069266248611 90.849532778 0.00000152214633334 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 evaluationsReflecting through centroid. Contracting by 0.5 in 1D Reflecting through centroid. Contracting by 0.5 in 1D1-D shrink unsuccessful 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 Appendix B. Configuration B.1. FDTD and TEMPESTB.1.1. tempest patches B.1.2. tempest limitations B.1.3. Advice common to all or most FDTD programsb. Understanding convergence c. Hints and KinksB.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 Airy pattern Antenna 1, 39, 40, 45, 51, 52 INDEXBASICSTEP 16, 18, 24, 54, 62 Boundary 7, 15, 19, 29, 31 C++ 7Periodic boundary conditions 31, 34