036307-6-T USER MANUAL FOR MR-TETRA.F: A MULTIRESOLUTION FINITE ELEMENT-BOUNDARY INTEGRAL (FE-BI) CODE WITH HIERARCHICAL TETRAHEDRAL ELEMENTS Lars S. Andersen John L. Volakis Mission Research Corp. 147 John Sims Parkway Valparaiso, FL. 32580 Air Force Office of Scientific Research AFOSR/NM Boiling AFB, DC 20332-0001 June 1999 36307-6-T = RL-2507

PROJECT INFORMATION PROJECT TITLE: REPORT TITLE: Focused Application Software for Design of Ferrite Patch Antennas USER MANUAL FOR MR-TETRA.F: A MULTIRESOLUTION FINITE ELEMENT-BOUNDARY INTEGRAL (FE-BI) CODE WITH HIERARCHICAL TETRAHEDRAL ELEMENTS U-M REPORT No.: 036307-6-T CONTRACT START DATE: END DATE: July 1997 June 1999 June 1999 DATE: SPONSOR: S. McCool and K. Trott Mission Research Corp. 147 John Sims Parkway Valparaiso, FL. 32580 trott@ fwb.gulf.net Phone: 904-729-4494 Fax: 904-729-2443 Dr. Arje Nachman Air Force Office of Scientific Research AFOSR/NM Bolling AFB, DC 20332-0001 nachman@afosr.af.mil SPONSOR CONTRACT No.: U-M PRINCIPAL INVESTIGATOR: SC-0022-97-001 (Air Force Prime Contract no. F49620-97-C-0022) John L. Volakis EECS Dept. University of Michigan 1301 Beal Ave Ann Arbor, MI 48109-2122 Phone: (313) 764-0500 FAX: (313) 747-2106 volakis @umich.edu http://www-personal.engin.umich.edu/-volakis/ CONTRIBUTORS TO THIS REPORT: Lars S. Andersen and John L. Volakis

User manual for MR-TETRA.f: A multi-resolution finite element - boundary integral (FE-BI) code with hierarchical tetrahedral elements Lars S. Andersen and John L. Volakis Radiation Laboratory Department of Electrical Engineering and Computer Science University of Michigan Ann Arbor, MI 48109-2122, USA Abstract The following is a manual for a multi-resolution finite element / boundary integral code for analysis of printed antennas backed by material-filled (dielectric or ferrite) metallic cavities recessed in infinite metallic ground planes. The main characteristic of the code is several attractive options for expanding the electric field within the cavity using lowest order, higher order or any combination of lowest and higher order hierarchical mixed-order tangential vector finite elements for tetrahedra. 1 Introduction The following is a manual for a multi-resolution finite element / boundary integral (FE/BI) code (referred to as MR-TETRA in the following) for analysis of printed metallic antennas backed by material-filled metallic cavities recessed in infinite metallic ground planes. It is based on a standard FE/BI formulation where the cavity volume is discretized into tetrahedral elements using a FE approach and the mesh is rigorously truncated at the aperture surface using a BI [1]. The electric field within each tetrahedral element is expanded using hierarchical mixed-order tangential vector finite elements (TVFEs) of order 0.5 or 1.5 (referred to as lowest and higher order TVFEs in the following) as presented by Andersen and Volakis [2]. The resulting linear equation system is solved using an iterative solver. The matrix-vector products within the iterative solver can be carried out either directly or via a two-dimensional discrete Fourier transform (DFT) in the spectral domain through application of the adaptive integral method (AIM) [3]. A quasi minimal residual (QMR) [4] or conjugate gradient squared (CGS) [4] solver can be used in case matrix-vector products are computed directly while a QMR or biconjugate gradient (BCG) [4] solver can be used in 1

case matrix-vector products are computed using a DFT. Upon solution of the resulting linear equation system, input impedance and, optionally, near and far fields (patterns, polarization characteristics) are computed. The purpose of this manual is to familiarize the reader with the steps involved in running MR-TETRA. This is achieved via a general description as well as a set of examples. It is not the purpose of this manual to demonstrate the merits of the approach on which the code is based. For such a demonstration, see for instance [5]. This manual is organized as follows. Section 2 provides a general description of how to run MR-TETRA. All relevant steps needed by the user will be presented. Section 3 outlines the construction of the input file for MR-TETRA. Section 4 offers a few examples of how to run MR-TETRA. Surface meshes, input files as well as results are given in order to demonstrate various capabilities of the code. Section 5 summarizes the manual. 2 General description of how to run MR-TETRA The root directory of MR-TETRA consists of two sub-directories. The sub-directory src contains all source files (.f), parameter files (.h) and a makefile MMR-TETRA. Upon compilation, it will also contain all object files (.o). The executable MR-TETRA will be in the sub-directory bin where the input file MR-TETRA. in and the universal file describing the mesh (see below) must also be placed. Upon execution, this sub-directory will also contain all output files. The sub-directory bin further contains a sub-directory plot with MATLAB codes for plotting various output files. These MATLAB codes are not documented in this manual. To run the code, the following steps must be followed: 1. A mesh in universal file format must be provided in the sub-directory bin. 2. In the sub-directory src, several dimensioning parameters in the file paraml.h must be initialized. They must be estimated conservatively, otherwise the code will either crash or produce erroneous results. After a first run for a given mesh, they can manually be set to their smallest possible values for that particular mesh. The parameters are Nnomax (max number of nodes), Nsnomax (max number of metallic nodes), Nedmax (max number of edges), Nsedmax (max number of metallic edges), Nfamax (max number of triangular faces), Nsfamax (max number of metallic triangular faces), Nmgmax (max number of material groups), Npredmax (max number of probe edges), Nhovolmax (max number of higher order sections), Nbifamax (max number of triangular BI faces), Nbiedmax (max number of BI edges), IAIMmax (max number of AIM grid points in the x-direction), JAIMmax (max number of AIM grid points in the y-direction), IJAIMmax (max of IAIMmax and JAIMmax) and Nnffamax (max number of near field triangular AIM faces). Note that the file parami.h contains several additional dimensioning parameters. These have been set conservatively so the code will run for fully higher order cavities without AIM. This, however, is a tremendous overkill if lowest order 2

TVFEs and AIM is applied. The user is highly advised to familiarize himself / herself with these parameters and set them optimally for each application. Specifically, BWmaxl and Nunkmax can be lowered from 100 and 2*Nedmax+2*Nfamax to 30 and Nedmax if only lowest order TVFEs are applied and Nmatmax can be lowered from BWmaxl*Nunkmax+Nbiedmax**2 to BWmaxl*Nunkmax if AIM is used. This will lead to significant memory savings without altering the results. 3. Build the executable MR-TETRA (type make -f MMR-TETRA in the sub-directory src 1). 4. Construct the input file MR-TETRA. in in the sub-directory bin (see below). 5. Run the code (type MR-TETRA in the sub-directory bin). 3 Construction of input file The input file for MR-TETRA can be broken into 12 different sections. The format of each of these 12 sections will be described below. "I" denotes an integer, "R" denotes a double precision real number and "C" denotes a double precision complex number. "Text" denotes a text line that allows the user to describe the following input parameter(s). Information in all such lines is irrelevant to the code but it makes the input file easier to read and hereby easier to modify. We note that some parts seem unnecessary or redundant. This is partly due to the fact that the code has more options than those described in this manual and partly a result of convenient choices made while the code was under development. We also note that the following description is very general and that specific examples follow later in the report. Section 1 This section describes the universal file that represents the mesh. Users unfamiliar with universal file formats are referred to SDRC I-DEAS manuals. A universal file section with the descriptor 2411 must give the node coordinates in cm (the ground plane and the metallic antenna must be in the plane z = 0 and the cavity must be in the half space z < 0). A universal file section with the descriptor 2412 must give the element connectivity and a material group for each element. Optionally, a universal file section with descriptor 2417 (used to describe groups of nodes) and name starting with "F" can be added. This section is used for describing the metallic triangular faces forming the printed antenna. Each node triplet describes a metallic face and hence the group will contain three times as many nodes as there are metallic faces. Note that this is the only way to uniquely describe the metallic antenna. For simple geometries, it can be done with nodes alone or with node doublets (edges) but in the general case node triplets (triangular faces) must be used. MR-TETRA allows easy specification of rectangular and circular patches via simple geometrical parameters (see Section 2). In these cases, a universal file section with descriptor 2417 is unnecessary. Note that the universal file can be generated using SDRC I-DEAS, the automatic mesher written at 1We note that the compilation of source files was tested on a SUN ULTRA30 work station. Certain parts of the makefile MMR-TETRA as well as the timing commands in the code may have to be modified on other platforms. 3

the University of Michigan or any other meshing package capable of producing a tetrahedral mesh described in universal file format. The format of this section is the following: Text "filename": Universal file Section 2 This section describes the metallic cavity and, optionally, the metallic patch. A cavity code is 1 / 2 for a rectangular / circular cavity. In the former case, (x, y, Z)start and (x, y, z)stop for the rectangular cavity as well as (x, y)start and (x, y)stop for a rectangular patch in the plane z = 0 is given. In the latter case, the radius rcavity and height hcavity of a circular cavity centered at (0, 0, -hcavity/2) as well as the radius rpatch of a circular patch centered at (0, 0, 0) is given. For cavities where the metallic faces are described in the universal file, all patch parameters should be set to zero. The format of this section is either of the following Text 1: Cavity code (rectangular cavity) Text R: xcavitystart R: ycavitystart R: zcavitystart R: xcavity.stop R: ycavitystop R: zcavitystop R: xpatchstart R: y-patch.start R: xpatchstop R: ypatchstop Text 2: Cavity code (circular cavity) Text R: r.cavity R: hcavity R: rpatch Section 3 This section describes the probe excitation. Each probe has unit magnitude and zero phase and must start and end at a node in the mesh. A probe code (must be 1) and the number of probes is given along with (x, y, Z)start and (x, y, z)stop for each probe. The format of this section is the following: 4

Text 1: Probe code Text I: Number of probes Text R: x-probe-start R: y-probe-start R: zprobestart Repeated for each probe R: x-probestop R: y-probe-stop R: z-probe-stop Section 4 This section describes the higher order sections in the mesh. A higher order code of 0 / 1 / 2 indicates that no higher order TVFEs are used / that higher order TVFEs are contained within rectangular brick sections of the mesh (for rectangular cavities) / that higher order TVFEs are contained within cylindrical shell sections of the mesh (for circular cavities). In case higher order TVFEs are used, the number of higher order sections is given along with a geometrical description of each higher order section. For rectangular brick sections, (X, y, Z)start and (x, y, Z)stop are given. The higher order section is then the brick having these points as opposite corners. For cylindrical shell sections, the inner and outer radii rstart and rstop and the azimuthal angles <start and tstop are given. Ostart and cs4top must be in degrees and be in the interval [0~, 360~). The higher order section is then the shell from rstart to Vstop radially, from (start to Ustop (counter-clockwise) azimuthally and from -hcavity to 0 in the z-direction. The format of this section is either of the following Text 0: Higher order code (no higher order TVFEs) Text 1 Higher order code (rectangular cavity) Text I: Number of higher order sections Text R: y_ho_start R z_ho_start Repeated for each higher order section R: x-ho-stop R: y-ho-stop R: z_ho_stop 5

Text 2: Higher order code (circular cavity) Text I: Number of higher order sections Text R: rho-start R rhostop > Repeated for each higher order section R: phi-ho-start R: phi-ho_stop Section 5 This section describes the TVFEs used for field expansion. A TVFE code 1 / 3 indicates that lowest order TVFEs / a combination of lowest and higher order TVFEs are used. Note that for a TVFE code of 1, the information about higher order TVFEs in Section 4 is not used. This allows the user to toggle between lowest order TVFEs and a combination of lowest and higher order TVFEs by changing only the TVFE code. The format of this section is the following: Text I: TVFE code Section 6 This section describes the different material groups. The number of material groups is given. For each material group, the relative permittivity and permeability tensor is described. The former is described by the nine elements e,,..., z and a conductivity r in S/cm that gives rise to a frequency dependent loss. The latter is either the unit tensor or a Polder tensor for a gyrotropic material. It is described by a bias code (0 means a unit tensor while 1 / 2 / 3 means a Polder tensor with x- / y- / z-bias) and the value of the precession and saturation frequencies fo and fm in Hz [6]. The format of this section is the following: Text I: Number of material groups Text C: epsilon-xx C: epsilon-xy C: epsilon-xz C: epsilon-yx C: epsilon-yy C epsilon.. yz C epsilon-yz Repeated for each material group C: epsilon-zx C: epsilon-zy C: epsilon-zz R: sigma I: Bias axis R: f_0 R: fm 6

Section 7 This section describes the orders of the Gauss-Legendre integrations for computing the FEM matrix elements, BI near zone elements, BI far zone elements and far fields. The first type of integration is a volume integration over a tetrahedron where an integer 1 / 2 / 3 / 4 denotes 1- / 4- / 5- / 11-point Gauss-Legendre integration. The last three types of integrations are surface integrations over triangles where an integer 1 / 2 / 3 / 4 denotes 1- / 3- / 4- / 7-point Gauss-Legendre integration. The integration order for FEM should be at least 2 for lowest order TVFEs and 4 for higher order TVFEs if exact results are desired. Since accurate evaluation of the near zone elements in the BI matrix is crucial for accurate analysis, the integration order for the BI near zone elements should be set to 4. The integration order for the far zone elements in the BI matrix is usually set to 1 but can be increased to give slightly more accurate results at the price of a significantly slower code. The integration order for far fields is usually set to 4. The format of this section is the following: Text I: Integration order for FEM I: Integration order for BI_near I: Integration order for BIfar I: Integration order for far field Section 8 This section describes the frequencies at which the geometry is analyzed. The start and stop frequency in Hz as well as the number of frequency points is given. The format of this section is the following: Text R: f-start R: fstop Text I: Number of frequency points Section 9 This section describes the iterative solver used for solving the resulting matrix equation system as well as the tolerance imposed on the relative residual as a stopping criterion. A solver code of 1 indicates a QMR solver. A solver code of 2 indicates a CGS solver for the case where matrix-vector products are computed directly and a BCG solver for the case where AIM is applied and matrix-vector products are computed using a DFT. The format of this section is the following: Text I: Solver code Text R: Solver tolerance 7

Section 10 This section describes whether AIM is used as well as all AIM parameters. An AIM code 0 / 1 indicates that AIM is not / is used. Regardless of the AIM code, (x, y)start and (x, y)stop for the AIM grid are given along with the number of AIM grid points in the x- and y-direction and the number of near zone AIM grid points in both directions. This information is not used when AIM is not used. Reading the information anyway allows the user to toggle between direct and DFT computation of matrix-vector products by changing only the AIM code. The format of this section is the following Text I: AIM code Text R x AIM start R y-AIM-start R: x-AIM-stop R: y-AIMstop I: N_AIM_x I: N-AIM-y I: NAIM-near Section 11 This section describes whether near and far field distributions are computed. For near fields, a near field code of 0 / 1 indicates that they are not / are computed. For far fields, a far field code of 0 / 1 / 2 / 3 indicates that nothing / patterns only / polarization characteristics only / patterns and polarization characteristics are computed. For near fields, the cavity is sampled in 33 x 33 points in 4 equidistant planes from the bottom of the cavity to the top, as illustrated in Fig. 1. The file xpyp.dat has 33. 33 = 1089 lines containing the values of (x, y) for the sample points in the plane z = 0. The file FxFy.dat has 33. 33 = 1089 lines containing the corresponding values of the real part of the x- and y-directed electric field as well as the imaginary part of the x- and y-directed electric field. This surface electric field is the one that is being radiated (as an equivalent magnetic current) to give the far field. The file xpypzp.dat has 33. 33 - 4 = 4356 lines containing the values of (x,y,z) for the sample points in the 4 planes. The file Fmgn.dat has 33. 33 - 4 = 4356 lines containing the corresponding magnitudes V/E|x2 + \Ey\2 + \Ez12 of the electric field. Note that an expert user can easily alter the routine nearfielddist.f to compute whatever near field quantity might be of interest or to change the sampling rate. For far fields, the half space above the antenna is sampled in 36 x 16 = 576 far field points (36 0-values in [0~,90~] and 16 q-values in [0~,3600] with 0 and q being the traditional polar and azimuthal angles of a spherical coordinate system related to the (x, y, z)-coordinate system). The file ff. dat has 36 - 16 lines containing the values of 0 and ( along with a quantity proportional to /Ee\2 + IE^\2 in the far field. Note that an expert user can easily alter the routine farfielddist.f to compute whatever far field quantity might be of interest or to change the sampling rate. The file pol.dat contains information about the polarization of the far field normal to the plane of the antenna. It gives the axial ratio (in dB) of the polarization ellipse, the phase 8

difference (in degrees) between the two components and the tilt angle (in degrees) of the polarization ellipse measured from the x-axis. The format of this section is the following Text I: Near field code Text I: Near field code 9h -------- - - - - - - - - - - - - - - - - - - ---- Figure 1: Illustration of cavity sampling. Section 12 This section describes file names of output files containing geometrical parameters, input impedance information (the real and imaginary part of the input impedance for each frequency), convergence information (the number of iterations and the final relative residual for each frequency) and timing information (the time for the FEM part, BI part and solver for each frequency as well as time for pre-processing and "the rest of the code"). The format of this section is the following: Text "filename": Geometry file Text "filename": Impedance file Text "filename": Convergence file Text "filename": Timing file 4 Examples Consider a square metallic patch antenna backed by a rectangular cavity recessed in an infinite metallic ground plane, as illustrated in Fig. 2 (side view) and Fig. 3 (top view). The cavity-backed patch antenna is situated in free space characterized by the permittivity Co and the permeability [o. The cavity is of dimensions 1.85 cm x 1.85 cm x 0.15 cm and filled with a dielectric material of permittivity 10 Eo and conductivity 0.0003 S/cm. The patch is of side length 0.925 cm and centered in the cavity aperture. It is fed by a vertical coaxial 9

line whose outer conductor is attached to the ground plane and whose inner conductor is attached to the patch at the mid point of an edge, as illustrated in Fig. 2 and Fig. 3. The coaxial feed will be modeled as a vertical probe of constant current. Co go 0.15 cm 10~o do 0.0003 S/cm Figure 2: Side view of square metallic patch antenna backed by a dielectric-filled rectangular cavity recessed in an infinite metallic ground plane. Figure 3: Top view of square metallic patch antenna backed by a dielectric-filled rectangular cavity recessed in an infinite metallic ground plane. A coarse surface mesh is given in Fig. 4. For mixed-order TVFEs of order 0.5 and 1.5, the particular mesh is too coarse to yield the exact resonant frequency of 4.43GHz as obtained by Schuster and Luebbers [7] and confirmed by Andersen and Volakis for finer meshes [5]. Nevertheless, the mesh is very useful for illustrating the capabilities of MR-TETRA. Let us analyze the above antenna at 21 frequency points in the interval [3.5GHz,4.5GHz]. We will solve resulting linear equation system using a QMR solver with tolerance 103. AIM will not be used. The input file is given below. Name of universal file "mesh. unv" PEC code 1 (x,y,z)_start, (x,y,z)_stop, (xpa,y_pa)_start, (x_pa,y_pa)_stop -0.4625d0 -0.4625d0 10

1.2 - _ _ _/ _ _ _ _ __ _ _ _ _ 0.8 - l l _l_ / l l/ /I 0.6 - 0.2 - _- -__\_ _ -__\ -0.2 - 1/l/l/l/l/l/l/l/l -0.4 -_____________________ -0.5 0 0.5 1 1.5 Figure 4: Top view of coarse surface mesh for square metallic patch antenna backed by a dielectric-filled rectangular cavity recessed in an infinite metallic ground plane. -0. 15d0 1.3875d0 1.3875d0 OdO OdO OdO 0.925d0 0.925d0 Probe code 1 Number of probes 1 (xpr,ypr,zpr)_start, (x.pr,ypr,zpr).stop for first probe OdO 0.4625d0 -0.15d0 OdO 0.4625d0 OdO HO code 0 TVFE code 1 Number of material groups 1 Parameters for the first group (10d0,-Od0) (OdO,-OdO) (OdO,-OdO) (OdO,-OdO) 11

(10dO,-OdO) (OdO,-OdO) (OdO,-OdO) (OdO,-OdO) (10dO,-OdO) 0.0003d0 0 OdO OdO Integration order: FEM, BInear, BI_far and farfield 4 4 1 4 Start and stop frequencies 3.50d9 4.50d9 Number of frequency points 21 Solvercode 1 Tolerance ld-3 AIM code 0 AIM parameters (xstart, ystart, x_stop, ystop, I, J, K) OdO OdO OdO OdO 0 0 0 Near field distribution code 0 Far field distribution code 0 Name of file with geometrical data "3D-FEM-BI.outl" Name of file with input impedances "3D-FEM-BI.impl" Name of file with convergence data "3D-FEM-BI.conl" Name of file with timing data "3D-FEM-BI.timl" 12

The output file containing the input impedance for each frequency is given below and also plotted in Fig 5. The antenna is seen to experience the expected resonant behavior. Also, the resonant frequency is seen to be around 3.97GHz which is much smaller than the true resonant frequency 4.43GHz. This is due to order TVFEs have been used. the coarse mesh and the fact that lowest 3500000000.0000 3550000000.0000 3600000000.0000 3650000000.0000 3700000000.0000 3750000000.0000 3800000000.0000 3850000000.0000 3900000000.0000 3950000000.0000 4000000000.0000 4050000000.0000 4100000000.0000 4150000000.0000 4200000000.0000 4250000000.0000 4300000000.0000 4350000000.0000 4400000000.0000 4450000000.0000 4500000000.0000 2.5588456293470 3.2429206914211 4.2286657863228 5.7163554782377 8.1003200176110 12.236883071922 20.294341985730 38.929881217035 94.338002704657 278.78777970206 248.05473459466 87.641136274513 39.664558360104 22.234424951724 14.241706452950 9.9540096442860 7.3944826938964 5.7400309760047 4.6063185864701 3.8071120800757 3.2156270709780 41.628451744725 45.778072533465 50.904979316993 57.440189012613 66.098467821245 78.132689480470 95.967863359927 124.44216366128 169.57056063719 152.06654344528 -138.59626984899 -133.27080997795 -92.517092077889 -66.678035812629 -50.108946578041 -38.749056086777 -30.493855847262 -24.197511427906 -19.211509191283 -15.235725559360 -11.848352913871 At the resonant frequency 3.97GHz, let us now look at the polarization of the far field normal to the antenna as well as the far field patterns. Hence, we change Section 8 and 11 of the input file to read Start and stop frequencies 3.97d9 3.97d9 Number of frequency points 1 Near field distribution code 0 Far field distribution code 3 The polarization of the far field normal to the antenna is described in the file pol. dat: 3970000000.0000 26.122202028302 -85.365890875233 0.23007593744868 13

Aflno -- TVFEcode= 1 300C- - - TVFE code=3........................ 0 0.............................................. 3.5 3.6 3.7 3.8 3.9 4 4.1 4.2 4.3 4.4 4.5 Frequency [GHz] 200 - 10 0......................... -i -100e|- TVFE code = 1 -200 -............ 3.5 3.6 3.7 3.8 3.9 4 4.1 4.2 4.3 4.4 4.5 Frequency [GHz] Figure 5: Real and imaginary part of the input impedance of the antenna in Fig. 2 and Fig. 3. This output implies that the polarization ellipse is tilted 0.23~ from the x-axis and has an axial ratio of 26.12dB. That is, we essentially have linear polarization along the x-axis which is what we expect for this patch. The patterns in the planes y = 0 and x = 0 at 3.97GHz are given in Fig. 6 (polar angles in [0~, 90~] on the plot corresponds to 4-values in [0(), 180~] while polar angles in [90~, 180~] on the plot corresponds to q-values in [1800, 3600]; also, the polar angle 90~ on the plot corresponds to 0 = 0~ and the polar angles 00 and 1800 correspond to 0 = 90~). These are consistent with the polarization of the far field normal to the antenna and exactly what we expect for this patch. 90 24.5054 90 24.5054 120 60 120 60 150/ 122527 30 150/ 122.27 2 30 180............. 0 180....... 0 210. ' " '330 210 ' /330 240 00 240 00 270 270 Figure 6: Patterns in the planes y = 0 and x = 0 at 3.97GHz for the antenna in Fig. 2 and Fig. 3. To demonstrate the use of higher order TVFEs for improved field modeling, we apply higher order TVFEs around the radiating edges of the patch where the field is known to be the strongest and experience the most variation. Hence, we change Section 4 and 5 of the input file to read 14

HO code 1 Number of higher order volumes 2 (x,y,z)_start, -0.23125d0 -0.23125d0 -0.15d0 0.23125d0 1.15625d0 OdO (x,y,z)-start, 0.69375d0 -0.23125d0 -0.15d0 1.15625d0 1.15625d0 OdO (x,y,z)_stop for first higher order volume (x,y,z)_stop for second higher order volume TVFE code 3 Let us consider the application of AIM and use a DFT to compute the matrix-vector products in the iterative solver. This is terribly inefficient for small problems like the one considered here and only serves the purpose of demonstrating how AIM is used. To invoke AIM, we change Section 10 of the input file to read AIM code 1 AIM parameters (xstart, ystart, xstop, ystop, I, J, K) -0.668055555555d0 -0.668055555555d0 1.5930555555555d0 1.5930555555555d0 23 23 8 The output file containing the input impedance for each frequency is given below and also plotted in Fig 5. The antenna is again seen to experience the expected resonant behavior. The resonant frequency is seen to be 4.32GHz which is much closer to the true resonant frequency 4.43GHz than the 3.97GHz found when lowest order TVFEs were used throughout the cavity. 3500000000.0000 3550000000.0000 3600000000.0000 1.0325245571199 1.1787602915226 1.3630059719665 37.824886278503 39.744351943887 41.855292671137 15

3650000000.0000 3700000000.0000 3750000000.0000 3800000000.0000 3850000000.0000 3900000000.0000 3950000000.0000 4000000000.0000 4050000000.0000 4100000000.0000 4150000000.0000 4200000000.0000 4250000000.0000 4300000000.0000 4350000000.0000 4400000000.0000 4450000000.0000 4500000000.0000 1.5908198292434 1.8827809279585 2.2511366594008 2.7452949855043 3.4128718362669 4.3376897449660 5.6748167344092 7.6934532336221 10.930495769666 16.607146734682 27.615471055343 52.159382628536 122.49555861383 320.78659000805 267.24188214807 103.160187785883 48.490909355324 27.621340718090 44.211407702366 46.855809984373 49.890077408460 53.404271011424 57.554844531757 62.542584004342 68.707629601590 76.555049192618 86.945891370685 101.292500931489 122.03417410905 154.29204802524 200.49906174576 152.29837557693 -137.93155070714 -137.62354804365 -96.045352685244 -67.895949125991 5 Summary This manual describes how to run MR-TETRA, a multi-resolution finite element / boundary integral code for analysis of printed antennas backed by material-filled (dielectric or ferrite) metallic cavities recessed in infinite metallic ground planes. Both a general description as well as examples are given. References [1] J.L. Volakis, A. Chatterjee, and L.C. Kempel. Finite element method for electromagnetics. IEEE Press, USA, 1998. [2] L.S. Andersen and J.L. Volakis. 'Hierarchical tangential vector finite elements for tetrahedra'. IEEE Microwave and Guided Wave Letters, vol. 8, pp. 127-129, March 1998. [3] E. Bleszynski, M. Bleszynski, and T. Jaroszewicz. 'AIM: Adaptive integral method for solving large-scale electromagnetic scattering and radiation problems'. Radio Sci., vol. 31, pp. 1225-1251, 1996. [4] Y. Saad. Iterative methods for sparse linear systems. PWS Publishing Company, 1996. [5] L.S. Andersen and J.L. Volakis. 'Accurate and efficient simulation of antennas using hierarchical mixed-order tangential vector finite elements for tetrahedra'. Accepted for publication in IEEE Trans. Antennas Propagat. [6] D.M. Pozar. Microwave engineering. Addison-Wesley Publishing Company, Inc., 1990. [7] J.W. Schuster and R. Luebbers. Private communication. 16